A Spalart-Allmaras turbulence model implementation in a discontinuous Galerkin solver for incompressible flows

In this paper the artificial compressibility flux Discontinuous Galerkin (DG) method for the solution of the incompressible Navier-Stokes equations has been extended to deal with the Reynolds-Averaged Navier-Stokes (RANS) equations coupled with the Spalart-Allmaras (SA) turbulence model. DG implementations of the RANS and SA equations for compressible flows have already been reported in the literature, including the description of limiting or stabilization techniques adopted in order to prevent the turbulent viscosity @[email protected]? from becoming negative. In this paper we introduce an SA model implementation that deals with negative @[email protected]? values by modifying the source and diffusion terms in the SA model equation only when the working variable or one of the model closure functions become negative. This results in an efficient high-order implementation where either stabilization terms or even additional equations are avoided. We remark that the proposed implementation is not DG specific and it is well suited for any numerical discretization of the RANS-SA governing equations. The reliability, robustness and accuracy of the proposed implementation have been assessed by computing several high Reynolds number turbulent test cases: the flow over a flat plate (Re=10^7), the flow past a backward-facing step (Re=37400) and the flow around a NACA 0012 airfoil at different angles of attack (@a=0^o,10^o,15^o) and Reynolds numbers (Re=2.88x10^6,6x10^6).

[1]  Marco Luciano Savini,et al.  Discontinuous Galerkin solution of the Reynolds-averaged Navier–Stokes and k–ω turbulence model equations , 2005 .

[2]  H. L. Seegmiller,et al.  Features of a reattaching turbulent shear layer in divergent channel flow , 1985 .

[3]  Guido Kanschat,et al.  Local Discontinuous Galerkin Methods for the Stokes System , 2002, SIAM J. Numer. Anal..

[4]  Andrea Crivellini,et al.  An implicit matrix-free Discontinuous Galerkin solver for viscous and turbulent aerodynamic simulations , 2011 .

[5]  Andrea Crivellini,et al.  An artificial compressibility flux for the discontinuous Galerkin solution of the incompressible Navier-Stokes equations , 2006, J. Comput. Phys..

[6]  C. Ross Ethier,et al.  A high-order discontinuous Galerkin method for the unsteady incompressible Navier-Stokes equations , 2007, J. Comput. Phys..

[7]  K Wieghardt,et al.  On the turbulent friction layer for rising pressure , 1951 .

[8]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[9]  David L. Darmofal,et al.  An unsteady adaptation algorithm for discontinuous Galerkin discretizations of the RANS equations , 2007 .

[10]  Ronald Cools,et al.  An encyclopaedia of cubature formulas , 2003, J. Complex..

[11]  Andrea Crivellini,et al.  A high-order discontinuous Galerkin method for natural convection problems , 2006 .

[12]  Daniele A. Di Pietro,et al.  A pressure-correction scheme for convection-dominated incompressible flows with discontinuous velocity and continuous pressure , 2011, J. Comput. Phys..

[13]  C. Rumsey.  Apparent Transition Behavior of Widely-Used Turbulence Models , 2006 .

[14]  Steven Allmaras,et al.  Multigrid for the 2-D compressible Navier-Stokes equations , 1999 .

[15]  Tuncer Cebeci,et al.  Analysis of turbulent flows , 2004 .

[16]  Guido Kanschat,et al.  A locally conservative LDG method for the incompressible Navier-Stokes equations , 2004, Math. Comput..

[17]  Nicholas J. Georgiadis,et al.  Implementation and Validation of the Chien k-epsilon Turbulence Model in the Wind Navier-Stokes Code , 1999 .

[18]  Ralf Hartmann,et al.  Discontinuous Galerkin methods for computational aerodynamics — 3D adaptive flow simulation with the DLR PADGE code , 2010 .

[19]  Per-Olof Persson,et al.  The Compact Discontinuous Galerkin (CDG) Method for Elliptic Problems , 2007, SIAM J. Sci. Comput..

[20]  Brian R. Smith,et al.  Description of a Website Resource for Turbulence Modeling Verification and Validation , 2010 .

[21]  Per-Olof Persson,et al.  RANS Solutions Using High Order Discontinuous Galerkin Methods , 2007 .

[22]  S. Rebay,et al.  An implicit high-order discontinuous Galerkin method for steady and unsteady incompressible flows , 2007 .

[23]  Dominik Obrist,et al.  High-order accurate solution of the incompressible Navier-Stokes equations on massively parallel computers , 2010, Journal of Computational Physics.

[24]  D. Wilcox.  Turbulence modeling for CFD , 1993 .

[25]  William G. Johnson,et al.  Pressure distributions from high Reynolds number transonic tests of an NACA 0012 airfoil in the Langley 0.3-meter transonic cryogenic tunnel , 1987 .

[26]  Ewald Krämer,et al.  A parallel, high-order discontinuous Galerkin code for laminar and turbulent flows , 2008 .

[27]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[28]  Chi-Wang Shu,et al.  A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows , 2000 .

[29]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[30]  Todd A. Oliver.  A High-Order, Adaptive, Discontinuous Galerkin Finite Element Method for the Reynolds-Averaged Navier-Stokes Equations , 2008 .

[31]  D. E. Coles,et al.  PROCEEDINGS: COMPUTATION OF TURBULENT BOUNDARY LAYERS - 1968 AFOSR-IFP-STANFORD CONFERENCE, STANFORD UNIV., CALIF., 19-24 AUGUST 1968. VOLUME II. COMPILED DATA, , 1969 .

[32]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[33]  Bernardo Cockburn,et al.  An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier-Stokes equations , 2011, J. Comput. Phys..

[34]  Guido Kanschat,et al.  The local discontinuous Galerkin method for linearized incompressible fluid flow: a review , 2005 .