MathLabs
TheoremProved

Existence of the solution by separation of variables

Statement

For the heat equation ut=αuxxu_t = \alpha u_{xx} on 0<x<L0<x<L with boundary conditions u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 and initial condition u(x,0)=f(x)u(x,0) = f(x) where ff is piecewise continuous, the series u(x,t)=∑n=1∞bnsin⁡ ⁣(nπxL)e−α(nπ/L)2tu(x,t) = \displaystyle\sum_{n=1}^{\infty} b_n \sin\!\left(\frac{n\pi x}{L}\right) e^{-\alpha (n\pi/L)^2 t} with coefficients bn=2L∫0Lf(x)sin⁡ ⁣(nπxL)dxb_n = \dfrac{2}{L}\displaystyle\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx converges (for t>0t>0) to a solution satisfying all three conditions.

Why is it true?

Sine functions with these particular frequencies are exactly the shapes that keep zero temperature at both ends while decaying independently in time; since every reasonable initial shape can be written as a sum of such sine waves (a Fourier sine series), the same decomposition immediately produces the solution at all later times.

Proof sketch

Seek a solution of the separated form u(x,t)=X(x)T(t)u(x,t)=X(x)T(t). Substituting into ut=αuxxu_t = \alpha u_{xx} gives X(x)T′(t)=αX′′(x)T(t)X(x)T'(t) = \alpha X''(x)T(t), and dividing by αX(x)T(t)\alpha X(x)T(t) separates the variables: T′(t)αT(t)=X′′(x)X(x)=−λ\dfrac{T'(t)}{\alpha T(t)} = \dfrac{X''(x)}{X(x)} = -\lambda, where λ\lambda must be a constant because the left side depends only on tt and the right side only on xx.

The boundary conditions u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 force X(0)=X(L)=0X(0)=X(L)=0. The eigenvalue problem X′′+λX=0X'' + \lambda X = 0 with these boundary conditions has nontrivial solutions only for λn=(nπ/L)2\lambda_n = (n\pi/L)^2, n=1,2,3,…n=1,2,3,\dots, with eigenfunctions sin⁡ ⁣(nπxL)\sin\!\left(\frac{n\pi x}{L}\right); any other value of λ\lambda forces X≡0X \equiv 0. Solving T′(t)=−αλnT(t)T'(t) = -\alpha\lambda_n T(t) then gives Tn(t)=e−α(nπ/L)2tT_n(t) = e^{-\alpha (n\pi/L)^2 t}.

Each product un(x,t)=sin⁡ ⁣(nπxL)e−α(nπ/L)2tu_n(x,t) = \sin\!\left(\frac{n\pi x}{L}\right)e^{-\alpha (n\pi/L)^2 t} solves the equation and boundary conditions. By linearity, so does any finite or (under mild convergence conditions) infinite sum u(x,t)=∑n=1∞bnsin⁡ ⁣(nπxL)e−α(nπ/L)2tu(x,t) = \displaystyle\sum_{n=1}^{\infty} b_n \sin\!\left(\frac{n\pi x}{L}\right) e^{-\alpha (n\pi/L)^2 t}. Setting t=0t=0 and using the orthogonality relation ∫0Lsin⁡(nπxL)sin⁡(mπxL) dx=0\int_0^L \sin(\frac{n\pi x}{L})\sin(\frac{m\pi x}{L})\,dx = 0 for n≠mn \ne m (and =L/2=L/2 for n=mn=m) to project f(x)f(x) onto each mode yields exactly the coefficients bn=2L∫0Lf(x)sin⁡ ⁣(nπxL)dxb_n = \dfrac{2}{L}\displaystyle\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx, so the initial condition u(x,0)=f(x)u(x,0) = f(x) is matched termwise.

Topics that use this theorem

Step-by-step proofs

No step-by-step proof yet for this theorem.

References

  1. Lawrence C. Evans (2010). Partial Differential Equations
  2. James Ward Brown, Ruel V. Churchill (2011). Fourier Series and Boundary Value Problems