Premium
An implicit mixed finite‐volume–finite‐element method for solving 3D turbulent compressible flows
Author(s) -
Hallo L.,
Le Ribault C.,
Buffat M.
Publication year - 1997
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/(sici)1097-0363(19971215)25:11<1241::aid-fld595>3.0.co;2-1
Subject(s) - finite volume method , turbulence , discretization , finite element method , mathematics , riemann solver , compressible flow , compressibility , reynolds averaged navier–stokes equations , solver , hexahedron , computational fluid dynamics , turbulence modeling , k epsilon turbulence model , mechanics , mathematical analysis , mathematical optimization , physics , thermodynamics
The development of new aeronautic projects require accurate and efficient simulations of compressible flows in complex geometries. It is well known that most flows of interest are at least locally turbulent and that the modelling of this turbulence is critical for the reliability of the computations. A turbulence closure model which is both cheap and reasonably accurate is an essential part of a compressible code. An implicit algorithm to solve the 2D and 3D compressible Navier–Stokes equations on unstructured triangular/tetrahedral grids has been extended to turbulent flows. This numerical scheme is based on second‐order finite element–finite volume discretization: the diffusive and source terms of the Navier–Stokes equations are computed using a finite element method, while the other terms are computed with a finite volume method. Finite volume cells are built around each node by means of the medians. The convective fluxes are evaluated with the approximate Riemann solver of Roe coupled with the van Albada limiter. The standard k –ϵ model has been introduced to take into account turbulence. Implicit integration schemes with efficient numerical methods (CGS, GMRES and various preconditioning techniques) have also been implemented. Our interest is to present the whole method and to demonstrate its limitations on some well‐known test cases in three‐dimensional geometries. © 1997 John Wiley & Sons, Ltd.