Initial conditions of perturbation for single‐mode RMI.
Abstract
Richtmyer‐Meshkov instability and turbulent mixing are fundamental problems of multi‐materials interface dynamics, which mainly focuses on the growth of perturbation on the interface and mixing of different materials. It is very important in many applications such as inertial confinement fusion, high‐speed combustion, supernova, etc. In this chapter, we will gain advances in understanding this problem by numerical investigations, including the numerical method and program we used, the verification and validation of numerical method and program, the growth laws and mechanics of turbulent mixing, the effects of initial conditions, the dynamic behavior, and some new phenomenon for Richtmyer‐Meshkov instability and turbulent mixing.
Keywords
- Richtmyer‐Meshkov instability
- turbulent mixing
- interface dynamics
- dynamic behavior
1. Introduction
When a shock wave accelerates a perturbed interface separating two different fluids, the Richtmyer‐Meshkov (RM) instability will occur. This phenomenon was theoretically analyzed by Richtmyer [1] for the first time in 1960, which was confirmed in experiment by Meshkov [2] in 1969. The main mechanism of the Richtmyer‐Meshkov instability (RMI) is the baroclinic vorticity deposition at the interface due to the misalignment of the pressure gradient across the shock and the local density gradient at the interface (∇
The Richtmyer‐Meshkov instability and induced turbulent mixing are very important in a variety of man‐made applications and natural phenomena such as inertial confinement fusion (ICF) [4], deflagration‐to‐detonation transition (DDT) [5], high‐speed combustion [6], and astrophysics (i.e., supernova explosions) [7]. For ICF, the ablative shell that encapsulates the deuterium‐tritium fuel becomes RM unstable as it is accelerated inward by the ablation of its outer surface by laser or secondary X‐ray radiation. The degree of compression achievable in laser fusion experiments is ultimately limited by Richtmyer‐Meshkov and Rayleigh‐Taylor instabilities. Thus, these instabilities represent the most significant barriers to attaining positive‐net‐yield fusion reactions in laser fusion facilities. For the supersonic combustion, the Richtmyer‐Meshkov instability caused by the interaction of a shock wave with a flame front can greatly promote the mixing of fuel and oxidant and enhance the burning rate. For the supernova explosions, the Richtmyer‐Meshkov instability was believed to occur when the outward propagating shock wave generated by the collapsing core of a dying star passes over the helium‐hydrogen interface. Observations of the optical output of the supernova 1987A suggest that the outer regions of the supernova were much more uniformly mixed than expected, and indicating significant Richtmyer‐Meshkov mixing had occurred [8]. The Richtmyer‐Meshkov instability and turbulent mixing are so important and have gained significant attention. However, the turbulent mixing is a complicated three‐dimensional unstable flow, which spans a wide range of time‐space scales.
2. Numerical method and program
By applying the piecewise parabolic method (PPM), volume of fluid (VOF), parallel technique, and so on, we have developed a series of Euler compressible multi‐fluid dynamic programs with three orders precision, such as MFPPP [9], MVPPM [10], and multi‐viscous‐flow and turbulence (MVFT) [11, 12]. MFPPM program is not considering the fluid viscosity, which only solves the Euler equations. MVPPM program is considering the molecular viscosity but it is not changed with temperature. Multi‐viscous‐flow and turbulence (MVFT) is a large‐eddy simulation (LES) program that has four choices of subgrid‐scale (SGS) stress models including the Smagorinsky model [13], Vreman model (VM) [14], dynamic Smagorinsky model (DSM) [15], and stretched‐vortex model (SVM) [16].
2.1. Governing equations
The governing equations of large‐eddy simulation are the Favre‐filtered compressible multi‐viscous‐flow Navier‐Stokes (NS) equations and are written as in tensor form:
here
where
For a multi‐material mixture case, the average quantities and physical properties of a mixed phase are supposed as a volume‐weighted sum over the set of components
For large‐eddy simulation, in the filtering operation, the unresolved‐scale motions identified as the “subgrid‐scale” are filtered, but the effects of subgrid‐scale motions on the resolved‐scale motions are retained in the governing equations in the form of SGS turbulence transport fluxes, which must be modeled to complete the closure of LES equations. The SGS stress tensor, the heat, and scalar transport flux are defined as
2.2. Subgrid‐scale stress models for LES
2.2.1. Smagorinsky model
The SGS turbulence behavior is assumed to be analogous to the molecular dissipative mechanism, so that the molecular viscosity and diffusion models can be used to simulate the SGS fluxes, and the SGS stress tensor, the heat, and the scalar transport flux are calculated as follows [13]
where the SGS turbulent viscosity
where the dimensionless coefficient
2.2.2. Vreman model
The Vreman SGS model is also an eddy viscosity model and is as follows [14]:
with
The model constant
2.2.3. Dynamic Smagorinsky model
In general, the turbulent kinetic energy is transferred from large scales to small scales in turbulent flows, which is called forward scatter of energy, and then it is dissipated by the viscous action. However, the reversed energy flow from small scales to large scales (the process called backscatter) may also occur. Although the backscatter is just a small range of local phenomenon, it has been shown to be of importance, especially in the transition regime. The physical mechanism of most of the widely used SGS models is the forward scatter of energy; in other words, the SGS models are absolutely dissipative, such as the eddy viscosity models (Smagorinsky and Vreman SGS models). In order to account for the backscatter of energy, several modifications to the eddy viscosity models have been proposed. An improvement is to calculate the eddy viscosity coefficient dynamically which is a function of space and time and can be negative in some regions [15].
The dynamic Smagorinsky model for compressible turbulence is as follows, the model coefficient is
where
the overbars of “¯” and “^” denote the grid filter and test filter, respectively.
2.2.4. Stretched‐vortex model
The stretched‐vortex model is based on an explicit structural modeling of small‐scale dynamics. It can simulate the multiscale compressible turbulence and allows the anisotropy of the flow to be extended to the dissipation scale. In the stretched‐vortex model, the flow within a computational grid cell is assumed to result from an ensemble of straight, nearly axisymmetric vortices aligned with the local resolved scale strain or vorticity. The resulting SGS flux terms are [16]
where
2.3. Numerical algorithm
For the convenience, in the following sections, the overbar and overtilde of variables are omitted. Operator splitting technique is used to solve the physical problems, described by Eq. (1), into three sub‐processes, i.e., the computation of inviscid flux, viscous flux, and heat flux. Eq. (1) can be split into two equations as follows:
and
For the inviscid flux, the three‐dimensional problem can be simplified into one‐dimensional (1D) problem in three directions of
And then the one‐dimensional Eq. (27) in each direction is resolved by the two‐step Lagrange‐Remapping algorithm. Also, one time step is divided into four substeps: (i) the piecewise parabolic interpolating of physical quantities, (ii) solving the Riemann problems approximately, (iii) marching of the Lagrange equations, and (iv) remapping the physical quantities back to the stationary Euler meshes. The orders of accuracy of the spatial and temporal schemes are the third and second, respectively, for smooth flows.
2.3.1. Lagrange step of finite volume method
The Lagrange equations in 1D for multi‐fluid can be written as
where
Because of the calculations are referred to as scalars, here the averaged physical quantities in a cell are written as a uniformity
2.3.2. Remap step
In Lagrange step, the computational cells distort to follow the material motion, and in Remap step, the averaged physical quantities at the distorted Lagrange cells are remapped back to the stationary Euler meshes. The piecewise parabolic interpolation and integral methods in this step are the same as the ones in Lagrange step.
After Remap step, we define the physical quantities in Euler mesh as
where
2.3.3. Viscous flux, heat flux, and scalar flux
The viscous flux, heat flux, and scalar flux of Eq. (26) are calculated based on the computed inviscid flux by using second‐order spatial center difference and two‐step Runge‐Kutta time marching. The first equation of Eq. (26) can be neglected, and which is written in conserved form
and
In the Cartesian coordinate system, the spatial derivative terms of Eq. (40) can be dispersed as follows:
Then, Eq. (40) is solved by two Runge‐Kutta time marching,
The detailed descriptions of numerical algorithm are referred to Ref. [17].
3. Verification and validation
In this section, the validity and reliability of our compressible multi‐fluid dynamic programs are to be verified and validated by comparisons with analytical solutions and hydrodynamic interface instability experiments in a shock tube.
3.1. Riemann problem of condensed matter
A copper pellet collides with a copper target with three velocities of 2, 4, and 8 km/s. Using the reduction form of Gruneisen EOS for copper in simulations, the one‐dimensional numerical results of postshock density (a), pressure (b), and velocity (c) compare with theoretical solutions, as shown in Figure 2, the solid lines corresponding to numerical results and the dot lines corresponding to the theoretical solutions.
3.2. Riemann problem of one‐dimensional gas/liquid
At initial time, the region [0, 1.0 cm] is filled with gas with high pressure 1.0 × 108 Pa and density 1.29 g/cm3. The region [1.0 cm, 5.0 cm] is filled with water with pressure 1.0 × 105 Pa and density 1.0 g/cm3. The gas and water are all described by Stiffen gas EOS. The left and right boundaries are flow. Figure 3 shows the distributions of the density (a), pressure (b), and velocity (c) at 20 μs. There are a forward shock wave and a backward rarefaction wave after interaction. The pressure and velocity around the interface are well continuous and have no nonphysical oscillation.
3.3. Single‐mode Richtmyer‐Meshkov instability
The two‐ and three‐dimensional single‐mode Richtmyer‐Meshkov instabilities are numerical simulated by MVPPM program, which also compare with the theoretical model [18]. The initial small perturbation is a sinusoidal one with wavelength 60 mm (global wavelength for a three‐dimensional case) and amplitude 1.0 mm. The incident air shock wave with Mach number 1.2 impacts the air/SF6 single‐mode interface. Figure 4 shows the comparisons of amplitude of single‐mode perturbed interface with the linear impulsive and nonlinear models, the left one corresponding to the two‐dimensional (2D) results, and the right one corresponding to the three‐dimensional results.
3.4. Simulations of shock tube experiments
The compressible multi‐fluid dynamic programs are used to simulated several shock tube experiments of interface instability for validation. These experiments include planar interface and gas cylinder shock tube experiment and planar and cylindrical jelly experiments.
Figure 5 shows the comparison of 2D calculated width of turbulent mixing zone (TMZ) [19] of Leinov’s shock tube experiment with reshock with experiment [20] in which the Mach number of incident air shock wave is 1.2, which impacts the air/SF6 interface, the grid size is 50 μm. Figure 6 shows the 3D calculated TMZ width [21] of Poggi’s multi‐mode shock tube experiment [22] with reshock in which the Mach number of incident SF6 shock wave is 1.453, which impacts the SF6/air interface. Figures 7 [10] and 8 [23] show the calculated and experimental interface images of AWE’s SF6 half‐height and double‐bump shock tube experiment [24, 25], the Mach number of incident air shock wave is 1.26, respectively. Figure 9 [26] shows the SF6 gas cylinder evolution at different times under the initial air shock wave with the Mach number 1.2, Figure 10 shows the width and height of gas cylinder at different times for experiment and numerical simulations.
For the problems of interface instability with a high density ratio, our hydrodynamic programs are also applicable. For the LLNL’s jelly ring experiments [27], the jelly ring only has a sinusoidal periodic initial perturbation at the outer interface and is driven by expansion of the explosion products of a gaseous mixture of C2H2 and O2. The jelly is mainly made of water. The thickness of ring is 15 mm, the mode number of perturbation is 6 and 36, respectively, and the initial amplitude is 1.0 mm. Figure 11 shows the evolutions of the jelly ring at different times from simulated from LLNL’ CALE program (a) and our MVPPM program (b), and the experimental image (c) at 600 μs [28]. Figure 12 shows the time histories of radius (a) and the amplitude (b) of the outer and inner interfaces of the jelly ring simulated from CALE and MVPPM programs [28]. Figure 13 shows the images of jelly ring at 600 μs for mode number 36 from experiment (a), CALE (b), and MVPPM (c) programs [28].
From the above comparisons between our numerical simulations and theoretical model, experiments for the problems of hydrodynamic interface instability from the density ratio low to the high density ratio, the numerical results agree well with theory and experiments, so the validity and reliability of our compressible multi‐fluid dynamic programs have been verified and validated.
4. Interface instability and turbulent mixing
As we know, the physical mechanism for the occurrence of Richtmyer‐Meshkov instability is the baroclinic vorticity deposition at the interface resulting from the misalignment of the pressure gradient across the shock front and the local density gradient across the interface. The evolution equation of vorticity is as follows:
where
In this section, we will introduce the growth laws of the Richtmyer‐Meshkov instability and the induced turbulent mixing and its dynamic behavior by numerical simulations.
4.1. Growth laws of Richtmyer‐Meshkov instability and turbulent mixing
Figure 14 shows the amplitude (a) and growth rate (b) of two‐dimensional single‐mode Richtmyer‐Meshkov instability without reshock for initial perturbation amplitude 1 and 12 mm and wavelength 60 mm. After initial shock, the perturbation enters the nonlinear stage quickly. The growth rate increases fast and reaches the highest peak, then will reduce owing to the reduction of the effect of compressibility and the dominant role of flow nonlinearity. At the late times, the amplitude is increasing linearly, and the growth rate remains a constant. For the larger initial amplitude, the perturbation enters the later linear growth stage earlier.
For the multi‐mode Richtmyer‐Meshkov instability and the induced turbulent mixing with reshock, the initial air shock wave with Mach number 1.2 impacts the air/SF6 interface. The multi‐mode initial perturbation is composed of eight dominant mode wavelengths of 0.8, 1.0, 1.25, 1.6, 2.0, 2.5, 3.2, and 4.0 mm superimposed with a very small random disturbance. The shock tube experiment can be referred to Ref. [29]. Figure 15 shows the schematic of shock tube and computational model. The transmitted wave rebounds between the interface and the end‐wall of shock tube and produces multiple shock‐interface interactions. At about 1.7 ms, the transmitted shock wave reflected back off the end‐wall impacts the interface. Figure 16 shows the time history of TMZ width, the black line denotes the numerically calculated width and the red line is obtained from fitting of the numerical results. As can be seen, after the initial shock, the TMZ width starts to grow as a power law ;
4.2. Evaluation of different subgrid‐scale stress models
In large‐eddy simulation, the effect of small scales on large‐scale motions is represented by the SGS stress model. Most of the commonly used SGS models assume that the main function of subgrid scales is to remove energy from the large scales and dissipate it through the action of the viscous forces. But, as we know, in fact the energy is also transferred from the small scales to the large scales (backscatter) in a small and local range. The SGS turbulent dissipation, which is the work of SGS stress, represents the energy transfer between resolved and subgrid scales,
If it is negative, the subgrid scales remove energy from the resolved scales (forward scatter); if it is positive, they release energy to the resolved scales (backscatter). It is easy to see that the eddy viscosity models such as Smagorinsky model, Vreman model, etc. are absolutely dissipative.
Figure 18 shows the distribution of calculated SGS turbulent dissipation in streamwise direction for Smagorinsky and Vreman models and for AWE’s SF6 half‐height shock tube experiment [24] at two times [11]. The SGS dissipation of the Smagorinsky model is much greater than the Vreman model over 1.5 times; therefore, the dissipation is too great for the Smagorinsky SGS model. The SGS turbulent dissipations of the Vreman model (VM), the dynamic Smagorinsky model (DSM), and the stretched‐vortex model (SVM) based on a planar Richtmyer‐Meshkov instability with incident Mach number 1.2 are shown in Figure 19 [30]. Before the interfacial flow has developed to be turbulent completely, the dynamic and stretched‐vortex models have all predicted the energy backscatter, but the energy backscatter predicted by the dynamic model is larger and its range is wider. After reshock, the turbulent fluctuations are stronger extremely the turbulent dissipation also increases. The dynamic model’s dissipation is the highest, then followed by the Vreman model, and the stretched‐vortex model’s dissipation is the lowest. At the late time, the SGS dissipation of the dynamic and Vreman models is the same basically, and the dynamic model is also to be absolutely dissipative, yet the stretched‐vortex model is still able to predict the local phenomenon of energy backscatter in a small range. So, the dynamic model is poor in representing the energy backscatter. The Vreman and stretched‐vortex models are all robust, but the former is absolutely dissipative.
4.3. Effects of the initial conditions on the growth of RMI and the turbulent mixing
4.3.1. Effects of the initial conditions on the growth of single‐mode RMI
First, we consider the effects of initial conditions of perturbation amplitude and wavelength on the growth of single‐mode Richtmyer‐Meshkov instability without reshock [18]. The incident shock wave with Mach number 1.2 impacts the air/SF6 interface. The initial perturbation amplitude and wavelength are listed in Table 1. The calculations are carried out in two dimensions, and the RMI does not develop into turbulent mixing completely. Figure 20 shows the effects of the initial perturbation amplitude on the growth of single‐mode RMI for the fixed initial perturbation wavelength
Wavelength (mm) | 60 | 30 | 40 | 50 | 70 | 80 | |||||||||
Amplitude (mm) | 1 | 2 | 3 | 4 | 5 | 6 | 8 | 10 | 12 | 14 | 3 |
4.3.2. Effects of the initial conditions on the growth of multi‐mode RMI and the induced turbulent mixing
For the effects of three‐dimensional initial multi‐mode conditions, the case of Richtmyer‐Meshkov instability is same as the above Section 4.1. Table 2 lists the initial condition of perturbations, the perturbation strength (
CASE1 | CASE2 | CASE3 | CASE4 | CASE5 | |
---|---|---|---|---|---|
0.07 | 0.14 | 0.28 | 0.56 | 1.12 | |
0.035 | .0.7 | 0.14 | 0.28 | 0.56 |
where
Figure 22 shows the growth history of TMZ width. Figures 23 and 24 show the time histories of the perk values of turbulent kinetic energy and enstrophy, respectively. For the larger initial perturbation strength, the TMZ grows faster, the turbulent kinetic energy is also larger or the turbulence strength is also stronger, the deposited vorticity is larger too. The development of turbulent mixing has a strong dependence on the initial conditions between the initial shock and the impingement of the first reflected rarefaction wave, after that the evolution of the turbulent mixing has lost the memory of the initial conditions.
4.4. Dynamic characters of multi‐mode RM instability and induced turbulent mixing
Figures 25–27 show the spatial profiles of the turbulent kinetic energy, dissipation rate and enstrophy along the motion direction of shock wave at different times before and after reshock individually. They all have a spatial profile similar to Gaussian distribution. The strongest turbulent intensity is located in the center of TMZ. Figures 28–30 show the temporal evolutions of the peak values of the turbulent kinetic energy, dissipation rate, and enstrophy in the spatial profiles, along with their fitted results, respectively. The turbulent kinetic energy, dissipation rate, and enstrophy decay gradually because of dissipation and diffusion. After the initial shock and before the reshock, the turbulent kinetic energy and enstrophy decay with time as a power law ;
4.5. Numerical study of the elliptic gas cylinder instability
Our group first performed the experimental and numerical investigations of the elliptic gas cylinder instability. As we know that the initial density distribution of the gas cylinder is hard to determine in the experiment, which can only give the one‐dimensional radial concentration distribution for circular gas cylinder as an approximate Gaussian function. Drawing to the experience on the one‐dimensional diffusive interfacial transition layer with finite thickness for circular gas cylinder, we constructed a two‐dimensional diffusive interfacial transition layer with finite thickness for elliptic gas cylinder through numerical simulation [32],
where
Figure 34 shows the vorticity at the center of the core and the distance between the two vortex cores of these two simulations. When the incident shock accelerates the elliptic gas cylinder along the major axis, the absolute vorticity in the vortex core |
4.6. Foundation of new RM instability and mechanism
At present, the investigation for the Richtmyer‐Meshkov instability is performed in a uniform flow field. We first study the Richtmyer‐Meshkov instability and turbulent mixing in a nonuniform flow field by experiment and numerical simulations and find some new phenomena. Figure 35 shows the test section schematic of shock tube. The air shock wave with Mach number 1.27 accelerates the two‐mode sinusoidal air/SF6 interface (amplitude
Figure 38 shows the evolution of the interface and the propagation of the transmitted shock wave, the calculated results agree well with the experiment. Due to the nonuniform flow field of SF6 gas, the propagating velocity of transmitted shock wave in the upper part of shock tube is faster than in the bottom of shock tube, and it forms an oblique shock wave front and the interface. Figure 39 shows the shock front locations between experiment and numerical simulations, Figure 40 shows the perturbation amplitude of experiment, numerical simulation, and theories, they are in good agreement. Figure 41 shows the calculated perturbation amplitude of RM instability for different modes in the initial uniform and nonuniform flow fields. As we can see that, at late times, the growth of small perturbation in a low‐density zone catches up and exceeds the large perturbation in a high‐density zone, which is opposite to the case of uniform flow field. Figure 42 shows the perturbation amplitudes of four different initial amplitudes 2.5, 5.0, 7.5, and 10 mm in low‐ and high‐density zones of nonuniform flow fields. It shows that the effect of initial amplitude on the growth of RM instability in the nonuniform flow field is different from the case of uniform one, which can be explained by the baroclinic vorticity shown in Figures 43 and 44, the baroclinic vorticity produced in the low‐density zone is larger than that in the high‐density zone.
For the above investigations of RM instability in the nonuniform flow field, they are all without reshock. In the following section, we will study the effect of nonuniformity of flow field and the reshock on the RM instability. Here, we consider two kinds of nonuniformity: the nonuniform coefficient
5. Prospect
Now, we investigated the Richtmyer‐Meshkov instability and the induced turbulent mixing in fluid flow by using large‐eddy simulation, but the turbulent mixing is a complicated three‐dimensional problem with multiple time‐space scales, and the more engineering applications involve the materials mixing with strength. So, the future work will be carried out in the following aspects: (a) the direct numerical simulations with high precision and high resolution on the platform of supercomputer, (b) the interface instability and turbulent materials mixing with strength, this may involve the material properties such as the deformation, fracture, melt, phase transition, material microstructures, etc., and it may make this problem to be more difficult.
Acknowledgments
Professor Jingsong Bai and Associate professor Tao Wang from the Institute of Fluid Physics, China Academy of Engineering Physics, People’s Republic of China and some other from the same institute should be thanked for the contributions and helps to these works, such as Prof. Ping Li, Dr. Liyong Zou, Kun Liu, Mr. Jinhong Liu, Min Zhong and Bing Wang, and Jiaxin Xiao.
The authors would like to thank the support from “Science Challenge Project” (no. TZ2016001), the National Natural Science Foundation of China (nos. 11372294, 11202195, 11072228, and 11532012), the Science Foundation of China Academy of Engineering Physics (nos. 2011B0202005, 2011A0201002, and 2008B0202011), and the Foundation of National Key Laboratory of Shock Wave and Detonation Physics (no. 9140C670301150C67290).
References
- 1.
Richtmyer RD. Taylor instability in a shock acceleration of compressible of fluids. Communications on Pure and Applied Mathematics. 1960; 13 (2):297–319 - 2.
Meshkov EE. Instability of the interface of two gases accelerated by a shock wave. Izv Akad Nauk SSSR Mekh Zhidk Gaza. 1969; 4 (5):151–157 - 3.
Zoldi CA. Simulation of a shock‐accelerated gas cylinder and comparison with experiment images and velocity fields. In: 8th International Workshop on the Physics of Compressible Turbulent Mixing; 9–14 December; Pasadena, California, USA, 2001. p. 12 - 4.
Amendt P, Colvin JD, Tipton RE, et al. Indirect‐drive noncryogenic double‐shell ignition targets for the National Ignition Facility: Design and analysis. Physics of Plasmas. 2002; 9 :2221–2233 - 5.
Khokhlov AM, Oran ES, Thomas GO. Numerical simulation of deflagration‐to‐detonation transition: The role of shock‐flame interaction in turbulence flames. Combustion and Flames. 1999; 117 :323–339 - 6.
Waitz IA, Marble FE, Zukoski EE. Investigation of a contoured wall injector for hypervelocity mixing augmentation. AIAA Journal. 1993; 31 :1014–1021 - 7.
Maran SP, Sonneborn G, Pun CSJ, et al. Physical conditions in circumstellar gas surrounding SN 1987A 12 years after outburst. Journal of Astrophysical. 2000; 545 :390–398 - 8.
Arnett WD, Bahcall JN, Kirshner RP, Woosley SE. Supernova 1987 A. Annual Review of Astronomy Astrophysics. 1989; 27 :629 - 9.
Bai JS, Li P, Tan DW. Simulations of the instabilities experiment s in St ratified cylindrical shells. Chinese Physics Letter. 2006; 23 (7):1850–1852 - 10.
Bai JS, Wang T, Li P, Zou LY, Liu CL. Numerical simulation of the hydrodynamic instability experiments and flow mixing. Science in China Series G: Physics, Mechanics & Astronomy. 2009; 52 (12):2027–2040 - 11.
Wang T, Bai JS, Li P, Liu K. Large‐eddy simulations of the Richtmyer‐Meshkov instability of rectangular interface accelerated by shock waves. Science in China Series G: Physics, Mechanics & Astronomy. 2010; 53 (5):905–914 - 12.
Bai JS, Liu JH, Wang T, Zou LY, Li P, Tan DW. Investigation of the Richtmyer‐Meshkov instability with double perturbation interface in nonuniform flows. Physical Review E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics. 2010; 81 (2):056302 - 13.
Smagorinsky J. General circulation experiments with the primitive equations. Monthly Weather Review. 1963; 91 :99–164 - 14.
Vreman AW. An eddy‐viscosity subgrid‐scale model for turbulent shear flow: Algebraic theory and applications. Physics of Fluids. 2004; 16 (10):3670–3681 - 15.
Moin P, Squires K, Cabot W, Lee S. A dynamic subgrid‐scale model for compressible turbulence and scalar transport. Physics of Fluids A (Fluid Dynamics). 1991; 3 (11): 2746–2757 - 16.
Hill DJ, Pantano C, Pullin DI. Large‐eddy simulation and multiscale modelling of a Richtmyer‐Meshkov instability with reshock. Journal of Fluid Mechanics. 2006; 557 :29–61 - 17.
Bai JS. High resolution numerical methods and adaptive mesh refinement algorithm for compressible multi‐fluid dynamics [PhD thesis]. China Academy of Engineering Physics; 2003. (in Chinese) - 18.
Wang T, Bai JS, Li P, et al. Two and three dimensional numerical investigations of the single‐mode Richtmyer‐Meshkov instability. Chinese Journal of High Pressure Physics. 2013; 27 (2):277–286 - 19.
Wang T, Bai JS, Li P, Zhong M. The numerical study of shock‐induced hydrodynamic instability and mixing. Chinese Physics B. 2009; 18 (3):1127–1135 - 20.
Malamud G, Elabz Y, Leinov E, et al. Two‐ and three‐dimensional numerical simulations of the interaction of a Richtmyer‐Meshkov instability induced turbulent mixing zone with several re‐shocks. In: 10th IWPCTM; 17–21 July 2006, Paris, France. 2006. - 21.
Bai JS, Wang T, Liu K, et al. Large‐Eddy simulation of the Three‐dimensional experiment on Richtmyer‐Meshkov instability induced turbulence. International Journal of Astronomy and Astrophysics. 2012; 2 :28–36 - 22.
Poggi F, Thorembey MH, Rodriguez G. Velocity measurements in turbulent gaseous mixtures induced by Richtmyer–Meshkov instability. Physics of Fluids. 1998; 10 :2698–2712 - 23.
Wang T, Li P, Bai JS, Tao G, Wang B, Zou LY. Large‐eddy simulation of the Richtmyer‐Meshkov instability. Canadian Journal of Physics. 2015; 93 (10):1124–1130 - 24.
Bates KR, Nikiforakis N, Holder D. Richtmyer‐Meshkov instability induced by the interaction of a shock wave with a rectangular block of SF6. Physics of Fluids. 2007; 19 :036101 - 25.
Holder DA, Smith AV, Barton CJ, Youngs DL. Shock‐tube experiments on Richtmyer‐Meshkov instability growth using an enlarged double‐bump perturbation. Laser and Particle Beams. 2003; 21 :411–418 - 26.
Li P, Bai JS, Wang T, et al. Large eddy simulation of a shocked gas cylinder instability induced turbulence. Science China Physics, Mechanics & Astronomy. 2010; 53 (2):262–268 - 27.
Mikaelian KO. Rayleigh‐Taylor and Richtmyer‐Meshkov instabilities and mixing in stratified spherical shells. Physical Review A (Atomic, Molecular, and Optical Physics). 1990; 42 (16):3400–3420 - 28.
Bai JS, Wang T, Zou LY, et al. Numerical simulation for shock tube and jelly interface instability experiments. Chinese Journal of Applied Mechanics. 2009; 26 (3):418–425 - 29.
Erez L, Sadot O, Oron D, et al. Study of the membrane effect on turbulent mixing measurements in shock tubes. Shock Waves. 2000; 10 :241–251 - 30.
Wang T, Tao G, Bai JS, et al. Numerical comparative analysis of Richtmyer‐Meshkov instability simulated by different SGS models. Canadian Journal of Physics. 2015; 93 (5):519–525 - 31.
Wang T, Tao G, Bai JS, et al. Analysis of dependence of multi‐mode Richtmyer‐Meshkov instability on initial conditions. Chinese Journal of High Pressure Physics. 2016; 30 (5):380–386 - 32.
Bai JS, Zou LY, Wang T, et al. Experimental and numerical study of the shock‐accelerated elliptic heavy gas cylinders. Physical Review E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics. 2010; 82 (5):056318 - 33.
Bai JS, Wang B, Wang T, Liu K. Numerical simulation of the Richtmyer‐Meshkov instability in initially nonuniform flows and mixing with reshock. Physical Review E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics. 2012; 86 (6):066319 - 34.
Xiao JX, Bai JS, Wang T. Numerical study of initial perturbation effects on Richtmyer‐Meshkov instability in nonuniform flows. Physical Review E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics. 2016; 94 (1):013112