advanced33 min studyLesson 411 in Mathematics

Characteristic surfaces and hyperbolicity

Read propagation directions from principal symbols while distinguishing scalar characteristic roots from a system's well-posedness structure.

On this page

What you’ll understand

  • Test a hypersurface using the principal symbol and its normal covector.
  • Compute directional wave speeds in a constant anisotropic equation.
  • Distinguish real characteristic roots from stronger system assumptions.

A propagating disturbance is organized by surfaces in space and time, not only by a solution curve at one point. A characteristic surface is one on which the highest-derivative part of a PDE fails to determine the corresponding normal derivative from tangential information. Its normal therefore reveals the equation's propagation geometry. This language connects transport rays, wave cones, and first-order systems, while keeping clear which structural assumptions are needed for a stable initial-value problem.

A symbol acts on a normal covector

For the constant-speed wave equation utt−c2Δu=0u_{tt}-c^2\Delta u=0, with c positive, choose principal symbol p(τ,ξ)=τ2−c2∣ξ∣2p(\tau,\xi)=\tau^2-c^2|\xi|^2. A smooth hypersurface ϕ(t,x)=0\phi(t,x)=0 is characteristic where its nonzero normal covector satisfies p(ϕt,∇ϕ)=0p(\phi_t,\nabla\phi)=0. The normal covector contains derivatives of the surface-defining function; it is not a point on the surface or a velocity vector. Multiplying the PDE or the defining function by a nonzero factor does not change this zero test.

A constant-time initial surface is noncharacteristic for this equation, so displacement and velocity can be prescribed there under the usual compatibility and regularity assumptions. In one spatial dimension the characteristic lines have slopes determined by the two traveling directions. A characteristic surface is instead adapted to transporting information along a front, and arbitrary independent Cauchy data cannot generally be prescribed on it in the same way. The noncharacteristic test is structural; global existence still requires a separate analysis of the equation and data.

p(ϕt,∇ϕ)=0,p(τ,ξ)=τ2−c2∣ξ∣2,τ=±c∣ξ∣p(\phi_t,\nabla\phi)=0,\qquad p(\tau,\xi)=\tau^2-c^2|\xi|^2,\qquad \tau=\pm c|\xi|
For every nonzero spatial covector, the wave symbol has two distinct real time roots. The characteristic condition concerns the normal to a space-time surface.↔ Scroll wide equations horizontally.

Directional speeds need not be isotropic

Replacing the same spatial coefficient in every direction by different positive coefficients changes the propagation cone. A plane front is tested by inserting its spatial normal and temporal frequency into the symbol. Dividing the temporal frequency by the Euclidean length of the spatial normal gives the front's normal phase speed. A non-unit spatial normal must therefore be normalized before its coefficient is described as a physical speed; otherwise rescaling the surface equation would appear to change propagation.

In an anisotropic medium, the direction of a plane front's normal need not coincide with the direction in which a localized wave packet's group velocity points. Both arise from the same dispersion relation, but they are different derivatives or ratios of it. The examples here calculate normal front speed only. This precision matters when converting a symbol calculation into a domain-of-dependence shape or comparing it with an experiment that tracks packets rather than an extended plane crest.

utt=4uxx+9uyy,ω(ξ,η)=4ξ2+9η2,cnormal=ωξ2+η2u_{tt}=4u_{xx}+9u_{yy},\qquad\omega(\xi,\eta)=\sqrt{4\xi^2+9\eta^2},\qquad c_{\mathrm{normal}}=\frac{\omega}{\sqrt{\xi^2+\eta^2}}
For a nonzero spatial normal, the normal phase speed lies between two and three. Anisotropy makes its value depend on direction.↔ Scroll wide equations horizontally.

Systems need a controlled characteristic basis

For a first-order constant-coefficient system with time derivative isolated, spatial Fourier transformation gives an ordinary differential equation involving a matrix symbol for each spatial direction. Real eigenvalues describe real characteristic speeds, but a defective matrix can still create frequency-dependent growth. A Jordan block produces a factor proportional to frequency times time in its exponential. High-frequency perturbations then lose the uniform estimate expected of a well-posed evolution in a fixed Sobolev norm.

Strong hyperbolicity supplies a uniformly controlled diagonalization or an equivalent suitable symmetrizer, rather than real eigenvalues alone. Symmetric hyperbolic systems offer a concrete sufficient structure: a uniformly positive time matrix and symmetric spatial matrices yield energy estimates. These distinctions are especially important for repeated characteristic speeds. Strictly separated real speeds in a smooth finite-dimensional constant-coefficient setting are a convenient case, but a scalar classification slogan should not be substituted for a complete system estimate.

det⁡ ⁣(τA0+∑j=1dξjAj)=0\det\!\left(\tau A^0+\sum_{j=1}^{d}\xi_jA^j\right)=0
For A⁰ invertible, this determinant locates system characteristic covectors. An energy or uniform-basis condition is still needed to justify a robust evolution estimate.↔ Scroll wide equations horizontally.

A tilted front in an anisotropic wave equation

For u_tt=4u_xx+9u_yy, determine when the plane φ=x+2y−at is characteristic and find its Euclidean normal propagation speed.

  1. The surface normal covector is (φ_t,φ_x,φ_y)=(−a,1,2). Insert it into the principal symbol to obtain a²−4−9·4=a²−40.

  2. The plane is characteristic precisely when a=±√40=±2√10. The coefficient a alone depends on the chosen normalization of the plane's spatial normal.

  3. The spatial normal has length √5, so the magnitude of the normal front speed is √40/√5=√8≈2.828427. This lies between the axial speeds two and three.

a2=40,∣cnormal∣=405=8a^2=40,\qquad |c_{\mathrm{normal}}|=\frac{\sqrt{40}}{\sqrt5}=\sqrt8
Normalizing the spatial covector removes the arbitrary scale of the surface-defining function.↔ Scroll wide equations horizontally.
The result

The two signs describe opposite normal propagation directions. The result is a front speed, not an assertion about the packet group-velocity direction.

Explore the model

Interactive lab

Interactive lab · waves

Interference changes a profile, not its fixed transport speed

Pause and Reset, then set phase difference 0° and amplitude ratio 1. Predict the sum at x=π/4,t=0, then Play to follow a crest. For a fresh comparison, Reset and set phase 180° and ratio 1. Predict the cancellation. Each component has wave number two and angular frequency one.

Interference changes a profile, not its fixed transport speedTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). 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 forThe first initial value is two; the second sum cancels. Nonzero profiles travel at speed one half in this fixed linear model. Changing phase does not change its principal symbol or simulate an anisotropic wave cone.

Check your understanding

Why can real eigenvalues of a system's directional symbol be insufficient for a uniform high-frequency evolution estimate?

Try it yourself

For u_tt−4u_xx=0 in one dimension, test φ=x−2t and ψ=x−3t as characteristic surfaces. Is the initial surface t=0 characteristic?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Insert each surface's temporal and spatial derivatives into τ²−4ξ².

Compare with the solution

For φ, the symbol is (−2)²−4·1²=0, so its zero set is characteristic. For ψ it is 9−4=5, so it is noncharacteristic. For the initial surface φ₀=t, the normal is (1,0) and the symbol is one, so t=0 is noncharacteristic. This permits the standard two initial wave data, subject to the appropriate function-space or smoothness assumptions.

Take it with you

  • Characteristic tests use the principal symbol evaluated on a nonzero surface normal.
  • Normalize spatial normals before interpreting a front coefficient as a speed.
  • Real characteristic roots alone do not replace a controlled system energy estimate.

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.