🔢 mathematics

A quasi-incompressible Cahn-Hilliard-Darcy model for two immiscible fluids in porous media

This paper derives a quasi-incompressible Cahn-Hilliard-Darcy model with logarithmic free energy for two-phase flows in porous media, demonstrates its convergence to the classical Muskat problem, and proves the global existence of weak solutions in both two and three dimensions.

Original authors: Daozhi Han, Sayantan Sarkar

Published 2026-07-07

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Daozhi Han, Sayantan Sarkar

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 watching two different liquids, like oil and water, trying to mix inside a sponge. In the real world, they don't snap apart instantly like a sharp line; instead, they have a fuzzy, blurry edge where they slowly blend into one another before separating again.

This paper is about creating a new mathematical "rulebook" to describe exactly how these two fluids move and mix inside that sponge, especially when the fluids are slightly squishy (compressible) rather than perfectly rigid.

Here is a breakdown of the paper's ideas using everyday analogies:

1. The Problem: The "Sharp Edge" is Too Fragile

Historically, scientists used a model called the "Muskat problem" to describe these fluids. Imagine drawing a razor-sharp line on a piece of paper to separate the oil from the water.

The Issue: In the real world, if that line gets too wiggly or unstable, the math breaks down. It's like trying to balance a pencil on its tip; eventually, it falls, and the math says the result is "infinity" or "undefined." This happens in real life too (like when a fluid interface collapses), but the old math couldn't handle it gracefully.

2. The Solution: The "Fuzzy Edge" (Phase-Field)

The authors propose a new model called qCHD (quasi-incompressible Cahn–Hilliard–Darcy).

The Analogy: Instead of a razor-sharp line, imagine the boundary between the fluids is a fuzzy, soft transition zone, like a gradient of color from blue to red.
Why it helps: This "fuzziness" acts like a safety net. Even if the fluids get chaotic, the math stays smooth and doesn't break. It treats the interface as a thin, but real, layer rather than a mathematical ghost line.

3. The Ingredients of the New Rulebook

The authors built this model by combining three main concepts:

The "Fuzzy" Rule (Cahn-Hilliard): This describes how the fluids want to separate. They use a specific "energy map" (called the Flory-Huggins potential) that acts like a valley with two deep pits. The fluids naturally roll down into these pits to stay pure, but the "fuzzy edge" keeps them from snapping apart too violently.
The "Sponge" Rule (Darcy's Law): This describes how fluids move through a porous material (like a sponge or soil). It's like water flowing through a coffee filter; the sponge slows the fluid down.
The "Squishy" Rule (Quasi-incompressible): Most old models assumed the fluids were perfectly hard to squeeze (incompressible). This new model admits that fluids can be slightly squishy. It accounts for the fact that when the fluids mix or separate, the total volume might shift just a tiny bit, which changes the pressure.

4. The "Energy Law" (The Battery)

One of the paper's biggest achievements is proving that this new system obeys a strict Energy Law.

The Analogy: Think of the system as a battery. The paper proves that as the fluids move and mix, the system can only lose energy (like a battery draining), never create it out of nowhere.
Why it matters: This "energy drain" is what keeps the math stable. It guarantees that the fluids won't suddenly start moving at infinite speeds or behaving in impossible ways.

5. The Big Result: "It Works!"

The authors spent a lot of time doing very advanced math to prove two things:

Existence: They proved that a solution to their new equations actually exists. In other words, the math doesn't break; you can always find an answer for how the fluids will behave.
Universality: They proved this works for both 2D (flat) and 3D (real-world) scenarios.

They also showed that if you make the "fuzzy edge" extremely thin (approaching zero), their new model turns back into the old, classic "Muskat" problem. This proves their new model is a more robust, upgraded version of the old one.

Summary

In short, the authors built a mathematical simulation engine for two fluids mixing in a sponge.

Old way: Used a sharp line that often broke the math.
New way: Uses a fuzzy, slightly squishy model that includes a "safety net" (energy dissipation).
The Win: They proved mathematically that this new way always works and never crashes, making it a reliable tool for understanding complex fluid flows in things like soil, oil reservoirs, or fuel cells.

Technical Summary: A Quasi-Incompressible Cahn–Hilliard–Darcy Model for Two Immiscible Fluids in Porous Media

Problem Statement
The study addresses the mathematical modeling of two-phase flows of immiscible fluids within porous media. While the classical Muskat problem effectively describes fluid interfaces as infinitely thin discontinuities, it is prone to finite-time singularities, such as unbounded growth of smooth profiles and instability even in seemingly stable regimes. To overcome these limitations, diffuse-interface (phase-field) formulations are employed to replace sharp discontinuities with narrow transition layers, thereby regularizing the dynamics.

Existing literature has extensively covered the Cahn–Hilliard–Navier–Stokes (CHNS) system for matched and unmatched densities, as well as Cahn–Hilliard–Hele–Shaw/Darcy (CHHS/CHD) models. However, a complete theoretical framework for a quasi-incompressible Cahn–Hilliard–Darcy (qCHD) model that incorporates a singular logarithmic Flory–Huggins potential and non-uniform fluid properties (specifically differing densities) remains an unresolved issue. Most existing Darcy-type models either assume strict incompressibility ( $\nabla \cdot u = 0$ ) or utilize polynomial free energy potentials, which lack the physical realism of the logarithmic potential near phase boundaries.

Methodology
The authors derive the qCHD model and establish the global existence of weak solutions through the following methodological steps:

Variational Derivation:
- The model is derived using Onsager's extremum principle (equivalent to Rayleigh's least energy dissipation principle).
- The system utilizes the volume fraction difference $\phi$ as the order parameter and the mass-averaged mixture velocity $u$ , which is non-solenoidal ( $\nabla \cdot u \neq 0$ ) due to density differences.
- The total energy functional $E(\phi)$ includes the Flory–Huggins homogeneous free energy $F(\phi)$ , a gradient term representing interfacial energy, and a gravitational potential term.
- By introducing a Lagrange multiplier $p$ to relax the continuity constraint and minimizing the Rayleighian, the authors derive the coupled system of equations governing momentum balance (Darcy's law), mass conservation, and the evolution of the order parameter.
Model Formulation:
- The resulting qCHD system consists of:
  - A modified Darcy equation relating velocity, pressure, chemical potential, and gravity.
  - A continuity equation for the mixture density.
  - A Cahn–Hilliard type equation for the order parameter $\phi$ , driven by the chemical potential $\mu_p = \mu + \alpha p$ , where $\alpha$ is the Atwood number.
- The free energy density $F(\phi)$ is the singular logarithmic function: $F(\phi) = \frac{\theta_0}{2}[(1+\phi)\ln(1+\phi) + (1-\phi)\ln(1-\phi)] - \frac{\theta_c}{2}\phi^2 + C$ .
Analytical Framework for Existence:
- Time-Discretization: A time-marching scheme is constructed using a convex-concave splitting of the free energy potential ( $F = F_v + F_n$ ).
- Fixed Point Argument: The existence of solutions to the time-discrete system is proven using the Browder-Minty lemma, relying on the strict monotonicity and coercivity of the associated operator.
- A Priori Estimates: The authors derive uniform estimates for the approximate solutions, leveraging the energy dissipation law and specific estimates on the pressure derived from Darcy's equations.
- Compactness: The Aubin-Lions-Simon lemma is employed to establish strong convergence of the order parameter, allowing the passage to the limit in the nonlinear terms, particularly the singular logarithmic potential.

Key Contributions and Results

Derivation of the qCHD Model: The paper presents a systematic derivation of a quasi-incompressible Cahn–Hilliard–Darcy model that satisfies local mass conservation and an energy-dissipation law. The model accounts for density differences via the Atwood number and uses a mass-averaged velocity.
Global Existence of Weak Solutions: The primary mathematical result is the proof of the global existence of weak solutions for the qCHD system in both 2D and 3D domains. This is achieved for the physically realistic Flory–Huggins logarithmic potential, which ensures the order parameter remains strictly within the physical bounds $(-1, 1)$ .
Sharp Interface Limit: Through formal matched asymptotic expansion, the authors demonstrate that the qCHD model converges to a modified Muskat problem in the sharp interface limit, confirming the model's consistency with classical interfacial dynamics.
Energy Stability: The derived system satisfies a rigorous energy-dissipation law, ensuring that the total energy decreases over time due to viscous dissipation and chemical diffusion.

Significance
The paper claims to provide the first study presenting a quasi-incompressible Cahn–Hilliard–Darcy system that simultaneously incorporates a singular logarithmic potential and non-uniform fluid properties (unequal densities).

The significance of these results is highlighted by the fact that global weak existence for quasi-incompressible Cahn–Hilliard–Navier–Stokes equations with such potentials is not currently known. By establishing well-posedness for the Darcy analogue, this work fills a critical gap in the mathematical theory of two-phase flows in porous media, offering a robust framework for applications ranging from unsaturated soil flow to geological CO2 storage and fuel cell water management. The use of the logarithmic potential ensures the model remains physically valid near pure phases, avoiding the unphysical behavior that can arise with polynomial approximations.