Heat Flow and Thermodynamics: From Fourier’s Law to Energy BalanceHeat flow is a fundamental concept that links microscopic particle motion with macroscopic energy transfer. It underpins technologies from building insulation to power plants, and it provides the bridge between classical thermodynamics and transport phenomena. This article explains heat flow mechanisms, derives and interprets Fourier’s law, places heat flow within thermodynamic frameworks, and applies these ideas to practical energy balance problems.
1. Heat and Heat Flow: Definitions and physical meaning
Heat is energy transferred between systems (or system and surroundings) because of a temperature difference. Unlike internal energy, heat is not a property of a system; it is a transient form of energy in transfer. Heat flow (or heat flux) quantifies the rate at which thermal energy crosses a surface.
- Heat (symbol Q) — energy transferred; units: joules (J).
- Heat transfer rate or power (Q̇) — energy per unit time; units: watts (W).
- Heat flux (q) — heat transfer per unit area; units: W/m².
- Temperature (T) — a scalar field describing thermal state; units: kelvin (K).
Heat flows spontaneously from regions of higher temperature to lower temperature until thermal equilibrium is achieved. The microscopic cause is the exchange of kinetic and potential energy between particles (molecules, electrons, phonons) across the interface or within a medium.
2. Mechanisms of heat transfer
There are three classical modes of heat transfer:
- Conduction: energy transfer within or between stationary bodies due to microscopic collisions and energy exchange. Dominant in solids and stagnant fluids.
- Convection: energy transfer by bulk motion of a fluid (natural or forced), combining conduction and advection.
- Radiation: transfer of energy by electromagnetic waves; no medium required.
Conduction is the primary focus when discussing Fourier’s law and many energy-balance derivations; convection and radiation require additional relations (Newton’s law of cooling, Stefan–Boltzmann law) and often couple with conduction in practical problems.
3. Fourier’s law of heat conduction
Fourier’s law provides a constitutive relation linking heat flux to the temperature gradient in a material. In its most common differential form for steady or transient conduction:
q = −k ∇T
where:
- q is the heat flux vector (W/m²),
- k is the thermal conductivity (W/m·K), a material property,
- ∇T is the temperature gradient (K/m),
- the negative sign indicates heat flows from high to low temperature.
In one dimension along x:
q_x = −k (dT/dx)
For isotropic materials, k is a scalar; for anisotropic materials k becomes a tensor, and Fourier’s law generalizes to qi = −k{ij} ∂T/∂x_j.
Fourier’s law is empirically based but consistent with microscopic transport theories (e.g., kinetic theory for gases, phonon transport in solids). Its validity depends on the continuum hypothesis and that local thermodynamic equilibrium holds over the scales of interest.
4. The heat equation: conservation + Fourier
Combine Fourier’s law with energy conservation to derive the heat (diffusion) equation. Consider a control volume with density ρ, specific heat capacity c_p (or c for constant-volume), and internal heat generation per unit volume Q̇_gen (W/m³). Conservation of energy yields:
ρ c ∂T/∂t = ∇·(k ∇T) + Q̇_gen
For constant thermal conductivity and no internal generation:
∂T/∂t = α ∇²T
where α = k / (ρ c) is the thermal diffusivity (m²/s). This parabolic partial differential equation governs transient conduction problems: diffusion of thermal energy in space and time.
Boundary conditions typically specify temperature (Dirichlet), heat flux (Neumann), or convective/radiative exchange (Robin/mixed).
5. Thermal resistance and steady-state conduction
For steady 1-D conduction through a plane wall of thickness L and area A, with constant k and temperatures T1 and T2 on opposite faces, the heat transfer rate is:
Q̇ = (k A / L) (T1 − T2)
This is analogous to Ohm’s law; define thermal resistance R_th = L / (k A), then Q̇ = (T1 − T2) / R_th. Composite walls, cylindrical shells (pipes), and spherical shells have corresponding resistance formulas—useful for engineering calculations and lumped parameter models.
6. Convection and radiation: coupling with conduction
Convection at a solid boundary is often modeled by Newton’s law of cooling:
q_conv = h (T_surface − T_∞)
where h is the convective heat transfer coefficient (W/m²K), dependent on flow conditions, fluid properties, and geometry. Combine with conduction in the solid using boundary conditions: −k (∂T/∂n)|_surface = h (T_surface − T_∞).
Radiation exchange between surfaces follows the Stefan–Boltzmann law for a blackbody: q_rad = σ T^4. For real surfaces and between surfaces, view factors and emissivity ε modify the net exchange. In many engineering problems, linearized radiation (approximation around an operating temperature) is used to combine radiation with convection in a mixed boundary condition.
7. Thermodynamics perspective: entropy and irreversible heat flow
From thermodynamics, heat flow is an irreversible process that increases total entropy. For a small heat δQ transferred reversibly at temperature T, the change in entropy is δS = δQ_rev / T. For irreversible heat transfer between two reservoirs at T_hot and T_cold, entropy production σ > 0:
σ = Q̇ (1/T_cold − 1/T_hot) ≥ 0
Local non-equilibrium thermodynamics links heat flux to thermodynamic forces: Fourier’s law can be seen as a linear phenomenological law where heat flux is proportional to the gradient of the intensive variable (temperature). The Onsager reciprocal relations extend these ideas to coupled transport processes (thermoelectric effects, thermal diffusion).
8. Energy balances: formulating practical problems
Energy balances are statements of conservation applied to a control mass or control volume. General form for a control volume in engineering:
dE_cv/dt = Q̇_in − Q̇_out + Ẇ_in − Ẇ_out + Σṁ_in h_in − Σṁ_out h_out + Q̇_gen
For pure heat conduction problems with negligible mass flow and work, simplify to:
ρ c ∂T/∂t = ∇·(k ∇T) + Q̇_gen
Common practical steps:
- Define system and control volume.
- Choose appropriate assumptions (steady vs transient, 1-D vs multi-D, constant properties).
- Write governing equations (heat equation, boundary/initial conditions).
- Solve analytically (separation of variables, similarity solutions, Laplace transforms) or numerically (finite difference, finite element, finite volume).
- Validate and interpret—compute heat rates, temperature profiles, and check energy conservation.
Examples:
- Cooling of a heated rod (transient 1‑D conduction).
- Heat loss through building walls (steady conduction with convection boundary conditions).
- Thermal design of heat exchangers (coupled convection–conduction analysis).
9. Lumped-capacitance model
When a solid object has Biot number Bi = hLc/k << 0.1 (Lc is characteristic length), internal temperature gradients are negligible and the object can be treated as isothermal. The lumped model yields:
dT/dt = −(h A)/(ρ c V) (T − T_∞) => T(t) − T_∞ = (T(0) − T_∞) e^(−t/τ)
where τ = ρ c V / (h A) is the thermal time constant. This simplifies transient analysis and gives quick engineering estimates.
10. Microscale and non-classical effects
At micro/nanoscale or in rarefied gases, Fourier’s law may fail. Examples:
- Ballistic transport and phonon mean free paths comparable to device dimensions cause non-local heat transport.
- Time-dependent heat conduction with finite propagation speed is modeled by hyperbolic heat equations (Cattaneo–Vernotte) addressing the paradox of infinite speed in Fourier’s parabolic equation: τ_q ∂q/∂t + q = −k ∇T
where τ_q is a relaxation time.
- Thermoelectric coupling (Seebeck, Peltier) couples heat and charge transport; constitutive relations include cross-coefficients.
11. Numerical methods and practical computation
Most real-world heat transfer problems require numerical solution. Common approaches:
- Finite difference method (FDM): simple grids, good for structured domains.
- Finite volume method (FVM): conserves fluxes; widely used in CFD and conjugate heat transfer.
- Finite element method (FEM): flexible for complex geometries and variable properties.
Key numerical concerns: mesh resolution near boundaries, stability/time-stepping for transient problems, treatment of nonlinearities (temperature-dependent k, radiation), and coupling between conduction, convection, and fluid flow.
12. Example problem (steady 1-D conduction through composite wall)
Consider two layers with thicknesses L1, L2 and conductivities k1, k2. With temperatures T_left and T_right, steady heat rate:
Q̇ = (T_left − T_right) / (R1 + R2), where R_i = L_i/(k_i A).
Temperature at interface T_int = T_left − Q̇ R1.
This simple model illustrates thermal resistances in series and how material choice and thickness influence heat flow.
13. Measurement and instrumentation
Heat flow and temperature measurement techniques:
- Thermocouples, RTDs, thermistors for temperature.
- Heat flux sensors (heat flux transducers) measure q directly using a known thermal resistance and thermopile.
- Calorimetry (differential scanning calorimetry for material properties).
- Infrared thermography for surface temperature fields and qualitative heat flow visualization.
Material properties (k, c, ρ) vary with temperature; accurate modeling often requires experimental characterization.
14. Applications and engineering relevance
- Building energy efficiency: insulation, thermal bridging, HVAC sizing.
- Electronics cooling: managing heat in chips, PCBs, and power electronics.
- Thermal management in aerospace: re-entry heating, thermal protection systems.
- Energy systems: heat exchangers, boilers, solar thermal collectors.
- Materials processing: heat treatment, welding, additive manufacturing.
Understanding heat flow enables design choices that minimize losses, avoid thermal stresses, and optimize performance.
15. Summary: connecting Fourier, thermodynamics, and energy balances
- Fourier’s law (q = −k ∇T) provides the local relation between temperature gradients and heat flux.
- Conservation of energy plus Fourier’s law yields the heat equation, which governs diffusion of thermal energy.
- Thermodynamic principles (entropy production, irreversible processes) frame heat flow as a driver of irreversibility and link to coupled transport phenomena.
- Practical analysis uses thermal resistances, lumped models, or numerical methods depending on geometry, time scales, and required accuracy.
Understanding the interplay of these elements—constitutive laws, conservation, boundary conditions, and material behavior—lets engineers and scientists predict and control heat flow across scales and applications.
Leave a Reply