Effect of Atwood number on convergent Richtmyer–Meshkov instability

Author(s):  
Jinggang Tang ◽  
Fu Zhang ◽  
Xisheng Luo ◽  
Zhigang Zhai
Keyword(s):  
2017 ◽  
Author(s):  
◽  
Constantine Gregory Avgoustopoulos

This paper investigates the experimental work in Shock Driven Multiphase Instabilities (SDMI). SDMIs occur when an interface consisting of a particle seeded gas is instantaneously accelerated and begins mixing. SDMIs have similar flow morphologies to the Richtmyer-Meshkov Instability (RMI), however, the driving force inducing this flow is very different. SDMIs occur when there is a relative velocity difference between surrounding gas and the moving particles. This results to a shear at the edges and ultimately leads to rollups that are similar to a RMI. To investigate this phenomena, a shock tube facility was designed, calibrated, and tested to perform experiments. The experimental data was qualitatively compared to simulations performed, as well as to literature of similar experiments. Quantitative data was analyzed using Particle Imaging Velocimetry (PIV) to understand the flow of the instability. The flow morphologies observed in experiments have similar behavior to those performed in simulations. Additionally, the qualitative observations of experiments performed in this lab are also in agreement with experimental literature. Two different effective Atwood numbers are investigated in this study. The first case looks at a gas cylinder interface with an effective Atwood number of -0.01 and a gas Atwood number of -0.02, shocked with a Mach 1.66 shock wave. The observations show a dominating instability resulting in the gas Atwood number. What ends up happening is the smaller particles are pulled into the vortex and the large particles separate and trail behind. The second case looks at the same gas cylinder perturbation but with an effective Atwood number of 0.03 and a gas Atwood number of 0, shocked at Mach 1.66. The higher Atwood number was achieved by modifying the experimental apparatus slightly to deliver a greater number of particles to the shock tube. The experiments observed show that there is agreement with literature and simulations. Certain unusual filaments begin forming at late times, 4.0ms after shock. This was thought to only appear in a pure RMI. In the case of a SDMI, these filaments are a result of colliding particles.


2021 ◽  
Vol 928 ◽  
Author(s):  
Xinliang Li ◽  
Yaowei Fu ◽  
Changping Yu ◽  
Li Li

In this paper, the Richtmyer–Meshkov instabilities in spherical and cylindrical converging geometries with a Mach number of approximately 1.5 are investigated by using the high resolution implicit large eddy simulation method, and the influence of the geometric effect on the turbulent mixing is investigated. The heavy fluid is sulphur hexafluoride (SF6), and the light fluid is nitrogen (N2). The shock wave converges from the heavy fluid into the light fluid. The Atwood number is 0.678. The total structured and uniform Cartesian grid node number in the main computational domain is 20483. In addition, to avoid the influence of boundary reflection, a sufficiently long sponge layer with 50 non-uniform coarse grids is added for each non-periodic boundary. Present numerical simulations have high and nonlinear initial perturbation levels, which rapidly lead to turbulent mixing in the mixing layers. Firstly, some physical-variable mean profiles, including mass fraction, Taylor Reynolds number, turbulent kinetic energy, enstrophy and helicity, are provided. Second, the mixing characteristics in the spherical and cylindrical turbulent mixing layers are investigated, such as molecular mixing fraction, efficiency Atwood number, turbulent mass-flux velocity and density self-correlation. Then, Reynolds stress and anisotropy are also investigated. Finally, the radial velocity, velocity divergence and enstrophy in the spherical and cylindrical turbulent mixing layers are studied using the method of conditional statistical analysis. Present numerical results show that the geometric effect has a great influence on the converging Richtmyer–Meshkov instability mixing layers.


Shock Waves ◽  
2012 ◽  
Vol 22 (4) ◽  
pp. 317-325 ◽  
Author(s):  
C. Weber ◽  
N. Haehn ◽  
J. Oakley ◽  
M. Anderson ◽  
R. Bonazza
Keyword(s):  

2016 ◽  
Vol 795 ◽  
pp. 313-355 ◽  
Author(s):  
Bhanesh Akula ◽  
Devesh Ranjan

Simultaneous density and velocity turbulence statistics for Rayleigh–Taylor-driven flows at a moderately high Atwood number ($A_{t}$) of $0.73\pm 0.02$ are obtained using a new convective type or statistically steady gas tunnel facility. Air and air–helium mixture are used as working fluids to create a density difference in this facility, with a thin splitter plate separating the two streams flowing parallel to each other at the same velocity ($U=3~\text{m}~\text{s}^{-1}$). At the end of the splitter plate, the two miscible fluids are allowed to mix and the instability develops. Visualization and Mie-scattering techniques are used to obtain structure shape, volume fraction profile and mixing height growth information. Particle image velocimetry (PIV) and hot-wire techniques are used to measure planar and point-wise velocity statistics in the developing mixing layer. Asymmetry is evident in the flow field from the Mie-scattering images, with the spike side showing a more gradual decline in volume fraction than the bubble side. The spike side of the mixing layer grows 50 % faster than the bubble side. PIV is implemented for the first time in these moderately high-Atwood-number experiments ($A_{t}>0.1$) to obtain root-mean-square velocities, anisotropy tensor components and Reynolds stresses across the mixing layer. Overall, the turbulence statistics measured have shown different scaling compared to small-Atwood-number experiments. However, the total probability density functions for the velocities and turbulent mass fluxes exhibit behaviour similar to small-Atwood-number experiments. Conditional statistics reveal different values for turbulence statistics for spikes and bubbles, unlike small-Atwood-number experiments.


2001 ◽  
Author(s):  
O. Sadot ◽  
A. Yosef-Hai ◽  
Dan Oron ◽  
A. Rikanati ◽  
D. Kartoon ◽  
...  

Author(s):  
Arindam Banerjee

The release of a mass of hydrogen fuel (gas) into the ambient atmosphere results in the transient formation of flammable mixture zones that represent potential fire, explosion and toxic hazards. The formation of mixing zones of air and hydrogen for this simple geometry follows the classical Rayleigh Taylor (R-T) instability, which is induced when a heavy fluid is placed over a light fluid in a gravitational field. Buoyancy driven mixing in such flow configurations is studied by using the Boussinesq approximation and considering the flow to be laminar. However, this approximation is valid only at low Atwood numbers (non-dimensional density differences). As Atwood number increases (>0.1, i.e. large density differences) the Boussinesq approximation is no longer valid and a distinct bubble and spike geometry of Rayleigh-Taylor buoyant plumes is formed. Aside from asymmetry in the flow the Atwood number also affects key parameters such as the growth constants and molecular mix. The effect of initial conditions on the growth rate of turbulent Rayleigh-Taylor (RT) mixing has been studied using carefully formulated numerical simulations. A monotone integrated large-eddy simulation (MILES) using a finite-volume technique was employed to solve the three-dimensional incompressible Euler equations with numerical dissipation for air and hydrogen mixing at Atwood number 0.875. The study also provides preliminary guidelines for reducing the fire and explosion hazards in enclosures where such situations are present.


Author(s):  
Arindam Banerjee ◽  
Malcolm J. Andrews

A novel gas channel experiment is used to study the non-equilibrium development of high Atwood number Rayleigh-Taylor mixing. Two gas streams, one containing air-helium mixture and the other air, flow parallel to each other separated by a thin splitter plate. The streams meet at the end of a splitter plate leading to the formation of an unstable interface and initiation of buoyancy driven mixing. This set up is statistically steady and allows for long data collection times. Here, we describe initial measurements to determine the density profile and growth rate along the mix at low density differences (At ~ 0.05). The facility is however designed capable of large Atwood number studies (At ~ 0.75). Diagnostics include high resolution digital image analysis, which is used to determine the density profile across the mix. The growth parameter (α) is also estimated by a “moving window” calculation. The results are then verified with measurements of α made by a Constant temperature (CT) hot-wire probe and with the growth parameter obtained from small Atwood number experiments (At ~ 0.001). However, there were some inherent errors in the density profile measurements because of non-uniformity in the concentration of smoke. To verify that these errors were indeed measurement errors and not as a result of lack of statistical convergence, a detailed statistical convergence test was performed. It showed that convergence was a direct consequence of the number of different large 3D structures that were averaged over the duration of the run.


2018 ◽  
Vol 838 ◽  
pp. 320-355 ◽  
Author(s):  
R. V. Morgan ◽  
W. H. Cabot ◽  
J. A. Greenough ◽  
J. W. Jacobs

Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Single-mode two-dimensional, and single-mode three-dimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a non-constant acceleration, and a time decreasing Atwood number, $A=(\unicode[STIX]{x1D70C}_{2}-\unicode[STIX]{x1D70C}_{1})/(\unicode[STIX]{x1D70C}_{2}+\unicode[STIX]{x1D70C}_{1})$, where $\unicode[STIX]{x1D70C}_{2}$ and $\unicode[STIX]{x1D70C}_{1}$ are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of $A=0.49$, $A=0.63$, $A=0.82$ and $A=0.94$. Nominally two-dimensional (2-D) experiments (initiated with nearly 2-D perturbations) and 2-D simulations are observed to approach an intermediate-time velocity plateau that is in disagreement with the late-time velocity obtained from the incompressible model of Goncharov (Phys. Rev. Lett., vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2-D bubbles in large wavenumber, $k=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D706}=0.247~\text{mm}^{-1}$, experiments and simulations, where $\unicode[STIX]{x1D706}$ is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These late-time velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian (Phys. Fluids., vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a free-fall like behaviour. Finally, experiments initiated with three-dimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.


Author(s):  
Mohammadreza Yavari ◽  
Elaheh Bagherizadeh ◽  
Majid Bazargan

Miscible displacement flow in a curved pipe with two fluids of equal viscosity and in the regime of low Atwood number is studied experimentally at different Reynolds and Atwood numbers. By using a curved pipe the effect of curvature on miscible displacement flow is studied. Curvature may be present in many displacement flow processes in nature or industry which underlines the necessity of studying its effect on displacement flow. As the flow is controlled by gravity as the main driving force, only low to moderate Reynolds numbers are considered and Atwood number is varied by adding NaCl salt to fresh water. This makes it possible to create different values of Atwood numbers which can be varied continuously. The position of the leading front is carefully recorded using a digital camera and is used as a measure of displacement efficiency. It is observed that the ratio of the front velocity to the mean velocity approaches a certain value as Reynolds number increases. The effect of Atwood number on flow dynamics is also studied based on experimental results and is interpreted following the conventions employed in some of the previous researches.


Sign in / Sign up

Export Citation Format

Share Document