Global weak solutions to a diffuse-interface model for quasi-incompressible two-phase flows with unmatched densities and singular potential
Pith reviewed 2026-05-07 11:37 UTC · model grok-4.3
The pith
Global weak solutions exist for a quasi-incompressible diffuse-interface model of two-phase flows with unmatched densities and singular potentials.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the initial-boundary value problem in T^3 with a class of physically relevant singular free energy densities, global-in-time weak solutions exist. The proof reduces the original system to a Korteweg-type fluid model, applies a two-layer approximation, and obtains delicate estimates for mass density and the phase-field variable inspired by the Bresch-Desjardins entropy. Capillarity at the free interface damps the density evolution, while tail estimates exclude concentrations of the singular potential even without a priori pressure integrability. This yields the first existence result for the Navier-Stokes/Cahn-Hilliard type system with unmatched densities and mass-averaged velocity in the
What carries the argument
The central mechanism is the reduction of the original Navier-Stokes/Cahn-Hilliard system to a Korteweg-type fluid model, combined with a two-layer approximation and Bresch-Desjardins entropy estimates that control the density and phase field while exploiting capillarity damping.
If this is right
- Global existence holds for the quasi-incompressible system with singular potentials that are physically relevant.
- The model admits weak solutions even though pressure integrability is unavailable a priori.
- Capillarity supplies a damping mechanism that prevents uncontrolled density growth.
- The result covers the unmatched-density case without requiring spatial regularization.
Where Pith is reading between the lines
- The techniques may extend to related diffuse-interface models that incorporate temperature or other conserved quantities.
- Numerical methods preserving the Bresch-Desjardins entropy structure could inherit global existence from the continuous level.
- The quasi-incompressible formulation might be compared quantitatively with fully compressible or incompressible approximations for the same physical regimes.
Load-bearing premise
The reduction to the Korteweg-type model together with the two-layer approximation converges to a weak solution of the original system without loss of the key physical structure.
What would settle it
An explicit initial datum on the torus for which the density or singular potential concentrates in finite time, violating the weak formulation or the energy inequality, would disprove global existence.
Figures
read the original abstract
We study a thermodynamically consistent diffuse-interface model that describes the motion of two macroscopically immiscible, incompressible, and viscous Newtonian fluids with unmatched densities. This model is compatible with continuum mixture theory. It adopts a mass-averaged (barycentric) velocity so that the two-phase flow is quasi-incompressible: the velocity is no longer divergence-free, and the pressure enters the equation of the chemical potential. For the initial-boundary value problem in $\mathbb{T}^3$ with a class of physically relevant singular free energy densities, we prove the existence of global-in-time weak solutions. The proof relies on a suitable reduction of the original system to a Korteweg-type fluid model combined with a two-layer approximation, together with delicate estimates for the mass density and the phase-field variable inspired by the celebrated Bresch-Desjardins entropy. A key observation is that capillarity at the free interface provides a damping effect on the density evolution. For the limiting procedure, we derive delicate tail estimates to exclude possible concentrations of the singular potential, since no integrability of the pressure is available \textit{a priori}. This work appears to be the first existence result for the Navier-Stokes/Cahn-Hilliard type system with unmatched densities and mass-averaged velocity without spatial regularization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proves the existence of global-in-time weak solutions to the initial-boundary value problem on the three-torus for a thermodynamically consistent diffuse-interface model of quasi-incompressible two-phase flows with unmatched densities and a class of singular free-energy densities. The proof proceeds by reducing the system to a Korteweg-type fluid model, introducing a two-layer approximation, deriving estimates for the mass density and phase field via a Bresch-Desjardins-type entropy, and obtaining tail controls on the singular potential that exploit capillarity damping; the limiting procedure is justified by excluding concentrations in the absence of a priori pressure integrability.
Significance. If the estimates hold, the result supplies the first global existence theorem for the Navier-Stokes/Cahn-Hilliard system with mass-averaged velocity and unmatched densities without spatial regularization. It therefore extends the existing theory for matched-density or regularized cases and supplies a mathematically rigorous foundation for continuum-mixture models of two-phase flows with physically relevant singular potentials.
major comments (2)
- [Limiting procedure (after the two-layer approximation)] The tail estimates that prevent concentrations of the singular potential (invoked because pressure integrability is unavailable a priori) are load-bearing for the passage to the limit; the manuscript should make explicit the quantitative damping rate furnished by the capillarity term and verify that it is uniform with respect to the approximation parameters.
- [Reduction step and entropy estimates] In the reduction to the Korteweg-type system, the additional terms arising from unmatched densities must be shown to remain controlled under the Bresch-Desjardins entropy; any extra commutator or lower-order contribution should be estimated explicitly to confirm that the entropy dissipation still yields the required a priori bounds.
minor comments (3)
- [Introduction] The claim that the result is the first of its kind would benefit from a short, explicit comparison (in the introduction) with the closest prior works on matched densities or regularized models.
- [Model formulation] Notation for the chemical potential and the pressure gradient terms should be unified between the original system and the reduced Korteweg formulation to avoid reader confusion.
- [Weak formulation] A few typographical inconsistencies appear in the statement of the weak formulation (e.g., the precise sense in which the chemical-potential equation is tested); these can be corrected without altering the argument.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. The positive assessment is appreciated, and we address the two major comments below by providing the requested clarifications and indicating the corresponding revisions.
read point-by-point responses
-
Referee: [Limiting procedure (after the two-layer approximation)] The tail estimates that prevent concentrations of the singular potential (invoked because pressure integrability is unavailable a priori) are load-bearing for the passage to the limit; the manuscript should make explicit the quantitative damping rate furnished by the capillarity term and verify that it is uniform with respect to the approximation parameters.
Authors: We agree that an explicit quantitative statement strengthens the argument. In the revised version we will add a dedicated lemma (placed immediately before the passage to the limit) that computes the damping rate induced by the capillarity term in the density equation. The lemma will show that the resulting integral bound is independent of both approximation parameters and yields a uniform control on the singular potential that rules out concentrations. The proof of uniformity follows directly from the structure of the two-layer approximation and the Bresch-Desjardins entropy already derived. revision: yes
-
Referee: [Reduction step and entropy estimates] In the reduction to the Korteweg-type system, the additional terms arising from unmatched densities must be shown to remain controlled under the Bresch-Desjardins entropy; any extra commutator or lower-order contribution should be estimated explicitly to confirm that the entropy dissipation still yields the required a priori bounds.
Authors: We thank the referee for this observation. While the reduction is carried out in Section 3 and the entropy estimates appear in Section 4, the commutator terms generated by the unmatched-density contributions were treated somewhat concisely. In the revision we will insert an explicit calculation (new Proposition 3.2) that bounds each commutator and lower-order term by quantities already controlled by the Bresch-Desjardins entropy. The resulting estimates confirm that the dissipation structure remains intact and that the a priori bounds are unaffected. revision: yes
Circularity Check
No circularity in existence proof
full rationale
The paper establishes global weak solutions for the quasi-incompressible NS/CH system via reduction to a Korteweg-type model, two-layer approximation, and estimates inspired by the external Bresch-Desjardins entropy. These steps are drawn from prior independent literature and do not reduce to self-definitions, fitted parameters renamed as predictions, or load-bearing self-citations. The tail controls on the singular potential and capillarity damping are derived within the proof without assuming the target existence result. The derivation remains self-contained against external mathematical benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Initial data satisfy sufficient regularity and compatibility conditions for the weak formulation to make sense.
- domain assumption The singular free energy densities belong to a physically relevant class permitting the Bresch-Desjardins-type estimates.
Reference graph
Works this paper leans on
-
[1]
Abels, Existence of weak solutions for a diffuse interface model for viscous, incompressible fluids with general densities, Comm
H. Abels, Existence of weak solutions for a diffuse interface model for viscous, incompressible fluids with general densities, Comm. Math. Phys.289(2009), no. 1, 45–73
2009
-
[2]
Abels, Strong well-posedness of a diffuse interface model for a viscous, quasi-incompressible two-phase flow, SIAM J
H. Abels, Strong well-posedness of a diffuse interface model for a viscous, quasi-incompressible two-phase flow, SIAM J. Math. Anal.44(2012), 316–340
2012
-
[3]
Abels, D
H. Abels, D. Depner and H. Garcke, Existence of weak solutions for a diffuse interface model for two-phase flows of incompressible fluids with different densities, J. Math. Fluid Mech.15(2013), no. 3, 453–480
2013
-
[4]
Abels, D
H. Abels, D. Depner and H. Garcke, On an incompressible Navier–Stokes/Cahn–Hilliard system with degen- erate mobility, Ann. Inst. H. Poincaré C Anal. Non Linéaire30(2013), no. 6, 1175–1190
2013
-
[5]
Abels, H
H. Abels, H. Garcke and G. Grün, Thermodynamically consistent, frame indifferent diffuse interface models for incompressible two-phase flows with different densities, Math. Models Methods Appl. Sci.22(2012), no. 3, 1150013, 40 pp
2012
-
[6]
Abels, H
H. Abels, H. Garcke and A. Poiatti, Mathematical analysis of a diffuse interface model for multi-phase flows of incompressible viscous fluids with different densities, J. Math. Fluid Mech.26(2024), no. 2, Paper No. 29
2024
-
[7]
Abels, H
H. Abels, H. Garcke and J. Wittmann, Diffuse interface models for two-phase flows with phase transition: modeling and existence of weak solutions, J. Math. Fluid Mech.28(2026), no. 7, 1–37
2026
-
[8]
Abels, H
H. Abels, H. Garcke and A. Giorgini, Global regularity and asymptotic stabilization for the incompressible Navier–Stokes–Cahn–Hilliard model with unmatched densities, Math. Ann.389(2024), no. 2, 1267–1321. 40 MINGWEN FEI, XIANG FEI, YADONG LIU, AND HAO WU
2024
-
[9]
Abels and J
H. Abels and J. Weber, Local well-posedness of a quasi-incompressible two-phase flow, J. Evol. Equ.21 (2021), 3477–3502
2021
-
[10]
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(2014), no. 5, 827–861
2014
-
[11]
D. M. Anderson, G. B. McFadden and A. A. Wheeler, Diffuse-interface methods in fluid mechanics, Ann. Rev. Fluid Mech.30(1998), 139–165
1998
-
[12]
Antonelli and S
P. Antonelli and S. Spirito, Global existence of weak solutions to the Navier–Stokes–Korteweg equations, Ann. Inst. H. Poincaré C Anal. Non Linéaire39(2022), no. 1, 171–200
2022
-
[13]
Bresch and B
D. Bresch and B. Desjardins, Existence of global weak solutions for a 2D viscous shallow water equations and convergence to the quasi-geostrophic model, Comm. Math. Phys.238(2003), no. 1-2, 211–223
2003
-
[14]
Bresch and B
D. Bresch and B. Desjardins, On the existence of global weak solutions to the Navier–Stokes equations for viscous compressible and heat conducting fluids, J. Math. Pures Appl. (9)87(2007), no. 1, 57–90
2007
-
[15]
Bresch, B
D. Bresch, B. Desjardins and C.-K. Lin, On some compressible fluid models: Korteweg, lubrication, and shallow water systems, Comm. Partial Differential Equations28(2003), no. 3-4, 843–868
2003
-
[16]
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. Fluids 31(2002), 41–68
2002
-
[17]
J. W. Cahn and J. E. Hilliard, Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys. 28(1958), 258–267
1958
-
[18]
N. F. Carnahan and K. E. Starling, Equation of state for nonattracting rigid spheres. J. Chem. Phys.51 (1969), 635–636
1969
-
[19]
H. Ding, P. D. M. Spelt and C. Shu, Diffuse interface model for incompressible two-phase flows with large density ratios, J. Comput. Phys.226(2007), 2078–2095
2007
-
[20]
21, (eds
Q.DuandX.-B.Feng, Chapter5–Thephasefieldmethodforgeometricmovinginterfacesandtheirnumerical approximations, in Handbook of Numerical Analysis, Vol. 21, (eds. A. Bonito and R. H. Nochetto), Elsevier, (2020), 425–508
2020
-
[21]
M. Fei, X. Fei, D. Han and Y. Liu, Local-in-time existence of strong solutions to a quasi-incompressible Cahn– Hilliard–Navier–Stokes system, to appear in Z. Angew. Math. Mech., 2026, DOI:10.1002/zamm.70443
- [22]
-
[23]
Feireisl, Y
E. Feireisl, Y. Lu and J. Málek, On PDE analysis of flows of quasi-incompressible fluids, ZAMM Z. Angew. Math. Mech.96(2016), no. 4, 491–508
2016
-
[24]
Feireisl, A
E. Feireisl, A. Novotný and H. Petzeltová, On the existence of globally defined weak solutions to the Navier– Stokes equations, J. Math. Fluid Mech.3(2001), 358–392
2001
-
[25]
Feireisl and H
E. Feireisl and H. Petzeltová, On integrability up to the boundary of the weak solutions of the Navier–Stokes equations of compressible flow, Comm. Partial Differential Equations25(2000), no. 3-4, 755–767
2000
-
[26]
Feireisl, H
E. Feireisl, H. Petzeltová, E. Rocca and G. Schimperna, Analysis of a phase-field model for two-phase com- pressible fluids, Math. Models Methods Appl. Sci.20(2010), no. 7, 1129–1160
2010
-
[27]
Feireisl and P
E. Feireisl and P. Zhang, Quasi-neutral limit for a model of viscous plasma, Arch. Ration. Mech. Anal.197 (2010), no. 1, 271–295
2010
-
[28]
Freistühler and M
H. Freistühler and M. Kotschote, Phase-field and Korteweg-type models for the time-dependent flow of compressible two-phase fluids, Arch. Ration. Mech. Anal.224(2017), no. 1, 1–20
2017
-
[29]
Freistühler and M
H. Freistühler and M. Kotschote, Phase-field descriptions of two-phase compressible fluid flow: interstitial working and a reduction to Korteweg theory, Quart. Appl. Math.77(2019) no. 3, 489–496
2019
-
[30]
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(2016), no. 10, 1955–1993
2016
-
[31]
Frigeri, On a nonlocal Cahn–Hilliard/Navier–Stokes system with degenerate mobility and singular potential for incompressible fluids with different densities, Ann
S. Frigeri, On a nonlocal Cahn–Hilliard/Navier–Stokes system with degenerate mobility and singular potential for incompressible fluids with different densities, Ann. Inst. H. Poincaré C Anal. Non Linéaire38(2021), no. 3, 647–687
2021
-
[32]
C. G. Gal, A. Giorgini, M. Grasselli and A. Poiatti, Global well-posedness and convergence to equilibrium for the Abels-Garcke-Grün model with nonlocal free energy, J. Math. Pures Appl. (9)178(2023), 46–109
2023
-
[33]
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(2019), no. 1, 1–56
2019
-
[34]
C. G. Gal, M.-Y. Lv and H. Wu, On a thermodynamically consistent diffuse interface model for two-phase incompressible flows with non-matched densities: Dynamics of moving contact lines, surface diffusion, and mass transfer, J. Differential Equations,465(2026), 114203, 70 pages
2026
-
[35]
C. G. Gal and A. Poiatti, Unified framework for the separation property in binary phase-segregation processes with singular entropy densities, European J. Appl. Math.36(2025), 40–67
2025
- [36]
-
[37]
Geissert, H
M. Geissert, H. Heck and M. Hieber, On the equationdivu=gand Bogovski˘ ı’s operator in Sobolev spaces of negative order, inPartial Differential Equations and Functional Analysis, 113–121, Oper. Theory Adv. Appl., 168, Birkhäuser, Basel, 2006. NSCH SYSTEM FOR QUASI-INCOMPRESSIBLE TWO-PHASE FLOWS 41
2006
-
[38]
Giorgini, Well-posedness of the two-dimensional Abels–Garcke–Grün model for two-phase flows with un- matched densities, Calc
A. Giorgini, Well-posedness of the two-dimensional Abels–Garcke–Grün model for two-phase flows with un- matched densities, Calc. Var. Partial Differ. Equ.60(2021), Paper No.100, 40 pp
2021
-
[39]
Giorgini, Existence and stability of strong solutions to the Abels–Garcke–Grün model in three dimensions, Interfaces Free Bound.24(2022), 565–608
A. Giorgini, Existence and stability of strong solutions to the Abels–Garcke–Grün model in three dimensions, Interfaces Free Bound.24(2022), 565–608
2022
-
[40]
Giorgini, M
A. Giorgini, M. Grasselli and H. Wu, The Cahn–Hilliard–Hele–Shaw system with singular potential, Ann. Inst. H. Poincaré C Anal. Non Linéaire.35(2018) no. 4, 1079–1118
2018
-
[41]
A. Giorgini, J.-N. He and H. Wu, On a thermodynamically consistent diffuse-interface model for incompress- ible two-phase flows with chemotaxis and mass transport, preprint, 2025, arXiv:2512.22585
-
[42]
Giorgini and P
A. Giorgini and P. Knopf, Two-phase flows with bulk–surface interaction: Thermodynamically consistent Navier–Stokes–Cahn–Hilliard models with dynamic boundary conditions, J. Math. Fluid Mech.25(2023), no. 3, Paper No. 65
2023
-
[43]
Gisclon and I
M. Gisclon and I. Lacroix-Violet, About the barotropic compressible quantum Navier–Stokes equations, Nonlinear Anal.128(2015), 106–121
2015
-
[44]
Gomez and K
H. Gomez and K. G. van der Zee, Computational phase-field modeling, Encyclopedia of Computational Mechanics, 2nd edn., 2016
2016
-
[45]
M. Grasselli and A. Poiatti, Convergence to equilibrium of weak solutions to the Cahn–Hilliard equation with non-degenerate mobility and singular potential, preprint, 2025, arXiv:2510.17296
-
[46]
P. C. Hohenberg and B. I. Halperin, Theory of dynamic critical phenomena, Rev. Mod. Phys.49(1977), 435–479
1977
-
[47]
Lam and H
K.-F. Lam and H. Wu, Thermodynamically consistent Navier–Stokes–Cahn–Hilliard models with mass trans- fer and chemotaxis, European J. Appl. Math.29(2018), no. 4, 595–644
2018
-
[48]
Liu and D
Y. Liu and D. Trautwein, On a diffuse interface model for incompressible viscoelastic two-phase flows, J. Nonlinear Sci.35(2025), no. 1, Paper No. 12, 71 pp
2025
-
[49]
Londen and H
S.-O. Londen and H. Petzeltová, Regularity and separation from potential barriers for the Cahn-Hilliard equation with singular potential, J. Evol. Equ.18(2018), no. 3, 1381–1393
2018
-
[50]
J. S. Lowengrub and L. Truskinovsky, Quasi-incompressible Cahn–Hilliard fluids and topological transitions, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci.454(1998), no. 1978, 2617–2654
1998
-
[51]
Miranville,The Cahn–Hilliard equation, Recent advances and applications
A. Miranville,The Cahn–Hilliard equation, Recent advances and applications. CBMS-NSF Regional Conf. Ser. in Appl. Math., 95, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, 2019
2019
-
[52]
Miranville and S
A. Miranville and S. Zelik, Robust exponential attractors for Cahn–Hilliard type equations with singular potentials, Math. Meth. Appl. Sci.27(2004) 545–582
2004
-
[53]
P. B. Mucha, M. Pokorný and E. Zatorska, Chemically reacting mixtures in terms of degenerated parabolic setting, J. Math. Phys.54(2013), no. 7, 071501, 17 pp
2013
-
[54]
P. B. Mucha, M. Szlenk and E. Zatorska, Construction of weak solutions to a model of pressureless viscous flow driven by nonlocal attraction-repulsion, J. Math. Pures Appl. (9)195(2025), Paper No. 103671, 59 pp
2025
-
[55]
Novotny and I
A. Novotny and I. Straškraba,Introduction to the Mathematical Theory of Compressible Flow, vol. 27 of Oxford Lecture Series in Mathematics and its Applications, Oxford University Press, Oxford, 2004
2004
-
[56]
J. Shen, X. Yang and Q. Wang, Mass and volume conservation in phase field models for binary fluids, Commun. Comput. Phys.13(2013), 1045–1065
2013
-
[57]
Shokrpour Roudbari, G
M. Shokrpour Roudbari, G. Şimşek, 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 Methods Appl. Sci.28(2018), no. 4, 733–770
2018
-
[58]
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(2023), no. 1, 175–221
2023
-
[59]
Truesdell,Rational Thermodynamics, Springer-Verlag, New York, 1984
C. Truesdell,Rational Thermodynamics, Springer-Verlag, New York, 1984
1984
-
[60]
J. D. van der Waals, The thermodynamic theory of capillarity under the hypothesis of a continuous variation of density, J. Stat. Phys.20(1979), 197–244
1979
-
[61]
A. F. Vasseur and C. Yu, Existence of global weak solutions for 3D degenerate compressible Navier–Stokes equations, Invent. Math.206(2016), no. 3, 935–974
2016
-
[62]
A. F. Vasseur and C. Yu, Global weak solutions to the compressible quantum Navier–Stokes equations with damping, SIAM J. Math. Anal.48(2016), no. 2, 1489–1511
2016
-
[63]
P. Yue, J. J. Feng, C. Liu and J. Shen, A diffuse-interface method for simulating two-phase flows of complex fluids, J. Fluid Mech.515(2004), 293–317. 42 MINGWEN FEI, XIANG FEI, YADONG LIU, AND HAO WU School of Mathematics and Statistics, Anhui Normal University, Wuhu 241002, P. R. China Email address:mwfei@ahnu.edu.cn School of Mathematics and Statistics...
2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.