REVIEW 1 major objections 5 minor 1 cited by
Global weak solutions to a compressible Navier--Stokes/Cahn--Hilliard system with singular entropy of mixing
T0 review · 1 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper establishes global-in-time weak solutions for the compressible Navier–Stokes/Cahn–Hilliard system with the singular Flory–Huggins entropy, under γ>3/2, with the phase variable confined to (−1,1) wherever the density is positive.
desk verdict First global weak solutions for compressible NS/Cahn-Hilliard with singular Flory-Huggins potential; the proof is convincing but Lemma 4.1 misstates the pressure integrability exponent. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the control of the density-weighted singular chemical-potential term through a conserved zero-mass identity. Because the relative mass M_r∈(−1,1) is conserved, the quantity ρ_ε(c_ε−M_r) has zero total mass for every approximating solution; choosing c_ε−M_r as a test function in the weak chemical-potential identity (3.8) yields a uniform $L^{2}$_t estimate for ∫ρ_εF'_ε(c_ε)(c_ε−M_r)dx (estimate (4.30)). Since F'_ε is odd, monotone, and behaves like a constant multiple of |F'_ε| outside the band |c_ε|≤K_0, choosing K_0 with |M_r|<K_0<1 converts this single estimate into uniform bounds on ∫ρ_ε|F'_ε(c_ε)|dx, then on μ_ε in $L^{2}$(0,T;$W^{{1,2}}$(Ω)), on √ρ_εF'_ε(c_ε) in $L^{2}$, and on ∇c_ε in $L^{2}$(0,T;$L^{{2p_4}}$) with p_4=3γ/(γ+3)>1. These bounds feed a pressure estimate obtained with a divergence-inverting operator and the renormalized-continuity compactness argument, yielding strong convergence of densities and concentrations; the integral lower-semicontinuity lemma then forces the strict inequality −1<c<1 on {ρ>0}.
What would settle it
Take admissible initial data with M_r=1, equivalently c0=1 on {ρ0>0}; the proof's choice of K_0 with |M_r|<K_0<1 becomes impossible, so the argument converting (4.30) into the $L^{2}$_t bound on ∫ρ_ε|F'_ε(c_ε)|dx fails. If the conclusion of Theorem 1.1 nonetheless held for such data, it would require a mechanism beyond the one constructed; conversely, exhibiting finite-energy initial data satisfying (1.20)–(1.22) for which an approximating sequence fails to keep c in [−1,1], or fails to converge to a weak solution, would refute the theorem.
Extended reading notes
Core claim
The central claim is Theorem 1.1: for any time horizon T, any bounded $C^{2}$ domain Ω⊂$ℝ^{3}$, and any adiabatic exponent γ>3/2, the Navier–Stokes/Cahn–Hilliard system (1.1)–(1.10) with free energy f(ρ,c)=ρ^γ/(γ−1)+F(c)−(θ0/2)$c^{2}$ and singular entropy F(c)=(θ/2)[(1+c)ln(1+c)+(1−c)ln(1−c)] admits a weak solution (ρ,u,c) with chemical potential μ whenever the initial data satisfy ρ0≥0, c0∈[−1,1], finite energy E0, total mass M>0, and relative mass M_r=(∫ρ0c0 dx)/M∈(−1,1). The solution obeys the renormalized continuity equation, balance of momentum, transport equation for the concentration, and the chemical-potential identity in weak form, together with the energy inequality. On the set where the density is positive, the phase variable satisfies the strict physical bound −1<c<1 almost everywhere. The proof approximates F by quadratic-growth potentials, derives uniform bounds for the approximating solutions, and passes to the limit with the standard compactness machinery for compressible Navier–Stokes equations; the key new step is a set of estimates for the density-dependent Cahn–Hilliard equation that control ∫ρF'(c) and the full $W^{{1,2}}$-norm of μ using only γ-integrability of ρ with γ>3/2.
Load-bearing premise
The proof's central assumption is that the initial relative mass M_r=(∫ρ0c0 dx)/(∫ρ0 dx) lies strictly between −1 and 1; if only one phase is present initially, so that c0≡±1 wherever ρ0>0, then no constant K_0 with |M_r|<K_0<1 exists and the key estimate controlling the singular chemical-potential term collapses.
Editorial extensions
If this is right
- For γ>3/2, the compressible model now has global weak solutions with the physically relevant logarithmic entropy, closing the gap between incompressible and compressible Navier–Stokes/Cahn–Hilliard theory for singular potentials.
- The strict bound −1<c<1 on {ρ>0} means solutions respect the interpretation of c as a difference of mass concentrations, so the singular potential can be used without separately enforcing the constraint.
- No positivity or boundedness of the initial density is needed, only finite energy and γ-integrability, so the admissible initial-data class is broader than in earlier results requiring ρ*≤ρ0≤ρ*.
- The same proof, with minor adjustments, covers two-dimensional domains, more general pressure laws, an additional H(c)lnρ mixing term, and any singular potential F∈C([−1,1])∩C^2(−1,1) with F'(±1)=±∞ and F''≥α>0, as stated in Remark 1.3.
- Because the time horizon T is arbitrary, the existence statement is global-in-time and permits arbitrarily large finite-energy initial data.
Reading between the lines
- Editorial inference: the method's reliance on K_0>|M_r| suggests the argument cannot be pushed to the pure-phase limits M_r→±1; extending existence to initial data with only one phase present would require a mechanism different from the zero-mass estimate (4.30).
- Editorial inference: the uniform bound on ∫ρ|F'(c)| may provide quantitative control of phase-separation dissipation, offering a route toward long-time convergence rates to equilibrium that the paper does not address.
- Editorial inference: the strategy of controlling a singular nonlinearity through a conserved relative-mass identity may apply to other density-dependent phase-field systems with logarithmic potentials, such as degenerate-mobility or nonlocal Cahn–Hilliard variants, although the paper does not treat those cases.
- Editorial inference: a numerical check of the regularized system for γ just above 3/2 and M_r close to 1 would test whether the constant 1/(K_0−|M_r|) in (4.31), which blows up as M_r→1, corresponds to an observable loss of uniform bounds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proves the existence of global weak solutions to a three-dimensional compressible Navier–Stokes/Cahn–Hilliard system with the singular Flory–Huggins entropy F(c)=θ/2[(1+c)ln(1+c)+(1−c)ln(1−c)]. The main theorem (Theorem 1.1) requires γ>3/2, finite-energy initial data with c0∈[−1,1], and a relative initial mass M_r=∫ρ0c0/∫ρ0 strictly between −1 and 1. The authors regularize F by the quadratic extensions F_ε of [31], use the existence theory of [4] for the approximate system, and prove uniform estimates for ρ_ε F'_ε(c_ε), μ_ε and c_ε. They then pass to the limit via the Lions–Feireisl strategy, prove strong convergence of densities and concentrations, and conclude that −1<c<1 almost everywhere on {ρ>0}.
Significance. If the exponent statements are corrected, this is a significant advance: it is the first existence result for the compressible phase-field model with the physical logarithmic potential, in the range γ>3/2, and it rigorously enforces the physical bounds on the order parameter. The proof is detailed, follows the standard Lions–Feireisl framework in a careful way, and the paper is transparent about the restriction M_r∈(−1,1) and about the approximation of F. The new chemical-potential and entropy estimates (4.30)–(4.32) are convincing and likely to be useful in later work on singular potentials in compressible phase-field models.
major comments (1)
- [§4.4, Lemma 4.1, (1.19)] The claimed pressure integrability q(γ)=min{5/3−1/γ, 3/2} is not established. The estimates in §4.4 yield q(γ)=5/3−1/γ for 3/2<γ<6 and q(γ)=4/3−1/(2γ) for γ>6, and the endpoint γ=6 is not covered by the improved estimate because the Sobolev exponent p=6γ/(7γ−6) equals 1 there, below the range allowed for Bogovskii's operator in Lemma 2.3. Since Lemma 4.1(4.13) and equation (1.19) assert the larger exponent 3/2 for γ>6 (and 1.5 at γ=6), they overclaim the result. Section A.3's use of 'estimate (4.40)' for the L^2-threshold γ≥9/5 is only correct with the improved piecewise exponent. Please correct the statements to the piecewise formula in Remark 4.2 and verify the downstream uses; the existence proof is unaffected because the derived q satisfies q>1 for all γ>3/2.
minor comments (5)
- [§4.3(v)] The displayed estimate contains the term θ0ϱεcε^2; based on (3.6) and (1.4) it should be θ0ϱεcε. Please correct the typo.
- [§4.4] The test function b(z)=z^ν is not bounded on [0,∞) as stated; please recall the standard truncation (e.g., b_k(z)=min{z^ν,k^ν}) that makes the Bogovskii-based identities legitimate.
- [§4.9] After (4.39), the phrase 'letting ψ→1' is informal because ψ∈C_c∞(0,T); the standard approximation by cutoffs should be mentioned.
- [§4.3(ii)] In the passage following (4.31), there is a missing variable in 'for any ∈(−1,1)'; it should read 'for any s∈(−1,1)'.
- [§1.2 / Remark 1.2] The paper should state explicitly that the case |M_r|=1 is left open; the proof of (4.28)–(4.31) requires a K0 with |M_r|<K0<1, so the current argument does not cover pure-phase initial data.
Assumptions & free parameters
assumptions (9)
- standard math Existence of Bogovskii operator (Lemma 2.3) on bounded Lipschitz domains
- standard math Aubin-Lions compactness lemma
- standard math Sobolev embeddings and Poincare inequalities (Lemma 2.1, 2.2, 2.4)
- standard math Feireisl's oscillation defect measure and renormalized continuity equation framework
- standard math Commutator lemma (from [20]) and Div-Curl lemma (from [19])
- domain assumption Domain Omega is bounded and of class C^2
- domain assumption Viscosity coefficients satisfy 0 < eta <= eta(c) <= bar eta, 0 <= lambda(c) <= bar lambda
- domain assumption Initial data satisfy (1.20)-(1.22), including M_r in (-1,1)
- domain assumption Thermodynamic parameters 0 < theta < theta0
Cite this review
Pith. "Pith review of Global weak solutions to a compressible Navier--Stokes/Cahn--Hilliard system with singular entropy of mixing." pith.science (2026). https://pith.science/paper/ZX7KE5ZS
@misc{pith2026250607835,
author = {Pith},
title = {Pith review of: Global weak solutions to a compressible Navier--Stokes/Cahn--Hilliard system with singular entropy of mixing},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZX7KE5ZS}},
note = {Machine review of arXiv:2506.07835}
}
abstract
We study a Navier-Stokes/Cahn-Hilliard system modeling the evolution of a compressible binary mixture of viscous fluids undergoing phase separation. The novelty of this work is a free energy potential including the physically relevant Flory-Huggins (logarithmic) entropy, as opposed to previous studies in the literature, which only consider regular potentials with polynomial growth. Our main result establishes the existence of global-in-time weak solutions in three-dimensional bounded domains for arbitrarily large initial data. The core contribution is the derivation of new estimates for the chemical potential and the Flory-Huggins entropy arising from a density-dependent Cahn-Hilliard equation under minimal assumptions: non-negative $\gamma$-integrable density with $\gamma>\frac32$. In addition, we prove that the phase variable, which represents the difference of the mass concentrations, takes value within the physical interval $(-1,1)$ almost everywhere on the set where the density is positive.
Forward citations
Cited by 1 Pith paper
-
Global strong solutions for 1D compressible Navier-Stokes/Cahn-Hilliard equations with vacuum
Global strong solutions exist and are unique for 1D compressible NS/CH with vacuum without compatibility conditions via singular-in-time weighted estimates.
Reference graph
Works this paper leans on
-
[4]
H. Abels and E. Feireisl,On a diffuse interface model for a two-phase flow of compressible viscous fluids, Indiana Univ. Math. J.57, 659–698, 2008
work page 2008
-
[31]
A. Giorgini and R. Temam,Weak and strong solutions to the nonhomogeneous incompressible Navier-Stokes- Cahn-Hilliard system, J. Math. Pures Appl.144, 194–249, 2020
work page 2020
-
[1]
H. Abels,Existence of weak solutions for a diffuse interface model for viscous, incompressible fluids with general densities, Commun. Math. Phys.289, 45–73, 2009
work page 2009
-
[2]
Abels,On a diffuse interface model for two-phase flows of viscous, incompressible fluids with matched densities, Arch
H. Abels,On a diffuse interface model for two-phase flows of viscous, incompressible fluids with matched densities, Arch. Ration. Mech. Anal.194, 463–506, 2009
2009
- [3]
- [5]
- [6]
- [7]
Show all 44 references
-
[8]
Abels, Y
H. Abels, Y. Liu and ˇS. Neˇ casov´ a,Low Mach number limit of a diffuse interface model for two-phase flows of compressible viscous fluids, GAMM-Mitteilungen.47, 2024, e202470008
2024
-
[9]
G.L. Aki, W. Dreyer, J. Giesselmann and C. Kraus,A quasi-incompressible diffuse interface model with phase transition, Math. Models Methods Appl. Sci.24, 827–861, 2014
2014
-
[10]
Boyer,Mathematical study of multi-phase flow under shear through order parameter formulation, Asympt
F. Boyer,Mathematical study of multi-phase flow under shear through order parameter formulation, Asympt. Anal.20, 175–212, 1999
1999
-
[11]
Boyer,Nonhomogeneous Cahn-Hilliard fluids, Ann
F. Boyer,Nonhomogeneous Cahn-Hilliard fluids, Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire18, 225–259, 2001
2001
-
[12]
Boyer,A theoretical and numerical model for the study of incompressible mixture flows, Comput
F. Boyer,A theoretical and numerical model for the study of incompressible mixture flows, Comput. Fluids31, 41–68, 2002
2002
-
[13]
Cherfils, E
L. Cherfils, E. Feireisl, M. Mich´ alek, A. Miranville, M. Petcu and D. Praˇ z´ ak,The compressible Navier-Stokes- Cahn-Hilliard equations with dynamic boundary conditions, Math. Models Methods Appl. Sci.29, 2557–2584, 2019
2019
-
[14]
Ding, P.D.M
H. Ding, P.D.M. Spelt and C. Shu,Diffuse interface model for incompressible two-phase flows with large density ratiosJ. Comput. Phys.226, 2078–2095, 2007
2007
-
[15]
DiPerna and P.L.-Lions,Ordinary differential equations, transport theory and Sobolev spaces, Invent
R.J. DiPerna and P.L.-Lions,Ordinary differential equations, transport theory and Sobolev spaces, Invent. Math. 98, 511–547, 1989
1989
-
[16]
ten Eikelder, E.H
M.F.P. ten Eikelder, E.H. van Brummelen and D. Schillinger,CompressibleN-phase fluid mixture models, arXiv preprint arXiv:2503.24225, 2025
2025 arXiv
-
[17]
ten Eikelder, K.G
M.F.P. ten Eikelder, K.G. van der Zee, I. Akkerman and D. Schillinger,A unified framework for Navier–Stokes Cahn–Hilliard models with non-matching densities, Math. Models Methods Appl. Sci.33, 175—221, 2023
2023
-
[18]
Elbar and A
C. Elbar and A. Poulain,Analysis and numerical simulation of a generalized compressible Cahn–Hilliard–Navier– Stokes model with friction effects, ESAIM: M2AN,58, 1989-2034, 2024
1989
-
[19]
Feireisl,Dynamics of viscous compressible fluids, Oxford University Press, Oxford; 2003
E. Feireisl,Dynamics of viscous compressible fluids, Oxford University Press, Oxford; 2003
2003
-
[20]
Feireisl,On the motion of a viscous, compressible, and heat conducting fluid, Indiana Univ
E. Feireisl,On the motion of a viscous, compressible, and heat conducting fluid, Indiana Univ. Math. J.53, 1707–1740, 2004
2004
-
[21]
Feireisl, H
E. Feireisl, H. Petzeltov´ a,On integrability up to the boundary of the weak solutions of the Navier-Stokes equations of compressible flow, Comm. Partial Differential Equations25, 755–767, 2000
2000
-
[22]
Feireisl, H
E. Feireisl, H. Petzeltov´ a, E. Rocca and G. Schimperna,Analysis of a phase-field model for two-phase compressible fluids, Math. Models Methods Appl. Sci.20, 1129–1160, 2010
2010
-
[23]
Freist¨ uhler,Phase transitions and traveling waves in compressible fluids, Arch
H. Freist¨ uhler,Phase transitions and traveling waves in compressible fluids, Arch. Ration. Mech. Anal.211, 189–204, 2014
2014
-
[24]
Frigeri,Global existence of weak solutions for a nonlocal model for two-phase flows of incompressible fluids with unmatched densities, Math
S. Frigeri,Global existence of weak solutions for a nonlocal model for two-phase flows of incompressible fluids with unmatched densities, Math. Models Methods Appl. Sci.26, 1955–1993, 2016
1955
-
[25]
Frigeri,On a nonlocal Cahn-Hilliard/Navier-Stokes system with degenerate mobility and singular potential for in-compressible fluids with different densities, Ann
S. Frigeri,On a nonlocal Cahn-Hilliard/Navier-Stokes system with degenerate mobility and singular potential for in-compressible fluids with different densities, Ann. Inst. Henri Poincar´ e Anal. Non Lin´ eaire38, 647–687, 2021
2021
-
[26]
C.G. Gal, A. Giorgini, M. Grasselli and A. Poiatti,Global well-posedness and convergence to equilibrium for the Abels-Garcke-Gr¨ un model with nonlocal free energy, J. Math. Pures Appl.178, 46–109, 2023
2023
-
[27]
Gal and M
C.G. Gal and M. Grasselli,Asymptotic behavior of a Cahn-Hilliard-Navier-Stokes in 2D, Ann. Inst. Henri Poincar´ e Anal. Non Lin´ eaire27, 401–436, 2010
2010
-
[28]
C.G. Gal, M. Grasselli and H. Wu,Global weak solutions to a diffuse interface model for incompressible two-phase flows with moving contact lines and different densities, Arch. Ration. Mech. Anal.234, 1–56, 2019
2019
-
[29]
Galdi,An introduction to the mathematical theory of the Navier-Stokes equations, I, Springer-Verlag, New- York; 1994
G.P. Galdi,An introduction to the mathematical theory of the Navier-Stokes equations, I, Springer-Verlag, New- York; 1994. 29
1994
-
[30]
Giorgini, A
A. Giorgini, A. Miranville and R. Temam,Uniqueness and Regularity for the Navier–Stokes–Cahn–Hilliard system, SIAM J. Math. Anal.51, 2535–2574, 2019
2019
-
[32]
Gomez and K.G
H. Gomez and K.G. van der Zee,Computational phase-field modeling, Encyclopedia of Computational Mechanics Second Edition, pages 1–35, 2018
2018
-
[33]
Gurtin, D
M.E. Gurtin, D. Polignone and J. Vi˜ nals,Two-phase binary fluids and immiscible fluids described by an order parameter, Math. Models Methods Appl. Sci.6, 815–831, 1996
1996
-
[34]
Heida, J
M. Heida, J. M´ alek and K.R. Rajagopal,On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework, Z. Angew. Math. Phys.63, 145–169, 2012
2012
-
[35]
Hohenberg and B.I
P.C. Hohenberg and B.I. Halperin,Theory of dynamic critical phenomena, Rev. Mod. Phys.49, 435–479, 1977
1977
-
[36]
C. Hurm, P. Knopf and A. Poiatti,Nonlocal-to-local convergence rates for strong solutions to a Navier-Stokes- Cahn-Hilliard system with singular potential, Comm. Partial Differential Equations49, 832–971, 2024
2024
-
[37]
Kotschote,Strong solutions of the Navier-Stokes equations for a compressible fluid of Allen-Cahn type, Arch
M. Kotschote,Strong solutions of the Navier-Stokes equations for a compressible fluid of Allen-Cahn type, Arch. Ration. Mech. Anal.206, 489–514, 2012
2012
-
[38]
Kotschote and R
M. Kotschote and R. Zacher,Strong solutions in the dynamical theory of compressible fluid mixtures, Math. Models Methods Appl. Sci.25, 1217–1256, 2015
2015
-
[39]
Liang and D
Z. Liang and D. Wang,Stationary Cahn–Hilliard–Navier–Stokes equations for the diffuse interface model of compressible flows, Math. Models Methods Appl. Sci.30, 2445–2486, 2020
2020
-
[40]
Liang and D
Z. Liang and D. Wang,Weak solutions to the stationary Cahn–Hilliard/Navier–Stokes equations for compressible fluids, J. Nonlinear Sci.32, 1–25, 2022
2022
-
[41]
2: Compressible Models, Oxford Science Publication, Oxford; 1998
P.L.-Lions,Mathematical Topics in Fluid Dynamics, Vol. 2: Compressible Models, Oxford Science Publication, Oxford; 1998
1998
-
[42]
Lowengrub and L
J. Lowengrub and L. Truskinovsky,Quasi–incompressible Cahn–Hilliard fluids and topological transitions, Proc. Royal Soc. London A: Math. Phys. Eng. Sci.454, 2617–2654, 1998
1998
-
[43]
J. Shen, X. Yang and Q. Wang,Mass and volume conservation in phase field models for binary fluids, Commun. Comput. Phys.13, 1045–1065, 2013
2013
-
[44]
Shokrpour Roudbari, G
M. Shokrpour Roudbari, G. S ¸im¸ sek, E.H. van Brummelen and K.G. van der Zee,Diffuse-interface two-phase flow models with different densities: A new quasi-incompressible form and a linear energy-stable method, Math. Models Meth. Appl. Sci.28, 733–770, 2017. 30
2017
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.