advanced32 min studyLesson 393 in Mathematics

Coercivity and energy estimates

Track the constants that control a weak solution, including lower-order terms and approach to resonance.

On this page

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.

λ∥∇u∥22≤F(u)≤∥F∥H−1∥∇u∥2,∥∇u∥2≤λ−1∥F∥H−1\lambda\|\nabla u\|_2^2\leq F(u)\leq\|F\|_{H^{-1}}\|\nabla u\|_2,\qquad \|\nabla u\|_2\leq\lambda^{-1}\|F\|_{H^{-1}}
H⁻¹ is the dual of H₀¹ equipped with gradient norm. The inequality assumes homogeneous Dirichlet data and no adverse lower-order contribution.↔ Scroll wide equations horizontally.

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.

a(v,v)≥(λ−κCP2)∥∇v∥22if c≥−κ,∥v∥2≤CP∥∇v∥2a(v,v)\geq\bigl(\lambda-\kappa C_P^2\bigr)\|\nabla v\|_2^2\quad\text{if }c\geq-\kappa,\qquad \|v\|_2\leq C_P\|\nabla v\|_2
A positive displayed coefficient gives coercivity in the gradient norm. A zero or negative coefficient leaves this sufficient test inconclusive; it is not itself a solvability classification.↔ Scroll wide equations horizontally.

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.

Watch the energy margin approach zero

On (0,π), solve −u′′−(1−ε)u=sin⁡x-u''-(1-\varepsilon)u=\sin x with zero endpoint values, for ε>0. Relate the answer to coercivity, and evaluate ε=1/4.

  1. 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.

  2. For 0<ε≤1, the form is ∫∣v′∣2−(1−ε)∣v∣2\int |v'|^2-(1-\varepsilon)|v|^2, 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.

  3. At ε=1/4, u=4 sin x and its gradient norm is 4π/24\sqrt{\pi/2}. 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 result

The large response near resonance is predicted by the shrinking coercivity margin. Exactly at resonance, compatibility must be checked rather than dividing by zero.

Explore the model

Interactive lab

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.

Use an energy mode with a known decay lawThe 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. Current readings: Current amplitude: 1; Spatial L² norm: 1.772; Decay rate · κn²: 3.6 /s. All parameters are adjustable with the labeled controls below.Periodic heat equation mode 3 at time 0: the amplitude is 1.0π/2π3π/22π−1−0.500.51x · periodic intervalTemperature perturbation uExact mode · n = 3
Current temperature modeInitial sineCurrent amplitude
Current amplitude
1
Spatial L² norm
1.772
Decay rate · κn²
3.6 /s
t = 0 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.

Check your understanding

A coercivity estimate gives λ−κC_P²=0. What follows from this estimate alone?

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 ∫01x2(1−x)2dx=1/30\int_0^1x^2(1-x)^2dx=1/30. The bound reads 1/30≤1/(3π2)1/30\leq1/(3\pi^2), equivalent to π²≤10, which is true. The energy identity gives ∫2u=1/3\int 2u=1/3 as well. The inequality is close but not exact because this polynomial is not the first Dirichlet sine eigenfunction.

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.