High Order Numerical Methods for the Dynamic SGS Model of Turbulent Flows with Shocks

Journal Article · 01 February 2016 · Communications in Computational Physics

DOI: https://doi.org/10.4208/cicp.211014.040915a · OSTI ID:1477827

Kotov, D. V. ^[1]; Yee, H. C. ^[2]; Wray, A. A. ^[2]; Hadjadj, A. ^[3]; Sjögreen, B. ^[4]

---

1. Bay Area Environmental Research Institute, Petaluma, CA (United States)
2. NASA-Ames Research Center, Moffett Field, CA (United States)
3. CORIA UMR 6614 & INSA de Rouen, St-Etienne du Rouvray (France)
4. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

Simulation of turbulent flows with shocks employing subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. This paper addresses the accuracy improvement of LES of turbulent flows in two ways: (a) from the SGS model standpoint and (b) from the numerical method improvement standpoint. In an internal report, Kotov et al. ( “High Order Numerical Methods for large eddy simulation (LES) of Turbulent Flows with Shocks”, CTR Tech Brief, Oct. 2014, Stanford University), we performed a preliminary comparative study of different approaches to reduce the loss of accuracy within the framework of the dynamic Germano SGS model. The high order low dissipative method of Yee & Sjögreen (2009) using local flow sensors to control the amount of numerical dissipation where needed is used for the LES simulation. The considered improved dynamics model approaches include applying the one-sided SGS test filter of Sagaut & Germano (2005) and/or disabling the SGS terms at the shock location. For Mach 1.5 and 3 canonical shock-turbulence interaction problems, both of these approaches show a similar accuracy improvement to that of the full use of the SGS terms. The present study focuses on a five levels of grid refinement study to obtain the reference direct numerical simulation (DNS) solution for additional LES SGS comparison and approaches. As a result, one of the numerical accuracy improvements included here applies Harten's subcell resolution procedure to locate and sharpen the shock, and uses a one-sided test filter at the grid points adjacent to the exact shock location.

• Research Organization: Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
• Sponsoring Organization: USDOE National Nuclear Security Administration (NNSA)
• Grant/Contract Number: AC52-07NA27344
• OSTI ID: 1477827
• Report Number(s): LLNL-JRNL-748257; applab; 933337
• Journal Information: Communications in Computational Physics, Vol. 19, Issue 02; ISSN 1815-2406
• Publisher: Global Science PressCopyright Statement
• Country of Publication: United States
• Language: English