Brinkman–Cahn–Hilliard Flow Models

• Brinkman–Cahn–Hilliard systems are diffuse-interface models that couple phase separation dynamics with Brinkman-type porous flow, applicable to binary mixtures and tumor growth.
• They employ variational formulations and energy dissipation laws to ensure mass conservation, well-posedness, and convergence of solutions.
• Advanced variants integrate nonlocal interactions, curvature effects, and optimal control, enabling robust numerical schemes and modeling of complex multiphase phenomena.

Searching arXiv for recent and foundational Brinkman–Cahn–Hilliard papers to ground the article. The Brinkman–Cahn–Hilliard system is a class of diffuse-interface partial differential equation models coupling a Cahn–Hilliard-type phase evolution for an order parameter with a Brinkman-type momentum balance for the velocity field. In these models, the phase variable describes composition or volume-fraction contrast in a binary mixture or related multiphase medium, while the Brinkman equation accounts for viscous flow in porous or partially porous environments through a combination of viscous stress and linear drag. The coupling is bidirectional: the flow advects the phase field, and the phase field feeds back into the momentum equation through a capillary or Korteweg force such as $\mu\nabla\phi$ or $(\mu+\chi\sigma)\nabla\varphi$. Across the literature, Brinkman–Cahn–Hilliard systems appear in incompressible binary-fluid flow, porous-medium phase separation, tumor-growth models, nonlocal interaction models, higher-order curvature-driven phase separation, optimal control, stochastic extensions, and vanishing-viscosity limits toward Darcy or Hele–Shaw regimes (Bosia et al., 2014, Porta et al., 2016, Ebenbeck et al., 2018, Colli et al., 24 Sep 2025, Brzeźniak et al., 10 Jan 2026).

1. Core mathematical structure

A standard local Brinkman–Cahn–Hilliard model couples

$$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$

with a Brinkman momentum balance

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$

or equivalent variants differing by normalization, coefficient dependence, and source terms (Bosia et al., 2014). In this formulation, $\phi$ is the order parameter, $\mu$ the chemical potential, $\mathbf u$ the averaged velocity, and $p$ the pressure. The term $\mu\nabla\phi$ is the Korteweg force, while $\eta\mathbf u$ is the porous-medium drag characteristic of Brinkman flow (Bosia et al., 2014, Chen et al., 8 Jun 2025).

This baseline structure admits several important generalizations. In nonlocal models, the chemical potential takes the form

$(\mu+\chi\sigma)\nabla\varphi$0

so the local gradient term is replaced by convolution with an interaction kernel $(\mu+\chi\sigma)\nabla\varphi$1 (Porta et al., 2016, Dharmatti et al., 2019, Dharmatti et al., 6 Feb 2025). In tumor-growth systems, additional fields such as nutrient concentration $(\mu+\chi\sigma)\nabla\varphi$2 enter the phase and momentum equations, and the velocity may be non-solenoidal because $(\mu+\chi\sigma)\nabla\varphi$3 is prescribed by solution-dependent source terms (Ebenbeck et al., 2018, Ebenbeck et al., 2019, Colli et al., 2022, Schimperna, 2024). In higher-order curvature models, the phase subsystem is reformulated using auxiliary variables and the chemical potential itself becomes fourth order in the phase field, yielding an overall sixth-order evolution (Colli et al., 24 Sep 2025, Bag, 8 Dec 2025, Colli et al., 9 Jan 2026).

The principal mechanism is always the same: advection by Brinkman flow modifies phase separation, while diffuse-interface thermodynamics modifies momentum transport. This makes the system an intermediate model between pure Cahn–Hilliard diffusion, Darcy-type porous-medium flow, and Stokes-like viscous flow (Bosia et al., 2014, Chen et al., 8 Jun 2025).

2. Variational and energetic formulation

A defining feature of Brinkman–Cahn–Hilliard systems is their energetic structure. For the classical local model, the free energy is

$(\mu+\chi\sigma)\nabla\varphi$4

and, under periodic or no-flux/incompressible settings, the continuous dissipation law takes the form

$(\mu+\chi\sigma)\nabla\varphi$5

in one common normalization (Chen et al., 8 Jun 2025). In the constant-coefficient CHB analysis, the corresponding identity is

$(\mu+\chi\sigma)\nabla\varphi$6

(Bosia et al., 2014).

For nonlocal systems, the energy functional is

$(\mu+\chi\sigma)\nabla\varphi$7

equivalently

$(\mu+\chi\sigma)\nabla\varphi$8

(Porta et al., 2016). This formulation is analytically advantageous because the nonlocal chemical potential remains second-order in space rather than fourth-order, which is central to the nonlocal CHHS limit theory (Porta et al., 2016).

In curvature-driven models, the energy is of Willmore/functionalized-Cahn–Hilliard type: $(\mu+\chi\sigma)\nabla\varphi$9 with

$$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$0

(Colli et al., 24 Sep 2025). This is the mechanism behind the “curvature-induced phase separation” terminology: the chemical potential is the first variation of a higher-order free energy, so the phase dynamics becomes sixth-order after elimination of $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$1 (Colli et al., 24 Sep 2025, Bag, 8 Dec 2025).

Source terms can break strict gradient-flow structure. In the higher-order Brinkman system with

$$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$2

the spatial mean of the phase is not conserved, and the mean of $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$3 must be estimated separately (Colli et al., 24 Sep 2025). In tumor models, source terms in both divergence and phase equations similarly destroy mass conservation and complicate control of chemical-potential averages (Ebenbeck et al., 2018, Ebenbeck et al., 2019).

3. Principal model variants

The Brinkman–Cahn–Hilliard literature is heterogeneous. The main variants differ in the order of the phase operator, local versus nonlocal interaction, compressibility, and auxiliary couplings such as nutrient or chemotaxis.

Variant	Characteristic feature	Representative papers
Classical local CHB	Fourth-order phase field, Brinkman drag and viscosity	(Bosia et al., 2014, Chen et al., 8 Jun 2025)
Nonlocal CHB	$$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$4	(Porta et al., 2016, Dharmatti et al., 2019, Dharmatti et al., 6 Feb 2025)
Tumor-growth CHB	Nutrient coupling, sources, often non-solenoidal velocity	(Ebenbeck et al., 2018, Ebenbeck et al., 2019, Colli et al., 2022, Schimperna, 2024)
Sixth-order Brinkman–CH	Curvature effects via $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$5	(Colli et al., 24 Sep 2025, Bag, 8 Dec 2025, Colli et al., 9 Jan 2026)
Stochastic CHB	Multiplicative Wiener noise, dynamic boundary conditions	(Brzeźniak et al., 10 Jan 2026)

The classical CHB system for incompressible binary flow in porous media treats $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$6 as divergence-free and typically imposes no-slip or traction-free velocity conditions together with Neumann conditions for $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$7 and $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$8 (Bosia et al., 2014, Porta et al., 2016). Nonlocal CHB replaces the local interfacial term with long-range interaction through $$\partial_t \phi + \nabla\cdot(\phi \mathbf u)=\Delta \mu,\qquad \mu=-\Delta\phi+f(\phi),$$9, which leads to stronger well-posedness results for the corresponding Hele–Shaw limit than are known in the local case (Porta et al., 2016).

Tumor-growth variants add nutrient concentration $$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$0, mass-production terms, chemotaxis, and compressibility. One model uses

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$1

with the phase equation containing both convection and source terms, and a quasi-static elliptic nutrient equation (Ebenbeck et al., 2018, Ebenbeck et al., 2019). In this setting, traction-free or frictionless velocity conditions are analytically important because Dirichlet conditions would impose the restrictive compatibility condition $$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$2 (Ebenbeck et al., 2018, Ebenbeck et al., 2019).

The higher-order sixth-order system studied in three dimensions introduces the unknowns $$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$3 and the coupled equations

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$4

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$5

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$6

(Colli et al., 24 Sep 2025). This formulation makes the sixth-order structure analytically tractable by reducing it to a mixed second-order system.

4. Existence, uniqueness, and regularity theory

The foundational local CHB analysis established global weak well-posedness for $$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$7 with constant coefficients under polynomial-growth assumptions on $$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$8 and showed uniqueness under stronger regularity assumptions on the potential (Bosia et al., 2014). In that framework, weak solutions satisfy

$$-\nabla\cdot(\nu D(\mathbf u))+\eta \mathbf u+\nabla p=\mu\nabla\phi,\qquad \nabla\cdot\mathbf u=0,$$9

with mass conservation and the energy identity as central a priori tools (Bosia et al., 2014). The same work proved eventual higher regularity, existence of a global attractor, and convergence of trajectories to single equilibria via the Łojasiewicz–Simon inequality (Bosia et al., 2014).

For the nonlocal CHB system,

$\phi$0

global weak existence follows under assumptions on $\phi$1, $\phi$2, viscosity, permeability, and forcing (Porta et al., 2016). With constant viscosity, weak solutions are unique, and bounded initial data yield $\phi$3 (Porta et al., 2016). In the two-dimensional regular-potential setting, this theory can be upgraded to strong solutions with

$\phi$4

(Dharmatti et al., 2019).

Singular-potential theory is substantially more delicate. For tumor-growth CHB systems with source terms and non-solenoidal velocity, weak and stationary solutions have been established for both double-obstacle and logarithmic potentials (Ebenbeck et al., 2019). A key difficulty there is that source terms destroy the usual conservation of the mean phase, so one must replace the standard mass-conservation argument by a proof that the spatial mean of $\phi$5 remains strictly inside $\phi$6; this then yields control of the mean of $\phi$7, which is essential for singular-potential analysis (Ebenbeck et al., 2019). A related tumor-growth model with Dirichlet boundary condition for the chemical potential was shown to admit weak solutions for regular, logarithmic, and double-obstacle potentials; the Dirichlet condition is analytically decisive because it allows direct control of $\phi$8 in $\phi$9 from $\mu$0 without requiring a separate estimate of its mean (Colli et al., 2022).

The sixth-order Brinkman–Cahn–Hilliard system with curvature effects admits weak solutions in three dimensions under variable shear viscosity, variable drag, and nonconservative source terms (Colli et al., 24 Sep 2025). In the divergence-free variational framework,

$\mu$1

$\mu$2

and the main existence theorem provides a uniform a priori bound of this full mixed system (Colli et al., 24 Sep 2025). Uniqueness and continuous dependence are proved only when mobility and shear viscosity are constants, precisely because difference estimates become tractable only in the absence of nonlinear coefficient variations (Colli et al., 24 Sep 2025).

5. Singular potentials, degenerate mobility, and separation phenomena

Singular free energies and degenerate mobilities play a special role because they enforce physically meaningful bounds on the phase and modify regularity mechanisms. In the nonlocal CHB system with logarithmic potential

$\mu$3

and degenerate mobility

$\mu$4

the main recent result is a separation property in two dimensions: if $\mu$5 and $\mu$6, then for any $\mu$7 there exists $\mu$8 such that

$\mu$9

(Dharmatti et al., 6 Feb 2025). Equivalently, after any positive waiting time the solution stays uniformly away from the pure phases $\mathbf u$0.

This result is nontrivial because the singularity of $\mathbf u$1 near $\mathbf u$2 and the vanishing of $\mathbf u$3 at $\mathbf u$4 create opposite effects: singularity repels the solution from pure phases, while degeneracy weakens diffusion there (Dharmatti et al., 6 Feb 2025). The proof uses a De Giorgi iteration adapted to the nonlocal Brinkman-coupled setting, together with a crucial uniform $\mathbf u$5 bound on $\mathbf u$6 and an entropy test built from a function $\mathbf u$7 satisfying $\mathbf u$8 (Dharmatti et al., 6 Feb 2025). The transport term caused by the Brinkman velocity is harmless in the iteration because incompressibility gives exact cancellation when truncation functions are used (Dharmatti et al., 6 Feb 2025).

Singular potentials also dominate tumor-growth CHB theory. The double-obstacle and logarithmic cases ensure $\mathbf u$9 remains in the physically relevant interval $p$0 or $p$1, which is important when $p$2 is interpreted as a difference of volume fractions (Ebenbeck et al., 2019, Colli et al., 2022). In modified chemotaxis-Brinkman-Cahn–Hilliard models, the Flory–Huggins logarithmic potential

$p$3

is used together with nonlinear chemotactic sensitivity; the singular potential yields a uniform bound on $p$4, which is essential for controlling the cross-diffusion structure (Schimperna, 2024).

6. Limits, control, computation, and extensions

A major structural theme is the vanishing-viscosity limit from Brinkman to Darcy or Hele–Shaw dynamics. In the classical local CHB system, letting the Brinkman viscosity parameter tend to zero yields the Cahn–Hilliard–Hele–Shaw system, and in two dimensions there is even a quantitative estimate comparing CHB and CHHS solutions (Bosia et al., 2014). For the nonlocal CHB system, weak solutions of CHHS are obtained as limits of CHB solutions under $p$5, and bounded weak solutions of the nonlocal CHHS problem are unique under constant permeability (Porta et al., 2016). In tumor-growth CHB with nutrient coupling, vanishing viscosities lead to a Cahn–Hilliard–Darcy system with pressure satisfying an elliptic problem and velocity determined algebraically from capillary and nutrient forcing (Ebenbeck et al., 2018). In the sixth-order curvature model, the Darcy limit is rigorously established by a sequence of viscosity coefficients $p$6 with $p$7, producing a reduced variational law

$p$8

(Colli et al., 24 Sep 2025).

Optimal control theory has recently become an active part of the Brinkman–Cahn–Hilliard literature. For the sixth-order system with constant mobility and viscosity, distributed forcing in the Brinkman equation has been treated as the control variable, and first-order necessary optimality conditions have been derived both for smooth quadratic costs and for sparsity-promoting functionals containing an $p$9 term (Bag, 8 Dec 2025). In that setting the control-to-state map is Fréchet differentiable, the adjoint system can be written explicitly, and the optimality condition takes the form of a variational inequality; for $\mu\nabla\phi$0, the nonsmooth term yields sparse controls through subdifferential conditions (Bag, 8 Dec 2025). A closely related higher-order optimal velocity control problem for a curvature-regularized Brinkman–Cahn–Hilliard system also derives adjoint-based first-order conditions and pointwise projection formulas for box-constrained controls (Colli et al., 9 Jan 2026). In two-dimensional nonlocal CHB with regular potential, distributed forcing control has likewise been analyzed using strong-solution theory, a linearized system, and an adjoint system (Dharmatti et al., 2019).

Numerically, recent work has focused on high-order, unconditionally energy-stable discretizations. A class of variable-step R-IMEX-BDF$\mu\nabla\phi$1 schemes based on a generalized scalar auxiliary variable approach with relaxation was developed for the CHB system with periodic boundary conditions (Chen et al., 8 Jun 2025). These schemes are proved unconditionally energy-stable in the modified scalar energy variable, remain stable without adjacent-step-ratio restrictions, and support adaptive and hybrid-order time-stepping strategies tuned to the multiscale dynamics of phase separation (Chen et al., 8 Jun 2025). The numerical experiments include manufactured-solution convergence tests, spinodal decomposition, and buoyancy-driven flow, and show that the hybrid adaptive algorithm effectively combines short-time third-order accuracy with long-time second-order robustness (Chen et al., 8 Jun 2025).

The model class has also expanded beyond deterministic bulk equations. A stochastic Cahn–Hilliard–Brinkman model with dynamic boundary conditions and independent multiplicative Wiener noises in bulk and surface phase equations has been formulated and weak martingale solutions have been proved in dimensions two and three (Brzeźniak et al., 10 Jan 2026). This stochastic extension preserves the deterministic Brinkman–Cahn–Hilliard coupling but destroys exact mass conservation, which introduces substantial new difficulties in controlling chemical-potential means and boundary phase variables (Brzeźniak et al., 10 Jan 2026).

Overall, the Brinkman–Cahn–Hilliard system is best understood not as a single PDE model but as a family of coupled diffuse-interface flow systems organized around a common mechanism: Cahn–Hilliard thermodynamics for phase segregation coupled to Brinkman momentum transport for porous or porous-like flow. The main directions of current research are the treatment of singular and higher-order energies, rigorous Brinkman-to-Darcy limits, nonlocal and chemotactic couplings, optimal control, robust energy-stable numerics, and stochastic generalizations (Bosia et al., 2014, Porta et al., 2016, Ebenbeck et al., 2018, Colli et al., 24 Sep 2025, Bag, 8 Dec 2025, Chen et al., 8 Jun 2025, Brzeźniak et al., 10 Jan 2026).