晶体相场方程的高精度 Fourier 拟谱逼近

摘要

针对能量稳定和质量守恒的晶体相场方程的周期边界条件,提出一种具有新的能量稳定性的高精度数值算法。算法对方程的离散在空间上采用Fourier拟谱逼近,在时间上采用三阶精度的后向微分形式,并增加Douglas-Dupont项,以保证数值格式的能量稳定性。理论上证明了数值解的存在唯一性以及格式的能量稳定性。数值仿真验证了算法的高精度和稳定性。