Finite elements for scalar convection-dominated equations and incompressible flow problems - A never ending story?

The contents of this paper is twofold. First, important recent results concerning finite element methods for convection-dominated problems and incompressible flow problems are described that illustrate the activities in these topics. Second, a number of, in our opinion, important problems in these fields are discussed

POD model order reduction with space-adapted snapshots for incompressible flows

We consider model order reduction based on proper orthogonal decomposition (POD) for unsteady incompressible Navier-Stokes problems, assuming that the snapshots are given by spatially adapted finite element solutions. We propose two approaches of deriving stable POD-Galerkin reduced-order models for this context. In the first approach, the pressure term and the continuity equation are eliminated by imposing a weak incompressibility constraint with respect to a pressure reference space. In the second approach, we derive an inf-sup stable velocity-pressure reduced-order model by enriching the velocity reduced space with supremizers computed on a velocity reference space. For problems with inhomogeneous Dirichlet conditions, we show how suitable lifting functions can be obtained from standard adaptive finite element computations. We provide a numerical comparison of the considered methods for a regularized lid-driven cavity problem

A posteriori error estimation and adaptivity based on VMS for the incompressible Navier–Stokes equations

In this work an explicit a posteriori error estimator for the steady incompressible Navier–Stokes equations is investigated. The error estimator is based on the variational multiscale theory, where the numerical solution is decomposed in resolved scales (FEM solution) and unresolved scales (FEM error). The error is estimated locally considering the residuals that emerge from the numerical solution and the error inverse-velocity scales, t’s, associated with each type of residual. These error scales are provided in this paper, which have been computed a-priori solving a set of local problems with unit residuals. Therefore, the computational effort to predict the error is small and its implementation in any FEM code is simple. As an application, a strategy to develop adaptive meshes with the aim of optimizing the computational effort is shown. Numerical examples are presented to test the behavior of the error estimator

Stabilized lowest order finite element approximation for linear three-field poroelasticity

A stabilized conforming mixed finite element method for the three-field (displacement, fluid flux and pressure) poroelasticity problem is developed and analyzed. We use the lowest possible approximation order, namely piecewise constant approximation for the pressure and piecewise linear continuous elements for the displacements and fluid flux. By applying a local pressure jump stabilization term to the mass conservation equation we ensure stability and avoid pressure oscillations. Importantly, the discretization leads to a symmetric linear system. For the fully discretized problem we prove existence and uniqueness, an energy estimate and an optimal a-priori error estimate, including an error estimate for the divergence of the fluid flux. Numerical experiments in 2D and 3D illustrate the convergence of the method, show the effectiveness of the method to overcome spurious pressure oscillations, and evaluate the added mass effect of the stabilization term.Comment: 25 page

Adaptive finite element simulation of incompressible flows by hybrid continuous-discontinuous Galerkin formulations

In this work we design hybrid continuous-discontinuous finite element spaces that permit discontinuities on nonmatching element interfaces of nonconforming meshes. Then we develop an equal-order stabilized finite element formulation for incompressible flows over these hybrid spaces, which combines the element interior stabilization of SUPG-type continuous Galerkin formulations and the jump stabilization of discontinuous Galerkin formulations. Optimal stability and convergence results are obtained. For the adaptive setting, we use a standard error estimator and marking strategy. Numerical experiments show the optimal accuracy of the hybrid algorithm for both uniformly and adaptively refined nonconforming meshes. The outcome of this work is a finite element formulation that can naturally be used on nonconforming meshes, as discontinuous Galerkin formulations, while keeping the much lower CPU cost of continuous Galerkin formulations

Variational Multiscale error estimators for solid mechanics adaptive simulations: an Orthogonal Subgrid Scale approach

In this work we present a general error estimator for the finite element solution of solid mechanics problems based on the Variational Multiscale method. The main idea is to consider a rich model for the subgrid scales as an error estimator. The subscales are considered to belong to a space orthogonal to the finite element space (Orthogonal Subgrid Scales) and we take into account their contribution both in the element interiors and on the element boundaries (Subscales on the Element Boundaries). A simple analysis shows that the upper bound for the obtained error estimator is sharper than in other error estimators based on the Variational Multiscale Method. Numerical examples show that the proposed error estimator is an accurate approximation for the energy norm error and can be used both in simple linear constitutive models and in more complex non-linear cases.Peer ReviewedPostprint (author's final draft

Variational Multiscale error estimators for solid mechanics adaptive simulations: an Orthogonal Subgrid Scale approach

In this work we present a general error estimator for the finite element solution of solid mechanics problems based on the Variational Multiscale method. The main idea is to consider a rich model for the subgrid scales as an error estimator. The subscales are considered to belong to a space orthogonal to the finite element space (Orthogonal Subgrid Scales) and we take into account their contribution both in the element interiors and on the element boundaries (Subscales on the Element Boundaries). A simple analysis shows that the upper bound for the obtained error estimator is sharper than in other error estimators based on the Variational Multiscale Method. Numerical examples show that the proposed error estimator is an accurate approximation for the energy norm error and can be used both in simple linear constitutive models and in more complex non-linear cases.Peer ReviewedPostprint (author's final draft

Augmented mixed finite element method for the Oseen problem: A priori and a posteriori error analyses

We propose a new augmented dual-mixed method for the Oseen problem based on the pseudostress–velocity formulation. The stabilized formulation is obtained by adding to the dual-mixed approach suitable least squares terms that arise from the constitutive and equilibrium equations. We prove that for appropriate values of the stabilization parameters, the new variational formulation and the corresponding Galerkin scheme are well-posed, and a Céa estimate holds for any finite element subspaces. We also provide the rate of convergence when each row of the pseudostress is approximated by Raviart–Thomas or Brezzi–Douglas–Marini elements and the velocity is approximated by continuous piecewise polynomials. Moreover, we derive a simple a posteriori error estimator of residual type that consists of two residual terms and prove that it is reliable and locally efficient. Finally, we include several numerical experiments that support the theoretical results.Dirección de Investigación of the Universidad Católica de la Santísima Concepción (Chile) y CONICYT-Chile FONDECYT; Ministerio de Ciencia e Innovación del Gobierno de España