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

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

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

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

Stokes问题有限元逼近中求导数的提取法

1.引言 我们知道Poisson方程和平面弹性问题的解的导数的近似值可以通过所谓提取公式得到,而不必对近似解直接求导数.这样我们可以得到具有与近似解本身同阶精度的导数的近似值.这一方法已被用于基于插值误差的后验误差估计及相应的自适应有限元方法中本文将这一方法应用于Stokes问题的有限元逼近,从Stokes方程的解的     

自适应有限元和后验误差估计——等价估计

近年来,自适应有限元计算得到了广泛的研究.后验误差估计是实现自适应计算的基础.70年代末,I. Babuska和 W.Rheinboldt给出了建立等价后验估计的一般原则和方法,在一维情形给出了完整的结果.80年代初,I.Babuska和A.Miller对平面弹性问题正方形双线性元,给出了一种后验估计的方法,并证明了其等价性和渐近准确     

基于自适应遗传算法的逐次超松驰迭代法

确定逐次超松驰迭代法中的最佳松驰因子,迄今,人们还没有给出一可行实用的方法.利用自适应遗传算法全局搜索性能、并行性及其遗传操作,构造出近似确定最佳松驰因子的一种自适应进化方法,并由此得到一近似确定ω功能的自适应逐次超松驰迭代算法.数值算例表明,该算法在求解线性方程组中是可行的,实用和快捷的.     

基于拟Shannon小波浅水长波近似方程组的数值解

本文研究了浅水长波近似方程组初边值问题的数值解.利用小波多尺度分析和区间拟Shannon小波,对浅水长波近似方程组空间导数实施空间离散,用时间步长自适应精细积分法对其变换所的非线性常微分方程组进行求解,得到了浅水长波近似方程组的数值解,并将此方法计算的结果与其解析解进行比较和验证.     

Wilson元插值误差渐近估计

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

杆件轴向受迫振动的Galerkin有限元EEP法自适应求解

基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.     

泊松方程及平面弹性问题有限元方法中求高阶导数的提取法

在许多有限元计算中经常在求得近似解后还要求得到近似的解的导数.如在弹性计算中,如何从计算得到的位移近似解较好地计算应力早已被研究多年.如果计算中包含直接对近似解求导数,必然会丧失部分精度,得不到满意的结果.特别,若近似解为分片常数函数,则根本无法从直接求导数得到应力的近似值.Babuska和 Miller提出了所谓“提取法”,即利用推导出来的提取公式来求解的导数的近似值,以得到与近似解本身同     

泊松方程及平面弹性问题有限元方法中求高阶导数的提取法

在许多有限元计算中经常在求得近似解后还要求得到近似的解的导数.如在弹性计算中,如何从计算得到的位移近似解较好地计算应力早已被研究多年.如果计算中包含直接对近似解求导数,必然会丧失部分精度,得不到满意的结果.特别,若近似解为分片常数函数,则根本无法从直接求导数得到应力的近似值.Babuska和 Miller提出了所谓“提取法”,即利用推导出来的提取公式来求解的导数的近似值,以得到与近似解本身同     

非线性积分方程迭代配置法的渐近展开及其外推

非线性积分方程迭代配置法的渐近展开及其外推韩国强（华南理工大学计算机工程与科学系）ＡＳＹＭＰＴＯＴＩＣＥＲＲＯＲＥＸＭＮＳＩＯＮＳＡＮＤＥＸＴＲＡＰＯＬＡＴＩＯＮＦＯＲＴＨＥＩＴＥＲＡＴＥＤＣＯＬＬＯＣＡＴＩＯＮＭＥＴＨＯＤＳＯＦＮＯＮＬＩＮＥＡＲＩ...     

Asymptotically exact a posteriori error estimator for biquadratic elements

This paper addresses the finite element method with a posteriori error estimation for elements of degree p = 1 and p = 2. It gives the formulae for the error indicators and error estimators. Basic mathematical characterizations of the estimators are given and it is shown that the estimators for p = 1 and p = 2 have different structures. Numerical examples show the effectivity of the approach and the high quality of the estimator.     

An Analysis of a Finite Element Method of Low Degree for the Navier-Stokes Problems

There are many papers in which approximate solution of Navier-Stokes problem is discussed by finite element method. Their error estimates are optimal, but degree of piecewise polynomials for pressure p or degree of piecewise polynomials for velocity u are not the lowest. In this papre a new element is given. Its degre for p and degree for u are the lowest and error estimates are optimal.     

用样条逼近Sobolev函数类W_p~1(R)的逼近误差

作者研究了定义在全实轴上的Sobolev函数类W_p~1(R)的逼近问题.以一次样条函数作为逼近工具,给出了p=1和p=∞时的逼近误差.     

Error estimates for 3- narrow finite elements

We obtain error estimates for finite element approximations of the lowest degree valid uniformly for a class of three-dimensional narrow elements. First, for the Lagrange interpolation we prove optimal error estimates, both in order and regularity, in for . For it is known that this result is not true. Applying extrapolation results we obtain an optimal order error estimate for functions sligthly more regular than . These results are valid both for tetrahedral and rectangular elements. Second, for the case of rectangular elements, we obtain optimal, in order and regularity, error estimates for an average interpolation valid for functions in with and .

     

On irreducible elements of semilattices

The concept of (join)-irreducible elements works well, especially for distributive lattices. Therefore our definition of elements of a given degree of irreducibility employs the notion of distributivity as much as possible, even if the irreducibility is defined for elements of a (meet)-semilattice. Via the lattice of hereditary subsets of the poset ofk-irreducible elements of a semilattice (wherek is a cardinal) we obtain a new construction of a D1_k-reflection (a sort of free distributive extension) of the semilattice, provided that there are sufficiently manyk-irreducible elements. The last property is satisfied, for example, if the original semilattice is the dual of an algebraic lattice [Dilworth and Crawley, 1960], but this condition is too restrictive for semilattices. It turns out that, under certain limitations, the D1_k-reflection of a semilattice both preserves and reflects the degree of irreducibility.Presented by R. Freese.     

-curved finite elements and applications to plate and shell problems

Many thin-plate and thin-shell problems are set on plane reference domains with a curved boundary. Their approximation by conforming finite-elements methods requires ¹-curved finite elements entirely compatible with the associated ¹-rectilinear finite elements. In this contribution we introduce a ¹-curved finite element compatible with the P₅-Argyris element, we study its approximation properties, and then, we use such an element to approximate the solution of thin-plate or thin-shell problems set on a plane-curved boundary domain. We prove the convergence and we get a priori asymptotic error estimates which show the very high degree of accuracy of the method. Moreover we obtain criteria to observe when choosing the numerical integration schemes in order to preserve the order of the error estimates obtained for exact integration.     

Residual-free bubbles for advection-diffusion problems: the general error analysis

Summary. We develop the general a priori error analysis of residual-free bubble finite element approximations to non-self-adjoint elliptic problems of the form subject to homogeneous Dirichlet boundary condition, where A is a symmetric second-order elliptic operator, C is a skew-symmetric first-order differential operator, and is a positive parameter. Optimal-order error bounds are derived in various norms, using piecewise polynomial finite elements of degree . Received October 1, 1998/ Revised version received April 6, 1999 / Published online January 27, 2000     

The convergence of spline collocation for strongly elliptic equations on curves

Summary Most boundary element methods for two-dimensional boundary value problems are based on point collocation on the boundary and the use of splines as trial functions. Here we present a unified asymptotic error analysis for even as well as for odd degree splines subordinate to uniform or smoothly graded meshes and prove asymptotic convergence of optimal order. The equations are collocated at the breakpoints for odd degree and the internodal midpoints for even degree splines. The crucial assumption for the generalized boundary integral and integro-differential operators is strong ellipticity. Our analysis is based on simple Fourier expansions. In particular, we extend results by J. Saranen and W.L. Wendland from constant to variable coefficient equations. Our results include the first convergence proof of midpoint collocation with piecewise constant functions, i.e., the panel method for solving systems of Cauchy singular integral equations.Dedicated to Prof. Dr. Dr. h.c. mult. Lothar Collatz on the occasion of his 75th birthdayThis work was begun at the Technische Hochschule Darmstadt where Professor Arnold was supported by a North Atlantic Treaty Organization Postdoctoral Fellowship. The work of Professor Arnold is supported by NSF grant BMS-8313247. The work of Professor Wendland was supported by the Stiftung Volkswagenwerk     

EFK方程一个新的低阶非协调混合有限元方法的高精度分析

对Extended Fisher-Kolmogorov(EFK)方程,利用EQ_1~(rot)元和零阶RaviartThomas(R-T)元建立了一个新的非协调混合元逼近格式.首先,证明了半离散格式逼近解的一个先验估计并证明了逼近解的存在唯一性.在半离散格式下,利用上述两种元的高精度分析结果以及这个先验估计,在不需要有限元解u_h属于L~∞的条件下,得到了原始变量u和中间变量v=-?u的H~1-模以及流量p=u的(L~2)~2-模意义下O(h~2)阶的超逼近性质.同时,借助插值后处理技术,证明了上述变量的具有O(h~2)阶的整体超收敛结果.其次,建立了一个新的线性化向后Euler全离散格式并证明了其逼近解的存在唯一性.另一方面,通过对相容误差和非线性项采取与传统误差分析不同的新的分裂技巧,分别导出了以往文献中尚未涉及的关于u和v在H~1-模以及p在(L~2)~2-模意义下具有O(h~2+τ)阶的超逼近性质,进一步地,借助插值后处理技术,得到了上述变量的整体超收敛结果.这里h和τ分别表示空间剖分参数和时间步长.最后,给出了一个数值算例,计算结果验证了理论分析的正确性.