论文部分内容阅读
瀑布型多重网格法是求解大型边值问题的一种有效迭代解法.其主要的优点是不要求粗网格校正,故又称单步多重网格法.Gisela Timmermann用瀑布型多重网格法对半线性椭圆问题进行了求解,在粗网格上用牛顿法(Newton)将由线性有限元离散而得到的非线性方程组线性化,在细网格上用瀑布型多重网格法解这个牛顿方程,并且提出了算法和对算法的收敛性进行了研究.该文将瀑布型多重网格法推广到半线性抛物问题,证明在能量范数下算法误差的最优阶,且有最优或有拟最优的计算工作量,并进行数值试验.该文以半线性二阶抛物型偏微分方程初边值问题为模型问题构造了瀑布型多重网格法(CMG),在最粗网格上用牛顿法(Newton)将由线性有限元离散而得到的非线性方程组线性化,在细网格上用瀑布型多重网格法解这个牛顿方程首先采用Richardson迭代法作为光滑子,我们证明了瀑布型多重网格法对二维半线性抛物型边值问题在能量范数下可获得最优收敛阶.然后又对采用共轭梯度法(CG)作为光滑子进行了研究,并得到了在此情况下,瀑布型多重网格法对二维半线性抛物型边值问题在能量范数下可获得最优收敛阶.同时对这两种情形,分析了计算工作量,得到了工作量的最优性或拟最优性,这说明在半线性情形的计算工作量与线性情形是大致相当的.数值实验也显示了该算法的有效性.