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

详细信息    查看全文

• 作者：Oriol Colom&eacute ; s^a ; ^{ocolomes@cimne.upc.edu" class="auth_mail" title="E-mail the corresponding author} ; Santiago Badia^a ; ^b ; ^{sbadia@cimne.upc.edu" class="auth_mail" title="E-mail the corresponding author} ; Javier Principe^a ; ^b ; ^{principe@cimne.upc.edu" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Large eddy simulation ; Turbulence ; Variational multiscale ; Block recursive preconditioning ; Grad&ndash ; div stabilization
• 刊名：Computer Methods in Applied Mechanics and Engineering
• 出版年：2016
• 出版时间：1 June 2016
• 年：2016
• 卷：304
• 期：Complete
• 页码：294-318
• 全文大小：1762 K

文摘

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=1600$Re=1600$, and the Turbulent Channel Flow at Re_τ=395$R{e}_{\tau }=395$ and Re_τ=590$R{e}_{\tau }=590$.