对流扩散特征值问题的一种多水平校正方法-计算数学专业论文.docxVIP

I万方数据 I 万方数据 摘 要 本文提出一种对流扩散特征值问题的多水平校正方法, 这是一种解对流扩散特征 值问题的新型, 有效的方法.它既适用于单特征值也适用于重特征值. 这种方法将在细网 格上解特征值问题归结为在粗网格上解特征值问题以及在加密后的网格上解边值问题. 文章证明了在每次校正步后特征对的逼近结果的精度能得到提高, 并且还给出了这个 方案能达到与解相应边值问题几乎相同的最优精度的理论证明, 同时, 这个方法是容易 在Matlab 软件下实现的, 并且实验的结果是令人满意的. 关键词: 对流扩散特征值问题; 有限元; 多水平校正; Matlab. II II 万方数据 Abstract In this paper, we propose a multilevel correction method for the convection-diffusion eigenvalue problems, and the method is a new type and efficiency to solve the eigenvalue problems, which is suitable for not only simple but also multiple eigenvalues. In this method, solving eigenvalue problem is transformed to a series of solutions of boundary value problem by the multilevel method on the finest meshes and a serious of eigenvalue problems on the coarsest finite element space. we prove that the accuracy of resulting eigenpair approximations can be improved after each correction step, and in this paper, we give the theoretical proof of this scheme that almost can achieve the optimal accuracy as same as solving the corresponding boundary value problem. Our scheme is easy to realize with Matlab, and numerical results are satisfactory . Key Words: convection-diffusion eigenvalue problems; finite element; multi- level correction; Matlab. PAGE 1 PAGE 1 万方数据 1 绪论 对流扩散特征值问题有着很重要的物理背景, 比如在流体力学和环境问题下的对流 扩散, 因此, 用有限元方法解这种问题变成了一个重要的课题, 并且这也吸引了力学以及 物理学领域众多学者的注意力. 参考文献[1] 讨论了一种后验误差估计, [2–7] 讨论了自适 应算法, 文献[8]讨论的是自适应同伦方案，文章[9]讨论了二水平算法, 文献[10]是讨论了 函数值还原算法, 文章[11, 12]讨论了外推法等. 这篇文章讨论基于林群和谢和虎的校正法 （见文章[13, 14]）的有限元多水平离散, 即将此方法应用于对流扩散特征值问题的求解. 在 文 章[13, 14]中, 林 群 和 谢 和 虎 引 入 了 一 种 新 型 的 多 水 平 校 正 方 案, 随 后, 这 种 校正方案得到更进一步的发展并被成功的应用于Steklov特征值问题（见文章[15]）, Helmholtz传输特征值问题（见文章[16]）等. 在这篇文章中, 作者应用林群和谢和虎的方法到对流扩散特征值问题, 得到一种多 水平方法. 这种方法描述如下： 第一步, 构造一个粗的有限元空间并在这空间上解原特 征值问题和其共轭特征值问题；第二步, 用上一步中所求得的特征值乘以相应的特征 函数作为边值问题右端, 然后在加密的空间里解两个相关的边值问题；第三步, 将第二 步所获得的近似特征函数和第一步中的最粗有限元空间相结合得到一个新的有限元空 间, 并在其上再一次解原特征值问题和其共轭的特征值问题, 然后转至第二步进行下一 个循环. 这种方法不仅适合于单特征值而且也适用于重特征值. 此外, 文章还给出了这个 方法能达到与解相应边值问题相同的最优收敛阶的理论