一阶双曲问题间断有限元的后验误差分析

一阶双曲问题的有限元后验误差估计至今没有得到很好的解决.本文对d维区域上一阶双曲问题的k次间断有限元逼近提出了一种新的后验误差分析方法, 进而建立了间断有限元解在DG范数下(强于L₂范数)基于误差余量型的后验误差估计. 数值计算验证了本文理论分析的有效性. 本文方法也适用于其他变分问题有限元逼近的后验误差分析.     

A UNIFIED A POSTERIORI ERROR ANALYSIS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF REACTIVE TRANSPORT EQUATIONS

Four primal discontinuous Galerkin methods are applied to solve reactive transportproblems, namely, Oden-Babuska-Baumann DG (OBB-DG), non-symmetric interior penaltyGalerkin (NIPG), symmetric interior penalty Galerkin (SIPG), and incomplete interiorpenalty Galerkin (IIPG). A unified a posteriori residual-type error estimation is derivedexplicitly for these methods. From the computed solution and given data, explicit esti-mators can be computed efficiently and directly, which can be used as error indicators foradaptation. Unlike in the reference [10], we obtain the error estimators in L~2 (L~2) norm byusing duality techniques instead of in L~2 (H~1) norm.     

积分算子本征值问题伽略金法的后验误差指示子

1引言 设G是R~n中有界域，积分算予Tx(s)=integral from n k(s.t)x(t)dt.(s∈G)是映L~2(G)到L~2(G)中的自共轭全连续算子。△={△}是G的拟一致部分。     

A POSTERIORI ENERGY-NORM ERROR ESTIMATES FOR ADVECTION-DIFFUSION EQUATIONS APPROXIMATED BY WEIGHTED INTERIOR PENALTY METHODS

We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms： a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.     

A POSTERIORI ERROR ESTIMATES FOR A MODIFIED WEAK GALERKIN FINITE ELEMENT APPROXIMATION OF SECOND ORDER ELLIPTIC PROBLEMS WITH DG NORM

In this paper,we derive a residual based a posteriori error estimator for a modified weak Galerkin formulation of second order elliptic problems.We prove that the error estimator used for interior penalty discontinuous Galerkin methods still gives both upper and lower bounds for the modified weak Galerkin method,though they have essentially different bilinear forms.More precisely,we prove its reliability and efficiency for the actual error measured in the standard DG norm.We further provide an improved a priori error estimate under minimal regularity assumptions on the exact solution.Numerical results are presented to verify the theoretical analysis.     

UNIFORM STABILITY AND ERROR ANALYSIS FOR SOME DISCONTINUOUS GALERKIN METHODS

In this paper,we provide a number of new estimates on the stability and convergence of both hybrid discontinuous Galerkin(HDG)and weak Galerkin(WG)methods.By using the standard Brezzi theory on mixed methods,we carefully define appropriate norms for the various discretization variables and then establish that the stability and error estimates hold uniformly with respect to stabilization and discretization parameters.As a result,by taking appropriate limit of the stabilization parameters,we show that the HDG method converges to a primal conforming method and the WG method converges to a mixed conforming method.     

SUPERCONVERGENCE AND A POSTERIORI ERROR ESTIMATES FOR BOUNDARY CONTROL GOVERNED BY STOKES EQUATIONS

In this paper, the superconvergence results are derived for a class of boundary con-trol problems governed by Stokes equations. We derive superconvergence results for boththe control and the state approximation. Base on superconvergence results, we obtainasymptotically exact a posteriori error estimates.     

A PRIORI AND A POSTERIORI ERROR ESTIMATES OF A WEAKLY OVER-PENALIZED INTERIOR PENALTY METHOD FOR NON-SELF-ADJOINT AND INDEFINITE PROBLEMS

In this paper, we study a weakly over-penalized interior penalty method for non-self- adjoint and indefinite problems. An optimal a priori error estimate in the energy norm is derived. In addition, we introduce a residual-based a posteriori error estimator, which is proved to be both reliable and efficient in the energy norm. Some numerical testes are presented to validate our theoretical analysis.     

稳态Poisson-Nernst-Planck方程的残量型后验误差估计

针对稳态的Poisson-Nernst-Planck方程研究了一种残量型的后验误差估计子,对方程的两个解-浓度和电势,都分别给出了上界和下界估计.数值实验表明,基于这种后验误差估计子构造的自适应有限元算法对于稳态的Poisson-Nernst-Planck方程是有效的.     

A POSTERIORI ERROR ESTIMATES OF A NON-CONFORMING FINITE ELEMENT METHOD FOR PROBLEMS WITH ARTIFICIAL BOUNDARY CONDITIONS

An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.     

A NEW APPROACH TO RECOVERY OF DISCONTINUOUS GALERKIN

A new recovery operator P ：Qn^disc（T）→Qn＋1^disc（M） for discontinuous Galerkin is derived. It is based on the idea of projecting a discontinuous, piecewise polynomial solution on a given mesh T into a higher order polynomial space on a macro mesh M. In order to do so, we define local degrees of freedom using polynomial moments and provide global degrees of freedom on the macro mesh. We prove consistency with respect to the local L2-projection, stability results in several norms and optimal anisotropic error estimates. As an example, we apply this new recovery technique to a stabilized solution of a singularly perturbed convection-diffusion problem using bilinear elements.     

SUPERCONVERGENCE OF DISCONTINUOUS GALERKIN METHOD FOR NONSTATIONARY HYPERBOLIC EQUATION

AbstractFor the first order nonstationary hyperbolic equation taking the piecewise linear discontinuous Galerkin solver, we prove that under the uniform rectangular partition, such a discontinuous solver, after postprossesing, can have two and half approximative order which is half order higher than the optimal estimate by Lesaint and Raviart under the rectangular partition.     

A NEW DIRECT DISCONTINUOUS GALERKIN METHOD WITH SYMMETRIC STRUCTURE FOR NONLINEAR DIFFUSION EQUATIONS

In this paper we continue the study of discontinuous Galerkin finite element methods for nonlinear diffusion equations following the direct discontinuous Galerkin （DDG） meth- ods for diffusion problems [17] and the direct discontinuous Galerkin （DDG） methods for diffusion with interface corrections [18]. We introduce a numerical flux for the test func- tion, and obtain a new direct discontinuous Galerkin method with symmetric structure. Second order derivative jump terms are included in the numerical flux formula and explicit guidelines for choosing the numerical flux are given. The constructed scheme has a sym- metric property and an optimal L2 （L2） error estimate is obtained. Numerical examples are carried out to demonstrate the optimal （k ＋ 1）th order of accuracy for the method with pk polynomial approximations for both linear and nonlinear problems, under one-dimensional and two-dimensional settings.     

MIXED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR ONE-DIMENSIONAL FULLY NONLINEAR SECOND ORDER ELLIPTIC AND PARABOLIC EQUATIONS

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations （PDEs）. In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin （IP-DG） methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative uxx of the solution u, three independent functions p1,p2 and p3 are introduced to represent numerical derivatives using various one-sided limits. The proposed DG frame- work, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative uxz. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator which is consistent with the differential operator and satisfies certain monotonicity （called g-monotonicity） properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG frame- work. The g-monotonicity gives the DG methods the ability to select the mathematically ＂correct＂ solution （i.e., the viscosity solution） among all possible solutions. Moreover, the g-monotonicity allows for the possible development of more efficient nonlinear solvers as the special nonlinearity of the algebraic systems can be explored to decouple the equations. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.     

ERROR ESTIMATES OF GALERKIN METHOD FOR HIGH DIMENSIONAL STEADY STATE KURAMOTO-SIVASHINSKY EQUATION

1　ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅｐａｐｅｒ,ｗｅｃｏｎｓｉｄｅｒｔｈｅｆｏｌｌｏｗｉｎｇＫ -Ｓｅｑｕａｔｉｏｎ (ｆｏｒｅｘａｍｐｌｅ,ｓｅｅ [1 ] -[3]ａｎｄｒｅｆｅｒｅｎｃｅｓｔｈｅｒｅｉｎ)ｉｎｓｐａｃｅｄｉｍｅｎｓｉｏｎ≥ 3:Δ2 ｕ+Δｕ+ 12 ｜ ｕ｜2 -12 Ｎｉ=1ｌｉ∫ＩＮ｜ ｕ｜2 ｄｘ=0 ,ｘ∈ＩＮ,(1 )ｗｉｔｈＬ ｐｅｒｉｏｄｉｃｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎｓ,Ｌ =(ｌ1,ｌ2 ,… ,ｌＮ)ｂｅｉｎｇｔｈｅｓｉｚｅｏｆａｔｙｐｉｃａｌｐａｔｔｅｒｎｃｅｌｌ (ａｌｓｏｂｉｆｕｒｃａｔｉｏｎｐａｒ…     

SUB-OPTIMAL CONVERGENCE OF DISCONTINUOUS GALERKIN METHODS WITH CENTRAL FLUXES FOR LINEAR HYPERBOLIC EQUATIONS WITH EVEN DEGREE POLYNOMIAL APPROXIMATIONS

In this paper,we theoretically and numerically verify that the discontinuous Galerkin(DG)methods with central fluxes for linear hyperbolic equations on non-uniform meshes have sub-optimal convergence properties when measured in the L2-norm for even degree polynomial approximations.On uniform meshes,the optimal error estimates are provided for arbitrary number of cells in one and multi-dimensions,improving previous results.The theoretical findings are found to be sharp and consistent with numerical results.