导数小片插值恢复技术与超收敛性

１．引言 有限元超收敛的研究自七十年代起至今方兴未艾．现有的研究工作基本遵循两种途径：一是找出有限元插值逼近的超收敛点,然后再利用插值弱估计等手段导出有限元解本身所具有的超收敛性质［１,２］;二是利用各种后处理技术,如平均技术,投影技术,外插技术和插值有限元技术等[3,6],来导出经过后处理的有限元解的超收敛性．近年来一种新的超收敛后处理技术,即所谓的“Ｚ-Ｚ导数小片恢复技术”,得到众多的研究［７-１１］,并被 Ｂａｂｕｓｋａ等人认为是用于渐进准确的后验误差估计效果最好的技术之一［１２］．这种技术是利用…     

Ultraconvergence of the patch recovery technique II

The ultraconvergence property of a gradient recovery technique proposed by Zienkiewicz and Zhu is analyzed for the Laplace equation in the two dimensional setting. Under the assumption that the pollution effect is not present or is properly controlled, it is shown that the convergence rate of the recovered gradient at an interior node is two orders higher than the optimal global convergence rate when even-order finite element spaces and local uniform rectangular meshes are used.

     

Ultraconvergence of the patch recovery technique

The ultraconvergence property of a derivative recovery technique recently proposed by Zienkiewicz and Zhu is analyzed for two-point boundary value problems. Under certain regularity assumptions on the exact solution, it is shown that the convergence rate of the recovered derivative at an internal nodal point is two orders higher than the optimal global convergence rate when even-order finite element spaces and local uniform meshes are used.

     

Stokes问题的一个非协调有限元

<正>1引言定常Stokes问题是流体力学中的一种重要问题,是标准的混合问题,速度与压力同时计算.关于该问题有限元求解的文章很多([1],[2],[4],[5],[6],[8],[9]),分析的难点在于单元必须满足离散的Babuska-Brezzi([2])条件.在([4])中提出了著名的非协调Crouzeix-Raviart     

鱼骨形网格上二阶方程混合元的超收敛

In this paper, superconvergence of the lowest order Raviart-Thomas mixed finite element approximation for second order Neumann boundary value problem on fishbone shape meshes is analyzed. The main term of the error between the exact solution and the finite element interpolating function is determined by Bramble-Hilbert lemma on the individual finite element. A part of the main term of the error on two adjacent finite elements can be cancelled along the special direction, and thus the higher order error estimate is obtained on the whole domain by summation. Compared with the general finite element error estimate,the convergence rate can be increased from order one to order two in L2-norm by postprocessing superconvergence technique.     

Stokes问题Q_2-P_1混合元外推方法

考虑Stokes问题的有限元解与精确解插值的Q2-P1混合元的渐近误差展开和分裂外推.首先利用积分恒等式技巧确定了微分方程精确解与有限元插值之间积分式的主项,其次再借助插值后处理和分裂外推技术,得到了比通常的误差估计提高两阶的收敛速度.     

Interior estimates for a low order finite element method for the Reissner–Mindlin plate model

Interior error estimates are obtained for a low order finite element introduced by Arnold and Falk for the Reissner–Mindlin plates. It is proved that the approximation error of the finite element solution in the interior domain is bounded above by two parts: one measures the local approximability of the exact solution by the finite element space and the other the global approximability of the finite element method. As an application, we show that for the soft simply supported plate, the Arnold–Falk element still achieves an almost optimal convergence rate in the energy norm away from the boundary layer, even though optimal order convergence cannot hold globally due to the boundary layer. Numerical results are given which support our conclusion. This revised version was published online in June 2006 with corrections to the Cover Date.     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

二阶椭圆问题的一类广义有限元法

本文提出了求解二阶椭圆问题的一类广义有限元方法,分析了广义有限元方法的优越性,证明了二阶椭圆问题的广义有限元方法具有比标准的Galerkin有限元方法更高阶的收敛速度,根据插值算子的性质,进一步证明了有限元解的亏量迭代校正收敛到广义有限元解,并用数值例子说明广义有限元方法是有效的.     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

Error analysis of finite element approximations of the stochastic Stokes equations

Numerical solutions of the stochastic Stokes equations driven by white noise perturbed forcing terms using finite element methods are considered. The discretization of the white noise and finite element approximation algorithms are studied. The rate of convergence of the finite element approximations is proved to be almost first order (h|ln h|) in two dimensions and one half order ( h^\frac12h^{\frac{1}{2}}) in three dimensions. Numerical results using the algorithms developed are also presented.     

Postprocessing and higher order convergence for the mixed finite element approximations of the eigenvalue problem

In this paper, we propose a method to improve the convergence rate of the lowest order Raviart-Thomas mixed finite element approximations for the second order elliptic eigenvalue problem. Here, we prove a supercloseness result for the eigenfunction approximations and use a type of finite element postprocessing operator to construct an auxiliary source problem. Then solving the auxiliary additional source problem on an augmented mixed finite element space constructed by refining the mesh or by using the same mesh but increasing the order of corresponding mixed finite element space, we can increase the convergence order of the eigenpair approximation. This postprocessing method costs less computation than solving the eigenvalue problem on the finer mesh directly. Some numerical results are used to confirm the theoretical analysis.     

Raviart-Thomas混合元的超收敛

考虑二阶椭圆方程Dirichlet边值问题在正则矩形网格上k阶RaviartThomas混合有限元的超收敛.对有限元解经插值处理后,与通常的有限元最优误差估计相比,收敛速度提高了两阶.     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

半线性Sobolev方程的H~1-Galerkin混合有限元方法

利用H~1-Galerkin混合有限元方法研究了一维半线性Sobolev方程,得到了半离散解的最优阶误差估计,优点是不需验证LBB相容性条件.     

SUPERCONVERGENCE ANALYSIS OF THE POLYNOMIAL PRESERVING RECOVERY FOR ELLIPTIC PROBLEMS WITH ROBIN BOUNDARY CONDITIONS

We analyze the superconvergence property of the linear finite element method based on the polynomial preserving recovery(PPR)for Robin boundary elliptic problems on triangulartions.First,we improve the convergence rate between the finite element solution and the linear interpolation under the H1-norm by introducing a class of meshes satisfying the Condition(α,σ,μ).Then we prove the superconvergence of the recovered gradients post-processed by PPR and define an asymptotically exact a posteriori error estimator.Finally,numerical tests are provided to verify the theoretical findings.     

Error estimates of H^1-Galerkin mixed finite element method for Schrodinger equation

An H1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.     

DISCONTINUOUS FINITE ELEMENT METHOD FOR CONVECTION-DIFFUSION EQUATIONS

1. IlltroductionThe finite element approximation of the convection--diffusin equations has been investigated using several different approaches (see e.g. [3] [4] and the references therein).Previous analysis in primal formulation of these problems was done for two types ofapproximation schemes: one which produces a continuous piecewise polynomial approximation and one which produces a piecewise polynomial approximation which arecontinuous for certain number of moments accross interelement edge…     

Parabolic finite volume element equations in nonconvex polygonal domains

We study spatially semidiscrete and fully discrete finite volume element approximations of the heat equation with homogeneous Dirichlet boundary conditions in a plane polygonal domain with one reentrant corner. We show that, as a result of the singularity in the solution near the reentrant corner, the convergence rate is reduced from optimal second order, similarly to what was shown for the finite element method in the earlier work 2 . Optimal order convergence may be restored by mesh refinement near the corners of the domain. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

加罚Navier—Stokes方程的最佳非线性Galerkin算法

该文提出了求解二维加罚Navier-Stokes方程的最佳非线性Galerkin算法．这个算法在于在粗网格有限元空间上求解一非线性子问题，在细网格增量有限元空间Wh上求解一线性子问题．如果线性有限元被使用及，则该算法具有和有限元Galerkin算法同阶的收敛速度．然而该文提出的算法可以节省可观的计算时间．