Introduction & Context

The calculation presented here estimates the steady-state surface temperature of a cylindrical mill (or similar rotating equipment) that generates heat internally due to power dissipation. In process engineering, excessive surface temperatures can lead to material degradation, safety hazards, and reduced equipment lifespan. Engineers use this analysis to verify that the mill’s operating conditions keep the surface temperature below a prescribed allowable limit, thereby ensuring reliable operation and compliance with safety standards.

This methodology is typically applied to:

  • Grinding or mixing mills where motor losses are converted to heat.
  • Rotary kilns and reactors with internal heat sources.
  • Any cylindrical component where lumped-capacitance assumptions are reasonable and heat is removed by convection and radiation to the surrounding air.

Methodology & Formulas

The analysis proceeds through a series of physics-based steps, each expressed algebraically.

1. Geometry and Material Properties

Define the following symbols:

  • \(L\) – Length of the cylinder (m)
  • \(r_1\) – Inner radius where heat is generated (m)
  • \(r_2\) – Outer radius exposed to the environment (m)
  • \(k\) – Thermal conductivity of the cylinder material (W·m\(^{-1}\)·K\(^{-1}\))
  • \(h\) – Convective heat-transfer coefficient (W·m\(^{-2}\)·K\(^{-1}\))
  • \(\varepsilon\) – Surface emissivity (dimensionless)
  • \(\sigma\) – Stefan-Boltzmann constant (W·m\(^{-2}\)·K\(^{-4}\))
  • \(T_{\text{air}}\) – Ambient temperature (°C)
  • \(T_{s,\max}\) – Allowable maximum surface temperature (°C)
  • \(P\) – Electrical power dissipated as heat (W)

2. Internal Heat Generation

The volume in which heat is generated is the inner cylinder:

\[ V_{\text{gen}} = \pi r_1^{2} L \]

The volumetric heat-generation rate is then:

\[ \dot{q}_{\text{gen}} = \frac{P}{V_{\text{gen}}} \]

3. External Surface Area

Heat leaves the outer surface (radius \(r_2\)) by convection and radiation. The relevant area is the lateral surface:

\[ A = 2\pi r_2 L \]

4. Energy Balance at the Surface

At steady state, the total heat generated equals the sum of convective and radiative heat losses:

\[ \dot{q}_{\text{gen}} V_{\text{gen}} = h A \bigl(T_s - T_{\text{air}}\bigr) + \varepsilon \sigma A \bigl[(T_s+273.15)^{4} - (T_{\text{air}}+273.15)^{4}\bigr] \]

Rearranging yields a root-finding function \(f(T_s)\):

\[ f(T_s) = h A \bigl(T_s - T_{\text{air}}\bigr) + \varepsilon \sigma A \bigl[(T_s+273.15)^{4} - (T_{\text{air}}+273.15)^{4}\bigr] - \dot{q}_{\text{gen}} V_{\text{gen}} \]

5. Newton–Raphson Solution

The temperature is obtained iteratively using:

\[ T_s^{(n+1)} = T_s^{(n)} - \frac{f\!\bigl(T_s^{(n)}\bigr)}{f'\!\bigl(T_s^{(n)}\bigr)} \]

where the derivative is:

\[ f'(T_s) = h A + 4 \varepsilon \sigma A \bigl(T_s+273.15\bigr)^{3} \]

6. Validity Checks (Lumped-Capacitance Assumption)

The lumped-capacitance model is acceptable when the Biot number is small:

\[ \text{Bi} = \frac{h\,(r_2 - r_1)}{k} \]

Additional engineering checks include the range of the convective coefficient and the admissible emissivity values.

CriterionExpressionAcceptable Range
Biot Number \(\displaystyle \text{Bi} = \frac{h\,(r_2 - r_1)}{k}\) \(\text{Bi} < 0.1\)
Convective Coefficient \(h\) \(5 \le h \le 25\) W·m\(^{-2}\)·K\(^{-1}\)
Surface Emissivity \(\varepsilon\) \(0.05 \le \varepsilon \le 0.95\)

7. Decision Logic

After solving for \(T_s\), compare it with the allowable limit \(T_{s,\max}\):

  • If \(T_s > T_{s,\max}\) → ACTION REQUIRED: redesign cooling, reduce power, or modify geometry.
  • If \(T_s \le T_{s,\max}\) → Ts within allowable limit: operation is acceptable under the assumed conditions.