Coercivity and energy estimates
Track the constants that control a weak solution, including lower-order terms and approach to resonance.
On this page
The question for today
What you’ll understand
- Derive a gradient stability estimate by testing with the solution.
- Use Poincaré and Young inequalities to handle adverse lower-order terms.
- Distinguish global coercivity from a local cutoff energy estimate.
An energy estimate is both a proof tool and a sensitivity statement. It describes how strongly the equation controls its solution and which constants deteriorate when the problem changes. Merely writing an ellipticity lower bound is not enough if a negative potential cancels the principal energy or if boundary conditions leave constants uncontrolled. We will keep the coefficient bounds, domain inequality, and source norm visible. The resulting estimates help diagnose why a solution can become very large even though a solver reports a small residual.
Spend ellipticity against the source
For the homogeneous Dirichlet equation minus divergence of A gradient u equals F on a bounded Lipschitz domain, assume A is real symmetric, bounded measurable, and uniformly elliptic. Testing with u gives a lower bound lambda times the squared gradient norm. Duality bounds the source pairing by its H-minus-one norm times the gradient norm. After dividing when the gradient is nonzero, the solution norm is bounded by the source norm divided by lambda; the zero-gradient case already satisfies the estimate.
When the source is an L-two function, Cauchy–Schwarz and Poincaré replace its dual norm by the L-two norm times a domain constant. That estimate is useful without knowing the solution explicitly. It also applies to perturbations: the difference of solutions with the same operator and boundary data solves an equation driven by the source difference. A stability constant therefore turns source uncertainty into a solution bound, provided the operator and norm convention are unchanged.
A negative potential can consume the margin
For a potential c bounded below by minus kappa, the form is at least lambda times the gradient energy minus kappa times the L-two energy. If the Poincaré inequality has constant C_P, the remaining gradient coefficient is lambda minus kappa C_P squared. A strictly positive remainder proves coercivity. This is a sufficient condition based on the chosen bounds; its failure alone does not prove that the boundary problem is unsolvable. It means that this particular coercivity argument has stopped providing control.
Bounded first-order terms can be treated with Cauchy–Schwarz and Young's inequality, paying a small part of the gradient energy and leaving an L-two remainder. This produces a Gårding-type lower bound: the form plus a sufficiently large multiple of the L-two norm squared is coercive. The shift changes the operator and gives an invertible reference problem; it does not assert that the unshifted operator has no kernel. Fredholm theory later identifies the remaining finite-dimensional compatibility obstruction.
Local energy needs a cutoff and a larger region
An interior estimate uses a compactly supported cutoff eta and tests a homogeneous equation with eta squared times u. The resulting cross term is controlled by Young's inequality, leaving the gradient energy on the smaller region bounded by the squared solution times the cutoff gradient on the surrounding region. This Caccioppoli mechanism needs only ellipticity and bounded coefficients for scalar divergence-form weak solutions. It avoids prescribing boundary data because the test function vanishes before reaching the outer boundary.
A typical concentric-ball estimate bounds the gradient L-two norm on radius r by a constant times the solution L-two norm on radius R divided by R−r, for r<R and the larger closed ball inside the domain. The constant depends on dimension and ellipticity bounds. The gap between the balls is essential: localization has a derivative cost. Such a local estimate is different from global uniqueness, and it cannot remove a constant-mode ambiguity in a pure Neumann problem without a normalization.
Work it through
Watch the energy margin approach zero
On (0,π), solve with zero endpoint values, for ε>0. Relate the answer to coercivity, and evaluate ε=1/4.
The first Dirichlet sine mode satisfies −(sin x)″=sin x. Substituting u=A sin x makes the left side equal εA sin x, so A=1/ε. The Poincaré constant on (0,π) is one.
For 0<ε≤1, the form is , at least ε times the squared gradient norm by Poincaré. On the first sine mode the inequality is exact, so the coercivity margin genuinely shrinks to zero.
At ε=1/4, u=4 sin x and its gradient norm is . The forcing has not changed, yet the response amplitude is four. At ε=0, testing against sin x would give zero on the operator side and π/2 on the source side, so this particular resonant forcing is incompatible.
The large response near resonance is predicted by the shrinking coercivity margin. Exactly at resonance, compatibility must be checked rather than dividing by zero.
The living notebook
Explore the model
Interactive lab · heat
Use an energy mode with a known decay law
Pause, Reset, and set Fourier mode n=1 and diffusivity κ=0.5. Predict that at time one the L² readout is √π·e⁻⁰·⁵. Repeat from Reset with both n=2 and κ=0.5; the time-one norm is √π·e⁻². This periodic analytic heat mode illustrates a gradient-energy scale through its decay; it is not a stationary Dirichlet solver or a simulation of the negative potential above.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forMode number changes the gradient-energy scale quadratically, but domain and boundary conditions determine which modes an elliptic estimate actually controls.
Pause & reason
Check your understanding
Your turn
Try it yourself
On (0,1), take u=x(1−x), which solves −u″=2 with zero endpoint values. Compute ∫|u′|² and ∫u², then check the sharp interval Poincaré bound ∫u²≤π⁻²∫|u′|².
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Expand the two polynomials and integrate term by term. The squared Poincaré constant is 1/π², not 1/π.
Compare with the solution
Since u′=1−2x, its squared integral is 1/3. Also . The bound reads , equivalent to π²≤10, which is true. The energy identity gives as well. The inequality is close but not exact because this polynomial is not the first Dirichlet sine eigenfunction.
Keep the idea
Take it with you
- Energy estimates quantify sensitivity through explicit coercivity constants.
- Adverse lower-order terms can consume energy control without changing principal ellipticity.
- Interior cutoff estimates pay for localization and do not replace global boundary conditions.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.