The manipulation of the Richtmyer–Meshkov instability growth at a heavy–light interface via successive shocks is theoretically analysed and experimentally realized in a specific shock-tube facility. An analytical model is developed to forecast the interface evolution before and after the second shock impact, and the possibilities for the amplitude evolution pattern are systematically discussed. Based on the model, the parameter conditions for each scenario are designed, and all possibilities are experimentally realized by altering the time interval between two shock impacts. These findings may enhance the understanding of how successive shocks influence hydrodynamic instabilities in practical applications.

An asymptotic matching modal model is established based on the singular perturbation method for predicting mode evolution in single- and dual-mode interfaces accelerated by a shock wave. The startup process is incorporated into the model to provide a complete description of the mode evolution after the shock impact. Through considering the feedback from high-order harmonic to the third-order harmonic, the model accuracy is improved and the model divergence is prevented. In addition, the model can evaluate the mutual-coupling effect on the amplitude variations of high-order harmonics besides the ‘beat modes’. To validate the model, experiments on both light–heavy and heavy–light interfaces subject to a shock wave are conducted, and both single- and dual-mode interfaces formed by the soap-film technique are involved. The interface profiles extracted from mode decomposition and predicted by the model show high consistency with the experimental counterparts. Good agreement of the mode amplitude growths between the experiments and theoretical predictions shows the superiority of the model, especially for the heavy–light interface.

The effects of reshock conditions, including the interface evolution state before reshock and the second shock intensity, on interface instability induced by two successive shocks propagating in the same direction are investigated via shock-tube experiments. It is observed that the reshock promotes the interface instability, and the post-reshock perturbation evolution relates to both the pre-reshock interface evolution state and second shock intensity. For the linear evolution of the twice-shocked interface, existing models perform poorly when either the pre-reshock interface shape effect or the secondary compression effect is pronounced, as current reduction factors fail to accurately describe these effects. Besides, the reshock-induced linear amplitude growth rate shows a non-monotonic dependence on the scaled pre-reshock amplitude, primarily due to the shape effect of the pre-reshock interface. For the post-reshock nonlinear evolution, the model proposed by Zhang & Guo (J. Fluid Mech., vol. 786, 2016, pp. 47–61) offers reasonable predictions when the second shock is weak. However, when the second shock is moderately strong, the model overestimates the bubble growth and underestimates the spike evolution under the influence of the significant secondary compression effect. Furthermore, empirical linear and nonlinear models capable of describing the dependence of the post-reshock evolution on reshock conditions are proposed based on the present experimental results and existing models.

The coupling of Richtmyer–Meshkov instability (RMI) and Kelvin–Helmholtz instability (KHI), referred to as RM-KHI, on a shock-accelerated inclined single-mode air–SF$_6$ interface is studied through shock-tube experiments, focusing on the evolution of the perturbation distributed along the inclined interface. To clearly capture the linear (overall linear to nonlinear) evolution of RM-KHI, a series of experiments with a weak (relatively strong) incident shock is conducted. For each series of experiments, various $\theta _{i}$ (angle between incident shock and equilibrium position of the initial interface) are considered. The nonlinear flow features manifest earlier and develop faster when $\theta _{i}$ is larger and/or shock is stronger. In addition, the interface with $\theta _{i}>0^{\circ }$ evolves obliquely along its equilibrium position under the effect of KHI. RMI dominates the early-time amplitude evolution regardless of $\theta _{i}$ and shock intensity, which arises from the discrepancy in the evolution laws between RMI and KHI. KHI promotes the post-early-stage amplitude growth and its contribution is related positively to $\theta _{i}$. An evident exponential-like amplitude evolution behaviour emerges in RM-KHI with a relatively strong shock and large $\theta _{i}$. The linear model proposed by Mikaelian (Phys. Fluids, vol. 6, 1994, pp. 1943–1945) is valid for RM-KHI within the linear period. In contrast, the adaptive vortex model (Sohn et al., Phys. Rev. E, vol. 82, 2010, p. 046711) can effectively predict both the interface morphology and overall amplitude evolutions from the linear to nonlinear regimes.

This paper characterizes the refraction of a triple-shock configuration at planar fast–slow gas interfaces. The primary objective is to reveal the wave configurations and elucidate the mechanisms governing circulation deposition and velocity perturbation on the interface caused by triple-shock refraction. The incident triple-shock configuration is generated by diffracting a planar shock around a rigid cylinder, and four interfaces with various $Z_{{r}}$ (i.e. acoustic impedance ratio across the interface) are considered. An analytical model describing the triple-shock refraction is developed, which accurately predicts both the wave configurations as well as circulation deposition and velocity perturbation. Depending on $Z_{{r}}$, three distinct patterns of transmitted waves can be anticipated: a triple-shock configuration; a four-shock configuration; a four-wave configuration. The underlying mechanism for the formation of these wave configurations is elucidated through shock polar analysis. A novel physical insight into the contribution of triple-shock refraction to the interface perturbation growth is provided. The results indicate that the reflected shock in an incident triple-shock configuration makes significant negative contribution to both circulation deposition and velocity perturbation. This investigation elucidates the underlying mechanism responsible for the relatively insignificant contribution of baroclinic circulation to the Richtmyer–Meshkov-like instability induced by a non-uniform shock, and provides an explanation for the decrease in growth rate of interface perturbation amplitude with increasing Atwood number.

Shock-tube experiments and theoretical studies have been performed to highlight mode-coupling in an air–SF$_6$–air fluid layer. Initially, the two interfaces of the layer are designed as single mode with different basic modes. It is found that as the two perturbed interfaces become closer, interface coupling induces a different mode from the basic mode on each interface. Then mode coupling further generates new modes. Based on the linear model (Jacobs et al., J. Fluid Mech., vol. 295, 1995, pp. 23–42), a modified model is established by considering the different accelerations of two interfaces and the waves’ effects in the layer, and provides good predictions to the linear growth rates of the basic modes and the modes generated by interface coupling. It is observed that interface coupling behaves differently to the nonlinear growth of the basic modes, which can be characterized generally by the existing or modified nonlinear model. Moreover, a new modal model is established to quantify the mode-coupling effect in the layer. The mode-coupling effect on the amplitude growth is negligible for the basic modes, but is significant for the interface-coupling modes when the initial wavenumber of one interface is twice the wavenumber of the other interface. Finally, amplitude freeze-out of the second single-mode interface is achieved theoretically and experimentally through interface coupling. These findings may be helpful for designing the target in inertial confinement fusion to suppress the hydrodynamic instabilities.

We report the first experiments on hydrodynamic instabilities of a single-mode light/heavy interface driven by co-directional rarefaction and shock waves. The experiments are conducted in a specially designed rarefaction-shock tube that enables the decoupling of interfacial instabilities caused by these co-directional waves. After the impacts of rarefaction and shock waves, the interface evolution transitions into Richtmyer–Meshkov unstable states from Rayleigh–Taylor (RT) stable states, which is different from the finding in the previous case with counter-directional rarefaction and shock waves. A scaling method is proposed, which effectively collapses the RT stable perturbation growths. An analytical theory for predicting the time-dependent acceleration and density induced by rarefaction waves is established. Based on the analytical theory, the model proposed by Mikaelian (Phys. Fluids, vol. 21, 2009, p. 024103) is revised to provide a good description of the dimensionless RT stable behaviour. Before the shock arrival, the unequal interface velocities, caused by rarefaction-induced uneven vorticity, result in a V-shape-like interface. The linear growth rate of the amplitude is insensitive to the pre-shock interface shape, and can be well predicted by the linear superposition of growth rates induced by rarefaction and shock waves. The nonlinear growth rate is higher than that of a pure single-mode case, which can be predicted by the nonlinear models (Sadot et al., Phys. Rev. Lett., vol. 80, 1998, pp. 1654–1657; Dimonte & Ramaprabhu, Phys. Fluids, vol. 22, 2010, p. 014104).

Richtmyer–Meshkov instability (RMI) at a light–heavy single-mode interface over a wide range of post-shock Atwood numbers $A_1$ is studied systematically through elaborate experiments. The interface generation and $A_1$ variation are achieved by the soap-film technology and gas-layer scheme, respectively. Qualitatively, the nonlinear interface evolution features, including spike, bubble and roll-up structures, are more significant in RMI with higher $A_1$. Quantitatively, both the impulsive model and an analytical linear model perform well in predicting the linear growth rate under a wide range of $A_1$ conditions. For the weakly nonlinear stage, the significant spike acceleration occurring when $A_1$ is high, which is observed experimentally for the first time, results in the evolution law of RMI with high $A_1$ being different from the counterpart with low or intermediate $A_1$. None of the considered nonlinear models is found to be applicable for RMI under all $A_1$ conditions, and the predictive capabilities of these models are analysed and summarized. Based on the present experimental results, an empirical nonlinear model is proposed for RMI over a wide range of $A_1$. Further, modal analysis shows that in RMI with high (low or intermediate) $A_1$, high-order harmonics evolve rapidly (slowly) and cannot (can) be ignored. Accordingly, for RMI with high (low or intermediate) $A_1$, the modal model proposed by Zhang & Sohn (Phys. Fluids, vol. 9, 1997, pp. 1106–1124) is less (more) accurate than the one proposed by Vandenboomgaerde et al. (Phys. Fluids, vol. 14, 2002, pp. 1111–1122), since the former ignores perturbation solutions higher than fourth order (the latter retains only terms with the highest power in time).

Attenuation and even freeze-out (amplitude growth stagnation) of the perturbation amplitude growth of a shocked SF$_6$–air interface are first realized in shock-tube experiments through reflected rarefaction waves, which produce reverse baroclinic vorticity offsetting the vorticity deposited by the shock. A theoretical model is constructed to predict the perturbation growth after the impact of rarefaction waves, and seven possibilities of amplitude growth are analysed. Experimentally, a planar air–helium interface is used to produce reflected rarefaction waves. Through changing the perturbation wavelength and the time interval of two impacts, five experiments with specific initial conditions are carried out, and three different possibilities of perturbation growth attenuation are realized.

We report the first shock-tube experiments on two-dimensional dual-mode air–SF$_6$ interfaces with different initial spectra subjected to a convergent shock wave. The convergent shock tube is specially designed with a tail opening to highlight the Bell–Plesset (BP) and mode-coupling effects on amplitude development of fundamental mode (FM). The results show that the BP effect promotes the occurrence of mode coupling, and the feedback of high-order modes to the FM also arises earlier in convergent geometry than that in its planar counterpart. Relatively, the amplitude growth of the FM with a higher mode number is inhibited by the feedback, and saturates earlier. The FM with a lower mode number is affected more heavily by the BP effect, and finally dominates the flow. A new model is proposed to well predict the amplitude growths of the FM and high-order modes in convergent geometry. In particular, for FM that reaches its saturation amplitude, the post-saturation relation is introduced in the model to achieve a better prediction.

‘Freeze-out’ of amplitude growth, i.e. the amplitude growth stagnation of a shocked helium–air interface, is realized through a reflected shock, which produces baroclinic vorticity of the opposite sign to that deposited by the first shock. Theoretically, a model is constructed to calculate the relations among the initial parameters for achieving freeze-out. In particular, if the amplitude growth is within the linear regime at the arrival of the reflected shock, the time interval between the impacts of two shock waves is linearly related to the initial perturbation wavelength, and is independent of the initial perturbation amplitude. Experimentally, an air–SF$_6$ (or air–argon) plane interface is adopted to produce a weak reflected shock. Seven experimental runs with specific initial conditions are examined. For all cases, freeze-out is achieved after the reflected shock impact under the designed conditions.

Richtmyer–Meshkov instability induced by two successive shock waves is experimentally studied in a specific shock tube. To create two successive shock waves synchronously, a driver section is added between the driver and driven sections of the standard shock tube, and an electronically controlled membrane rupture equipment is adopted. The shock-tube flow after the membranes rupture is well described by combining the shock relations, isentropic wave relations with compatibility relations across the contact surface (region). The new shock tube is capable of generating two successive shock waves with controllable strengths and time interval, and provides a relatively ‘clean’ wave system. Then the developments of single-mode light–heavy interfaces with different initial conditions induced by two successive shock waves are investigated. The initial amplitudes are all small enough such that the first-shocked interface is within the linear growth regime at the arrival of the second shock. The results show that if the pre-second-shock perturbation amplitude is small, the linear, nonlinear and modal evolutions of the double-shocked interface can be reasonably predicted by the existing models proposed for predicting the perturbation growth induced by a single shock. For the double-shocked interface, the second shock provides an additional perturbation velocity field to the original one introduced by the first shock impact. The validity of the linear superposition model indicates that the linear superposition of these two perturbation velocity fields is satisfied. Therefore, a double-shocked interface evolves similarly to a single-shocked interface provided that their postshock amplitudes and linear growth rates are the same.

The two-dimensional (2-D) unsteady shock reflection over a single wedge is studied theoretically and numerically, and the emphasis is placed on the trajectory of the triple point (TP). Skews’ relation and the three-shock theory are, respectively, used for determining the trajectory angles of the corner-generated disturbance and the TP, and, subsequently, a model capable of predicting the TP trajectory is established for 2-D unsteady shock reflections over a single wedge. Then, a systematically numerical study is carried out on the 2-D unsteady shock reflection over a single wedge, including five types of shock reflection with the wedge angle increased and five types of shock reflection with the wedge angle decreased. It is found that the Mach stem is always slightly concavely (convexly) curved for reflections with the wedge angle increased (decreased), which should be caused by corner-generated compression (rarefaction) waves propagating along the Mach stem. The new model reasonably predicts the TP trajectory in all the 2-D unsteady shock reflections, and its performance is related not only to the variation trend of the wedge angle (increase or decrease), but also to the type of shock and the initial wedge angle. Specifically, for the shock reflection with the wedge angle increased (decreased), the model generally provides a slightly better (worse) prediction if the TP trajectory angle calculated by the three-shock theory is compulsively modified. The shock-shock/expansion polar analysis presented can partially explain the model performance for predicting the TP trajectory.

Shock-tube experiments on eight kinds of two-dimensional multi-mode air–SF$_6$ interface with controllable initial conditions are performed to examine the dependence of perturbation growth on initial spectra. We deduce and demonstrate experimentally that the amplitude development of each mode is influenced by the mode-competition effect from quasi-linear stages. It is confirmed that the mode-competition effect is closely related to initial spectra, including the wavenumber, the phase and the initial amplitude of constituent modes. By considering both the mode-competition effect and the high-order harmonics effect, a nonlinear model is established based on initial spectra to predict the amplitude growth of each individual mode. The nonlinear model is validated by the present experiments and data in the literature by considering diverse initial spectra, shock intensities and density ratios. Moreover, the nonlinear model is successfully extended based on the superposition principle to predict the growths of the total perturbation width and the bubble/spike width from quasi-linear to nonlinear stages.

The convergent Richtmyer–Meshkov instability (RMI) of an $\textrm {SF}_6$ layer with a uniform outer surface and a sinusoidal inner surface surrounded by air (generated by a novel soap film technique) is studied in a semiannular convergent shock tube using high-speed schlieren photography. The outer interface initially suffers only a slight deformation over a long period of time, but distorts quickly at late stages when the inner interface is close to it and produces strong coupling effects. The development of the inner interface can be divided into three stages. At stage I, the interface amplitude first drops suddenly to a lower value due to shock compression, then decreases gradually to zero (phase reversal) and later increases sustainedly in the negative direction. After the reshock (stage II), the perturbation amplitude exhibits long-term quasi-linear growth with time. The quasi-linear growth rate depends weakly on the pre-reshock amplitude and wavelength, but strongly on the pre-reshock growth rate. An empirical model for the growth of convergent RMI under reshock is proposed, which reasonably predicts the present results and those in the literature. At stage III, the perturbation growth is promoted by the Rayleigh–Taylor instability caused by a rarefaction wave reflected from the outer interface. It is found that both the Rayleigh–Taylor effect and the interface coupling depend heavily on the layer thickness. Therefore, controlling the layer thickness is an effective way to modulate the late-stage instability growth, which may be useful for the target design.

We report the first experiments on divergent shock-driven Richtmyer–Meshkov instability (RMI) at well-controlled single-mode interfaces. These experiments are performed in a novel divergent shock tube designed by shock dynamics theory. Generally, the perturbation growth can be divided into three successive stages: linear growth, quick reduction in growth rate and instability freeze-out. It is observed that the growth rate at each stage is far lower than its counterpart in planar or convergent geometry due to geometric divergence. We also found that nonlinearity is much weaker than that in planar or convergent RMI, and has a negligible influence on the overall amplitude growth even at late stages when it has become strong. This weak nonlinear effect is because the growth of the third harmonic counteracts its feedback to the fundamental mode. As a consequence, the linear theory of Bell (report no. LA-1321) accounting for geometric divergence and Rayleigh–Taylor (RT) stabilization caused by flow deceleration can reasonably predict the present results from early to late stages. The instability freeze-out at late times is ascribed to the negative growth induced by geometric divergence and RT stabilization, and is also well reproduced by the linear theory.

Richtmyer–Meshkov instability of the SF6 gas layer surrounded by air is experimentally investigated. Using the soap film technique, five kinds of gas layer with two sharp interfaces are generated such that the development of each individual interface is highlighted. The flow patterns are determined by the amplitudes and phases of two corrugated interfaces. For a layer with both interfaces planar, the interface velocity shows that the reflected rarefaction waves from the second interface accelerate the first interface motion. For a layer with the second interface corrugated but the first interface planar, the reflected rarefaction waves make the first interface develop with the same phase as the second interface. For a layer with the first interface corrugated but the second interface planar, the rippled shock seeded from the first interface makes the second interface develop with the same phase as the first interface and the layer evolves into an ‘upstream mushroom’ shape. For two interfaces corrugated with opposite (the same) phase but a larger amplitude for the first interface, the layer evolves into ‘sinuous’ shape (‘bow and arrow’ shape, which has never been observed previously). For the interface amplitude growth in the linear stage, the waves’ effects are considered in the model to give a better prediction. In the nonlinear stage, the effect of the rarefaction waves on the first interface evolution is quantitatively evaluated, and the nonlinear growth is well predicted. It is the first time in experiments to quantify the interfacial instability induced by the rarefaction waves inside the heavy gas layer.

The interaction of rarefaction waves and a heavy/light interface is investigated using numerical simulations by solving the compressible Euler equations. An upwind space–time conservation element and solution element (CE/SE) scheme with second-order accuracy in both space and time is adopted. Rarefaction waves are generated by simulating the shock-tube problem. In this work, the SF6/air interface evolution under different conditions is considered. First, the gas physical parameters before and after the rarefaction waves impact the interface are calculated using one-dimensional gas dynamics theory. Then, the interaction between the rarefaction waves and a single-mode perturbation interface is investigated, and both the interface evolution and the wave patterns are obtained. Afterwards, the amplitude growth of the interface over time is compared between cases, considering the effects of the interaction period and the strength of the rarefaction waves. During the interaction of the rarefaction waves with the interface, the Rayleigh–Taylor instability induced by the rarefaction waves is well predicted by modifying the nonlinear model proposed by Zhang & Guo (J. Fluid Mech., vol. 786, 2016, pp. 47–61), considering the variable acceleration. After the rarefaction waves leave the interface, the equivalent Richtmyer–Meshkov instability is well depicted by the nonlinear model proposed by Zhang et al. (Phys. Rev. Lett., vol. 121(17), 2018, 174502), considering the growth rate transition from Rayleigh–Taylor instability to Richtmyer–Meshkov instability. The differences in the heavy/light interface amplitude growth under the rarefaction wave condition and the shock wave condition are compared. The interface perturbation is shown to be more unstable under rarefaction waves than under a shock wave.

The regular reflection to Mach reflection ($\text{RR}\rightarrow \text{MR}$) transition in inviscid perfect air for shock reflection over convex and straight wedges is investigated. Provided that the variation of shock intensity only has a second-order effect on the wave transition, the possible cases for the occurrence of the $\text{RR}\rightarrow \text{MR}$ transition for curved shock reflection over a wedge are discussed. For a planar shock reflecting over a convex wedge, four different flow regions are classified and the mechanism of the disturbance propagation is interpreted. It is found that the flow-induced rarefaction waves exist between disturbances generated from neighbouring positions and isolate them. For a curved shock reflecting over a convex wedge, although the distributions of the flow regions are different from those in planar shock reflection, the analysis for the planar shock case can be extended to curved shock cases as long as the wedge is convex. When a diverging shock reflects off a straight wedge, the flow-induced rarefaction waves are absent. However, the disturbances generated earlier cannot overtake the reflection point before the pseudo-steady criterion is satisfied. In the cases considered, the flow in the vicinity of the reflection point will not be influenced by the unsteady flow caused by the shock intensity and the wedge angle variations. This is clearly a local property of the shock–wall interaction, no matter what the history of the shock trajectory is. For validation, extensive inviscid numerical simulations are performed, and the numerical results show the reliability of the pseudo-steady criterion for predicting the $\text{RR}\rightarrow \text{MR}$ transition on a convex wedge.

The evolution of an $\text{SF}_{6}$ layer surrounded by air is experimentally studied in a semi-annular convergent shock tube by high-speed schlieren photography. The gas layer with a sinusoidal outer interface and a circular inner interface is realized by the soap-film technique such that the initial condition is well controlled. Results show that the thicker the gas layer, the weaker the interface–coupling effect and the slower the evolution of the outer interface. Induced by the distorted transmitted shock and the interface coupling, the inner interface exhibits a slow perturbation growth which can be largely suppressed by reducing the layer thickness. After the reshock, the inner perturbation increases linearly at a growth rate independent of the initial layer thickness as well as of the outer perturbation amplitude and wavelength, and the growth rate can be well predicted by the model of Mikaelian (Physica D, vol. 36, 1989, pp. 343–357) with an empirical coefficient of 0.31. After the linear stage, the growth rate decreases continuously, and finally the perturbation freezes at a constant amplitude caused by the successive stagnation of spikes and bubbles. The convergent geometry constraint as well as the very weak compressibility at late stages are responsible for this instability freeze-out.