Table of Contents

Heat Equation Forms

1D: ∂u/∂t = α ∂²u/∂x²
2D: ∂u/∂t = α (∂²u/∂x² + ∂²u/∂y²)
3D: ∂u/∂t = α ∇²u
Where α = k/(ρc)

Introduction to the Heat Equation

The heat equation is one of the most important partial differential equations in mathematics, physics, and engineering. It describes how heat (or temperature) diffuses through a given region over time, making it fundamental to understanding thermal processes in everything from engineering design to climate science.

Historical Significance:

  • First studied by Joseph Fourier in the early 19th century
  • Pioneered the mathematical study of heat conduction
  • Led to the development of Fourier series and transforms
  • Foundation for modern partial differential equation theory
  • Applications span physics, engineering, finance, and image processing

In this comprehensive guide, we'll explore the heat equation from its mathematical foundations to practical applications, with interactive tools to help you visualize and understand heat diffusion processes.

What is the Heat Equation?

The heat equation is a parabolic partial differential equation that describes the distribution of heat (or variation in temperature) in a given region over time. It's based on Fourier's law of heat conduction and the principle of conservation of energy.

∂u/∂t = α ∇²u

Where:

  • u(x,t) is the temperature at position x and time t
  • α is the thermal diffusivity (α = k/(ρc))
  • k is the thermal conductivity
  • ρ is the density
  • c is the specific heat capacity
  • ∇² is the Laplace operator (second spatial derivatives)

Common Forms:

1D Heat Equation: ∂u/∂t = α ∂²u/∂x²

2D Heat Equation: ∂u/∂t = α (∂²u/∂x² + ∂²u/∂y²)

3D Heat Equation: ∂u/∂t = α (∂²u/∂x² + ∂²u/∂y² + ∂²u/∂z²)

Key Properties
  • Parabolic: Solutions smooth out over time (diffusion)
  • Linear: Principle of superposition applies
  • Maximum Principle: Maximum temperature decreases over time
  • Smoothing: Irregular initial conditions become smooth
  • Energy Conservation: Total heat energy is conserved

Track your progress by practicing with the PDE calculator.

Mathematical Derivation

The heat equation can be derived from basic physical principles: Fourier's law of heat conduction and conservation of energy.

1️⃣

Step 1: Fourier's Law

Heat Flux: q = -k ∇u

Heat flows from hot to cold regions

Rate proportional to temperature gradient

k is thermal conductivity (material property)

2️⃣

Step 2: Energy Conservation

Rate of Change: ∂/∂t ∫ ρc u dV = -∫ q·n dA

Heat change in volume = Net heat flux through surface

ρ is density, c is specific heat capacity

3️⃣

Step 3: Divergence Theorem

Convert Surface to Volume: ∫ q·n dA = ∫ ∇·q dV

Apply divergence theorem

Surface integral becomes volume integral

∇·q = -k ∇²u

4️⃣

Step 4: Final Equation

Combine Results: ρc ∂u/∂t = k ∇²u

Simplify: ∂u/∂t = α ∇²u

Where α = k/(ρc) is thermal diffusivity

This is the heat equation!

Thermal Properties Calculator

Enter material properties and click "Calculate"

Default values are for aluminum

Analytical Solutions

The heat equation can be solved analytically for simple geometries and boundary conditions using various mathematical techniques:

📈

Separation of Variables

Assume: u(x,t) = X(x)T(t)

Leads to ordinary differential equations

X''/X = T'/(αT) = -λ (separation constant)

Solve eigenvalue problem for X(x)

∞

Fourier Series

Solution: u(x,t) = Σ Aₙ sin(nπx/L) e^{-α(nπ/L)²t}

Aₙ determined from initial condition

Works for finite domains with boundary conditions

Each term decays exponentially

∫

Fourier Transform

Infinite Domain: Use Fourier transform

û(k,t) = ∫ u(x,t) e^{-ikx} dx

Transformed equation: ∂û/∂t = -αk² û

Solution: û(k,t) = û(k,0) e^{-αk²t}

🔍

Fundamental Solution

Heat Kernel: K(x,t) = (4παt)^{-1/2} e^{-x²/(4αt)}

Solution for delta function initial condition

General solution: u(x,t) = ∫ K(x-y,t) u₀(y) dy

Shows smoothing property

Example: Rod with Fixed Ends

Consider a rod of length L with ends held at 0°C:

∂u/∂t = α ∂²u/∂x², 0 < x < L, t > 0
u(0,t) = u(L,t) = 0
u(x,0) = f(x)

Solution using separation of variables:

u(x,t) = Σ_{n=1}∞ Aₙ sin(nπx/L) e^{-α(nπ/L)²t}
Aₙ = (2/L) ∫₀ᴸ f(x) sin(nπx/L) dx

Each mode decays exponentially, with higher frequencies decaying faster.

Engage in hands-on learning and sharpen your skills with the PDE calculator.

Numerical Methods

For complex geometries and boundary conditions, numerical methods are essential for solving the heat equation:

🔢

Finite Difference Method

Discretize: Replace derivatives with finite differences

∂u/∂t ≈ (uᵢʲ⁺¹ - uᵢʲ)/Δt

∂²u/∂x² ≈ (uᵢ₊₁ʲ - 2uᵢʲ + uᵢ₋₁ʲ)/Δx²

Explicit and implicit schemes available

🔺

Finite Element Method

Weak Form: Multiply by test function and integrate

Works for complex geometries

Adaptive mesh refinement possible

Industry standard for engineering problems

⚡

Explicit Scheme

Forward Euler: uᵢʲ⁺¹ = uᵢʲ + r(uᵢ₊₁ʲ - 2uᵢʲ + uᵢ₋₁ʲ)

r = αΔt/Δx² must be ≤ 0.5 for stability

Simple to implement

Conditionally stable

🔄

Implicit Scheme

Crank-Nicolson: Unconditionally stable

Requires solving linear system each time step

More accurate than explicit methods

Standard choice for many applications

Finite Difference Calculator

20
100

Real-World Applications

The heat equation has numerous applications across science and engineering:

🏗️

Building Design

Thermal Insulation: Optimize wall thickness and materials

HVAC Systems: Design efficient heating and cooling

Energy Efficiency: Calculate heat loss through buildings

Passive Solar: Design buildings for natural heating

🚗

Automotive Engineering

Engine Cooling: Design radiator and cooling systems

Brake Systems: Analyze heat dissipation in brakes

Exhaust Systems: Manage heat in exhaust components

Battery Thermal Management: Critical for electric vehicles

💻

Electronics Cooling

CPU/GPU Cooling: Design heat sinks and fans

Circuit Board Design: Prevent overheating components

Thermal Management: Essential for high-performance computing

Heat Spreader Design: Distribute heat evenly

🌡️

Medical Applications

Hyperthermia Treatment: Heat cancer cells

Cryotherapy: Freeze tissue for treatment

Thermal Imaging: Analyze heat patterns in body

Drug Delivery: Temperature-sensitive release systems

Material Thermal Properties

Different materials have vastly different thermal properties:

Material k (W/m·K) ρ (kg/m³) c (J/kg·K) α (m²/s)
Copper 401 8960 385 1.16×10⁻⁴
Aluminum 237 2700 900 9.75×10⁻⁵
Steel 50 7850 460 1.38×10⁻⁵
Glass 1.05 2500 840 5.00×10⁻⁷
Wood 0.15 700 2300 9.32×10⁻⁸
Air 0.026 1.2 1005 2.16×10⁻⁵

Engage in hands-on learning and sharpen your skills with the PDE calculator.

Interactive Heat Simulation

2D Heat Equation Simulation

Visualize heat diffusion in a 2D plate with interactive controls.

Simulation Parameters:

0.01
5

Boundary Conditions

Boundary conditions are essential for solving the heat equation. Different conditions lead to different physical situations:

Dirichlet Condition

u = g on boundary

Temperature specified at boundary

Example: Rod ends held at fixed temperature

Neumann Condition

∂u/∂n = h on boundary

Heat flux specified at boundary

Example: Insulated boundary (h=0)

Robin Condition

a u + b ∂u/∂n = c

Mixed condition (convection)

Example: Newton's law of cooling

Periodic Condition

u(0,t) = u(L,t)

Solution repeats periodically

Useful for infinite periodic structures

Example: Insulated Rod

Consider a rod with insulated ends (no heat flux):

∂u/∂t = α ∂²u/∂x², 0 < x < L, t > 0
∂u/∂x(0,t) = ∂u/∂x(L,t) = 0 (Neumann conditions)
u(x,0) = f(x)

Solution:

u(x,t) = A₀/2 + Σ_{n=1}∞ Aₙ cos(nπx/L) e^{-α(nπ/L)²t}
Aₙ = (2/L) ∫₀ᴸ f(x) cos(nπx/L) dx

Note: As t→∞, temperature becomes uniform: u(x,∞) = A₀/2

Turn theory into practice with real-world problems using the PDE calculator.

Advanced Topics

Beyond the basic heat equation, several advanced topics extend its applicability:

Nonlinear Heat Equation

When thermal properties depend on temperature:

ρc(u) ∂u/∂t = ∇·(k(u) ∇u)

Requires numerical solution

Important for phase change problems

Anisotropic Materials

Heat conduction varies with direction:

ρc ∂u/∂t = ∇·(K ∇u)

K is conductivity tensor

Important for composite materials

Heat Equation with Sources

Includes heat generation/absorption:

∂u/∂t = α ∇²u + f(x,t)

f(x,t) represents heat sources

Example: electrical heating

Fractional Heat Equation

Uses fractional derivatives:

∂ᵅu/∂tᵅ = α ∇²u

Models anomalous diffusion

Applications in complex media

Measure your understanding of partial differential equations by using the PDE calculator.

Practice Problems

Problem 1: A copper rod of length 1m has its ends maintained at 0°C. The initial temperature distribution is u(x,0) = 100sin(πx). Find the temperature distribution u(x,t) for t > 0. (α for copper = 1.16×10⁻⁴ m²/s)

Solution:

1. The problem has Dirichlet boundary conditions: u(0,t) = u(1,t) = 0

2. Using separation of variables, the solution is:

u(x,t) = Σ_{n=1}∞ Aₙ sin(nπx) e^{-α(nπ)²t}

3. The initial condition is already a sine function: u(x,0) = 100sin(πx)

4. Comparing with the series, we see that only n=1 term is nonzero:

A₁ = 100, Aₙ = 0 for n ≠ 1

5. Therefore, the solution is:

u(x,t) = 100 sin(πx) e^{-απ²t}

6. Substitute α = 1.16×10⁻⁴:

u(x,t) = 100 sin(πx) e^{-1.16×10⁻⁴ π² t}
Problem 2: A steel plate (α = 1.38×10⁻⁵ m²/s) of thickness 0.1m is initially at 20°C. Suddenly, both surfaces are raised to 100°C. How long will it take for the center to reach 50°C?

Solution:

1. This is a symmetric heating problem with Dirichlet conditions.

2. The solution for a plate with sudden surface temperature change is:

u(x,t) = T_s + (T_i - T_s) Σ_{n=0}∞ (-1)ⁿ erfc((2n+1)L-x)/(2√(αt)) + erfc((2n+1)L+x)/(2√(αt))

3. For the center (x = L/2 = 0.05m), we can use the first term approximation:

(T - T_s)/(T_i - T_s) ≈ erfc(L/(2√(αt)))

4. Plug in values: T = 50°C, T_s = 100°C, T_i = 20°C

(50-100)/(20-100) = 50/80 = 0.625 = erfc(0.05/(2√(1.38×10⁻⁵ t)))

5. Solve: erfc(z) = 0.625 ⇒ z ≈ 0.31 (from erfc table)

0.05/(2√(1.38×10⁻⁵ t)) = 0.31

6. Solve for t: t ≈ 184 seconds ≈ 3 minutes

Continue Your Mathematical Journey

Explore more advanced concepts in partial differential equations and their applications:

What are Partial Differential Equations?

Understand the fundamentals of PDEs and how they model complex systems in physics and engineering.

Heat Equation Explained

Learn how heat flow and diffusion are described using one of the most important PDEs.

Wave Equation Applications

Explore how PDEs model waves, vibrations, and signals across different fields.

Numerical Methods for PDEs

Discover computational approaches used to approximate solutions to complex PDEs.