Similarity Solutions of Shock Formation for First-order Strictly Hyperbolic Systems

Basic Information

• Paper ID: 2511.00672
• Title: Similarity Solutions of Shock Formation for First-order Strictly Hyperbolic Systems
• Authors: Jun Eshima, Luc Deike, Howard A. Stone (Princeton University)
• Classification: math.AP (Analysis of PDEs), math-ph (Mathematical Physics), math.MP (Mathematical Physics)
• Submission Date: November 1, 2025
• Paper Link: https://arxiv.org/abs/2511.00672

Abstract

Shock phenomena arising from hyperbolic partial differential equations are ubiquitous in mathematics and science. A classical example is shock formation in the inviscid Burgers equation $\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}=0$. Previous research has demonstrated that when shocks form in the inviscid Burgers equation, the dynamics near the shock singularity in space and time are locally self-similar and universal—that is, equivalent regardless of initial conditions. This paper proves that, in fact, shock formation in general first-order strictly hyperbolic partial differential equations in one-dimensional space is also self-similar and universal, with self-similarity analogous to the inviscid Burgers equation. The paper derives and verifies analytical formulas for the self-similar universal solutions.

Research Background and Motivation

Problem Background

1. Universality of Shock Phenomena: Shock phenomena arising from hyperbolic partial differential equations appear in numerous physical domains, including fluid mechanics (shallow water equations, compressible gas equations), soft matter physics (nonlinear elasticity), magnetohydrodynamics (astrophysics and plasma physics), and traffic flow.
2. Shock Formation Mechanism: Traditional mathematical research has focused primarily on fundamental questions regarding shock uniqueness, existence, and well-posedness, but has given less attention to the question of how shocks form. Understanding the shock formation mechanism has practical value for understanding subsequent regularization.
3. Leading Role of Burgers Equation: Research by Pomeau and by Eggers and Fontelos on the inviscid Burgers equation has demonstrated that shock formation possesses a locally self-similar universal solution structure, providing a foundation for understanding shock formation in more general systems.

Research Motivation

1. Universality Question: Although shock formation in the Burgers equation is known to be universal, whether all strictly hyperbolic systems possess similar universality remains an open problem.
2. Theoretical Refinement: There is a need to establish a general theoretical framework that extends the results from the Burgers equation to broader hyperbolic systems.
3. Practical Applications: To provide a unified mathematical description tool for shock phenomena in various physical systems.

Core Contributions

1. Theoretical Universality Proof: Proves that shock formation in first-order strictly hyperbolic partial differential equation systems in one-dimensional space possesses local self-similarity and universality.
2. Analytical Formula Derivation: Derives analytical expressions for self-similar universal solutions, with forms analogous to the Burgers equation solution.
3. Methodological Innovation: Provides a systematic method for analyzing shock formation in general hyperbolic systems.
4. Numerical Verification: Validates theoretical predictions through numerical examples using the shallow water equations.

Detailed Methodology

Problem Formulation

Consider a first-order strictly hyperbolic partial differential equation system: $\frac{\partial \mathbf{f}}{\partial t} = \mathbf{M}(\mathbf{f}) \cdot \frac{\partial \mathbf{f}}{\partial x}$

where $\mathbf{f}(x,t) = (f_1(x,t), \ldots, f_N(x,t))$ represents N dependent variables, and $\mathbf{M}(\mathbf{f})$ is an $N \times N$ matrix. The system is called strictly hyperbolic when the matrix $\mathbf{M}$ is diagonalizable and possesses distinct real eigenvalues.

Theoretical Framework

1. Leading-Order Expansion Analysis

Near the shock formation point $(x_*, t_*)$, introduce local variables: $(x', \tau, \mathbf{f}') = (x-x_*, t_*-t, \mathbf{f}-\mathbf{f}_*)$

The leading-order linearized equation is: $\frac{\partial \mathbf{f}'}{\partial \tau} + \mathbf{M}|_{\mathbf{f}=\mathbf{f}_*} \frac{\partial \mathbf{f}'}{\partial x'} = 0$

This is a linear advection equation that does not produce shocks; therefore, higher-order terms must be considered.

2. Higher-Order Nonlinear Analysis

The equation including the next-order nonlinear terms is: $\frac{\partial \mathbf{f}'_i}{\partial \tau} + \sum_{j=1}^N M_{ij}|_{\mathbf{f}=\mathbf{f}_*} \frac{\partial \mathbf{f}'_j}{\partial x'} + \sum_{j=1}^N \sum_{k=1}^N M_{ij,k}|_{\mathbf{f}=\mathbf{f}_*} f'_k \frac{\partial \mathbf{f}'_j}{\partial x'} = 0$

3. Self-Similar Transformation

Through balance analysis, determine the scaling relations: $x = O(\tau^\alpha)$, $\mathbf{f}' = O(\tau^{\alpha-1})$, where $\alpha > 1$.

In the shock reference frame $x = x' - \lambda\tau$ (where $\lambda$ is a particular eigenvalue), the solution expansion takes the form: $\mathbf{f}'(x,\tau) = g(x,\tau)\mathbf{e} + \mathbf{q}\tau^{\alpha-1} + \mathbf{h}(x,\tau)$

where $\mathbf{e}$ is the eigenvector corresponding to eigenvalue $\lambda$.

4. Reduction to Burgers Equation

Through projection by the left eigenvector $\mathbf{e}^L$, the system reduces to: $\frac{\partial g}{\partial \tau} + c_1\tau^{\alpha-1}\frac{\partial g}{\partial x} - c_2 g\frac{\partial g}{\partial x} = 0$

After appropriate coordinate transformation, the standard Burgers equation is obtained: $\frac{\partial g}{\partial \tau} - g\frac{\partial g}{\partial x_s} = 0$

Main Results

The local self-similar solution for shock formation is: $\mathbf{f}(x,t) = \mathbf{f}_* + (t_*-t)^{1/2} F\left(\frac{x-x_*-\lambda(t_*-t)}{c(t_*-t)^{3/2}}\right) \mathbf{e}$

where the function $F(\xi)$ satisfies: $-\xi = F + KF^3$

The constant $c$ can be analytically determined from the matrix $\mathbf{M}$ and its derivatives: $c = -\frac{\sum_{i,j,k} M_{ij,k}|_{\mathbf{f}=\mathbf{f}_*} e^L_i e_j e_k}{\sum_i e^L_i e_i}$

Technical Innovations

1. Systematic Methodology: Systematically extends the analytical method for the Burgers equation to general hyperbolic systems.
2. Eigenvector Direction: Proves that shocks form along the eigenvector direction, a key application of strict hyperbolicity.
3. Matching Conditions: Determines undetermined constants in the solution by requiring matching between local and far-field solutions.
4. Analytical Expressions: Provides explicit analytical expressions for the coefficient $c$, avoiding numerical fitting.

Experimental Setup

Verification Example: Shallow Water Equations

The one-dimensional shallow water equations are selected as the verification example: $\frac{\partial u}{\partial t} = -u\frac{\partial u}{\partial x} - \frac{\partial \eta}{\partial x}$$\frac{\partial \eta}{\partial t} = -\eta\frac{\partial u}{\partial x} - u\frac{\partial \eta}{\partial x}$

where $u(x,t)$ is velocity and $\eta(x,t)$ is water surface elevation.

Initial and Boundary Conditions

• Initial conditions: $(u(x,0), \eta(x,0)) = (\sin(2\pi x), 1)$
• Computational domain: $x \in [0,1]$ with periodic boundary conditions
• Shock formation time: $t_* \approx 0.196$

Numerical Method

High-precision numerical methods are employed to solve the shallow water equations, tracking the shock formation process to precision levels of $t_* - t = 10^{-4}$.

Experimental Results

Eigenvalue and Eigenvector Verification

The matrix for the shallow water equations is: $\mathbf{M} = \begin{pmatrix} -u & -1 \\ -\eta & -u \end{pmatrix}$

The eigenvalues and eigenvectors are: $\lambda^{(\pm)} = -u_* \pm \sqrt{\eta_*}, \quad \mathbf{e}^{(\pm)} = (1, \mp\sqrt{\eta_*})$

Numerical results show that the shock forms along the $\mathbf{e}^{(-)}$ direction.

Power-Law Verification

The theoretically predicted divergence behavior of first-order derivatives: $\max\left|\frac{\partial f_i}{\partial x}\right| = \frac{|e_i c^{-1}|}{t_* - t}$

The divergence behavior of second-order derivatives: $\max\left|\frac{\partial^2 f_i}{\partial x^2}\right| = \frac{|e_i c^{-2}|}{(t_* - t)^{5/2}} \frac{25\sqrt{15}}{108} K^{1/2}$

Numerical results show perfect agreement with the power-law exponents predicted by theory.

Self-Similarity Verification

By rescaling solutions at different times according to self-similar coordinates, the self-similar convergence of the solution is verified. The fitted constant $K \approx 0.14$ is consistent with theoretical predictions.

Key Findings

1. Universality Confirmed: Shock formation in the shallow water equations indeed follows the same self-similar laws as the Burgers equation.
2. Analytical Predictions Accurate: The power-law exponents and coefficients given by theory match numerical results precisely.
3. Parameter Determination: The unique unknown parameter $K$ can be determined through numerical fitting.

Related Work

Foundations of Shock Theory

• Pioneering work by Challis, Stokes, Riemann, and others established the foundations of shock theory
• Dafermos's monograph systematically presents the theory of hyperbolic conservation laws

Self-Similar Solution Research

• Pioneering research on self-similar solutions of the Burgers equation by Pomeau and by Eggers and Fontelos
• Similar self-similarity discovered in specific hyperbolic systems (such as thin-film Marangoni flows)

Singularity Formation Theory

• John's theory on the inevitable formation of shocks in strictly hyperbolic systems
• Applications of finite-time singularity formation in problems such as the Navier-Stokes equations

Conclusions and Discussion

Main Conclusions

1. Universality Theorem: Shock formation in first-order strictly hyperbolic partial differential equation systems possesses universal self-similarity.
2. Burgers-Type Behavior: Shock formation in all such systems is analogous to the Burgers equation.
3. Analytical Solvability: Self-similar solutions can be completely determined through analytical methods (except for one constant).

Limitations

1. Dimensional Restriction: The current theory applies only to one-dimensional space.
2. Strict Hyperbolicity Assumption: Requires the condition that eigenvalues are strictly distinct.
3. Smoothness Requirements: Requires the matrix function $\mathbf{M}(\mathbf{f})$ to be smooth at the shock point.
4. Higher-Order Corrections: When certain components of eigenvectors vanish, higher-order correction terms must be considered.

Future Directions

1. Repeated Eigenvalue Cases: Extend to hyperbolic systems with repeated eigenvalues.
2. Multidimensional Generalization: Study self-similar shock formation in multidimensional space.
3. Application Extension: Apply the method to more physical systems.

In-Depth Evaluation

Strengths

1. Theoretical Depth: Provides rigorous mathematical proofs that generalize classical results from the Burgers equation to the general case.
2. Systematic Methodology: Establishes a complete framework for analyzing shock formation in general hyperbolic systems.
3. Practical Value: Analytical formulas are convenient for use in various physical applications.
4. Comprehensive Numerical Verification: Theoretical predictions are thoroughly validated through the shallow water equations.

Weaknesses

1. Limited Applicability: Restricted to one-dimensional strictly hyperbolic systems.
2. Numerical Precision Limitations: Verification experiments are limited by numerical method precision, achieving only $t_* - t = 10^{-4}$.
3. Degenerate Cases: Limited discussion of degenerate cases where $c_2 = 0$.

Impact

1. Theoretical Contribution: Provides important supplementation to singularity theory for hyperbolic partial differential equations.
2. Application Prospects: Provides a unified tool for shock research in fluid mechanics, plasma physics, and other fields.
3. Methodological Value: Demonstrates how to systematically analyze singularity behavior in nonlinear partial differential equations.

Applicable Scenarios

1. Fluid Mechanics: Shock analysis in compressible flows, shallow water waves, and similar systems.
2. Plasma Physics: Shock phenomena in magnetohydrodynamics.
3. Traffic Flow Theory: Congestion formation in nonlinear traffic flow models.
4. Numerical Methods: Provides theoretical foundation for shock-capturing algorithms.

References

The paper cites important literature in the field, including:

• Dafermos's monograph on hyperbolic conservation laws
• Classical work by Eggers and Fontelos on singularities in partial differential equations
• John's pioneering research on singularity formation in nonlinear wave propagation
• Recent findings in specific physical systems

This paper possesses significant theoretical value in mathematical physics, providing a unified mathematical framework for understanding shock phenomena in a wide range of physical systems. Its rigorous theoretical derivations and detailed numerical verification make it an important contribution to the field.