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>

A Parallel Domain Decomposition-Based Implicit Method for the Cahn-Hilliard-Cook Phase-Field Equation in 3D

The high-order exponential semi-implicit scalar auxiliary variable approach for nonlocal Cahn-Hilliard equation

Parallel Energy-Stable Phase Field Crystal Simulations Based on Domain Decomposition Methods.

Efficient Numerical Simulation of Cahn-Hilliard Type Models by a Dimension Splitting Method

An Implicit Parallel Multigrid Computing Scheme to Solve Coupled Thermal-Solute Phase-Field Equations for Dendrite Evolution

An Efficient Finite Difference-based Implicit Solver for Phase-Field Equations with Spatially and Temporally Varying Parameters

A Parallel 3D Dendritic Growth Simulator Using the Phase-Field Method

A Scalable Implicit Solver for Phase Field Crystal Simulations

Parallel and Fully Implicit Simulations of the Thermal Convection in the Earth’s Outer Core

NKS Method for the Implicit Solution of a Coupled Allen-Cahn/Cahn-Hilliard System

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

Parallel Domain Decomposition Methods with Mixed Order Discretization for Fully Implicit Solution of Tracer Transport Problems on the Cubed-Sphere.

A parallel finite volume procedure for phase-field simulation of solidification

Simulation of Cavity Flows by an Implicit Domain Decomposition Algorithm for the Lattice Boltzmann Equations

Parallel Numerical Simulation of the Thermal Convection in the Earth'S Outer Core on the Cubed-Sphere

An efficient phase-field method for turbulent multiphase flows

Development of a parallel adaptive multigrid algorithm for solving the multi-scale thermal-solute 3D phase-field problems

On the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

Projection Method for Steady states of Cahn-Hilliard Equation with the Dynamic Boundary Condition

Parallel Energy-Stable Solver for a Coupled Allen--Cahn and Cahn--Hilliard System

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