Systems of conservation laws and wave families
Decompose coupled perturbations into characteristic families and state the local hypotheses behind nonlinear Riemann constructions.
On this page
The question for today
What you’ll understand
- Classify characteristic fields by the variation of their speed along an eigenvector.
- Solve a two-family linear acoustic Riemann problem.
- Distinguish local small-wave systems theory from the bounded scalar entropy theorem.
A vector conservation law couples several quantities, so one scalar shock or fan is usually insufficient to connect arbitrary neighboring states. Each eigenvalue of the flux Jacobian identifies a characteristic family, and the associated eigenvectors determine which state changes belong to it. Nonlinear families can steepen into shocks or spread into fans; linearly degenerate families carry different kinds of discontinuities. The local wave decomposition is a powerful construction, but its assumptions should remain visible when moving beyond a scalar equation.
Eigenvectors describe permitted infinitesimal state changes
For a smooth flux , the classical system is . Strict hyperbolicity in a neighborhood means that its Jacobian has m distinct real eigenvalues, ordered increasingly there. The right eigenvector of a family points along its infinitesimal wave curve. Across a finite discontinuity, however, the full vector Rankine–Hugoniot equation is required; replacing a finite jump by a local eigenvector relation is only a small-amplitude approximation.
A field is genuinely nonlinear when its characteristic speed changes nontrivially along its eigenvector direction throughout the neighborhood. A field is linearly degenerate when that directional derivative vanishes identically there. The classification concerns variation in state space, not whether a particular spatial solution happens to look like a straight line. Nonlinear gas models can contain both types. Loss of strict hyperbolicity or a field changing type introduces complications that the elementary local construction does not automatically resolve.
A local Riemann solver orders several families
For a sufficiently smooth strictly hyperbolic one-dimensional flux, with each field genuinely nonlinear or linearly degenerate near a reference state, the standard local Riemann construction connects sufficiently nearby left and right states by one elementary wave from each family. These waves are ordered by their separated speeds, with intermediate constant states between them. Genuinely nonlinear fields use admissible shocks or rarefactions; linearly degenerate fields use their corresponding characteristic discontinuities. The statement is local in wave amplitude and state space.
This construction is not a theorem that arbitrary large systems data have a unique bounded entropy solution. Interactions between families, possible vacuum or degeneracy, and the available entropy structure all matter. Small-total-variation global theory in one space dimension requires its own hypotheses and estimates. The scalar Kružkov contraction does not transfer component by component to a nonlinear vector system. A convex system entropy, when available, constrains admissibility but should not be treated as a universal replacement for the complete systems theory.
Linear acoustics gives an exact two-family model
For constant acoustic density and bulk modulus, the pressure-velocity flux matrix has eigenvalues minus c and plus c. Both speeds are constant in state, so both fields are linearly degenerate. A Riemann jump decomposes exactly into a left-moving and a right-moving linear wave with an intermediate pressure and velocity. The intermediate state is constrained by preserving the appropriate characteristic invariant across each wave, rather than by separately averaging pressure and velocity without impedance weights.
The left-moving wave changes the minus characteristic variable while the plus variable retains its left-data value. The right-moving wave changes the plus variable while the minus variable retains its right-data value. Solving those two relations yields the impedance-weighted formulas below. Different impedances on the two sides create an interface-medium problem requiring modified matching conditions; the formulas here assume one homogeneous medium. This exact linear example makes family geometry concrete without claiming to simulate nonlinear gas shocks.
Work it through
A pressure step launches two families
Take K=8,rho=2, left state (p,v)=(4,0), and right state (0,0). Find the intermediate state and verify both wave jumps using eigenvectors.
The speed is c=2 and impedance Z=4. The intermediate formulas give p*=2 and v*=(4−0)/(2·4)=0.5.
For variables (p,v), choose r₋=(−4,1) and r₊=(4,1). The left-to-middle jump is (−2,0.5)=0.5r₋, traveling at speed −2.
The middle-to-right jump is (−2,−0.5)=−0.5r₊, traveling at speed +2. Both jumps satisfy the vector equation s[U]=A[U], with A=[[0,8],[0.5,0]].
Between x=−2t and x=2t the state is (2,0.5). Outside, the original left and right states remain until other interactions are introduced.
The living notebook
Explore the model
Interactive lab · vectors
One perturbation can excite two characteristic families
Pause and Reset, then set angle 0° and magnitude 2. Interpret the coordinates as nondimensional (p/Z,v). Predict the coefficients in U=α₊(1,1)+α₋(−1,1). For a fresh comparison, Reset and set angle 90° and magnitude 2, and predict the new two-family coefficients.
- x component
- 1.393
- y component
- 0.975
- Magnitude
- 1.7
Playback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. 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 pressure-axis perturbation has coefficients (1,−1), while the velocity-axis perturbation has (1,1). The instrument shows the state vector's coordinate components; it does not move the two families or solve their Riemann problem.
Pause & reason
Check your understanding
Your turn
Try it yourself
In the same medium with Z=4 and c=2, take left state (p,v)=(2,0.5) and right state (2,−0.5). Find the intermediate state and verify the two invariants inherited from the initial sides.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The pressure difference is zero, but the velocity difference contributes to the intermediate pressure.
Compare with the solution
The formulas give p*=2+(4/2)(0.5−(−0.5))=4 and v*=0. The inherited plus invariant is pL+ZvL=2+2=4, matching p*+Zv*=4. The inherited minus invariant is pR−ZvR=2−4(−0.5)=4, matching p*−Zv*=4. Thus the converging velocity perturbations produce an intermediate pressure increase even though the initial pressures agree.
Keep the idea
Take it with you
- A characteristic family is organized by an eigenvalue and its state-space eigenvector.
- Standard nonlinear Riemann constructions require local strict hyperbolicity and field hypotheses.
- Systems have interacting families and do not inherit scalar contraction componentwise.
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.