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>

Development of a numerical phase optimized upwinding combined compact difference scheme for solving the Camassa–Holm equation with different initial solitary waves

Development of Two-Dimensional Non-Hydrostatic Wave Model Based on Central-Upwind Scheme

Numerical Modeling of Lock-Exchange Gravity/turbidity Currents by a High-Order Upwinding Combined Compact Difference Scheme

Long-time asymptotic solution structure of Camassa-Holm equation subject to an initial condition with non-zero reflection coefficient of the scattering data

Up-wind Scheme Based on the Undetermined Coefficient Method

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

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

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation

Error estimates of invariant-preserving difference schemes for the rotation-two-component Camassa--Holm system with small energy

Conservative second-order finite difference method for Camassa–Holm equation with periodic boundary condition

Fully Conservative Difference Schemes for the Rotation-Two-Component Camassa-Holm System with Smooth/Nonsmooth Initial Data

Well-balanced fifth-order finite difference Hermite WENO scheme for the shallow water equations

Error estimates of invariant-preserving difference schemes for the rotation-two-component Camassa–Holm system with small energy

Dynamics of the modified Rotation-Camassa-Holm Equation with backward diffusion term

On the long-time asymptotics of the modified Camassa-Holm equation with step-like initial data

A Riemann-Hilbert approach to the two-component modified Camassa-Holm equation

Fourth order finite difference methods for the wave equation with mesh refinement interfaces

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

Error estimates for Galerkin finite element methods for the Camassa–Holm equation

A Less Time-Consuming Upwind Compact Difference Method with Adjusted Dissipation Property for Solving the Unsteady Incompressible Navier-Stokes Equations

Multiquadric Solution For Shallow Water Equations