The Wave Equation

Part 14, Chapter 14: Differential Equations

Learning objectives

  • State the wave equation and identify the wave speed cc
  • Write d'Alembert's solution on the infinite line as two travelling waves plus an integral over initial velocity
  • Solve the finite-string problem by separation of variables and recognize the resulting standing waves
  • Contrast wave propagation (finite speed, energy conservation) with heat diffusion (infinite speed, smoothing)

The wave equation governs everything from the vibrating guitar string to electromagnetic radiation to seismic shock waves. Its solutions split cleanly into travelling waves (on the line) or standing waves (on a finite string), and unlike the heat equation it preserves energy, nothing decays, nothing smooths. The two pillars are d'Alembert's formula (which gives the explicit solution on mathbbR\mathbb{R} as two waves travelling in opposite directions) and separation of variables (which on a finite string produces the normal modes that are the vibrating-string Fourier basis).

The equation

The one-dimensional wave equation:

dfracpartial2upartialt2=c2dfracpartial2upartialx2\dfrac{\partial^2 u}{\partial t^2} = c^2 \dfrac{\partial^2 u}{\partial x^2}

c>0c > 0 is the wave speed, the speed at which disturbances propagate. The PDE is hyperbolic, linear, homogeneous, and second order in both xx and tt. Because tt enters at second order, you need two initial conditions: u(x,0)=f(x)u(x, 0) = f(x) (initial displacement) and ut(x,0)=g(x)u_t(x, 0) = g(x)t(x,0)=g(x) (initial velocity).

D'Alembert's solution on the infinite line

On all of mathbbR\mathbb{R}, the wave equation has a closed-form solution. Substituting xi=xct\xi = x - ct and eta=x+ct\eta = x + ct turns uttc2uxx=0u_{tt} - c^2 u_{xx} = 0xx=0 into uxieta=0u_{\xi \eta} = 0xieta=0, whose general solution is

u(x,t)=F(xct)+G(x+ct)u(x, t) = F(x - ct) + G(x + ct)

for arbitrary functions F,GF, G. The first term F(xct)F(x - ct) is a wave travelling rightward at speed cc; the second is a wave travelling leftward at speed cc. Matching the initial conditions u(x,0)=f(x)u(x, 0) = f(x) and ut(x,0)=g(x)u_t(x, 0) = g(x)t(x,0)=g(x) gives d'Alembert's formula:

u(x, t) = \tfrac{1}{2}\bigl[f(x - ct) + f(x + ct)\bigr] + \dfrac{1}{2c} \displaystyle\int_{x - ct}^{x + ct} g(s) \, ds

If the initial velocity gequiv0g \equiv 0, this reduces to two half-amplitude copies of the initial profile travelling in opposite directions, an exact bilateral split.

The finite-string problem: standing waves

For a string of length LL fixed at both ends, u(0,t)=u(L,t)=0u(0, t) = u(L, t) = 0, separation of variables produces the same spatial eigenfunctions as the heat equation, sin(npix/L)\sin(n\pi x / L), but the time ODE is now second order and oscillatory:

Tn(t)+(cnpi/L)2Tn(t)=0T_n''(t) + (cn\pi/L)^2 T_n(t) = 0n(t)=0

so Tn(t)=Ancos(omegant)+Bnsin(omegant)T_n(t) = A_n \cos(\omega_n t) + B_n \sin(\omega_n t)n(t)=Ancos(omegant)+Bnsin(omegant) with natural frequencies omegan=cnpi/L\omega_n = cn\pi/Ln=cnpi/L. The general solution is

u(x, t) = \displaystyle\sum_{n=1}^{\infty} \bigl[A_n \cos(\omega_n t) + B_n \sin(\omega_n t)\bigr] \sin(n\pi x / L)

Each term sin(npix/L)cos(omegant)\sin(n\pi x / L) \cos(\omega_n t)nt) is a standing wave: a fixed spatial shape multiplied by a time-varying amplitude. The points where sin(npix/L)=0\sin(n\pi x / L) = 0 are nodes, they stay still. The maxima between nodes are antinodes, oscillating at the natural frequency.

Plot the first standing wave u(x,t)=sin(pix/L)cos(cpit/L)u(x, t) = \sin(\pi x / L) \cos(c\pi t / L) for several values of tt. The shape stays sin(pix/L)\sin(\pi x / L); only the overall amplitude oscillates. Now plot the second mode sin(2pix/L)cos(2cpit/L)\sin(2\pi x / L) \cos(2c\pi t / L): notice the new node at x=L/2x = L/2, fixed forever in time.

Energy conservation

The total mechanical energy of the vibrating string is

E(t) = \tfrac{1}{2} \displaystyle\int_0^L \bigl[(u_t)^2 + c^2 (u_x)^2\bigr] \, dx

Differentiating EE in time and using utt=c2uxxu_{tt} = c^2 u_{xx}xx with integration by parts shows dE/dt=0dE/dt = 0: energy is conserved exactly. This is in sharp contrast to the heat equation, where energy (in the sense of L2L^2 norm) decays exponentially to zero. The wave equation is time-reversible; the heat equation is not.

Where this shows up
  • Acoustics and musical instruments: A vibrating guitar string, organ pipe, or drum head satisfies the wave equation. The natural frequencies omegan=cnpi/L\omega_n = cn\pi/Ln=cnpi/L are exactly the harmonics you hear, tuning a guitar string changes cc (tension over linear density), which scales the entire harmonic series.
  • Electromagnetism: Maxwell's equations in vacuum reduce to the 3D wave equation nabla2E=(1/c2)Ett\nabla^2 E = (1/c^2) E_{tt}tt with c=1/sqrtvarepsilon0mu0c = 1/\sqrt{\varepsilon_0 \mu_0}, the speed of light. Every light wave, radio wave, and microwave you have ever experienced satisfies this PDE.
  • Seismic waves: P-waves and S-waves through the Earth satisfy wave equations with different speeds. The difference in arrival times between P and S at a seismograph is what lets seismologists triangulate earthquake epicenters.
  • Quantum field theory: Free particles in QFT (Klein-Gordon, Dirac) satisfy wave equations or their relativistic generalisations. The mass term adds m2um^2 u to the right side, but the propagation-at-finite-speed structure is identical.
  • Pause and think: Where does the factor 1/(2c)1/(2c) in d'Alembert's formula come from? Hint: differentiate u(x,t)=F(xct)+G(x+ct)u(x, t) = F(x - ct) + G(x + ct) in tt and set t=0t = 0 to recover the initial-velocity condition.

    Try it

    • Identify the wave speed in utt=9uxxu_{tt} = 9 u_{xx}xx.
    • A string of length L=1L = 1 has wave speed c=2c = 2. Compute the first three natural frequencies omegan\omega_nn.
    • Use d'Alembert with f(x)=ex2f(x) = e^{-x^2} and gequiv0g \equiv 0, c=1c = 1. Sketch the solution at t=0,t=1,t=2t = 0, t = 1, t = 2.
    • For the nn-th standing-wave mode un(x,t)=sin(npix/L)cos(omegant)u_n(x, t) = \sin(n\pi x / L) \cos(\omega_n t)nt), list the positions of the nodes inside [0, L].
    • Show that the energy E(t)E(t) is conserved: differentiate, substitute utt=c2uxxu_{tt} = c^2 u_{xx}xx, and integrate by parts using the fixed boundary conditions.
    • A trap to watch for

      The wave equation has finite propagation speed cc: an initial disturbance at x_0x_0 can only affect points xx with xx0leqct|x - x_0| \leq ct by time tt. The triangular region (x,t):tgeq0,xx0leqct\{(x, t) : t \geq 0, |x - x_0| \leq ct\} is the future light cone of (x0,0)(x_0, 0); outside this cone uu is unaffected by the initial datum at x0x_0. Students who pattern-match the heat equation's "anything affects everything instantly" intuition onto the wave equation will get the qualitative behavior wrong. Heat: instant smearing. Wave: travelling pulse.

      What you now know

      You can write down d'Alembert's solution to the wave equation on mathbbR\mathbb{R}, solve the finite-string problem as a Fourier series of standing waves with natural frequencies cnpi/Lcn\pi/L, identify nodes and antinodes, and prove that total mechanical energy is exactly conserved. With this chapter complete, you have the basic Fourier-and-PDE toolkit that underlies all of mathematical physics, sufficient to read advanced texts on quantum mechanics, fluid dynamics, signal processing, or numerical methods.

      Mark section complete →

      References

      • Garrity, T. (2002). All the Mathematics You Missed: But Need to Know for Graduate School. Cambridge University Press, ch. 14.
      • Evans, L. C. (2010). Partial Differential Equations (2nd ed.). AMS, ch. 2 (wave equation, d'Alembert).
      • Strauss, W. A. (2007). Partial Differential Equations: An Introduction (2nd ed.). Wiley, ch. 2 (wave equation).
      • John, F. (1991). Partial Differential Equations (4th ed.). Springer, ch. 2-3.
      • Jackson, J. D. (1999). Classical Electrodynamics (3rd ed.). Wiley, ch. 6-7 (electromagnetic waves).

      This page is prerendered for SEO and accessibility. The interactive widgets above hydrate on JavaScript load.