Abstract:A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment for the singular logarithmic terms and the surface diffusion term, and an explicit update for the expansive concave term. The convective term for the phase variable, as well as the coupled term in the Stokes equation, are approximated in a semi-implicit manner. In the spatial discretization, the marker and cell (MAC) difference method is applied, which evaluates the velocity components, the pressure and the phase variable at different cell locations. Such an approach ensures the divergence-free feature of the discrete velocity, and this property plays an important role in the analysis. The positivity-preserving property and the unique solvability of the proposed numerical scheme are theoretically justified, utilizing the singular nature of the logarithmic term as the phase variable approaches the singular limit values. An unconditional energy stability analysis is standard, as an outcome of the convex-concave decomposition technique. A convergence analysis with accompanying error estimate is provided for the proposed numerical scheme. In particular, a higher order consistency analysis, accomplished by supplementary functions, is performed to ensure the separation properties of numerical solution. In turn, using the approach of rough and refined error (RRE) estimates, we are able to derive an optimal rate convergence. To conclude, several numerical experiments are presented to validate the theoretical analysis.

A Conservative Numerical Method for the Cahn-Hilliard Equation with Generalized Mobilities on Curved Surfaces in Three-Dimensional Space

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

A Conservative Finite Difference Scheme for the N-component Cahn–Hilliard System on Curved Surfaces in 3D

Explicit Hybrid Numerical Method for the Allen-Cahn Type Equations on Curved Surfaces

Direct Discretization Method for the Cahn---Hilliard Equation on an Evolving Surface

A Compact Fourth-Order Finite Difference Scheme for the Three-Dimensional Cahn–Hilliard Equation

An Efficient Numerical Method for Reaction-Diffusion Equation on the General Curved Surfaces.

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

A Finite Volume / Discontinuous Galerkin Method for the Advective Cahn–Hilliard Equation with Degenerate Mobility on Porous Domains Stemming from Micro-Ct Imaging

Two linear energy stable lumped mass finite element schemes for the viscous Cahn-Hilliard equation on curved surfaces in 3D

A Practical Finite Difference Scheme for the Navier-Stokes Equation on Curved Surfaces in R3

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

A fourth-order spatial accurate and practically stable compact scheme for the Cahn–Hilliard equation

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

First- and Second-Order Unconditionally Stable Direct Discretization Methods for Multi-Component Cahn–Hilliard System on Surfaces

A Simple Benchmark Problem for the Numerical Methods of the Cahn–Hilliard Equation

A Curved Boundary Treatment for Discontinuous Galerkin Method Applied to Euler Equations on Triangular and Tetrahedral Grids.

A Stable Second-Order BDF Scheme for the Three-Dimensional Cahn–Hilliard–Hele–Shaw System

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

Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential

Structure preserving finite element schemes for the Navier-Stokes-Cahn-Hilliard system with degenerate mobility