Stabilized immersed isogeometric analysis for the Navier-Stokes-Cahn-Hilliard equations, with applications to binary-fluid flow through porous media

Stein K.F. Stoter,Tom B. van Sluijs,Tristan H.B. Demont,E. Harald van Brummelen,Clemens V. Verhoosel

2023-05-31

Abstract:Binary-fluid flows can be modeled using the Navier-Stokes-Cahn-Hilliard equations, which represent the boundary between the fluid constituents by a diffuse interface. The diffuse-interface model allows for complex geometries and topological changes of the binary-fluid interface. In this work, we propose an immersed isogeometric analysis framework to solve the Navier-Stokes-Cahn-Hilliard equations on domains with geometrically complex external binary-fluid boundaries. The use of optimal-regularity B-splines results in a computationally efficient higher-order method. The key features of the proposed framework are a generalized Navier-slip boundary condition for the tangential velocity components, Nitsche's method for the convective impermeability boundary condition, and skeleton- and ghost-penalties to guarantee stability. A binary-fluid Taylor-Couette flow is considered for benchmarking. Porous medium simulations demonstrate the ability of the immersed isogeometric analysis framework to model complex binary-fluid flow phenomena such as break-up and coalescence in complex geometries.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to simulate and analyze the flow phenomena of binary fluids in porous media, especially in cases involving complex topological changes such as interface rupture and fusion. Specifically, the authors propose a new immersed isogeometric analysis framework for solving the Navier - Stokes - Cahn - Hilliard equations that describe the flow of binary fluids.

### Main problems
1. **Complex physical behaviors**: The flow of binary fluids involves complex physical phenomena, such as the rupture and fusion of the interface layer. These phenomena not only change the fluid shape but also the topological structure.
2. **Contact line dynamics**: The dynamic behavior of the contact line (i.e., the intersection of the interface layer and the solid boundary) is crucial to fluid flow and must be accurately simulated by the computational model.
3. **Complex geometries**: The geometric shape of the binary fluid flow region can be very complex (for example, scan - based microstructures), which makes automatic high - quality mesh generation almost impossible.

### Solutions
To address the above challenges, the authors propose the following solutions:
- **Diffusion interface model**: Adopt a diffusion interface model to describe the boundary between binary fluids, allowing for complex geometries and topological changes.
- **Immersed isogeometric analysis**: Combine the optimal - regularity B - spline to propose an efficient high - order method to solve the Navier - Stokes - Cahn - Hilliard equations for the complex external binary fluid boundary.
- **Generalized Navier slip boundary condition**: For the tangential velocity component, use the generalized Navier slip boundary condition to avoid traction singularity.
- **Nitsche method**: For the convective impermeable boundary condition, use the Nitsche method.
- **Skeleton and ghost penalty**: Introduce skeleton and ghost penalty terms to ensure stability.

### Application examples
The paper demonstrates the effectiveness of this method in handling complex fluid phenomena (such as rupture and fusion) through Taylor - Couette flow and binary fluid flow simulations in porous media.

### Key formulas
The main formulas involved in the paper include:

- **Navier - Stokes - Cahn - Hilliard equations**:
\[
\begin{aligned}
&\partial_t(\rho \mathbf{u})+\nabla\cdot(\rho \mathbf{u}\otimes\mathbf{u})+\nabla\cdot(\mathbf{u}\otimes\mathbf{J})-\nabla\cdot\boldsymbol{\tau}-\nabla\cdot\boldsymbol{\zeta}+\nabla p = 0, \\
&\nabla\cdot\mathbf{u} = 0, \\
&\partial_t\varphi+\nabla\cdot(\varphi \mathbf{u})-\nabla\cdot(m\nabla\mu)= 0.
\end{aligned}
\]

- **Chemical potential**:
\[
\mu=-\sigma\epsilon\Delta\varphi+\frac{\sigma}{\epsilon}\Psi'(\varphi).
\]

- **Generalized Navier boundary condition**:
\[
P_{\Gamma}(\alpha_{GN}(\mathbf{u}-\mathbf{u}_{\text{wall}})+(\boldsymbol{\tau}\mathbf{n}+\boldsymbol{\zeta}\mathbf{n}))-\nabla_{\Gamma}\sigma_{sf}(\varphi)= 0.
\]

Through these methods and techniques, the authors have successfully developed a method that can be efficient and accurate.