吉林大学学报(理学版) ›› 2022, Vol. 60 ›› Issue (3): 721-728.

• • 上一篇    下一篇

一类相场方程的能量稳定性分析及数值模拟

霍俊蓉1, 刘昊1, 温学兵1, 张荣培2, 蔚喜军3   

  1. 1. 沈阳师范大学 数学与系统科学学院, 沈阳 110034; 2. 广东工业大学 应用数学学院, 广州 510006; 3. 北京应用物理与计算数学研究所, 北京 100088
  • 收稿日期:2021-06-21 出版日期:2022-05-26 发布日期:2022-05-26
  • 通讯作者: 张荣培 E-mail:rongpeizhang@163.com

Energy Stability Analysis and Numerical Simulation of a Class of Phase-Field Equations

HUO Junrong1, LIU Hao1, WEN Xuebing1, ZHANG Rongpei2, WEI Xijun3   

  1. 1. College of Mathematics and Systems Science, Shenyang Normal University, Shenyang 110034, China; 2. College of Applied Mathematics, Guangdong University of Technology, Guangzhou 510006, China; 3. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
  • Received:2021-06-21 Online:2022-05-26 Published:2022-05-26

摘要: 提出一种快速、 稳定的数值方法求解具有恒定迁移率的二维Cahn-Hilliard方程. 在空间离散上采用二阶有限差分方法, 在时间离散上采用Crank-Nicolson方法, 从理论上证明离散能量随时间发展具有耗散性质. 针对全离散格式下的非线性代数方程组, 应用不动点迭代方法求解, 并利用快速离散余弦变换(FDCT)以提高计算效率. 数值实验结果表明, 离散自由能关于时间是非递增的, 该方法具有稳定性好、 存储量小、 计算速度快等优点.

关键词: Cahn-Hilliard 方程, 有限差分方法, Crank-Nicolson方法, 能量耗散

Abstract: We proposed a fast and stable numerical method to solve the two-dimensional Cahn-Hilliard equation with constant mobility. The second order finite difference method was used in spatial discretization and Crank-Nicolson method was used in time discretization. We proved theoretically that the discrete energy had the property of dissipation with time evolving. The fixed point iteration method was used to solve the nonlinear algebraic equations in the fully discrete scheme, and the fast discrete cosine transform (FDCT) was used to improve the computational efficiency. The numerical results show that the discrete free energy is non increasing with respect to time, and the method has the advantaes of good stability, small storage and fast computation speed.

Key words: Cahn-Hilliard equation, finite difference method, Crank-Nicolson method, energy dissipation

中图分类号: 

  • O411.1