Implementation of a stabilized finite element formulation for the incompressible Navier-Stokes equations based on a pressure gradient projection

We discuss in this paper some implementation aspects of a finite element formulation for the incompressible Navier-Stokes equations which allows the use of equal order velocity-pressure interpolations. The method consists in introducing the projection of the pressure gradient and adding the differen...

Descripción completa

Detalles Bibliográficos
Autores: Codina, Ramon|||0000-0002-7412-778X, Blasco Lorente, Jorge, Buscaglia, G C, Huerta, Antonio|||0000-0003-4198-3798
Tipo de recurso: artículo
Fecha de publicación:2001
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/8531
Acceso en línea:https://hdl.handle.net/2117/8531
https://dx.doi.org/10.1002/fld.182
Access Level:acceso abierto
Palabra clave:Navier-Stokes equations--Numerical solutions
finite element
incompressible flow
Navier-Stokes equations
pressure gradient projection
Equacions de Navier-Stokes -- Mètodes numèrics
Àrees temàtiques de la UPC::Matemàtiques i estadística::Anàlisi numèrica::Mètodes en elements finits
Descripción
Sumario:We discuss in this paper some implementation aspects of a finite element formulation for the incompressible Navier-Stokes equations which allows the use of equal order velocity-pressure interpolations. The method consists in introducing the projection of the pressure gradient and adding the difference between the pressure Laplacian and the divergence of this new field to the incompressibility equation, both multiplied by suitable algorithmic parameters. The main purpose of this paper is to discuss how to deal with the new variable in the implementation of the algorithm. Obviously, it could be treated as one extra unknown, either explicitly or as a condensed variable. However, we take for granted that the only way for the algorithm to be efficient is to uncouple it from the velocity-pressure calculation in one way or another. Here we discuss some iterative schemes to perform this uncoupling of the pressure gradient projection (PGP) from the calculation of the velocity and the pressure, both for the stationary and the transient Navier-Stokes equations. In the first case, the strategies analyzed refer to the interaction of the linearization loop and the iterative segregation of the PGP, whereas in the second the main dilemma concerns the explicit or implicit treatment of the PGP.