High accuracy analysis of a new nonconforming mixed finite element scheme for Sobolev equations

A nonconforming mixed finite element scheme is proposed for Sobolev equations based on a new mixed variational form under semi-discrete and Euler fully-discrete schemes. The corresponding optimal convergence error estimates and superclose property are obtained without using Ritz projection, which are the same as the traditional mixed finite elements. Furthemore, the global superconvergence is obtained through interpolation postprocessing technique. The numerical results show the validity of the theoretical analysis.     

Sobolev方程一个新的H~1-Galerkin混合有限元分析

研究了Sobolev方程的H~1-Galerkin混合有限元方法.利用不完全双二次元Q_2~-和一阶BDFM元,建立了一个新的混合元模式,通过Bramble-Hilbert引理,证明了单元对应的插值算子具有的高精度结果.进一步,对于半离散和向后欧拉全离散格式,分别导出了原始变量u在H~1-模和中间变量p在H(div)-模意义下的超逼近性质.     

Analysis of least-squares mixed finite element methods for nonlinear nonstationary convection-diffusion problems

Some least-squares mixed finite element methods for convection-diffusion problems, steady or nonstationary, are formulated, and convergence of these schemes is analyzed. The main results are that a new optimal a priori error estimate of a least-squares mixed finite element method for a steady convection-diffusion problem is developed and that four fully-discrete least-squares mixed finite element schemes for an initial-boundary value problem of a nonlinear nonstationary convection-diffusion equation are formulated. Also, some systematic theories on convergence of these schemes are established.

     

Time discontinuous Galerkin space-time finite element method for nonlinear Sobolev equations

This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introducing the equivalent integral equation of the primal equation. The proposed scheme does not explicitly include the jump terms in time, which represent the discontinuity characteristics of approximate solution. And then the complexity of the theoretical analysis is reduced. The existence and uniqueness of the approximate solution and the stability of the scheme are proved. The optimalorder error estimates in L²(H¹) and L²(L²) norms are derived. These estimates are valid under weak restrictions on the space-time mesh, namely, without the condition k_n≥ch², which is necessary in traditional space-time discontinuous Galerkin methods. Numerical experiments are presented to verify the theoretical results.     

Nonconforming elements in least-squares mixed finite element methods

In this paper we analyze the finite element discretization for the first-order system least squares mixed model for the second-order elliptic problem by means of using nonconforming and conforming elements to approximate displacement and stress, respectively. Moreover, on arbitrary regular quadrilaterals, we propose new variants of both the rotated nonconforming element and the lowest-order Raviart-Thomas element.

     

A splitting positive definite mixed element method for second‐order hyperbolic equations

In this article, we establish a new mixed finite element procedure, in which the mixed element system is symmetric positive definite, to solve the second‐order hyperbolic equations. The convergence of the mixed element methods with continuous‐ and discrete‐time scheme is proved. And the corresponding error estimates are given. Finally some numerical results are presented. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A reduced-order extrapolating space–time continuous finite element method for the 2D Sobolev equation

This paper mainly concerns with the order reduction to the coefficient vectors of the classical space–time continuous finite element (STCFE) solutions for a two-dimensional Sobolev equation. The classical STCFE model is first constructed for the governing equation, and the theoretical results of the existence, stability, and convergence are provided for the STCFE solutions. We then employ a proper orthogonal decomposition to develop a reduced-order extrapolating STCFE (ROESTCFE) vector model with the lower dimension, and demonstrate the existence, stability, and convergence for the ROESTCFE solutions by the matrix means, resulting in the very concise and flexible theoretical analysis. Lastly, we examine the effectiveness of the developed ROESTCFE model by several numerical tests. It is shown that the ROESTCFE method is computationally very cheap in actual applications.     

Convergence and stability of two-level penalty mixed finite element method for stationary Navier-Stokes equations

The two-level penalty mixed finite element method for the stationary Navier-Stokes equations based on Taylor-Hood element is considered in this paper. Two algorithms are proposed and analyzed. Moreover, the optimal stability analysis and error estimate for these two algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms.     

Error Estimates for Mixed Finite Element Methods for Sobolev Equation

1 IntroductionLet fl be a bounded domain in R2 with Lipschitz continuous boundaxy 0fl. For thed0 < T < co, we consider the fo1lowing initial-boun'lar}-ralue problem for thc Sobolevequation:where ut denotes the time derivative of the function (1. Vu denotes the gradient of thefunction u, and divv denotes the divergence of the vect{Jr tulued function v, a1 b1, f, anduo are known functions.The standard finite element method for (1.1) (1.3) llas received considerable attentionand is well studied…     

New mixed finite elements for plane elasticity and Stokes equations

We consider mixed finite elements for the plane elasticity system and the Stokes equation. For the unmodified Hellinger-Reissner formulation of elasticity in which the stress and displacement fields are the primary unknowns, we derive two new nonconforming mixed finite elements of triangle type. Both elements use piecewise rigid motions to approximate the displacement and piecewise polynomial functions to approximate the stress, where no vertex degrees of freedom are involved. The two stress finite element ...     

Crank–Nicolson least-squares Galerkin procedures for parabolic integro-differential equations

Two Crank–Nicolson least-squares Galerkin finite element schemes are formulated to solve parabolic integro-differential equations. The advantage of this method is that it is not subject to the LBB condition. The convergence analysis shows that the methods yield the approximate solutions with optimal accuracy in H(div; Ω) × H¹(Ω) and (L²(Ω))² × L²(Ω), respectively. Moreover, the two methods both get the approximate solutions with second-order accuracy in time increment.     

A remark on split least-squares mixed element procedures for pseudo-parabolic equations

In this paper, we introduce two novel split least-squares mixed element procedures for pseudo-parabolic equations. By selecting the least-squares functional properly, each procedure can be split into two independent symmetric positive definite sub-procedures. One of sub-procedures is for the primitive unknown variable u, which is the same as the standard Galerkin finite element procedure and the other is for the introduced flux variable σ. Optimal order error estimates are developed. A numerical example is given to show the efficiency of the introduced schemes.     

Least‐squares mixed finite element methods for the RLW equations

A least‐squares mixed finite element (LSMFE) schemes are formulated to solve the 1D regularized long wave (RLW) equations and the convergence is discussed. The L² error estimates of LSMFE methods for RLW equations under the standard regularity assumption on the finite element partition are given.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A mixed nonconforming finite element for linear elasticity

This article considers a mixed finite element method for linear elasticity. It is based on a modified mixed formulation that enforces the continuity of the stress weakly by adding a jump term of the approximated stress on interior edges. The symmetric stress are approximated by nonconforming linear elements and the displacement by piecewise constants. We establish ??(h) error bound in the (broken) L² norm for the divergence of the stress and ??(h) error bound in the L² norm for both the displacement and the stress tensor. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

A Newton linearized compact finite difference scheme for one class of Sobolev equations

In this article, a Newton linearized compact finite difference scheme is proposed to numerically solve a class of Sobolev equations. The unique solvability, convergence, and stability of the proposed scheme are proved. It is shown that the proposed method is of order 2 in temporal direction and order 4 in spatial direction. Moreover, compare to the classical extrapolated Crank‐Nicolson method or the second‐order multistep implicit–explicit methods, the proposed scheme is easier to be implemented as it only requires one starting value. Finally, numerical experiments on one and two‐dimensional problems are presented to illustrate our theoretical results.     

Error estimates for mixed finite element approximations of the viscoelasticity wave equation

This paper studies mixed finite element approximations to the solution of the viscoelasticity wave equation. Two new transformations are introduced and a corresponding system of first‐order differential‐integral equations is derived. The semi‐discrete and full‐discrete mixed finite element methods are then proposed for the problem based on the Raviart–Thomas–Nedelec spaces. The optimal error estimates in L²‐norm are obtained for the semi‐discrete and full‐discrete mixed approximations of the general viscoelasticity wave equation. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

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 posteriori error estimates with the finite element method of lines for a Sobolev equation

A posteriori error estimates for semidiscrete finite element methods for a nonlinear Sobolev equation are considered. The error estimates are obtained by solving local nonlinear or linear pseudo‐parabolic equations for corrections to the solution on each element. The ratios of these estimates and the true errors are proved to converge to 1, implying that the estimates can be used as indicators in adaptive schemes for the problem. Numerical results underline our theoretical results. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005