A Thermodynamically Consistent Free Boundary Model for Two-Phase Flows in an Evolving Domain with Bulk-Surface Interaction
Pith reviewed 2026-05-19 06:40 UTC · model grok-4.3
The pith
The paper derives a thermodynamically consistent free boundary model for two-phase flows that includes bulk-surface material interactions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We derive from local mass balance a system of Navier-Stokes equations in the bulk coupled to convective Cahn-Hilliard equations on the bulk and surface. The model is thermodynamically consistent, satisfies local mass balance, allows material transfer between bulk and surface, produces variable contact angles, and recovers previous models by dropping dynamic boundary conditions or fixing the domain boundary.
What carries the argument
The coupled bulk-surface convective Cahn-Hilliard system with interaction terms designed for thermodynamic consistency, together with the generalized Navier slip boundary condition.
If this is right
- Local mass balance is satisfied for the components.
- Material can transfer between the bulk fluid and the surface.
- Contact angles at the boundary vary dynamically.
- Earlier models without dynamic boundaries or with fixed domains emerge as special cases.
Where Pith is reading between the lines
- Similar derivation techniques may apply to other systems involving interfaces with mass exchange, such as in multiphase porous media.
- Testing the model against experiments with measurable surface adsorption rates could validate the transfer terms.
Load-bearing premise
The interaction terms between bulk and surface phase fields are assumed to take a form that guarantees the required energy dissipation while maintaining consistency with mass balance.
What would settle it
A calculation or simulation demonstrating that the energy dissipation law fails to hold under material transfer between bulk and surface would disprove the thermodynamic consistency.
Figures
read the original abstract
We derive a thermodynamically consistent model, which describes the time evolution of a two-phase flow in an evolving domain. The movement of the free boundary of the domain is driven by the velocity field of the mixture in the bulk, which is determined by a Navier--Stokes equation. In order to take interactions between bulk and boundary into account, we further consider two materials on the boundary, which may be the same or different materials as those in the bulk. The bulk and the surface materials are represented by respective phase-fields, whose time evolution is described by a bulk-surface convective Cahn--Hilliard equation. This approach allows for a transfer of material between bulk and surface as well as variable contact angles between the diffuse interface in the bulk and the boundary of the domain. To provide a more accurate description of the corresponding contact line motion, we include a generalized Navier slip boundary condition on the velocity field. Based on local mass balance laws, we derive our model from scratch in two different ways: by the Lagrange Multiplier Approach and (in the case of matched densities and no mass flux between bulk and surface) by the Energetic Variational Approach. We further show that our model generalizes previous models from the literature, which can be recovered from our system by either dropping the dynamic boundary conditions or assuming a static boundary of the domain.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives a thermodynamically consistent free-boundary model for two-phase incompressible flows in an evolving domain. Bulk flow is governed by Navier–Stokes with a generalized Navier slip condition; materials in the bulk and on the surface are represented by phase fields whose evolution is given by a coupled bulk-surface convective Cahn–Hilliard system that permits mass transfer across the interface. The model is obtained from local mass balance via a Lagrange-multiplier formulation (general case) and an energetic variational formulation (matched densities, no mass flux). It recovers several earlier models by specializing the boundary conditions or fixing the domain.
Significance. If the energy-dissipation identity is established for the general case with nonzero mass flux and distinct densities, the work supplies a unified, first-principles framework that simultaneously incorporates material exchange, variable contact angles, and dynamic contact-line motion. Such a model is of clear interest for the mathematical analysis of diffuse-interface flows with evolving boundaries.
major comments (2)
- [§3] §3 (Lagrange-multiplier derivation): the integration-by-parts identities that close the energy dissipation law when the surface velocity differs from the bulk trace and the domain is moving rely on the precise functional form of the bulk-surface coupling terms. Please exhibit the explicit cancellation of the cross terms for arbitrary densities and nonzero mass flux; the current argument appears to invoke the matched-density simplification used in the energetic-variational section.
- [§4] §4 (energy law): the dissipation identity is stated to be non-positive, yet the proof sketch is given only under the assumption of matched densities. A separate verification for the general case (distinct densities, active mass transfer) is required, because the curvature and transport terms introduced by the evolving domain may otherwise produce uncontrolled contributions.
minor comments (2)
- [Introduction] The notation for the surface phase field and the precise definition of the mass-transfer flux should be introduced with a single, self-contained table or diagram in the introduction.
- [Introduction] Several references to prior phase-field models with dynamic boundary conditions are cited only by author-year; please add the specific equation numbers from those works that are recovered as special cases.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address the two major comments below and will revise the manuscript to strengthen the presentation of the energy-dissipation law.
read point-by-point responses
-
Referee: [§3] §3 (Lagrange-multiplier derivation): the integration-by-parts identities that close the energy dissipation law when the surface velocity differs from the bulk trace and the domain is moving rely on the precise functional form of the bulk-surface coupling terms. Please exhibit the explicit cancellation of the cross terms for arbitrary densities and nonzero mass flux; the current argument appears to invoke the matched-density simplification used in the energetic-variational section.
Authors: We appreciate the referee highlighting this point. The Lagrange-multiplier derivation in §3 is constructed for the general case with arbitrary densities and nonzero mass flux between bulk and surface. The bulk-surface coupling terms are chosen precisely so that the cross terms cancel upon integration by parts, even when the surface velocity differs from the bulk trace and the domain evolves. To improve clarity and address the concern directly, we will add an explicit step-by-step computation of these cancellations in the revised manuscript, without invoking the matched-density assumption used later in the energetic-variational section. revision: yes
-
Referee: [§4] §4 (energy law): the dissipation identity is stated to be non-positive, yet the proof sketch is given only under the assumption of matched densities. A separate verification for the general case (distinct densities, active mass transfer) is required, because the curvature and transport terms introduced by the evolving domain may otherwise produce uncontrolled contributions.
Authors: We agree that the current sketch of the energy law focuses on the matched-density case. Although the Lagrange-multiplier formulation guarantees thermodynamic consistency for the general system, an explicit verification for distinct densities and active mass transfer is indeed desirable to control the additional curvature and transport terms arising from the moving domain. In the revised version we will supply a complete, self-contained proof of the dissipation identity in the general case, showing that all such terms cancel or are absorbed into the non-positive dissipation. revision: yes
Circularity Check
Derivation from local mass balance laws via Lagrange multipliers and energetic variational approach is self-contained
full rationale
The paper derives its system from scratch starting from local mass balance laws, employing two independent methods (Lagrange multiplier approach for the general case and energetic variational approach for matched densities with no mass flux). The bulk-surface coupling terms are constructed within this derivation to enforce thermodynamic consistency and mass conservation, rather than being fitted or renamed from the target result. The model generalizes prior literature by dropping dynamic boundary conditions or fixing the domain, providing independent content. Self-citations to phase-field models supply background but do not load-bear the central claims or force uniqueness via internal theorems.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Local mass balance laws hold separately in the bulk and on the surface, with possible exchange terms between them.
- domain assumption The total energy functional includes bulk and surface contributions whose time derivative yields a non-positive dissipation.
invented entities (1)
-
Bulk-surface interaction terms in the phase-field equations
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We derive our model from scratch by considering local mass balance and energy dissipation laws. Finally, the derivation is completed via the Langrange multiplier approach.
-
IndisputableMonolith/Foundation/ArrowOfTime.leanentropy_monotone echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
Sufficiently regular solutions of system (1.1) satisfy the energy-dissipation law (1.10)
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]
- [2]
-
[3]
H. Abels and H. Garcke. Weak solutions and diffuse interface models for incompressible two-phase flows. In Handbook of mathematical analysis in mechanics of viscous fluids , pages 1267–1327. Springer, Cham, 2018
work page 2018
- [4]
- [5]
- [6]
- [7]
- [8]
-
[9]
H. Abels and J. Weber. Local well-posedness of a quasi-incompressible two-phase flow. J. Evol. Equ. , 21(3):3477–3502, 2021
work page 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(5):827–861, 2014
work page 2014
-
[11]
E. Bachini, V. Krause, I. Nitschke, and A. Voigt. Derivation and simulation of a two-phase fluid deformable surface model. J. Fluid Mech. , 977:Paper No. A41, 37, 2023
work page 2023
-
[12]
E. Bachini, V. Krause, and A. Voigt. The interplay of geometry and coarsening in multicomponent lipid vesicles under the influence of hydrodynamics. Physics of Fluids , 35(4):042102, 04 2023
work page 2023
-
[13]
E. B¨ ansch, K. Deckelnick, H. Garcke, and P. Pozzi.Interfaces: Modeling, Analysis, Numerics. Oberwolfach Seminars. Birkh¨ auser Cham, 2023
work page 2023
- [14]
- [15]
-
[16]
J. W. Barrett, H. Garcke, and R. N¨ urnberg. Finite element approximation for the dynamics of fluidic two-phase biomembranes. ESAIM Math. Model. Numer. Anal. , 51(6):2319–2366, 2017
work page 2017
-
[17]
J. W. Barrett, H. Garcke, and R. N¨ urnberg. Parametric finite element approximations of curvature-driven interface evolutions. In Geometric partial differential equations. Part I , volume 21 of Handb. Numer. Anal., pages 275–423. Elsevier/North-Holland, Amsterdam, [2020] ©2020
work page 2020
-
[18]
K. Binder and H. L. Frisch. Dynamics of surface enrichment: A theory based on the Kawasaki spin- exchange model in the presence of a wall. Z. Phys. B , 84:403—-418, 1991
work page 1991
-
[19]
F. Boyer. A theoretical and numerical model for the study of incompressible mixture flows. Comput. & Fluids, 31(1):41–68, 2002
work page 2002
-
[20]
N. Bullerjahn. Error estimates for full discretization by an almost mass conservation technique for Cahn– Hilliard systems with dynamic boundary conditions. Preprint: arXiv:2502.03847 [math.AP], 2025
-
[21]
N. Bullerjahn and B. Kov´ acs. Error estimates for full discretization of cahn–hilliard equation with dynamic boundary conditions. IMA Journal of Numerical Analysis , page draf009, 04 2025
work page 2025
-
[22]
D. Caetano and C. M. Elliott. Cahn-Hilliard equations on an evolving surface. European J. Appl. Math. , 32(5):937–1000, 2021
work page 2021
-
[23]
D. Caetano, C. M. Elliott, M. Grasselli, and A. Poiatti. Regularization and separation for evolving surface Cahn-Hilliard equations. SIAM J. Math. Anal. , 55(6):6625–6675, 2023
work page 2023
-
[24]
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(14):2557–2584, 2019
work page 2019
- [25]
- [26]
- [27]
- [28]
- [29]
- [30]
-
[31]
H. Ding, P. D. Spelt, and C. Shu. Diffuse interface model for incompressible two-phase flows with large density ratios. J. Comput. Phys. , 226(2):2078–2095, 2007
work page 2078
- [32]
-
[33]
C. M. Elliott and T. Sales. Navier–Stokes–Cahn–Hilliard equations on evolving surfaces. Interfaces Free Bound., online first, 2024
work page 2024
-
[34]
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(7):1129–1160, 2010
work page 2010
-
[35]
H. Fischer, P. Maass, and W. Dieterich. Novel Surface Modes in Spinodal Decomposition. Phys. Rev. Lett., 79:893–896, 1997
work page 1997
-
[36]
H. Fischer, J. Reinhard, W. Dieterich, J.-F. Gouyet, P. Maass, A. Majhofer, and D. Reinel. Time- dependent density functional theory and the kinetics of lattice gas systems in contact with a wall. J. Chem. Phys., 108(7):3028–3037, 1998
work page 1998
-
[37]
H. Freist¨ uhler and M. Kotschote. Phase-field and Korteweg-type models for the time-dependent flow of compressible two-phase fluids. Archive for Rational Mechanics and Analysis , 224(1):1–20, 2017
work page 2017
-
[38]
T. Fukao and H. Wu. Separation property and convergence to equilibrium for the equation and dynamic boundary condition of Cahn–Hilliard type with singular potential. Asymptot. Anal. , 124(3-4):303–341, 2021
work page 2021
-
[39]
C. G. Gal. A Cahn–Hilliard model in bounded domains with permeable walls. Math. Methods Appl. Sci. , 29(17):2009–2036, 2006
work page 2009
-
[40]
C. G. Gal, M. Grasselli, and A. Miranville. Cahn–Hilliard–Navier–Stokes systems with moving contact lines. Calc. Var. Partial Differential Equations , 55(3):Art. 50, 47, 2016
work page 2016
-
[41]
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):1–56, 2019
work page 2019
-
[42]
C. G. Gal, M. 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. Preprint on ResearchGate: https://doi.org/10.13140/RG.2.2.27875.73764, 2024
-
[43]
H. Garcke and P. Knopf. Weak solutions of the Cahn–Hilliard system with dynamic boundary conditions: a gradient flow approach. SIAM J. Math. Anal. , 52(1):340–369, 2020
work page 2020
- [44]
- [45]
-
[46]
M.-H. Giga, A. Kirshtein, and C. Liu. Variational modeling and complex fluids. In Handbook of mathe- matical analysis in mechanics of viscous fluids , pages 73–113. Springer, Cham, 2018
work page 2018
-
[47]
G. Gilardi and J. Sprekels. Asymptotic limits and optimal control for the Cahn–Hilliard system with convection and dynamic boundary conditions. Nonlinear Anal., 178:1–31, 2019
work page 2019
- [48]
- [49]
-
[50]
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(3):Pa- per No. 65, 44, 2023
work page 2023
-
[51]
A. Giorgini, P. Knopf, and J. Stange. Regularity of a bulk-surface Cahn-Hilliard model driven by Leray velocity fields. Preprint: arXiv:2506.18617 [math.AP]
-
[52]
G. R. Goldstein, A. Miranville, and G. Schimperna. A Cahn–Hilliard model in a domain with non- permeable walls. Phys. D, 240(8):754–766, 2011
work page 2011
- [53]
-
[54]
P. Harder and B. Kov´ acs. Error estimates for the Cahn–Hilliard equation with dynamic boundary condi- tions. IMA J. Numer. Anal. , 42(3):2589–2620, 2022
work page 2022
- [55]
-
[56]
P. C. Hohenberg and B. I. Halperin. Theory of dynamic critical phenomena. Rev. Mod. Phys., 49:435–479, 1977
work page 1977
-
[57]
R. Kenzler, F. Eurich, P. Maass, B. Rinn, J. Schropp, E. Bohl, and W. Dietrich. Phase separation in confined geometries: Solving the Cahn–Hilliard equation with generic boundary conditions. Comp. Phys. Comm., 133:139–157, 2001
work page 2001
-
[58]
P. Knopf and K. F. Lam. Convergence of a Robin boundary approximation for a Cahn–Hilliard system with dynamic boundary conditions. Nonlinearity, 33(8):4191–4235, 2020
work page 2020
- [59]
-
[60]
P. Knopf and A. Signori. On the nonlocal Cahn–Hilliard equation with nonlocal dynamic boundary condition and boundary penalization. J. Differential Equations , 280:236–291, 2021
work page 2021
-
[61]
P. Knopf and J. Stange. Well-posedness of a bulk-surface convective Cahn–Hilliard system with dynamic boundary conditions. NoDEA Nonlinear Differential Equations Appl. , 31(5):Paper No. 82, 48, 2024
work page 2024
-
[62]
P. Knopf and J. Stange. Strong well-posedness and separation properties for a bulk-surface convective Cahn–Hilliard system with singular potentials. J. Differential Equations , 439:113408, 2025
work page 2025
- [63]
-
[64]
I.-S. Liu. Continuum mechanics. Advanced Texts in Physics. Springer-Verlag, Berlin, 2002
work page 2002
-
[65]
J. Lowengrub and L. Truskinovsky. Quasi-incompressible Cahn–Hilliard fluids and topological transitions. R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. , 454(1978):2617–2654, 1998
work page 1978
- [66]
-
[67]
X. Meng, X. Bao, and Z. Zhang. Second order stabilized semi-implicit scheme for the Cahn–Hilliard model with dynamic boundary conditions. J. Comput. Appl. Math. , 428:Paper No. 115145, 22, 2023
work page 2023
-
[68]
S. Metzger. An efficient and convergent finite element scheme for Cahn–Hilliard equations with dynamic boundary conditions. SIAM J. Numer. Anal. , 59(1):219–248, 2021
work page 2021
-
[69]
S. Metzger. A convergent SAV scheme for Cahn–Hilliard equations with dynamic boundary conditions. IMA Journal of Numerical Analysis , page drac078, 01 2023
work page 2023
-
[70]
A. Miranville. The Cahn–Hilliard equation: Recent advances and applications , volume 95 of CBMS-NSF Regional Conference Series in Applied Mathematics . Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2019
work page 2019
-
[71]
A. Miranville and H. Wu. Long-time behavior of the Cahn–Hilliard equation with dynamic boundary condition. J. Elliptic Parabol. Equ. , 6(1):283–309, 2020
work page 2020
-
[72]
J. Pr¨ uss and G. Simonett.Moving interfaces and quasilinear parabolic evolution equations , volume 105 of Monographs in Mathematics. Birkh¨ auser/Springer, [Cham], 2016
work page 2016
-
[73]
T. Qian, X.-P. Wang, and P. Sheng. A variational approach to moving contact line hydrodynamics. J. Fluid Mech., 564:333–360, 2006
work page 2006
- [74]
-
[75]
J. Shen, X. Yang, and Q. Wang. Mass and volume conservation in phase field models for binary fluids. Communications in Computational Physics , 13(4):1045–1065, 2013
work page 2013
-
[76]
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 Methods Appl. Sci. , 28(4):733–770, 2018
work page 2018
-
[77]
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(1):175– 221, 2023
work page 2023
-
[78]
H. Wu. A review on the Cahn–Hilliard equation: classical results and recent advances in dynamic boundary conditions. Electron. Res. Arch., 30(8):2788–2832, 2022
work page 2022
-
[79]
J. Yao, Z. Zhang, and W. Ren. Modelling moving contact lines on inextensible elastic sheets in two dimensions. J. Fluid Mech. , 955:Paper No. A25, 29, 2023. 23
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.