高波数Helmholtz方程的超收敛分析

本文考虑求解Helmholtz方程的有限元方法的超逼近性质以及基于PPR后处理方法的超收敛性质.我们首先给出了矩形网格上的p-次元在收敛条件k(kh)~(2p+1)≤C_0下的有限元解和基于Lobatto点的有限元插值之间的超逼近以及重构的有限元梯度和精确解之间的超收敛分析.然后我们给出了四边形网格上的线性有限元方法的分析.这些估计都给出了与波数k和网格尺寸h的依赖关系.同时我们回顾了三角形网格上的线性有限元的超收敛结果.最后我们给出了数值实验并且结合Richardson外推进一步减少了误差.     

高波数问题的超收敛性

 本文考虑求解Helmholtz方程的有限元方法的超逼近性质以及基于PPR后处理方法的超收敛性质.我们首先给出了矩形网格上的p-次元在收敛条件k（kh）^2p+1≤C₀下的有限元解和基于Lobatto点的有限元插值之间的超逼近以及重构的有限元梯度和精确解之间的超收敛分析.然后我们给出了四边形网格上的线性有限元方法的分析.这些估计都给出了与波数k和网格尺寸h的依赖关系.同时我们回顾了三角形网格上的线性有限元的超收敛结果.最后我们给出了数值实验并且结合Richardson外推进一步减少了误差.     

曲边区域上的多边形网格间断有限元离散及其多重网格算法

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.     

四阶障碍问题的稳定化混合有限元方法

本文讨论了四阶障碍问题的稳定化混合有限元方法.首先,引入网格依赖范数,通过加罚方法得到了与四阶障碍问题的等价的混合变分形式.随后给出了基于C~0协调有限元空间(W_h,V_h)的混合有限元逼近,例如P_k-P_k三角形有限元.在网格依赖范数下,(W_h,V_h)满足离散的inf-sup条件.最后,我们在不同的假设下,得到了一些误差估计.     

拟等级网格下非线性延迟微分方程间断有限元法

本文利用间断有限元法求解非线性延迟微分方程,在拟等级网格下.给出非线性延迟微分方程间断有限元解的整体收敛阶和局部超收敛阶,数值实验验证了理论结果的正确性.     

反应扩散方程的连续时空有限元方法

本文研究了反应扩散方程的连续时空有限元方法.首先建立了其连续时空有限元格式并证明了有限元解的存在唯一性及稳定性.然后通过引入时空投影算子在没有时空网格限制的条件下给出其近似解在节点处的L~2,H~1最优范数估计以及全局L~2(L~2),L~2(H~1)最优范数估计.最后给出两个数值算例来验证方法的有效性与灵活性并说明结论的正确性.     

各向异性EQrot1非协调元高精度分析的一般格式

本文在各向异性网格下讨论了一般二阶椭圆方程的EQrot1非协调有限元逼近.利用Taylor展开,积分恒等式和平均值技巧导出了一些关于该元新的高精度估计.再结合该元所具有的二个特殊性质:(a)当精确解属于H3(Ω)时,其相容误差为O(h2)阶比它的插值误差高一阶;(b)插值算子与Ritz投影算子等价,得到了在能量模意义下O(h2)阶的超逼近性质.进而,借助于插值后处理技术给出了整体超收敛的一般估计式.     

可混溶驱动问题的超收敛性

本文讨论多孔介质中两相可混溶渗流驱动问题的有限元方法，采用一致网格剖分、指标为k的Raviart－Thomas空间对压力作混合有限元逼近，用正则剖分、逼近阶为l的标准有限元方法处理浓度方程，通过核函数对有限元解作卷积进行局部平均确定非线性项的系数，得到了浓度误差H1范数的超收敛估计，经高阶插值，得到了整体高精度的逼近．     

变度量各向异性网格生成算法和各向异性椭圆偏微分方程匹配网格

对于各向异性问题而言,在进行数值求解时,一个相匹配的各向异性网格至关重要.本文针对二阶椭圆偏微分方程,当其系数矩阵为各向异性时,给出了系数矩阵的逆作为相匹配的各向异性度量,在此度量下生成的各向异性网格即为匹配网格.本文还提出了新的变度量各向异性网格生成算法,该算法通过在各向异性背景网格上结合各向异性Delaunay准则和基于力平衡的结点优化函数生成高质量的各向异性三角化网格.通过数值算例可知,在匹配的各向异性三角化网格上,不仅有条件数小的代数离散系统,还有高精度的数值解,更甚者,在匹配网格结点上的l2范数误差有超收敛现象.     

变系数四维椭圆方程张量积矩形有限元的高精度分析

<正>1引言有限元超收敛性对于有限元实际计算的精度提高、后验误差估计和自适应计算等都有着很重要的作用,因而一直备受人们的关注.目前,关于一维、二维和三维问题的有限元超收敛研究已有许多文献[1-12].由于区域复杂性等原因,三维有限元超收敛研究相比一、二维困难很多,且目前所获结果大都是L2范数意义下的平均超收敛结果.本文将研     

Computer-based proof of the existence of superconvergence points in the finite element method; superconvergence of the derivatives in finite element solutions of Laplace's,Poisson's,and the elasticity equations

In this article we address the problem of the existence of superconvergence points for finite element solutions of systems of linear elliptic equations. Our approach is quite different from all other studies of superconvergence. We prove that the existence of superconvergence points can be guaranteed by a numerical algorithm, which employs a finite number of operations (provided that there is no roundoff-error). By employing this approach, we can reproduce all known results on superconvergence of finite element solutions for linear elliptic problems and we can obtain many new results. Here, in particular, we address the problem of the superconvergence points for the gradient of finite element solutions of Laplace's and Poisson's equations and we show that the sets of superconvergence points are very different for these two cases. We also study the superconvergence of the components of the gradient of the displacement, the strain and stress for finite element solutions of the equations of elasticity. For Laplace's and Poisson's equations (resp. the equations of elasticity), we consider meshes of triangular as well as square elements of degree p, 1 ? p ? 7 (resp. 1 ? p ? 4). For the meshes of triangular elements we investigate the effect of the geometry of the mesh by considering four mesh patterns that typically occur in practical meshes, while in the case of square elements, we study the effect of the element type (tensor-product, serendipity, or other). © 1996 John Wiley & Sons, Inc.     

高次Galerkin有限元法的超收敛性

通过局部误差估计,对具有光滑边界的二阶常系数椭圆型方程,给出了高次Ｇａｌｅｒｋｉｎ 有限元法的超收敛性． 运用对称技巧和积分恒等式技巧,在局部对称矩形网格或三角形网格上,我们得到了改进的超收敛性（提高１－ ３ 阶）．     

Convergence and superconvergence analysis of a new quadratic Hermite-type triangular element on anisotropic meshes

A new quadratic Hermite-type triangular finite element is conceived to solve a class of two-dimensional second-order elliptic boundary value problems. Its error estimates on anisotropic meshes are developed. Furthermore, we verify that some conditions set to the meshes contribute to the proof of its superconvergence properties, which can improve the approximation results. Numerical examples are given to confirm our theoretical analysis.     

Numerical study of natural superconvergence in least-squares finite element methods for elliptic problems

Natural superconvergence of the least-squares finite element method is surveyed for the one-and two-dimensional Poisson equation. For two-dimensional problems, both the families of Lagrange elements and Raviart-Thomas elements have been considered on uniform triangular and rectangular meshes. Numerical experiments reveal that many superconvergence properties of the standard Galerkin method are preserved by the least-squares finite element method. The second author was supported in part by the US National Science Foundation under Grant DMS-0612908.     

Superconvergence analysis of FEMs for the Stokes–Darcy system

We consider a superconvergence analysis for quadratic finite element approximations of the Stokes–Darcy system. The superclose property of an extra half order is proven for uniform triangular meshes. Based on the result of the superclose property, global superconvergence is derived by applying a postprocessing technique. In addition, some numerical examples are presented to demonstrate our theoretical results. Copyright © 2010 John Wiley & Sons, Ltd.     

Recovery type a posteriori estimates and superconvergence for nonconforming FEM of eigenvalue problems

The main goal of this paper is to present recovery type a posteriori error estimators and superconvergence for the nonconforming finite element eigenvalue approximation of self-adjoint elliptic equations by projection methods. Based on the superconvergence results of nonconforming finite element for the eigenfunction we derive superconvergence and recovery type a posteriori error estimates of the eigenvalue. The results are based on some regularity assumption for the elliptic problem and are applicable to the lowest order nonconforming finite element approximations of self-adjoint elliptic eigenvalue problems with quasi-regular partitions. Therefore, the results of this paper can be employed to provide useful a posteriori error estimators in practical computing under unstructured meshes.     

多项时间分数阶扩散方程各向异性线性三角元的高精度分析

在各向异性网格下,针对具有Caputo导数的二维多项时间分数阶扩散方程,给出了线性三角形元的高精度分析.首先,基于线性三角形元和改进的L1格式,建立了一个全离散逼近格式,并证明了其无条件稳定性;其次,利用有限元插值算子与Riesz投影算子之间的关系及相关的高精度结果,导出了超逼近性质.进而,借助于插值后处理技术得到了超收敛估计.值得指出的是,单独利用插值算子或Riesz投影都无法得到上述超逼近和超收敛结果.最后,利用数值算例验证了理论分析的正确性.此外,对一些常见的有限单元在该方程的数值逼近方面,作了进一步探讨.     

非线性sine-Gordon方程的各向异性线性元高精度分析新模式

在各向异性网格下,针对一类非线性sine-Gordon方程提出了线性三角形元新的高精度分析模式.基于该元的积分恒等式结果,导出了插值与Riesz投影之间的误差估计,再借助于插值后处理技术得到了在半离散和全离散格式下单独利用插值或Riesz投影所无法得到的超逼近和超收敛结果.最后,对一些常见的单元作了进一步探讨.     

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.