Pseudostress-Based Mixed Finite Element Methods for the Stokes Problem in ℝnwith Dirichlet Boundary Conditions. I: A Priori Error Analysis

Abstract

AbstractWe consider a non-standard mixed method for the Stokes problem in ℝn,nЄ {2,3}, with Dirichlet boundary conditions, in which, after using the incompressibility condition to eliminate the pressure, the pseudostress tensor σ and the velocity vectorubecome the only unknowns. Then, we apply the Babuška-Brezzi theory to prove the well-posedness of the corresponding continuous and discrete formulations. In particular, we show that Raviart-Thomas elements of orderk≥0 for σ and piecewise polynomials of degreekforuensure unique solvability and stability of the associated Galerkin scheme. In addition, we introduce and analyze an augmented approach for our pseudostress-velocity formulation. The methodology employed is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and the Dirichlet boundary condition for the velocity, all of them multiplied by suitable stabilization parameters. We show that these parameters can be 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. For instance, Raviart-Thomas elements of orderk≥0 for σ and continuous piecewise polynomials of degreek+1 forubecome a feasible choice in this case. Finally, extensive numerical experiments illustrating the good performance of the methods and comparing them with other procedures available in the literature, are provided.