🔢 mathematics

Local-in-time existence of strong solutions to a quasi-incompressible Cahn--Hilliard--Navier--Stokes system

This paper establishes the local existence and uniqueness of strong solutions for a quasi-incompressible Cahn–Hilliard–Navier–Stokes system modeling two-phase flows with unequal densities using the Banach fixed point theorem and maximal regularity theory.

Original authors: Mingwen Fei, Xiang Fei, Daozhi Han, Yadong Liu

Published 2026-04-28

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Mingwen Fei, Xiang Fei, Daozhi Han, Yadong Liu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to model how two different liquids—like oil and water—mix and move inside a container. This isn't just a simple stir; it’s a complex dance where the liquids have different weights (densities), they push against each other, and they are constantly trying to separate or blend at their boundary.

This mathematical paper tackles a very specific, high-level version of this problem. Here is the breakdown of what they did, using everyday analogies.

1. The Problem: The "Heavy vs. Light" Dance

Most scientific models assume that the two liquids are roughly the same weight. But in the real world, oil is much lighter than water. This difference in weight creates a "tug-of-war" effect.

The researchers are looking at a "Quasi-Incompressible" system.

The Analogy: Imagine a crowded subway car. If everyone is the same size, it’s easy to predict how the crowd moves. But if you suddenly mix giant sumo wrestlers with small children, the way the crowd shifts becomes much more chaotic. The "pressure" doesn't just push people around; it changes how the "density" of the crowd shifts. In this paper, the math has to account for the fact that the pressure and the chemical makeup of the liquids are "talking" to each other constantly.

2. The "Order Parameter": The Boundary Line

In these models, scientists use something called an "order parameter" ( $\phi$ ).

The Analogy: Think of a coloring book. Instead of a sharp black line separating the blue section from the red section, imagine a "gradient" where the colors bleed into each other. The order parameter is the mathematical way of saying, "At this exact spot, it is 70% oil and 30% water." The paper focuses on how this "blurry line" moves and evolves over time.

3. The Challenge: The "Mathematical Friction"

The authors mention a "higher-order term" that makes this much harder than previous models.

The Analogy: Imagine you are trying to write a poem while riding a unicycle on a tightrope. Previous scientists had already figured out how to write the poem on solid ground. But this new model adds a "wobble" (the higher-order term) that makes the pen slip. If the pen slips too much, the poem becomes gibberish (the math breaks).

To fix this, the researchers had to add a "smallness assumption." They basically said, "We can solve this beautiful, complex poem, provided the initial wobble isn't too violent."

4. The Solution: The "Fixed Point" Strategy

How do you solve an equation where everything depends on everything else? (The velocity depends on the density, which depends on the pressure, which depends on the velocity...)

They used something called the Banach Fixed Point Theorem.

The Analogy: Imagine you are trying to find the exact center of a spinning merry-go-round. You can't just jump to the center. Instead, you make a guess, see how far off you are, adjust your position slightly, and repeat. If your adjustments are smart enough, you will eventually "settle" into the perfect center. The researchers proved that their mathematical "adjustments" would always settle into a single, unique, and stable answer.

Summary: What did they actually achieve?

They proved that for a very complex, realistic model of two liquids with different weights, a solution actually exists and is unique for at least a short period of time.

They provided the "mathematical permission" for other scientists and engineers to use these equations in computer simulations to predict how chemicals, fuels, or biological fluids will behave in the real world. They proved that the math isn't just a fantasy—it's a stable, reliable map of reality.

Technical Summary: Local-in-Time Existence of Strong Solutions to a Quasi-Incompressible Cahn–Hilliard–Navier–Stokes System

1. Problem Statement

The paper investigates the mathematical well-posedness of a quasi-incompressible Cahn–Hilliard–Navier–Stokes (qCHNS) system. This system is a diffuse interface model designed to describe two-phase flows where the two fluids have mismatched (unmatched) densities.

Unlike standard incompressible models (where the velocity field is divergence-free), this model utilizes a mass-averaged (barycentric) velocity. Because the densities of the two phases differ, the mass-averaged velocity is not divergence-free within the diffusive interface, leading to a "quasi-incompressible" nature. A critical feature of this model is the strong coupling between the fluid mechanics and the phase separation: the pressure $p$ enters the equation for the chemical potential $\mu$ , which in turn influences the evolution of the order parameter $\phi$ (the volume fraction difference).

The authors consider the system in a smooth bounded domain $\Omega \subset \mathbb{R}^3$ with Navier-slip boundary conditions for the velocity and specific boundary conditions for the order parameter and chemical potential.

2. Methodology

The authors employ a rigorous analytical framework based on maximal regularity theory and the Banach fixed-point theorem. The proof strategy is divided into several sophisticated steps:

Reformulation and Linearization: To handle the complex coupling, the authors eliminate the pressure $p$ and chemical potential $\mu$ by expressing them in terms of the velocity $u$ and order parameter $\phi$ using the weak Neumann-Laplace operator. The system is then rewritten as a nonlinear evolution equation of the form $L(\phi) \begin{pmatrix} u \\ \phi' \end{pmatrix} = F(u, \phi)$ , where $L(\phi)$ is a linear operator and $F$ contains the nonlinearities.
Maximal Regularity and Semigroup Theory: The authors prove that the linearized operator $L(\phi_0)$ is an isomorphism between specific anisotropic fractional Sobolev spaces (Hilbert spaces). This is achieved by decomposing the velocity field using the Helmholtz decomposition and applying the theory of analytic semigroups to the resulting Stokes-like and Cahn-Hilliard-like operators.
Contraction Mapping Principle: By establishing that the nonlinear terms $F(u, \phi)$ are Lipschitz continuous with respect to the solution norms (with a constant that vanishes as the time $T \to 0$ ), the authors apply the Banach fixed-point theorem to guarantee a unique local-in-time solution.
Handling Higher-Order Terms: A significant technical challenge identified is the presence of a higher-order term $\text{div}(\phi u)$ in the order parameter equation. To overcome the difficulty this poses to the standard analysis framework, the authors impose a smallness assumption on the initial data $\phi_0$ (ensuring $|\phi_0|$ is sufficiently small), which keeps the density $\rho$ bounded away from zero.

3. Key Contributions

Mathematical Treatment of Quasi-Incompressibility: The paper provides a rigorous treatment of the pressure-dependent chemical potential, a feature that distinguishes this model from the more common "matched density" Model H.
Advanced Boundary Condition Analysis: The authors successfully incorporate Navier-slip boundary conditions into the existence proof, utilizing complex trace estimates in Besov and fractional Sobolev spaces.
Overcoming Nonlinearity Challenges: The paper identifies and solves the analytical difficulty posed by the $\text{div}(\phi u)$ term, which is more singular than the terms found in previous models (like the Lowengrub–Truskinovsky model).
Unified Framework: The work contributes to the broader effort of unifying various diffuse interface models by providing a rigorous foundation for the mass-averaged velocity approach.

4. Main Results

The central result is Theorem 1.1, which states:
Given initial data $u_0 \in H^1(\Omega)$ and $\phi_0 \in H^2(\Omega)$ (with $|\phi_0|$ sufficiently small), there exists a unique strong solution $(u, \phi)$ on a short time interval $[0, T]$ .

The regularity of the solution is established as:

$u \in H^1(0, T; L^2(\Omega)) \cap L^2(0, T; H^2(\Omega))$
$\phi \in L^2(0, T; H^3(\Omega)) \cap H^1(0, T; H^1(\Omega))$
The solution maintains physical bounds: $|\phi| \leq 1/2$ and $\rho \geq 1/4$ .

5. Significance

This research is significant for both the mathematical and computational fluid dynamics communities.

Mathematically, it advances the theory of non-homogeneous Navier-Stokes equations coupled with phase-field models, pushing the boundaries of what can be proven regarding regularity and well-posedness in the presence of density mismatches.
Physically/Computationally, the qCHNS model is highly relevant for simulating real-world two-phase flows (such as oil and water) where density differences are non-negligible. The authors note that their specific formulation (using volume fraction difference) is particularly well-suited for numerical implementation, making this theoretical result a vital foundation for stable and accurate numerical simulations.