1. INTRODUCTION
The interaction of a shock wave with an initially perturbed interface separating two fluids with different properties is generally referred to as the Richtmyer-Meshkov (RM) instability. The pioneering theoretical and numerical analysis on this shock-induced instability was performed by Richtmyer (Reference Richtmyer1960) and then confirmed by Meshkov (Reference Meshkov1969) in a shock tube. From then on, a large number of theoretical, numerical, and experimental studies on the RM instability have been carried out because of its wide physical applications such as inertial confinement fusion (Lindl et al., Reference Lindl, Mccrory and Campbell1992), turbulent mixing in scramjet (Yang et al., Reference Yang, Kubota and Zukoski1993), and collapse in supernova (Arnett et al., Reference Arnett, Bahcall, Kirshner and Woosley1989) and its great academic significance in interfacial instability, vortex dynamics, and formation of compressible turbulence. Several comprehensive reviews have been made to summarize the progress of studies on the RM instability in the past years (Zabusky, Reference Zabusky1999; Brouillette, Reference Brouillette2002; Ranjan et al., Reference Ranjan, Oakley and Bonazza2011). However, most previous researches mainly focused on planar shock cases and the studies on the RM instability induced by converging shock waves are still desirable for many practical applications. Exceptional studies include numerical simulations (Mikaelian,Reference Mikaelian1995; Zhang & Graham, Reference Zhang and Graham1998; Glimm et al., Reference Glimm, Grove, Zhang and Dutta2002) and experiments (Lanier et al., Reference Lanier, Barnes, Batha, Dar, Magelssen, Scott, Dunne, Parker and Rothman2003; Hosseini & Takayama, Reference Hosseini and Takayama2005a; Fincke et al., Reference Fincke, Lanier, Batha, Hueckstaedt, Magelssen, Rothman, Parker and Horsfield2005), which have reported that the convergence can cause the apparent postponement of saturation and the late-time increase in growth rate. Furthermore, the converging shock can result in extreme conditions near the focal point. It is therefore significant to explore the characteristics of the RM instability with respect to converging shock cases.
Shock tube is one of the most used facilities to produce a shock wave in laboratory situations. In recent decades, attempts have been made to generate converging shock waves in shock tube environments. Perry and Kantrowitz (Reference Perry and Kantrowitz1951) established a horizontal annular coaxial shock tube for the first time to observe the cylindrical converging shock. Later, Takayama et al. (Reference Takayama, Kleine and Gronig1987) built a horizontal converging shock tube with struts supporting the inner core. Subsequently, a vertical coaxial diaphragmless shock tube was constructed to produce uniform cylindrical converging shock waves (Hosseini et al., Reference Hosseini, Onodera and Takayama2000). Cylindrical converging shock wave was also generated in a shock tube with a shock-shaping end section (Apazidis & Lesser, Reference Apazidis and Lesser1996; Apazidis et al., Reference Apazidis, Lesser, Tillmark and Johansson2002; Kjellander et al., Reference Kjellander, Tillmark and Apazidis2011). Besides, Dimotakis and Samtaney (Reference Dimotakis and Samtaney2005) reported a gas lens technique in a two-dimensional wedge geometry and Hosseini and Takayama (Reference Hosseini and Takayama2005b) designed an aspheric lens-shaped transparent test section to produce converging shock waves. It must be mentioned that so far the RM instability experiments have seldom been performed in such shock tubes probably due to the difficulties in settling the initial interfaces and measuring relevant quantities. Recently, a simple but effective technique for generating cylindrical converging shock waves in an ordinary shock tube was developed based on shock dynamics theory in our group (Zhai et al., Reference Zhai, Liu, Qin, Yang and Luo2010; Luo et al., Reference Luo, Si, Yang and Zhai2014). The present experimental study is motivated by the method, aiming at more potential of this experimental method in study of the RM instability.
In experimental aspect of the RM instability study, the generation of initially perturbed interface and the employment of diagnostic methods are very important. Interfaces used in previous experiments can mainly be classified into two types, i.e., with and without membranes. Interfaces with membrane consist of soap film interface (Haas & Sturtevant, Reference Haas and Sturtevant1987; Ranjan et al., Reference Ranjan, Anderson, Oakley and Bonazza2005; Layes et al., Reference Layes, Jourdan and Houas2009; Zhai et al., Reference Zhai, Si, Luo and Yang2011; Si et al., Reference Si, Zhai, Luo and Yang2012; Wang et al., Reference Wang, Si and Luo2013; Luo et al., Reference Luo, Wang and Si2013) and microfilm interface (Bates & Nikiforakisb, Reference Bates and Nikiforakisb2007; Mariani et al., Reference Mariani, Vandenboomgaerde, Jourdan, Souffland and Houas2008). In these methods, the shape of interface with sharp boundaries can accurately be determined. The other type of interface is considered to be membraneless interface, which includes gas cylinders (Jacobs, Reference Jacobs1993; Tomkins et al., Reference Tomkins, Kumar, Orlicz and Prestridge2008; Zou et al., Reference Zou, Liu, Tan, Huang and Luo2010), gas curtains (Jacobs et al., Reference Jacobs, Klein, Jenkins and Benjamin1992; Orlicz et al., Reference Orlicz, Balakumar, Tomkins and Prestridge2009; Balasubramanian et al., Reference Balasubramanian, Orlicz, Balakumar and Prestridge2012), single-mode interfaces by vibration (Sheeley & Jcobs, Reference Sheeley and Jacobs1995; Jones & Jacobs, Reference Jones and Jacobs1997; Jacobs & Krivets, Reference Jacobs and Krivets2005; Long et al., Reference Long, Krivets, Greenough and Jacobs2009), interfaces by retracting metal plate initially separating two different gases (Puranik et al., Reference Puranik, Oakley, Anderson and Bonazza2004) and some others. In these methods, the influences of membrane and support would be avoided, but gaseous diffusion often happens initially between two fluids and the exact initial conditions have to be determined afterward. From the above-mentioned methods, different shapes of gaseous interface can be obtained, including single and multi-mode interfaces (Mariani et al., Reference Mariani, Vandenboomgaerde, Jourdan, Souffland and Houas2008; Long et al., Reference Long, Krivets, Greenough and Jacobs2009; Balasubramanian et al., Reference Balasubramanian, Orlicz, Balakumar and Prestridge2012), spherical bubbles (Ranjan et al., Reference Ranjan, Oakley and Bonazza2011), circular cylinders (Hass & Sturtevant, Reference Haas and Sturtevant1987; Hosseini & Takayama, Reference Hosseini and Takayama2005a; Tomkins et al., Reference Tomkins, Kumar, Orlicz and Prestridge2008), elliptical cylinders (Zou et al., Reference Zou, Liu, Tan, Huang and Luo2010), and polygonal blocks (Bates & Nikiforakisb, Reference Bates and Nikiforakisb2007; Wang et al., Reference Wang, Si and Luo2013). As to the diagnostics, various methods such as shadowgraphy (Hass & Sturtevant, Reference Haas and Sturtevant1987; Layes et al., Reference Layes, Jourdan and Houas2009), schlieren (Zhai et al., Reference Zhai, Si, Luo and Yang2011), double-exposure holographic interferometry (Hosseini & Takayama, Reference Hosseini and Takayama2005a), and laser sheet imaging (Brouillette, Reference Brouillette2002; Ranjan et al., Reference Ranjan, Oakley and Bonazza2011; Orlicz et al., Reference Orlicz, Balakumar and Prestridge2013) have been used. As a result, rich phenomena and relevant quantities of the distorted interface can be obtained for understanding the physical rules and mechanisms of the RM instability.
When a shock wave passes through an initially perturbed interface, baroclinic vorticity will be generated due to the misalignment between the density and pressure gradients. The development behaviors of gaseous interfaces are closely dependent on the amplitude and distribution of baroclinic vorticity, which is attributed to the shock strength and the initial shapes of both interface and shock wave. In our previous study (Si et al., Reference Si, Zhai, Luo and Yang2014), some preliminary experiments were performed based on the interaction of a cylindrical shock with a heavy-gas cylinder and it was found that the results are different from ones with respect to the planar shock. In this work, we will revisit two types of interface used in planar shock cases, i.e., gas bubble and gas cylinder, and settle them in the cylindrical converging shock tube to investigate the evolution of such interfaces accelerated by the converging shock and the reflected shock formed from the center of convergence.
2. EXPERIMENTAL METHODS
2.1. Generation of Cylindrical Shock Waves
Here we will restrict ourselves to the RM instability experiments and briefly describe the generation of cylindrical shock waves. For the detailed mathematical formulism and validation, the reader is referred to our previous work (Zhai et al., Reference Zhai, Liu, Qin, Yang and Luo2010; Reference Zhai, Si, Luo, Yang, Liu, Tan and Zou2012; Luo et al., Reference Luo, Si, Yang and Zhai2014). Figure 1 illustrates the schematic of the test section which can directly convert an incident planar shock wave into a cylindrical one. The test section mainly consists of a straight wall (OAA′O′) and a curved wall (OBB′O′) designed based on the shock dynamics theory (Zhai et al., Reference Zhai, Liu, Qin, Yang and Luo2010). The wall profile (curved line from point B (B′) to point Q (Q′)) can be obtained for a given group of controllable parameters including the converging angle θ0, the incident planar shock Mach number M 0 and the shock tube height h. In order to optimize the experimental observation and measurement, a sufficient converging part (OPQ-O′P′Q′) is needed. Generally, large M 0, large h, and moderate θ0 will be favorable for performing RM instability experiments (Zhai et al., Reference Zhai, Si, Luo, Yang, Liu, Tan and Zou2012).
Experiments are conducted in two sets of horizontal rectangular shock tubes with different h, respectively. The first shock tube (ST-1) consists of a 1.7 m driver section and a 2.5 m driven section with cross-sectional area of 70 mm × 40 mm (i.e., h = 70 mm or 40 mm), and the second shock tube (ST-2) consists of a 2.0 m driver section and a 7.0 m driven section with cross-sectional area of 95 mm × 95 mm (i.e., h = 95 mm). For each shock tube, the controllable parameters M 0, h, and θ0 are chosen to calculate the wall profile and the test section is then manufactured accurately by the linear cutting technique. Initially, both the driver and driven sections of the ST-1 or ST-2 are full of air at atmospheric pressure and the driver section is further filled with nitrogen to generate the incident planar shock by the rupture of a polypropylene diaphragm initially separating two sections. The velocity (or the Mach number) of the incident shock is determined by the signals of two piezoelectric pressure transducers mounted on the shock tube side wall.
2.2. Formation of Initial Interfaces
The experimental facility is capable of settling most above-mentioned interfaces with arbitrary shapes. In this work, we mainly focus on spherical and cylindrical gaseous interfaces making use of heavy tested gas (i.e., SF6), as sketched in Figure 2. Soap film is used to create the spherical gaseous interface. Soap film component is very thin (about 0.25–1 μm), flimsy, and durable so that its presence may have very little influence on experiments. The spherical gaseous interface is formed by filling a soap bubble with SF6 and the bubble is hung on a top support in the test section. The gas cylinder flowing vertically under gravity into the test section is created in a careful procedure similar to previous studies (Tomkins et al., Reference Tomkins, Kumar, Orlicz and Prestridge2008; Si et al., Reference Si, Zhai, Luo and Yang2014). Prior to experiments glycol droplets produced by a commercial theatrical fog generator are well mixed with the tested gas and served as tracer. The accurate track of glycol droplets in SF6 has been confirmed by Prestridge et al., (Reference Prestridge, Rightley, Vorobieff, Benjamin and Kurnit2000). The heavy gas cylinder flows slowly (~0.1 m/s) through a round nozzle on the top wall and is sucked mildly through a plenum on the bottom wall by a vacuum pump. In this way the gas cylinder moves steadily with smooth edges.
2.3. Flow Diagnostics
A high-speed video camera (FASTCAM SA5, Photron Limited) is utilized to record the evolution of interfaces based on two optical systems. Figure 3 presents the schematics of the shock tube and flow diagnostic system. The first one is the schlieren photography as shown in Figure 3a, in which a Z-fold schlieren system made up of a knife-edge, a slit, two convex lenses, and two concave mirrors is adopted and the flow field is illuminated by a DC regulated light source (DCR III, SCHOTT North America, Inc., 200 W). Another optical system is the planar Mie scattering photography as sketched in Figure 3b, in which a laser-sheet illuminating the flow is provided by a continuous laser (SDL-532-15000T, Shanghai Dream Lasers technology Co. Ltd., 15 W, 532 nm) combined with a cylindrical concave mirror and a convex lens. An inclined flat mirror located in front of the test section is used to easily adjust the position of the laser-sheet and make sure that the laser-sheet is perpendicular to the axis of gas cylinder (see Fig. 2b). The timing and triggering system involves a four channel delay generator (DG645, Stanford Research Systems), two piezoelectric pressure transducers, a charge amplifier, an oscilloscope, and some accessories. The schlieren photography records the shock propagation and interface deformation in an integrated view, while the planar Mie scattering photography monitors the flow based on the additional tracer in a cross-sectional view. Depending on the initial conditions in experiments, different optical systems are employed.
3. RESULTS AND DISCUSSION
3.1. Gas Bubble
The shock-bubble interaction is a fundamental configuration for studying the RM instability and the planar shock case has been well understood (Ranjan et al., Reference Ranjan, Oakley and Bonazza2011). In this work, the experiment is conducted in the ST-1 with chosen parameters M 0 = 1.2, θ0 = 15° and h = 70 mm. The initial Mach number at point Q and the radius of the converging part are calculated from the shock dynamics theory as M q = 1.29 and R q = 106 mm, respectively. The bubble diameter D 0 = 16 mm and the distance from the bubble center to the center of convergence L 0 = 86 mm are measured, respectively. A pair of optical windows of 70 mm × 150 mm are located at the position where the converging part is in sight, as shown in Figure 2b. The flow is monitored by the high-speed schlieren photography, as indicated in Figure 3a. The high-speed video camera with a frame rate of 50,000 fps corresponding to the time interval between two consecutive frames of 20 μs, a spatial resolution of 192 × 704 and a very fast shutter of down to 1 μs is adopted.
The evolution process of the shock propagation and the bubble deformation is captured in a single test run as presented in Figure 4. At the beginning the incident shock moves into the converging part and the shock front bends to a circular arc (frame 0). When the cylindrical shock passes across the SF6 bubble, a reflected shock moves back and a refracted shock propagates inside the gas volume more slowly than the incident shock due to the larger acoustic impedance inside the bubble (frames 1 to 3). Then typical phenomena including the shock focusing and the jet formation appear (frame 4). In this process the refracted shock converges to a very small zone inside the SF6 bubble, resulting in very high pressure that promotes the interface deformation. As time proceeds, the SF6 jet is elongated and the interface body develops into a vortex ring (frames 5 to 9). Later on, as the shock reflects from the center of convergence and interacts with the evolving interface (frames 14 to 22), the evolution of the interface is intensified and the flow quickly resembles turbulent mixing.
Note that the phenomenon presented here is similar to that in the planar shock case (Zhai et al., Reference Zhai, Si, Luo and Yang2011; Si et al., Reference Si, Zhai, Luo and Yang2012), but a great difference between them exists. On one hand, the shock focusing coming from the cylindrical shock acceleration is much stronger because of the gradually enhanced shock strength, thus the SF6 jet moves more quickly for the same incident shock strength. On the other hand, the shock structure seems more complex in the converging situation. For example, the transmitted and refracted shock waves propagating toward the center of convergence impinge on the wedge walls and the reflected waves will impact with the distorted interface again. These additional waves will certainly complicate the interaction process.
The displacements and structural scales of the interface are further measured. Figure 5 gives the variations of the interface displacements for the left and the right interfaces (assuming the incident shock moves from left to right as inserted in the figure) and structural dimensions including the interface length and the vortex ring length. The positions of measurement always refer to the leftmost or rightmost portion of the interface. At the beginning, the left interface first obtains a velocity while the right interface is still stationary. Then the left interface keeps moving with nearly the same velocity and the right interface moves quickly because of the SF6 jet formation. In this stage, the velocities of the left and right interfaces are estimated to be about 67.7 m/s and 156.6 m/s, respectively. As the reflected shock arrives and interacts with the evolving interface, the SF6 jet is first compressed, obtains a sudden opposite velocity, and finally diffuses gradually. In the process, the left interface first stops moving and then moves back with a very small velocity (about 8.8 m/s). Due to the variation of the interface displacements, the structural dimensions of the interface also change with time, as shown in Figure 5b. After the incident shock passage, the interface length decreases because of the compression. As the SF6 jet and the vortex ring are formed, all interface structures increase in size. When the reflected shock arrives, the interface length and the vortex ring length start decreasing. At the late stage the intensified mixing between two fluids makes the overall length increase.
3.2. Gas Cylinder
The gas cylinder is one of the typical membraneless interfaces widely used in past RM instability experiments with respect to a planar shock wave (Jacobs, Reference Jacobs1993; Tomkins et al., Reference Tomkins, Kumar, Orlicz and Prestridge2008). In this work, the shock tube (ST-2) with M 0 = 1.2, θ0 = 15°, and h = 95 mm which are obtained theoretically (Si et al., Reference Si, Zhai, Luo and Yang2014) is used. The SF6 cylinder well mixed with glycol droplets is positioned within the converging part. The diameter D 0 of the gas cylinder and the distance L 0 from the gas cylinder axis to the center of convergence can be adjusted. It should be pointed out that besides the feature of the gas cylinder (e.g., D 0 and L 0), the initial conditions in experiments cover the shape of the test section (e.g., M 0, θ0, and h) and the category of the gas (e.g., SF6, krypton, argon, and others). In present work, we mainly focus on the influence of L 0 and more studies are the scope of our next work. Here the diameter D 0 = 5 mm and two different distances (L 0 = 93.3 mm and L 0 = 115.5 mm ) are adopted. The flow is monitored by the continuous laser-sheet imaging method, as sketched in Figure 3b. Different frame rates of the high-speed video camera (30,000 fps and 35,000 fps) are used for different temporal and spatial resolutions (the corresponding time intervals are 33.3 μs and 28.6 μs and spatial resolutions are 256 × 800 and 192 × 856, respectively). The shutter of the camera is down to 1 μs. The width of the laser-sheet is carefully adjusted in order to gain an acceptable scattering light intensity when the flow field is illuminated.
Figures 6 and 7 present the sequences of cross-sectional view of the SF6 cylinder accelerated by the cylindrical converging shock and the reflected shock from the center of convergence (assuming the incident shock moves from bottom to top as shown in the figure). Note that the initial conditions have been confirmed by the direct observation of the flow prior to each run. For L 0 = 93.3 mm, the shock Mach number is about 1.35 when the cylindrical shock impacts with the SF6 cylinder. At the beginning, the SF6 cylinder develops into a crescent shape quickly. As time proceeds, a counter-rotating vortex pair is formed due to the baroclinic vorticity initially deposited on the boundary of the gas cylinder. As the evolving interface undergoes the reflected shock, the vortex pair gathers most of the interface mass at its top boundaries and continues rotating in the same directions. The bottom arcuate interface becomes flat and then breaks up into two parts, moving together with the vortex pair. At late times, the vortex pair diffuses gradually due to the mixing of two gases. For L 0 = 115.5 mm, the shock Mach number is about 1.32 when the cylindrical shock arrives at the position of the SF6 cylinder. In the early stage, the evolution of the interface is similar to that described above. As the reflected shock arrives, the evolving interface suddenly stops moving onward and the top vortex pair rolls up with opposite rotating directions. Later on, the bottom circular interface grows in size and gradually diffuses, and the top vortex pair enhances the mixing of the SF6 gas with the ambient air.
It is interesting that there is a reversal of the rotating direction of the re-accelerated vortex pair when the value of L 0 increases. The previous experimental results for L 0 = 181.5 mm (Si et al., Reference Si, Zhai, Luo and Yang2014) also demonstrated this observation. The differences between these cases are mainly attributed to the Mach number effect, or rather, to the generation and distribution of baroclinic vorticity due to the misalignment between the density and pressure gradients. On one hand, the initial shock Mach numbers are different when the cylindrical shock impacts with the SF6 cylinder for different L 0. A larger L 0 corresponds to a smaller initial shock Mach number. On the other hand, although the initial gas cylinder diameter and the cylindrical shock shape are nearly the same for these cases, the evolving interfaces between them at the time when the reflected shock arrives are different due to the different time period of evolution. As a whole, the rotating direction of the vortex pair is dependent on the positive and negative values of the superposed baroclinic vorticity deposited on the interface.
Figure 8 presents the experimental data for both L 0 = 93.3 mm and L 0 = 115.5 mm showing the variation with time of the interface displacements for the left and the right interfaces (assuming the incident shock moves from left to right as inserted in the figure) and structural dimensions including the interface length and the vortex pair core distance. Before the arrival of the reflected shock, both the left and the right interfaces move toward the center of convergence and the velocities of the interface for L 0 = 93.3 mm are larger than those for L 0= 115.5 mm. The difference mainly comes from their different shock Mach numbers as the cylindrical shock wave interacts with the SF6 cylinder. In this period, the right interface moves more quickly than the left interface due to the development of the vortex pair. When the reflected shock propagates through the evolving bubble, each interface suddenly stops moving and then obtains a negative velocity. Although the rotating directions of the main vortex pair are opposite for L 0 = 93.3 mm and L 0 = 115.5 mm, the tendencies of the interface structures changing with time between them have the similar regularity except when the evolving interface for L 0 = 115.5 mm develops into a reversed vortex pair and a circular interface at the late stage. The overall lengths of the interfaces for both L 0 = 93.3 mm and L 0 = 115.5 mm are plotted in Figure 8b, where the vortex pair core distance is also given. Before the arrival of the reflected shock, the interface length first decreases because of the incident shock passage and then gradually increases after the generation of the vortex pair. In this process, the transmitted and refracted shock waves propagating toward the center of convergence reflect from the wedge walls and further compress the interface, making the vortex pair core distance be reduced. After the reflected shock passes through the distorted interface, different phenomena for L 0 = 93.3 mm and L 0 = 115.5 mm result in different tendencies of the structural dimensions changing with time. For L 0 = 93.3 mm, as the counter-rotating vortex pair develops still in the same direction, the overall interface length first decreases in a very short time period and then slightly increases, and the vortex pair core distance first keeps nearly constant and then gradually increases. Later on, both of the structural scales seldom change accompanied with the diffusion of the interface. For L 0 = 115.5 mm, the overall interface length always increases in the process except when the reflected shock just impacts with the evolving interface, and the vortex pair core distance seems decreasing all the time. One can find oscillations appearing within the experimental data, which may be caused by the interaction of the interface with the shock reflected from the wedge walls of the shock tube test section. As a whole, the complete process of the interface evolution can be revealed by the variations of the interface structures with time.
4. CONCLUSIONS
Experiments have been carried out to explore the potential of experimental method in study of RM instability induced by cylindrical converging shock waves. The test section of the shock tube was well designed based on shock dynamics theory to directly convert the incident planar shock into the cylindrical one. Spherical and cylindrical gaseous interfaces were employed in two sets of shock tubes with different control parameters. High-speed imaging combined alternately with two optical systems including schlieren and planar Mie scattering photography was utilized to capture the evolution of flow structures in a single test run for each case. Rich phenomena have been observed in different cases. The typical shock focusing and jet formation are found corresponding to the gas bubble. The counter-rotating vortex pair is formed in case of gas cylinder and a reversal is obtained when the position of the gas cylinder in the test section changes. The experiments also indicate that the phenomena of the interface accelerated by the converging shock and the reflected shock formed from the center of convergence are a little similar to those found in the planar shock cases, but differences also exist between them because of the shock curvature and Mach number effects in the converging cases, which can be attributed to the baroclinic mechanism due to the misalignment of density and pressure gradients. The experimental facility reported here provides a significant tool in study of RM instability with respectto converging shock waves. It must be mentioned that more experimental results and quantitative analysis are planned to be performed in the future. The effects of initial conditions on the growth rate of perturbations and the turbulent mixing of two gases should be studied thoroughly in order to go deep into the physical mechanism of interface evolution and compressible turbulence.
ACKNOWLEDGMENTS
The authors would like to thank Xiansheng Wang and Fu Zhang for their help during the experiments. The work was supported by the National Natural Science Foundation of China (11272308 and 11302219), the Knowledge Innovation Program of the Chinese Academy of Sciences (GY2011053006), and the Fundamental Research Funds for Central Universities (WK2090050020).