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

Title: Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves

Journal Article · · Journal of Computational Physics
 [1];  [1];  [2];  [3];  [1];  [2];  [4];  [2];  [3];  [5];  [4];  [2]
  1. Center for Turbulence Research, Stanford University, Stanford, CA 94305 (United States)
  2. Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305 (United States)
  3. Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States)
  4. Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095 (United States)
  5. NASA Ames Research Center, Moffett Field, CA 94035 (United States)

Flows in which shock waves and turbulence are present and interact dynamically occur in a wide range of applications, including inertial confinement fusion, supernovae explosion, and scramjet propulsion. Accurate simulations of such problems are challenging because of the contradictory requirements of numerical methods used to simulate turbulence, which must minimize any numerical dissipation that would otherwise overwhelm the small scales, and shock-capturing schemes, which introduce numerical dissipation to stabilize the solution. The objective of the present work is to evaluate the performance of several numerical methods capable of simultaneously handling turbulence and shock waves. A comprehensive range of high-resolution methods (WENO, hybrid WENO/central difference, artificial diffusivity, adaptive characteristic-based filter, and shock fitting) and suite of test cases (Taylor-Green vortex, Shu-Osher problem, shock-vorticity/entropy wave interaction, Noh problem, compressible isotropic turbulence) relevant to problems with shocks and turbulence are considered. The results indicate that the WENO methods provide sharp shock profiles, but overwhelm the physical dissipation. The hybrid method is minimally dissipative and leads to sharp shocks and well-resolved broadband turbulence, but relies on an appropriate shock sensor. Artificial diffusivity methods in which the artificial bulk viscosity is based on the magnitude of the strain-rate tensor resolve vortical structures well but damp dilatational modes in compressible turbulence; dilatation-based artificial bulk viscosity methods significantly improve this behavior. For well-defined shocks, the shock fitting approach yields good results.

OSTI ID:
21333930
Journal Information:
Journal of Computational Physics, Vol. 229, Issue 4; Other Information: DOI: 10.1016/j.jcp.2009.10.028; PII: S0021-9991(09)00580-4; Copyright (c) 2009 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
Country of Publication:
United States
Language:
English