针对扩散反应方程在各向同性和各向异性网格上的USFEM方法和段氏方法.pdfVIP

第一置绪论 第一章 绪论 1．1稳定方法的简介 本节我们简单地回顾稳定方法。稳定方法就是在标准Galerkin方法的基 础上加上与网格单元有关的、相容的、能够使格式保持数值稳定的变分项 所得的方法。此类方法的优点在于通过付出较小的代价，保证方法数值 稳定性以及收敛性。在稳定方法的引进和推广方面，tJ．R．Hughes教授做出 主要的贡献。早在1982年的文章[33】中，针对不可压缩的Navier．Stokes方程， formulationsJ。 uDwind／Petrov-Gaterkin 他给出了著名的SUPG格式(Streamline 法(Gelarkin—Least．Sqilaresmethod)，并将其应用到几种结构和流体问题上去。 在八十年代，工J．R．Hughes教授做出不凡的努力，在许多领域都有开创性 的工作。在有关声学的工作[381139]中，稳定方法付出较小代价，而得到方法 的收敛性。以前这个领域中，有限元方法要得到一定精度所付出的代价很高。 所以稳定方法有很大的应用前景。80年代有关可压缩的Navier．Stokes方程的工 作[4J0—“1，也为九十年代一直到现在有关稳定方法的许多发展提供了动力和来 源。 Brezzi． 进入九十年代，对稳定方法的认识有很大地提高。在【45]中，F M．O．Bristeau和LPFranca做了很重要的工作，他们讨论了稳定方法和 Galerkinmetllodwithbubble 带Bubble函数的Galerkin方法(the functions)之间的关 系，他们认为通过静凝聚方法(static 法可以导出一种稳定方法。这是后面寻求理想的bubble functions的工作[46— 系．这样对稳定方法的来源的认识又提高到一个新阶段。同时针对不同 的模型的稳定方法也有很大的发展，例如，针对对流扩散方程，反应扩 散方程，以及反应对流扩散方程．有一系列的工作，看[4】[5][50．57]。在文 方法，本文中的许多资料摘于此。 ’ 除了上面所提到的SUPG和GLS这两种稳定方法以外．还有其他多种稳定 第一章绪论 法(theUnusualStabilizedFiniteElement Method)和段火元在[8]给出的一种稳定方 法。对稳定方法来说，稳定参数的选择是一件重要而且困难的事．它对稳定 定参数改进了[414中针对对流扩散反应方程的USFEM方法。通过USFEM方法和 with 带Bubble函数的Galerkia方法(theGalerkinmethodhubMe缸蜮j∞s)之间的关 系，我们来确定在各向异性网格上稳定参数的选择，进而将它们推广到各向异 性网格上去。 第一量绪论 1．2自适应方法的简介 自适应算法是根据中间计算结果自动控制计算过程的求偏微分方程的方 法。它主要利用中间计算结果自动计算所需的网格，选取最佳的离散方式，从 而逐步对误差自动地作适当调节以达到所需精度。该方法具有较高的识别能力 和选择最优参数的能力，一尽量少的计算量达到所要求的精度。有限元自适应 算法在最近十几年里发展最快。由于它能够自动判断在结构误差最大处加密网 格，提高了计算效率．并且通过误差估计可以告诉用户计算结果的误差范围。 因此，它在工程界受到了广泛的欢迎。 最近十几年所进行的有限元自适应方法的研究取得了重要进展，自适应分 析过程主要包括两个方面一可靠的误差估计方法和功能强大的网格自动剖分技 术。