[Paper Review] Remarks on the formation and decay of multidimensional shock waves

В. Г. Данилов|ArXiv.org|Dec 27, 2005

Advanced Mathematical Physics Problems2 references5 citations

TL;DR

This paper generalizes the one-dimensional shock wave formation and decay mechanism to multidimensional scalar conservation laws with convex nonlinearity. By introducing a level-set formulation based on characteristics and a generalized implicit equation for the shock front, it constructs piecewise-smooth solutions where shock formation occurs when trajectories converge at finite time, with conditions for absolute nonstability derived from second derivatives of flux functions.

ABSTRACT

In this paper, we present a formula describing the formation and decay of shock wave type solutions in some special cases.

Motivation & Objective

• To extend the one-dimensional shock wave formation and decay framework to multidimensional scalar conservation laws with convex flux functions.
• To derive a multidimensional analog of the implicit equation governing shock formation in 1D, based on characteristic trajectories and level surfaces.
• To establish conditions under which a piecewise-constant initial condition evolves into a shock wave, including criteria for absolute nonstability.
• To formalize the transition from initial data with a discontinuous step to a shock wave via smooth intermediate profiles defined on characteristic bundles.
• To provide a foundation for constructing weak asymptotic solutions in higher dimensions, analogous to those in 1D.

Proposed method

• Introduces a generalized implicit equation (5) involving second derivatives of flux functions and a scalar function K(s) to describe shock formation in nD.
• Uses characteristic systems (6) to trace trajectories from two disjoint hypersurfaces Γ₁ and Γ₂, with u₁ evolving linearly along them.
• Defines the shock time t₀(s) = 1/K(s) as the time when characteristics from Γ₁ intersect those from Γ₂, forming a shock front.
• Constructs initial data (8) as a three-region function: constant U before Γ₁, u₀₀ after Γ₂, and u₁(x) in between, using characteristic functions H⁻, H⁰, H⁺.
• Applies Lemma 1 to show that the Jacobian determinant J vanishes at t = 1/K(s), signaling wave breaking and shock formation.
• Derives conditions for absolute nonstability (13) based on the sign of the difference in second derivatives of flux functions at the jump states.

Experimental results

Research questions

• RQ1 How can the 1D shock formation mechanism based on implicit equations be generalized to multidimensional scalar conservation laws?
• RQ2 What conditions ensure that a nonstable step function evolves into a shock wave in nD, and how is the shock front defined?
• RQ3 What role do the second derivatives of the flux functions play in determining the stability of the shock front?
• RQ4 How do characteristic trajectories from two disjoint hypersurfaces determine the time and location of shock formation?
• RQ5 Under what geometric and analytic conditions does the solution transition from a smooth intermediate profile to a shock wave?

Key findings

• Shock formation occurs at time t₀(s) = 1/K(s) when characteristic trajectories from Γ₁ and Γ₂ intersect, with K(s) > 0 ensuring convergence.
• The solution evolves into a shock wave of the form u = U + H(S(x,t))(u₀₀ - U), where S(x,t) = 0 defines the shock front.
• The jump across the shock is absolutely nonstable if (f′′uu(u₀₀))limₓ→x̄u(x,t) − (f′′uu(U))limₓ→x̄u(x,t) > 0, indicating instability under perturbations.
• The intermediate solution u₁(x) satisfies a first-order PDE (5) that generalizes the 1D implicit equation f′(u₁(x)) = −Kx + b.
• The construction reduces the nD problem to a 1D problem along characteristic curves, with explicit integration of system (6) yielding u₁ = U − K(s)τ and Xᵢ(τ,s) = χᵢ¹(s) + (1/K)(f′ᵢ(U) − f′ᵢ(U − Kτ)).
• For t = 1/K(s), the shock front forms at the intersection of trajectories, and the time difference between arrival at Γ₁ and Γ₂ is governed by the inner product expression involving f′(U) and f′(u₀₀).