一维有限元的h-p方法的后验误差估计

有限元的h-p方法,是指在增加有限元空间的维数时,既加密某些单元的网格,同时也增加某些单元的次数.对h-p方法,人们希望得到O(h~mp~(-n))(m,n>0)形状的误差估计.这种误差估计的结果包括了对传统的h方法以及p方法的结果.关于h-p方法的     

Wilson元插值误差渐近估计

Wilson元是工程界常用的一种有限元计算方法,但在理论分析中插值误差估计的常数只知道存在,不知道具体值.本文给出了在L~2、H~1范数意义下Wilson元在参考单元和一般单元上插值误差渐近估计,导出了主要常数.这种精确的估计为有限元后验误差估计和自适应计算提供保障.     

双奇次有限元的渐近准确误差估计

§1.引言 近年来自适应有限元方法无论在数学理论还是在实际应用方面都已得到迅速发展.I.Babuska 等首先提出了双线性单元(p=1)的h型自适应方法.此后作者与Babuska又发展了双二次单元(p=2)的h型自适应方法并进行了一系列数值计算.这些成果已被应用于美国马里兰大学的自适应有限元程序FEARS中.自适应方法的基础在于对有限元近似解作后验误差估计,这些估计应是便于计算的.作者在[5]中已对任     

细菌模型的非协调有限元分析

将非协调元应用于描述细菌传播的反应扩散方程组的初边值问题.借助单元的一些特性和非协调误差估计技巧,分别在半离散和全离散有限元格式下,研究了其数值解与精确解的误差估计,得到了最优的误差估计以及超逼近结果.     

一个新非协调单元对扩散对流反应方程的应用

利用最近提出的一个新型非协调双参数单元,将流线扩散有限元方法成功地应用于对流占优的扩散对流反应方程,并且得到流线扩散模意义下的误差估计结果.     

非协调有限元逼近的梯度恢复型后验误差估计（英文）

本文给出二阶椭圆型方程的非协调有限元的梯度恢复型后验误差估计.后验误差估计是在Crouzeix-Raviart非协调有限单元上得到的,并且给出误差的上下界,更进一步可以证明所得的后验误差估计在拟一致网格上是渐近精确的,所以误差估计是可行的、有效的.上界证明过程依赖于"Helmholtz分解",下界证明主要依赖"bubble函数".数值结果验证了理论的正确性.     

非线性粘弹性方程的EQrot1非协调有限元分析

针对非线性粘弹性方程,在半离散和全离散格式下给出EQrot1非协调有限元逼近.由于该单元的相容误差(O(h2)阶)比插值误差(O(h)阶)高一阶,可得到在H1模意义下的O(h2)阶超逼近结果,并利用插值后处理技术导出整体超收敛.进而,基于该单元的渐近展开式,构造新的插值后处理算子和外推格式,给出O(h4)阶的外推结果.最后,运用与以往文献不同的方法得到全离散逼近格式的最优误差估计.     

非线性粘弹性方程的EQ1rot非协调有限元分析

针对非线性粘弹性方程,在半离散和全离散格式下给出EQ1rot非协调有限元逼近.由于该单元的相容误差 (O(h2)阶)比插值误差 (O(h)阶)高一阶,可得到在H1模意义下的O(h2)阶超逼近结果,并利用插值后处理技术导出整体超收敛.进而,基于该单元的渐近展开式,构造新的插值后处理算子和外推格式,给出O(h4)阶的外推结果.最后,运用与以往文献不同的方法得到全离散逼近格式的最优误差估计.     

平面弹性问题自适应有限元方法的收敛性分析

针对平面弹性问题,首先采用基于最新顶点二分法的网格加密方法,给出一种不需要标记振荡项和加密单元、不需要满足"内节点"性质的自适应有限元方法.其次,通过对各层网格上解函数和误差指示子的分析,利用相邻网格层上解函数的正交性、解函数和真解函数的能量误差的上界估计、相邻网格层上误差指示子的近似压缩性等结果,从理论上严格证明了该自适应有限元方法是收敛的.最后数值实验验证了该自适应有限元方法是收敛的和鲁棒的.     

高次三角形有限元的超收敛问题

关于二维区域二阶线性椭圆问题的有限元求解,[1,2]各自独立地对低次奇妙族矩形元采用单元合并技巧,获得能量的近似正交性(或称插值误差的第一弱估计),从而获得应力佳点定理.若获得更佳形式的能量正交性(或称插值误差的第二弱估计),则可获得位移佳点定理.运用以上方法,[1—8]解决了奇妙族矩形任意次元及三角形线元、二次元     

DOMAIN DECOMPOSITION METHODS WITH NONMATCHING GRIDS FOR THE UNILATERAL PROBLEM

1. IntroductionDomain decomposition methods (DDMs) with nonmatching grids, which have been deveLoped in recent years3 are a quite new class of nonconforming DDMs. As this kind of DDMscan be aPplied to solving many practical problems which can't be handled…     

Sharp upper global a posteriori error estimates for nonlinear elliptic variational problems

The paper is devoted to the problem of verification of accuracy of approximate solutions obtained in computer simulations. This problem is strongly related to a posteriori error estimates, giving computable bounds for computational errors and detecting zones in the solution domain where such errors are too large and certain mesh refinements should be performed. A mathematical model embracing nonlinear elliptic variational problems is considered in this work. Based on functional type estimates developed on an abstract level, we present a general technology for constructing computable sharp upper bounds for the global error for various particular classes of elliptic problems. Here the global error is understood as a suitable energy type difference between the true and computed solutions. The estimates obtained are completely independent of the numerical technique used to obtain approximate solutions, and are sharp in the sense that they can be, in principle, made as close to the true error as resources of the used computer allow. The latter can be achieved by suitably tuning the auxiliary parameter functions, involved in the proposed upper error bounds, in the course of the calculations.     

THE SHIFTED-INVERSE POWER WEAK GALERKIN METHOD FOR EIGENVALUE PROBLEMS

This paper proposes and analyzes a new weak Galerkin method for the eigenvalue problem by using the shifted-inverse power technique. A high order lower bound can be obtained at a relatively low cost via the proposed method. The error estimates for both eigenvalue and eigenfunction are provided and asymptotic lower bounds are shown as well under some conditions. Numerical examples are presented to validate the theoretical analysis.     

ERROR ANALYSIS OF A STABILIZED FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES PROBLEM

This paper is devoted to the establishment of sharper $a$ $priori$stability and error estimates of a stabilized finite element method proposed by Barrenechea and Valentin for solving the generalized Stokes problem, which involves a viscosity $\nu$ and a reaction constant $\sigma$. With the establishment of sharper stability estimates and the help of $ad$ $hoc$finite element projections, we can explicitly establish the dependence of error bounds of velocity and pressure on the viscosity $\nu$, the reaction constant $\sigma$, and the mesh size $h$. Our analysis reveals that the viscosity $\nu$ and the reaction constant $\sigma$ respectively act in the numerator position and the denominator position in the error estimates of velocity and pressure in standard norms without any weights. Consequently, the stabilization method is indeed suitable for the generalized Stokes problem with a small viscosity $\nu$ and a large reaction constant $\sigma$. The sharper error estimates agree very well with the numerical results.     

加权总体最小二乘问题的分析

总体最小二乘问题由Golub和Van Loan首先进行数学的分析,随后人们对于总体最小二乘问题的算法、解的各种形式、总体最小二乘解和最小二乘解的关系、总体最小二乘解的扰动理论以及数值试验作了大量的研究工作。近来,[10]中给出了总体最小二乘问题(TLS)较一般地讨论。另一方面,Golub和Van Loan研究了总体最小二乘问题的特殊均加权形式。本文试图在[10,11]的基础上讨论最一般的总体最小二     

THE DUAL MIXED METHOD FOR AN UNILATERAL PROBLEM

In this paper, the dual mixed method for an unilateral problem, which is the simplified modelling of scalar function for the friction-free contact problem, is considered. The dual mixed problem is introduced, the existence and uniqeness of the solution of the problem are presented, and error bounds O(h^3/4 ) and O(h^3/2 ) are obtained for the dual mixed finite element approximations of Raviart-Thomas elements for k= 0 and k= 1 respectively.     

A posteriori error estimates for a compressible Stokes system

A linearized compressible viscous Stokes system is considered. The a posteriori error estimates are defined and compared with the true error. They are shown to be globally upper and locally lower bounds for the true error of the finite element solution. Some numerical examples are given, showing an efficiency of the estimator. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 412–431, 2004.     

Two-sided a posteriori error estimates for linear elliptic problems with mixed boundary conditions

The paper is devoted to verification of accuracy of approximate solutions obtained in computer simulations. This problem is strongly related to a posteriori error estimates, giving computable bounds for computational errors and detecting zones in the solution domain where such errors are too large and certain mesh refinements should be performed. A mathematical model consisting of a linear elliptic (reaction-diffusion) equation with a mixed Dirichlet/Neumann/Robin boundary condition is considered in this work. On the base of this model, we present simple technologies for straightforward constructing computable upper and lower bounds for the error, which is understood as the difference between the exact solution of the model and its approximation measured in the corresponding energy norm. The estimates obtained are completely independent of the numerical technique used to obtain approximate solutions and are “flexible” in the sense that they can be, in principle, made as close to the true error as the resources of the used computer allow. This work was supported by the Academy Research Fellowship No. 208628 from the Academy of Finland.     

非线性Sobolev方程的特征－差分方法

非线性Ｓｏｂｏｌｅｖ方程的特征－差分方法由同顺（南开大学）ＴＨＥＣＨＡＲＡＣＴＥＲＩＳＴＩＣＤＩＦＦＥＲＥＮＣＥＭＥＴＨＯＤＦＯＲＡＮＯＮＬＩＮＥＡＲＳＯＢＯＬＥＶＥＱＵＡＴＩＯＮ￥ＹｏｕＴｏｎｇ－ｓｈｕｎ（ＮａｎｋａｉＵｎｉｖｅｒｓｉｔｙ）Ａｂｓｔ...     

Error bounds for anisotropic RBF interpolation

We present error bounds for the interpolation with anisotropically transformed radial basis functions for both a function and its partial derivatives. The bounds rely on a growth function and do not contain unknown constants. For polyharmonic basic functions in ${\mathbb{R}}^2$, we show that the anisotropic estimates predict a significant improvement of the approximation error if both the target function and the placement of the centers are anisotropic, and this improvement is confirmed numerically.