9.1 Conduction, Convection, and Radiation Fundamentals

Key Takeaways

  • Fourier's law of conduction states that heat flux is proportional to the negative temperature gradient: q = -k dT/dx
  • Newton's law of cooling governs convection: Q = h A (T_s - T_\infty), where h is flow-dependent and not a fluid property
  • Key convection dimensionless numbers: Nusselt (convective to conductive ratio), Prandtl (momentum to thermal diffusivity), and Grashof (buoyancy to viscous force)
  • Stefan-Boltzmann law defines blackbody emission: E_b = \sigma T^4, requiring absolute temperatures in Kelvin or Rankine
  • View factors (F_ij) represent geometry fractions for radiation exchange and obey reciprocity (A_i F_ij = A_j F_ji) and summation rules
Last updated: July 2026

9.1 Conduction, Convection, and Radiation Fundamentals

Heat transfer is thermal energy in transit due to a spatial temperature difference. The study of heat transfer differs from thermodynamics: while thermodynamics deals with the end states and the quantity of energy required to transition between them, heat transfer predicts the rate at which energy is exchanged and the temperature distribution within a system. For chemical engineers, mastering the three fundamental modes of heat transfer—conduction, convection, and radiation—is essential for designing reactors, separators, and utility systems, and is a major focus of the FE Chemical exam.

Conduction & Fourier's Law

Conduction is the transfer of energy from more energetic particles of a substance to adjacent, less energetic particles due to interactions between them. In gases and liquids, conduction is due to the collisions and diffusion of molecules during their random motion. In solids, conduction is attributed to lattice vibrations (phonons) and the translation of free electrons.

The fundamental rate equation for conduction is Fourier’s Law. In its general vector form, it states that the heat flux is proportional to the negative gradient of temperature:

q=kT\vec{q}'' = -k \nabla T

where $\vec{q}''$ is the heat flux vector ($\text{W/m}^2$), and $k$ is the thermal conductivity ($\text{W/(m·K)}$), which is a transport property of the material. For steady-state, one-dimensional heat conduction in a Cartesian coordinate system (e.g., through a plane wall of thickness $L$ and surface area $A$), Fourier's Law integrates to:

Q˙=kA(T1T2)L\dot{Q} = \frac{kA(T_1 - T_2)}{L}

where $\dot{Q}$ is the heat transfer rate ($\text{W}$), and $T_1$ and $T_2$ are the temperatures at the boundaries ($T_1 > T_2$).

Cylindrical Geometries

For radial conduction through a cylindrical wall (such as pipe walls or pipe insulation), the heat transfer area changes with radius ($A = 2\pi r L$). Integrating Fourier's Law in cylindrical coordinates yields:

Q˙=2πkL(T1T2)ln(r2/r1)\dot{Q} = \frac{2\pi k L (T_1 - T_2)}{\ln(r_2/r_1)}

where $L$ is the length of the cylinder, and $r_1$ and $r_2$ are the inner and outer radii, respectively.

Thermal Resistance

By analogy to Ohm's law ($I = V/R$), we can define a thermal resistance for conduction ($R_{cond}$):

  • Plane wall conduction resistance: $R_{cond} = \frac{L}{kA}$
  • Cylindrical shell conduction resistance: $R_{cond, cyl} = \frac{\ln(r_2/r_1)}{2\pi k L}$

Convection & Newton's Law of Cooling

Convection heat transfer is comprised of two mechanisms: energy transfer due to random molecular motion (diffusion) and energy transfer due to the bulk, or macroscopic, motion of the fluid (advection). When a fluid moves over a surface at a different temperature, the heat transfer is governed by Newton's Law of Cooling:

Q˙=hA(TsT)\dot{Q} = h A (T_s - T_\infty)

where $h$ is the convective heat transfer coefficient ($\text{W/(m}^2\text{·K)}$), $T_s$ is the surface temperature, and $T_\infty$ is the bulk fluid temperature. The coefficient $h$ is not a thermodynamic property of the fluid; instead, it depends on surface geometry, the nature of fluid motion, fluid properties, and bulk fluid velocity.

Convection is broadly classified into two categories:

  1. Forced Convection: The flow is caused by external means, such as a pump, fan, or atmospheric winds.
  2. Natural (Free) Convection: The flow is induced by buoyancy forces within the fluid, which arise from density differences caused by temperature variations.

Convective Dimensionless Numbers

To characterize convective heat transfer and scale experimental data, chemical engineers rely on several key dimensionless groups:

Dimensionless NumberSymbolFormulaPhysical Interpretation
Nusselt Number$Nu$$\frac{h L_c}{k_f}$Ratio of convective heat transfer to conductive heat transfer within the fluid boundary layer.
Prandtl Number$Pr$$\frac{C_p \mu}{k} = \frac{\nu}{\alpha}$Ratio of momentum diffusivity ($\nu$) to thermal diffusivity ($\alpha$). Describes the relative thickness of velocity and thermal boundary layers.
Grashof Number$Gr$$\frac{g \beta (T_s - T_\infty) L_c^3}{\nu^2}$Ratio of buoyancy forces to viscous forces in the fluid. Governs natural convection.
Rayleigh Number$Ra$$Gr \cdot Pr$Combined buoyancy and transport parameter. Determines the transition from laminar to turbulent flow in natural convection.
Reynolds Number$Re$$\frac{\rho v L_c}{\mu}$Ratio of inertial forces to viscous forces. Determines flow regime (laminar vs. turbulent) in forced convection.

For forced convection, correlations typically take the form $Nu = f(Re, Pr)$. For natural convection, correlations take the form $Nu = f(Gr, Pr)$ or $Nu = f(Ra)$.

Radiation & Stefan-Boltzmann Law

Radiation is energy emitted by matter that is at a non-zero temperature. Unlike conduction and convection, radiation does not require the presence of a material medium and propagates at the speed of light via electromagnetic waves.

A blackbody is an idealized physical body that absorbs all incident electromagnetic radiation and acts as a perfect emitter. The emissive power of a blackbody is given by the Stefan-Boltzmann Law:

Eb=σT4E_b = \sigma T^4

where $T$ is the absolute temperature ($\text{K}$ or $R$) and $\sigma$ is the Stefan-Boltzmann constant:

  • $\sigma = 5.670 \times 10^{-8} \text{ W/(m}^2\text{·K}^4\text{)}$ in SI units
  • $\sigma = 0.1714 \times 10^{-8} \text{ Btu/(h·ft}^2\text{·R}^4\text{)}$ in USCS units

Real surfaces emit less energy than a blackbody at the same temperature. This is accounted for by the emissivity ($\epsilon$), a dimensionless surface property ($0 \le \epsilon \le 1$):

E=ϵσT4E = \epsilon \sigma T^4

View Factors

When analyzing radiation exchange between multiple surfaces, we must account for their orientation relative to one another. The view factor (or shape factor), $F_{ij}$, is defined as the fraction of the radiation leaving surface $i$ that is directly intercepted by surface $j$.

Two fundamental rules govern view factors in an enclosure of $N$ surfaces:

  1. Reciprocity Relation: AiFij=AjFjiA_i F_{ij} = A_j F_{ji}
  2. Summation Rule: j=1NFij=1\sum_{j=1}^N F_{ij} = 1

For a flat or convex surface, none of the radiation leaving the surface can strike itself, so $F_{ii} = 0$. For a concave surface, some radiation strikes the surface itself, so $F_{ii} > 0$.

Worked Examples

Worked Example 1: Cylindrical Heat Conduction

A carbon steel pipe ($k = 50 \text{ W/(m·K)}$) with an inner radius $r_1 = 0.04 \text{ m}$ and an outer radius $r_2 = 0.05 \text{ m}$ is covered with a layer of fiberglass insulation ($k = 0.04 \text{ W/(m·K)}$) of thickness $0.03 \text{ m}$ ($r_3 = 0.08 \text{ m}$). The inner surface of the steel pipe is maintained at $220^\circ\text{C}$, and the outer surface of the insulation is at $40^\circ\text{C}$. Calculate the rate of heat loss per meter of pipe length.

Step 1: Identify the thermal resistances in series.

The heat must conduct radially through the pipe wall and then through the insulation. We calculate the resistances per unit length ($R'$):

Rcond,pipe=ln(r2/r1)2πkpipe=ln(0.05/0.04)2π×50=0.2231314.160.00071 mK/WR'_{cond, pipe} = \frac{\ln(r_2/r_1)}{2\pi k_{pipe}} = \frac{\ln(0.05 / 0.04)}{2\pi \times 50} = \frac{0.2231}{314.16} \approx 0.00071 \text{ m}\cdot\text{K/W} Rcond,ins=ln(r3/r2)2πkins=ln(0.08/0.05)2π×0.04=0.47000.25131.870 mK/WR'_{cond, ins} = \frac{\ln(r_3/r_2)}{2\pi k_{ins}} = \frac{\ln(0.08 / 0.05)}{2\pi \times 0.04} = \frac{0.4700}{0.2513} \approx 1.870 \text{ m}\cdot\text{K/W}

Step 2: Sum the resistances to find the total thermal resistance per unit length.

Rtotal=Rcond,pipe+Rcond,ins=0.00071+1.870=1.8707 mK/WR'_{total} = R'_{cond, pipe} + R'_{cond, ins} = 0.00071 + 1.870 = 1.8707 \text{ m}\cdot\text{K/W}

(Note that the thermal resistance of the metal pipe is negligible compared to the insulation layer, which is typical in industrial piping).

Step 3: Apply the radial heat rate equation.

q=TinnerTouterRtotal=220401.870796.2 W/mq' = \frac{T_{inner} - T_{outer}}{R'_{total}} = \frac{220 - 40}{1.8707} \approx 96.2 \text{ W/m}

The rate of heat loss per meter of pipe length is $96.2 \text{ W/m}$.

Worked Example 2: View Factor Summation and Reciprocity

A long, triangular duct has an equilateral cross-section. Surface 1 is the bottom flat plate, Surface 2 is the left plate, and Surface 3 is the right plate. Find all nine view factors $F_{ij}$ assuming the duct is infinitely long.

Step 1: Identify self-view factors.

Since all three surfaces are flat plates, none of them can see themselves:

F11=F22=F33=0F_{11} = F_{22} = F_{33} = 0

Step 2: Apply the summation rule for Surface 1.

F11+F12+F13=1    F12+F13=1F_{11} + F_{12} + F_{13} = 1 \implies F_{12} + F_{13} = 1

By symmetry of the equilateral shape, the fraction of radiation going to Surface 2 must equal that going to Surface 3:

F12=F13=0.5F_{12} = F_{13} = 0.5

Step 3: Determine the remaining view factors by symmetry.

Because the cross-section is an equilateral triangle and the duct is symmetric:

F21=F23=0.5F_{21} = F_{23} = 0.5 F31=F32=0.5F_{31} = F_{32} = 0.5

Thus, the complete matrix of view factors is:

F=[0&0.5&0.50.5&0&0.50.5&0.5&0]F = \begin{bmatrix} 0 \& 0.5 \& 0.5 \\ 0.5 \& 0 \& 0.5 \\ 0.5 \& 0.5 \& 0 \end{bmatrix}

Test Your Knowledge

The Prandtl number (Pr) is a dimensionless parameter that represents the ratio of momentum diffusivity to thermal diffusivity. Which of the following fluids would be expected to have the lowest Prandtl number at room temperature?

A
B
C
D
Test Your Knowledge

A hemispherical dome (surface 1) of radius R is placed on a flat circular base (surface 2) of the same radius. What is the view factor F11 (the fraction of radiation leaving the dome that is intercepted by the dome itself)?

A
B
C
D
Test Your Knowledge

A plane wall of thickness L = 0.1 m has a thermal conductivity of k = 1.5 W/(m·K). The left surface is maintained at 120°C and the right surface is exposed to convection with a fluid at 20°C and a convective heat transfer coefficient of h = 15 W/(m²·K). Under steady-state conditions, what is the temperature of the right surface of the wall?

A
B
C
D