skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Journal Article · · Communications in Computational Physics
 [1];  [2];  [2];  [3];  [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
Citation Metrics:
Cited by: 14 works
Citation information provided by
Web of Science

References (27)

Efficient implementation of essentially non-oscillatory shock-capturing schemes, II journal July 1989
Development of low dissipative high order filter schemes for multiscale Navier–Stokes/MHD systems journal July 2007
Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers journal November 2012
A proposed modification of the Germano subgrid‐scale closure method journal March 1992
Toward the large-eddy simulation of compressible turbulent flows journal May 1992
Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters journal March 1999
Large-Eddy Simulation of the Shock/Turbulence Interaction journal July 1999
Atwood ratio dependence of Richtmyer–Meshkov flows under reshock conditions using large-eddy simulations journal February 2011
Interaction of isotropic turbulence with shock waves: effect of shock strength journal June 1997
Simulation of Richtmyer–Meshkov instability by sixth-order filter methods journal September 2007
Summation by parts, projections, and stability. I journal September 1995
High-Order Fluxes for Conservative Skew-Symmetric-like Schemes in Structured Meshes: Application to Compressible Flows journal June 2000
Entropy Splitting and Numerical Dissipation journal July 2000
Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for High Order Methods journal April 2004
Prediction of turbulent separation over a backward-facing smooth ramp journal January 2005
LES of temporally evolving mixing layers by an eighth-order filter scheme: LES OF TEMPORALLY EVOLVING MIXING LAYERS BY EIGHTH-ORDER FILTER SCHEME journal January 2012
A dynamic subgrid‐scale eddy viscosity model journal July 1991
Entropy Splitting for High-Order Numerical Simulation of Compressible Turbulence journal May 2002
Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves journal February 2010
A Subgrid-Scale Deconvolution Approach for Shock Capturing journal May 2002
ENO schemes with subcell resolution journal July 1989
High order finite difference methods with subcell resolution for advection equations with stiff source terms journal January 2012
Higher order compact implicit schemes for the wave equation journal January 1975
Comparison of three large-eddy simulations of shock-induced turbulent separation bubbles journal August 2009
Direct numerical simulation of canonical shock/turbulence interaction journal December 2009
Consistent initial conditions for the DNS of compressible turbulence journal January 1997
On the filtering paradigm for LES of flows with discontinuities journal January 2005

Cited By (1)

Skew-Symmetric Splitting and Stability of High Order Central Schemes journal May 2017