View source for Ball Mill (Overfilling)

== Description ==

This article describes several methods for estimating the maximum volumetric flow capacity of an overflow ball mill, including the '''Shi''' (2016) and '''Arbiter''' (1990) approaches.{{Shi (2016)}}{{Arbiter (1990)}}

== Model theory ==

=== Shi (2016) ===

[[File:BallMillOverfilling1.png|thumb|450px|Figure 1. Overflow discharge tumbling mill profile showing Shi's assumed charge and slurry pool areas.]]

Shi estimated the residence time of an overflow ball mill by considering the volume of slurry present in the grinding media interstices and slurry pool.

This was accomplished by adopting Morrell's simplified charge geometry and calculating the volume of slurry resident in each of the areas A-C in Figure 1:{{Morrell (1996a)}}

:A. In the grinding media interstices above the slurry pool level
:B. In the grinding media interstices below the slurry pool level
:C. In the slurry pool

The total volume of slurry in the mill is therefore the sum of slurry in areas A-C.

Estimating the slurry volume in areas A-C requires definition of the simplified charge geometry and slurry pool level for the mill in question.

==== Charge position ====

The simplified charge geometry is defined by three parameters:
* the angular position of the ''toe'' of the charge, <math>\theta_t</math> (rad)
* the angular position of the ''shoulder'' charge, <math>\theta_s</math> (rad)
* the radius of the ''inner charge surface'', <math>R_i</math> (m)

Morrell's (1996) equations are used to define the charge position, i.e.{{Morrell (1996a)}}

The position of the '''toe''' of the charge is defined as: 

:<math>\theta_t = 2.5307 \left ( 1.2796 - J_t \right ) \left (1 - {\rm e}^{-19.42(\phi_c - \phi)} \right ) + \frac{\pi}{2}</math>

where:
* <math>J_{t}</math> is the volumetric fraction of the mill occupied by balls and coarse rock (including void space and interstitial slurry) (v/v)
* <math>\phi</math> (frac) is the theoretical [[Tumbling Mill (Speed)|fraction critical speed]]

The fraction of critical speed at which centrifuging actually occurs, <math>\phi_c</math> (frac), is:

:<math> \phi_c =
    \begin{cases}
      \phi & \phi>0.35(3.364-J_t)\\
      0.35(3.364-J_t)& \phi \leq 0.35(3.364-J_t)\\
    \end{cases}
</math>

The position of the '''shoulder''' of the charge is: 

:<math>\theta_s = \frac{\pi}{2} - \left ( \theta_t - \frac{\pi}{2} \right ) \big[ \left ( 0.336 + 0.1041 \phi \right ) + \left (1.54 - 2.5673 \phi \right ) J_t \big]</math>

The inner '''charge surface radius''' is:

:<math>R_i = r_m \left ( 1 - \frac{2 \pi \beta J_t}{2 \pi + \theta_s - \theta_t} \right )^{0.5}</math>

where the fraction of total charge in the active region, <math>\beta</math> (frac), is:

:<math>\beta = \frac{t_c}{t_f + t_c}</math>

The time taken to travel between the toe and shoulder of the charge during one revolution, <math>t_c</math> (s), is:

:<math>t_c = \frac{2 \pi - \theta_t + \theta_s}{2 \pi \bar N}</math>

where the mean rotational rate, <math>\bar N</math> (rev/s), is:

:<math>\bar N \approx \frac{N_m}{2}</math>

The time taken to travel between the shoulder and toe of the charge in free flight during one revolution, <math>t_f</math> (s), is:

:<math>t_f \approx \left ( \frac{2 \bar r \left (\sin \theta_s - \sin \theta_t \right )}{g} \right )^{0.5}</math>

where the mean radial position <math>\bar r</math> (m), is:
:<math>\bar r = \frac{r_m}{2} \left [ 1 + \left ( 1 - \frac{2 \pi J_t}{2 \pi + \theta_s - \theta_t} \right )^{0.5} \right ]</math>

==== Slurry hold-up below pool level ====

The volume of slurry held up below the pool level in areas B and C, <math>V_1</math> (m<sup>3</sup>), is:

:<math>V_1 = L_m (S_{pool} - S_{net\,ball})</math>

where:
* <math>L_m</math> is effective mill length (m)
* <math>S_{pool}</math> is cross-sectional area of the slurry pool (m<sup>2</sup>)
* <math>S_{net\,ball}</math> is cross-sectional area occupied by the grinding media, excluding balls, below the slurry pool level(m<sup>2</sup>)

This implementation of the Shi model computes <math>L_m</math> as an ''effective'' mill length from:

:<math>L_m = \dfrac{V_m}{\pi {R_m}^2}</math>

where <math>V_m</math> is the volume of the mill (m<sup>3</sup>) and <math>R_m</math> is the radius of the mill (m).

The volume of the mill is calculated as the sum of a cylinder and two right  circular frustums:{{Gupta and Yan (2016)}}

:<math>V_m = \pi {R_m}^2L + 2 \cdot \bigg[ \dfrac{\pi}{3} (R_m - R_t) \cdot \tan \alpha_c \cdot \left ( {R_{m}}^{2} + R_{m} R_{t} + {R_{t}}^{2} \right) \bigg]</math>

wher:
* <math>L</math> is the length of the cylindrical (belly) section of the mill (m)
* <math>R_t</math> is the radius of the discharge trunnion (m)
* <math>\alpha_c</math> is the cone angle, measured as the angular displacement of the cone surface from the vertical direction (rad)

<math>S_{pool}</math> is computed from:

:<math>S_{pool} = 0.5 {R_m}^2 (\theta_p - \sin \theta_p)</math>

:<math>\theta_p = \pi - 2 \alpha</math>

:<math>\alpha = \arcsin \left (\dfrac{R_t - \Delta h}{R_m} \right )</math>

where <math>R_t</math> is the radius of the discharge trunnion (m), and <math>\Delta h</math> is the height of the slurry pool level over the trunnion lip (m).

The value of <math>\Delta h</math> may be estimated  as <math>0.67 h</math>, where <math>h</math> is the slurry height above the trunnion level inside the mill (m):{{Morrell (2016)}}

:<math>h = k_{\mu} \cdot \frac{1}{0.67} \left ( \frac{2Q}{\pi R_t g^{0.5}} \right )^{\frac{2}{3}}</math>

and:
* <math>k_{\mu}</math> is an empirical coefficient related to the ''viscidity'' of the slurry (Morrell suggests <math>k_{\mu} \approx 2</math>, Shi applies <math>k_{\mu} = 1</math>)
* <math>Q</math> is the volumetric discharge rate of pulp from the mill (m<sup>3</sup>/s)
* <math>g</math> is acceleration due to gravity (m/s<sup>2</sup>)

<math>S_{net\,ball}</math> is computed as:

:<math>S_{net\,ball} = 0.5 \theta_b (1 - \epsilon_d)({R_m}^2 - {R_t}^2)</math>

:<math>\theta_b = 2 \pi - \alpha - \theta_t</math>

where the ''dynamic charge porosity'', <math>\epsilon_d</math> (v/v), is:

:<math>\epsilon_d = 0.4 + 0.228 \left [ 1 - \exp \left (-0.315\dfrac{\phi}{J_t} \right) \right ]</math>

==== Slurry hold up above pool level ====

The volume of slurry held up in the grinding media interstices above the slurry pool level in area A, <math>V_2</math> (m<sup>3</sup>), is:

:<math>V_2 = L_m . \epsilon_d . S_{ball\,above\,pool}</math>

where:

:<math>S_{ball\,above\,pool} = 0.5 (\theta_s + \alpha)({R_m}^2 - {R_t}^2)</math>

==== Residence time and overfilling ====

The total volume of slurry hold up in the mill is <math>V_{total}=V_1+V_2</math> (m<sup>3</sup>).

The residence time of slurry in the mill, <math>t_{Res}</math> (s), is then:

:<math>t_{Res} = \dfrac{V_{total}}{Q}</math>

The ''axial velocity'' of slurry through the mill, <math>v_{ax}</math> (m/s), is:

:<math>v_{ax} = \dfrac{L_m}{t_{Res}}</math>

Shi suggests the following residence time limits for overflow discharge ball mills (s):

:<math>
t_{Limit} = 
\begin{cases}
120\text{ s} & D_m<5.9\text{ m}\\
60\text{ s} & D_m\geq5.9\text{ m}\\
\end{cases}
</math>

where <math>D_m</math> is the diameter of the mill (m), i.e. <math>D_m = 2 R_m</math>.

The maximum volumetric flow rate of the mill, <math>Q_{Max}</math> (m<sup>3</sup>/s), at the residence time limit <math>t_{Limit}</math> may be back-calculated using the above relations. However, as the height of slurry above the trunnion lip (<math>h</math>), and hence residence time (<math>t_{Res}</math>), is a function of flow rate, no analytical solution is available and a numerical method is required for computation.

=== Arbiter (1990) ===

Arbiter describes two approaches for identifying the overfilling of an overflow discharge ball mill:
#. the ''critical flow number'' method
#. the ''axial flow velocity'' method

Arbiter's limits are rearranged to yield the following relations for the maximum volumetric discharge rate of slurry from a mill:

:<math>Q_{Max}^{N_q} = 0.0217 \dfrac{\pi^2}{4} {D_m}^3 N_{eq} (0.5 - 0.66 J_t) \cdot 60</math>

:<math>Q_{Max}^{v_{ax}} = 0.0769 \dfrac{\pi}{4} {D_m}^{2} (0.5 - 0.66 J_t) \cdot 3600</math>

where:
* <math>Q_{Max}^{N_q}</math> is the critical flow number volumetric slurry discharge rate limit (m<sup>3</sup>/h)
* <math>Q_{Max}^{v_{ax}}</math> is the axial velocity volumetric slurry discharge rate limit (m/s)
* <math>D_m</math> is mill diameter (m)
* <math>N_{eq}</math> is mill rotational speed (rev/min)
* <math>J_t</math> is the volumetric charge fraction in the mill (v/v)

Overfilling may be a risk if the actual volumetric flow rate from the mill approaches or exceeds either of <math>Q_{Max}^{N_q}</math> or <math>Q_{Max}^{v_{ax}}</math>.

== Additional notes ==

Shi estimated the residence times for a database of 121 overflow ball mills. These residence times are presented as cumulative frequency distributions in Figure 2, allowing the performance a given mill to be ranked against the database.

:{|
|-
|style="padding: 10px" | [[File:BallMillOverfilling2.png|thumb|none|650px|Figure 2. Cumulative frequency distribution of the volume-based residence time of 121 mills in Shi's database (after Shi, 2016).{{Shi (2016)}}]]
|}

Arbiter's relations were developed from a limited database of five mills, the largest of which was 5.3 m diameter x 6.4 m length.

== Excel ==

=== Shi ===

The Shi overflow discharge ball mill overfilling model may be invoked from the Excel formula bar with the following function call:

<syntaxhighlight lang="vb">=mdMillOverfilling_Shi(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_m\text{ (m)}\\
L\text{ (m)}\\
R_t\text{ (m)}\\
\alpha_{c}\text{ (degrees)}\\
\phi\text{ (frac)}\\
J_t\text{ (v/v)}\\
Q\text{ (m}^{3}\text{/h)}\\
k_{\mu}\\
\end{bmatrix},\;\;\;\;\;\;

mdMillOverfilling\_Shi=
\begin{bmatrix}
v_{ax}\text{ (m/s)}\\
t_{Res}\text{ (s)}\\
t_{Limit}\text{ (s)}\\
Q_{Max}\text{ (m}^{3}\text{/h)}\\
V_{1}\text{ (m}^{3}\text{)}\\
V_{2}\text{ (m}^{3}\text{)}\\
V_{total}\text{ (m}^{3}\text{)}\\
\end{bmatrix}\;\;\;\;\;\;
</math>

| [[File:BallMillOverfilling3.png|frame|Figure 3. Example showing the selection of the '''Parameters''' (blue frame), and '''Results''' (light blue frame) arrays in Excel.]]
|}

=== Arbiter ===

The Arbiter overflow discharge ball mill overfilling model may be invoked from the Excel formula bar with the following function call:

<syntaxhighlight lang="vb">=mdMillOverfilling_Arbiter(Parameters as Range)(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_m\text{ (m)}\\
\phi\text{ (frac)}\\
J_t\text{ (v/v)}\\
\end{bmatrix},\;\;\;\;\;\;

mdMillOverfilling\_Arbiter=
\begin{bmatrix}
Q_{Max}^{N_q}\text{ (m}^{3}\text{/h)}\\
Q_{Max}^{v_{ax}}\text{ (m}^{3}\text{/h)}\\
\end{bmatrix}\;\;\;\;\;\;
</math>

| [[File:BallMillOverfilling4.png|frame|Figure 4. Example showing the selection of the '''Parameters''' (blue frame), and '''Results''' (light blue frame) arrays in Excel.]]
|}

== SysCAD ==

The Ball Mill Overfilling Indicator model is an optional calculation for tumbling mill units. If selected, the input and display parameters below are shown.

{{SysCAD_Table_Header}}
|-
! colspan="3" style="text-align:left;" |''Shi''
|-
|MillDiameter
|style="background: #eaecf0" | Input/Display
|Diameter of the mill (inside liners).
|-
|BellyLength
|style="background: #eaecf0" | Input/Display
|Length of the cylindrical section (belly) of the mill (inside liners).
|-
|TrunnionDiameter
|style="background: #eaecf0" | Input/Display
|Diameter of the trunnion (inside liners).
|-
|ConeAngle
|style="background: #eaecf0" | Input/Display
|Angular displacement of the cone surface from the vertical direction.
|-
|FracCS
|style="background: #eaecf0" | Input/Display
|Fraction critical speed of the mill.
|-
|Jt
|style="background: #eaecf0" | Input/Display
|Volumetric fraction of the mill occupied by balls and coarse rock (including voids).
|-
|ViscidityCoeff
| Input
| Coefficient of viscidity.
|-
|V1
|style="background: #eaecf0" | Display
|Volume of slurry below slurry pool level.
|-
|V2
|style="background: #eaecf0" | Display
|Volume of slurry above slurry pool level.
|-
|V1
|style="background: #eaecf0" | Display
|Total volume of slurry in charge and pool.
|-
|AxialVelocity
|style="background: #eaecf0" | Display
|Axial velocity of slurry flow through the charge and pool.
|-
|tRes
|style="background: #eaecf0" | Display
|Residence time of slurry in mill.
|-
|tLimit
|style="background: #eaecf0" | Display
|Shi's suggested lower limit of residence time for mill. Based on mill diameter.
|-
|Feed.SLQv
|style="background: #eaecf0" | Display
|Volumetric flow rate of slurry (solids + liquids) in mill feed.
|-
|Feed.SLQv.Limit
|style="background: #eaecf0" | Display
|Estimated volumetric flow rate of slurry (solids + liquids) in mill feed at tLimit.
|-
|Overfilled
|style="background: #eaecf0" | True/False
|Indicates if tRes is less than tLimit, i.e. mill is overfilled.
|-
! colspan="3" style="text-align:left;" |''Arbiter''
|-
|MillDiameter
|style="background: #eaecf0" | Input/Display
|Diameter of the mill (inside liners).
|-
|FracCS
|style="background: #eaecf0" | Input/Display
|Fraction critical speed of the mill.
|-
|Jt
|style="background: #eaecf0" | Input/Display
|Volumetric fraction of the mill occupied by balls and coarse rock (including voids).
|-
|Feed.SLQv
|style="background: #eaecf0" | Display
|Volumetric flow rate of slurry (solids + liquids) in mill feed.
|-
|Feed.SLQv.NqLimit
|style="background: #eaecf0" | Display
|Estimated volumetric flow rate of slurry (solids + liquids) at Arbiters critical flow number limit.
|-
|Feed.SLQv.VmaxLimit
|style="background: #eaecf0" | Display
|Estimated volumetric flow rate of slurry (solids + liquids) at Arbiters axial velocity limit.
|}

== See also ==

* [[Tumbling Mill (Power, Morrell Continuum)]]
* [[Ball Mill (Perfect Mixing)|Perfect Mixing ball mill model]]
* [[Mill (Herbst-Fuerstenau)]]

== References ==

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