A Second Order Accurate in Time, Energy Stable Finite Element Scheme for the Flory-Huggins-Cahn-Hilliard Equation

A Second Order Accurate in Time, Energy Stable Finite Element Scheme for the Flory-Huggins-Cahn-Hilliard Equation

Year:    2022

Author:    Maoqin Yuan, Wenbin Chen, Cheng Wang, Steven M. Wise, Zhengru Zhang

Advances in Applied Mathematics and Mechanics, Vol. 14 (2022), Iss. 6 : pp. 1477–1508

Abstract

In this paper, we propose and analyze a second order accurate in time, mass lumped mixed finite element scheme for the Cahn-Hilliard equation with a logarithmic Flory-Huggins energy potential. The standard backward differentiation formula (BDF) stencil is applied in the temporal discretization. In the chemical potential approximation, both the logarithmic singular terms and the surface diffusion term are treated implicitly, while the expansive term is explicitly updated via a second-order Adams-Bashforth extrapolation formula, following the idea of the convex-concave decomposition of the energy functional. In addition, an artificial Douglas-Dupont regularization term is added to ensure the energy dissipativity. In the spatial discretization, the mass lumped finite element method is adopted. We provide a theoretical justification of the unique solvability of the mass lumped finite element scheme, using a piecewise linear element. In particular, the positivity is always preserved for the logarithmic arguments in the sense that the phase variable is always located between −1 and 1. In fact, the singular nature of the implicit terms and the mass lumped approach play an essential role in the positivity preservation in the discrete setting. Subsequently, an unconditional energy stability is proven for the proposed numerical scheme. In addition, the convergence analysis and error estimate of the numerical scheme are also presented. Two numerical experiments are carried out to verify the theoretical properties.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/aamm.OA-2021-0331

Advances in Applied Mathematics and Mechanics, Vol. 14 (2022), Iss. 6 : pp. 1477–1508

Published online:    2022-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    32

Keywords:    Cahn-Hilliard equations Flory Huggins energy potential mass lumped FEM convex-concave decomposition energy stability positivity preserving.

Author Details

Maoqin Yuan

Wenbin Chen

Cheng Wang

Steven M. Wise

Zhengru Zhang

  1. Explicit P1 Finite Element Solution of the Maxwell-Wave Equation Coupling Problem with Absorbing b. c.

    Beilina, Larisa | Ruas, Vitoriano

    Mathematics, Vol. 12 (2024), Iss. 7 P.936

    https://doi.org/10.3390/math12070936 [Citations: 0]
  2. Numerical simulation of acoustic wave propagation by finite element method based on optimized matrices

    Li, Lei | Wen, Xiaotao | Tang, Chao | Zhou, Dongyong | Zhang, Songgen

    Journal of Geophysics and Engineering, Vol. 21 (2024), Iss. 3 P.1027

    https://doi.org/10.1093/jge/gxae055 [Citations: 0]
  3. Numerical Analysis of a Convex-Splitting BDF2 Method with Variable Time-Steps for the Cahn–Hilliard Model

    Hu, Xiuling | Cheng, Lulu

    Journal of Scientific Computing, Vol. 98 (2024), Iss. 1

    https://doi.org/10.1007/s10915-023-02400-5 [Citations: 0]
  4. Energy dissipative and positivity preserving schemes for large-convection ion transport with steric and solvation effects

    Ding, Jie | Wang, Zhongming | Zhou, Shenggao

    Journal of Computational Physics, Vol. 488 (2023), Iss. P.112206

    https://doi.org/10.1016/j.jcp.2023.112206 [Citations: 3]
  5. Convergence analysis of a second order numerical scheme for the Flory–Huggins–Cahn–Hilliard–Navier–Stokes system

    Chen, Wenbin | Jing, Jianyu | Liu, Qianqian | Wang, Cheng | Wang, Xiaoming

    Journal of Computational and Applied Mathematics, Vol. 450 (2024), Iss. P.115981

    https://doi.org/10.1016/j.cam.2024.115981 [Citations: 3]
  6. A third-order positivity-preserving and energy stable numerical scheme for the Cahn-Hilliard equation with logarithmic potential

    Yuhuan, Li | Jianyu, Jing | Qianqian, Liu | Cheng, Wang | Wenbin, Chen

    SCIENTIA SINICA Mathematica, Vol. (2024), Iss.

    https://doi.org/10.1360/SSM-20223-0014 [Citations: 0]
  7. Convergence and stability analysis of energy stable and bound‐preserving numerical schemes for binary fluid‐surfactant phase‐field equations

    Duan, Jiayi | Li, Xiao | Qiao, Zhonghua

    Numerical Methods for Partial Differential Equations, Vol. 40 (2024), Iss. 6

    https://doi.org/10.1002/num.23125 [Citations: 0]
  8. Fully discrete scheme for a time-dependent Ginzburg-Landau equation in macromolecular microsphere composite hydrogels

    Hou, Bingrui | Yuan, Maoqin | Huang, Pengzhan

    Computers & Mathematics with Applications, Vol. 151 (2023), Iss. P.127

    https://doi.org/10.1016/j.camwa.2023.09.038 [Citations: 0]
  9. EnVarA-FEM for the flux-limited porous medium equation

    Liu, Qianqian | Duan, Chenghua | Chen, Wenbin

    Journal of Computational Physics, Vol. 493 (2023), Iss. P.112432

    https://doi.org/10.1016/j.jcp.2023.112432 [Citations: 0]
  10. Optimal rate convergence analysis of a numerical scheme for the ternary Cahn–Hilliard system with a Flory–Huggins–deGennes energy potential

    Dong, Lixiu | Wang, Cheng | Wise, Steven M. | Zhang, Zhengru

    Journal of Computational and Applied Mathematics, Vol. 415 (2022), Iss. P.114474

    https://doi.org/10.1016/j.cam.2022.114474 [Citations: 4]
  11. An unconditionally energy stable linear scheme for Poisson–Nernst–Planck equations

    Qiao, Tian | Qiao, Zhonghua | Sun, Shuyu | Zhou, Shenggao

    Journal of Computational and Applied Mathematics, Vol. 443 (2024), Iss. P.115759

    https://doi.org/10.1016/j.cam.2024.115759 [Citations: 3]
  12. A Positivity-Preserving, Energy Stable BDF2 Scheme with Variable Steps for the Cahn–Hilliard Equation with Logarithmic Potential

    Liu, Qianqian | Jing, Jianyu | Yuan, Maoqin | Chen, Wenbin

    Journal of Scientific Computing, Vol. 95 (2023), Iss. 2

    https://doi.org/10.1007/s10915-023-02163-z [Citations: 5]
  13. Maximum bound principle and non-negativity preserving ETD schemes for a phase field model of prostate cancer growth with treatment

    Huang, Qiumei | Qiao, Zhonghua | Yang, Huiting

    Computer Methods in Applied Mechanics and Engineering, Vol. 426 (2024), Iss. P.116981

    https://doi.org/10.1016/j.cma.2024.116981 [Citations: 1]
  14. A second-order, mass-conservative, unconditionally stable and bound-preserving finite element method for the quasi-incompressible Cahn-Hilliard-Darcy system

    Gao, Yali | Han, Daozhi | Wang, Xiaoming

    Journal of Computational Physics, Vol. 518 (2024), Iss. P.113340

    https://doi.org/10.1016/j.jcp.2024.113340 [Citations: 0]