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

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

作　　者：霍俊蓉 刘昊 温学兵[1] 张荣培 蔚喜军[3] HUO Junrong;LIU Hao;WEN Xuebing;ZHANG Rongpei;WEI Xijun(College of Mathematics and Systems Science,Shenyang Normal University,Shenyang 110034,China;College of Applied Mathematics,Guangdong University of Technology,Guangzhou 510006,China;Institute of Applied Physics and Computational Mathematics,Beijing 100088,China)

机构地区：[1]沈阳师范大学数学与系统科学学院,沈阳110034 [2]广东工业大学应用数学学院,广州510006 [3]北京应用物理与计算数学研究所,北京100088

出　　处：《吉林大学学报（理学版）》2022年第3期721-728,共8页Journal of Jilin University:Science Edition

基　　金：辽宁省自然科学基金(批准号:20180550996);国防科技重点实验室基金(批准号:6142A0502020717).

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

关 键 词：CAHN-HILLIARD方程 有限差分方法 Crank-Nicolson方法 能量耗散 

分 类 号：O411.1[理学—理论物理]