高次有限元外推的进一步研究

对高次矩形元,我们给出了一个外推算法.利用离散格林函数权模估计和渐进不等式展开,证明了外推算法可以提高泊松方程有限元逼近解的精度.     

二次有限元解的渐近展开与外推

在Poisson方程的求解域Ω存在一致的三角剖分,并且相邻两初始单元构成平行四边形的假设下,证明了若Poisson方程的解u属于H6(Ω),那么二次有限元的误差有h4的渐近展开.基于误差的渐近展开,可以利用h4-Richardson外推进一步提高数值解的精度阶,并且能够得到一个后验误差估计.最后,一个数值算例验证了理论分析.     

数值积分对特征值有限元外推的影响

§1.引言在分片一致三角形剖分下,用线性有限元法解特征值问题求得近似特征值λ、λ~(h/2).[1]证明了对λ~h.λ~(h/2)作外推可提高收敛阶:     

板振动问题Ciarlet-RaViart混合有限元渐近展开式与Richardson外推

1引言 Richardson外推是数值计算各个方面广泛使用的一种高效算法,它的理论基础是渐近展开式.1983年文[6]对三角形网格首先证明了二阶椭圆问题协调有限元渐近展开式,这以后有限元外推迅速发展成一个国际性研究课题,文[5],[7]对有限元外推作了系统、深入的论述.板弯曲和板振动是数学力学界关注的课题.对简支板弯曲问题,文[11]证明了Ciarlet-Raviart(简记为C-R)混合有限元渐近展开式.     

关于本征值问题有限元外推的一个注记

§1.引言考虑二阶椭圆微分算子本征值问题: 求,使其中G是R~N(N=1,2)中边界分段光滑的有界域,a(·,·)是定义在H_0~1(G)上的对称、连续,H_0~1(G)椭圆双线性形式,(·,·)是L_2(G)中内积,||·||_0是L_2(G)中范数。假设满足正则性条件:对fL_2(G),“a(w,v)=(f,v),vH_0~1(G)”的解     

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

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

n维矩形域上椭圆问题有限元单方向外推

1 引言 Richardson外推应用于椭圆偏微方程边值问题有限元法始于1978年(见[1],并于1983年在理论研究方面取得突破性进展(见[2]).自那以后有限元外推得到迅速发展,成为一个富于竞争的国际性研究课题(见[3],[4],[5]及其所列参考文献).但是通常的有限元外推需同时在每一个方向上分半加密网格,因此,对n维问题,细网格的结点数是粗网格的2~n倍,结果当n较大时(高维问题),细网格上的计算工作量十分庞大.为了克服这个缺点,发展了有限元单方向外推.对Poisson方程边值问题,[6]研究了2维矩形域上双线性有限元解的单方向外推,[7]研究了3维矩形域上三线性有限元单方向外推必须的插值渐近展开式,[8]研究了n维矩形域上n线性有限元解的区域分裂外推.本文旨在研究n维矩形域上Poisson方程边值问题及其对应的本征值问题n线性有限元解的单方向外推.始终假设本文出现的函数u是连续的.     

粘弹性方程ACM有限元的超收敛分析和外推

本文通过积分恒等式技巧和插值后处理技术,得到了粘弹性方程ACM有限元的超逼近和超收敛性质.进一步利用Bramble-Hilbert引理构造了一个合适的外推格式,得到了比以往文献高一阶的外推结果.     

四阶方程两点边值问题Hermite有限元解的渐近展式与外推

1引言有限元解的渐近展式是提高微分方程数值解精度的重要工具,比如亏量校正和外推就是建立在有限元解的渐近展式的基础之上．许多作者对此进行了大量的研究(见[1]-[4]),特别是文[1],提出了在研究有限元解的渐近展式中十分有用的能量嵌入技巧．本文利用能量嵌入定理得到了四阶方程两点边值问题Hermite有限元解及其二阶平均导数的渐近展式,进一步我们还讨论了它们的Richardson外推公式．考虑四阶方程两点边值问题     

自适应有限元和后验误差估计——渐近准确估计

在[7]中,作者讨论了有限元误差的1-模等价估计.本文是[7]的继续,给出一种自适应有限元计算中误差的1-模渐近准确估计,即对于误差的1-模||e||_1,Ω给出可计算的估计量?,当||e||_1,Ω→0时,成立?/||e||_1,Ω→1. 本文将沿用[7]中的定义及符号.     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.     

有限元插值误差的下界估计

基于一维区域上的拟一致剖分,证明了线性元插值误差的最优下界估计.基于此并利用超收敛理论,我们得到了有限元离散误差的上、下界.     

Finite Field Elements with Specified Order and Traces

Results of F. Özbudak on the existence of elements in a finite extension of a finite field with specified order and at which given rational functions have specified traces over the ground field are obtained by the method of character sums. There are significant applications to the construction of ordered arrays.AMS classification: 11T30, 11T24     

Analysis of Linear Triangular Elements for Convection-diffusion Problems by Streamline Diffusion Finite Element Methods

This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems.In [8],under the condition thatε≤h~2 the optimal finite element error estimate was obtained in L~2-norm.In the present paper,however,the same error estimate result is gained under the weaker condition thatε≤h.     

A Stabilized Low Order Finite Element Method for Three Dimensional Elasticity Problems

We introduce a low order finite element method for three dimensional elasticity problems. We extend Kouhia-Stenberg element [12] by using two nonconforming components and one conforming component, adding stabilizing terms to the associated bilinear form to ensure the discrete Korn's inequality. Using the second Strang's lemma, we show that our scheme has optimal convergence rates in $L^2$ and piecewise $H^1$-norms even when Poisson ratio $\nu$ approaches $1/2$. Even though some efforts have been made to design a low order method for three dimensional problems in [11,16], their method uses some higher degree basis functions. Our scheme is the first true low order method. We provide three numerical examples which support our analysis. We compute two examples having analytic solutions. We observe the optimal $L^2$ and $H^1$ errors for many different choices of Poisson ratios including the nearly incompressible cases. In the last example, we simulate the driven cavity problem. Our scheme shows non-locking phenomena for the driven cavity problems also.     

Development and Analysis of Higher Order Finite Volume Methods over Rectangles for Elliptic Equations

Currently used finite volume methods are essentially low order methods. In this paper, we present a systematic way to derive higher order finite volume schemes from higher order mixed finite element methods. Mostly for convenience but sometimes from necessity, our procedure starts from the hybridization of the mixed method. It then approximates the inner product of vector functions by an appropriate, critical quadrature rule; this allows the elimination of the flux and Lagrange multiplier parameters so as to obtain equations in the scalar variable, which will define the finite volume method. Following this derivation with different mixed finite element spaces leads to a variety of finite volume schemes. In particular, we restrict ourselves to finite volume methods posed over rectangular partitions and begin by studying an efficient second-order finite volume method based on the Brezzi–Douglas–Fortin–Marini space of index two. Then, we present a general global analysis of the difference between the solution of the underlying mixed finite element method and its related finite volume method. Then, we derive finite volume methods of all orders from the Raviart–Thomas two-dimensional rectangular elements; we also find finite volume methods to associate with BDFM ₂ three-dimensional rectangles. In each case, we obtain optimal error estimates for both the scalar variable and the recovered flux.     

Stokes问题低阶非协调混合元的改进加罚算法(英文)

通过对由经典加罚算法得到的两个解进行线性组合,研究Stokes方程低阶非协调混合元的改进加罚算法.该方法利用较大的罚参数能得到同使用较小参数的经典加罚方法一样的收敛阶.此外,基于单元的特性和插值后处理技巧,得到一些超收敛结果,从而改进以往的文献结果.