外推多网格法——-大规模科学计算的新算法.pdf

摘 要 快速求解由差分法或有限元法导出的大型线性方程组是大规模科 学和工程计算中的重要问题．人们的目标是对Ⅳ阶方程组用O(N)次 乘除运算得到所需精度的解．而多网格法(MG)第一次实现了这个目 标，成为求解大规模问题最有效的方法．经典MG的理论方法已经比 较成熟，由于使用了3种运算：插值，限制和迭代，求解时需要在多重 粗细网格上反复迭代，所以程序较为复杂．1996年德国Bornemann等人 提出瀑布式多网格法(CMG)，即从粗网到细网的单向计算，只采用了 插值与迭代两种运算，程序容易实现令人注目． 本论文研究一类新的瀑布式多重网格法，主要工作和创新点如下． 首次提出了外推瀑布式多网格法(EXCMG)．本论文基于有限元的 渐近展开式，导出了新的外推公式．新算法沿用CMG的思想，但将粗 网上的线性插值改为新外推与二次插值，为密网提供更好的初值，本 质上减小了初始误差，对加速收敛起着关键的作用．新算法对函数和 导数都收敛，且都有高精度．本论文在PC机上用EXCMG求解了400 万未知数的问题，对函数和导数都达到8-10位精度，仅费时lo分钟，进 一步证实了这些优势． 首次证明了共轭梯度法(CG)按f2模的有界性和收敛性．对CG法 至今只有按能量模的收敛估计，因此这是一个新结果．本文还发现当 第i层网格的迭代次数取幂指数形式(mLp仁t)时，基于CG迭代的多 重网格法，存在着—个网格层数的“阀值”io，即当ito时，第i层上的 误差快速衰减CG迭代压缩效果非常明显；而当i硒时，CG主要起 磨光作用．此特性对证明EXCMG按离散L2模的收敛性起重要作用． 对正则性较弱的解u∈H3，首次证明了双线性元的渐近展开在节 点集上按离散Lz模有高精度O(h3)，因此EXCMG仍有效．已有的渐近 展开式研究主要是对光滑解(如牡∈C。)，并只得到了在逐点意义下的 结果．上述结果拓宽了EXCMG的应用范围．数值试验表明，EXCMG 对非光滑解u∈H3甚至H2仍有很高精度． 关键词：有限元；新外推公式；二次插值；共轭梯度法；瀑布式多 重网格法 ABSTRACT solvethe of derivedfinitedifference Howto linear systemequationsby in methodorfiniteelementmethod isan issuelarge-scale rapidlyimportant and istoobtainthesolution scientific engineeringcomputing．People’Sgoal with ofN-order system multiplication requiredaccuracy throughO(N)times anddivision this forthefirst operations．Multi-gridmethod(MC)achievesgoal tosolve becomesthemosteffectivemethod problems． time，and large-scale Since methodhasbeen asystematic 1970’S，themulti-grid