Premium
An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow
Author(s) -
Tu Shuangzhang,
Aliabadi Shahrouz,
Patel Reena,
Watts Marvin
Publication year - 2008
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/fld.1864
Subject(s) - solver , turbulence , turbulence modeling , finite volume method , compressibility , pressure correction method , incompressible flow , finite element method , mathematics , physics , classical mechanics , mechanics , mathematical optimization , thermodynamics
In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.