a-Contraction with Shifts in Shock Stability

• a-Contraction with Shifts is a weighted relative-entropy method for analyzing shock stability in hyperbolic conservation laws.
• It employs a shift function h(t) and weight a to modulate a pseudo-distance that compensates for asymmetry and neutral translation of shock profiles.
• The approach yields necessary and sufficient conditions for stability in both hyperbolic and viscous-dispersive shock settings.

In the analysis of shock stability for systems of conservation laws, $a$-contraction with shifts denotes a weighted relative-entropy method in which an entropic Rankine–Hugoniot shock $(u_L,u_R,\sigma_{LR})$ is compared against a solution $u(t,x)$ through a pseudo-distance that is both asymmetric and translation-modulated. For a weight $a>0$ and a shift $h(t)$, the central functional is

$$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$

where $\eta(u\mid v)$ is the relative entropy associated with a strictly convex entropy–entropy flux pair $(\eta,q)$. The goal is to choose $h(t)$ so that $E(t)$ is nonincreasing, thereby proving orbital $(u_L,u_R,\sigma_{LR})$0 stability up to translation; recent work develops necessary and sufficient conditions for this mechanism in the small-shock regime for general systems, especially for interior characteristic families (Faile, 4 Apr 2025).

1. Hyperbolic–entropy formulation

The underlying PDE is the one-dimensional system of conservation laws

$(u_L,u_R,\sigma_{LR})$1

with $(u_L,u_R,\sigma_{LR})$2, $(u_L,u_R,\sigma_{LR})$3 open and connected, and $(u_L,u_R,\sigma_{LR})$4. The flux Jacobian is $(u_L,u_R,\sigma_{LR})$5. The system is assumed to admit a convex entropy–entropy flux pair $(u_L,u_R,\sigma_{LR})$6 satisfying

$(u_L,u_R,\sigma_{LR})$7

together with the entropy inequality

$(u_L,u_R,\sigma_{LR})$8

in the sense of distributions, and strict convexity $(u_L,u_R,\sigma_{LR})$9 for all $u(t,x)$0 (Faile, 4 Apr 2025).

For any fixed $u(t,x)$1, the relative entropy and relative entropy flux are

$u(t,x)$2

$u(t,x)$3

The weighted pseudo-distance with shift,

$u(t,x)$4

is equivalent to a weighted $u(t,x)$5 norm on compact subsets: for any compact $u(t,x)$6 there exist $u(t,x)$7 such that

$u(t,x)$8

The shift $u(t,x)$9 is chosen to follow a generalized characteristic in continuous regions or the shock path at discontinuities. Along a smooth path $a>0$0, one has for almost all $a>0$1,

$a>0$2

and

$a>0$3

When no jump is present and $a>0$4 is chosen as a generalized characteristic of the $a>0$5-th family, $a>0$6 (Faile, 4 Apr 2025).

A standard abbreviation is

$a>0$7

With strong traces, the basic differential inequality is

$a>0$8

This is the analytic core of the method: the weight $a>0$9 compensates for the asymmetry of the two states, while the shift removes the neutral translational mode.

2. Dissipation mechanism and $h(t)$0-relative entropy stability

The method separates dissipation into a continuous contribution and a Rankine–Hugoniot contribution. When $h(t)$1, the continuous dissipation is

$h(t)$2

When $h(t)$3 crosses a discontinuity and $h(t)$4, the discontinuous dissipation is

$h(t)$5

An entropic shock $h(t)$6 is called $h(t)$7-relative entropy stable if the following hold:

$h(t)$8

$h(t)$9

A theorem attributed in the paper to Kang–Vasseur states that if, for any $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$0, there exists a Lipschitz $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$1 such that $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$2 is nonincreasing, then the shock is $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$3-relative entropy stable, hence satisfies $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$4–$$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$5 (Faile, 4 Apr 2025).

In the small-shock local-attractor regime, the target estimate becomes quantitative: $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$6 This formulation is stronger than mere nonincrease of $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$7: it identifies a coercive mechanism tied to the left trace $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$8, and it is precisely this strengthened dissipation that recent work characterizes for general systems.

3. Necessary and sufficient conditions for small shocks

For small shocks of a genuinely nonlinear $$E(t) = \int_{-\infty}^{h(t)} a\, \eta(u(t,x)\mid u_L)\,dx + \int_{h(t)}^{\infty} \eta(u(t,x)\mid u_R)\,dx,$$9-th family, the decisive criterion is a matrix inequality on the complement of the $\eta(u\mid v)$0-th eigendirection. Let $\eta(u\mid v)$1, with Rankine–Hugoniot speed $\eta(u\mid v)$2, and let $\eta(u\mid v)$3 and $\eta(u\mid v)$4 denote the right eigenvector and characteristic speed of $\eta(u\mid v)$5. A sufficient condition is the existence of $\eta(u\mid v)$6 such that

$\eta(u\mid v)$7

is negative definite on $\eta(u\mid v)$8. Under this hypothesis, there exists $\eta(u\mid v)$9 such that for all $(\eta,q)$0 one can choose

$(\eta,q)$1

and $(\eta,q)$2 so that for perturbations $(\eta,q)$3,

$(\eta,q)$4

Hence the shock is a local attractor (Faile, 4 Apr 2025).

The corresponding necessary condition has the same geometric content. If for all small $(\eta,q)$5 the shock $(\eta,q)$6 is a local attractor with weight $(\eta,q)$7, then

$(\eta,q)$8

is negative semidefinite on $(\eta,q)$9 (Faile, 4 Apr 2025).

This criterion is obtained through the analysis of a “maximal shock” $h(t)$0 and of the associated functional $h(t)$1, together with the asymptotic Hessian

$h(t)$2

on the transverse eigenspace. A plausible implication is that $h(t)$3-contraction is governed by a genuinely transverse second-order compatibility between entropy geometry, shock family, and nonlinear flux curvature.

4. Extremal families, interior families, and examples

The theory sharply distinguishes extremal and interior characteristic families. For extremal families $h(t)$4 or $h(t)$5, the operator $h(t)$6 is sign-definite on $h(t)$7, so the sufficient condition is always satisfiable. This yields $h(t)$8-contraction with appropriate weights even for large perturbations in $h(t)$9. For interior families $E(t)$0, the necessary matrix condition can fail; when it is not negative semidefinite, the paper states that for all small $E(t)$1 and any $E(t)$2, there exist nearby states $E(t)$3 with $E(t)$4 and $E(t)$5 but

$E(t)$6

so no local-attractor $E(t)$7-contraction holds (Faile, 4 Apr 2025).

The article also places the single-entropy, constant-weight theory in context. In “rich systems,” Serre–Vasseur proved contraction with shifts for any Liu–Majda stable shock by constructing two distinct entropies $E(t)$8 depending on the shock, including for intermediate families. They also treated the Euler contact discontinuities using a single entropy in a special contact structure case. By contrast, the recent small-shock theory works in general systems with a single entropy and constant weights $E(t)$9, and it shows explicitly that some interior shocks cannot be local attractors under this framework (Faile, 4 Apr 2025).

A model $(u_L,u_R,\sigma_{LR})$00 example illustrates positive results. For the system with state $(u_L,u_R,\sigma_{LR})$01, entropy $(u_L,u_R,\sigma_{LR})$02, and intermediate family $(u_L,u_R,\sigma_{LR})$03, the critical matrix becomes

$(u_L,u_R,\sigma_{LR})$04

which is negative definite on $(u_L,u_R,\sigma_{LR})$05 for $(u_L,u_R,\sigma_{LR})$06 and $(u_L,u_R,\sigma_{LR})$07. Hence small $(u_L,u_R,\sigma_{LR})$08 shocks are local attractors with $(u_L,u_R,\sigma_{LR})$09. By contrast, for 2D isentropic MHD and interior family $(u_L,u_R,\sigma_{LR})$10, the analogous matrix fails to be negative semidefinite on $(u_L,u_R,\sigma_{LR})$11 for many states, including large specific volume $(u_L,u_R,\sigma_{LR})$12, yielding cases where small interior shocks are not local attractors for any constant weight $(u_L,u_R,\sigma_{LR})$13 (Faile, 4 Apr 2025).

A common misconception is that the weight in $(u_L,u_R,\sigma_{LR})$14-contraction is always a spatial profile attached to the shock. In the small-shock hyperbolic theory just described, the weight is constant in space and depends only on the shock amplitude $(u_L,u_R,\sigma_{LR})$15; the paper states explicitly that there is no ODE for $(u_L,u_R,\sigma_{LR})$16 along the shock profile.

5. Viscous and viscous–dispersive shock applications

For the one-dimensional barotropic Navier–Stokes system in Lagrangian coordinates,

$(u_L,u_R,\sigma_{LR})$17

the method is adapted to viscous shock profiles $(u_L,u_R,\sigma_{LR})$18 by introducing a shock-aligned weight

$(u_L,u_R,\sigma_{LR})$19

together with a shifted weight $(u_L,u_R,\sigma_{LR})$20 and the weighted relative entropy

$(u_L,u_R,\sigma_{LR})$21

A central point of the 2024 analysis is that it does not employ the effective velocity variable $(u_L,u_R,\sigma_{LR})$22 even for higher order estimates. The resulting contraction estimate has the form

$(u_L,u_R,\sigma_{LR})$23

and yields global existence, convergence in sup norm toward the shifted viscous shock, and $(u_L,u_R,\sigma_{LR})$24 (Han et al., 2024).

For the Navier–Stokes–Korteweg system, the 2026 extension reintroduces an effective velocity

$(u_L,u_R,\sigma_{LR})$25

which transforms the system into two degenerate parabolic equations. The paper then applies the method of $(u_L,u_R,\sigma_{LR})$26-contraction with shifts to prove the contraction property of any large solution perturbed from a viscous-dispersive shock wave. The contraction property does not depend on the strengths of viscosity and capillarity. This uniformity yields global existence for large perturbations and zero viscosity–capillarity limits, on which Riemann shocks are unique and stable up to shifts (Eun et al., 14 Feb 2026).

Taken together, these results show that $(u_L,u_R,\sigma_{LR})$27-contraction with shifts is not merely a hyperbolic relative-entropy criterion; it is a transportable stability architecture for viscous and viscous–dispersive regularizations, with the shift absorbing translation and the $(u_L,u_R,\sigma_{LR})$28-weight encoding the asymmetry of the shock.

6. Related meanings of contraction with shifts

The phrase “contraction with shifts” also appears in other branches of analysis, but with different technical content. In nonlinear dynamical systems, contraction “with shifts” is also called transverse contraction, horizontal contraction, or contraction modulo time/phase reparametrization. There, contraction is required only in directions transverse to the flow $(u_L,u_R,\sigma_{LR})$29, while tangential directions along $(u_L,u_R,\sigma_{LR})$30 are neutral. Two trajectories are considered equivalent if there exists a reparameterization $(u_L,u_R,\sigma_{LR})$31 such that

$(u_L,u_R,\sigma_{LR})$32

and uniform transverse contraction on a compact forward-invariant set without equilibria implies existence of a unique periodic orbit together with exponential orbital convergence (Giesl et al., 2022).

In symbolic dynamics, the analogous structure is a contraction homotopy

$(u_L,u_R,\sigma_{LR})$33

satisfying

$(u_L,u_R,\sigma_{LR})$34

Here the “time parameter” is itself a binary full shift, and homotopies must be block maps commuting with the group action. The notion is used to characterize retracts of full shifts, dense periodic points under finite periodic asymptotic dimension, and several map-extension properties (Poirier et al., 2024).

A further, unrelated fixed-point usage employs “shifting distance functions” $(u_L,u_R,\sigma_{LR})$35 satisfying

$(u_L,u_R,\sigma_{LR})$36

which generalizes the Banach contraction principle on complete metric spaces (Berzig, 2013).

These usages are mathematically distinct. What they share is the replacement of literal pointwise contraction by contraction after a symmetry reduction, a shift, or a transformed comparison. In the conservation-law setting, that symmetry is spatial translation of a shock profile, and the resulting $(u_L,u_R,\sigma_{LR})$37-contraction with shifts has become a precise tool for proving orbital $(u_L,u_R,\sigma_{LR})$38 stability and for identifying when shocks are, or are not, local attractors.