8.1 Steady-State Conduction in Planar, Cylindrical, and Spherical Geometries

8.1 Steady-State Conduction in Planar, Cylindrical, and Spherical Geometries

Conduction is the mode of heat transfer in which energy transfers from more energetic particles of a substance to adjacent less energetic particles as a result of microscopic interactions (molecular collisions, lattice vibrations, and free electron transport). In chemical processing plants—spanning catalytic fixed-bed reactors, cryogenic liquid nitrogen vessels, refractory-lined blast furnaces, and multi-mile steam distribution headers—steady-state conduction governs thermal barrier sizing, structural safety, energy conservation, and condensation prevention.

On the NCEES PE Chemical Exam, conduction problems test your mastery of Fourier's law, the thermal resistance network analogy (Ohm's law analogy), coordinate system transformations (planar, cylindrical, spherical), and the subtle, counterintuitive physics of the critical radius of insulation.

---

1. Fundamentals of Thermal Conduction & Fourier's Law

Thermal conduction is governed phenomenologically by Fourier's Law of Heat Conduction, which states that the local heat flux vector $\mathbf{q}''$ is directly proportional to the negative gradient of the temperature field:

$$\mathbf{q}'' = -k \nabla T$$

Where:

• $\mathbf{q}''$ = heat flux vector ($\text{W/m}^2$ or $\text{Btu/(hr}\cdot\text{ft}^2)$).
• $k$ = thermal conductivity of the medium ($\text{W/(m}\cdot\text{K)}$ or $\text{Btu/(hr}\cdot\text{ft}\cdot^\circ\text{F)}$).
• $\nabla T$ = spatial temperature gradient vector ($\text{K/m}$ or $^\circ\text{F/ft}$).
• The negative sign enforces the Second Law of Thermodynamics: heat flows spontaneously in the direction of decreasing temperature.

For one-dimensional, steady-state conduction without internal heat generation in the Cartesian $x$-direction:

$$q_x = -k A \frac{dT}{dx}$$

Where $q_x$ is the total heat transfer rate ($\text{W}$ or $\text{Btu/hr}$) and $A$ is the cross-sectional area perpendicular to the direction of heat flow ($\text{m}^2$ or $\text{ft}^2$).

Temperature-Dependent Thermal Conductivity

In most process engineering materials, thermal conductivity varies with temperature. Over moderate operating ranges, this dependence is modeled linearly:

$$k(T) = k_0 (1 + \beta T)$$

Integrating Fourier's law over a slab of thickness $L$ with boundary temperatures $T_1$ and $T_2$ yields:

$$q = \frac{A}{L} \int_{T_2}^{T_1} k(T) dT = \frac{k_m A (T_1 - T_2)}{L}$$

Where $k_m$ is evaluated precisely at the arithmetic mean temperature of the slab:

$$k_m = k(T_m) = k_0 \left( 1 + \beta \frac{T_1 + T_2}{2} \right)$$

---

2. The Thermal Resistance Analogy (Ohm's Law Analogy)

Heat transfer calculations for complex composite equipment are greatly simplified by treating heat flow analogously to electrical current in direct-current (DC) electric circuits:

$$\text{Thermal Current } (q) = \frac{\text{Thermal Driving Force } (\Delta T)}{\text{Thermal Resistance } (R_{th})}$$

Where:

• $q$ corresponds to electric current ($I$).
• $\Delta T = T_{hot} - T_{cold}$ corresponds to electrical potential difference ($\Delta V$).
• $R_{th}$ is the thermal resistance ($\text{K/W}$ or $^\circ\text{C/W}$ or $\text{hr}\cdot^\circ\text{F/Btu}$).

Boundary Convection Resistance

When a solid boundary is exposed to a moving fluid stream at bulk temperature $T_\infty$, heat transfers from the surface at $T_s$ via convection, governed by Newton's Law of Cooling:

$$q = h A (T_s - T_\infty) \implies R_{conv} = \frac{1}{h A}$$

Where $h$ is the convective heat transfer coefficient ($\text{W/(m}^2\cdot\text{K)}$ or $\text{Btu/(hr}\cdot\text{ft}^2\cdot^\circ\text{F)}$).

Overall Heat Transfer Coefficient ($U$)

For a composite system with total thermal resistance $R_{total}$:

$$q = U A \Delta T_{overall} = \frac{\Delta T_{overall}}{R_{total}} \implies U A = \frac{1}{R_{total}}$$

Where $U$ is the overall heat transfer coefficient referenced to a designated area $A$.

---

3. Planar Conduction & Composite Multi-Layer Walls

For a flat, homogeneous solid wall of thickness $L$, cross-sectional area $A$, and constant thermal conductivity $k$:

$$q = -k A \frac{dT}{dx} \implies q \int_0^L dx = -k A \int_{T_1}^{T_2} dT \implies q = \frac{k A (T_1 - T_2)}{L}$$

The conductive thermal resistance of a flat plate is:

$$R_{cond,wall} = \frac{L}{k A}$$

The specific thermal resistance per unit area (area-normalized resistance) is:

$$R''_{cond} = R_{cond} \cdot A = \frac{L}{k} \quad \left[\text{m}^2\cdot\text{K/W}\right]$$

Series Composite Planar Wall

Consider a multi-layer refractory wall (e.g., an industrial reformer or furnace) consisting of $n$ distinct material layers in series, exposed to a hot process gas ($T_{\infty,1}, h_1$) on the inside and ambient air ($T_{\infty,2}, h_2$) on the outside.

At steady state, the rate of heat transfer $q$ through each layer is identical:

$$q = h_1 A (T_{\infty,1} - T_1) = \frac{k_1 A}{L_1}(T_1 - T_2) = \frac{k_2 A}{L_2}(T_2 - T_3) = \dots = h_2 A (T_n - T_{\infty,2})$$

Summing all temperature drops yields the total thermal resistance:

$$R_{total} = \frac{1}{h_1 A} + \sum_{i=1}^n \frac{L_i}{k_i A} + \frac{1}{h_2 A}$$ $$q = \frac{T_{\infty,1} - T_{\infty,2}}{R_{total}}$$ $$\frac{1}{U} = \frac{1}{h_1} + \sum_{i=1}^n \frac{L_i}{k_i} + \frac{1}{h_2}$$

Hot Gas   | Layer 1 | Layer 2 | Layer 3 | Ambient Air
T_inf,1   |  (k_1)  |  (k_2)  |  (k_3)  |   T_inf,2
  --->    |         |         |         |    --->
   h_1    |   L_1   |   L_2   |   L_3   |     h_2

Interface Temperature Determination

To verify whether a refractory brick layer will exceed its maximum allowable continuous operating temperature, the interface temperature $T_{interface}$ is computed directly by isolating individual resistances:

$$T_{i} - T_{i+1} = q \cdot R_{th,i} = q \left( \frac{L_i}{k_i A} \right)$$

---

4. Cylindrical Conduction & The Log-Mean Area

In circular pipes, tubes, shell-and-tube exchangers, and tubular reactors, heat flows radially. Because the surface area perpendicular to the heat flux increases linearly with radius ($A(r) = 2 \pi r L$), the heat flux $q'' = q / (2 \pi r L)$ decreases outward, creating a logarithmic temperature profile.

Derivation of Cylindrical Conduction

Applying Fourier's law at radius $r$:

$$q = -k A(r) \frac{dT}{dr} = -k (2 \pi r L) \frac{dT}{dr}$$

Separating variables and integrating from inner radius $r_1$ (temperature $T_1$) to outer radius $r_2$ (temperature $T_2$):

$$q \int_{r_1}^{r_2} \frac{dr}{r} = -2 \pi k L \int_{T_1}^{T_2} dT$$ $$q \ln\left(\frac{r_2}{r_1}\right) = 2 \pi k L (T_1 - T_2)$$ $$q = \frac{2 \pi k L (T_1 - T_2)}{\ln(r_2/r_1)}$$

The conductive thermal resistance for a hollow cylindrical shell is:

$$R_{cond,cyl} = \frac{\ln(r_2/r_1)}{2 \pi k L}$$

The resistance per unit length ($R'_{cond}$ in $\text{K}\cdot\text{m/W}$) is:

$$R'_{cond,cyl} = \frac{\ln(r_2/r_1)}{2 \pi k}$$

Radial Temperature Distribution

Integrating from $r_1$ to an arbitrary radius $r$ ($r_1 \le r \le r_2$):

$$T(r) = T_1 - (T_1 - T_2) \frac{\ln(r/r_1)}{\ln(r_2/r_1)}$$

The Log-Mean Area ($A_{lm}$)

To express cylindrical conduction in the familiar planar form ($q = k A \Delta T / \Delta r$):

$$q = \frac{k A_{lm} (T_1 - T_2)}{r_2 - r_1}$$

Equating this to the exact logarithmic equation yields the logarithmic mean area:

$$A_{lm} = \frac{A_2 - A_1}{\ln(A_2 / A_1)} = 2 \pi r_{lm} L$$

Where the log-mean radius is:

$$r_{lm} = \frac{r_2 - r_1}{\ln(r_2 / r_1)}$$

Composite Cylindrical Shells with Convection

For an inner fluid at $T_{\infty,i}$ flowing inside a pipe covered with insulation and exposed to outer ambient air at $T_{\infty,o}$:

$$R'_{total} = \frac{1}{2 \pi r_1 h_i} + \frac{\ln(r_2/r_1)}{2 \pi k_{pipe}} + \frac{\ln(r_3/r_2)}{2 \pi k_{ins}} + \frac{1}{2 \pi r_3 h_o}$$ $$q' = \frac{q}{L} = \frac{T_{\infty,i} - T_{\infty,o}}{R'_{total}}$$

---

5. Spherical Conduction & The Geometric-Mean Area

Spherical geometries are common in chemical storage for pressurized liquefied petroleum gases (LPG, propane, butane), spherical LNG cryogenic tanks, and microencapsulated catalyst pellets. In spherical coordinates, the area through which heat conducts expands quadratically with radius: $A(r) = 4 \pi r^2$.

Derivation of Spherical Conduction

Applying Fourier's law at radius $r$:

$$q = -k A(r) \frac{dT}{dr} = -k (4 \pi r^2) \frac{dT}{dr}$$

Separating variables and integrating from $r_1$ to $r_2$:

$$q \int_{r_1}^{r_2} \frac{dr}{r^2} = -4 \pi k \int_{T_1}^{T_2} dT$$ $$q \left[ -\frac{1}{r} \right]_{r_1}^{r_2} = q \left( \frac{1}{r_1} - \frac{1}{r_2} \right) = 4 \pi k (T_1 - T_2)$$ $$q = \frac{4 \pi k r_1 r_2 (T_1 - T_2)}{r_2 - r_1} = \frac{T_1 - T_2}{\frac{1}{4 \pi k} \left( \frac{1}{r_1} - \frac{1}{r_2} \right)}$$

The conductive thermal resistance of a hollow spherical shell is:

$$R_{cond,sph} = \frac{r_2 - r_1}{4 \pi k r_1 r_2} = \frac{1}{4 \pi k} \left( \frac{1}{r_1} - \frac{1}{r_2} \right)$$

Geometric Mean Area ($A_g$)

Writing spherical conduction in the equivalent planar format:

$$q = \frac{k A_g (T_1 - T_2)}{r_2 - r_1}$$

Equating to the exact equation reveals that the effective area is the geometric mean area:

$$A_g = \sqrt{A_1 A_2} = \sqrt{(4 \pi r_1^2)(4 \pi r_2^2)} = 4 \pi r_1 r_2$$

---

6. The Critical Radius of Insulation Phenomenon

Adding insulation to a flat planar wall always increases total thermal resistance and decreases heat transfer rate. In radial systems (cylinders and spheres), however, adding insulation introduces two competing physical effects:

1. Conductive resistance increases as the insulation thickness increases ($\ln(r/r_o)$ or $1/r_o - 1/r$ grows).
2. Convective resistance decreases because the outer exposed surface area ($A = 2 \pi r L$ or $4 \pi r^2$) expands, enhancing convective dissipation to the ambient fluid.

Analytical Derivation for a Cylinder

Consider a pipe of outer radius $r_o$ covered by an insulation layer of outer radius $r$ ($r \ge r_o$), thermal conductivity $k_{ins}$, and external convective coefficient $h_o$. The total thermal resistance per unit length as a function of outer insulation radius $r$ is:

$$R'_{total}(r) = R'_{internal} + \frac{\ln(r/r_o)}{2 \pi k_{ins}} + \frac{1}{2 \pi r h_o}$$

To find the extremum, differentiate $R'_{total}$ with respect to $r$ and set the derivative to zero:

$$\frac{d R'_{total}}{dr} = \frac{1}{2 \pi k_{ins} r} - \frac{1}{2 \pi h_o r^2} = 0$$ $$\frac{1}{k_{ins} r} = \frac{1}{h_o r^2} \implies r_{crit,cyl} = \frac{k_{ins}}{h_o}$$

Evaluating the second derivative:

$$\frac{d^2 R'_{total}}{dr^2} = -\frac{1}{2 \pi k_{ins} r^2} + \frac{1}{\pi h_o r^3}$$

At $r = r_{crit} = k_{ins}/h_o$:

$$\left. \frac{d^2 R'_{total}}{dr^2} \right|_{r = r_{crit}} = -\frac{h_o^2}{2 \pi k_{ins}^3} + \frac{h_o^2}{\pi k_{ins}^3} = +\frac{h_o^2}{2 \pi k_{ins}^3} > 0$$

Because the second derivative is strictly positive, the thermal resistance is at an absolute MINIMUM at $r = r_{crit}$. Consequently, the heat transfer rate $q$ is at an absolute MAXIMUM at $r = r_{crit}$!

Critical Radius for Spheres

Performing the identical differentiation for a spherical shell of radius $r$:

$$R_{total}(r) = R_{internal} + \frac{1}{4 \pi k_{ins}}\left( \frac{1}{r_o} - \frac{1}{r} \right) + \frac{1}{4 \pi r^2 h_o}$$ $$\frac{d R_{total}}{dr} = \frac{1}{4 \pi k_{ins} r^2} - \frac{2}{4 \pi h_o r^3} = 0 \implies r_{crit,sph} = \frac{2 k_{ins}}{h_o}$$

Practical Engineering Implications

Heat Loss (q)
     ^
     |             Peak Heat Loss at r = r_crit
     |                 /---\
     |                /     \
     |               /       \      Decreasing heat loss
     |  Increasing  /         \     (r > r_crit)
     |  heat loss  /           \---
     |  (r < r_crit)
     +------------+-----------+---------------------> Outer Radius (r)
                 r_o        r_crit

1. If $r_o < r_{crit}$: Adding insulation increases the outer surface area faster than it adds conductive resistance. Heat loss increases until the outer radius reaches $r_{crit}$. Any insulation thickness where $r_o < r < r_{crit}$ backfires if the goal is thermal insulation (though it is exploited intentionally to cool electrical transmission wires and microelectronic leads).
2. If $r_o \ge r_{crit}$: Adding any amount of insulation monotonically decreases heat loss. In industrial chemical plants, process pipes typically have radii $r_o \ge 25\text{ mm}$ ($1\text{ inch}$), while standard insulations (mineral wool, calcium silicate) have $k_{ins} \approx 0.05\text{ W/(m}\cdot\text{K)}$ in ambient air ($h_o \approx 10\text{ W/(m}^2\cdot\text{K)}$), giving $r_{crit} = 0.05 / 10 = 0.005\text{ m} = 5\text{ mm}$. Because $r_o \gg r_{crit}$, industrial process pipe insulation always reduces heat loss.

---

7. Summary Table: Conduction Across Common Geometries

---

8. Comprehensive Worked Numerical Example: High-Pressure Steam Header Insulation Design

Problem Statement

A saturated steam header in a petrochemical plant is fabricated from Schedule 40 carbon steel ($k_{steel} = 45.0\text{ W/(m}\cdot\text{K)}$) with an inner radius $r_1 = 50.0\text{ mm}$ ($0.050\text{ m}$) and outer radius $r_2 = 55.0\text{ mm}$ ($0.055\text{ m}$). High-pressure saturated steam condenses inside the pipe at $T_{steam} = 220.0^\circ\text{C}$ ($493.15\text{ K}$) with an internal condensation heat transfer coefficient $h_i = 3,000\text{ W/(m}^2\cdot\text{K)}$.

The pipe is covered with a $35.0\text{ mm}$ thick jacket of calcium silicate insulation ($k_{ins} = 0.065\text{ W/(m}\cdot\text{K)}$), bringing the outer radius to $r_3 = 55.0 + 35.0 = 90.0\text{ mm}$ ($0.090\text{ m}$). The external ambient air is at $T_\infty = 20.0^\circ\text{C}$ ($293.15\text{ K}$) with an external combined convection-radiation coefficient $h_o = 12.0\text{ W/(m}^2\cdot\text{K)}$. Total pipe length is $L = 10.0\text{ m}$.

Calculate:

1. The critical radius of insulation $r_{crit}$ and confirm whether the insulation reduces heat loss.
2. The individual thermal resistances per unit length ($R'$ in $\text{K}\cdot\text{m/W}$) for internal convection, steel pipe wall, insulation layer, and external convection.
3. The steady-state heat loss rate per meter ($q'$) and the total heat loss rate ($q$) for the $10\text{ m}$ pipe.
4. The temperature at the outer surface of the insulation ($T_{s,o}$) to check compliance with OSHA personnel protection limits ($T_{touch} \le 60^\circ\text{C}$).
5. The heat loss rate per meter if the pipe were completely uninsulated ($q'_{unins}$), and the percentage energy savings achieved by the insulation.

---

Step 1: Critical Radius Check

$$r_{crit} = \frac{k_{ins}}{h_o} = \frac{0.065\text{ W/(m}\cdot\text{K)}}{12.0\text{ W/(m}^2\cdot\text{K)}} = 0.00542\text{ m} = \mathbf{5.42\text{ mm}}$$

Since the outer bare pipe radius $r_2 = 55.0\text{ mm} \gg r_{crit} = 5.42\text{ mm}$, the pipe is well above the critical radius threshold. Adding any insulation will strictly decrease heat loss.

---

Step 2: Thermal Resistances Per Unit Length

1. Internal Convective Resistance ($R'_{conv,i}$):

   $$R'_{conv,i} = \frac{1}{2 \pi r_1 h_i} = \frac{1}{2 \pi (0.050\text{ m})(3,000\text{ W/(m}^2\cdot\text{K)})} = \frac{1}{942.48} = \mathbf{0.001061\text{ K}\cdot\text{m/W}}$$

2. Steel Pipe Wall Conductive Resistance ($R'_{pipe}$):

   $$R'_{pipe} = \frac{\ln(r_2/r_1)}{2 \pi k_{steel}} = \frac{\ln(0.055 / 0.050)}{2 \pi (45.0\text{ W/(m}\cdot\text{K)})} = \frac{\ln(1.100)}{282.743} = \frac{0.095310}{282.743} = \mathbf{0.000337\text{ K}\cdot\text{m/W}}$$

3. Calcium Silicate Insulation Resistance ($R'_{ins}$):

   $$R'_{ins} = \frac{\ln(r_3/r_2)}{2 \pi k_{ins}} = \frac{\ln(0.090 / 0.055)}{2 \pi (0.065\text{ W/(m}\cdot\text{K)})} = \frac{\ln(1.63636)}{0.408407} = \frac{0.492476}{0.408407} = \mathbf{1.20585\text{ K}\cdot\text{m/W}}$$

4. External Convective Resistance ($R'_{conv,o}$):

   $$R'_{conv,o} = \frac{1}{2 \pi r_3 h_o} = \frac{1}{2 \pi (0.090\text{ m})(12.0\text{ W/(m}^2\cdot\text{K)})} = \frac{1}{6.78584} = \mathbf{0.147365\text{ K}\cdot\text{m/W}}$$

5. Total Thermal Resistance Per Unit Length ($R'_{total}$):

   $$R'_{total} = R'_{conv,i} + R'_{pipe} + R'_{ins} + R'_{conv,o}$$ $$R'_{total} = 0.001061 + 0.000337 + 1.20585 + 0.147365 = \mathbf{1.35461\text{ K}\cdot\text{m/W}}$$

Notice that the insulation layer contributes $1.20585 / 1.35461 = 89.0\%$ of the total resistance, while the steel pipe contributes less than $0.03\%$.

---

Step 3: Steady-State Heat Loss Rate

Overall driving temperature difference:

$$\Delta T = T_{steam} - T_\infty = 220.0^\circ\text{C} - 20.0^\circ\text{C} = 200.0\text{ K}$$

Heat loss per unit length ($q'$):

$$q' = \frac{\Delta T}{R'_{total}} = \frac{200.0\text{ K}}{1.35461\text{ K}\cdot\text{m/W}} = \mathbf{147.64\text{ W/m}}$$

Total heat loss for the $10.0\text{ m}$ pipe:

$$q = q' \cdot L = 147.64\text{ W/m} \times 10.0\text{ m} = \mathbf{1,476.4\text{ W}} = \mathbf{1.476\text{ kW}}$$

---

Step 4: Outer Insulation Surface Temperature ($T_{s,o}$)

Using the external convective resistance from the outer surface to ambient air:

$$q' = \frac{T_{s,o} - T_\infty}{R'_{conv,o}} \implies T_{s,o} = T_\infty + q' \cdot R'_{conv,o}$$ $$T_{s,o} = 20.0^\circ\text{C} + (147.64\text{ W/m} \times 0.147365\text{ K}\cdot\text{m/W}) = 20.0^\circ\text{C} + 21.76^\circ\text{C} = \mathbf{41.76^\circ\text{C}}$$

Because $41.8^\circ\text{C} < 60.0^\circ\text{C}$, the insulation thickness satisfies safety requirements to prevent skin burn injury upon incidental contact.

---

Step 5: Uninsulated Heat Loss Comparison

For the bare uninsulated pipe ($r = r_2 = 0.055\text{ m}$):

$$R'_{conv,o,unins} = \frac{1}{2 \pi r_2 h_o} = \frac{1}{2 \pi (0.055\text{ m})(12.0\text{ W/(m}^2\cdot\text{K)})} = \frac{1}{4.14690} = 0.241144\text{ K}\cdot\text{m/W}$$ $$R'_{total,unins} = R'_{conv,i} + R'_{pipe} + R'_{conv,o,unins} = 0.001061 + 0.000337 + 0.241144 = 0.242542\text{ K}\cdot\text{m/W}$$

Uninsulated heat loss per unit length:

$$q'_{unins} = \frac{200.0\text{ K}}{0.242542\text{ K}\cdot\text{m/W}} = \mathbf{824.60\text{ W/m}}$$

Total uninsulated heat loss for $10\text{ m}$:

$$q_{unins} = 8,246.0\text{ W} = 8.246\text{ kW}$$ $$\text{Energy Savings} = \frac{q_{unins} - q}{q_{unins}} \times 100\% = \frac{8,246.0 - 1,476.4}{8,246.0} \times 100\% = \mathbf{82.09\%}$$

---

9. Critical PE Exam Traps & Pitfalls

Trap 1: Confusing Heat Rate ($q$) with Heat Flux ($q''$) in Radial Geometries
In steady-state planar conduction, both heat rate $q$ and heat flux $q'' = q/A$ are constant across every layer. In cylindrical and spherical systems, only the total heat flow rate $q$ is constant. Heat flux $q''(r) = q / A(r)$ drops inversely with $r$ in cylinders and with $r^2$ in spheres. Never equate heat fluxes across concentric radial layers!

Trap 2: The Critical Radius Fallacy for Process Piping
Candidates often memorize $r_{crit} = k/h$ and mistakenly assert that adding the first layer of insulation to a large process pipe will increase heat loss. Remember: $r_o$ of standard process pipes ($25\text{ to }150\text{ mm}$) is almost always much larger than $r_{crit}$ ($5\text{ to }10\text{ mm}$). Heat loss increases only when the bare outer radius is smaller than $r_{crit}$ (thin instrument lines, capillary tubes, and electrical wiring).

Trap 3: Using Arithmetic Mean Area When Radius Ratios are Large
Using $A_{am} = (A_1 + A_2)/2$ instead of $A_{lm} = (A_2 - A_1)/\ln(A_2/A_1)$ introduces significant error when $r_2 / r_1 > 1.5$. On the PE exam, always default to the exact logarithmic expression $\ln(r_2/r_1)/(2\pi k L)$ unless an explicitly thin-walled tube is specified.

Trap 4: Forgetting Contact Resistance at Metal Interfaces
In composite metal walls subjected to high heat flux, microscopic surface roughness creates microscopic air pockets with high thermal resistance. If a problem states a thermal contact resistance $R''_{t,c}$, treat it as an additional series resistance: $R_{t,c} = R''_{t,c} / A$.