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

Une méthode hybride multi-échelles pour les problèmes de Darcy utilisant des solveurs locaux à éléments finis mixtes

International audienceMultiscaled Hybrid Mixed (MHM) method refers to a numerical technique targeted to approximate systems of differential equations with strongly varying solutions. For fluid flows, normal fluxes (multiplier) over macro element boundaries, and coarse piecewise constant potential approximations in each macro element are computed (upscaling). Then, small details are resolved by local problems, using fine representations inside the macro elements, setting the multiplier as Neumann boundary conditions (downscaling). In this work a variant of the method is developed, denoted by MHM-H(div), adopting mixed finite elements at the dowscaling stage, instead of continuous finite elements used in all previous publications of the method. Thus, this alternative MHM method inherits improvements typical of mixed methods, as better flux accuracy, and local mass conservation at the mi-cro scale level inside the macro elements, which are important properties for multi-phase flows in rough heterogeneous media. Different two-scale stable space settings are considered. Vector face functions are supposed to have normal components restricted to a given finite dimensional trace space defined over the macro element boundaries. In each macro element, the internal flux components, with vanishing normal traces, and the potential approximations, may be enriched in different extents: with respect to internal mesh size, internal polynomial degree, or both, the choice being determined by the problem at hands. A unified general error analysis of the MHM-H(div) method is presented for all these two-scale space scenarios. Both MHM versions are compared for 2D test problems, with smooth solutions, for convergence rates verification, and for a Darcy's flow in heterogeneous media. MHM-H(div) 3D simulations are presented for a known singular Darcy's solution, using adaptive macro partitions, and for an oscillatory permeability scenario

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

Enhanced error estimator based on a nearly equilibrated moving least squares recovery technique for FEM and XFEM

In this paper a new technique aimed to obtain accurate estimates of the error in energy norm using a moving least squares (MLS) recovery-based procedure is presented. We explore the capabilities of a recovery technique based on an enhanced MLS fitting, which directly provides continuous interpolated fields, to obtain estimates of the error in energy norm as an alternative to the superconvergent patch recovery (SPR). Boundary equilibrium is enforced using a nearest point approach that modifies the MLS functional. Lagrange multipliers are used to impose a nearly exact satisfaction of the internal equilibrium equation. The numerical results show the high accuracy of the proposed error estimator

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