Mathematics
Fourier Series for Solving Heat Equations
Quick fact
Joseph Fourier introduced his series while studying heat conduction, and his work was initially rejected by mathematicians because he claimed that even discontinuous functions could be represented by sums of smooth sine and cosine curves—a revolutionary idea that later became a cornerstone of modern analysis.
Why this is interesting
Why does a hot metal rod cool down in that smooth, predictable way? And what do musical tones have to do with it?
Read the full explanation
Understanding Fourier Series for Solving Heat Equations
Imagine you have a metal rod with a temperature distribution that varies along its length. The heat equation describes how that temperature changes over time. Fourier's key insight was to look at the initial temperature pattern as a sum of simpler 'wave-like' patterns—sine and cosine curves. Each curve has a certain 'wavelength' (roughly, how bumpy it is). By breaking the initial state into these simple waves, we can solve the heat equation for each wave individually. This works because heat diffusion is linear: the total effect is just the sum of the effects of each wave. This is the principle of superposition. So, if we can figure out how each wave decays, we can add up the decaying waves to get the temperature at any later time.
A deeper explanation
The heat equation is a partial differential equation (PDE) that relates the rate of change of temperature with respect to time to the curvature of temperature in space. Mathematically, it is ∂u/∂t = α ∂²u/∂x², where u(x, t) is temperature. To solve it, we assume the solution can be separated: u(x, t) = X(x)T(t). Substituting this into the equation yields two ordinary differential equations, linked by a separation constant. For appropriate boundary conditions (e.g., fixed temperature at the ends), the spatial equation has solutions that are sine and cosine functions, the Fourier series. The key property is that these sinusoidal functions are eigenfunctions of the second derivative operator: taking the second derivative yields a constant multiplied by the same function. This reduces the PDE to a simple ordinary differential equation for T(t) that has exponential decay. The decay rate depends on the frequency (wavelength) of the wave: high-frequency waves (short wavelength) decay quickly, while low-frequency waves (long wavelength) persist longer. Thus, the temperature distribution smooths out over time, and the final state is a slow-varying pattern. The Fourier series coefficients are determined by the initial temperature distribution, using orthogonality properties of sine and cosine functions. This method is powerful because it works for many other linear PDEs, such as the wave equation, and it has led to the development of Fourier analysis in signal processing, image compression, and solving differential equations across physics and engineering.