Local projection FEM stabilization for the time‐dependent incompressible Navier–Stokes problem

We consider conforming finite element (FE) approximations of the time‐dependent, incompressible Navier–Stokes problem with inf‐sup stable approximation of velocity and pressure. In case of high Reynolds numbers, a local projection stabilization method is considered. In particular, the idea of streamline upwinding is combined with stabilization of the divergence‐free constraint. For the arising nonlinear semidiscrete problem, a stability and convergence analysis is given. Our approach improves some results of a recent paper by Matthies and Tobiska (IMA J. Numer. Anal., to appear) for the linearized model and takes partly advantage of the analysis in Burman and Fernández, Numer. Math. 107 (2007), 39–77 for edge‐stabilized FE approximation of the Navier–Stokes problem. Some numerical experiments complement the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1224–1250, 2015     

On the relationship between finite volume and finite element methods applied to the Stokes equations

We investigate the relationship between finite volume and finite element approximations for the lower‐order elements, both conforming and nonconforming for the Stokes equations. These elements include conforming, linear velocity‐constant pressure on triangles, conforming bilinear velocity‐constant pressure on rectangles and their macro‐element versions, and nonconforming linear velocity‐constant pressure on triangles and nonconforming rotated bilinear velocity‐constant pressure on rectangles. By applying the relationship between the two methods, we obtain the convergence finite volume solutions for the Stokes equations. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 440–453, 2001.     

On the convergence of nonnested multigrid methods with nested spaces on coarse grids

Multigrid methods for discretized partial differential problems using nonnested conforming and nonconforming finite elements are here defined in the general setting. The coarse‐grid corrections of these multigrid methods make use of different finite element spaces from those on the finest grid. In general, the finite element spaces on the finest grid are nonnested, while the spaces are nested on the coarse grids. An abstract convergence theory is developed for these multigrid methods for differential problems without full elliptic regularity. This theory applies to multigrid methods of nonnested conforming and nonconforming finite elements with the coarse‐grid corrections established on nested conforming finite element spaces. Uniform convergence rates (independent of the number of grid levels) are obtained for both the V and W‐cycle methods with one smoothing on all coarse grids and with a sufficiently large number of smoothings solely on the finest grid. In some cases, these uniform rates are attained even with one smoothing on all grids. The present theory also applies to multigrid methods for discretized partial differential problems using mixed finite element methods. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 265–284, 2000     

Conforming and nonconforming harmonic finite elements

ABSTRACT

Instead of using the full polynomial space, a conforming and a nonconforming finite element methods are designed where only harmonic polynomials (a much smaller space) are employed in the computation. The conforming quadratic harmonic polynomial finite element is defined only on a special triangular grid. The nonconforming quadratic harmonic finite element is defined on general triangular grids. The optimal order of convergence is proved for both finite element methods, and confirmed by numerical computations. In addition, numerical comparisons with the standard conforming and nonconforming finite elements are presented.     

非定常对流扩散问题的非协调局部投影有限元方法

本文将近年来基于协调有限元逼近提出的涡旋粘性法推广应用到非协调有限元逼近,对非定常的对流占优扩散问题,空间采用非协调Crouzeix-Raviart元逼近,时间用Crank-Nicolson差分离散格式,提出了Crank-Nicolson差分-局部投影法稳定化有限元格式,我们对稳定性和误差估计给出了详细的分析,得出了最...     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Error analysis of a fully discrete finite element variational multiscale method for time‐dependent incompressible Navier–Stokes equations

A finite element variational multiscale method based on two local Gauss integrations is applied to solve numerically the time‐dependent incompressible Navier–Stokes equations. A significant feature of the method is that the definition of the stabilization term is derived via two local Guass integrations at element level, making it more efficient than the usual projection‐based variational multiscale methods. It is computationally cheap and gives an accurate approximation to the quantities sought. Based on backward Euler and Crank–Nicolson schemes for temporal discretization, we derive error bounds of the fully discrete solution which are first and second order in time, respectively. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

热传导方程一个新的H~1-Galerkin非协调混合有限元格式

对热传导方程提出了一个新的H~1-Galerkin非协调混合有限元格式,其逼近空间不需满足LBB相容性条件,且在不引进传统的Rutz投影的情况下,得到了与以往协调有限元方法相同的L~2-模和H~1-模的误差估计.     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Stabilized Crouzeix–Raviart element for Darcy–Forchheimer model

In this article, a stabilized mixed finite element method for steady Darcy–Forchheimer flow is introduced, in which the velocity and pressure are approximated by nonconforming Crouzeix–Raviart element and piecewise constant, respectively. A discrete inf‐sup condition and a priori error estimates are derived. An iterative scheme is given for practical computation. Finally, some numerical examples are carried out to verify the theoretical analysis and a comparison between two discretizations is given to demonstrate that one of the discretizations has better properties. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1568–1588, 2015     

A unified framework for a posteriori error estimation for the Stokes problem

In this paper, a unified framework for a posteriori error estimation for the Stokes problem is developed. It is based on $[H^1_0(\Omega )]^d$ -conforming velocity reconstruction and $\underline{\varvec{H}}(\mathrm{div},\Omega )$ -conforming, locally conservative flux (stress) reconstruction. It?gives guaranteed, fully computable global upper bounds as well as local lower bounds on the energy error. In order to apply this framework to a given numerical method, two simple conditions need to be checked. We show how to do this for various conforming and conforming stabilized finite element methods, the discontinuous Galerkin method, the Crouzeix–Raviart nonconforming finite element method, the mixed finite element method, and a general class of finite volume methods. The tools developed and used include a new simple equilibration on dual meshes and the solution of local Poisson-type Neumann problems by the mixed finite element method. Numerical experiments illustrate the theoretical developments.     

A combined finite volume–finite element scheme for the discretization of strongly nonlinear convection–diffusion–reaction problems on nonmatching grids

We propose and analyze in this paper a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the time evolution, convection, reaction, and source terms on a given grid, which can be nonmatching and can contain nonconvex elements, by means of the cell‐centered finite volume method. To discretize the diffusion term, we construct a conforming simplicial mesh with the vertices given by the original grid and use the conforming piecewise linear finite element method. In this way, the scheme is fully consistent and the discrete solution is naturally continuous across the interfaces between the subdomains with nonmatching grids, without introducing any supplementary equations and unknowns or using any interpolation at the interfaces. We allow for general inhomogeneous and anisotropic diffusion–dispersion tensors, propose two variants corresponding respectively to arithmetic and harmonic averaging, and use the local Péclet upstream weighting in order to only add the minimal numerical diffusion necessary to avoid spurious oscillations in the convection‐dominated case. The scheme is robust, efficient since it leads to positive definite matrices and one unknown per element, locally conservative, and satisfies the discrete maximum principle under the conditions on the simplicial mesh and the diffusion tensor usual in the finite element method. We prove its convergence using a priori estimates and the Kolmogorov relative compactness theorem and illustrate its behavior on a numerical experiment. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

An adaptive least squares mixed finite element method for the stress‐displacement formulation of linear elasticity

A least‐squares mixed finite element method for linear elasticity, based on a stress‐displacement formulation, is investigated in terms of computational efficiency. For the stress approximation quadratic Raviart‐Thomas elements are used and these are coupled with the quadratic nonconforming finite element spaces of Fortin and Soulie for approximating the displacement. The local evaluation of the least‐squares functional serves as an a posteriori error estimator to be used in an adaptive refinement algorithm. We present computational results for a benchmark test problem of planar elasticity including nearly incompressible material parameters in order to verify the effectiveness of our adaptive strategy. For comparison, conforming quadratic finite elements are also used for the displacement approximation showing convergence orders similar to the nonconforming case, which are, however, not independent of the Lamé parameters. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Superconvergence of a two‐grid method for fourth‐order dispersive‐dissipative wave equations

In this paper, the superconvergence analysis of a two‐grid method (TGM) with low‐order finite elements is presented for the fourth‐order dispersive‐dissipative wave equations for a second order fully discrete scheme. The superclose estimates in the H¹‐norm on the two grids are obtained by the combination technique of the interpolation and Ritz projection. Then, with the help of the interpolated postprocessing technique, the global superconvergence properties are deduced. Finally, numerical results are provided to show the performance of the proposed TGM for conforming bilinear element and nonconforming element, respectively. It shows that the TGM is an effective method to the problem considered of our paper compared with the traditional Galerkin finite element method (FEM).     

Local projection stabilization for advection–diffusion–reaction problems: One-level vs. two-level approach

Local projection stabilization (LPS) of finite element methods is a new technique for the numerical solution of transport-dominated problems. The main aim of this paper is a critical discussion and comparison of the one- and two-level approaches to LPS for the linear advection–diffusion–reaction problem. Moreover, the paper contains several other novel contributions to the theory of LPS. In particular, we derive an error estimate showing not only the usual error dependence on the mesh width but also on the polynomial degree of the finite element space. Based on this error estimate, we propose a definition of the stabilization parameter depending on the data of the solved problem. Unlike other papers on LPS methods, we observe that the consistency error may deteriorate the convergence order. Finally, we explain the relation between the LPS method and residual-based stabilization techniques for simplicial finite elements.     

A posteriori error estimation for the Stokes–Darcy coupled problem on anisotropic discretization

This paper presents an a posteriori error analysis for the stationary Stokes–Darcy coupled problem approximated by finite element methods on anisotropic meshes in or 3. Korn's inequality for piecewise linear vector fields on anisotropic meshes is established and is applied to non‐conforming finite element method. Then the existence and uniqueness of the approximation solution are deduced for non‐conforming case. With the obtained finite element solutions, the error estimators are constructed and based on the residual of model equations plus the stabilization terms. The lower error bound is proved by means of bubble functions and the corresponding anisotropic inverse inequalities. In order to prove the upper error bound, it is vital that an anisotropic mesh corresponds to the anisotropic function under consideration. To measure this correspondence, a so‐called matching function is defined, and its discussion shows it to be useful tool. With its help, the upper error bound is shown by means of the corresponding anisotropic interpolation estimates and a special Helmholtz decomposition in both media. Copyright © 2016 John Wiley & Sons, Ltd.     

Stabilized spectral element approximation for the Navier–Stokes equations

The conforming spectral element methods are applied to solve the linearized Navier–Stokes equations by the help of stabilization techniques like those applied for finite elements. The stability and convergence analysis is carried out and essential numerical results are presented demonstrating the high accuracy of the method as well as its robustness. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 115–141, 1998     

Analysis of a finite volume element method for the Stokes problem

In this paper we propose a stabilized conforming finite volume element method for the Stokes equations. On stating the convergence of the method, optimal a priori error estimates in different norms are obtained by establishing the adequate connection between the finite volume and stabilized finite element formulations. A superconvergence result is also derived by using a postprocessing projection method. In particular, the stabilization of the continuous lowest equal order pair finite volume element discretization is achieved by enriching the velocity space with local functions that do not necessarily vanish on the element boundaries. Finally, some numerical experiments that confirm the predicted behavior of the method are provided.     

An immersed finite element space and its approximation capability

This article discusses an immersed finite element (IFE) space introduced for solving a second‐order elliptic boundary value problem with discontinuous coefficients (interface problem). The IFE space is nonconforming and its partition can be independent of the interface. The error estimates for the interpolation of a function in the usual Sobolev space indicate that this IFE space has an approximation capability similar to that of the standard conforming linear finite element space based on body‐fit partitions. Numerical examples of the related finite element method based on this IFE space are provided. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 338–367, 2004     

Two-level stabilized nonconforming finite element method for the Stokes equations

In this article, we present a new two-level stabilized nonconforming finite elements method for the two dimensional Stokes problem. This method is based on a local Gauss integration technique and the mixed nonconforming finite element of the NCP ₁ — P ₁ pair (nonconforming linear element for the velocity, conforming linear element for the pressure). The two-level stabilized finite element method involves solving a small stabilized Stokes problem on a coarse mesh with mesh size H and a large stabilized Stokes problem on a fine mesh size h = H/3. Numerical results are presented to show the convergence performance of this combined algorithm.