Abstract:A fully decoupled, linearized, and unconditionally stable finite element method is developed to solve the Cahn–Hilliard–Navier–Stokes–Darcy model in the coupled free fluid region and porous medium region. By introducing two auxiliary energy variables, we derive the equivalent system that is consistent with the original system. The energy dissipation law of the proposed equivalent model is proven. To lay a solid foundation, we first present a coupled linearized time-stepping method for the reformulated system, and prove its unconditionally energy stability. In order to further improve the computational efficiency, special treatment for the interface conditions and the artificial compression approach are utilized to decouple the two subdomains and the Navier–Stokes equation. Therefore, with the discretization techniques of two existing auxiliary variable approaches, a fully decoupled and linearized numerical scheme can be developed, under the framework of a semi-implicit semi-explicit scheme for temporal discretization and Galerkin finite element method for spatial discretization. The grad-div stabilization is also employed to further improve the stability of auxiliary variable algorithm. The full discretization obeys the desired energy dissipation law without any temporal restriction. Moreover, the implementation process is discussed, including the adaptive mesh strategy to accurately capture the diffuse interface. Ample numerical experiments are performed to validate the typical features of developed numerical schemes, such as the accuracy, energy stability without restriction for time step size, and adaptive mesh refinement in space. Furthermore, we apply the proposed numerical method to simulate the shape relaxation and the Buoyancy-driven flows, which demonstrate the applicability of the proposed method.

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

STABILIZED FINITE ELEMENT METHODS BASED ON MULTISCALE ENRICHMENT FOR ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS

A Linear Unconditionally Stable Scheme for the Incompressible Cahn–Hilliard–Navier–Stokes Phase-Field Model

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

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

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

Discontinuous Finite Volume Element Method for a Coupled Navier-Stokes-Cahn-Hilliard 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 the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

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

Decoupled, Energy Stable Scheme for Hydrodynamic Allen-Cahn Phase Field Moving Contact Line Model

Stabilized Crank-Nicolson/Adams-Bashforth Schemes for Phase Field Models

Analysis of the Energy Stability for Stabilized Semi-implicit Schemes of the Functionalized Cahn-Hilliard Mass-conserving Gradient Flow Equation.

Fully Discrete Energy Stable Scheme for a Phase-Field Moving Contact Line Model with Variable Densities and Viscosities

Decoupled, linear, unconditionally energy stable and charge-conservative finite element method for an inductionless magnetohydrodynamic phase-field model

An Energy-Stable Finite Element Method for Incompressible Magnetohydrodynamic-Cahn-Hilliard Coupled Model

Error Analysis of a Decoupled, Linear Stabilization Scheme for the Cahn–Hilliard Model of Two-Phase Incompressible Flows

Efficient Fully Discrete Finite-Element Numerical Scheme with Second-Order Temporal Accuracy for the Phase-Field Crystal Model

Efficient and Energy Stable Scheme for the Hydrodynamically Coupled Three Components Cahn-Hilliard Phase-Field Model Using the Stabilized-Invariant Energy Quadratization (S-IEQ) Approach.

A Fully Decoupled Numerical Method for Cahn–Hilliard–Navier–Stokes–Darcy Equations Based on Auxiliary Variable Approaches

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