Stabilized finite element methods based on multiscale enrichment for the Stokes problem

This work concerns the development of stabilized finite element methods for the Stokes problem considering nonstable different (or equal) order of velocity and pressure interpolations. The approach is based on the enrichment of the standard polynomial space for the velocity component with multiscale functions which no longer vanish on the element boundary. On the other hand, since the test function space is enriched with bubble-like functions, a Petrov--Galerkin approach is employed. We use such a strategy to propose stable variational formulations for continuous piecewise linear in velocity and pressure and for piecewise linear/piecewise constant interpolation pairs. Optimal order convergence results are derived and numerical tests validate the proposed methods

Stabilization arising from PGEM : a review and further developments

The aim of this paper is twofold. First, we review the recent Petrov-Galerkin enriched method (PGEM) to stabilize numerical solutions of BVP's in primal and mixed forms. Then, we extend such enrichment technique to a mixed singularly perturbed problem, namely, the generalized Stokes problem, and focus on a stabilized finite element method arising in a natural way after performing static condensation. The resulting stabilized method is shown to lead to optimal convergences, and afterward, it is numerically validated

On the stability of bubble functions and a stabilized mixed finite element formulation for the Stokes problem

In this paper we investigate the relationship between stabilized and enriched finite element formulations for the Stokes problem. We also present a new stabilized mixed formulation for which the stability parameter is derived purely by the method of weighted residuals. This new formulation allows equal order interpolation for the velocity and pressure fields. Finally, we show by counterexample that a direct equivalence between subgrid-based stabilized finite element methods and Galerkin methods enriched by bubble functions cannot be constructed for quadrilateral and hexahedral elements using standard bubble functions.Comment: 25 pages, 13 figures (The previous version was compiled by mistake with the wrong style file, the current one uses amsart, and there is no difference in the text or the figures

Consistent local projection stabilized finite element methods

This work establishes a formal derivation of local projection stabilized methods as a result of an enriched Petrov-Galerkin strategy for the Stokes problem. Both velocity and pressure finite element spaces are enhanced with solutions of residual-based local problems, and then the static condensation procedure is applied to derive new methods. The approach keeps degrees of freedom unchanged while gives rise to new stable and consistent methods for continuous and discontinuous approximation spaces for the pressure. The resulting methods do not need the use of a macro-element grid structure and are parameter-free. The numerical analysis is carried out showing optimal convergence in natural norms, and moreover, two ways of rendering the velocity field locally mass conservative are proposed. Some numerics validate the theoretical results

A symmetric nodal conservative finite element method for the Darcy equation

This work introduces and analyzes novel stable Petrov-Galerkin EnrichedMethods (PGEM) for the Darcy problem based on the simplest but unstable continuous P1/P0 pair. Stability is recovered inside a Petrov-Galerkin framework where element-wise dependent residual functions, named multi-scale functions, enrich both velocity and pressure trial spaces. Unlike the velocity test space that is augmented with bubble-like functions, multi-scale functions correct edge residuals as well. The multi-scale functions turn out to be the well-known lowest order Raviart-Thomas basis functions for the velocity and discontinuous quadratics polynomial functions for the pressure. The enrichment strategy suggests the way to recover the local mass conservation property for nodal-based interpolation spaces. We prove that the method and its symmetric version are well-posed and achieve optimal error estimates in natural norms. Numerical validations confirm claimed theoretical results

Local projection finite element stabilization for the generalized Stokes problem

We analyze pressure stabilized finite element methods for the solution of the generalized Stokes problem and investigate their stability and convergence properties. An important feature of the method is that the pressure gradient unknowns can be eliminated locally thus leading to a decoupled system of equations. Although stability of the method has been established, for the homogeneous Stokes equations, the proof given here is based on the existence of a special interpolant with additional orthogonal property with respect to the projection space. This, makes it a lot simpler and more attractive. The resulting stabilized method is shown to lead to optimal rates of convergence for both velocity and pressure approximations

A two-level enriched finite element method for a mixed problem

The simplest pair of spaces is made inf-sup stable for the mixed form of the Darcy equation. The key ingredient is to enhance the finite element spaces inside a Petrov-Galerkin framework with functions satisfying element-wise local Darcy problems with right hand sides depending on the residuals over elements and edges. The enriched method is symmetric, locally mass conservative and keeps the degrees of freedom of the original interpolation spaces. First, we assume local enrichments exactly computed and we prove uniqueness and optimal error estimates in natural norms. Then, a low cost two-level finite element method is proposed to effectively obtain enhancing basis functions. The approach lays on a two-scale numerical analysis and shows that well-posedness and optimality is kept, despite the second level numerical approximation. Several numerical experiments validate the theoretical results and compares (favourably in some cases) our results with the classical Raviart-Thomas elemen