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

doi:10.1016/S0252-9602(10)60064-1

Abstract

Abstract:

The dynamics of the reshocked multi-mode Richtmyer--Meshkov instability is investigated using 513×257² 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/SF₆ 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

CLC Number:

35LXX

Oleg Schilling, Marco Latini. HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER--MESHKOV INSTABILITY TO LATE TIMES: DYNAMICS, DEPENDENCE ON INITIAL CONDITIONS, AND COMPARISONS TO EXPERIMENTAL DATA[J].Acta mathematica scientia,Series B, 2010, 30(2): 595-620.

Trendmd

References

[1] Balsara D S, Shu C-W. Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J Comput Phys, 2000, 160: 405--452

[2] Boris J P, Grinstein F F, Oran E S, Kolbe R L. New insights into large eddy simulation. Fluid Dyn Res, 1992, 10: 199--228

[3] Brouillette M. The Richtmyer--Meshkov instability. Annu Rev Fluid Mech, 2002, 34: 445--468

[4] Cohen R H, Dannevik W P, Dimits A M, Eliason D E, Mirin A A, Zhou Y, Porter D H, Woodward P R. Three-dimensional simulation of
a Richtmyer--Meshkov instability with a two-scale initial perturbation. Phys Fluids, 2002, 14: 3692--3709}

[5] Colella P, Woodward P R. The piecewise parabolic method (PPM) for gas-dynamical simulations. J Comput Phys, 1984, 54: 174--201

[6] Collins B D, Jacobs J W. PLIF flow visualization and measurements of the Richtmyer--Meshkov instability of an air/SF₆interface. J Fluid Mech, 2002, 464: 113--136

[7] Cook R C, McEachern R L, Stephens R B. Representative surface profile power spectra from capsules used in Nova and Omega implosion experiments. Fusion Tech, 1999, 35: 224--228

[8] Drikakis D. Embedded turbulence model in numerical methods for hyperbolic conservation laws. Int J Numer Methods in Fluids, 2002, 39: 763--781

[9] Drikakis D. Advances in turbulent flow computations using high-resolution methods. Prog Aerosp Sci, 2003, 39: 405--424

[10] Drikakis D, Hahn M, Mosedale A, Thornber B. Large eddy simulation using high-resolution and high-order methods. Phil Trans Royal Soc London A, 2009, 367: 2985--2997

[11] Hill D J, Pantano C, Pullin D I. Large-eddy simulation and multiscale modelling of a Richtmyer--Meshkov instability with reshock. J Fluid Mech, 2006, 557: 29--61

[12] Jameson L. High order schemes for resolving waves: number of points per wavelength. J Sci Comput, 2000, 15: 417--439

[13] Kreiss H O, Oliger J. Comparison of accurate methods for the integration of hyperbolic equations. Tellus, 1972, 24: 199--215

[14] Latini M, Schilling O, Don W S. Effects of {WENO} flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer--Meshkov instability. J Comput Phys, 2007, 221: 805--836

[15] Latini M, Schilling O, Don W S. High-resolution simulations and modeling of reshocked single-mode Richtmyer--Meshkov instability: comparison to experimental data and to amplitude growth model predictions. Phys Fluids, 2007, 19: 024104-1--024104-19

[16] Leinov E, Sadot O, Formoza A, Malamud G, Elbaz Y, Levin L A, Ben-Dor G, Shvarts D. Investigation of the Richtmyer--Meshkov instability under re-shock conditions. Phys Scripta T, 2008, 132: 014014-1--014014-5

[17] Leinov E, Malamud G, Elbaz Y, Levin L A, Ben-Dor G, Shvarts D, Sadot O. Experimental and numerical investigation of the Richtmyer--Meshkov instability under re-shock conditions. J Fluid Mech, 2009, 626: 449--475

[18] Meloon M R. Models of Richtmyer--Meshkov instability in continuously stratified fluids
[D]. California Institute of Technology, 1998

[19] Mikaelian K O. Turbulent mixing generated by Rayleigh--Taylor and Richtmyer--Meshkov instabilities. Physica D, 1989, 36: 343--357

[20] Quirk J J, Karni S. On the dynamics of a shock--bubble interaction. J Fluid Mech, 1996, 318: 129--163

[21] Read K I. Experimental investigation of turbulent mixing by Rayleigh--Taylor instability. Physica D, 1984, 12: 45--58

[22] Roe P L. Approximate Riemann solvers, parameter vectors, and difference schemes. J Comput Phys, 1981, 43: 357--382

[23] Schilling O, Jacobs J W. Richtmyer--Meshkov instability and re-accelerated inhomogeneous flows. Scholarpedia, 2008, 3: 6090

[24] Schilling O, Latini M, Don W S. Physics of reshock and mixing in single-mode Richtmyer--Meshkov instability. Phys Rev E, 2007, 76: 026319-1--026319-28

[25] Shu C-W. High order weighted essentially nonoscillatory schemes for convection dominated problems. SIAM Rev, 2009, 51: 82--126

[26] Ukai S, Genin F, Srinivasan S, Menon S. Large eddy simulation of re-shocked Richtmyer--Meshkov instability. AIAA Paper, 2009, 2009-944

[27] Vetter M, Sturtevant B. Experiments on the Richtmyer--Meshkov instability of an air/SF₆ interface. Shock Waves, 1995, 4: 247--252

[28] von Kármán T. Progress in the statistical theory of turbulence. Proc Nat Acad Sci, 1948, 34: 530--539

[29] Youngs D L. Numerical simulation of turbulent mixing by Rayleigh--Taylor instability. Physica D, 1984, 12: 32--44

[30] Zhang Y-T, Shi J, Shu C-W, Zhou Y. Numerical viscosity and resolution of high-order weighted essentially nonoscillatory schemes for compressible flows with high Reynolds numbers. Phys Rev E, 2003, 68: 046709-1--046709-16