Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions

Basic Information

• Paper ID: 2511.06847
• Title: Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions
• Author: Jonas Stange (Universität Regensburg, Germany)
• Classification: math.AP (Analysis of PDEs), math-ph (Mathematical Physics), math.MP (Mathematical Physics)
• Submission Date: November 10, 2025
• Paper Link: https://arxiv.org/abs/2511.06847v1

Abstract

This paper investigates a thermodynamically consistent diffuse interface model for describing bulk-surface viscous fluid mixtures. The model couples the Navier-Stokes-Cahn-Hilliard equations in the bulk domain with a surface Navier-Stokes-Cahn-Hilliard system on the boundary. The author proves global well-posedness of strong solutions in the two-dimensional case, covering physically meaningful non-degenerate mobility functions. Furthermore, uniqueness of strong solutions and continuous dependence on initial data are established. The approach is based on novel well-posedness and regularity theory for the convective bulk-surface Cahn-Hilliard equation with non-degenerate mobilities, as well as regularity theory for bulk-surface Stokes equations with non-constant coefficients.

Research Background and Motivation

Problem Background

1. Physical Phenomenon Modeling: The research involves diffuse interface models for two-phase fluid mixtures with different densities, where viscous dynamics occur both in the bulk domain and on the boundary. Such models have important applications in materials science, biomembrane dynamics, and related fields.
2. Mathematical Challenges: The model couples:
   • Navier-Stokes-Cahn-Hilliard equations in the bulk domain (describing fluid motion and phase separation)
   • Surface Navier-Stokes-Cahn-Hilliard system on the boundary (describing surface dynamics)
   • Dynamic boundary conditions (allowing mass exchange between bulk and surface)
3. Limitations of Existing Theory:
   • For diffuse interface models with dynamic boundary conditions, existing literature primarily focuses on weak solution existence
   • Strong solution existence is less studied, particularly for non-degenerate mobility functions
   • Global strong solution theory in three dimensions remains incomplete

Research Motivation

1. Physical Reasonableness: Non-degenerate mobility functions (i.e., mobilities with positive upper and lower bounds throughout the domain) are more physically realistic
2. Mathematical Rigor: Strong solutions provide finer regularity estimates, which are crucial for understanding long-time behavior of solutions
3. Theoretical Completeness: Filling the gap in strong solution theory for Navier-Stokes-Cahn-Hilliard systems with dynamic boundary conditions

Core Contributions

The main contributions of this paper include:

1. Global Strong Solution Existence: In the two-dimensional case, the author proves existence of global strong solutions for the bulk-surface Navier-Stokes-Cahn-Hilliard system with non-degenerate mobilities (Theorem 3.1)
2. Uniqueness and Continuous Dependence: Uniqueness of strong solutions and continuous dependence on initial data are established (Theorem 3.4), subject to the restriction $L \in (0,\infty]$
3. Novel Regularity Theory: Two key auxiliary theories are established:
   • Well-posedness and regularity theory for the convective bulk-surface Cahn-Hilliard equation with non-degenerate mobilities (Appendix A)
   • Regularity theory for bulk-surface Stokes equations with non-constant coefficients (Appendix B)
4. Energy Inequality: Solutions satisfy the energy dissipation law: $\frac{d}{dt}E_{tot}(v,w,\phi,\psi) + \text{dissipative terms} = 0$
5. Separation Property: Phase field functions strictly remain in the physical interval $(-1,1)$: $|\phi| < 1$ a.e. in $Q$, $|\psi| < 1$ a.e. on $\Sigma$

Methodology Details

Problem Formulation

Consider a two-dimensional bounded domain $\Omega \subset \mathbb{R}^2$ with boundary $\Gamma = \partial\Omega$. Set $Q = \Omega \times (0,\infty)$, $\Sigma = \Gamma \times (0,\infty)$.

Governing System (Equation 1.1) includes:

1. Bulk Momentum Equation: $\partial_t(\rho(\phi)v) + \text{div}(v \otimes (\rho(\phi)v + J)) = \text{div} T, \quad \text{div} v = 0 \text{ in } Q$
2. Surface Momentum Equation: $\partial_t(\sigma(\psi)w) + \text{div}_\Gamma(w \otimes (\sigma(\psi)w + K)) = \text{div}_\Gamma T_\Gamma + Z, \quad \text{div}_\Gamma w = 0 \text{ on } \Sigma$
3. Bulk Cahn-Hilliard Equation: $\partial_t\phi + \text{div}(\phi v) = \text{div}(m_\Omega(\phi)\nabla\mu), \quad \mu = -\Delta\phi + F'(\phi) \text{ in } Q$
4. Surface Cahn-Hilliard Equation: $\partial_t\psi + \text{div}_\Gamma(\psi w) = \text{div}_\Gamma(m_\Gamma(\psi)\nabla_\Gamma\theta) - \beta m_\Omega(\phi)\partial_n\mu \text{ on } \Sigma$$\theta = -\Delta_\Gamma\psi + G'(\psi) + \alpha\partial_n\phi \text{ on } \Sigma$
5. Boundary Conditions: $w = v|_\Gamma, \quad v \cdot n = 0 \text{ on } \Sigma$$K\partial_n\phi = \alpha\psi - \phi, \quad Lm_\Omega(\phi)\partial_n\mu = \beta\theta - \mu \text{ on } \Sigma$

Where:

• $v, w$ are bulk and surface velocity fields
• $\phi, \psi$ are bulk and surface phase field variables
• $\mu, \theta$ are chemical potentials
• $\rho(\phi), \sigma(\psi)$ are phase-dependent density functions
• $m_\Omega, m_\Gamma$ are mobility functions (satisfying $0 < m_* \leq m_\Omega, m_\Gamma \leq m^*$)

Model Architecture

Proof Strategy employs a semi-Galerkin approximation method:

Step One: Approximate Problem Construction

1. Introduce eigenfunctions $\{(\tilde{v}_j, \tilde{w}_j)\}_{j\in\mathbb{N}}$ of the bulk-surface Stokes operator $A$
2. Define finite-dimensional subspace $V_k = \text{span}\{(\tilde{v}_1, \tilde{w}_1), \ldots, (\tilde{v}_k, \tilde{w}_k)\}$
3. Construct approximate solutions $(v_k, w_k, \phi_k, \psi_k, \mu_k, \theta_k)$ satisfying Galerkin weak formulations

Step Two: Fixed Point Argument

For fixed velocity field $(v_*, w_*)$, solve:

• Convective bulk-surface Cahn-Hilliard equation (4.8), obtaining $(\phi_k, \psi_k, \mu_k, \theta_k)$
• Galerkin approximation of Navier-Stokes equation (4.12), obtaining $(v_k, w_k)$

By Schauder's fixed point theorem, the mapping $\Lambda: (v_*, w_*) \mapsto (v_k, w_k)$ admits a fixed point.

Step Three: A Priori Estimates

Establish key uniform estimates (independent of $k$):

1. Energy Estimates (Equations 4.35-4.37): $\|(v_k, w_k)\|_{L^\infty(0,\infty; L^2)} + \|(v_k, w_k)\|_{L^2(0,\infty; H^1)} \leq C$$\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; H^1)} \leq C$
2. Higher-Order Estimates (Equations 4.43-4.45): $\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; W^{2,p})} + \|(F'(\phi_k), G'(\psi_k))\|_{L^\infty(0,\infty; L^p)} \leq K_{1,p}$$\|(\mu_k, \theta_k)\|_{L^\infty(0,\infty; H^1)} + \|(\mu_k, \theta_k)\|_{L^2_{uloc}(0,\infty; H^3)} \leq K_2$
3. Velocity Field Higher-Order Estimates (Equation 4.61): $\sup_{t\geq 0} \int_t^{t+1} \|( v_k, w_k)\|_{H^2}^2 + \|(\partial_t v_k, \partial_t w_k)\|_{L^2}^2 + \|(p_k, q_k)\|_{H^1}^2 \, ds \leq C$

Step Four: Limiting Process

Using compactness arguments (Aubin-Lions lemma) and monotonicity techniques, extract convergent subsequences as $k \to \infty$ to obtain strong solutions.

Technical Innovations

1. Handling Non-Degenerate Mobilities: Unlike the constant mobility case, non-degenerate mobilities introduce additional nonlinear terms. The key technique is using elliptic regularity theory (Equations 2.13-2.15) to control these terms.
2. Higher-Order Energy Estimates: By differentiating the energy functional (Equation 4.46) and combining with bulk-surface Stokes regularity theory (Theorem B.2), establish $H^2$ estimates for velocity fields: $\frac{d}{dt}\text{(energy)} + \text{dissipation} \leq \varpi\|(v_k, w_k)\|_{H^2}^2 + \text{lower-order terms}$ Choose appropriate $\varpi$ to absorb the $H^2$ norm.
3. Preservation of Separation Property: Using singularity of logarithmic potentials and variants of maximum principles, prove that phase field functions strictly remain in $(-1, 1)$.
4. Uniqueness Proof: For the difference $(v, w, \phi, \psi) = (v_1 - v_2, w_1 - w_2, \phi_1 - \phi_2, \psi_1 - \psi_2)$ of two solutions, establish differential inequalities (Equations 5.2-5.3): $\frac{d}{dt}\left(\|(v, w)\|_{L^2}^2 + \|(\phi, \psi)\|_{K,\alpha}^2\right) + \text{dissipation} \leq P(t) \times \text{(norm of difference)}$ where $P(t) \in L^1(0,T)$. Apply Gronwall's lemma to obtain uniqueness.

Experimental Setup

Note: This is a pure theoretical mathematics paper without numerical experiments. The main results are mathematical proofs of existence and uniqueness.

Theoretical Verification Framework

The paper verifies theory through the following mathematical tools:

1. Function Space Setup:
   • Velocity field space: $H^1_{0,\text{div}} = \{(v,w) \in H^1_0: \text{div} v = 0, \text{div}_\Gamma w = 0\}$
   • Phase field space: $V^1_{L,\beta}$ (depending on coupling parameters $L, \beta$)
   • Chemical potential space: $H^1_{L,\beta}$
2. Regularity Requirements:
   • Initial data: $(v_0, w_0) \in H^1_{\text{div}}$, $(\phi_0, \psi_0) \in H^1$ satisfying $\|\phi_0\|_{L^\infty} \leq 1$, $\|\psi_0\|_{L^\infty} \leq 1$
   • Compatibility conditions: There exist $(\mu_0, \theta_0) \in H^1_{L,\beta}$ satisfying specific integral relations
3. Parameter Restrictions:
   • Spatial dimension: $d = 2$
   • Coupling parameters: $K \in (0, \infty)$, $L \in [0, \infty]$ (uniqueness requires $L \in (0, \infty]$)
   • Mobilities: $0 < m_* \leq m_\Omega, m_\Gamma \leq m^* < \infty$

Key Assumptions

(A1-A7) include:

• Domain regularity: $\Omega \subset \mathbb{R}^2$ with $C^3$ boundary
• Potential functions: $F(s) = F_0(s) - \frac{c_F}{2}s^2$, where $F_0$ is strongly convex and singular as $s \to \pm 1$
• Growth conditions: $F_0''(s) \geq \Theta > 0$ or specific growth rate restrictions (A7.1)-(A7.2)

Experimental Results

Main Theoretical Results

Theorem 3.1 (Global Strong Solution Existence): Under assumptions (A1)-(A7), for $K \in (0, \infty)$, $L \in [0, \infty]$, there exists a global strong solution $(v, w, \phi, \psi, \mu, \theta)$ satisfying:

1. Regularity:
   • $(v, w) \in BC([0,\infty); H^1_{0,\text{div}}) \cap L^2_{uloc}([0,\infty); H^2) \cap H^1_{uloc}(0,\infty; L^2_{\text{div}})$
   • $(\phi, \psi) \in L^\infty(0,\infty; H^3)$
   • $(F'(\phi), G'(\psi)), (F''(\phi), G''(\psi)) \in L^\infty(0,\infty; L^p)$ for all $2 \leq p < \infty$
2. Separation Property: $|\phi| < 1$ a.e. in $Q$, $|\psi| < 1$ a.e. on $\Sigma$
3. Energy Inequality: $E_{tot}(v(t), w(t), \phi(t), \psi(t)) + \int_0^t \text{(dissipation)} \, ds = E_{tot}(v_0, w_0, \phi_0, \psi_0)$

Theorem 3.4 (Uniqueness): Under the assumptions of Theorem 3.1, if $L \in (0, \infty]$, then the strong solution is unique.

Key Estimates

The paper establishes the following quantitative estimates (Equations A.10-A.12):

$\|(\phi, \psi)\|_{L^\infty(0,\infty; H^1)}^2 + \|(\partial_t\phi, \partial_t\psi)\|_{L^2(0,\infty; (H^1_{L,\beta})')}^2 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^2_{uloc}([0,\infty); W^{2,p})}^2 + \|(F'(\phi), G'(\psi))\|_{L^2_{uloc}([0,\infty); L^p)}^2 \leq C_p(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^4_{uloc}([0,\infty); H^2)}^4 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)^2$

Overcoming Technical Difficulties

1. Limitation of the $L = 0$ Case: When $L = 0$, boundary condition (1.1h)$_2$ does not allow reformulation of the normal derivative $\partial_n\mu$, causing the uniqueness proof to fail. The author notes this requires additional structural assumptions (Remark 3.5).
2. Challenges in Three Dimensions: The paper only addresses the two-dimensional case because in three dimensions:
   • The separation property is lost
   • More refined approximation techniques are needed to control nonlinear terms
   • Only local strong solutions rather than global ones can be expected

Related Work

Main Research Directions

1. Models without Dynamic Boundary Conditions:
   • Abels-Garcke-Grün (AGG) model: 2, 4, 16, 17, 21, 22
   • Mature strong solution theory already exists
2. Models with Dynamic Boundary Conditions:
   • Primarily focus on weak solution existence: 9, 14, 15, 19, 30
   • Quasi-strong solutions: Allen-Cahn-Navier-Stokes-Voigt system 13
3. Regularity of Cahn-Hilliard Equations:
   • Constant mobility case: 20, 28, 29
   • Non-degenerate mobilities: 11, 39 (prior work by the author)

Novelty of This Work

Compared to existing work, this paper:

• First addresses strong solutions for Navier-Stokes-Cahn-Hilliard systems with dynamic boundary conditions
• First covers the non-degenerate mobility case
• Establishes novel bulk-surface Stokes equation regularity theory (Appendix B)
• Develops higher-order estimate techniques for convective Cahn-Hilliard equations (Appendix A)

Conclusions and Discussion

Main Conclusions

1. In the two-dimensional case, the existence and uniqueness of global strong solutions for bulk-surface Navier-Stokes-Cahn-Hilliard systems with non-degenerate mobilities are proved (when $L \in (0, \infty]$)
2. Solutions satisfy strict separation properties and energy dissipation laws
3. Two important auxiliary theoretical frameworks are established, which can be independently applied to other problems

Limitations

1. Dimensional Restriction: Results apply only to the two-dimensional case. The three-dimensional case requires:
   • More refined nonlinear estimates
   • Possibly only local strong solutions
   • For constant mobility case, weak solution uniqueness and strong solution existence are known
2. Parameter Restrictions:
   • $K = 0$ is not feasible (non-degenerate mobility case)
   • Uniqueness for $L = 0$ remains unresolved without additional structural assumptions
3. Compatibility Conditions: When $L = 0$, density functions must satisfy the special relation $\beta(\sigma_2 - \sigma_1) = \rho_2 - \rho_1$ (Equation 3.2)

Future Directions

1. Three-Dimensional Extension: Develop techniques applicable to three dimensions, possibly requiring:
   • Local existence theory
   • Global results for small initial data or special structures
   • Constant mobility assumptions
2. The $L = 0$ Case: Seek new methods to handle this case or determine its infeasibility
3. Long-Time Behavior: Study asymptotic behavior of solutions and existence of attractors
4. Numerical Verification: Develop efficient numerical schemes to verify theoretical predictions

In-Depth Evaluation

Strengths

1. Theoretical Rigor:
   • Complete and detailed proofs with standard technical treatment
   • Clear explanations of various technical difficulties
   • Precise statement of lemmas and theorems
2. Novelty:
   • First strong solution theory for this class of systems
   • Non-degenerate mobility handling is challenging
   • Auxiliary theories in appendices have independent value
3. Clear Structure:
   • Main results clearly stated in Section 3
   • Proof strategy presented hierarchically
   • Technical details in appendices without affecting main narrative
4. Physical Reasonableness:
   • Thermodynamically consistent model
   • Non-degenerate mobilities more realistic
   • Energy dissipation and mass conservation guaranteed

Weaknesses

1. Dimensional Limitation: Two-dimensional results have limited practical application; three-dimensional case is more important but unresolved
2. Parameter Restrictions: Restrictions $K \in (0, \infty)$ and $L \in (0, \infty]$ exclude certain physical scenarios
3. Lack of Numerical Verification: As a pure theoretical work, no numerical examples support the theoretical results
4. Technical Complexity: Highly technical proofs requiring extensive background knowledge, potentially limiting readership
5. Constant Dependence: Many estimates involve constants $C$ depending on multiple parameters without explicit expressions

Impact

1. Theoretical Contribution:
   • Fills important gap in strong solution theory for Navier-Stokes-Cahn-Hilliard systems with dynamic boundary conditions
   • Provides technical tools and theoretical framework for subsequent research
2. Methodological Value:
   • Semi-Galerkin method combined with fixed point arguments applicable to other coupled systems
   • Higher-order estimate techniques are general
3. Application Potential:
   • Provides theoretical support for numerical methods
   • Applicable to modeling in materials science, biomembrane dynamics, and related fields

Applicable Scenarios

1. Materials Science: Phase separation and interfacial dynamics of two-phase fluid materials
2. Biophysics: Fluid mosaic model of cell membranes, lipid bilayer dynamics
3. Contact Line Dynamics: Motion and wetting of droplets on solid surfaces
4. Theoretical Research:
   • Regularity theory for nonlinear PDEs
   • Well-posedness analysis of coupled systems
   • Mathematical theory of dynamic boundary conditions

Selected References

2 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., 2024.

11 M. Conti, P. Galimberti, S. Gatti, and A. Giorgini. New results for the Cahn-Hilliard equation with non-degenerate mobility. Calc. Var. PDE, 2025.

19 A. Giorgini and P. Knopf. Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models. J. Math. Fluid Mech., 2023.

30 P. Knopf and J. Stange. A thermodynamically consistent model for bulk-surface viscous fluid mixtures. Preprint, 2025.

39 J. Stange. Well-posedness and long-time behavior of a bulk-surface Cahn-Hilliard model with non-degenerate mobility. Preprint, 2025.

Overall Assessment: This is a high-quality theoretical mathematics paper making significant contributions to PDE analysis. Despite dimensional and parameter restrictions, its technical innovations and theoretical depth are commendable. For mathematicians and applied scientists working in related fields, this paper provides valuable theoretical tools and research paradigms.