On weak solutions for the stationary Cahn-Hillard-Navier-Stokes equations with singular potential
Pith reviewed 2026-05-08 18:53 UTC · model grok-4.3
The pith
Weak solutions exist for the stationary Navier-Stokes-Cahn-Hilliard system with singular logarithmic free energy and vacuum states.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We prove the existence of weak solutions in a three-dimensional bounded domain under structural assumptions on the adiabatic exponent. The stationary setting poses two main mathematical challenges: the absence of an energy inequality driven by the evolution process to control the singular potential, and the degeneracy of the density near the vacuum. To address these issues, we introduce a specialized regularization of the logarithmic term that eliminates the quadratic growth induced by anti-diffusion, thereby restoring compactness. Uniform estimates are obtained through a special choice of artificial pressure and an interpolation argument that controls the desired norm of the density. A two-
What carries the argument
Specialized regularization of the logarithmic term that eliminates quadratic growth from anti-diffusion, together with artificial pressure and interpolation for uniform estimates, followed by a two-level limiting process.
If this is right
- Weak solutions exist that satisfy the stationary system in the distributional sense.
- The mass fraction remains in the interval between zero and one almost everywhere on the support of the density.
- Vacuum states are admissible within the class of weak solutions.
- The two-level limiting procedure yields solutions that respect the singular free-energy structure.
- The result provides the first existence theory for this steady system that includes both singularity and vacuum.
Where Pith is reading between the lines
- The regularization technique may extend to related time-dependent compressible mixture models with other singular potentials.
- The uniform estimates could inform the design of structure-preserving numerical schemes for phase-separating flows.
- Similar compactness arguments might apply to stationary models in unbounded domains once suitable decay conditions are imposed.
Load-bearing premise
The structural assumptions on the adiabatic exponent permit the specialized regularization to cancel quadratic growth from anti-diffusion and thereby restore compactness.
What would settle it
A counterexample or high-resolution numerical computation exhibiting non-existence or violation of the physical bounds for an adiabatic exponent outside the assumed structural class would falsify the result.
Figures
read the original abstract
The stationary Navier--Stokes--Cahn--Hilliard equations are considered, governing the motion of a compressible, two-phase fluid mixture with a diffuse interface. The free energy density in this paper has a singular logarithmic (Flory-uggins) form, ensuring that the mass fraction remains in the physical range and allowing for vacuum states. We prove the existence of weak solutions in a three-dimensional bounded domain under structural assumptions on the adiabatic exponent. The stationary setting poses two main mathematical challenges: the absence of an energy inequality driven by the evolution process to control the singular potential, and the degeneracy of the density near the vacuum. To address these issues, we introduce a specialized regularization of the logarithmic term that eliminates the quadratic growth induced by anti-diffusion, thereby restoring compactness. Uniform estimates are obtained through a special choice of artificial pressure and an interpolation argument that controls the desired norm of the density. A two-level limiting process then yields a weak solution that satisfies the physical bounds almost everywhere on the support of the density.To our knowledge, this is the first existence result for the steady compressible Navier--Stokes--Cahn--Hilliard system that incorporates both a singular free energy and vacuum regions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proves existence of weak solutions to the stationary compressible Navier-Stokes-Cahn-Hilliard system in a bounded 3D domain, with singular logarithmic (Flory-Huggins) free energy that permits vacuum states. The strategy combines a specialized regularization of the log term to remove quadratic anti-diffusion growth, an artificial-pressure term, an interpolation argument to obtain uniform density estimates, and a two-level limiting process (regularization parameter followed by artificial pressure) to pass to the limit while preserving physical bounds a.e. on the support of the density. Structural assumptions on the adiabatic exponent are imposed to close the estimates.
Significance. If the limiting arguments can be made rigorous, the result would constitute the first existence theorem for the steady compressible NS-CH system that simultaneously incorporates a singular free energy and vacuum regions. It directly addresses the two principal difficulties of the stationary setting: absence of an evolutionary energy inequality to control the singular potential, and degeneracy of the density near vacuum. The approach of tailored log regularization plus artificial pressure offers a potentially reusable technique for other stationary diffuse-interface models with vacuum.
major comments (2)
- [Abstract (two-level limiting process) and the section containing the a-priori estimates] The two-level limiting process (regularization parameter then artificial pressure) is load-bearing for the entire existence claim. The abstract states that uniform estimates are obtained via artificial pressure and interpolation, yet it is not clear from the outline whether the constants in the key interpolation or singular-potential estimates remain independent of the artificial-pressure parameter. If the bound on the regularized logarithmic term deteriorates as the artificial pressure tends to zero, compactness restored by the regularization may be lost before the vacuum degeneracy is controlled, exactly as flagged in the stress-test note.
- [Assumptions and the interpolation argument] The structural assumptions on the adiabatic exponent are invoked to close the interpolation argument that controls the desired norm of the density. The precise range (e.g., lower bound on gamma) must be stated explicitly and the dependence of all constants on this exponent tracked through the estimates; otherwise it is impossible to verify that the assumptions are sufficient for the 3D compactness step.
minor comments (2)
- [Regularization step] Clarify the precise definition of the specialized regularization of the logarithmic term (e.g., the cutoff or mollification parameter) and verify that it indeed eliminates quadratic growth without introducing new singularities.
- [Introduction] Add a short comparison paragraph with existing existence results for compressible NS-CH systems (even those without singular potentials or vacuum) to situate the novelty.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the constructive comments, which help clarify the presentation of the limiting arguments and assumptions. We address each major comment below.
read point-by-point responses
-
Referee: [Abstract (two-level limiting process) and the section containing the a-priori estimates] The two-level limiting process (regularization parameter then artificial pressure) is load-bearing for the entire existence claim. The abstract states that uniform estimates are obtained via artificial pressure and interpolation, yet it is not clear from the outline whether the constants in the key interpolation or singular-potential estimates remain independent of the artificial-pressure parameter. If the bound on the regularized logarithmic term deteriorates as the artificial pressure tends to zero, compactness restored by the regularization may be lost before the vacuum degeneracy is controlled, exactly as flagged in the stress-test note.
Authors: We appreciate the referee highlighting the need for greater clarity on the uniformity of estimates. The specialized regularization of the logarithmic term is constructed precisely so that the resulting bounds on the singular potential remain independent of the artificial-pressure parameter ε. The interpolation argument is applied after deriving ε-independent controls on the density and velocity fields, ensuring that compactness is retained before passing to the limit ε → 0. The two-level process proceeds by first removing the regularization for fixed ε and then letting ε vanish. We will revise the abstract and the a-priori estimates section to state explicitly that all constants are independent of ε and to detail the order of limits. revision: partial
-
Referee: [Assumptions and the interpolation argument] The structural assumptions on the adiabatic exponent are invoked to close the interpolation argument that controls the desired norm of the density. The precise range (e.g., lower bound on gamma) must be stated explicitly and the dependence of all constants on this exponent tracked through the estimates; otherwise it is impossible to verify that the assumptions are sufficient for the 3D compactness step.
Authors: We agree that the structural assumptions require more explicit statement. The analysis imposes a lower bound on the adiabatic exponent γ to close the interpolation in three dimensions. We will revise the introduction, the statement of the main theorem, and the estimates section to specify the precise range (γ > 3/2) and to include remarks tracking the dependence of all constants appearing in the estimates on γ. revision: yes
Circularity Check
No circularity: self-contained existence proof via a priori estimates and limits
full rationale
The paper derives existence of weak solutions for the stationary Navier-Stokes-Cahn-Hilliard system with singular logarithmic potential through regularization of the free energy, introduction of artificial pressure, interpolation-based uniform estimates, and a two-level limiting process. These steps rely on structural assumptions on the adiabatic exponent and standard PDE techniques to obtain compactness and pass to the limit while preserving physical bounds. No step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation; the central claim is obtained from the equations and estimates without circular reduction to inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- adiabatic exponent
axioms (1)
- standard math Standard Sobolev embeddings, interpolation inequalities, and compactness theorems in bounded domains
Lean theorems connected to this paper
-
IndisputableMonolith.Cost (Jcost = ½(x+x⁻¹)−1)washburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
F(c) = (θ_0/2)((1+c)ln(1+c) + (1−c)ln(1−c)) − (θ_c/2)c², c ∈ (−1,1).
-
Foundation/AxiomDischargePlan, Cost/FunctionalEquationn/a — no ratio-symmetric or J-cost structure invoked unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We introduce a specialized regularization of the logarithmic term that eliminates the quadratic growth induced by anti-diffusion, thereby restoring compactness. Uniform estimates are obtained through a special choice of artificial pressure (ln(1/δ))⁻¹ ρ¹¹ and an interpolation argument.
-
Constants (parameter-free derivations)n/a — adiabatic exponent is a free structural parameter, not derived unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
γ > 3/2 if ∇×g₁ = 0, γ > 5/3 if ∇×g₁ ≠ 0
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
H. Abels, On a diffuse interface model for two-phase flows of viscus, incompressible fluids with matched densities, Arch. Ration. Mech. Anal., 194, 463-506, 2009
work page 2009
-
[2]
H. Abels, Existence of weak solutions for a diffuse interface model for viscous, incompressible fluids with general densities, Comm. Math. Phys., 289, 45-73, 2009
work page 2009
- [3]
- [4]
-
[5]
H. Abels, J. Weber, Stationary Solutions for a Navier-Stokes/Cahn-Hilliard System with Singular Free En- ergies. In: Amann, H., Giga, Y., Kozono, H., Okamoto, H., Yamazaki, M. (eds) Recent Developments of Mathematical Fluid Mechanics. Advances in Mathematical Fluid Mechanics. Birkhauser, Basel, 2016
work page 2016
- [6]
- [7]
- [8]
-
[9]
D. Basaric, A. Giorgini, Global weak solutions to a compressible Navier-Stokes/Cahn-Hilliard system with singular entropy of mixing, arXiv:2506.07835v1, 2025
- [10]
-
[11]
L. Evans, Partial Differential Equations, 2nd ed., Graduate Studies in Mathematics 19, American Mathe- matical Society, Providence, RI, 2010
work page 2010
-
[12]
Evans, Weak Convergence Methods for Nonlinear Partial Differential Equations, Vol
L. Evans, Weak Convergence Methods for Nonlinear Partial Differential Equations, Vol. 74, American Math- ematical Society, 1990
work page 1990
-
[13]
E. Feireisl, A. Novotny, Singular limits in thermodynamics of viscous fluids, Advances in Mathematical Fluid Mechanics, Birkhauser, 2009
work page 2009
-
[14]
E. Feireisl, A. Novotny, H. Petzeltova, On the existence of globally defined weak solutions to the Navier-Stokes equations, J. Math. Fluid Mech., 358-392, 2001
work page 2001
- [15]
-
[16]
A. Mirnaville, R. Temam, On the Cahn-Hilliard-Oono-Navier-Stokes equations with singular potentials, Appl. Anal., 95(12), 2609-2624, 2015
work page 2015
-
[17]
S. Frigeri, On a nonlocal Cahn-Hilliard/Navier-Stokes system with degenerate mobility and singular potential for incompressible fluids with different densities, Ann. Inst. Henri Poincare Anal. Non Lin’eaire., 38(3), 647- 687, 2021
work page 2021
-
[18]
D. Gilbarg, N. Trudinger, Elliptic Partial Differential Equations of Second Order, 2nd edn., Grundlehren der Mathematischen Wissenschaften, Vol. 224, SpringerVerlag, 1983
work page 1983
-
[19]
C. Gal, A. Giorgini, M. Grasselli, A. Poiatti, Global well posedness and convergence to equilibrium for the Abels-Garcke-Grun model with nonlocal free energy, J. Math. Pures Appl., 178, 46-109, 2023
work page 2023
-
[20]
A. Giorgini, A. Miranville, R. Temam, Uniqueness and regularity for the Navier-Stokes-Cahn-Hilliard system, SIAM J. Math. Anal., 51(3), 2535-2574, 2007
work page 2007
-
[21]
A. Giorgini, R. Temam, Weak and strong solutions to the non-homogeneous incompressible Navier-Stokes- Cahn-Hilliard system, J. Math. Pures Appl., 144, 194-249, 2020
work page 2020
-
[22]
P. Hohenberg, B. Halperin, Theory of dynamic critical phenomena, Rev. Mod. Phys., 49, 435-479, 1977
work page 1977
- [23]
-
[24]
S. Ko, P. Pustejovska, E. Suli, Finite element approximation of an incompressible chemically reacting non- Newtonian fluid, Math. Mod. Numerical Appl., 52(2), 509-541, 2018
work page 2018
-
[25]
S. Ko, E. Suli, Finite element approximation of steady flows of generalized Newtonian fluids with concentration-dependent power-law index, Math. Comp., 88(317), 1061-1090, 2019
work page 2019
- [26]
- [27]
-
[28]
Lions, Mathematical topics in fluid mechanics, Vol
P. Lions, Mathematical topics in fluid mechanics, Vol. 2. Compressible models. Oxford Lecture Series in Mathematics and its Applications, Vol. 10, Oxford Science Publications/The Clarendon Press, 1998
work page 1998
- [29]
- [30]
-
[31]
A. Novotny, I. Straskraba, Introduction to the Mathematical Theory of Compressible Flow, Oxford Lecture Series in Mathematics and its Applications, Vol. 27, Oxford Univ, Press, 2004
work page 2004
-
[32]
C. Hurm, P. Knopf, A. Poiatti, Nonlocal-to-local convergence rates for strong solutions to a Navier-Stokes- CahnHilliard system with singular potential, Comm. Partial Differential Equations., 49(2), 1-40, 2024
work page 2024
-
[33]
J. Lowengrub, L. Truskinovsky, Quasi-incompressible Cahn-Hilliard fluids and topological transitions, Proc. R. Soc. Lond. Ser. A., 454, 2617-2654, 1998
work page 1998
-
[34]
P. Plotnikov, W. Weigant, Steady 3D viscous compressible flows with adiabatic exponentγ∈(1,∞), J. Math. Pures Appl., 104, 58-82, 2015
work page 2015
-
[35]
E. Stein, Singular integrals and differentiability properties of functions, Bulletin of the London Mathematical Society 5, 1973
work page 1973
-
[36]
W. Ziemer, Weakly Differentiable Functions: Sobolev Spaces and Functions of Bounded Variation, Springer, 1989 32 Z. LIANG, S. LIU, J. SHUAI, AND D. WANG School of Mathematics, Southwestern University of Finance and Economics, Chengdu 611130, China. Email address:liangzl@swufe.edu.cn School of Mathematics, Southwestern University of Finance and Economics, ...
work page 1989
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.