Nonlinear stability of compressible vortex sheets in three-dimensional elastodynamics

Technical Summary: Nonlinear Stability of Compressible Vortex Sheets in Three-Dimensional Elastodynamics

Problem Statement
This paper investigates the nonlinear stability and local existence of compressible vortex sheet solutions for three-dimensional (3D) isentropic elastic flows. The governing system is a nonlinear hyperbolic system of conservation laws describing inviscid compressible elastodynamics, where the fluid interface is a characteristic free boundary (a contact discontinuity). Across this interface, the tangential velocity exhibits a jump, while the normal velocity and pressure remain continuous.

The authors address the transition from two-dimensional (2D) to fully three-dimensional settings. While 2D compressible vortex sheets in Euler flows are known to be violently unstable for Mach numbers $M < \sqrt{2}$ and stable only in the supersonic regime, the inclusion of elasticity introduces new physical mechanisms. Previous work established weakly linear stability for 3D elastic vortex sheets under specific conditions. However, proving nonlinear stability in 3D presents significantly more intricate challenges due to the additional spatial dimension, which complicates frequency interactions between elasticity and velocity, and introduces higher-order degeneracies in the stability analysis that are absent in 2D.

Methodology
The analysis employs a combination of microlocal analysis, paradifferential calculus, and the Nash-Moser iteration scheme. The core of the methodology involves the following steps:

1. Linearization and Variable Transformation: The problem is first linearized around a piecewise constant background state (a planar vortex sheet). To handle the moving free boundary, the authors introduce a coordinate transformation (flattening the interface) and utilize Alinhac's "good unknowns" to eliminate the loss of derivatives associated with the free boundary term.
2. Paralinearization: The variable-coefficient linearized system is treated using Bony's paradifferential calculus. This allows the system to be approximated by a pseudo-differential operator system, facilitating energy estimates in weighted Sobolev spaces ($H^s_\gamma$).
3. Microlocal Analysis of the Lopatinski˘i Condition: A critical component is the detailed study of the roots of the Lopatinski˘i determinant. The authors identify specific geometric stability conditions associated with the deformation gradient tensor. They analyze the interaction between "poles" (where the system cannot be reduced to a non-characteristic form) and "roots" (where the Lopatinski˘i determinant vanishes).
4. Upper Triangularization and Diagonalization: To handle the variable coefficients and the loss of derivatives, the authors employ an upper triangularization technique. This isolates outgoing modes into a closed system where they appear only at the leading order. Furthermore, they refine Coulombel's diagonalization framework to construct suitable transformations that reduce the degeneracy order of the Lopatinski˘i matrix, particularly in regions where double roots coincide with double poles.
5. Energy Estimates: Using localized G˚arding-type inequalities and symmetrizers, the authors derive energy estimates that control the characteristic components despite the loss of tangential derivatives.
6. Nash-Moser Iteration: Finally, the local existence and nonlinear stability are established using a Nash-Moser iteration scheme, which is necessary to overcome the loss of regularity inherent in the linearized problem.

Key Contributions and Novelty
The paper highlights two major novelties in addressing the challenges of the 3D elastic vortex sheet problem:

1. Restoration of Ellipticity via Elasticity: In 3D compressible Euler vortex sheets, the front symbol exhibits degenerate ellipticity in certain frequency directions, making standard energy estimates insufficient to control the regularity of the free interface. The authors demonstrate that the non-parallel structure of the deformation gradient tensor ($F_1 \times F_2 \neq 0$) provides an ellipticity enhancement mechanism. This geometric condition allows for the recovery of one derivative in the regularity estimates of the free interface, a property not guaranteed in the purely Eulerian case.
2. Resolution of Higher-Order Singularities: Unlike the 2D case, the 3D setting features a co-dimension one set in frequency space where a double root of the Lopatinski˘i determinant coincides with a double pole. This "root-pole collision" creates a higher-order degeneracy that prevents the direct application of standard diagonalization techniques. The authors resolve this by constructing specific invertible mappings (symbols of type $\Gamma^2_0$) that transform the Lopatinski˘i matrix into a form where the degeneracy is reduced. This enables the use of localized G˚arding estimates to control the characteristic components.

Main Results
The authors establish the local existence and nonlinear stability of 3D compressible elastic vortex sheets under small initial perturbations. Specifically, they prove that if the background state satisfies:

• A geometric condition: $F_1 \times F_2 \neq 0$ (non-parallel deformation gradients).
• Stability conditions: $\bar{v}^2 < |\Pi^\perp_{F_2}(F_1)|^2/4$ and $\bar{v}^2 < G(F_1, F_2)$, where $G$ is a function derived from the deformation gradient and sound speed.

Then, for sufficiently small initial data in appropriate Sobolev spaces, there exists a unique solution to the vortex sheet problem on a time interval $[0, T]$. The solution maintains the regularity required for the Nash-Moser scheme, with the free interface and fluid variables belonging to specific weighted Sobolev spaces.

Significance
The paper claims that its results extend the understanding of multidimensional hyperbolic systems with characteristic free boundaries. By rigorously establishing the nonlinear stability of 3D elastic vortex sheets, the work demonstrates that elasticity can act as a stabilizing mechanism against the Kelvin-Helmholtz instability, provided specific geometric conditions on the deformation gradient are met. The techniques developed, particularly the refined diagonalization for handling double root-pole collisions and the exploitation of elasticity to recover ellipticity, are presented as tools potentially applicable to other free-boundary models with similar algebraic structures. The authors note that while their conditions are sufficient for nonlinear stability, they are more restrictive than the necessary and sufficient conditions for linear stability derived in their previous work, and extending the nonlinear result to the full linear stability domain remains an open challenge.