Tumbling Mill (Power, Hilden and Powell)

From Met Dynamics
Revision as of 13:21, 2 March 2023 by Scott Munro (talk | contribs) (1 revision imported)
Jump to navigation Jump to search

Description

This article describes the Hilden and Powell (2018) approach for estimating the power draw of a tumbling mill.[1]

Hilden and Powell extended Morrell's (1996, 2016) approaches to include the following features:[2][3]

  • Partial filling of the charge void space with slurry in two directions, from the charge shoulder towards the toe and from the outer shell towards the mill centre
  • Estimation of the fraction of a slurry pool which settles across the centre of the mill profile and thus does not contribute to power draw
  • No-load power based on mill drive type and size

Model theory

Hilden and Powell's approach adopts Morrell's (1996) expression for tumbling mill power draw:[2]

where is a lumped calibration parameter accounting for power losses.

Hilden and Powell expand Morrell's expression to:

where:

  • is the gross power draw of the mill (kW)
  • is the no-load power draw of the mill (kW)
  • is the power drawn by the charge component in the cylindrical section of the mill (kW)
  • is the power drawn by the slurry component in the cylindrical section of the mill (kW)
  • is the power drawn by the charge component in the cone end sections of the mill (kW)
  • is the power drawn by the slurry component in the cone end sections of the mill (kW)

Hilden and Powell correct Morrell's value of of 1.26 to a value of 1.25 based on analysis of an updated data set with the model relations outlined below.

Sub-components of the gross power equation are described below.

Charge motion power

Figure 1. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when , along with key dimensions.
Figure 2. Schematic of a tumbling mill showing key dimensions.

Using Morrell's (1996) equations, the power drawn by the charge component in the cylindrical and cone end sections of the mill, and , are:[2]

where:

  • is length of the cylindrical section (belly) of the mill inside liners (m)
  • is length of the cone end (m)
  • is acceleration due to gravity (m/s2)
  • is the rotational rate of the mill (rev/s)
  • is the density of the total charge (t/m3)
  • is the mill radius inside liners (m), which is half the mill diameter , i.e.
  • is the trunnion radius inside liners (m), which is half the trunnion diameter , i.e.
  • is the inner surface radius of the charge (m)
  • is the angle of the charge shoulder (radians)
  • is the angle of the charge toe (radians)
  • is an empirical term that relates the rotational rate of mass at radial positions within the charge to the rotational rate of the mill

and all angles are unit circle positive.

The length of the cone end, (m), is:

where is the angular displacement of the cone surface from the vertical direction (rad).

The parameter is:

where is the total mill filling of balls plus rocks (v/v).

The power drawn by the slurry component in the cylindrical and cone end sections, and , is computed by substituting the slurry component properties , , , and for the equivalent charge properties , , , and in the same equations for and , where:

  • is the density of the slurry component in the mill (t/m3)
  • is the inner surface radius of the slurry component (m)
  • is an empirical term that relates the rotational rate of mass at radial positions within the slurry component to the rotational rate of the mill
  • is the angle of the toe of the slurry component (radians)

The parameter is:

where is the apparent filling for slurry (v/v), i.e.:

and is defined as the volume fraction of the charge void space occupied by slurry (v/v).

The remaining parameters of the charge motion power equations are described below.

Charge position

The charge shoulder and toe positions and inner surface radius are calculated using Morrell's (1996) method.[2]

The position of the toe of the charge, (rad), is defined as:

where (frac) is the theoretical fraction critical speed, and the fraction of critical speed at which centrifuging actually occurs, (frac), is:

The position of the shoulder of the charge, (rad), is:

The inner surface radius of the charge, (m) is:

where the active load fraction, (frac), is:

Slurry position

The position of the shoulder of the slurry component is the same as the charge shoulder position.

The position of the toe of the slurry component varies based on the filling of the charge void space, :

  • if then the charge void space is partially or exactly filled with slurry and
  • if then there is an excess of slurry beyond the available charge void space, which manifests as a slurry pool, and .

Slurry toe when U≤1

When the charge void space is partially or exactly filled with slurry (), the position of the toe of the slurry component, (rad), is:

where is the profile shape factor (frac).

The inner surface radius of the slurry component, (m), is:

Slurry toe when U>1

Figure 3. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when , along with key dimensions and areas.
Figure 4. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when , and , along with key dimensions and areas.

When the charge void space is overfilled (), a slurry pool forms.

Hilden and Powell conceptualised the slurry pool as occupying two areas within the mill, as illustrated in Figure 3:

  • in the annular sector at the toe of the charge (), and
  • in the circular segment settled across the centre of the mill ().

Only the slurry in the annular sector () contributes to power draw.

The inner surface radius of the slurry component filling the charge void space and extending into the annular sector is equal to the charge radius, i.e. .

Locating the position of the toe of the slurry component, , then becomes a problem of determining the distribution of slurry between the two areas, and .

The total area of excess slurry, (m2), in the pool is:

The area of the annular segment, , is:

where the parameters and , as shown in Figure 4, are:

The area of the circular sector, , is:

The slurry toe parameter cannot be algebraically isolated in the above equations. Therefore, an iterative numerical method is used to solve the value of that yields the required value of (in this case, Ridders' method).

Charge density

The masses of balls, (t), and coarse rocks, (t), in the charge are:

where:

  • is the ball filling (v/v)
  • is the volume of of the cylindrical section of the mill (m3)
  • is the charge void fraction, typically 0.4 (v/v)
  • is the Specific Gravity or density of the balls (- or t/m3)
  • is the Specific Gravity or density of the rock in the charge (- or t/m3)

The volume of of the cylindrical section of the mill () is:

The charge density, (t/m3), excluding slurry, is then:

The apparent density of slurry within the charge, (t/m3), is:

where is the density of slurry (t/m3). The value of is computed as:

where is the density of liquids (t/m3) and is the volume fraction of solids in the slurry (v/v). is calculated from the mass fraction of solids in the slurry, (w/w), as:

No-load power

The no-load power () is estimated as:

where is the maximum power draw of the installed drive motor.

The motor loss factor may be obtained from drive vendor specifications or selected from the table below:

Table 1. Typical motor loss factors (after Hilden and Powell (2018))[1]
Drive Type Motor loss factor
Range (%) Average (%)
Gearless GMD 3.2 - 4.9 3.9
VSD LS synchronous 4.36 - 5.26 4.8
VSD with SER 5.79 - 7.11 6.4
VSD HS asynchronous 5.6 - 6.62 6.1

Excel

The Hilden and Powell mill power model may be invoked from the Excel formula bar with the following function call:

=mdMillPower_HildenPowell(Parameters as Range)

Invoking the function with no arguments will print Help text associated with the model, including a link to this page.

The Parameters array and model results are defined below in matrix notation, along with an example image showing the selection of the same arrays in the Excel interface:


where:

Figure 5. Example showing the selection of the Parameters (blue frame), and Results (light blue frame) arrays in Excel.

SysCAD

The Hilden and Powell power model is an optional calculation for tumbling mill units. If selected, the input and display parameters below are shown.

Tag (Long/Short) Input / Display Description/Calculated Variables/Options
HildenPowell
HelpLink ButtonModelHelp.png Opens a link to this page using the system default web browser. Note: Internet access is required.
MillDiameter Display Diameter of the mill (inside liners).
BellyLength Display Length of the cylindrical section (belly) of the mill (inside liners).
TrunnionDiameter Display Diameter of the trunnion (inside liners).
ConeAngle Display Angular displacement of the cone surface from the vertical direction.
Jt Display Volumetric fraction of the mill occupied by balls and coarse rock (including voids).
BallLoadVol Display Volumetric fraction of the mill occupied by balls (including voids).
SolidsSG Display Specific Gravity or density of solids.
BallSG Display Specific Gravity or density of balls.
LiquidsSG Display Specific Gravity or density of liquids.
Voidage Input/Display Volumetric fraction of interstitial void space in the charge. Usually 0.4.
Note: The Dynamic Perfect Mixing ball mill model uses Shi's dynamic charge porosity estimation to align power draw predictions with the adopted slurry filling approach.
U Display Volumetric fraction of interstitial grinding media voidage occupied by slurry.
Cw Display Mass fraction of solids in discharge slurry.
FracCS Display Fraction critical speed of the mill.
MillDriveLosses Input Mill drive motor power loss factor.
MotorInstalledPower Input Maximum power draw of installed mill motor.
NetPowerAdjust Input Lumped calibration parameter accounting for power losses. Found to be a value of 1.26 for industrial mills.
y Input Profile shape factor.
ThetaShoulder Display Angle of the charge shoulder.
ThetaToe Display Angle of the charge toe.
ThetaSlurryToe Display Angle of the toes of the slurry component
ChargeSurfaceRadius Display Radius of the inner surface of the charge.
SlurrySurfaceRadius Display Radius of the inner surface of the slurry component.
NoLoadPower Display Power input to the motor when the mill is empty (no balls, rocks or slurry).
NetPower Display Charge motion power, including losses.
GrossPower Display Power input to the motor.

See also

References

  1. 1.0 1.1 Hilden, M.M. and Powell, M.S., 2018. A model of SAG mill power incorporating slurry filling and its application to real-time mill performance analysis. Paper presented at Comminution '18, Cape Town, South Africa, April 16-19, 2018.
  2. 2.0 2.1 2.2 2.3 Morrell, S., 1996. Power draw of wet tumbling mills and its relationship to charge dynamics. Pt. 1: a continuum approach to mathematical modelling of mill power draw. Transactions of the Institution of Mining and Metallurgy. Section C. Mineral Processing and Extractive Metallurgy, 105.
  3. Morrell, S., 2016. Modelling the influence on power draw of the slurry phase in Autogenous (AG), Semi-autogenous (SAG) and ball mills. Minerals Engineering, 89, pp.148-156.