An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

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.

[1]  J. Abedi,et al.  Comparison of finite element and pendulum models for simulation of sloshing , 2003 .

[2]  M. Strelets.  Detached eddy simulation of massively separated flows , 2001 .

[3]  R. Cummings,et al.  Detached-eddy simulation with compressibility corrections applied to a supersonic axisymmetric base flow , 2002 .

[4]  M Delanaye,et al.  An accurate finite volume scheme for Euler and Navier-Stokes equations on unstructured adaptive grids , 1995 .

[5]  P. Spalart.  Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach , 1997 .

[6]  Shahrouz Aliabadi,et al.  Integrated high performance computational tools for simulations of transport and diffusion of contaminants in urban areas , 2006 .

[7]  Frans N. van de Vosse,et al.  An approximate projec-tion scheme for incompressible ow using spectral elements , 1996 .

[8]  Franck Nicoud,et al.  An approach to wall modeling in large-eddy simulations , 2000 .

[9]  Scott Percival,et al.  Viscous Drag Calculations for Ship Hull Geometry , 1995 .

[10]  Shuangzhang Tu,et al.  Development and Performance of Camel_Aero, a Truly Matrix-Free, Parallel and Vectorized Unstructured Finte Volume Solver for Compressible Flow , 2006 .

[11]  Subrata Roy,et al.  Detached-Eddy Simulation over a Reference Ahmed Car Model , 2003 .

[12]  Russell M. Cummings,et al.  DES GRID RESOLUTION ISSUES FOR VORTICAL FLOWS ON A DELTA WING AND AN F-18C , 2003 .

[13]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[14]  Shahrouz Aliabadi,et al.  Development of a hybrid finite volume/element solver for incompressible flows , 2007 .

[15]  P. Buning,et al.  Solution of the two-dimensional Euler equations with generalized coordinate transformation using flux vector splitting , 1982 .