advanced36 min studyLesson 419 in Mathematics

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

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 F:Rm→RmF:\mathbb R^m\to\mathbb R^m, the classical system is Ut+DF(U)Ux=0U_t+DF(U)U_x=0. 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.

DF(U)ri(U)=λi(U)ri(U),∇λi⋅ri≠0 (genuinely nonlinear),∇λi⋅ri≡0 (linearly degenerate)DF(U)r_i(U)=\lambda_i(U)r_i(U),\qquad\nabla\lambda_i\cdot r_i\ne0\ \text{(genuinely nonlinear)},\qquad\nabla\lambda_i\cdot r_i\equiv0\ \text{(linearly degenerate)}
The properties are imposed in a specified state neighborhood, not checked only at an isolated state. Eigenvector orientation can be chosen locally for a genuinely nonlinear field.↔ Scroll wide equations horizontally.

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.

UL=U(0)⟶U(1)⟶⋯⟶U(m)=URU_L=U^{(0)}\longrightarrow U^{(1)}\longrightarrow\cdots\longrightarrow U^{(m)}=U_R
Under the stated local strict-hyperbolicity and field assumptions, consecutive states are connected by an elementary wave of the corresponding family. This is not an arbitrary-large-data assertion.↔ Scroll wide equations horizontally.

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.

p∗=pL+pR2+Z2(vL−vR),v∗=vL+vR2+pL−pR2Zp_* =\frac{p_L+p_R}{2}+\frac Z2(v_L-v_R),\qquad v_* =\frac{v_L+v_R}{2}+\frac{p_L-p_R}{2Z}
In one constant acoustic medium, the intermediate state satisfies p*+Zv*=pL+ZvL and p*−Zv*=pR−ZvR. Wave speeds are −c and +c.↔ Scroll wide equations horizontally.

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.

  1. The speed is c=2 and impedance Z=4. The intermediate formulas give p*=2 and v*=(4−0)/(2·4)=0.5.

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

  3. 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]].

U∗=(2,0.5),U∗−UL=12(−4,1),UR−U∗=−12(4,1)U_*=(2,0.5),\qquad U_*-U_L=\tfrac12(-4,1),\qquad U_R-U_*=-\tfrac12(4,1)
The signs and different family vectors determine two opposite propagation directions from a single pressure jump.↔ Scroll wide equations horizontally.
The result

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.

Explore the model

Interactive lab

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.

One perturbation can excite two characteristic familiesPlayback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. Current readings: x component: 1.393; y component: 0.975; Magnitude: 1.7. All parameters are adjustable with the labeled controls below.−7.46−3.7303.737.46−3.5−1.7501.753.5xyvvₓ = 1.39
Vectorx / y projectionsTrace marker
x component
1.393
y component
0.975
Magnitude
1.7
t = 0 s

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.

Check your understanding

What distinguishes a linearly degenerate characteristic field in the neighborhood used for the local Riemann construction?

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.

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.