数学物理学报(英文版) ›› 2010, Vol. 30 ›› Issue (2): 595-620.doi: 10.1016/S0252-9602(10)60064-1

• 论文 • 上一篇    下一篇

HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER--MESHKOV INSTABILITY TO LATE TIMES: DYNAMICS, DEPENDENCE ON INITIAL CONDITIONS, AND COMPARISONS TO EXPERIMENTAL DATA

 Oleg Schilling, Marco Latini   

  1. Lawrence Livermore National Laboratory, Livermore, California 94551, USA
    Applied and Computational Mathematics, California Institute of Technology, Pasadena, California 91125, USA
  • 收稿日期:2010-01-12 出版日期:2010-03-20 发布日期:2010-03-20
  • 基金资助:

    This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER--MESHKOV INSTABILITY TO LATE TIMES: DYNAMICS, DEPENDENCE ON INITIAL CONDITIONS, AND COMPARISONS TO EXPERIMENTAL DATA

 Oleg Schilling, Marco Latini   

  1. Lawrence Livermore National Laboratory, Livermore, California 94551, USA
    Applied and Computational Mathematics, California Institute of Technology, Pasadena, California 91125, USA
  • Received:2010-01-12 Online:2010-03-20 Published:2010-03-20
  • Supported by:

    This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

摘要:

The dynamics of the reshocked multi-mode Richtmyer--Meshkov instability is investigated using 513×2572 three-dimensional ninth-order weighted essentially nonoscillatory shock-capturing simulations. A two-mode initial perturbation with superposed random noise is used to model the Mach 1.5 air/SF6 Vetter--Sturtevant shock tube experiment. The mass fraction and enstrophy isosurfaces, and density cross-sections are utilized to show the detailed flow structure before, during, and after reshock. It is shown that the mixing layer growth agrees well with the experimentally measured growth rate before and after reshock. The post-reshock growth rate is also in good agreement with the prediction of the Mikaelian model. A parametric study of the sensitivity of the layer growth to the choice of amplitudes of the short and long wavelength initial interfacial perturbation is also presented. Finally, the amplification effects of reshock are quantified using the evolution of the turbulent kinetic energy and turbulent enstrophy spectra, as well as the evolution of the baroclinic enstrophy production, buoyancy production, and shear production terms in the enstrophy and turbulent kinetic transport equations.

关键词: Richtmyer--Meshkov instability, reshock, mixing properties, weighted essentially nonoscillatory (WENO) method

Abstract:

The dynamics of the reshocked multi-mode Richtmyer--Meshkov instability is investigated using 513×2572 three-dimensional ninth-order weighted essentially nonoscillatory shock-capturing simulations. A two-mode initial perturbation with superposed random noise is used to model the Mach 1.5 air/SF6 Vetter--Sturtevant shock tube experiment. The mass fraction and enstrophy isosurfaces, and density cross-sections are utilized to show the detailed flow structure before, during, and after reshock. It is shown that the mixing layer growth agrees well with the experimentally measured growth rate before and after reshock. The post-reshock growth rate is also in good agreement with the prediction of the Mikaelian model. A parametric study of the sensitivity of the layer growth to the choice of amplitudes of the short and long wavelength initial interfacial perturbation is also presented. Finally, the amplification effects of reshock are quantified using the evolution of the turbulent kinetic energy and turbulent enstrophy spectra, as well as the evolution of the baroclinic enstrophy production, buoyancy production, and shear production terms in the enstrophy and turbulent kinetic transport equations.

Key words: Richtmyer--Meshkov instability, reshock, mixing properties, weighted essentially nonoscillatory (WENO) method

中图分类号: 

  • 35LXX