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>

Unconditional energy stable IEQ-FEMs for the Cahn-Hilliard-Navier-Stokes equations

Unconditionally energy stable IEQ-FEMs for the Cahn-Hilliard equation and Allen-Cahn equation

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

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

The subdivision-based IGA-EIEQ numerical scheme for the Navier–Stokes equations coupled with Cahn–Hilliard phase-field model of two-phase incompressible flow on complex curved surfaces

An Efficient Unconditional Energy-Stable Finite Element Method for the Electro-Hydrodynamic Equations

Semi-implicit, Unconditionally Energy Stable, Stabilized Finite Element Method Based on Multiscale Enrichment for the Cahn-Hilliard-Navier-Stokes Phase-Field Model

Efficient Numerical Schemes with Unconditional Energy Stabilities for the Modified Phase Field Crystal Equation

Unconditionally stable exponential integrator schemes for the 2D Cahn-Hilliard equation

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

Two novel numerical methods for gradient flows: generalizations of the Invariant Energy Quadratization method

Highly efficient, robust and unconditionally energy stable second order schemes for approximating the Cahn-Hilliard-Brinkman system

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

An Unconditionally Energy Stable Method for Binary Incompressible Heat Conductive Fluids Based on the Phase-Field Model

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

Gauge-Uzawa-based, highly efficient decoupled schemes for the diffuse interface model of two-phase magnetohydrodynamic

Global error estimates of high-order fully decoupled schemes for the Cahn-Hilliard-Navier-Stokes model of Two-Phase Incompressible Flows

An exponential-free Runge--Kutta framework for developing third-order unconditionally energy stable schemes for the Cahn--Hilliard equation

Error estimates for the highly efficient and energy stable schemes for the 2D/3D two-phase MHD

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