Abstract:<p>The Cahn-Hilliard-Navier-Stokes (CHNS) equations represent the fundamental building blocks of hydrodynamic phase-field models for multiphase fluid flow dynamics. Due to the coupling between the Navier-Stokes equation and the Cahn-Hilliard equation, the CHNS system is non-trivial to solve numerically. Traditionally, a numerical extrapolation for the coupling terms is used. However, such brute-force extrapolation usually destroys the intrinsic thermodynamic structures of this CHNS system. This paper proposes a new strategy to reformulate the CHNS system into a constraint gradient flow formation. Under the new formulation, the reversible and irreversible structures are clearly revealed. This guides us to propose operator splitting schemes that have several advantageous properties. First of all, the proposed schemes lead to several decoupled systems in smaller sizes to be solved at each time marching step. This significantly reduces computational costs. Secondly, the proposed schemes still guarantee the thermodynamic laws of the CHNS system at the discrete level. In addition, unlike the recently populated IEQ or SAV approach using auxiliary variables, our resulting energy laws are formulated in the original variables. This is a significant improvement, as the modified energy laws with auxiliary variables sometimes deviate from the original energy law. Our proposed framework lays a foundation for designing decoupled and energy stable numerical algorithms for hydrodynamic phase-field models. Furthermore, various numerical algorithms can be obtained given different splitting steps, making this framework rather general. The proposed numerical algorithms are implemented. Their second-order temporal and spatial accuracy are verified numerically. Some numerical examples and benchmark problems are calculated to verify the effectiveness of the proposed schemes.</p>

Fully decoupled unconditionally stable Crank–Nicolson leapfrog numerical methods for the Cahn–Hilliard–Darcy system

A new decoupled unconditionally stable scheme and its optimal error analysis for the Cahn-Hilliard-Navier-Stokes equations

Unconditionally Energy Stable Second-Order Numerical Schemes for the Functionalized Cahn-Hilliard Gradient Flow Equation Based on the SAV Approach.

Second‐order, Fully Decoupled, Linearized, and Unconditionally Stable Scalar Auxiliary Variable Schemes for Cahn–Hilliard–Darcy System

A Second-Order Decoupled Energy Stable Numerical Scheme for Cahn-Hilliard-Hele-Shaw System

Linear and Energy-Stable Method with Enhanced Consistency for the Incompressible Cahn–Hilliard–Navier–Stokes Two-Phase Flow Model

A Family of Second-Order Energy-Stable Schemes for Cahn–Hilliard Type Equations

Decoupled, Linear, and Energy Stable Finite Element Method for the Cahn-Hilliard-Navier-Stokes-Darcy Phase Field Model.

Highly efficient and unconditionally energy stable semi-discrete time-marching numerical scheme for the two-phase incompressible flow phase-field system with variable-density and viscosity

On Efficient Second Order Stabilized Semi-implicit Schemes for the Cahn–Hilliard Phase-Field Equation

A high-order accurate unconditionally stable bound-preserving numerical scheme for the Cahn-Hilliard-Navier-Stokes equations

Second-order decoupled energy-stable schemes for Cahn-Hilliard-Navier-Stokes equations

Crank–Nicolson Leap-Frog Time Stepping Decoupled Scheme for the Fluid–Fluid Interaction Problems

Convergence analysis and error estimates for a second order accurate finite element method for the Cahn–Hilliard–Navier–Stokes system

A fully-decoupled second-order-in-time and unconditionally energy stable scheme for a phase-field model of two phase flow with variable density

A Conservative Numerical Method for the Cahn-Hilliard Equation with Dirichlet Boundary Conditions in Complex Domains

Energy Stable Numerical Schemes for the Fractional-in-space Cahn–Hilliard Equation

An Unconditionally Energy-Stable Second-Order Time-Accurate Numerical Scheme For The Coupled Cahn-Hilliard System In Copolymer/Homopolymer Mixtures

Double Stabilizations and Convergence Analysis of a Second-Order Linear Numerical Scheme for the Nonlocal Cahn-Hilliard Equation

A novel fully decoupled scheme with second‐order time accuracy and unconditional energy stability for the Navier‐Stokes equations coupled with mass‐conserved Allen‐Cahn phase‐field model of two‐phase incompressible flow

Convergence and Stability Analysis of Energy Stable and Bound‐preserving Numerical Schemes for Binary Fluid‐surfactant Phase‐field Equations