Superconvergence of H 1-Galerkin Mixed Finite Element Methods for Second-Order Elliptic Equations

We derive superconvergence result for H ¹-Galerkin mixed finite element method for second-order elliptic equations over rectangular partitions. Compared to standard mixed finite element procedure, the method is not subject to the Ladyzhenskaya–Bab?ska–Brezzi (LBB) condition and the approximating finite element spaces are allowed to be of different polynomial degrees. Superconvergence estimate of order 𝒪(h ^k+3), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas elements, is established for the flux via a postprocessing technique.     

强阻尼波动方程的H1-Galerkin混合有限元超收敛分析

研究了强阻尼波动方程的H¹-Galerkin混合有限元方法的超收敛性. 借助于协调线性三角形元已有的分析估计式, 直接利用插值算子代替原始变量 u 的 Ritz 投影和应力变量 p 的 Ritz-Volterra 投影,对半离散和全离散格式, 得到了u在 H¹(Ω) 模和 p 在 H(div;Ω) 模意义下比以往文献高一阶的超逼近和超收敛结果.     

Further results on ultraconvergence derivative recovery for odd-order rectangular finite elements

For rectangular finite element, we give a superconvergence method by SPR technique based on the generalization of a new ultraconvergence record and the sharp Green function estimates, by which we prove that the derivative has ultra-convergence of order O(h ^k+3) (k ⩾ 3 being odd) and displacement has order of O(h ^k+4) (k ⩾ 4 being even) at the locally symmetry points.      

Nonconforming <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-Galerkin mixed FEM for Sobolev equations on anisotropic meshes

A nonconforming H^1-Calerkin mixed finite element method is analyzed for Sobolev equations on anisotropic meshes. The error estimates are obtained without using Ritz-Volterra projection.     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

Lagrange interpolation and finite element superconvergence

We consider the finite element approximation of the Laplacian operator with the homogeneous Dirichlet boundary condition, and study the corresponding Lagrange interpolation in the context of finite element superconvergence. For d‐dimensional Q_k‐type elements with d ≥ 1 and k ≥ 1, we prove that the interpolation points must be the Lobatto points if the Lagrange interpolation and the finite element solution are superclose in H¹ norm. For d‐dimensional P_k‐type elements, we consider the standard Lagrange interpolation—the Lagrange interpolation with interpolation points being the principle lattice points of simplicial elements. We prove for d ≥ 2 and k ≥ d + 1 that such interpolation and the finite element solution are not superclose in both H¹ and L² norms and that not all such interpolation points are superconvergence points for the finite element approximation. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 33–59, 2004.     

Superconvergence of triangular mixed finite elements for optimal control problems with an integral constraint

In this paper, we shall investigate the superconvergence property of quadratic elliptical optimal control problems by triangular mixed finite element methods. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite elements and the control is discretized by piecewise constant functions. We prove the superconvergence error estimate of h² in L²-norm between the approximated solution and the interpolation of the exact control variable. Moreover, by postprocessing technique, we find that the projection of the discrete adjoint state is superclose (in order h²) to the exact control variable.     

Error estimates of H 1-Galerkin mixed finite element method for Schrödinger equation

An H^1-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.     

Numerical Solutions for second-kind Volterra integral equations by Galerkin methods

In this paper, we study the global convergence for the numerical solutions of nonlinear Volterra integral equations of the second kind by means of Galerkin finite element methods. Global superconvergence properties are discussed by iterated finite element methods and interpolated finite element methods. Local superconvergence and iterative correction schemes are also considered by iterated finite element methods. We improve the corresponding results obtained by collocation methods in the recent papers [6] and [9] by H. Brunner, Q. Lin and N. Yan. Moreover, using an interpolation post-processing technique, we obtain a global superconvergence of the O(h ^2r )-convergence rate in the piecewise-polynomial space of degree not exceeding (r–1). As a by-product of our results, all these higher order numerical methods can also provide an a posteriori error estimator, which gives critical and useful information in the code development.     

A two-grid method with expanded mixed element for nonlinear reaction-diffusion equations

Expanded mixed finite element approximation of nonlinear reaction-diffusion equations is discussed. The equations considered here are used to model the hydrologic and bio-geochemical phenomena. To linearize the mixed-method equations, we use a two-grid method involving a small nonlinear system on a coarse gird of size H and a linear system on a fine grid of size h. Error estimates are derived which demonstrate that the error is O(△t + h k+1 + H 2k+2 d/2 ) (k ≥ 1), where k is the degree of the approximating space for the primary variable and d is the spatial dimension. The above estimates are useful for determining an appropriate H for the coarse grid problems.     

High accuracy analysis of the lowest order <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-Galerkin mixed finite element method for nonlinear sine-Gordon equations

The lowest order H¹-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart-Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h²)/O(h² + τ²) in H¹-norm and H(div;Ω)-norm are deduced for the semi-discrete and the fully-discrete schemes, where h, τ denote the mesh size and the time step, respectively, which improve the results in the previous literature.     

Superconvergence of Solution Derivatives of the Shortley–Weller Difference Approximation to Elliptic Equations with Singularities Involving the Mixed Type of Boundary Conditions

This paper presents a superconvergence analysis for the Shortley–Weller finite difference approximation of second-order self-adjoint elliptic equations with unbounded derivatives on a polygonal domain with the mixed type of boundary conditions. In this analysis, we first formulate the method as a special finite element/volume method. We then analyze the convergence of the method in a finite element framework. An O(h ^1.5)-order superconvergence of the solution derivatives in a discrete H ¹ norm is obtained. Finally, numerical experiments are provided to support the theoretical convergence rate obtained.     

Two order superconvergence of finite element methods for Sobolev equations

We consider finite element methods applied to a class of Sobolev equations inR ^d(d ≥ 1). Global strong superconvergence, which only requires that partitions are quais-uniform, is investigated for the error between the approximate solution and the Ritz-Sobolev projection of the exact solution. Two order superconvergence results are demonstrated inW ^1,p (Ω) andL _p(Ω) for 2 ≤p < ∞.     

Two-grid Methods for Finite Volume Element Approximations of Nonlinear Sobolev Equations

In this article, two-grid methods are studied for solving nonlinear Sobolev equation using the finite volume element method. The methods are based on one coarse grid space and one fine grid space. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H), and the fine grid solution (with grid size h) can be obtained in a single symmetric and linear step. The optimal H¹ error estimates are presented for the proposed methods, which show that the two-grid methods achieve optimal approximation as long as the mesh sizes satisfy h = 𝒪(H³|ln H|). As a result, solving such a large class of nonlinear Sobolev equations will not be much more difficult than solving one linearized equation.     

Two temporal second-order <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-Galerkin mixed finite element schemes for distributed-order fractional sub-diffusion equations

In this paper, a special point is found for the interpolation approximation of the distributed order fractional derivatives to achieve at least second-order accuracy. Then, two H¹-Galerkin mixed finite element schemes combined with the higher accurate interpolation approximation are introduced and analyzed to solve the distributed order fractional sub-diffusion equations. The stable results, which just depend on initial value and source item, are derived. Some a priori estimates with optimal order of convergence both for the unknown function and its flux are established rigorously. It is shown that the H¹-Galerkin mixed finite element approximations have the same rates of the convergence as in the classical mixed finite element method, but without LBB consistency condition and quasiuniformity requirement on the finite element mesh. Finally, some numerical experiments are presented to show the efficiency and accuracy of H¹-Galerkin mixed finite element schemes.     

A Priori Error Estimates for Finite Element Methods for H (2,1)-Elliptic Equations

The convergence of finite element methods for linear elliptic boundary value problems of second and forth order is well understood. In this article, we introduce finite element approximations of some linear semi-elliptic boundary value problem of mixed order on a two-dimensional rectangular domain Q. The equation is of second order in one direction and forth order in the other and appears in the optimal control of parabolic partial differential equations if one eliminates the control and the state (or the adjoint state) in the first order optimality conditions. We establish a regularity result and estimate for the finite element error of conforming approximations of this equation. The finite elements in use have a tensor product structure, in one dimension we use linear, quadratic or cubic Lagrange elements in the other dimension cubic Hermite elements. For these elements, we prove the error bound O(h ² + τ ^k ) in the energy norm and O((h ² + τ ^k )(h ² + τ)) in the L ²(Q)-norm.     

A note onC o Galerkin methods for two-point boundary problems

Summary As is known [4]. theC ^o Galerkin solution of a two-point boundary problem using piecewise polynomial functions, hasO(h ^2k ) convergence at the knots, wherek is the degree of the finite element space. Also, it can be proved [5] that at specific interior points, the Gauss-Legendre points the gradient hasO(h ^k+1) convergence, instead ofO(h ^k ). In this note, it is proved that on any segment there arek–1 interior points where the Galerkin solution is ofO(h ^k+2), one order better than the global order of convergence. These points are the Lobatto points.     

Spline solution of the generalized Burgers'-Fisher equation

A numerical method based on exponential spline and finite difference approximations is developed to solve the generalized Burgers'-Fisher equation. The error analysis, stability and convergence properties of the method are studied via energy method. The method is shown to be unconditionally stable and accurate of orders 𝒪(k?+?kh?+?h ²) and 𝒪(k?+?kh?+?h ⁴). Some test problems are given to demonstrate the applicability of the purposed method numerically. Numerical results verify the theoretical behaviour of convergence properties. The main superiority of the presented scheme is its simplicity and applicability in comparison with the existing well-known methods.     

Superconvergence of the split least‐squares method for second‐order hyperbolic equations

This article studies superconvergence phenomena of the split least‐squares mixed finite element method for second‐order hyperbolic equations. By selecting the least‐squares functional properly, the procedure can be split into two independent symmetric positive definite subprocedures, one of which is for the primitive unknown and the other is for the flux. Based on interpolation operators and an auxiliary projection, superconvergent H¹ error estimates for the primary variable u and L² error estimates for the introduced flux variable σ are obtained under the standard quasiuniform assumptions on finite element partition. A numerical example is given to show the performance of the introduced scheme. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 222‐238, 2014