View source for Tumbling Mill (Power, Hilden and Powell)

== Description ==

This article describes the '''Hilden and Powell''' (2018) approach for estimating the power draw of a tumbling mill.{{Hilden_and_Powell_(2018)}}

Hilden and Powell extended Morrell's ([[Tumbling Mill (Power, Morrell Continuum)|1996]], 2016) approaches to include the following features:{{Morrell (1996a)}}{{Morrell (2016)}}

* 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:{{Morrell (1996a)}}

:<math>\text{Gross power (kW)} = \text{No-load power} + (k \times \text{Charge motion power})</math>

where <math>k</math> is a lumped calibration parameter accounting for power losses.

Hilden and Powell expand Morrell's expression to:

:<math>P = P_{\rm N} + k.(P_{\rm L} + P'_{\rm L} + P_{\rm C} + P'_{\rm C})</math>

where:
* <math>P</math> is the gross power draw of the mill (kW)
* <math>P_{\rm N}</math> is the no-load power draw of the mill (kW)
* <math>P_{\rm L}</math> is the power drawn by the charge component in the cylindrical section of the mill (kW)
* <math>P'_{\rm L}</math> is the power drawn by the slurry component in the cylindrical section of the mill (kW)
* <math>P_{\rm C}</math> is the power drawn by the charge component in the cone end sections of the mill (kW)
* <math>P'_{\rm C}</math> 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 <math>k</math> 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 ===

[[File:TumblingMillPower8.png|thumb|450px|Figure 1. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when <math>U<1</math>, along with key dimensions.]]

[[File:TumblingMillDimensions3.png|thumb|450px|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, <math>P_{\rm L}</math> and <math>P_{\rm C}</math>, are:{{Morrell (1996a)}}

:<math>P_{\rm L} = \frac{\pi g L \rho_{\rm c} N_{\rm m} R}{3(R - zr_{\rm i})} \left [ 2R^3 -3zR^2r_{\rm i} + r_{\rm i}^3(3z-2) \right ].(\sin \theta_{\rm S} - \sin \theta_{\rm T}) + L \rho_{\rm c} \left ( \frac{\pi N_{\rm m} R}{R - zr_{\rm i}} \right )^3 \left [ (R-zr_{\rm i})^4 - r_{\rm i}^4(z-1)^4 \right ] </math>

:<math>P_{\rm C} = \frac{\pi g L_{\rm d} \rho_{\rm c} N_{\rm m}}{3(R-r_{\rm T})}(R^4 - 4Rr_{\rm i}^3 + 3r_{\rm i}^4)(\sin \theta_{\rm S} - \sin \theta_{\rm T}) + \frac{2 \pi^3 N_{\rm m}^3 L_{\rm d} \rho_{\rm c}}{5(R-r_{\rm T})}(R^5 - 5Rr_{\rm i}^4 + 4r_{\rm i}^5)</math>

where:
* <math>L</math> is length of the cylindrical section (belly) of the mill inside liners (m)
* <math>L_{\rm d}</math> is length of the cone end (m)
* <math>g</math> is acceleration due to gravity (m/s<sup>2</sup>)
* <math>N_{\rm m}</math> is the rotational rate of the mill (rev/s)
* <math>\rho_{\rm c}</math> is the density of the total charge (t/m<sup>3</sup>)
* <math>R</math> is the mill radius inside liners (m), which is half the mill diameter <math>D</math>, i.e. <math>R=0.5D</math>
* <math>r_{\rm T}</math> is the trunnion radius inside liners (m), which is half the trunnion diameter <math>d_{\rm T}</math>, i.e. <math>r_{\rm T}=0.5d_{\rm T}</math>
* <math>r_{\rm i}</math> is the inner surface radius of the charge (m)
* <math>\theta_{\rm S}</math> is the angle of the charge shoulder (radians)
* <math>\theta_{\rm T}</math> is the angle of the charge toe (radians)
* <math>z</math> 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, <math>L_{\rm d}</math> (m), is:

:<math>L_{\rm d} = (R - r_{\rm T}) \tan \alpha_{c}</math>

where <math>\alpha_{\rm c}</math> is the angular displacement of the cone surface from the vertical direction (rad).

The parameter <math>z</math> is:

:<math>z = \left (  1 - J_{\rm T} \right )^{0.4532}</math>

where <math>J_{\rm T}</math> 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, <math>P'_{\rm L}</math> and <math>P'_{\rm C}</math>, is computed by substituting the slurry component properties <math>\rho'_c</math>, <math>r'_{\rm i}</math>, <math>z'</math>, and <math>\theta_{\rm L}</math> for the equivalent charge properties <math>\rho_{\rm c}</math>, <math>r_{\rm i}</math>, <math>z</math>, and <math>\theta_{\rm T}</math> in the same equations for <math>P_{\rm L}</math> and <math>P_{\rm C}</math>, where:
* <math>\rho'_c</math> is the density of the slurry component in the mill (t/m<sup>3</sup>)
* <math>r'_{\rm i}</math> is the inner surface radius of the slurry component (m)
* <math>z'</math> 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
* <math>\theta_{\rm L}</math> is the angle of the toe of the slurry component (radians)

The parameter <math>z'</math> is:

:<math>z' = \left (  1 - J'_{\rm T} \right )^{0.4532}</math>

where <math>J'_{\rm T}</math> is the apparent filling for slurry (v/v), i.e.:

:<math>
J'_{\rm T} =
\begin{cases}
U.J_{\rm T} & U \leq 1\\
J_{\rm T} + (U - 1) . U . \varepsilon & U>1
\end{cases}
</math>

and <math>U</math> 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.{{Morrell (1996a)}}

The position of the '''toe of the charge''', <math>\theta_{\rm T}</math> (rad), is defined as: 

:<math>\theta_{\rm T} = 2.5307 \left ( 1.2796 - J_{\rm t} \right ) \left (1 - {\rm e}^{-19.42(\phi_{\rm c} - \phi)} \right ) + \frac{\pi}{2}</math>

where <math>\phi</math> (frac) is the theoretical [[Tumbling Mill (Speed)|fraction critical speed]], and the fraction of critical speed at which centrifuging actually occurs, <math>\phi_{\rm c}</math> (frac), is:

:<math> \phi_{\rm c} =
    \begin{cases}
      \phi & \phi>0.35(3.364-J_{\rm T}\\
      0.35(3.364-J_{\rm T}& \phi \leq 0.35(3.364-J_{\rm T}\\
    \end{cases}
</math>

The position of the '''shoulder of the charge''', <math>\theta_{\rm S}</math> (rad), is: 

:<math>\theta_{\rm S} = \frac{\pi}{2} - \left ( \theta_{\rm T} - \frac{\pi}{2} \right ) \big[ \left ( 0.336 + 0.1041 \phi \right ) + \left (1.54 - 2.5673 \phi \right ) J_{\rm T} \big]</math>

The inner surface '''radius of the charge''', <math>r_{\rm i}</math> (m) is:

:<math>r_{\rm i} = R \sqrt{1 - \frac{2 \pi \beta J_{\rm T}}{2 \pi + \theta_{\rm S} - \theta_{\rm T}}}</math>

where the active load fraction, <math>\beta</math> (frac), is:

:<math>\beta = \frac{\left ( \dfrac{2 \pi + \theta_{\rm S} - \theta_{\rm T}}{\pi N_{\rm m}} \right )}{\sqrt{\dfrac{R}{g} \left ( 1 + \sqrt{1 - \dfrac{2 \pi J_{\rm T}}{2 \pi + \theta_{\rm S} - \theta_{\rm T}} } \right ) (\sin \theta_{\rm S} - \sin \theta_{\rm T})} + \left( \dfrac{2 \pi + \theta_{\rm S} - \theta_{\rm T}}{\pi N_{\rm m}} \right )}</math>

=== 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, <math>U</math>:
* if <math>U\leq1</math> then the charge void space is partially or exactly filled with slurry and <math>\theta_{\rm L} \geq \theta_{\rm T}</math>
* if <math>U>1</math> then there is an excess of slurry beyond the available charge void space, which manifests as a slurry pool, and <math>\theta_{\rm L} < \theta_{\rm T}</math>.

==== Slurry toe when U≤1 ====

When the charge void space is partially or exactly filled with slurry (<math>U\leq1</math>), the position of the '''toe of the slurry component''', <math>\theta_{\rm L}</math> (rad), is:

:<math>\theta_{\rm L} = \theta_{\rm T} + (2 \pi + \theta_{\rm S} - \theta_{\rm T}).(1 - U).y</math>

where <math>y</math> is the profile shape factor (frac).

The inner surface '''radius of the slurry component''', <math>r'_{\rm i}</math> (m), is:

:<math>r'_{\rm i} = \sqrt{R^2 - \frac{U(R^2-r_{\rm i}^2).(2\pi+\theta_{\rm S}-\theta_{\rm T})}{2\pi+\theta_{\rm S}-\theta_{\rm L}}}</math>

==== Slurry toe when U>1 ====

[[File:TumblingMillPower9.png|thumb|450px|Figure 3. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when <math>U>1</math>, along with key dimensions and areas.]]

[[File:TumblingMillPower10.png|thumb|450px|Figure 4. Tumbling mill profile showing Hilden and Powell's assumed charge and slurry component shapes when <math>U>1</math>, <math>h<r_{\rm i}</math> and <math>h_2<r_{\rm i}</math>, along with key dimensions and areas.]]

When the charge void space is overfilled (<math>U>1</math>), 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 (<math>A_{\rm Sect}</math>), and
* in the circular segment settled across the centre of the mill (<math>A_{\rm Seg}</math>).

Only the slurry in the annular sector (<math>A_{\rm Sect}</math>) 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. <math>r'_{\rm i} = r_{\rm i}</math>.

Locating the position of the '''toe of the slurry component''', <math>\theta_{\rm L}</math>, then becomes a problem of determining the distribution of slurry between the two areas, <math>A_{\rm Seg}</math> and <math>A_{\rm Sect}</math>.

The total area of excess slurry, <math>A_x</math> (m<sup>2</sup>), in the pool is:

:<math>A_x = A_{\rm Seg} + A_{\rm Sect} = (U - 1) . \frac{1}{2} (R^2 - r_{\rm i}^2).(2 \pi + \theta_{\rm S} - \theta_{\rm T})</math>

The area of the annular segment, <math>A_{\rm Seg}</math>, is:

:<math>A_{\rm Seg} = 
\begin{cases}
r_{\rm i}^2 . \arccos \left ( \frac{h}{r_{\rm i}} \right ) - h . \sqrt{r_{\rm i}^2 - h^2} & \text{if }h<r_{\rm i}\\
\left ( r_{\rm i}^2.\arccos \left ( \frac{h}{r_{\rm i}} \right ) - h.\sqrt{r_{\rm i}^2 - h^2} \right ) - \left ( r_{\rm i}^2.\arccos \left ( \frac{h_2}{r_{\rm i}} \right ) - h_2.\sqrt{r_{\rm i}^2 - h_2^2} \right ) & \text{if }h<r_{\rm i} \text{ and }h_2<r_{\rm i}
\end{cases}
</math>

where the parameters <math>h</math> and <math>h_2</math>, as shown in Figure 4, are:

:<math>h = - \left ( \frac{R + r_{\rm i}}{2} \right ) . \sin \theta_{\rm L}</math>

:<math>h_2 = - \left ( \frac{R + r_{\rm i}}{2} \right ) . \sin \theta_{\rm T}</math>

The area of the circular sector, <math>A_{\rm Sect}</math>, is:

:<math>A_{\rm Sect} = \frac{(\theta_{\rm T} - \theta_{\rm L}}{2} . (R^2 - r_{\rm i}^2)</math>

The slurry toe parameter <math>\theta_{\rm L}</math> cannot be algebraically isolated in the above equations. Therefore, an iterative numerical method is used to solve the value of <math>\theta_{\rm L}</math> that yields the required value of <math>A_x</math>.

=== Charge density ===

The masses of balls, <math>m_{\rm B}</math> (t), and coarse rocks, <math>m_{\rm R}</math> (t), in the charge are:

:<math>m_{\rm B} =J_{\rm B}.V_{\rm C}.(1 - \varepsilon). \rho_{\rm B}</math>

:<math>m_{\rm R} = \rho_{\rm R}.V_{\rm C}. \big (J_{\rm T}(1 - \varepsilon) - J_{\rm B}(1 - \varepsilon) \big )</math>

where:
* <math>J_{\rm B}</math> is the ball filling (v/v)
* <math>V_{\rm C}</math> is the volume of of the cylindrical section of the mill (m<sup>3</sup>)
* <math>\varepsilon</math> is the charge void fraction, typically 0.4 (v/v)
* <math>\rho_{\rm B}</math> is the Specific Gravity or density of the balls (- or t/m<sup>3</sup>)
* <math>\rho_{\rm R}</math> is the Specific Gravity or density of the rock in the charge (- or t/m<sup>3</sup>) 

The volume of of the cylindrical section of the mill (<math>V_{\rm C}</math>) is:

:<math>V_{\rm C} = \pi R^2 L</math>

The charge density, <math>\rho_{\rm c}</math> (t/m<sup>3</sup>), excluding slurry, is then:

:<math>\rho_{\rm c} = \frac{m_{\rm B} + m_{\rm R}}{J_{\rm T} V_{\rm C}}</math>

The apparent density of slurry within the charge, <math>\rho'_c</math> (t/m<sup>3</sup>), is:

:<math>\rho'_c = \varepsilon.\rho_{\rm s}</math>

where <math>\rho_{\rm s}</math> is the density of slurry (t/m<sup>3</sup>). The value of <math>\rho_{\rm s}</math> is computed as:

:<math>\rho_{\rm s} = S_{\rm V}.(\rho_{\rm R} - 1) + \rho_{\rm L}</math>

where <math>\rho_{\rm L}</math> is the density of liquids (t/m<sup>3</sup>) and <math>S_{\rm V}</math> is the volume fraction of solids in the slurry (v/v). <math>S_{\rm V}</math> is calculated from the mass fraction of solids in the slurry, <math>S</math> (w/w), as:

:<math>S_{\rm V} = \frac{S}{\rho_{\rm R} \left ( \dfrac{S}{\rho_{\rm R}} + \dfrac{1 - S}{\rho_{\rm L}} \right )}</math>

=== No-load power ===

The no-load power (<math>P_{\rm N}</math>) is estimated as:

:<math>P_{\rm N} = \text{Motor loss factor (kW/kW)} \times P_{\rm max}</math>

where <math>P_{\rm max}</math> 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:

:{| class="wikitable" style="width: 20%;"
|+ Table 1. Typical motor loss factors (after Hilden and Powell (2018)){{Hilden and Powell (2018)}}
|-
! rowspan="2" | Drive Type
! colspan="2" | 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:

<syntaxhighlight lang="vb">=mdMillPower_HildenPowell(Parameters as Range)</syntaxhighlight>

{{Excel (Text, Help, No Arguments)}}

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:

{| 
|- style="vertical-align:top;"
|
:<math>Parameters=
\begin{bmatrix}
D\text{ (m)}\\
L\text{ (m)}\\
d_{\rm T}\text{ (m)}\\
\alpha_{c}\text{ (degrees)}\\
J_{\rm T}\text{ (v/v)}\\
J_{\rm B}\text{ (v/v)}\\
\rho_{\rm R}\text{ (t/m}^{\text{3}}\text{)}\\
\rho_{\rm L}\text{ (t/m}^{\text{3}}\text{)}\\
\rho_{\rm B}\text{ (t/m}^{\text{3}}\text{)}\\
\varepsilon\text{ (v/v)}\\
U\text{ (v/v)}\\
S\text{ (}\% \text{ w/w)}\\
\phi\text{ (frac)}\\
\text{Motor loss factor (kW/kW)}\\
\text{Installed motor power (kW)}\\
k\text{ (kW/kW)}\\
y\text{ (frac)}\\
\end{bmatrix},\;\;\;\;\;\;

mdMillPower\_HildenPowell=
\begin{bmatrix}
P_{\rm N}\text{ (kW)}\\
P_{\rm Net}\text{ (kW)}\\
P\text{ (kW)}\\
\theta_{S} \text{ (rad)}\\
\theta_{\rm T} \text{ (rad)}\\
\theta_{L} \text{ (rad)}\\
r_{\rm i}\text{ (m)}\\
r'_{\rm i}\text{ (m)}\\
\end{bmatrix}\;\;\;\;\;\;
</math>


where:

:<math>P_{\rm Net} = k.(P_{\rm L} + P'_{\rm L} + P_{\rm C} + P'_{\rm C})</math>

| [[File:TumblingMillPower11.png|frame|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.

{{SysCAD (Text, Table Header)}}

|-
! colspan="3" style="text-align:left;" |''HildenPowell''

{{SysCAD (Text, Help Link)}}

|-
|MillDiameter
|style="background: #eaecf0" | Display
|Diameter of the mill (inside liners).
|-
|BellyLength
|style="background: #eaecf0" | Display
|Length of the cylindrical section (belly) of the mill (inside liners).
|-
|TrunnionDiameter
|style="background: #eaecf0" | Display
|Diameter of the trunnion (inside liners).
|-
|ConeAngle
|style="background: #eaecf0" | Display
|Angular displacement of the cone surface from the vertical direction.
|-
|Jt
|style="background: #eaecf0" | Display
|Volumetric fraction of the mill occupied by balls and coarse rock (including voids).
|-
|BallLoadVol
|style="background: #eaecf0" | Display
|Volumetric fraction of the mill occupied by balls (including voids).
|-
|SolidsSG
|style="background: #eaecf0" | Display
|Specific Gravity or density of solids.
|-
|BallSG
|style="background: #eaecf0" | Display
|Specific Gravity or density of balls.
|-
|LiquidsSG
|style="background: #eaecf0" | Display
|Specific Gravity or density of liquids.
|- style="vertical-align:top;"
|Voidage
|style="background: #eaecf0" | Input/Display
|Volumetric fraction of interstitial void space in the charge. Usually 0.4. <br> Note: The [[Ball Mill (Perfect Mixing, Dynamic)|Dynamic Perfect Mixing ball mill model]] uses [[Ball Mill (Overfilling)#Slurry hold-up below pool level|Shi's dynamic charge porosity estimation]] to align power draw predictions with the adopted slurry filling approach.
|-
|U
|style="background: #eaecf0" | Display
|Volumetric fraction of interstitial grinding media voidage occupied by slurry.
|-
|Cw
|style="background: #eaecf0" | Display
|Mass fraction of solids in discharge slurry.
|-
|FracCS
|style="background: #eaecf0" | 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
|style="background: #eaecf0" | Display
|Angle of the charge shoulder.
|-
|ThetaToe
|style="background: #eaecf0" | Display
|Angle of the charge toe.
|-
|ThetaSlurryToe
|style="background: #eaecf0" | Display
|Angle of the toes of the slurry component
|-
|ChargeSurfaceRadius
|style="background: #eaecf0" | Display
|Radius of the inner surface of the charge.
|-
|SlurrySurfaceRadius
|style="background: #eaecf0" | Display
|Radius of the inner surface of the slurry component.
|-
|NoLoadPower
|style="background: #eaecf0" | Display
|Power input to the motor when the mill is empty (no balls, rocks or slurry).
|-
|NetPower
|style="background: #eaecf0" | Display
|Charge motion power, including losses.
|-
|GrossPower
|style="background: #eaecf0" | Display
|Power input to the motor.
|}

== See also ==

* [[Tumbling Mill (Power, Morrell Continuum)|Morrell Continuum model]]

== References ==

[[Category:Excel]]
[[Category:SysCAD]]