Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential

被引：18

作者：

Chen W. ^{[1
,4
]}

Wang C. ^{[2
]}

Wang X. ^{[3
,4
,5
]}

Wise S.M. ^{[6
]}

机构：

[1] Shanghai Key Laboratory for Contemporary Applied Mathematics, School of Mathematical Sciences, Fudan University, Shanghai

[2] Mathematics Department, University of Massachusetts, North Dartmouth, 02747, MA

[3] Department of Mathematics, Southern University of Science and Technology, Shenzhen

[4] Fudan University, Shanghai

[5] Florida State University, Tallahassee, 32306, FL

[6] Mathematics Department, University of Tennessee, Knoxville, 37996, TN

来源：

Journal of Computational Physics: X | 2019年 / 3卷

基金：

美国国家科学基金会; 中国国家自然科学基金;

关键词：

Cahn-Hilliard equation; Energy stability; Logarithmic Flory Huggins energy potential; Optimal rate convergence analysis; Positivity preserving; Second order BDF scheme;

D O I：

10.1016/j.jcpx.2019.100031

中图分类号：

学科分类号：

摘要：

In this paper we present and analyze finite difference numerical schemes for the Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential. Both first and second order accurate temporal algorithms are considered. In the first order scheme, we treat the nonlinear logarithmic terms and the surface diffusion term implicitly, and update the linear expansive term and the mobility explicitly. We provide a theoretical justification that this numerical algorithm has a unique solution, such that the positivity is always preserved for the logarithmic arguments, i.e., the phase variable is always between −1 and 1, at a point-wise level. In particular, our analysis reveals a subtle fact: the singular nature of the logarithmic term around the values of −1 and 1 prevents the numerical solution reaching these singular values, so that the numerical scheme is always well-defined as long as the numerical solution stays similarly bounded at the previous time step. Furthermore, an unconditional energy stability of the numerical scheme is derived, without any restriction for the time step size. Such an analysis technique can also be applied to a second order numerical scheme in which the BDF temporal stencil is applied, the expansive term is updated by a second order Adams-Bashforth explicit extrapolation formula, and an artificial Douglas-Dupont regularization term is added to ensure the energy dissipativity. The unique solvability and the positivity-preserving property for the second order scheme are proved using similar ideas, namely, the singular nature of the logarithmic term plays an essential role. For both the first and second order accurate schemes, we are able to derive an optimal rate convergence analysis. The case with a non-constant mobility is analyzed as well. We also describe a practical and efficient multigrid solver for the proposed numerical schemes, and present some numerical results, which demonstrate the robustness of the numerical schemes. © 2019

引用

共 67 条

[1]

Abels H., On a diffuse interface model for two-phase flows of viscous, incompressible fluids with matched densities, Arch. Ration. Mech. Anal., 194, 2, pp. 463-506, (2009)

[2]

Abels H., Wilke M., Convergence to equilibrium for the Cahn-Hilliard equation with a logarithmic free energy, Nonlinear Anal., 67, pp. 3176-3193, (2007)

[3]

Allen S.M., Cahn J.W., A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metall., 27, (1979)

[4]

Barrett J., Blowey J., An error bound for the finite element approximation of the Cahn-Hilliard equation with logarithmic free energy, Numer. Math., 72, pp. 1-20, (1995)

[5]

Barrett J., Blowey J., An error bound for the finite element approximation of a model for phase separation of a multi-component alloy, IMA J. Numer. Anal., 16, pp. 257-287, (1996)

[6]

Barrett J., Blowey J., Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy, Numer. Math., 77, pp. 1-34, (1997)

[7]

Barrett J., Blowey J., Finite element approximation of a model for phase separation of a multi-component alloy with a concentration-dependent mobility matrix, IMA J. Numer. Anal., 18, pp. 287-328, (1998)

[8]

Barrett J., Blowey J., Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility, Math. Comput., 68, pp. 487-517, (1999)

[9]

Barrett J., Blowey J., An improved error bound for a finite element approximation of a model for phase separation of a multi-component alloy with a concentration dependent mobility matrix, Numer. Math., 88, pp. 255-297, (2001)

[10]

Baskaran A., Hu Z., Lowengrub J., Wang C., Wise S.M., Zhou P., Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation, J. Comput. Phys., 250, pp. 270-292, (2013)

← 1 2 3 4 5 6 7 →