A term by term stabilization algorithm for finite element solution of incompressible flow problems

This paper introduces a stabilization technique for Finite Element numerical solution of 2D and 3D incompressible flow problems. It may be applied to stabilize the discretization of the pressure gradient, and also of any individual operator term such as the convection, curl or divergence operators, with specific levels of numerical diffusion for each one of them. Its computational complexity is reduced with respect to usual (residual-based) stabilization techniques. We consider piecewise affine Finite Elements, for which we obtain optimal error bounds for steady Navier-Stokes and also for generalized Stokes equations (including convection). We include some numerical experiment in well known 2D test cases, that show its good performances.

Authors and Affiliations

1. Departamento de Ecuaciones Diferenciales y Análisis Numérico, Universidad de Sevilla. C/ Tarfia, s/n. E-41080 Sevilla, Spain; e-mail: chacon@numer.us.es , , , , , , ES

   Tomás Chacón Rebollo

Received March 15, 1996 / Revised version received January 17, 1997

Reprints and permissions