Stirred Mill (Perfect Mixing): Difference between revisions
Scott Munro (talk | contribs) m (1 revision imported) |
imported>Scott.Munro m (→Results) |
||
Line 50: | Line 50: | ||
Adopting the same scaling approach and nomenclature as the [[Ball Mill (Perfect Mixing)|Perfect Mixing ball mill]], the <math>R_i/D_i^*</math> rate is scaled by the following relation: | Adopting the same scaling approach and nomenclature as the [[Ball Mill (Perfect Mixing)|Perfect Mixing ball mill]], the <math>R_i/D_i^*</math> rate is scaled by the following relation: | ||
:<math>\left(\frac{R}{D^*}\right)_{\rm Sim} = \left(\frac{R}{D^*}\right)_{\rm Orig} \cdot | :<math>\left(\frac{R}{D^*}\right)_{\rm Sim} = \left(\frac{R}{D^*}\right)_{\rm Orig} \cdot f_{\rm Ds} \cdot f_{\rm Hb} \cdot f_{\rm Ns} \cdot f_{\rm S} \cdot f_{\rm T} \cdot f_{\rm Db} \cdot f_{\rm WI}</math> | ||
{{Model theory (Text, Stirred Mill, Perfect Mixing, Breakage Scaling)}} | {{Model theory (Text, Stirred Mill, Perfect Mixing, Breakage Scaling)}} | ||
Line 60: | Line 60: | ||
=== Internal mesh series === | === Internal mesh series === | ||
{{Model theory (Text, Ball Mill, Perfect Mixing, Internal mesh)}} | {{Model theory (Text, Ball Mill, Perfect Mixing, Internal mesh)|Perfect Mixing}} | ||
=== Multi-component modelling === | === Multi-component modelling === | ||
{{Model theory (Text, Ball Mill, Perfect Mixing, Multi-component)}} | {{Model theory (Text, Ball Mill, Perfect Mixing, Multi-component)|Perfect Mixing|breakage rate}} | ||
== Excel == | == Excel == | ||
Line 176: | Line 176: | ||
\tau\text{ (min)}\\ | \tau\text{ (min)}\\ | ||
\text{R/D* factor (frac)}\\ | \text{R/D* factor (frac)}\\ | ||
f_{\rm Ds}\text{ (-)}\\ | |||
f_{\rm Hb}\text{ (-)}\\ | |||
f_{\rm Ns}\text{ (-)}\\ | |||
f_{\rm S}\text{ (-)}\\ | |||
f_{\rm T}\text{ (-)}\\ | |||
f_{\rm Db}\text{ (-)}\\ | |||
f_{\rm WI}\text{ (-)}\\ | |||
\end{bmatrix} | \end{bmatrix} | ||
Revision as of 06:57, 14 June 2023
Description
This article describes an implementation of the Perfect Mixing vertical stirred mill model outlined by Napier-Munn et al. (1996).[1]
The model described here is for steady-state simulation. For dynamic simulation, see Stirred Mill (Perfect Mixing, Dynamic).
Model theory
The Perfect Mixing model is based on a population balance of particles entering the mill, breaking into smaller sizes, and discharging as product. For a mill operating in steady-state, the diagram in Figure 1 below represents the balance for a given size fraction:
The steady-state population balance is formulated mathematically as:
where:
- is the index of the size interval, , is the number of size intervals
- is the mass flow rate of solids of size interval in the mill feed
- is the mass flow rate of solids of size interval in the mill product
- is the mass of solids of size interval in the mill load
- is the breakage rate of solids of size interval in the mill load
- is the appearance function, the distribution of particle mass arising from the breakage of a parent particle in size interval into progeny of size interval
As the mill is perfectly mixed, the product is related to the mill contents and discharge rate as:
where is the rate of discharge of solids in size interval from the mill.
Substitution and rearrangement of the above equations leads to:
The mill product can therefore be computed provided a feed rate, appearance function and breakage rate per discharge rate, , is available. Alternatively, the function can be determined from the feed rate, product rate and appearance function.
Discharge rate
The discharge rate, , is related to the residence time of pulp in the mill and can be corrected for different volumetric feed rates and mill volumes by first normalising to
where:
- is the residence time normalised discharge rate per size interval
- is the volume of slurry in the mill
- is the volumetric flow rate of slurry in mill feed (solids plus liquids)
- is the mean residence time of slurry in the mill
The volume of slurry in the mill, (m3), is the internal volume of the mill, (m3), minus the volume of the media, (m3),:
where:
- is the mill diameter (m)
- is the mill height, the distance from the feed inlet to the product outlet (m)
- is the height of the media charge (m)
- is the volume occupied by the stirrer mechanism between the feed inlet to the product outlet (m3)
- is volumetric fraction of interstitial void space in the charge, assumed here to be 0.4 (v/v)
The volume occupied by the stirrer mechanism is ignored as a simplification.
The actual discharge rate for a given mill and pulp feed rate is subsequently scaled as:
with becoming for mill model calibration and simulation.
Breakage rate
The rate is a function of particle size and is typically a smooth curve. In order to reduce the number of model parameters, the rates are specified only at three or four regularly spaced sizes (knots). Cubic spline interpolation is then used to reconstruct the complete set of rates for all size intervals.
The breakage rate, , is affected by mill operating conditions as follows:
- is the stirrer diameter (m)
- is the stirrer speed (rpm)
- is the number of starts of the screw stirrer
- is the number of turns per start along the full height of the screw stirrer ()
- is the mean media diameter (mm)
Adopting the same scaling approach and nomenclature as the Perfect Mixing ball mill, the rate is scaled by the following relation:
The scaling factors are defined as:
where:
- the subscript refers to the original mill from which was derived
- the subscript refers to the mill being simulated (scaled)
- are scaling exponents, , , and
The Work Index scaling factor, , is retained from the Perfect Mixing ball mill model, where is the Bond Ball Work Index value of the ore (kWh/t).
The model user may optionally change the values of to suit specific test work results or operating data.
Appearance function
The Appearance function describes the mass-by-size distribution of progeny particles resulting from the breakage of parent particles.
The Appearance function may be specified for a particular ore.
Morrell et al. (1993) found the appearance functions generated for the Perfect Mixing ball mill model were too fine for stirred mills. They recommended either generating new appearance functions at one-tenth of the standard breakage test input energy level, or using the Kelsall breakage function.[2]
Internal mesh series
The Perfect Mixing mill model is formulated internally with a geometric progression of 31 mesh sizes at intervals. Feed and product size fractions are automatically converted to and from the internal mesh series during model computation. The size intervals allow the appearance function to be specified as a one-dimensional matrix, rather than the two dimensional form defined above, since
when the intervals are so spaced.
Multi-component modelling
The original Perfect Mixing model formulation only considered the properties of a single ore type.
This implementation applies different appearance functions and breakage rate scaling factors to separate population balance computations for each ore type in the feed.
Excel
The Perfect Mixing stirred mill model may be invoked from the Excel formula bar with the following function call:
=mdUnit_StirredMill_PerfectMixing(Parameters as Range, Size as Range, MillFeed as Range, OreSG As Range, Appearance as Range, WorkIndexSim As Range, R/D*KnotPositions as Range, R/D*KnotsOrig as Range)
Invoking the function with no arguments will print Help text associated with the model, including a link to this page.
Inputs
The required inputs are defined below in matrix notation with elements corresponding to cells in Excel row () x column () format:
where:
- is the mass flow feed rate of liquids into the mill (t/h)
- is the Specific Gravity or density of the media in the mill (- or t/m3)
- is the Specific Gravity or density of liquids in the feed (- or t/m3)
- is the number of ore types
- is the number of breakage rate per discharge rate knots
- is the size of the square mesh interval that feed mass is retained on (mm)
- , i.e. descending size order from top size () to sub mesh ( mm)
- is the mass flow rate of particles in the feed (t/h)
- is the Specific Gravity or density of solids (- or t/m3)
Results
The results are displayed in Excel as an array corresponding to the matrix notation below:
where:
- is the media load fraction (v/v), defined as the fraction of mill volume occupied by media, including void spaces, i.e.
- is the discharge rate scaling factor
- is the mass flow rate of particles in the mill product (t/h)
- is the geometric mean size of the internal mesh series interval that mass is retained on (mm)
Example
The images below show the selection of input arrays and output results in the Excel interface.
SysCAD
The SysCAD interface for steady-state (ProBal) mode is described below. For SysCAD Dynamic, see Stirred Mill (Perfect Mixing, Dynamic).
MD_Mill page
The first tab page in the access window will have this name.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
Tag | Display | This name tag may be modified with the change tag option. |
Condition | Display | OK if no errors/warnings, otherwise lists errors/warnings. |
ConditionCount | Display | The current number of errors/warnings. If condition is OK, returns 0. |
GeneralDescription / GenDesc | Display | This is an automatically generated description for the unit. If the user has entered text in the 'EqpDesc' field on the Info tab (see below), this will be displayed here.
If this field is blank, then SysCAD will display the unit class ID. |
Requirements | ||
On | CheckBox | This enables the unit. If this box is not checked, then the material will pass straight through the mill with no change to the size distribution. |
Method | Fixed Discharge | The discharge particle size distribution is user defined. Different distributions can be used for different solids. |
AG/SAG (Variable Rates) | The Variable Rates AG/SAG mill model is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
Rod Mill (Lynch) | The Lynch rod mill model is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
Ball (Perfect Mixing) | The Perfect Mixing ball mill model (steady-state or dynamic) is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
Stirred (Perfect Mixing) | The Perfect Mixing stirred mill model (steady-state or dynamic) is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
Mill (Herbst-Fuerstenau) | The Herbst-Fuerstenau model is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
PowerModels | CheckBox | Show alternative mill power model calculations on the Power page. |
MediaStrings | CheckBox | Show media size distributions at recharge equilibrium on the MediaStrings page. |
Options | ||
ShowQFeed | CheckBox | QFeed and associated tab pages (eg Sp) will become visible, showing the properties of the combined feed stream. |
ShowQProd | CheckBox | QProd and associated tab pages (eg Sp) will become visible, showing the properties of the products. |
SizeForPassingFracCalc | Input | Size fraction for % Passing calculation. The size fraction input here will be shown in the Stream Summary section. |
FracForPassingSizeCalc | Input | Fraction passing for Size calculation. The fraction input here will be shown in the Stream Summary section. |
Stream Summary | ||
MassFlow / Qm | Display | The total mass flow in each stream. |
SolidMassFlow / SQm | Display | The Solids mass flow in each stream. |
LiquidMassFlow / LQm | Display | The Liquid mass flow in each stream. |
VolFlow / Qv | Display | The total Volume flow in each stream. |
Temperature / T | Display | The Temperature of each stream. |
Density / Rho | Display | The Density of each stream. |
SolidFrac / Sf | Display | The Solid Fraction in each stream. |
LiquidFrac / Lf | Display | The Liquid Fraction in each stream. |
Passing | Display | The mass fraction passing the user-specified size (in the field SizeForPassingFracCalc) in each stream. |
Passes | Display | The user-specified (in the field FracForPassesSizeCalc) fraction of material in each stream will pass this size fraction. |
Mill page
The Mill page is used to specify the input parameters for the mill model.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
PerfectMixing | ||
HelpLink | Opens a link to this page using the system default web browser. Note: Internet access is required. | |
Requirements | ||
NumParallelUnits | Input | The number of parallel, identical units to simulate:
|
MediaStringsP50 | CheckBox |
|
Stirred | ||
UserExponents | Checkbox | Indicates whether to use user-defined values for the scaling exponents. |
MillDiameter | Input | The inside liner diameter of the simulated mill. |
MillHeight | Input | The height of the simulated mill. |
MediaHeight | Input | Height of the media charge of the original and simulated mills. |
MediaDiameter | Input | Median diameter of the media in the original and simulated mills. |
StirrerDiameter | Input | Diameter of the stirrers of the original and simulated mills. |
StirrerSpeed | Input | Rotational speed of the stirrers of the original and simulated mills. |
StirrerStarts | Input | Number of starts of the (screw) stirrers of the original and simulated mills. |
StirrerTurns | Input | Number of turns per start of the (screw) stirrers of the original and simulated mills. |
WorkIndex | Input | Bond Ball Work Index of ore in the original mill. |
Exponent | Input |
|
R/DFunction | ||
NumSplineKnots | Input | Number of spline knots for the function. |
Size | Input | Spline knot size positions. |
Ln(R/D*) | Input | Values of at each spline knot position. |
Results | ||
MillVolume | Display | Volume inside the mill. |
MediaVolume | Display | Volume occupied by the media in the mill, excluding void space. |
SlurryVolume | Display | Volume of slurry in the mill, including slurry in media void space and above media charge. |
LoadFraction | Display | Volume of media charge, including void space and excluding stirrer volume, as a fraction of mill volume, excluding stirrer volume. |
ResidenceTime | Display | Mean residence time of slurry in the mill. |
Ore page
This page is used to define the comminution properties of SysCAD species with the size distribution quality in the project.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
Appearance | ||
DefaultAppearance | Sets all species to the the default Kelsall Appearance function. | |
OreSpecific | CheckBox |
|
Appearance | Input | User-specified Appearance function data for all species with size distribution property. |
WorkIndex | ||
WorkIndex.Sim | Input | Bond Ball Work Index data for all species with size distribution property. |
Ri/Di page
This page displays the scaling factors and breakage rate per discharge rate for each size interval computed by the Perfect Mixing model.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
Scaling | ||
RDStar/RD | Display | Value of the factor for discharge rate scaling. |
StirrerDiameter | Display | Value of the stirrer diameter factor for rate scaling. |
MediaHeight | Display | Value of the media height factor for rate scaling. |
StirrerSpeed | Display | Value of the stirrer speed factor for rate scaling. |
StirrerStarts | Display | Value of the stirrer starts factor for rate scaling. |
StirrerTurns | Display | Value of the stirrer turns per start factor for rate scaling. |
MediaSize | Display | Value of the media size factor for rate scaling. |
WorkIndex | Display | Value of the Work Index factor of each ore species for rate scaling. |
Ri/DiStar | ||
Size | Display | Size of each interval in internal mesh series. |
MeanSize | Display | Geometric mean size of each interval in internal mesh series. |
Ri/DiStar | Display | Value of normalised rate for each size interval, for each ore species. |
Ri/Di | ||
Size | Display | Size of each interval in internal mesh series. |
MeanSize | Display | Geometric mean size of each interval in internal mesh series. |
Ri/Di | Display | Value of rate for each size interval, for each ore species. |
Power page
This optional page displays the inputs and results for alternative mill power models. The page is only visible if PowerModels is selected on the MetDynamics*Mill page.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
Power | ||
Nitta | CheckBox | Shows inputs and results for stirred mill power calculations using the Nitta method. |
Heath | CheckBox | Shows inputs and results for stirred mill power calculations using the Heath method. |
MediaStrings page
This page displays the inputs and results for grinding mill media string calculations. The page is only visible if MediaStrings is selected on the MetDynamics*Mill page.
About page
This page is provides product and licensing information about the Met Dynamics Models SysCAD Add-On.
Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
---|---|---|
About | ||
HelpLink | Opens a link to the Installation and Licensing page using the system default web browser. Note: Internet access is required. | |
Information | Copies Product and License information to the Windows clipboard. | |
Product | ||
Name | Display | Met Dynamics software product name |
Version | Display | Met Dynamics software product version number. |
BuildDate | Display | Build date and time of the Met Dynamics Models SysCAD Add-On. |
License | ||
File | This is used to locate a Met Dynamics software license file. | |
Location | Display | Type of Met Dynamics software license or file name and path of license file. |
SiteCode | Display | Unique machine identifier for license authorisation. |
ReqdAuth | Display | Authorisation level required, MD-SysCAD Full or MD-SysCAD Runtime. |
Status | Display | License status, LICENSE_OK indicates a valid license, other messages report licensing errors. |
IssuedTo | Display | Only visible if Met Dynamics license file is used. Name of organisation/seat the license is authorised to. |
ExpiryDate | Display | Only visible if Met Dynamics license file is used. License expiry date. |
DaysLeft | Display | Only visible if Met Dynamics license file is used. Days left before the license expires. |
See also
References
- ↑ 1.0 1.1 Napier-Munn, T.J., Morrell, S., Morrison, R.D. and Kojovic, T., 1996. Mineral comminution circuits: their operation and optimisation. Julius Kruttschnitt Mineral Research Centre, Indooroopilly, QLD.
- ↑ 2.0 2.1 Morrell, S., Sterns, U.J. and Weller, K.R., 1993. The application of population balance models to very fine grinding in tower mills. 18th International Mineral Processing Congress; 23-28 May, 1993; Sydney, NSW. 1993, pp. 61-66.