Hostname: page-component-7b7567c8c-7pwfn Total loading time: 0 Render date: 2026-11-14T06:04:56.916Z Has data issue: false hasContentIssue false

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

Published online by Cambridge University Press:  20 August 2015

Gabriel N. Gatica*
Affiliation:
CI2MA and Departamento de Ingeniería Matemática, Universidad de Concepción, Casilla 160-C, Concepción, Chile
Antonio Márquez*
Affiliation:
CI2MA and Departamento de Ingeniería Matemática, Universidad de Concepción, Casilla 160-C, Concepción, Chile
Manuel A. Sánchez*
Affiliation:
CI2MA and Departamento de Ingeniería Matemática, Universidad de Concepción, Casilla 160-C, Concepción, Chile
*
Get access

Abstract

We 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 vector u become 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 order k≥0 for σ 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 σ 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.

Information

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Article purchase

Temporarily unavailable