Heat Conduction Simulation: Fourier's Law & Thermal Diffusivity

simulator intermediate ~8 min
Loading simulation...
Steel bar: α ≈ 1.28 × 10⁻⁵ m²/s

A steel bar (k=50, ρ=7800, cₚ=500) has thermal diffusivity α ≈ 1.28×10⁻⁵ m²/s. Heat from the 500°C end gradually diffuses toward the cold end following the diffusion equation.

Formula

q = -k dT/dx (Fourier's law)
α = k / (ρ cₚ) (thermal diffusivity)
∂T/∂t = α ∂²T/∂x² (1D heat equation)

Heat Flow Through Solids

When one end of a metal bar is heated, the temperature change doesn't appear instantly at the other end. Heat conducts through the bar atom by atom — energetic atoms vibrate and transfer kinetic energy to their neighbors. In metals, free electrons carry most of the thermal energy, which is why good electrical conductors are also good thermal conductors. This simulation visualizes the 1D heat equation as temperature evolves in real time.

Fourier's Law

Joseph Fourier established in 1822 that heat flux q (energy per unit area per unit time) is proportional to the temperature gradient: q = -k(dT/dx). The proportionality constant k is thermal conductivity — a material property ranging from ~0.02 W/(m·K) for aerogel to ~400 W/(m·K) for copper. The negative sign means heat flows from hot to cold, satisfying the second law of thermodynamics.

Thermal Diffusivity

While conductivity k tells you how much heat flows, thermal diffusivity α = k/(ρcₚ) tells you how fast temperature changes propagate. A material can have high conductivity but also high heat capacity, making it slow to change temperature. Diffusivity captures this balance. High-α materials like copper reach equilibrium quickly; low-α materials like concrete or brick store heat for hours — which is why masonry buildings stay cool in summer.

Engineering Applications

Thermal conduction analysis is essential in electronics cooling (preventing chip overheating), building insulation design, metallurgical heat treatment (quenching rate determines microstructure), and spacecraft thermal management. The heat equation, solved numerically with finite differences or finite elements, is one of the most computed PDEs in engineering practice.

FAQ

What is thermal diffusivity?

Thermal diffusivity α = k/(ρcₚ) measures how quickly temperature changes propagate through a material. High α means the material reaches thermal equilibrium fast. Copper (α ≈ 1.1×10⁻⁴) equilibrates ~8× faster than steel (α ≈ 1.3×10⁻⁵).

What is Fourier's law of heat conduction?

Fourier's law states that heat flux is proportional to the negative temperature gradient: q = -k(dT/dx). The thermal conductivity k determines how much heat flows per unit area per unit temperature gradient. It's the conduction equivalent of Ohm's law.

Why does copper feel colder than wood at room temperature?

Copper has much higher thermal conductivity (~400 W/m·K vs ~0.2 W/m·K for wood). When you touch copper, heat flows rapidly from your hand into the material, making it feel cold. Wood conducts heat so slowly that your skin stays warm.

What is the heat equation?

The heat equation ∂T/∂t = α ∂²T/∂x² is a partial differential equation that describes how temperature evolves in space and time. The thermal diffusivity α determines the rate. It was first studied by Joseph Fourier in 1822.

Sources

Embed

<iframe src="https://homo-deus.com/lab/materials-science/heat-conduction/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub