ERROR ANALYSIS OF A STABILIZED FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES PROBLEM

This paper is devoted to the establishment of sharper $a$ $priori$stability and error estimates of a stabilized finite element method proposed by Barrenechea and Valentin for solving the generalized Stokes problem, which involves a viscosity $\nu$ and a reaction constant $\sigma$. With the establishment of sharper stability estimates and the help of $ad$ $hoc$finite element projections, we can explicitly establish the dependence of error bounds of velocity and pressure on the viscosity $\nu$, the reaction constant $\sigma$, and the mesh size $h$. Our analysis reveals that the viscosity $\nu$ and the reaction constant $\sigma$ respectively act in the numerator position and the denominator position in the error estimates of velocity and pressure in standard norms without any weights. Consequently, the stabilization method is indeed suitable for the generalized Stokes problem with a small viscosity $\nu$ and a large reaction constant $\sigma$. The sharper error estimates agree very well with the numerical results.     

The Finite Element Analysis of the Controlled-Source Electromagnetic Induction Problems by Fractional-Step Projection Method

This paper provides a convergence analysis of a fractional-step projection method for the controlled-source electromagnetic induction problems in heterogenous electrically conduting media by means of finite element approximations. Error estimates in finite time are given. And it is verified that provided the time step $\tau$ is sufficiently small, the proposed algorithm yields for finite time $T$ an error of $\mathcal{O}(h^s+\tau)$) in the $L^2$-norm for the magnetic field $\boldsymbol{H},$ where $h$ is the mesh size and $1/2 < s≤1$.     

Derivative superconvergent points in finite element solutions of Poisson's equation for the serendipity and intermediate families - a theoretical justification

Finite element derivative superconvergent points for the Poisson equation under local rectangular mesh (in the two dimensional case) and local brick mesh (in the three dimensional situation) are investigated. All superconvergent points for the finite element space of any order that is contained in the tensor-product space and contains the intermediate family can be predicted. In case of the serendipity family, the results are given for finite element spaces of order below 7. Any finite element space that contains the complete polynomial space will have at least all superconvergent points of the related serendipity family.

     

半正定两相驱动问题的多步有限体积方法及其理论分析

考虑多维半正定两相驱动方程的初边值问题,在非结构网格上构造多步的迎风有限体积格式,利用微分方程先验估计理论证明了格式的离散模形式的误差估计为D(△t~2 h),其中△t和h分别表示时空步长．数值算例进一步验证了格式的有效性．     

Graded Meshes for Higher Order FEM

A singularly perturbed one-dimensional convection-diffusion problem is solved numerically by the finite element method based on higher order polynomials. Numerical solutions are obtained using S-type meshes with special emphasis on meshes which are graded (based on a mesh generating function) in the fine mesh region. Error estimates in the ε-weighted energy norm are proved. We derive an 'optimal' mesh generating function in order to minimize the constant in the error estimate. Two layer-adapted meshes defined by a recursive formulae in the fine mesh region are also considered and a new technique for proving error estimates for these meshes is presented. The aim of the paper is to emphasize the importance of using optimal meshes for higher order finite element methods. Numerical experiments support all theoretical results.     

CONVERGENCE OF AN ALTERNATING A-φ SCHEME FOR QUASI-MAGNETOSTATIC EDDY CURRENT PROBLEM

We propose in this paper an alternating A-$\phi$ method for the quasi-magnetostatic eddy current problem by means of finite element approximations. Bounds for continuous and discrete error in finite time are given. And it is verified that provided the time step $\tau$ is sufficiently small, the proposed algorithm yields for finite time $T$ an error of $O(h+\tau^{1/2})$ in the $L^2$-norm for the magnetic field $H(= \mu^{-1} \nabla \times A)$, where $h$ is the mesh size, $\mu$ the magnetic permeability.     

A POSTERIORI ERROR ESTIMATES IN ADINI FINITE ELEMENT FOR EIGENVALUE PROBLEMS

In this paper, we discuss a posteriori error estimates of the eigenvalue $\lambda_h$ given by Adini nonconforming finite element. We give an asymptotically exact error estimator of the $\lambda_h$. We prove that the order of convergence of the $\lambda_h$ is just 2 and the $\lambda_h$ converge from below for sufficiently small $h$.     

An Improved Two-Grid Technique for the Nonlinear Time-Fractional Parabolic Equation Based on the Block-Centered Finite Difference Method

A combined scheme of the improved two-grid technique with the block-centered finite difference method is constructed and analyzed to solve the nonlinear time-fractional parabolic equation. This method is considered where the nonlinear problem is solved only on a coarse grid of size $H$ and two linear problems based on the coarse-grid solutions and one Newton iteration is considered on a fine grid of size $h$. We provide the rigorous error estimate, which demonstrates that our scheme converges with order $\mathcal{O}(\Delta t^{2-\alpha}+h^2+H^4)$ on non-uniform rectangular grid. This result indicates that the improved two-grid method can obtain asymptotically optimal approximation as long as the mesh sizes satisfy $h=\mathcal{O}(H^2).$ Finally, numerical tests confirm the theoretical results of the presented method.     

Superconvergent Cluster Recovery Method for the Crouzeix-Raviart Element

In this paper, we propose and numerically investigate a superconvergent cluster recovery (SCR) method for the Crouzeix-Raviart (CR) element. The proposed recovery method reconstructs a $C^0$ linear gradient. A linear polynomial approximation is obtained by a least square fitting to the CR element approximation at certain sample points, and then taken derivatives to obtain the recovered gradient. The SCR recovery operator is superconvergent on uniform mesh of four patterns. Numerical examples show that SCR can produce a superconvergent gradient approximation for the CR element, and provide an asymptotically exact error estimator in the adaptive CR finite element method.     

Derivative superconvergent points in finite element solutions of harmonic functions-- A theoretical justification

Finite element derivative superconvergent points for harmonic functions under local rectangular mesh are investigated. All superconvergent points for the finite element space of any order that is contained in the tensor-product space and contains the intermediate family can be predicted. In the case of the serendipity family, results are given for finite element spaces of order below 6. The results justify the computer findings of Babuska, et al.

     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

TWO-GRID CHARACTERISTIC FINITE VOLUME METHODS FOR NONLINEAR PARABOLIC PROBLEMS＊

In this work, two-grid characteristic finite volume schemes for the nonlinear parabolic problem are considered. In our algorithms, the diffusion term is discretized by the finite volume method, while the temporal differentiation and advection terms are treated by the characteristic scheme. Under some conditions about the coefficients and exact solution, optimal error estimates for the numerical solution are obtained. Furthermore, the two- grid characteristic finite volume methods involve solving a nonlinear equation on coarse mesh with mesh size H, a large linear problem for the Oseen two-grid characteristic finite volume method on a fine mesh with mesh size h = O（H2） or a large linear problem for the Newton two-grid characteristic finite volume method on a fine mesh with mesh size h = 0（I log hll/2H3）. These methods we studied provide the same convergence rate as that of the characteristic finite volume method, which involves solving one large nonlinear problem on a fine mesh with mesh size h. Some numerical results are presented to demonstrate the efficiency of the proposed methods.     

Polynomial preserving recovery for quadratic elements on anisotropic meshes

Polynomial preserving gradient recovery technique under anisotropic meshes is further studied for quadratic elements. The analysis is performed for highly anisotropic meshes where the aspect ratios of element sides are unbounded. When the mesh is adapted to the solution that has significant changes in one direction but very little, if any, in another direction, the recovered gradient can be superconvergent. The results further explain why recovery type error estimator is robust even under nonstandard and highly distorted meshes. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Error estimates on anisotropic elements for functions in weighted Sobolev spaces

In this paper we prove error estimates for a piecewise average interpolation on anisotropic rectangular elements, i.e., rectangles with sides of different orders, in two and three dimensions.

Our error estimates are valid under the condition that neighboring elements have comparable size. This is a very mild assumption that includes more general meshes than those allowed in previous papers. In particular, strong anisotropic meshes arising naturally in the approximation of problems with boundary layers fall under our hypotheses.

Moreover, we generalize the error estimates allowing on the right-hand side some weighted Sobolev norms. This extension is of interest in singularly perturbed problems.

Finally, we consider the approximation of functions vanishing on the boundary by finite element functions with the same property, a point that was not considered in previous papers on average interpolations for anisotropic elements.

As an application we consider the approximation of a singularly perturbed reaction-diffusion equation and show that, as a consequence of our results, almost optimal order error estimates in the energy norm, valid uniformly in the perturbation parameter, can be obtained.

     

GLOBAL SUPERCONVERGENCE OF THE MIXED FINITEELEMENT METHODS FOR 2-D MAXWELL EQUATIONS

Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.     

Numerical Analysis of the Allen-Cahn Equation with Coarse Meshes

In this paper, we consider the finite difference semi-discretization of the Allen-Cahn equation with the diffuse interface parameter $\varepsilon$. While it is natural to make the mesh size parameter $h$ smaller than $\varepsilon$, it is desirable that $h$ is as big as possible in view of computational costs. In fact, when $h$ is bigger than $\varepsilon$ (i.e., the mesh is relatively coarse), it is observed that the numerical solution does not move at all. The purpose of this paper is to clarify the mechanism of this phenomenon. We will prove that the numerical solution converges to that of the ordinary equation without the diffusion term if $h$ is bigger than $\varepsilon$. Numerical examples are presented to support the result.     

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

The numerical approximation by a lower‐order anisotropic nonconforming finite element on appropriately graded meshes are considered for solving semisingular 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. Copyright © 2007 John Wiley & Sons, Ltd.     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS II： TIME DISCRETIZATION

In this article we consider the fully discrete two-level finite element Galerkin method for the two-dimensional nonstationary incompressible Navier-Stokes equations. This method consists in dealing with the fully discrete nonlinear Navier-Stokes problem on a coarse mesh with width $H$ and the fully discrete linear generalized Stokes problem on a fine mesh with width $h << H$. Our results show that if we choose $H=O(h^{1/2}$) this method is as the same stability and convergence as the fully discrete standard finite element Galerkin method which needs dealing with the fully discrete nonlinear Navier-Stokes problem on a fine mesh with width $h$. However, our method is cheaper than the standard fully discrete finite element Galerkin method.     

Mathematical analysis of Zienkiewicz—Zhu's derivative patch recovery technique

Zienkiewicz-Zhu's derivative patch recovery technique is analyzed for general quadrilateral finite elements. Under certain regular conditions on the meshes, the arithmetic mean of the absolute error of the recovered gradient at the nodal points is superconvergent for the second-order elliptic operators. For rectangular meshes and the Laplacian, the recovered gradient is superconvergent in the maximum norm at the nodal points. Furthermore, it is proved for a model two-point boundary-value problem that the recovery technique results in an “ultra-convergent” derivative recovery at the nodal points for quadratic finite elements when uniform meshes are used. © 1996 John Wiley & Sons, Inc.     

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.