5.1 Conduction Heat Transfer & Thermal Insulation Design

5.1 Conduction Heat Transfer & Thermal Insulation Design

Thermal conduction is the transfer of energy from more energetic particles of a substance to adjacent less energetic particles as a result of interactions between particles. In solids, conduction is due to the combination of lattice vibrations (phonons) and free electron transport. In liquids and gases, conduction is due to collisions and diffusion of molecules during their random motion.

---

1. Fundamentals of Thermal Conduction & Fourier's Law

The fundamental law governing rate of heat conduction is Fourier's Law of Heat Conduction, formulated by Joseph Fourier in 1822. For one-dimensional steady-state heat conduction, the heat transfer rate $q$ (expressed in Watts, $\text{W}$) is directly proportional to the temperature gradient $\frac{dT}{dx}$ and the heat transfer area $A$ perpendicular to the direction of heat flow:

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

Where:

• $q$ = Heat conduction rate ($\text{W}$ or $\text{J/s}$)
• $k$ = Thermal conductivity of the material ($\text{W/m}\cdot\text{K}$ or $\text{W/m}\cdot^\circ\text{C}$)
• $A$ = Cross-sectional area normal to heat flow ($\text{m}^2$)
• $\frac{dT}{dx}$ = Temperature gradient in the direction of heat flow ($\text{K/m}$)

The negative sign in Fourier's law enforces the Second Law of Thermodynamics: heat flows naturally in the direction of decreasing temperature (down the temperature gradient, where $\frac{dT}{dx} < 0$).

The conduction heat flux $q''$ (heat rate per unit surface area) is defined as:

$$q'' = \frac{q}{A} = -k \frac{dT}{dx} \quad [\text{W/m}^2]$$

Thermal Conductivity ($k$)

Thermal conductivity $k$ represents a material's capability to conduct heat. It varies significantly across material phases and structures:

---

2. One-Dimensional Steady-State Conduction & Thermal Resistance Networks

Under steady-state conditions without internal heat generation ($q_g = 0$), heat flow through a body is constant. By integrating Fourier's law across simple geometries, we derive standard expressions for heat rate and thermal resistance ($R_{th}$), analogous to Ohm's Law in electrical circuits ($I = \frac{\Delta V}{R_{elec}} \implies q = \frac{\Delta T}{R_{th}}$).

Plane Wall Geometries

For a single plane wall of thickness $L$, thermal conductivity $k$, and surface temperatures $T_{s1}$ and $T_{s2}$ ($T_{s1} > T_{s2}$):

$$q = k A \frac{T_{s1} - T_{s2}}{L} = \frac{T_{s1} - T_{s2}}{R_{cond,plane}}$$

Where conductive thermal resistance $R_{cond,plane}$ is:

$$R_{cond,plane} = \frac{L}{k A} \quad [\text{K/W}]$$

When a wall experiences surface convection from fluids at temperatures $T_{\infty,1}$ and $T_{\infty,2}$ with convective heat transfer coefficients $h_1$ and $h_2$, the convective thermal resistance is:

$$R_{conv} = \frac{1}{h A} \quad [\text{K/W}]$$

For a composite wall consisting of $N$ layers in series exposed to convection on both boundaries, the total thermal resistance $R_{total}$ and heat transfer rate $q$ are:

$$R_{total} = R_{conv,1} + \sum_{i=1}^{N} R_{cond,i} + R_{conv,2} = \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}} = U A (T_{\infty,1} - T_{\infty,2})$$

Where $U = \frac{1}{R_{total} A}$ is the Overall Heat Transfer Coefficient ($\text{W/m}^2\cdot\text{K}$).

Cylindrical Geometries (Pipes and Tubes)

For 1D radial heat conduction through a hollow cylinder of inner radius $r_1$, outer radius $r_2$, length $L$, and thermal conductivity $k$:

$$q = \frac{2 \pi k L (T_{s1} - T_{s2})}{\ln(r_2 / r_1)} = \frac{T_{s1} - T_{s2}}{R_{cond,cyl}}$$ $$R_{cond,cyl} = \frac{\ln(r_2 / r_1)}{2 \pi k L} \quad [\text{K/W}]$$

For a composite pipe insulated on the exterior ($r_1 \to r_2 \to r_3$) exposed to internal fluid $T_{\infty,1}$ ($h_1$) and external fluid $T_{\infty,2}$ ($h_2$):

$$R_{total} = \frac{1}{h_1 (2\pi r_1 L)} + \frac{\ln(r_2/r_1)}{2\pi k_1 L} + \frac{\ln(r_3/r_2)}{2\pi k_2 L} + \frac{1}{h_2 (2\pi r_3 L)}$$

Spherical Geometries

For 1D radial conduction through a hollow sphere of radii $r_1$ and $r_2$:

$$q = \frac{4 \pi k r_1 r_2 (T_{s1} - T_{s2})}{r_2 - r_1} = \frac{T_{s1} - T_{s2}}{R_{cond,sph}}$$ $$R_{cond,sph} = \frac{r_2 - r_1}{4 \pi k r_1 r_2} \quad [\text{K/W}]$$

---

3. Thermal Insulation Design & Critical Radius of Insulation ($r_{cr}$)

Adding insulation layer to a plane wall always increases conductive resistance and decreases heat transfer rate. However, for cylindrical or spherical geometries (e.g., bare electrical wires or small-diameter steam pipes), adding outer insulation increases conductive resistance ($R_{cond} \propto \ln(r_o/r_i)$) but simultaneously increases outer surface area ($A_o = 2 \pi r_o L$), which decreases convective resistance ($R_{conv} = \frac{1}{h A_o}$). This creates a critical insulation radius.

Derivation for Cylindrical Geometry

Consider a pipe of outer radius $r_i$ at temperature $T_i$, covered with insulation of outer radius $r$ and conductivity $k$, exposed to ambient air at $T_\infty$ with convective coefficient $h$. The total thermal resistance per unit length is:

$$R'_{total} = R'_{cond} + R'_{conv} = \frac{\ln(r / r_i)}{2 \pi k} + \frac{1}{h (2 \pi r)}$$

To find the radius $r$ that minimizes total resistance (maximizing heat loss), set $\frac{d R'_{total}}{dr} = 0$:

$$\frac{d R'_{total}}{dr} = \frac{1}{2 \pi k r} - \frac{1}{2 \pi h r^2} = 0 \implies \frac{1}{k} = \frac{1}{h r}$$ $$r_{cr,cyl} = \frac{k}{h} \quad [\text{m}]$$

For spherical geometries, a similar derivation yields:

$$r_{cr,sph} = \frac{2 k}{h} \quad [\text{m}]$$

Engineering Rules for Thermal Insulation Design

1. If $r_i < r_{cr}$: Adding insulation up to $r = r_{cr}$ increases heat transfer. The maximum heat loss occurs precisely at $r = r_{cr}$. Adding insulation beyond $r_{cr}$ reduces heat transfer below the maximum, but heat rate drops below uninsulated level only after outer radius exceeds a value significantly larger than $r_{cr}$.
2. If $r_i \ge r_{cr}$: Adding any amount of insulation immediately decreases heat transfer rate. (Typical industrial steam pipes have $r_i \gg r_{cr}$, so insulation always conserves energy).
3. Electrical Cables: Bare electrical wires are intentionally coated with insulation up to $r_o = r_{cr}$ to maximize heat dissipation from internal $I^2 R$ electrical resistance heating, maintaining lower conductor operating temperatures.

---

4. Extended Surfaces (Fins) & Performance Metrics

Extended surfaces (fins) are primary passive enhancements used to increase heat transfer from a surface by increasing the total effective surface area. They are widely applied in air-cooled engine cylinders, condenser tubes, electronic heatsinks, and refrigeration coils.

Governing Differential Equation for Uniform Fin

Applying an energy balance on a differential element of a pin or straight fin with cross-sectional area $A_c$, perimeter $P$, thermal conductivity $k$, and convective coefficient $h$:

$$\frac{d^2 \theta}{dx^2} - m^2 \theta = 0$$

Where:

• $\theta(x) = T(x) - T_\infty$ = Excess temperature at position $x$
• $m = \sqrt{\frac{h P}{k A_c}}$ = Fin parameter ($\text{m}^{-1}$)

Common Fin Tip Boundary Conditions and Heat Transfer Equations

Note: An adiabatic tip solution can be used for convective tip fins by substituting corrected fin length $L_c = L + \frac{t}{2}$ for rectangular fins or $L_c = L + \frac{D}{4}$ for pin fins.

Fin Performance Quantities

• Fin Effectiveness ($\epsilon_f$): Ratio of fin heat transfer rate to heat transfer rate from the base area without fin:

$$\epsilon_f = \frac{q_f}{h A_{c,b} (T_b - T_\infty)}$$

For an infinitely long fin: $\epsilon_f = \frac{\sqrt{h P k A_c}}{h A_c} = \sqrt{\frac{k P}{h A_c}}$. To justify the cost of adding fins, $\epsilon_f \ge 2$ is required as a general design criterion.

• Fin Efficiency ($\eta_f$): Ratio of actual fin heat transfer rate to maximum possible heat transfer rate if the entire fin surface were maintained at the base temperature $T_b$:

$$\eta_f = \frac{q_f}{q_{max}} = \frac{q_f}{h A_f (T_b - T_\infty)}$$

For an adiabatic tip fin of uniform cross-section: $\eta_f = \frac{\tanh(m L)}{m L}$.

---

5. Transient Heat Conduction & Lumped Capacitance Method

When a solid body experiences a sudden change in its thermal environment (e.g., quenching a hot metal billet in water), its temperature changes with time until thermal equilibrium is established. If internal thermal resistance within the body is negligible compared to external convective resistance across its boundary, the temperature within the body remains spatially uniform ($T(x,y,z,t) \approx T(t)$).

The Biot Number ($Bi$)

The validity of assuming uniform internal temperature is determined by the dimensionless Biot Number:

$$Bi = \frac{\text{Internal Resistance to Conduction}}{\text{External Resistance to Convection}} = \frac{L_c / k}{1 / h} = \frac{h L_c}{k}$$

Where Characteristic Length $L_c$ is:

$$L_c = \frac{V}{A_s}$$

• Plane Wall of thickness $2L$ (cooled from both sides): $L_c = L$
• Long Cylinder of radius $r_o$: $L_c = \frac{\pi r_o^2 L}{2 \pi r_o L} = \frac{r_o}{2}$
• Sphere of radius $r_o$: $L_c = \frac{\frac{4}{3} \pi r_o^3}{4 \pi r_o^2} = \frac{r_o}{3}$

Criterion for Lumped Capacitance Validity: The Lumped Capacitance Method is applicable if and only if $Bi < 0.1$. Errors are typically under 5%.

Lumped Capacitance Temperature Equations

Applying an energy balance on a body of mass $m$, volume $V$, density $\rho$, specific heat $c_p$, initial temperature $T_i$, immersed at $t=0$ into a fluid at $T_\infty$ with convective coefficient $h$:

$$\rho V c_p \frac{dT}{dt} = -h A_s (T - T_\infty)$$

Integrating from $t=0$ ($T=T_i$) to time $t$ ($T=T(t)$):

$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp\left( -\frac{h A_s}{\rho V c_p} t \right) = \exp\left( -\frac{t}{\tau_t} \right)$$

Where $\tau_t$ is the Thermal Time Constant of the system:

$$\tau_t = \left( \frac{1}{h A_s} \right) (\rho V c_p) = R_{conv} C_{th} \quad [\text{s}]$$

The dimensionless form using Fourier Number ($Fo = \frac{\alpha t}{L_c^2}$, where thermal diffusivity $\alpha = \frac{k}{\rho c_p}$) is:

$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp(-Bi \cdot Fo)$$

The total heat transfer $Q(t)$ transferred between $t=0$ and time $t$ is:

$$Q(t) = \rho V c_p (T_i - T(t)) = Q_{max} \left[ 1 - \exp\left( -\frac{t}{\tau_t} \right) \right]$$

---

6. Step-by-Step Worked Engineering Calculation (MELE Board Exam Style)

Problem Statement

A high-pressure steam pipe with an outer diameter of $D_1 = 60\text{ mm}$ ($r_1 = 30\text{ mm} = 0.030\text{ m}$) operates at a surface temperature of $T_{s1} = 250^\circ\text{C}$. The pipe is to be insulated with fiberglass insulation ($k = 0.040\text{ W/m}\cdot\text{K}$). The insulated pipe is exposed to ambient atmospheric air at $T_\infty = 25^\circ\text{C}$ with a convective heat transfer coefficient of $h = 8.0\text{ W/m}^2\cdot\text{K}$.

1. Calculate the critical radius of insulation $r_{cr}$.
2. Determine the heat loss per meter of pipe length ($q/L$) for the uninsulated pipe.
3. Calculate the heat loss per meter ($q/L$) when an insulation layer of thickness $t = 20\text{ mm}$ is added.
4. Determine the percentage reduction in heat loss achieved by adding this insulation layer.

Solution

Step 1: Calculate Critical Radius of Insulation

$$r_{cr} = \frac{k}{h} = \frac{0.040\text{ W/m}\cdot\text{K}}{8.0\text{ W/m}^2\cdot\text{K}} = 0.0050\text{ m} = 5.0\text{ mm}$$

Analysis: Since initial pipe radius $r_1 = 30.0\text{ mm}$ is greater than $r_{cr} = 5.0\text{ mm}$, adding insulation of any thickness will immediately decrease the rate of heat loss.

Step 2: Uninsulated Pipe Heat Loss per Unit Length

For uninsulated pipe, thermal resistance consists only of outer surface convection:

$$R'_{conv,un} = \frac{1}{h (2 \pi r_1)} = \frac{1}{8.0 \times (2 \pi \times 0.030)} = \frac{1}{1.50796} = 0.66315\text{ K}\cdot\text{m/W}$$ $$\left(\frac{q}{L}\right)_{uninsulated} = \frac{T_{s1} - T_\infty}{R'_{conv,un}} = \frac{250 - 25}{0.66315} = \frac{225}{0.66315} = 339.29\text{ W/m}$$

Step 3: Insulated Pipe Heat Loss per Unit Length

With insulation thickness $t = 20\text{ mm}$, outer insulation radius $r_2 = r_1 + t = 30 + 20 = 50\text{ mm} = 0.050\text{ m}$.

Calculate insulation conductive resistance per unit length:

$$R'_{cond,ins} = \frac{\ln(r_2 / r_1)}{2 \pi k} = \frac{\ln(0.050 / 0.030)}{2 \pi \times 0.040} = \frac{\ln(1.6667)}{0.25133} = \frac{0.51083}{0.25133} = 2.03251\text{ K}\cdot\text{m/W}$$

Calculate new outer convective resistance per unit length:

$$R'_{conv,ins} = \frac{1}{h (2 \pi r_2)} = \frac{1}{8.0 \times (2 \pi \times 0.050)} = \frac{1}{2.51327} = 0.39789\text{ K}\cdot\text{m/W}$$

Calculate total thermal resistance per unit length:

$$R'_{total} = R'_{cond,ins} + R'_{conv,ins} = 2.03251 + 0.39789 = 2.43040\text{ K}\cdot\text{m/W}$$

Calculate insulated heat loss rate:

$$\left(\frac{q}{L}\right)_{insulated} = \frac{T_{s1} - T_\infty}{R'_{total}} = \frac{250 - 25}{2.43040} = \frac{225}{2.43040} = 92.58\text{ W/m}$$

Step 4: Percentage Reduction in Heat Loss

$$\text{Reduction (\%)} = \frac{(q/L)_{uninsulated} - (q/L)_{insulated}}{(q/L)_{uninsulated}} \times 100\%$$ $$\text{Reduction (\%)} = \frac{339.29 - 92.58}{339.29} \times 100\% = \frac{246.71}{339.29} \times 100\% = 72.71\%$$

Adding 20 mm of fiberglass insulation reduces steam pipe thermal loss by 72.71%.