Principia Orthogona  ·  Book 6  ·  Differential Equations
Elliptic / Poisson Foundations Heat Equation Lifting to a Box Domain The Wave Equation Helix Toy Model The Cosmic No-Go A Nonlinear Reaction-Diffusion Fold
d/dt ⌁½(ut²+ux²)dx = 0
energy conserved, not decayed

The Wave Equation

the hyperbolic leg: a different toolkit, and the exact mirror of the heat equation's own H-theorem

The first three chapters of this corpus share one toolkit: a static bilinear form, coercivity, Lax–Milgram. That toolkit is elliptic machinery, and it does not extend to this chapter. The wave equation is genuinely second-order in time — it needs both an initial position and an initial velocity to be well-posed — and its natural well-posedness argument is the energy method: multiply by the time derivative, integrate by parts, and track a conserved quantity directly, rather than inverting a stationary operator. Completing the elliptic/parabolic/hyperbolic triad this deliberately is also what makes the comparison to the heat equation's own Part V exact rather than rhetorical, which is the real payoff of this chapter.

The model

Take the 1D wave equation on (0,1) with homogeneous Dirichlet boundary conditions and both initial position and initial velocity data:

utt = uxx,   u(0,t)=u(1,t)=0,   u(x,0)=u₀(x),  ut(x,0)=u₁(x).

Unlike the heat equation, there is no mass term and no dissipation anywhere in this equation — nothing to make a stationary bilinear form coercive, because there is no stationary bilinear form to begin with. What replaces it is a quantity that is exactly conserved along trajectories.

The energy method

Multiply the equation by ut and integrate over (0,1):

∫uttutdx = ∫uxxutdx.

The left side is ½ d/dt ∫ut²dx. Integrating the right side by parts gives a boundary term [uxut]₀¹, which vanishes exactly: the Dirichlet condition u(0,t)=u(1,t)=0 for all t forces ut(0,t)=ut(1,t)=0 as well, so both endpoints kill the term regardless of what ux does there — the same vanishing-boundary-term step every earlier chapter in this corpus has used, applied here to the time derivative of the boundary condition rather than to a test function. What remains is −½ d/dt ∫ux²dx. Together:

d/dt ⌃E(t) = 0,    E(t) := ½∫₀¹(ut²+ux²)dx.

E(t) is exactly conserved, not merely bounded. This single identity is the whole well-posedness argument for this class of problems in outline: a Galerkin approximation (project onto finitely many Dirichlet eigenmodes, solve the resulting finite-dimensional ODE system exactly, and pass to the limit) inherits the same energy identity at every finite level, which supplies the a priori bound needed to extract a convergent subsequence. No coercive bilinear form and no Lax–Milgram appear anywhere in that argument — energy conservation does the entire job elliptic coercivity did in the earlier chapters.

The exact mirror of the heat equation's H-theorem

The heat equation chapter's Part V (Prop. 20.1) proves H(t) := ½∫u²dx is strictly decreasing along solutions — the same multiply-and-integrate-by-parts technique, applied to u itself rather than ut, because the heat equation is only first-order in time. That decay is irreversible: Corollary 20.3 there notes the forward solution map loses information, since two different initial states can decay toward indistinguishable states as t → ∞, and nothing recovers the initial data from a late-time snapshot. This chapter's E(t) is the structurally identical computation with the opposite outcome: because the equation is second-order in time, the boundary term that killed the mass term's contribution here instead preserves both the kinetic (ut²) and potential (ux²) pieces exactly, and the solution map (u₀,u₁) ↦ (u(·,T),ut(·,T)) is invertible — run time backward from any snapshot and the original data returns exactly, because u(x,−t) solves the same equation whenever u(x,t) does. Heat loses information; waves do not. Same corpus, same integration-by-parts step, opposite conclusion — which is exactly what a second-order-in-time versus first-order-in-time distinction should produce, not a coincidence needing further explanation.

Verified this session (scipy.integrate.quad, exact standing-wave solution u(x,t) = cos(pi t) sin(pi x) on (0,1), Dirichlet BC, u0=sin(pi x), u1=0): PDE check u_tt = u_xx at 4 sample points: exact agreement (diff = 0.0) Boundary check u(0,t)=u(1,t)=0 at t=1.23: confirmed (~1e-17) Energy E(t) = (1/2) int (u_t^2 + u_x^2) dx, sampled at t = 0, 0.3, 0.77, 1.5, pi, 10.0: every value = 2.4674011003 (= pi^2/4, exact, to 10 digits) Conserved across more than 3 full periods -- not just at t=0.

Finite propagation speed

The other qualitative break from the heat equation is how information travels. The heat equation's Green's function is a Gaussian supported on all of for any t>0 — a localized disturbance is felt everywhere instantly, however faintly. The wave equation on the whole line admits the classical d'Alembert solution

u(x,t) = ½[u₀(x−t)+u₀(x+t)] + ½∫x−tx+tu₁(s)ds,

which shows the value at (x,t) depends only on initial data within the interval [x−t,x+t] — a finite domain of dependence, propagation speed exactly 1 in these units. A compactly supported disturbance stays compactly supported, expanding at a fixed rate rather than spreading into a tail everywhere at once.

Verified this session (numpy, d'Alembert's formula, u0 a smooth bump on [0.4,0.6], u1=0, sampled on a fine grid on the line): t=0.0: support = [0.4000, 0.6000] predicted [0.4000, 0.6000] t=0.5: support = [-0.1000, 1.1000] predicted [-0.1000, 1.1000] t=1.0: support = [-0.6000, 1.6000] predicted [-0.6000, 1.6000] t=2.0: support = [-1.6000, 2.6000] predicted [-1.6000, 2.6000] Point x=2.0 at t=1.0 (strictly outside the light cone [-0.6,1.6]): u = 0.0 exactly -- not yet reached, not a small residual.
StatementStatusNotes
Energy identity d/dt E(t) = 0derived + numerically checkedVerified against the exact standing-wave solution over 3+ periods, not just at t=0
Time-reversibility / invertible solution mapderivedFollows from u(x,-t) solving the same equation; contrasted directly against the heat equation's Cor. 20.3
Finite propagation speed (d'Alembert)derived + numerically checkedExact light-cone support confirmed on a fine grid, whole-line case
Galerkin / energy-method well-posedness on (0,1)sorryOutlined above; the finite-dimensional ODE step and the passage to the limit are not yet formalized
Lean / AXLE formalizationnot startedNo scaffold file yet; would need a genuinely new object (u_t as an independent field) rather than reusing H10 alone

What this chapter shows: the energy identity and the finite-propagation result, both derived honestly and both checked against an exact solution or exact formula, not merely asserted. What it does not show: a formalized existence-uniqueness theorem on the bounded domain (0,1) — the Galerkin argument is described, not carried out in Lean, and no scaffold file exists yet for this chapter, unlike the heat equation's HeatEquation_Step1.lean. The energy method outlined here is the standard textbook route (Evans, Partial Differential Equations, §7.2), not a novel result; the contribution of this chapter is placing it inside this corpus's own honest-ledger convention and making its contrast with the heat equation's H-theorem explicit rather than left implicit.

← Lifting to a Box Domain Book 6 index →