Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

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; 933337
Journal Information:
Communications in Computational Physics, Journal Name: Communications in Computational Physics Journal Issue: 02 Vol. 19; ISSN 1815-2406; ISSN applab
Publisher:
Global Science PressCopyright Statement
Country of Publication:
United States
Language:
English

References (27)

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
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
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
Entropy Splitting for High-Order Numerical Simulation of Compressible Turbulence journal May 2002
A Subgrid-Scale Deconvolution Approach for Shock Capturing journal May 2002
Simulation of Richtmyer–Meshkov instability by sixth-order filter methods journal September 2007
Comparison of three large-eddy simulations of shock-induced turbulent separation bubbles journal August 2009
Efficient implementation of essentially non-oscillatory shock-capturing schemes, II journal July 1989
ENO schemes with subcell resolution journal July 1989
Development of low dissipative high order filter schemes for multiscale Navier–Stokes/MHD systems journal July 2007
Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves journal February 2010
High order finite difference methods with subcell resolution for advection equations with stiff source terms journal January 2012
Atwood ratio dependence of Richtmyer–Meshkov flows under reshock conditions using large-eddy simulations journal February 2011
Toward the large-eddy simulation of compressible turbulent flows journal May 1992
Interaction of isotropic turbulence with shock waves: effect of shock strength journal June 1997
Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for High Order Methods journal April 2004
Direct numerical simulation of canonical shock/turbulence interaction journal December 2009
A dynamic subgrid‐scale eddy viscosity model journal July 1991
A proposed modification of the Germano subgrid‐scale closure method journal March 1992
Consistent initial conditions for the DNS of compressible turbulence journal January 1997
Prediction of turbulent separation over a backward-facing smooth ramp journal January 2005
On the filtering paradigm for LES of flows with discontinuities journal January 2005
Higher order compact implicit schemes for the wave equation journal January 1975
Summation by parts, projections, and stability. I journal September 1995
Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers journal November 2012

Cited By (1)

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

Similar Records

Construction of low dissipative high-order well-balanced filter schemes for non-equilibrium flows
Journal Article · Fri May 20 00:00:00 EDT 2011 · Journal of Computational Physics · OSTI ID:21499746

Eulerian-Eulerian large eddy simulation of two-phase dilute bubbly flows
Journal Article · Thu Aug 08 20:00:00 EDT 2019 · Chemical Engineering Science · OSTI ID:1799155

Three dimensional simulations of Richtmyer-Meshkov instabilities in gas-curtain shock-tube experiments
Conference · Wed Dec 31 23:00:00 EST 2008 · OSTI ID:989837