Cahn-Hilliard方程多重网格求解器收敛性分析

Cahn-Hilliard(CH)方程是相场模型中的一个基本的非线性方程，通常使用数值方法进行分析。在对CH方程进行数值离散后会得到一个非线性的方程组，全逼近格式(Full Approximation Storage, FAS)是求解这类非线性方程组的一个高效多重网格迭代格式。目前众多的求解CH方程主要关注数值格式的收敛性，而没有论证求解器的可靠性。文中给出了求解CH方程离散得到的非线性方程组的多重网格算法的收敛性证明，从理论上保证了计算过程的可靠性。针对CH方程的时间二阶全离散差分数值格式，利用快速子空间下降(Fast Subspace Descent, FASD)框架给出其FAS格式多重网格求解器的收敛常数估计。为了完成这一目标，首先将原本的差分问题转化为完全等价的有限元问题，再论证有限元问题来自一个凸泛函能量形式的极小化，然后验证能量形式及空间分解满足FASD框架假设，最终得到原多重网格算法的收敛系数估计。结果显示，在非线性情形下，CH方程中的参数ε对网格尺度添加了限制，太小的参数会导致数值计算过程不收敛。最后通过数值实验验证了收敛系数与方程参数及网格尺度的依赖关系。