Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer-Meshkov instability

  • Authors:

  • Marco Latini;Oleg Schilling;Wai Sun Don

  • Affiliations:

  • Applied and Computational Mathematics, California Institute of Technology, Pasadena, CA 91125, United States;University of California, Lawrence Livermore National Laboratory, Livermore, CA 94551, United States;Division of Applied Mathematics, Brown University, Providence, RI 02912, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.49

Visualization

Abstract

Weighted essentially non-oscillatory (WENO) simulations of the reshocked two-dimensional single-mode Richtmyer-Meshkov instability using third-, fifth- and ninth-order spatial flux reconstruction and uniform grid resolutions corresponding to 128, 256 and 512 points per initial perturbation wavelength are presented. The dependence of the density, vorticity, simulated density Schlieren and baroclinic production fields, mixing layer width, circulation deposition, mixing profiles, production and mixing fractions, energy spectra, statistics, probability distribution functions, numerical turbulent kinetic energy and enstrophy production/dissipation rates, numerical Reynolds numbers, and numerical viscosity on the order and resolution is investigated to long evolution times. The results are interpreted using the implicit numerical dissipation in the characteristic projection-based, finite-difference WENO method. It is shown that higher-order higher-resolution simulations have lower numerical dissipation. The sensitivity of the quantities considered to the order and resolution is further amplified following reshock, when the energy deposition by the second shock-interface interaction induces the formation of small-scale structures. Lower-order lower-resolution simulations preserve large-scale structures and flow symmetry to late times, while higher-order higher-resolution simulations exhibit fragmentation of the structures, symmetry breaking and increased mixing. Similar flow features are qualitatively and quantitatively captured by either approximately doubling the order or the resolution. Additionally, the computational scaling shows that increasing the order is more advantageous than increasing the resolution for the flow considered here. The present investigation suggests that the ninth-order WENO method is well-suited for the simulation and analysis of complex multi-scale flows and mixing generated by shock-induced hydrodynamic instabilities.