Hyperbolic PDE Theory (Wave Propagation and Characteristics)

Core Idea

Hyperbolic PDEs describe wave propagation phenomena with finite speed of information travel. The wave equation u_tt = c²Δu is the prototype, with solutions determined by initial data through D'Alembert's formula (1D) and Kirchhoff's formula (3D). Key features include: finite domain of dependence (the solution at a point depends only on data within the backward characteristic cone), preservation of singularities (unlike parabolic equations, waves carry discontinuities), and energy conservation. The theory of symmetric hyperbolic systems extends these results to Maxwell's equations, elasticity, and general relativity.

Explainer

Hyperbolic PDEs are the mathematical description of wave phenomena—sound, light, elastic vibrations, gravitational waves. Their defining characteristic is finite propagation speed: information travels along characteristic surfaces at a definite speed, creating domains of dependence and influence. This is physically natural (nothing travels faster than light) but mathematically distinctive: it means hyperbolic equations behave very differently from elliptic and parabolic ones.

For the wave equation u_tt = c²Δu in ℝⁿ, the solution with initial data u(x,0) = f(x), u_t(x,0) = g(x) is given by explicit formulas that depend on the dimension. In 1D, D'Alembert's formula u = ½[f(x-ct) + f(x+ct)] + (1/2c)∫g(s)ds shows the solution is a superposition of left- and right-traveling waves. In 3D, Kirchhoff's formula involves an average of the data over a sphere of radius ct—this is Huygens' principle, and it implies a sharp wavefront with no residual signal. In 2D, Hadamard's method of descent from 3D gives a formula involving an integral over a disk, producing wave tails.

The energy method is the primary tool for establishing well-posedness of hyperbolic problems. For the wave equation, E(t) = ½∫[u_t² + c²|∇u|²]dx is conserved, which gives uniqueness (if initial data is zero, energy is zero, so u ≡ 0) and continuous dependence. For general symmetric hyperbolic systems, the energy E = ½∫u^T A₀ u dx satisfies dE/dt ≤ CE for a constant depending on the coefficients, and Gronwall's inequality gives exponential bounds. These estimates are sharp: hyperbolic equations conserve regularity, neither smoothing initial data nor creating new singularities.

The theory of singularities for hyperbolic equations is much richer than for elliptic or parabolic equations. Singularities in the initial data propagate along characteristics: a discontinuity in the initial data traces out a characteristic surface in space-time, creating a wavefront. Diffraction, reflection, and focusing of singularities at boundaries and caustics produce complex patterns described by microlocal analysis and geometric optics. For nonlinear hyperbolic equations, new singularities can form spontaneously—shock waves in gas dynamics, caustics in nonlinear optics—even from smooth initial data. The global existence theory for nonlinear hyperbolic equations remains one of the great open challenges in PDE theory.