Introduction & Context

The calculation presented estimates the mechanical power required by a roller mill that processes a visco-elastic material (e.g., chocolate, polymer melt, or food slurry). In process engineering, the power consumption of a refining or grinding mill is a key design parameter because it directly influences motor sizing, energy cost, heat generation, and product quality. The model is applicable to continuous roll-refining operations where two cylindrical rolls compress a thin film of material, creating a nip zone in which shear and normal stresses develop.

Methodology & Formulas

The procedure follows the sequence of the Python script, expressed here with algebraic symbols.

1. Geometry and Kinematic Conversions

Convert practical dimensions to SI units using the appropriate conversion factors:

\[ D = D_{mm}\,\alpha_{mm\to m},\qquad R = \frac{D}{2},\qquad h_{in} = h_{in,mm}\,\alpha_{mm\to m},\qquad h_{out} = h_{out,mm}\,\alpha_{mm\to m} \]

Roll surface speed:

\[ v = \frac{\pi D N}{60} \]

where \(N\) is the rotational speed in revolutions per minute.

2. Shear Rate in the Gap

The effective gap is taken as the outlet film thickness \(h_{out}\). The shear rate is therefore

\[ \dot{\gamma} = \frac{v}{\max(h_{out},\varepsilon)} \]

with \(\varepsilon\) representing a very small positive number to avoid division by zero.

3. Apparent Viscosity (Power-Law Model)

The material is described by a power-law (Ostwald-de Waele) relationship:

\[ \mu_{app} = \frac{\tau_{0}}{\max(\dot{\gamma},\varepsilon)} + K\,\dot{\gamma}^{\,n-1} \]

Conversion to centipoise (cP) uses the factor \(\beta_{Pa\to cP}\):

\[ \mu_{app}^{cP} = \frac{\mu_{app}}{\beta_{Pa\to cP}} \]

4. Nip Contact Length

The nip length is approximated by the square-root of the product of roll radius and the film-thickness reduction:

\[ L_{nip} = \sqrt{R\,\Delta h},\qquad \Delta h = h_{in} - h_{out} \]

5. Rolling Force per Unit Width

The normal force per unit width of roll is the sum of the yield component and the viscous component integrated over the nip length:

\[ F' = \tau_{0}\,L_{nip} + K\left(\frac{v}{\max(h_{out},\varepsilon)}\right)^{n} L_{nip} \]

6. Power Consumption

The power per unit width is the product of rolling force and surface speed:

\[ P' = F'\,v \]

The total power for a mill of effective width \(W\) is

\[ P = P'\,W,\qquad P_{kW} = \frac{P}{10^{3}} \]

7. Validity Checks

The empirical model is reliable only when the operating conditions fall within established ranges. These ranges are summarized in the table below.

Parameter	Acceptable Symbolic Range	Comment
Power-law exponent \(n\)	\(n_{min} \le n \le n_{max}\)	Typical for shear-thinning foods
Yield stress \(\tau_{0}\) (Pa)	\(\tau_{0,\,min} \le \tau_{0} \le \tau_{0,\,max}\)	Ensures valid yield-stress regime
Consistency index \(K\) (Pa·s\(^n\))	\(K_{min} \le K \le K_{max}\)	Empirical calibration limits
Thickness ratio \(h_{in}/h_{out}\)	\(R_{ratio,\,min} \le \frac{h_{in}}{h_{out}} \le R_{ratio,\,max}\)	Ensures sufficient compression
Roll surface speed \(v\) (m·s\(^{-1}\))	\(v_{min} \le v \le v_{max}\)	Maintains laminar shear regime
Shear rate \(\dot{\gamma}\) (s\(^{-1}\))	\(\dot{\gamma}_{min} \le \dot{\gamma} \le \dot{\gamma}_{max}\)	Within rheometer-validated region
Apparent viscosity \(\mu_{app}^{cP}\)	\(\mu_{min}^{cP} \le \mu_{app}^{cP} \le \mu_{max}^{cP}\)	Typical for chocolate-refining
Nip-length ratio \(L_{nip}/R\)	\(\frac{L_{nip}}{R} < \lambda_{max}\)	Ensures validity of nip-theory approximation

If any of the above conditions are violated, the user should treat the resulting power estimate with caution and consider recalibrating the material model or adjusting operating parameters.

Throughput rate and material grindability (Bond Work Index).
Roll surface friction and profile condition.
Target product fineness (specific surface area or % passing 45 µm).
Hydraulic grinding pressure and roll gap setting.
Internal recirculation load and classifier speed.
Moisture content and feed size distribution.

Calculate the specific power (kWh t⁻¹) and compare it with OEM guarantees or historical best achieved.
Check the power ratio: P_actual / P_rated; sustained values > 1.05 indicate overload.
Monitor the power trend versus feed rate; a steepening slope shows diminishing efficiency.
Compare with peer plants using the same material and product fineness; ±10 % is typical scatter.

Main motor current and instantaneous kW from the VFD.
Differential pressure across the grinding zone (surrogate for material bed thickness).
Vibration RMS at mill gearbox and classifier bearings (prevents hidden overload).
Recirculation bucket-elevator kW (reveals internal circulation build-up).
Near-infrared or microwave moisture at feed belt (moisture spikes raise power).

Re-groove or re-texture rolls when power rises > 5 % for the same throughput.
Replace worn edge seals to stop false air; every 1 % false air can add 0.3 kWh t⁻¹.
Calibrate hydraulic accumulators to the correct pre-charge pressure; under-charged systems force higher hydraulic pressure and power.
Clean classifier blades and check tip speed; a 2 % drop in classifier efficiency can recycle extra material and raise mill kWh t⁻¹ by 3–4 %.

Worked Example – Power Consumption of a Roller Mill

Scenario: A food-processing plant uses a single-roller mill to grind a viscous paste. The roller has a diameter of 400 mm, a width of 250 mm, and rotates at 15 rpm. The material enters the nip with a thickness of 2 mm and exits with a thickness of 0.2 mm, giving a gap of 0.2 mm. The paste behaves as a Herschel-Bulkley fluid with a yield stress \(\tau_0 = 250\;{\rm Pa}\), consistency index \(K = 800\;{\rm Pa\;s^{n}}\), and flow index \(n = 0.45\). Determine the mechanical power required to operate the mill.

Knowns

Roller diameter, \(D = 0.400\;{\rm m}\)
Roller radius, \(R = 0.200\;{\rm m}\)
Roller width, \(W = 0.250\;{\rm m}\)
Rotational speed, \(N = 15\;{\rm rpm}\)
Conversion factor, \(1\;{\rm rpm} = 0.10472\;{\rm rad\;s^{-1}}\)
Inlet thickness, \(h_{\text{in}} = 2.0\;{\rm mm} = 0.002\;{\rm m}\)
Outlet thickness, \(h_{\text{out}} = 0.2\;{\rm mm} = 0.0002\;{\rm m}\)
Gap height, \(\Delta h = h_{\text{in}} - h_{\text{out}} = 0.0018\;{\rm m}\)
Yield stress, \(\tau_0 = 250\;{\rm Pa}\)
Consistency index, \(K = 800\;{\rm Pa\;s^{n}}\)
Flow index, \(n = 0.45\)

Step-by-Step Calculation

Convert rotational speed to angular velocity:
\[ \omega = N \times 0.10472 = 15 \times 0.10472 = 1.571\;{\rm rad\;s^{-1}} \]
Calculate the linear surface speed of the roller:
\[ v = \omega R = 1.571 \times 0.200 = 0.314\;{\rm m\;s^{-1}} \]
Determine the shear rate in the nip (gap height \(\Delta h = 0.0002\;{\rm m}\)):
\[ \gamma = \frac{v}{\Delta h} = \frac{0.314}{0.0002} = 1570.796\;{\rm s^{-1}} \]
Evaluate the Herschel-Bulkley shear stress:
\[ \tau = \tau_0 + K\,\gamma^{\,n} \] First compute \(\gamma^{\,n}\): \[ \gamma^{\,n} = 1570.796^{0.45} = 27.440 \] Then \[ \tau = 250 + 800 \times 27.440 = 250 + 21955.200 = 22205.200\;{\rm Pa} \]
Apparent viscosity (optional check):
\[ \mu_{\text{app}} = \frac{\tau}{\gamma} = \frac{22205.200}{1570.796} = 14.130\;{\rm Pa\;s} \]
Compute the volume of material in the nip. The nip length is approximated by
\[ L_{\text{nip}} = \sqrt{R\,\Delta h} = \sqrt{0.200 \times 0.0018} = 0.018974\;{\rm m} \] Volume: \[ V = L_{\text{nip}} \times W \times \Delta h = 0.018974 \times 0.250 \times 0.0002 = 9.487 \times 10^{-7}\;{\rm m^{3}} \]
Power density (shear stress × shear rate):
\[ \dot{q} = \tau \,\gamma = 22205.200 \times 1570.796 = 3.488 \times 10^{7}\;{\rm W\;m^{-3}} \]
Total mechanical power required:
\[ P = \dot{q}\,V = 3.488 \times 10^{7} \times 9.487 \times 10^{-7} = 33.075\;{\rm W} \] In kilowatts: \[ P_{\text{kW}} = \frac{33.075}{1000} = 0.033\;{\rm kW} \]

Final Answer

The roller mill consumes 33.075 W (≈ 0.033 kW) of mechanical power under the stated operating conditions.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle