Mixed finite element methods with convection stabilization for the large eddy simulation of incompressible turbulent flows

The variational multiscale method thought as an implicit large eddy simulation model for turbulent flows has been shown to be an alternative to the widely used physical-based models. This method is traditionally combined with equal-order velocity–pressure pairs, since it provides pressure stabilizat...

Descripción completa

Detalles Bibliográficos
Autores: Colomés Gené, Oriol, Badia, Santiago|||0000-0003-2391-4086, Principe, Ricardo Javier|||0000-0002-1478-2651
Tipo de recurso: artículo
Fecha de publicación:2016
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/87783
Acceso en línea:https://hdl.handle.net/2117/87783
https://dx.doi.org/10.1016/j.cma.2016.02.026
Access Level:acceso abierto
Palabra clave:Turbulence--Mathematical models
Large eddy simulation
Turbulence
Variational multiscale
Block recursive preconditioning
Grad–div stabilization
Turbulència -- Models matemàtics
Àrees temàtiques de la UPC::Matemàtiques i estadística::Anàlisi numèrica::Mètodes en elements finits
Àrees temàtiques de la UPC::Física::Física de fluids
Descripción
Sumario:The variational multiscale method thought as an implicit large eddy simulation model for turbulent flows has been shown to be an alternative to the widely used physical-based models. This method is traditionally combined with equal-order velocity–pressure pairs, since it provides pressure stabilization. In this work, we consider a different approach, based on inf–sup stable elements and convection-only stabilization. In order to do so, we consider a symmetric projection stabilization of the convective term using an orthogonal subscale decomposition. The accuracy and efficiency of this method compared with residual-based algebraic subgrid scales and orthogonal subscales methods for equal-order interpolation is assessed in this paper. Moreover, when inf–sup stable elements are used, the grad–div stabilization term has been shown to be essential to guarantee accurate solutions. Hence, a study of the influence of such term in the large eddy simulation of turbulent incompressible flows is also performed. Furthermore, a recursive block preconditioning strategy has been considered for the resolution of the problem with an implicit treatment of the projection terms. Two different benchmark tests have been solved: the Taylor–Green Vortex flow with Re=1600Re=1600, and the Turbulent Channel Flow at Ret=395Ret=395 and Ret=590Ret=590.