粘性Cahn-Hilliard方程的二阶稳定Crank-Nicolson-Leapfrog格式

本文研究了粘性Cahn-Hilliard方程的一个二阶逼近。对于双阱势函数,本文通过引入拉格朗日乘子得到一个等价形式。其次,使用Crank-Nicolson-Leapfrog格式进行时间离散,使用有限元方法进行空间离散,从而得出了一个二阶线性无条件稳定的数值格式。然后,证明了数值格式的无条件能量稳定性和误差分析。最后,给出了几个数值模拟,验证了格式的数值精度。