Abstract:We propose and analyze a first-order finite difference scheme for the functionalized Cahn-Hilliard (FCH) equation with a logarithmic Flory-Huggins potential. The semi-implicit numerical scheme is designed based on a suitable convex-concave decomposition of the FCH free energy. We prove unique solvability of the numerical algorithm and verify its unconditional energy stability without any restriction on the time step size. Thanks to the singular nature of the logarithmic part in the Flory-Huggins potential near the pure states $\pm 1$, we establish the so-called positivity-preserving property for the phase function at a theoretic level. As a consequence, the numerical solutions will never reach the singular values $\pm 1$ in the point-wise sense and the fully discrete scheme is well defined at each time step. Next, we present a detailed optimal rate convergence analysis and derive error estimates in $l^{\infty}(0,T;L_h^2)\cap l^2(0,T;H^3_h)$ under a linear refinement requirement $\Delta t\leq C_1 h$. To achieve the goal, a higher order asymptotic expansion (up to the second order temporal and spatial accuracy) based on the Fourier projection is utilized to control the discrete maximum norm of solutions to the numerical scheme. We show that if the exact solution to the continuous problem is strictly separated from the pure states $\pm 1$, then the numerical solutions can be kept away from $\pm 1$ by a positive distance that is uniform with respect to the size of the time step and the grid. Finally, a few numerical experiments are presented. Convergence test is performed to demonstrate the accuracy and robustness of the proposed numerical scheme. Pearling bifurcation, meandering instability and spinodal decomposition are observed in the numerical simulations.

A two‐grid finite element method for the Allen‐Cahn equation with the logarithmic potential

A two-level finite element method for the Allen–Cahn equation

An efficient two-grid method for the Cahn–Hilliard equation with the concentration-dependent mobility and the logarithmic Flory-Huggins bulk potential

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

A Two-Grid Finite Element Method for a Second-Order Nonlinear Hyperbolic Equation

Lp error estimate of nonlinear Schrödinger equation using a two‐grid finite element method

Maximum-Principle-Preserving Local Discontinuous Galerkin Methods for Allen-Cahn Equations

Two-grid method for two-dimensional nonlinear Schrödinger equation by mixed finite element method

Superconvergence analysis of two-grid finite element method for nonlinear Benjamin-Bona-Mahony equation

Adaptive discontinuous Galerkin finite element methods for the Allen-Cahn equation on polygonal meshes

Hybrid numerical method for the Allen–Cahn equation on nonuniform grids

Two-grid Methods for Finite Volume Element Approximations of Nonlinear Sobolev Equations

Optimal H1 spatial convergence of a fully discrete finite element method for the time-fractional Allen-Cahn equation

Two‐grid finite element method on grade meshes for time‐fractional nonlinear Schrödinger equation

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

A two-grid finite element method for nonlinear parabolic integro-differential equations

Adaptive linear second-order energy stable schemes for time-fractional Allen-Cahn equation with volume constraint

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

A two-grid method based on Newton iteration for the Navier-Stokes equations

Analysis of two-grid method for second-order hyperbolic equation by expanded mixed finite element methods

Parallel Energy-Stable Solver for a Coupled Allen--Cahn and Cahn--Hilliard System