A Modified Crank-Nicolson Numerical Scheme for the Flory-Huggins Cahn-HilliardModel
In this paper we propose and analyze a second order accurate numerical scheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. a modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order Adams- Bashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. a nonlinear artificial regularization term is added in the numerical scheme, which ensures the positivity-preserving property, i.e., the numerical value of the phase variable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special form of the logarithmic approximation term. in addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearized stability analysis. a few numerical results, including both the constant-mobility and solution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme.
W. Chen et al., "A Modified Crank-Nicolson Numerical Scheme for the Flory-Huggins Cahn-HilliardModel," Communications in Computational Physics, vol. 31, no. 1, pp. 60 - 93, Global Science Press, Jan 2021.
The definitive version is available at https://doi.org/10.4208/cicp.OA-2021-0074
Mathematics and Statistics
Keywords and Phrases
Cahn-Hilliard Equation; Energy Stability; Flory Huggins Energy Potential; Optimal Rate Convergence Estimate; Positivity Preserving; Second Order Accuracy
International Standard Serial Number (ISSN)
Article - Journal
© 2023 Global Science Press, All rights reserved.
01 Jan 2021