A Regularized Framework and Admissible Solutions for Liquid-Vapor Phase Transitions in Steady Compressible Flows

1. Problem Statement

The paper investigates the well-posedness and solution structure of the steady-state compressible isentropic Navier-Stokes system under the van der Waals equation of state.

• Core Difficulty: Unlike ideal gases, the van der Waals equation introduces a non-monotonic pressure-density relationship at subcritical temperatures ($0 < \theta < \theta_c$). This non-monotonicity creates a region of mechanical instability (the "hump" in the $p-v$ diagram) where pressure increases with volume.
• Physical Context: This instability drives liquid-vapor phase transitions. In the inviscid limit, the mathematical problem becomes ill-posed (mixed hyperbolic-elliptic type), leading to non-uniqueness of solutions and physical instabilities.
• Key Question: Among the infinite set of mathematically valid weak solutions (piecewise constant functions), which ones represent physically admissible phase transitions? Specifically, how does the system select a stable two-phase state with distinct interfaces?

2. Methodology

The authors employ a regularization approach combined with rigorous mathematical analysis to resolve the non-uniqueness and identify physical solutions.

• Artificial Viscosity: They introduce an artificial viscosity term ($\epsilon^2 v_{xx}$) into the steady-state equations. This transforms the ill-posed inviscid problem into a well-posed elliptic/parabolic system, allowing for the construction of smooth approximate solutions.
• Lagrangian Coordinates: The system is reformulated in Lagrangian mass coordinates to decouple density and velocity, simplifying the analysis to a single second-order ordinary differential equation (ODE) for specific volume $v$.
• Variational Framework: The problem is cast as a minimization of an energy functional:
  $$\min \int_{-1}^{1} \left( W(v) + \frac{\epsilon^2}{2} v_x^2 \right) dx$$
  subject to a constraint on the integral average of the specific volume ($\bar{v}$). Here, $W(v)$ is the potential energy derived from the pressure function.
• Analytical Tools:
  • Phase-Plane Analysis: Used to study the dynamics of the ODE and identify periodic orbits corresponding to phase transitions.
  • Calculus of Variations: Used to prove the existence of minimizers and analyze the second variation to determine stability.
  • Level-Set and Anti-Derivative Techniques: Employed to handle the singular limits and construct the relationship between the average volume and the solution structure.
  • Singular Limit Analysis: Investigating the behavior as the viscosity coefficient $\epsilon \to 0$ to recover the sharp-interface solution.

3. Key Contributions

A. Classification of Solution Regimes

The paper rigorously establishes that the existence and uniqueness of solutions depend entirely on the integral average of the specific volume ($\bar{v}$) relative to the Maxwell region ($\Omega_{Maxwell} = (\alpha_0, \beta_0)$), defined by the Maxwell equal-area rule:

1. Stable Region ($\bar{v} \notin \Omega_{Maxwell}$): If $\bar{v}$ lies outside the coexistence region, the unique solution is a constant single-phase state ($v(x) \equiv \bar{v}$). No phase transition occurs.
2. Maxwell Region ($\bar{v} \in \Omega_{Maxwell}$): If $\bar{v}$ lies inside the coexistence region, the inviscid system admits multiple piecewise constant solutions (phase transitions). The solution is a piecewise constant function taking only the two values $\alpha_0$ (liquid) and $\beta_0$ (vapor).

B. Definition of Admissible Solutions

To resolve the non-uniqueness in the Maxwell region, the authors define Admissible Solutions:

• They identify that among the infinite mathematical solutions, only those with a single-peak or single-valley profile (representing exactly two sharp interfaces) are physically meaningful.
• These solutions correspond to the global energy minimizers of the regularized system.
• Solutions with more than two interfaces (e.g., $2N$ interfaces) or constant solutions are shown to have higher energy and are unstable.

C. Convergence and Interface Formation

The paper proves that as the artificial viscosity $\epsilon \to 0$:

• The smooth approximate solutions (single-peak or single-valley) converge pointwise to the piecewise constant admissible solution.
• The diffuse interfaces in the approximate solution sharpen into discontinuities (shocks) at specific locations determined by the conservation of mass (the average $\bar{v}$).
• The width of the transition layer scales exponentially with $\epsilon$ (i.e., $O(e^{-C/\epsilon})$).

4. Main Results (Theorems)

• Theorem 1.1 (Existence/Uniqueness in Inviscid Limit):
  • Unique constant solution if $\bar{v} \in \Omega_{stable}$.
  • Multiple piecewise constant solutions if $\bar{v} \in \Omega_{Maxwell}$.
• Theorem 1.2 (Well-Posedness of Regularized System):
  • For large viscosity, the solution is unique and trivial.
  • For small viscosity and $\bar{v} \in \Omega_{unstable}$, non-trivial smooth solutions exist.
• Theorem 1.3 (Convergence to Admissible Solutions):
  • There exists a unique family of two-smooth-interface solutions (single-peak or single-valley) for the regularized system.
  • As $\epsilon \to 0$, these converge to the unique two-sharp-interface admissible solution of the inviscid system.
  • This establishes a nucleation mechanism: The non-monotonicity of pressure combined with the average volume condition triggers the formation of a stable two-phase state with two interfaces.

5. Significance

• Mathematical Rigor: The paper provides a rigorous mathematical justification for the Maxwell construction in the context of steady compressible flows, moving beyond heuristic arguments to a proof based on energy minimization and singular limits.
• Resolution of Non-Uniqueness: It offers a clear criterion (energy minimization via artificial viscosity) to select the physically relevant solution from the infinite set of weak solutions, a long-standing issue in phase transition problems.
• Physical Insight: It confirms that the instability of the pressure function acts as a nucleation mechanism. The system naturally evolves to a state with exactly two interfaces (one liquid-vapor, one vapor-liquid) to minimize energy, explaining the "double-interface" structures observed in numerical simulations.
• Framework for Future Work: The proposed regularized framework and definition of admissible solutions provide a robust foundation for analyzing more complex phase transition problems, including multi-dimensional flows and dynamic (time-dependent) cases.

In summary, the paper successfully bridges the gap between the mathematical ill-posedness of non-monotonic fluid models and physical reality by demonstrating that energy minimization selects a unique, stable phase transition profile characterized by two sharp interfaces.