吉林大学学报(理学版) ›› 2023, Vol. 61 ›› Issue (5): 1063-1072.

• • 上一篇    下一篇

黏性Cahn-Hilliard方程的二阶BDF数值格式

郭媛, 王旦霞, 张建文   

  1. 太原理工大学 数学学院, 太原 030024
  • 收稿日期:2022-11-18 出版日期:2023-09-26 发布日期:2023-09-26
  • 通讯作者: 王旦霞 E-mail: 2621259544@qq.com

Second Order BDF Numerical Scheme for Viscous Cahn-Hilliard Equation

GUO Yuan, WANG Danxia, ZHANG Jianwen   

  1. College of Mathematics, Taiyuan University of Technology, Taiyuan 030024, China
  • Received:2022-11-18 Online:2023-09-26 Published:2023-09-26

摘要: 采用有限元方法对黏性Cahn-Hilliard方程进行数值求解. 首先, 引入辅助变量Lagrange乘子r, 得到黏性Cahn-Hilliard方程的等价形式; 其次, 在空间上采用混合有限元逼近, 时间上采用隐式向后差分公式(BDF)进行离散, 给出黏性Cahn-Hilliard方程的二阶线性有限元数值格式, 并分析所给格式的无条件能量稳定性和误差估计;最后, 通过一系列数值算例验证所给格式的精确性和有效性. 结果表明, 该数值格式是理想的, 并具有同时满足线性、 无条件能量稳定和二阶精度的特点.

关键词: 黏性Cahn-Hilliard方程, Lagrange乘子, 向后差分公式(BDF), 无条件能量稳定

Abstract: We used  finite element method to numerically solve the viscous Cahn-Hilliard equation. Firstly, the equivalent form of the viscous Cahn-Hilliard equation was obtained by introducing the Lagrange multiplier r of the auxiliary variable. Secondly, the second order linear finite element numerical scheme for the viscous Cahn-Hilliard equation was given by using the mixed finite element approximation  in space and the implicit backward differentiation formula (BDF)  for discretization in time, and the unconditional stability in energy and error estimation of the given scheme were analyzed in detail. Finally, a series of numerical examples were used to verify the accuracy and effectiveness of the given scheme. The results show that the proposed numerical scheme is ideal and has the characteristics of simultaneously satisfying linear, unconditional stability in energy and second order accuracy.

Key words: viscous Cahn-Hilliard equation, Lagrange multiplier, backward differentiation formula (BDF), unconditional stability in energy

中图分类号: 

  • O221.6