解抛物问题的一类新的瀑布型多重网格法.pdfVIP

重网格法． 目前它已越来越受到科学计算工作者的重视，被广泛应用于求解线性和非线 性偏微分方程问题．与其它迭代法相比，多重网格法的最大优点是它的计算工作量 与未知量的总数同阶．因此，对大型问题它能成千上万倍地减少计算工作量，大大 地加快收敛速度． 瀑布型多重网格法是一类新的多重网格法，不要求粗网格校正，故又称单步 多重网格法．它产生于二十世纪九十年代中期，并得到迅速发展．Bornemann和 Richardson迭代法作为光滑子，则该法是最优的，但对二维问题仅是拟最优的．若 采用共轭斜量法作为光滑子(见文献[201)，则对二维和三维问题都是最优的．石钟 慈和许学军把上述结论推广到更一般的协调和非协调有限元(见文献[21—221)，考 察了板问题的瀑布型多重网格法(见文献[231)，证明了此时若采用共轭斜量法作 为光滑子，该法是拟最优的，而传统的光滑子不能使用，还讨论了抛物问题(见 G等将此方法应用于求解非线性椭圆边值问题(见文献 【251)．Timmermann 【26。27】)．曾金平和姜英军将此方法应用于求解非对称椭圆边值问题(见文献[28】)． 他们提供的数值实验结果表明，瀑布型多重网格法应用于上述不同问题都是有效 的． 从多重网格法的理论可知，粗网格上定义的有限元空间仅起到校正的作用， 对一些工程上常见的复杂有限元离散，没有必要在粗网格上使用同样复杂的有限 元空间，可以使用比较简单的有限元空间作校正，从而减少计算工作量．石钟慈和 许学军对椭圆问题提出了一类新的瀑布型多重网格法(见文献[29】)，该方法的要 点是在粗网格和最细网格上使用不同的有限元空间，利用粗网格上简单的有限元 空间去处理最细网格上复杂的有限元离散．将这一类新的瀑布型多重网格法推广 到抛物问题有一定的理论和实践意义． 本文推广石钟慈、许学军提出的一类新的瀑布型多重网格法到抛物问题，建 立相应的理论结果．主要研究这一新算法的最优性，以及新算法对一些二阶问题的 应用． 本文以下内容组织如下：第1章，介绍我们所讨论的模型问题及其离散形式， 构造可直接用于计算机编程的算法，并建立相应的收敛性理论和误差估计，讨论 算法的最优性，这是对算法计算性态的理论研究；第2章，将这一新方法应用于 性，得出了相应的理论结果；第3章，给出数值实验结果，以验证算法的有效性． 第1章问题与算法 本章我们简要介绍所讨论的模型问题及兵离散形式，构造可直接用于计算机 编程的新的瀑布型多重网格算法(NCMG)；对网格转移算子、磨光算子作了适当的 假设，讨论新算法的最优性，得到了关于收敛阶和工作量估计的一般性定理． 1．1模型问题及其有限元逼近 本文考虑如下模型问题： f峨一△“；，G，rl 0，f)∈Qx[o，r】 {“G，t)．0， 0，f)∈aQ×【o，r】 (1．1．1) IHb，o)一UoGl z∈Q “∈日；(Q)，“G，0)=U00)使 0。，v)+n0，v)t(，，vlVveH：晒lf∈[o，T】 (1．1．2) 其中 Ⅱ㈧2P眦， (，，y)。11f触 取t。=0 离散格式：求““∈H：心)，使 ． (1．1．3) 巧1仁n，v)+口-”，v)一(g．，y)，Vv∈月：(Ql刀-1,…，N． 其中 g，一f，+f?n“。． 于是求解(1．1．3)的问题化为求解带参数f的椭圆问蕊：求w∈日：(Q)，使 f一1(w，v)+口(w，v)-(g，v)，Vp∈日：(Q)． (1．1．4) 作拟一致剖分L，不妨设^，=hH／2．例如，对L。三角形单元，连接三边中点 即得四个1-i三角形单元；对rf_l矩形单元，连接对边中点即得四个0矩形单元．设巧 为r，对应的有限元空间，巧