Separation of Variables and the Heat Equation

30 min
0/4 practice checks

The heat equation ut=kuxxu_t = k\,u_{xx} models temperature u(x,t)u(x,t) diffusing along a rod (k>0k > 0 is the thermal diffusivity). Separation of variables looks for solutions of the special product form

u(x,t)=X(x)T(t).u(x, t) = X(x)\,T(t).

Substituting and dividing by kXTk\,X\,T gives

TkT=XX.\frac{T'}{k\,T} = \frac{X''}{X}.

The left side depends only on tt, the right only on xx; the only way they can be equal for all x,tx, t is if both equal the same constant λ-\lambda. This splits one PDE into two ordinary ones.

Worked example. From the split, T=λkTT' = -\lambda k\,T gives a decaying time factor T(t)=eλktT(t) = e^{-\lambda k t}, while X=λXX'' = -\lambda X gives oscillating X(x)=sin(λx)X(x) = \sin(\sqrt\lambda\, x) (with suitable boundary conditions). So each mode looks like u=eλktsin(λx)u = e^{-\lambda k t}\sin(\sqrt\lambda\, x). For k=4k = 4 and spatial part sin(2x)\sin(2x) (so λ=2\sqrt\lambda = 2, λ=4\lambda = 4), the decay rate is λk=4×4=16\lambda k = 4\times 4 = 16: u=e16tsin(2x)u = e^{-16t}\sin(2x). As tt \to \infty, every such mode decays to 00.