Abstract:In this paper we propose and analyze a finite difference numerical scheme for the Flory-Huggins-Cahn-Hilliard equation with dynamical boundary condition. The singular logarithmic potential is included in the Flory-Huggins energy expansion. Meanwhile, a dynamical evolution equation for the boundary profile corresponds to a lower-dimensional singular energy potential. In turn, a theoretical analysis for the coupled system becomes very challenging, since it contains nonlinear and singular energy potentials for both the interior region and on the boundary. In the numerical design, a convex splitting approach is applied to the chemical potential associated with the energy both at the interior region and on the boundary: implicit treatments for the singular and logarithmic terms, as well as the surface diffusion terms, combined with an explicit treatment for the concave expansive term. In addition, the discrete boundary condition for the phase variable is coupled with the evolutionary equation of the boundary profile. The resulting numerical system turns out to be highly nonlinear, singular and coupled. A careful finite difference approximation and convexity analysis reveals that such a numerical system could be represented as a minimization of a discrete numerical energy functional, which contains both the interior and boundary integrals. More importantly, all the singular terms correspond to a discrete convex functional. As a result, a unique solvability and positivity-preserving analysis could be theoretically justified, based on the subtle fact that the singular nature of the logarithmic terms around the singular limit values prevent the numerical solutions reaching these values. The total energy stability analysis could be established by a careful estimate over the finite difference inner product. Some numerical results are presented in this article.

A fully discrete evolving surface finite element method for the Cahn-Hilliard equation with a regular potential

An evolving surface finite element method for the Cahn-Hilliard equation with a logarithmic potential

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

Navier-Stokes-Cahn-Hilliard equations on evolving surfaces

An Unconditionally Energy-Stable Second-Order Time-Accurate Scheme for the Cahn–Hilliard Equation on Surfaces

An integral equation method for the Cahn-Hilliard equation in the wetting problem

Error estimates for full discretization of Cahn--Hilliard equation with dynamic boundary conditions

A Nonconforming Finite Element Method for the Cahn-Hilliard Equation

The evolving surface Cahn-Hilliard equation with a degenerate mobility

A structure-preserving discontinuous Galerkin scheme for the Cahn-Hilliard equation including time adaptivity

Error Estimates of Fully Discrete Finite Element Solutions for the 2D Cahn–Hilliard Equation with Infinite Time Horizon

Continuous Linear Finite Element Method for Biharmonic Problems on Surfaces

A higher order nonconforming virtual element method for the Cahn-Hilliard equation

An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation

Local Discontinuous Galerkin Methods for Solving Convection-Diffusion and Cahn-Hilliard Equations on Surfaces

Efficient and stable time integration of Cahn-Hilliard equations: explicit, implicit and explicit iterative schemes

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

A New Mixed Finite Element Method For The Cahn-Hilliard Equation

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

A uniquely solvable and positivity-preserving finite difference scheme for the Flory-Huggins-Cahn-Hilliard equation with dynamical boundary condition

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