边值问题离散系统数值稳定性的新度量

对于用有限元法或差分法求解微分方程边值问题,流行着这样的观点:当网格剖分出现小角度的三角形或窄长的矩形时,离散系统的数值稳定性就差.这种观点的根据是由于把系数矩阵条件数作为稳定性度量.     

凹角域应力强度因子的外推加速及MG算法

§1.引言 [1]最早讨论将外推用于嵌套迭代,[2]-[4]则讨论外推用于多重网格法,两者都没有涉及凹角域的情况.在凸域上,有渐近展式(例如[5]): u~h(x)=u~I(x)+d_1(x)h~2+O(h~τ),x∈Ω,(1.1)其中,τ> 2,u~h和u~I分别为椭圆边值问题解u的线性有限元逼近和线性插值函数.而     

矩形域低光滑解的渐近展开

研究外推法用于偏微分方程边值问题数值解的效果及理论,已有[1—6].[5,6]对有限元任意初始剖分证明了渐近展开,使外推理论向前迈出了新的一步.进一步的问题是:就实际情况而言,已有理论对光滑性要求仍然过高.在高光滑解的情况下,高次元可达到     

有限元方法的后验误差估计

有限元的误差分析理论已有丰富的成果,但一般以先验估计为主。这些理论对算法评价和指导计算无疑具有重要意义。但实际计算中常要求给出误差的定量估计,先验估计一般就做不到这一点,因为其中会有如微分方程解的导数等不可计算的量。后验误差估计的研究近年来已开始受到注意。在多重网格技术中基于误差渐近展开的外推可以给出逐点意义下的后验估计。(例如见[1],[2])另一条途径是通过局部化进行能量模意义下的后验误差估计。其中对一维边值问题的研究已较成熟。[7]就常系数