Augmented Mixed Finite Element Method Research Articles

An augmented mixed FEM for the convective Brinkman–Forchheimer problem: a priori and a posteriori error analysis

We propose and analyze an augmented mixed finite element method for the pseudostress–velocity formulation of the stationary convective Brinkman–Forchheimer problem in Rd,d∈{2,3}. Since the convective and Forchheimer terms forces the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin type terms. The resulting augmented scheme is written equivalently as a fixed point equation, so that the well-known Schauder and Banach theorems, combined with the Lax–Milgram theorem, allow to prove the unique solvability of the continuous problem. The finite element discretization involves Raviart–Thomas spaces of order k≥0 for the pseudostress tensor and continuous piecewise polynomials of degree ≤k+1 for the velocity. Stability, convergence, and a priori error estimates for the associated Galerkin scheme are obtained. In addition, we derive two reliable and efficient residual-based a posteriori error estimators for this problem on arbitrary polygonal and polyhedral regions. The reliability of the proposed estimators draws mainly upon the uniform ellipticity of the form involved, a suitable assumption on the data, a stable Helmholtz decomposition, and the local approximation properties of the Clément and Raviart–Thomas operators. In turn, inverse inequalities, the localization technique based on bubble functions, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, some numerical examples illustrating the performance of the mixed finite element method, confirming the theoretical rate of convergence and the properties of the estimators, and showing the behavior of the associated adaptive algorithms, are reported. In particular, the case of flow through a 2D porous media with fracture networks is considered.

On augmented finite element formulations for the Navier–Stokes equations with vorticity and variable viscosity

In this work, we propose and analyse an augmented mixed finite element method for solving the Navier–Stokes equations describing the motion of incompressible fluid. The model is written in terms of velocity, vorticity, and pressure, and takes into account non-constant viscosity and no-slip boundary conditions. The weak formulation of the method includes least-squares terms that arise from the constitutive equation and the incompressibility condition. We discuss the theoretical and practical implications of using augmentation in detail. Additionally, we use fixed–point strategies to show the existence and uniqueness of continuous and discrete solutions under the assumption of sufficiently small data. The method is constructed using any compatible finite element pair for velocity and pressure, as dictated by Stokes inf-sup stability, while for vorticity, any generic discrete space of arbitrary order can be used. We establish optimal a priori error estimates and provide a set of numerical tests in 2D and 3D to illustrate the behaviour of the discretisations and verify their theoretical convergence rates. Overall, this method provides an efficient and accurate solution for simulating fluid flow in a wide range of scenarios.

Velocity-vorticity-pressure formulation for the Oseen problem with variable viscosity

We propose and analyse an augmented mixed finite element method for the Oseen equations written in terms of velocity, vorticity, and pressure with non-constant viscosity and homogeneous Dirichlet boundary condition for the velocity. The weak formulation includes least-squares terms arising from the constitutive equation and from the incompressibility condition, and we show that it satisfies the hypotheses of the Babu\vska-Brezzi theory. Repeating the arguments of the continuous analysis, the stability and solvability of the discrete problem are established. The method is suited for any Stokes inf-sup stable finite element pair for velocity and pressure, while for vorticity any generic discrete space (of arbitrary order) can be used. A priori and a posteriori error estimates are derived using two specific families of discrete subspaces. Finally, we provide a set of numerical tests illustrating the behaviour of the scheme, verifying the theoretical convergence rates, and showing the performance of the adaptive algorithm guided by residual a posteriori error estimation.

New a posteriori error estimator for an stabilized mixed method applied to incompressible fluid flows

We consider an augmented mixed finite element method for incompressible fluid flows and develop a simple a posteriori error analysis. We obtain an a posteriori error estimator that is reliable and locally efficient. We provide numerical experiments that illustrate the performance of the corresponding adaptive algorithm and support its use in practice.

A priori and a posteriori error analysis of an augmented mixed-FEM for the Navier–Stokes–Brinkman problem

We introduce and analyze an augmented mixed finite element method for the Navier–Stokes–Brinkman problem with nonsolenoidal velocity. We employ a technique previously applied to the stationary Navier–Stokes equation, which consists of the introduction of a modified pseudostress tensor relating the gradient of the velocity and the pressure with the convective term, and propose an augmented pseudostress–velocity formulation for the model problem. The resulting augmented scheme is then written equivalently as a fixed point equation, so that the well-known Banach fixed point theorem, combined with the Lax–Milgram lemma, are applied to prove the unique solvability of the continuous and discrete systems. We point out that no discrete inf–sup conditions are required for the solvability analysis, and hence, in particular for the Galerkin scheme, arbitrary finite element subspaces of the respective continuous spaces can be utilized. For instance, given an integer k≥0, the Raviart–Thomas spaces of order k and continuous piecewise polynomials of degree ≤k+1 constitute feasible choices of discrete spaces for the pseudostress and the velocity, respectively, yielding optimal convergence. We also emphasize that, since the Dirichlet boundary condition becomes a natural condition, the analysis for both the continuous an discrete problems can be derived without introducing any lifting of the velocity boundary datum. In addition, we derive a reliable and efficient residual-based a posteriori error estimator for the augmented mixed method. The proof of reliability makes use of a global inf–sup condition, a Helmholtz decomposition, and local approximation properties of the Clément interpolant and Raviart–Thomas operator. On the other hand, inverse inequalities, the localization technique based on element-bubble and edge-bubble functions, approximation properties of the L2-orthogonal projector, and known results from previous works, are the main tools for proving the efficiency of the estimator. Finally, some numerical results illustrating the performance of the augmented mixed method, confirming the theoretical rate of convergence and properties of the estimator, and showing the behavior of the associated adaptive algorithms, are reported.

A posteriori error analysis of an augmented mixed method for the Navier–Stokes equations with nonlinear viscosity

In this work we develop the a posteriori error analysis of an augmented mixed finite element method for the 2D and 3D versions of the Navier–Stokes equations when the viscosity depends nonlinearly on the module of the velocity gradient. Two different reliable and efficient residual-based a posteriori error estimators for this problem on arbitrary (convex or non-convex) polygonal and polyhedral regions are derived. Our analysis of reliability of the proposed estimators draws mainly upon the global inf–sup condition satisfied by a suitable linearisation of the continuous formulation, an application of Helmholtz decomposition, and the local approximation properties of the Raviart–Thomas and Clément interpolation operators. In addition, differently from previous approaches for augmented mixed formulations, the boundedness of the Clément operator plays now an interesting role in the reliability estimate. On the other hand, inverse and discrete inequalities, and the localisation technique based on triangle-bubble and edge-bubble functions are utilised to show their efficiency. Finally, several numerical results are provided to illustrate the good performance of the augmented mixed method, to confirm the aforementioned properties of the a posteriori error estimators, and to show the behaviour of the associated adaptive algorithm.

Analysis of an augmented mixed-FEM for the Navier-Stokes problem

In this paper we propose and analyze a new augmented mixed finite element method for the Navier-Stokes problem. Our approach is based on the introduction of a “nonlinear-pseudostress” tensor linking the pseudostress tensor with the convective term, which leads to a mixed formulation with the nonlinearpseudostress tensor and the velocity as the main unknowns of the system. Further variables of interest, such as the fluid pressure, the fluid vorticity and the fluid velocity gradient, can be easily approximated as a simple postprocess of the finite element solutions with the same rate of convergence. The resulting mixed formulation is augmented by introducing Galerkin least-squares type terms arising from the constitutive and equilibrium equations of the Navier-Stokes equation and from the Dirichlet boundary condition, which are multiplied by stabilization parameters that are chosen in such a way that the resulting continuous formulation becomes well-posed. Then, the classical Banach's fixed point Theorem and Lax-Milgram's Lemma are applied to prove well-posedness of the continuous problem. Similarly, we establish wellposedness and the corresponding Cea's estimate of the associated Galerkin scheme considering any conforming finite element subspace for each unknown. In particular, the associated Galerkin scheme can be defined by employing Raviart-Thomas elements of degree k for the nonlinear-pseudostress tensor, and continuous piecewise polynomial elements of degree k+1 for the velocity, which leads to an optimal convergent scheme. In addition, we provide two iterative methods to solve the corresponding nonlinear system of equations and analyze their convergence. Finally, several numerical results illustrating the good performance of the method are provided. This contribution is based on joint work with Ricardo Oyarzua (Universidad del Bio-Bio) and Giordano Tierra (Charles University).

A fully-mixed finite element method for the Navier–Stokes/Darcy coupled problem with nonlinear viscosity

AbstractWe propose and analyze an augmented mixed finite element method for the coupling of fluid flow with porous media flow. The flows are governed by a class of nonlinear Navier–Stokes and linear Darcy equations, respectively, and the transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We apply dual-mixed formulations in both domains, and the nonlinearity involved in the Navier–Stokes region is handled by setting the strain and vorticity tensors as auxiliary unknowns. In turn, since the transmission conditions become essential, they are imposed weakly, which yields the introduction of the traces of the porous media pressure and the fluid velocity as the associated Lagrange multipliers. Furthermore, since the convective term in the fluid forces the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin type terms arising from the constitutive and equilibrium equations of the Navier–Stokes equations, and the relation defining the strain and vorticity tensors. The resulting augmented scheme is then written equivalently as a fixed point equation, so that the well-known Schauder and Banach theorems, combined with classical results on bijective monotone operators, are applied to prove the unique solvability of the continuous and discrete systems. In particular, given an integer

An augmented velocity–vorticity–pressure formulation for the Brinkman equations

SummaryThis paper deals with the analysis of a new augmented mixed finite element method in terms of vorticity, velocity, and pressure, for the Brinkman problem with nonstandard boundary conditions. The approach is based on the introduction of Galerkin least‐squares terms arising from the constitutive equation relating the aforementioned unknowns and from the incompressibility condition. We show that the resulting augmented bilinear form is continuous and elliptic, which, thanks to the Lax–Milgram theorem, and besides proving the well‐posedness of the continuous formulation, ensures the solvability and stability of the Galerkin scheme with any finite element subspace of the continuous space. In particular, Raviart–Thomas elements of any order for the velocity field, and piecewise continuous polynomials of degree k + 1 for both the vorticity and the pressure, can be utilized. A priori error estimates and the corresponding rates of convergence are also given here. Next, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem. The ellipticity of the bilinear form together with the local approximation properties of the Clément interpolation operator are the main tools for showing the reliability. In turn, inverse inequalities and the localization technique based on triangle‐bubble and edge‐bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of the method, confirming the properties of the estimators and showing the behavior of the associated adaptive algorithms, are reported. Copyright © 2015 John Wiley & Sons, Ltd.

Analysis of an augmented pseudostress-based mixed formulation for a nonlinear Brinkman model of porous media flow

Abstract In this paper we introduce and analyze an augmented mixed finite element method for the two-dimensional nonlinear Brinkman model of porous media flow with mixed boundary conditions. More precisely, we extend a previous approach for the respective linear model to the present nonlinear case, and employ a dual-mixed formulation in which the main unknowns are given by the gradient of the velocity and the pseudostress. In this way, and similarly as before, the original velocity and pressure unknowns are easily recovered through a simple postprocessing. In addition, since the Neumann boundary condition becomes essential, we impose it in a weak sense, which yields the introduction of the trace of the fluid velocity over the Neumann boundary as the associated Lagrange multiplier. We apply known results from nonlinear functional analysis to prove that the corresponding continuous and discrete schemes are well-posed. In particular, a feasible choice of finite element subspaces is given by Raviart–Thomas elements of order k ≥ 0 for the pseudostress, piecewise polynomials of degree ≤ k for the gradient of the velocity, and continuous piecewise polynomials of degree ≤ k + 1 for the Lagrange multiplier. We also derive a reliable and efficient residual-based a posteriori error estimator for this problem. Finally, several numerical results illustrating the performance and the robustness of the method, confirming the theoretical properties of the estimator, and showing the behavior of the associated adaptive algorithm, are provided.

A posteriori error analysis of an augmented mixed finite element method for Darcy flow

We develop an a posteriori error analysis of residual type of a stabilized mixed finite element method for Darcy flow. The stabilized formulation is obtained by adding to the standard dual-mixed approach suitable residual type terms arising from Darcy’s law and the mass conservation equation. We derive sufficient conditions on the stabilization parameters that guarantee that the augmented variational formulation and the corresponding Galerkin scheme are well-posed. Then, we obtain a simple a posteriori error estimator and prove that it is reliable and locally efficient. Finally, we provide several numerical experiments that illustrate the theoretical results and support the use of the corresponding adaptive algorithm in practice.

Low cost a posteriori error estimators for an augmented mixed FEM in linear elasticity

We consider an augmented mixed finite element method applied to the linear elasticity problem and derive a posteriori error estimators that are simpler and easier to implement than the ones available in the literature. In the case of homogeneous Dirichlet boundary conditions, the new a posteriori error estimator is reliable and locally efficient, whereas for non-homogeneous Dirichlet boundary conditions, we derive an a posteriori error estimator that is reliable and satisfies a quasi-efficiency bound. Numerical experiments illustrate the performance of the corresponding adaptive algorithms and support the theoretical results.

Analysis of an augmented fully-mixed approach for the coupling of quasi-Newtonian fluids and porous media

In this paper we introduce and analyze an augmented mixed finite element method for the coupling of quasi-Newtonian fluids and porous media. The flows are governed by a class of nonlinear Stokes and linear Darcy equations, respectively, and the transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We apply dual-mixed formulations in both domains, and, in order to handle the nonlinearity involved in the Stokes region, we set the strain and vorticity tensors as auxiliary unknowns. In turn, since the transmission conditions become essential, they are imposed weakly, which yields the introduction of the traces of the porous media pressure and the fluid velocity as the associated Lagrange multipliers. Moreover, in order to facilitate the analysis, we augment the formulation in the fluid by incorporating a redundant Galerkin-type term arising from the quasi-Newtonian constitutive law multiplied by a suitable stabilization parameter. In this way, under a suitable and explicit choice of this parameter, a generalization of the Babuska–Brezzi theory is utilized to show the well-posedness of the continuous and discrete formulations and to derive the corresponding a priori error estimate. In particular, the feasible finite element subspaces include PEERS and Arnold–Falk–Winther elements for the stress, velocity and vorticity in the fluid, Raviart–Thomas elements and piecewise constants for the velocity and pressure in the porous medium, together with piecewise constant Stokes strain tensor and continuous piecewise linear elements for the traces. Next, we employ classical approaches, which include linearization techniques, Clément’s interpolator and Helmholtz’s decomposition, together with known efficiency estimates, to derive a reliable and efficient residual-based a posteriori error estimator for the coupled problem. Finally, several numerical results confirming the good performance of the method and the properties of the a posteriori error estimator, and illustrating the capability of the corresponding adaptive algorithm to identify the singular regions of the solution, are reported.

An augmented mixed finite element method for the vorticity–velocity–pressure formulation of the Stokes equations

This paper deals with the numerical approximation of the stationary two-dimensional Stokes equations, formulated in terms of vorticity, velocity and pressure, with non-standard boundary conditions. Here, by introducing a Galerkin least-squares term, we end up with a stabilized variational formulation that can be recast as a twofold saddle point problem. We propose two families of mixed finite elements to solve the discrete problem, in the first family, the unknowns are approximated by piecewise continuous and quadratic elements, Brezzi–Douglas–Marini, and piecewise constant finite elements, respectively, while in the second family, the unknowns are approximated by piecewise linear and continuous, Raviart–Thomas, and piecewise constant finite elements, respectively. The wellposedness of the resulting continuous and discrete variational problems are studied employing an extension of the Babuška–Brezzi theory. We establish a priori error estimates in the natural norms, and we finally report some numerical experiments illustrating the behavior of the proposed schemes and confirming our theoretical findings on structured and unstructured meshes. Additional examples of cases not covered by our theory are also presented.

An augmented mixed finite element method for 3D linear elasticity problems

In this paper we introduce and analyze a new augmented mixed finite element method for linear elasticity problems in 3D. Our approach is an extension of a technique developed recently for plane elasticity, which is based on the introduction of consistent terms of Galerkin least-squares type. We consider non-homogeneous and homogeneous Dirichlet boundary conditions and prove that the resulting augmented variational formulations lead to strongly coercive bilinear forms. In this way, the associated Galerkin schemes become well posed for arbitrary choices of the corresponding finite element subspaces. In particular, Raviart–Thomas spaces of order 0 for the stress tensor, continuous piecewise linear elements for the displacement, and piecewise constants for the rotation can be utilized. Moreover, we show that in this case the number of unknowns behaves approximately as 9.5 times the number of elements (tetrahedrons) of the triangulation, which is cheaper, by a factor of 3, than the classical P E E R S in 3D. Several numerical results illustrating the good performance of the augmented schemes are provided.

A priori and a posteriori error analysis of an augmented mixed finite element method for incompressible fluid flows

In this paper we extend recent results on the a priori and a posteriori error analysis of an augmented mixed finite element method for the linear elasticity problem, to the case of incompressible fluid flows with symmetric stress tensor. Similarly as before, the present approach is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and from the relations defining the pressure in terms of the stress tensor and the rotation in terms of the displacement, all of them multiplied by stabilization parameters. We show that these parameters can be suitably chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well-posed for any choice of finite element subspaces. Next, we present a reliable and efficient residual-based a posteriori error estimator for the augmented mixed finite element scheme. Finally, several numerical results confirming the theoretical properties of this estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are reported.

A residual basedA POSTERIORIerror estimator for an augmented mixed finite element method in linear elasticity

In this paper we develop a residual based a posteriori error analysis for an augmented mixed finite element method applied to the problem of linear elasticity in the plane. More precisely, we derive a reliable and efficient a posteriori error estimator for the case of pure Dirichlet boundary condi- tions. In addition, several numerical experiments confirming the theoretical properties of the estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are also reported.

An augmented mixed finite element method with Lagrange multipliers: A priori and a posteriori error analyses

In this paper, we provide a priori and a posteriori error analyses of an augmented mixed finite element method with Lagrange multipliers applied to elliptic equations in divergence form with mixed boundary conditions. The augmented scheme is obtained by including the Galerkin least-squares terms arising from the constitutive and equilibrium equations. We use the classical Babuška–Brezzi theory to show that the resulting dual-mixed variational formulation and its Galerkin scheme defined with Raviart–Thomas spaces are well posed, and also to derive the corresponding a priori error estimates and rates of convergence. Then, we develop a reliable and efficient residual-based a posteriori error estimate and a reliable and quasi-efficient Ritz projection-based one, as well. Finally, several numerical results illustrating the performance of the augmented scheme and the associated adaptive algorithms are reported.

A mixed finite element method for the arch problem

We consider a general loaded arch problem. To approximate the solution of such a problem, an augmented mixed finite element method which takes into account an approximation of the middle line of the arch is considered. Some accuracy error estimates, which depend on the thickness of the arch, are then established.