Energy methods establish properties of PDE solutions by studying integral quantities (energies) rather than the solutions themselves. For the wave equation, the total energy E(t) = ½∫[u_t² + c²|∇u|²]dx is conserved, which immediately proves uniqueness and continuous dependence on data. For the heat equation, the L² norm ∫u²dx is non-increasing, reflecting dissipation. By choosing appropriate energy functionals and showing they satisfy differential inequalities (Gronwall-type arguments), energy methods prove existence, uniqueness, and stability results for a wide variety of linear and nonlinear PDEs without requiring explicit solution formulas.
Energy methods are among the most versatile and powerful techniques in PDE theory, applicable to problems where explicit solution formulas are unavailable. The basic strategy is: multiply the PDE by the solution u (or by u_t, or some other carefully chosen multiplier), integrate over the domain, use integration by parts to move derivatives around, and arrive at an identity or inequality for an energy-like integral. This integral captures the essential dynamics of the system without requiring detailed knowledge of the solution.
For the wave equation u_tt = c²Δu with homogeneous boundary conditions, multiplying by u_t and integrating yields the energy conservation law dE/dt = 0, where E = ½∫[u_t² + c²|∇u|²]dx is the sum of kinetic and potential energy. Conservation of energy immediately implies uniqueness: if two solutions share the same initial data, their difference has zero energy at t = 0 and hence zero energy for all time, which forces the difference to be identically zero. It also implies continuous dependence: if the initial data are close in the energy norm, the solutions remain close for all time.
For dissipative equations like the heat equation, energy methods yield decay rather than conservation. The L² norm ∫u²dx satisfies d/dt ∫u²dx = -2k∫|∇u|²dx ≤ 0, showing that the total "thermal energy" decreases monotonically. Using the Poincaré inequality ∫|∇u|²dx ≥ λ₁∫u²dx (where λ₁ is the first eigenvalue of -Δ), one obtains exponential decay: ∫u²dx ≤ e^(-2kλ₁t)∫u₀²dx. This quantitative decay rate depends on the geometry of the domain through λ₁.
The power of energy methods extends far beyond linear constant-coefficient equations. For nonlinear equations, carefully chosen energy functionals can establish global existence (the solution does not blow up), stability of equilibria, and asymptotic behavior. For systems, energy methods handle coupling between components naturally. In the variational approach to existence theory, the solution is found as a minimizer of an energy functional, and the energy method provides the a priori estimates needed to pass to limits. Energy methods also underlie the stability analysis of numerical schemes: a numerical method is stable if it preserves a discrete analogue of the continuous energy estimate.