EQ 1 rot nonconforming finite element approximation to Signorini problem

In this paper, we apply EQ rot 1 nonconforming finite element to approximate Signorini problem. If the exact solution u∈H5/2(Ω), the error estimate of order O(h) about the broken energy norm is obtained for quadrilateral meshes satisfying regularity assumption and bi-section condition. Furthermore, the superconvergence results of order O(h3/2) are derived for rectangular meshes. Numerical results are presented to confirm the considered theory.     

Superconvergence analysis and extrapolation of quasi-Wilson nonconforming finite element method for nonlinear Sobolev equations

Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equa-tions is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches. One is that (▽h ( u-Ihu )1, ▽hvh) h may be estimated as order O ( h2 ) when u ∈ H3 (Ω), where Ihu denotes the bilinear interpolation of u , vh is a polynomial belongs to quasi-Wilson finite element space and ▽h denotes the piecewise defined gradient operator, h is the mesh size tending to zero. The other is that the consistency error of this element is of order O ( h2 ) /O ( h3 ) in broken H 1-norm, which is one/two order higher than its interpolation error when u ∈ H3 (Ω) /H4 (Ω). Then we derive the optimal order error estimate and su- perclose property via mean-value method and the known high accuracy result of bilinear element. Furthermore, we deduce the global superconvergence through interpolation post processing technique. At last, an extrapola- tion result of order O ( h3 ), two order higher than traditional error estimate, is obtained by constructing a new suitable extrapolation scheme.     

A NONCONFORMING QUADRILATERAL FINITE ELEMENT APPROXIMATION TO NONLINEAR SCHRDINGER EQUATION

In this article, a nonconforming quadrilateral element(named modified quasiWilson element) is applied to solve the nonlinear schr¨odinger equation(NLSE). On the basis of a special character of this element, that is, its consistency error is of order O(h~3) for broken H1-norm on arbitrary quadrilateral meshes, which is two order higher than its interpolation error, the optimal order error estimate and superclose property are obtained. Moreover,the global superconvergence result is deduced with the help of interpolation postprocessing technique. Finally, some numerical results are provided to verify the theoretical analysis.     

Extrapolation of the Nédélec element for the Maxwell equations by the mixed finite element method

In this paper, we use the integral-identity argument to obtain asymptotic error expansions for the mixed finite element approximation of the Maxwell equations on a rectangular mesh. The extrapolation method is applied to improve the accuracy of the approximation via an interpolation postprocessing technique. With the extrapolation, the approximation accuracy can be improved from O(h) to O(h ⁴) in the L ²-norm. Illustrative numerical results are given to demonstrate the higher order accuracy of the extrapolation method. This research was supported by the National Natural Science Foundation of China (No.10471103), Social Science Foundation of the Ministry of Education of China (06JA630047), Tianjin Natural Science Foundation (07JCYBJC14300).     

Nonconforming finite element approximations of the Steklov eigenvalue problem and its lower bound approximations

The paper deals with error estimates and lower bound approximations of the Steklov eigenvalue problems on convex or concave domains by nonconforming finite element methods. We consider four types of nonconforming finite elements: Crouzeix-Raviart, Q ₁ ^rot , EQ ₁ ^rot and enriched Crouzeix-Raviart. We first derive error estimates for the nonconforming finite element approximations of the Steklov eigenvalue problem and then give the analysis of lower bound approximations. Some numerical results are presented to validate our theoretical results.     

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.     

A low order anisotropic nonconforming characteristic finite element method for a convection-dominated transport problem

A low order anisotropic nonconforming rectangular finite element method for the convection-diffusion problem with a modified characteristic finite element scheme is studied in this paper. The O(h²) order error estimate in L²-norm with respect to the space, one order higher than the expanded characteristic-mixed finite element scheme with order O(h), and the same as the conforming case for a modified characteristic finite element scheme under regular meshes, is obtained by use of some distinct properties of the interpolation operator and the mean value technique, instead of the so-called elliptic projection, which is an indispensable tool in the convergence analysis of the previous literature. Lastly, some numerical results of the element are provided to verify our theoretical analysis.     

Convergence of H 1-Galerkin mixed finite element method for parabolic problems with reduced regularity on initial data

We study the convergence of H ¹-Galerkin mixed finite element method for parabolic problems in one space dimension. Both semi-discrete and fully discrete schemes are analyzed assuming less regularity on initial data. More precisely, for the spatially discrete scheme, error estimates of order \(\mathcal{O}\) (h ² t ^?1/2) for positive time are established assuming the initial function p ₀ ∈ H ²(Ω) ∩ H ₀ ¹ (Ω). Further, we use energy technique together with parabolic duality argument to derive error estimates of order \(\mathcal{O}\) (h ² t ^?1) when p ₀ is only in H ₀ ¹ (Ω). A discrete-in-time backward Euler method is analyzed and almost optimal order error bounds are established.     

A second order nonconforming rectangular finite element method for approximating Maxwell’s equations

The main objective of this paper is to present a new rectangular nonconforming finite element scheme with the second order convergence behavior for approximation of Maxwell’s equations.Then the corresponding optimal error estimates are derived.The difficulty in construction of this finite element scheme is how to choose a compatible pair of degrees of freedom and shape function space so as to make the consistency error due to the nonconformity of the element being of order O(h 3 ) ,properly one order higher than that of its interpolation error O(h 2 ) in the broken energy norm,where h is the subdivision parameter tending to zero.     

Spatial High Accuracy Analysis of FEM for Two-dimensional Multi-term Time-fractional Diffusion-wave Equations

In this paper, high-order numerical analysis of finite element method (FEM) is presented for two-dimensional multi-term time-fractional diffusion-wave equation (TFDWE). First of all, a fully-discrete approximate scheme for multi-term TFDWE is established, which is based on bilinear FEM in spatial direction and Crank-Nicolson approximation in temporal direction, respectively. Then the proposed scheme is proved to be unconditionally stable and convergent. And then, rigorous proofs are given here for superclose properties in H¹?norm and temporal convergence in L²-norm with order \(O(h^{2}+\tau^{3-\alpha})\), where h and τ are the spatial size and time step, respectively. At the same time, theoretical analysis of global superconvergence in H¹-norm is derived by interpolation postprocessing technique. At last, numerical example is provided to demonstrate the theoretical analysis.     

Computing the lower and upper bounds of Laplace eigenvalue problem: by combining conforming and nonconforming finite element methods

We introduce some ways to compute the lower and upper bounds of the Laplace eigenvalue problem.By using the special nonconforming finite elements,i.e.,enriched Crouzeix-Raviart element and extended Q1ro t,we get the lower bound of the eigenvalue.Additionally,we use conforming finite elements to do the postprocessing to get the upper bound of the eigenvalue,which only needs to solve the corresponding source problems and a small eigenvalue problem if higher order postprocessing method is implemented.Thus,we can obtain the lower and upper bounds of the eigenvalues simultaneously by solving eigenvalue problem only once.Some numerical results are also presented to demonstrate our theoretical analysis.     

Nonconforming quadrilateral finite element method for Ginzburg–Landau equation

Nonconforming finite element method is studied for a linearized backward fully-discrete scheme of the Ginzburg–Landau equation with the quadrilateral element. The unconditional convergent result of order O(h + τ) in the broken H¹-norm is deduced rigorously based on a splitting technique, by which the ratio between the subdivision parameter h and the time step τ is removed. Furthermore, numerical results are provided to confirm the theoretical analysis. The analysis developed herein can be regarded as a framework to deal with the unconditional convergent analysis of the Ginzburg–Landau equation for other known low order nonconforming elements.     

Multi-neighboring grids schemes for solving PDE eigen-problems

Instead of most existing postprocessing schemes,a new preprocessing approach,called multineighboring grids(MNG),is proposed for solving PDE eigen-problems on an existing grid G(Δ).The linear or multi-linear element,based on box-splines,are taken as the frst stage Kh1Uh=λh1Mh1Uh.In this paper,the j-th stage neighboring-grid scheme is defned asKh jUh=λh j Mh jUh,where Kh j:=Mh j 1Kh1and Mh jUh is to be found as a better mass distribution over the j-th stage neighboring-gridG(Δ),and Kh jcan be seen as an expansion of Kh1on the j-th neighboring-grid with respect to the(j 1)-th mass distribution Mh j 1.It is shown that for an ODE model eigen-problem,the j-th stage scheme with 2j-th order B-spline basis can reach2j-th order accuracy and even(2j+2)-th order accuracy by perturbing the mass matrix.The argument can be extended to high dimensions with separable variable cases.For Laplace eigen-problems with some 2-D and 3-D structured uniform grids,some 2j-th order schemes are presented for j 3.     

Superconvergence and recovery type a posteriori error estimation for hybrid stress finite element method

Superconvergence and recovery type a posteriori error estimators are analyzed for Pian and Sumihara's 4-node hybrid stress quadrilateral finite element method for linear elasticity problems. Superconvergence of order O(h~(1+min){α,1}) is established for both the displacement approximation in H~1-norm and the stress approximation in L~2-norm under a mesh assumption, where α 0 is a parameter characterizing the distortion of meshes from parallelograms to quadrilaterals. Recovery type approximations for the displacement gradients and the stress tensor are constructed, and a posteriori error estimators based on the recovered quantities are shown to be asymptotically exact. Numerical experiments confirm the theoretical results.     

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.     

Mixed finite element methods for general quadrilateral grids

We study a new mixed finite element of lowest order for general quadrilateral grids which gives optimal order error in the H(div)-norm. This new element is designed so that the H(div)-projection Π_h satisfies ∇ · Π_h = P_hdiv. A rigorous optimal order error estimate is carried out by proving a modified version of the Bramble-Hilbert lemma for vector variables. We show that a local H(div)-projection reproducing certain polynomials suffices to yield an optimal L²-error estimate for the velocity and hence our approach also provides an improved error estimate for original Raviart-Thomas element of lowest order. Numerical experiments are presented to verify our theory.     

Extrapolation and superconvergence of the Steklov eigenvalue problem

On the basis of a transform lemma, an asymptotic expansion of the bilinear finite element is derived over graded meshes for the Steklov eigenvalue problem, such that the Richardson extrapolation can be applied to increase the accuracy of the approximation, from which the approximation of O(h ^3.5) is obtained. In addition, by means of the Rayleigh quotient acceleration technique and an interpolation postprocessing method, the superconvergence of the bilinear finite element is presented over graded meshes for the Steklov eigenvalue problem, and the approximation of O(h ³) is gained. Finally, numerical experiments are provided to demonstrate the theoretical results.     

Approximation and eigenvalue extrapolation of Stokes eigenvalue problem by nonconforming finite element methods

In this paper we analyze the stream function-vorticity-pressure method for the Stokes eigenvalue problem. Further, we obtain full order convergence rate of the eigenvalue approximations for the Stokes eigenvalue problem based on asymptotic error expansions for two nonconforming finite elements, Q ₁^rot and EQ ₁^rot. Using the technique of eigenvalue error expansion, the technique of integral identities and the extrapolation method, we can improve the accuracy of the eigenvalue approximations. This project is supported in part by the National Natural Science Foundation of China (10471103) and is subsidized by the National Basic Research Program of China under the grant 2005CB321701.     

Unconditional superconvergent analysis of a new mixed finite element method for Ginzburg–Landau equation

In this article, unconditional superconvergent analysis of a linearized fully discrete mixed finite element method is presented for a class of Ginzburg–Landau equation based on the bilinear element and zero‐order Nédélec's element pair (Q₁₁/Q₀₁ × Q₁₀). First, a time‐discrete system is introduced to split the error into temporal error and spatial error, and the corresponding error estimates are deduced rigorously. Second, the unconditional superclose and optimal estimate of order O(h² + τ) for u in H¹‐norm and p = ?u in L²‐norm are derived respectively without the restrictions on the ratio between h and τ, where h is the subdivision parameter and τ, the time step. Third, the global superconvergent results are obtained by interpolated postprocessing technique. Finally, some numerical results are carried out to confirm the theoretical analysis.     

Convergence and superconvergence analysis of an anisotropic nonconforming finite element methods for singularly perturbed reaction-diffusion problems

The numerical approximation by a lower order anisotropic nonconforming finite element on appropriately graded meshes are considered for solving singular perturbation problems. The quasi-optimal order error estimates are proved in the ε-weighted H¹-norm valid uniformly, up to a logarithmic factor, in the singular perturbation parameter. By using the interpolation postprocessing technique, the global superconvergent error estimates in ε-weighted H¹-norm are obtained. Numerical experiments are given to demonstrate validity of our theoretical analysis.