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 Positivity-Preserving, Energy Stable BDF2 Scheme with Variable Steps for the Cahn–Hilliard Equation with Logarithmic Potential

A Uniquely Solvable, Positivity-Preserving and Unconditionally Energy Stable Numerical Scheme for the Functionalized Cahn-Hilliard Equation with Logarithmic Potential

Efficient Variable Steps BDF2 Scheme for the Two-Dimensional Space Fractional Cahn-Hilliard Model

On energy stable, maximum-principle preserving, second order BDF scheme with variable steps for the Allen-Cahn equation

Convergence analysis of variable steps BDF2 method for the space fractional Cahn-Hilliard model

Asymptotically compatible energy of variable-step fractional BDF2 formula for time-fractional Cahn-Hilliard model

Analysis of the second order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection

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

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

Structure-preserving weighted BDF2 methods for Anisotropic Cahn-Hilliard model: uniform/variable-time-steps

Unconditionally optimal error estimate of a linearized variable-time-step BDF2 scheme for nonlinear parabolic equations

A Second-order Energy Stable Fourier Spectral Method for the Fractional Cahn-Hilliard Equation

Variable-step numerical schemes and energy dissipation laws for time fractional Cahn–Hilliard model

Energy Stable Numerical Schemes for the Fractional-in-space Cahn–Hilliard Equation

Convergence analysis of the maximum principle preserving BDF2 scheme with variable time-steps for the space fractional Allen-Cahn equation

Analysis of adaptive BDF2 scheme for diffusion equations

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

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

A Bound Preserving Energy Stable Scheme for a Nonlocal Cahn-Hilliard Equation

Stabilized energy factorization approach for Allen-Cahn equation with logarithmic Flory-Huggins potential

Convergence Analysis of a Preconditioned Steepest Descent Solver for the Cahn-Hilliard Equation with Logarithmic Potential