Thermal conduction

Thermal conduction transfers heat through microscopic interactions without bulk transport of matter. To predict the heating load through a bar, combine Fourier's law with conservation of energy and the imposed end temperatures. In the isotropic model with positive thermal conductivity $k$, the conductive heat flux points from high to low temperature, opposite to the temperature gradient. A spoon warming along its handle and heat escaping through a wall are familiar examples.

Consider a bar of length $L=0.10\mathrm{m}$, constant cross-sectional area $A=0.010{\mathrm{m}}^2$ and constant thermal conductivity $k=2\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. Maintain its ends at $T_h=40{}^{\circ }\mathrm{C}$ and $T_c=20{}^{\circ }\mathrm{C}$, with $x$ increasing toward the colder end. Assume a steady state, one-dimensional thermal conduction, no internal sources and no heat transfer through the sides. The temperature falls by $20\mathrm{K}$ over $0.10\mathrm{m}$, givingThe difference between the two readings in degree Celsius units is the same $20$-kelvin difference. The positive heat flux carries heat toward increasing $x$.

The general solution under these assumptions isA proof outline is enough to explain the straight line: conservation of energy in any small slice, with no accumulation, sources or side losses, gives $d\dot{Q}/dx=0$. Constant cross-sectional area gives $dq/dx=0$; constant thermal conductivity and Fourier's law then give $d^{2}T/d{x}^2=0$. Two integrations and the end boundary conditions determine the displayed affine function. This is why a visually simple temperature profile can encode a continuous flow of energy. ^[ref]

Doubling the length halves the rate, while doubling the cross-sectional area doubles it, at fixed end temperatures and thermal conductivity. This explains why a long narrow support can limit unwanted heat transfer, while a short broad path can help remove heat. Joining materials leads to thermal resistance in series. During initial warming, internal energy can accumulate and the temperature profile need not be a straight line; the time dependence is described by the heat equation and thermal diffusivity.

Now join two layers, each of length $L_1=L_2=0.10\mathrm{m}$, with constant cross-sectional area $A=0.010{\mathrm{m}}^2$ and stipulated constant thermal conductivities $k_1=1$ and $k_2=4\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. Keep the outer ends at $40$ and $20{}^{\circ }\mathrm{C}$. Assume a steady state, one-dimensional thermal conduction, no internal sources or side losses, and perfect contact: the temperature is continuous at the interface, with no extra contact drop. For each layer, use the area-specific thermal resistance $r_i=L_i/k_i$ and $R_i=r_i/A$. The series calculation gives

The heat flux is continuous: otherwise energy would accumulate at the interface. Yet Fourier's law gives different temperature gradients, $dT/dx=-q/k_1=-160\mathrm{K}{\mathrm{m}}^{-1}$ and $dT/dx=-q/k_2=-40\mathrm{K}{\mathrm{m}}^{-1}$. Lower thermal conductivity requires a steeper temperature gradient to carry the same heat flux; equal lengths therefore give drops in the ratio $4:1$. A kink in the temperature profile records the change of material while perfect contact keeps the temperature itself continuous.

Thermal resistance measures the temperature drop required to drive a given rate of heat transfer. For the steady one-dimensional bar model in thermal conduction,Here $r$ relates the temperature drop to heat flux, while $R_{\mathrm{th}}$ relates it to the rate of heat transfer. The single bar above has $r=0.050{\mathrm{m}}^2\mathrm{K}{\mathrm{W}}^{-1}$ and $R_{\mathrm{th}}=5\mathrm{K}{\mathrm{W}}^{-1}$, so its $20\mathrm{K}$ drop drives $4\mathrm{W}$.

To see why series thermal resistances add, apply Fourier's law in each layer to obtain $\Delta T_i=\dot{Q}{R}_i$. With no accumulation or side losses the rate of heat transfer is the same through both, so adding the drops gives $T_h-T_c=\dot{Q}(R_1+R_2)$. This lets us size layered walls without solving a new differential equation for every layer. The structure matches Ohm's law: temperature difference plays the role of electric potential difference, the rate of heat transfer plays the role of electric current, and thermal resistance plays the role of electrical resistance. ^[ref]

In the absence of heat advection, an insulating sphere of radius $a$ in a uniform remote temperature gradient $\mathbf{H}$ has the displayed exterior harmonic function as its temperature field. The reflected dipole term makes ${\partial }_{r}T=0$ at $r=a$. The surface temperature is $T_0+(3/2)a\mathbf{H}\cdot \mathbf{n}$, so the tangential gradient is amplified by $3/2$. If the sphere has thermal conductivity ratio $K$ relative to its surroundings, temperature and heat-flux continuity instead give surface gradient factor $3/(2+K)$.

With fixed pressure $p$, let $A=\mu p/\mathcal{R}$, so $\rho =A/T$. A material mass coordinate has ${\partial }_x=\rho {\partial }_m$ and the material time derivative becomes ${\partial }_t{|}_m$. The energy and continuity equations then give the displayed equation with $\tau =Ct$ and $C=(\gamma -1)(\mu /\mathcal{R}{)}^{2}p/\gamma$. A mass integral from a fixed origin is material only when the mass flux there vanishes; otherwise subtract its time-integrated flux to label parcels. Uniform pressure is a low-Mach approximation when sound crossing is faster than thermal evolution.

For conductivity $\lambda ={\lambda }_{0}T$, no cooling in the hot phase and a zero-temperature boundary at $m=0$, the mass-coordinate equation is linear diffusion. Its similarity profile solves $f^{\prime \prime }+\xi f^{\prime }/2=0$, $f(0)=0$, $f(\infty )=1$, giving the displayed error function. The thermal layer broadens as the square root of time. Conducted energy enters the cold radiating reservoir at rate per area $L={\lambda }_{0}A{T}_m(0)=A{T}_0\sqrt{{\lambda }_0}/\sqrt{\pi \tau }$, hence $L\propto t^{-1/2}$.

The temperature gradient is the gradient $\nabla T$ of the temperature field. It points in the direction of fastest temperature increase.

Fourier's law predicts the local conductive heat flux from the local temperature gradient. In a continuum approximation with isotropic thermal conductivity $k>0$, its vector form and one-dimensional form areThe heat flux $q$ has units $\mathrm{W}{\mathrm{m}}^{-2}$, the temperature gradient has units $\mathrm{K}{\mathrm{m}}^{-1}$ and the thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. The rate of heat transfer $\dot{Q}$ has units $\mathrm{W}$. These units distinguish a local material response from the total power crossing a surface.

The minus sign makes heat flow toward lower temperature: for $k>0$, $\mathbf{q}\cdot \nabla T=-k|\nabla T{|}^2\le 0$. In the worked bar in thermal conduction, $dT/dx=-200\mathrm{K}{\mathrm{m}}^{-1}$ therefore gives $q=+400\mathrm{W}{\mathrm{m}}^{-2}$. Reversing the end temperatures reverses the heat flux. A zero temperature gradient gives zero conductive heat flux in this model.

Fourier's law is a constitutive relation expressing measured material behavior; it is not a consequence of conservation of energy alone. It applies locally and does not require a steady state. The straight-line bar solution additionally assumes constant thermal conductivity, constant cross-sectional area, no sources or side losses, and fixed end temperatures. Anisotropic materials require the directional response described by the thermal conductivity tensor. The scalar $k$ used in the examples expresses the stated isotropy assumption. ^[ref]

Historically, Joseph Fourier's Théorie analytique de la chaleur was published in 1822. In Chapter II, article 112, he writes the heat crossing a spherical surface as the product of thermal conductivity, elapsed time, surface area and the negative radial temperature gradient; article 113 balances the incoming and outgoing heat to obtain the heat equation for a sphere. The enduring idea is this connection between a local flow rule and a global prediction of changing temperature. ^{[ref] [ref]}

Thermal conductivity is the material response coefficient in Fourier's law. For an isotropic material it is a scalar $k$, with units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$; larger positive $k$ means a larger conductive heat flux for the same temperature gradient. It can depend on temperature and composition. Direction-dependent response in an anisotropic material is described by a thermal conductivity tensor. The two worked examples use stipulated positive constants, rather than measured values for named materials.

Thermal conductivity describes material response, whereas thermal resistance also includes geometry: for the steady bar, $R_{\mathrm{th}}=L/(kA)$. Doubling $k$ at fixed length, cross-sectional area and end temperatures halves thermal resistance and doubles the rate of heat transfer. Doubling length changes thermal resistance without changing the material's thermal conductivity. This distinction lets a designer improve cooling either by choosing a higher thermal conductivity or by shortening and widening the path. It also explains why the less conductive layer in thermal resistance takes the greater temperature drop.

Thermal diffusivity answers a different question: how rapidly a temperature disturbance spreads. Its usual model $\kappa =k/(\rho c_p)$ also involves mass density $\rho$ and specific heat capacity $c_p$, so equal thermal conductivities need not imply equal rates of warming.

The thermal conductivity tensor maps a temperature gradient to the negative heat-flux density. In passive reciprocal conduction it is symmetric positive definite. This permits direction-dependent transport while ensuring positive dissipated energy.

Thermal diffusivity measures how rapidly thermal conduction smooths a temperature difference inside a material. Thermal conductivity measures the heat flux driven by a specified temperature gradient; thermal diffusivity also accounts for the heat capacity that must be warmed. A material can carry substantial heat yet change temperature slowly because it stores substantial internal energy per kelvin.

Write $\alpha =\kappa$ for thermal diffusivity, reserving $k$ for thermal conductivity. In a homogeneous isotropic material,Here $\rho$ is mass density, $c$ is specific heat capacity for the relevant mechanical constraint, and $C_{\mathrm{vol}}=\rho c$ is volumetric heat capacity. The familiar $k/(\rho c_p)$ uses specific heat capacity at constant pressure. For ordinary solids with negligible thermal expansion, the difference between the two constrained specific heat capacities is often neglected. A rigid fixed-volume model instead uses specific heat capacity at constant volume; a compressible moving fluid requires the accompanying work and transport terms. Thus the constraint is part of the approximation, not a choice made solely from the units. ^[ref]

Assume a stationary material in the continuum approximation, local thermal equilibration on times shorter than those resolved, no internal heating, no phase transition, and negligible convection and mechanical work. Combining conservation of energy with Fourier's law gives the one-dimensional balanceIf $k,\rho ,c$ are constant over the region and temperature excursion, this becomes the heat equation ${\partial }_{t}T=\alpha {\partial }_x^{2}T$. If $k$ varies, retain the differentiated product; replacing it with $k{\partial }_x^{2}T$ loses a term. Anisotropic materials need a thermal conductivity tensor. The scalar constant-$\alpha$ examples below deliberately exclude these complications. ^[ref]

For an original comparison, stipulate two model materials with equal thermal conductivity $k=2\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$ and equal mass density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, but specific heat capacities $c_A=1000$ and $c_B=4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. ThenAt the same temperature gradient, both have the same instantaneous heat flux by Fourier's law. At the same initial temperature profile and boundary conditions, material B takes four times as long to smooth a given spatial pattern. With fixed end temperatures, identical geometry and no sources, both eventually have the same steady profile and heat flux. The larger heat capacity changes the transient, not this steady state. These are stipulated parameters, not measurements of named substances.

The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity. It compares the rates at which momentum and temperature variations diffuse. In thermal-diffusion units the Boussinesq equations have momentum inertia multiplied by ${\mathrm{Pr}}^{-1}$; large, rather than small, Prandtl number supplies the usual instantaneous Stokes flow limit.

The thermal diffusion time $t_D=L^2/\alpha$, with $\alpha =\kappa$ the thermal diffusivity, is a characteristic scale for smoothing a temperature variation over length $L$. It is the thermal instance of the diffusive time scale. Rescaling $x=L\hat{x}$ and $t=(L^2/\alpha )\hat{t}$ removes $L$ and $\alpha$ from the constant-coefficient heat equation. Doubling length therefore multiplies comparable internal equilibration times by four; increasing thermal diffusivity by four divides them by four.

This scale is not a universal time of complete thermal equilibrium. A nonzero Fourier mode decays exponentially and never becomes exactly zero at finite time. The geometry, boundary conditions, initial pattern and accepted residual determine the numerical factor. For example, an insulated bar's slowest nonuniform Fourier mode has decay time $L^2/({\pi }^2\alpha )$, and reducing its amplitude to a fraction $\epsilon$ requiresExternal heat transfer can impose a longer cooling time even when internal temperature differences smooth quickly. Thermal diffusivity alone does not determine cooling to the room's temperature.

With thermal diffusivity $\xi$, a Fourier mode of angular wavenumber $k$ decays at rate $\beta =\xi k^2$. To avoid confusing this $k$ with thermal conductivity, use $\alpha =\xi$ and angular wavenumber $K$ below: the same rate is $\beta =\alpha K^2$, with units ${\mathrm{s}}^{-1}$.

Consider an insulated bar $0\le x\le L$ with constant thermal diffusivity, no sources and no heat flux through either end or the sides. An original initial temperature profile and its exact evolution areThe cosine has zero temperature gradient at both ends, satisfying the insulated Neumann boundary conditions. Its spatial mean is zero, so conservation of energy leaves the mean temperature $\overline{T}$ unchanged while the nonuniform part disappears. A short derivation substitutes $T=\overline{T}+A(t)\cos (K_{n}x)$ into the heat equation: the second spatial derivative supplies $-K_n^2$, hence $A^{\prime }=-\alpha K_n^{2}A$. The exponential decay time is ${\tau }_n=1/(\alpha K_n^2)$; doubling angular wavenumber quarters this time. If cycles-per-length spatial frequency $f_x$ is used instead, $K=2\pi f_x$ and the rate is $4{\pi }^2\alpha f_x^2$.

For the two model materials in thermal diffusivity and $L=0.010\mathrm{m}$, ${\tau }_{1,A}\approx 5.07\mathrm{s}$ and ${\tau }_{1,B}\approx 20.3\mathrm{s}$. After $10\mathrm{s}$, the remaining amplitude fractions are approximately $0.139$ and $0.610$. A narrow hot-and-cold pattern disappears faster than a broad one because its temperature curvature, and hence the imbalance in conductive heat flux, is larger.

The thermal diffusion length $\ell (t)=\sqrt{\alpha t}$ is the characteristic length over which thermal conduction changes a temperature distribution during time $t$ in the constant-thermal diffusivity model. It is a scale, not a sharp front or distance travelled at a constant speed: the heat equation has infinite propagation speed mathematically. Different definitions of penetration depth introduce numerical factors. In particular, the normalized one-dimensional heat kernel has spatial variance $2\alpha t$, so its root-mean-square width is $\sqrt{2\alpha t}$.

At $t=10\mathrm{s}$, the model materials have thermal diffusion lengths $4.47$ and $2.24\mathrm{mm}$. For $L=10\mathrm{mm}$ their thermal diffusion times are $50$ and $200\mathrm{s}$. These values describe the same heat equation as the thermal diffusion rate of a Fourier mode; the different numerical factors come from the definition of the scale.

Laser flash analysis measures thermal diffusivity by heating one face of a specimen briefly and recording the opposite face's temperature rise. The original 1961 experiment used a flash tube rather than a laser: a blackened thin metal specimen sat in an opaque ceramic holder, with thermocouple wires pressed against its rear face. An amplifier, oscilloscope and camera recorded the signal. Uniform front illumination and the black coating promoted uniform absorption; the opaque low-thermal conductivity holder blocked direct rear illumination and reduced edge losses. The delayed rear response reveals internal thermal conduction rather than direct heating of the sensor. ^[ref]

A modern version uses a laser or flash lamp, a specimen of measured thickness $L$, controlled initial temperature, and a rear-facing infrared detector with a synchronized recorder. The readout is a temperature rise versus time, starting from a pre-pulse baseline. Thermal conductivity additionally requires mass density and specific heat capacity; it is not obtained from timing alone. ^[ref]

Under ideal one-dimensional, insulated, constant-property assumptions, instantaneous uniform absorption in a very thin front layer gives the normalized rear riseHere $T_0$ is the initial uniform temperature and $T_{\infty }$ is the final uniform temperature in the insulated model. Fitting this curve estimates thermal diffusivity; each term has the same $\alpha K_n^2$ decay as the insulated-bar Fourier modes. Equivalently, the half-rise time $t_{1/2}$ solves $R(t_{1/2})=1/2$. Writing $s=\alpha t/L^2$, the series gives $s_{1/2}\approx 0.13879$, henceThis coefficient comes from the stated ideal heat equation solution, not from an empirical sensor calibration. It has been independently checked by solving the series equation. Absolute absorbed pulse energy cancels from the normalized curve; determining heat capacity from its absolute height needs additional information. ^[ref]

Parker and colleagues identify finite flash duration, heat loss, nonuniform illumination, recording bandwidth and electrical interference as practical limits. A specimen too thin responds during the pulse; one too thick allows losses before the rear rise finishes. They shielded the electrical chain and recorded a baseline; their reported tests were preliminary metal measurements near ambient and moderately elevated temperatures. ^[ref]

Useful additional controls are to measure thickness at several points, verify the recorder's time scale and detector response independently, record a no-pulse baseline, and repeat at different pulse energies while avoiding phase transitions. Compare different thicknesses of the same homogeneous material: ideal half-rise times scale as $L^2$. Failure of these checks can expose finite-pulse effects, sensor delay, optical transmission, coating storage or edge losses; a thickness error is especially consequential because the inferred thermal diffusivity depends on $L^2$. These are proposed diagnostic controls, not invented measurements. When the pulse is not short or the sample exchanges appreciable heat with its surroundings, fit a model that includes those effects instead of applying the ideal half-rise coefficient unchanged.