论文部分内容阅读
自适应有限元方法自上世纪七十年代以来一直都是科学与工程计算领域的研究热点,是一种以常规有限元方法为基础通过后验误差估计和网格自动调整来不断提高逼近解精度的数值计算方法。它能以尽量少的计算量达到所要求的逼近精度,因此是一种具有较高识别能力,应用范围更广,高效可靠的计算方法。 其中,后验误差估计是自适应有限元方法的关键和核心,其作用就是利用有限元逼近解的信息和已知量计算逼近解在求解域当前剖分网格上的逼近误差,得到逼近解误差的分布情况,进而对其逼近精度做出可靠性分析。后验误差估计方法的好坏直接影响着自适应有限元方法的效率和可靠性。网格自动调整指的是,根据误差均分原则,利用得到的后验误差估计量及其分布,对当前的网格实施剖分加密或放粗的调整,进而给出下一步有限元计算的网格。 本文针对变系数椭圆问题,给出了一种新型的后验误差估计方法。它的主要思想是,运用有限元方法在求解域当前的剖分网格Τh求解问题,得到问题的有限元逼近解uh;接着,全局一致加密Τh得到辅助网格Τh/2,在辅助网格上,首先利用网格Τh上的单元刚度矩阵信息快速形成Τh/2上的总刚度矩阵,然后,以uh在辅助网格对应的有限元空间上的插值为初值,利用光滑迭代法计算辅助网格上的有限元方程,得到问题的一个近似解uh/2,m;最后,以uh和uh/2,之差的能量范数作为后验误差估计量。基于此后验误差估计方法,提出了相应的自适应有限元算法。另外,对于自适应计算过程中的网格调整,运用的是能保证网格正则性和协调性,且易于实现二分法网格加密技术。 实践中,进行的多种类型的数值实验结果证实,本文给出的后验误差估计和自适应有限元算法是有效可行的。