论文部分内容阅读
本文的工作分两部分。第一部分主要研究了一类基于局部多项式逼近空间的单位分解法的最优误差估计。近年来,无网格方法被大量应用到科学与工程计算中。与经典有限元方法相比较,这类方法的共同特征是不再需要网格结构,它们在处理具有复杂域的问题或区域在求解过程中变化的问题时非常有效。单位分解法是非常重要的无网格方法之一,其两大主要特点是:一方面,它允许使用支集不依赖于网格或依赖于不与问题域一致的简单网格的单位分解函数(例如Shepard函数[4]),在此意义下,PUM是一种无网格方法,这个特征免去了网格生成;另一方面,局部逼近空间可以包含非多项式函数,从而很好地局部逼近未知解。这两大特征使单位分解法得到了迅速发展。虽然有关单位分解法的文献较多,但大部分都侧重于工程应用,只有极少数的数学理论分析,主要以I.Babu(?)ka和他的合作者为代表做了很多奠基性工作。但I.Babu(?)ka等现有的关于单位分解法的插值误差估计还没有获得最优阶,本文第一部分的主要工作就是:通过构造一种特殊的局部多项式近似空间,以获得最优阶插值误差估计。为此,作者从有限元方法的误差收敛阶入手,针对一类特殊的单位分解方法(取通常的有限元基函数作单位分解)进行分析,构造了一个特殊的局部多项式逼近空间,给出了一维下高次单位分解插值格式和二维下低次单位分解插值格式,推导了相应的最优阶插值误差,并研究了一维下Galerkin解的误差估计。第二部分主要研究了求解一类椭圆型变分不等式的修正自适应代数多重网格法及其并行化。根据离散的椭圆型变分不等方程所具有的线性互补性质,提出了一个基于积极集策略之上的修正代数多重网格解法,求解具有对称二阶椭圆算子的变分不等式的有限元离散问题。数值实验表明了该算法在一致网格和h-自适应网格上的计算有效性和健壮性。为了减少计算时间,本文还根据该修正算法内在的并行度,提出了一个并行计算格式,数值结果给出了该并行的加速比和效率。