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

We consider a non-standard mixed method for the Stokes problem in R-n, n is an element of {2,3}, with Dirichlet boundary conditions, in which, after using the incompressibility condition to eliminate the pressure, the pseudostress tensor sigma and the velocity vector u become the only unknowns. Then, we apply the Babuska-Brezzi theory to prove the well-posedness of the corresponding continuous and discrete formulations. In particular, we show that Raviart-Thomas elements of order k >= 0 for s and piecewise polynomials of degree k for u ensure 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 order k >= 0 for s and continuous piecewise polynomials of degree k+1 for u become 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.