Atwood Number Research Articles

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

The Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)] are extended to construct models for reshocked Richtmyer–Meshkov mixing. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

Rayleigh-Taylor instability in magnetohydrodynamics with finite resistivity in a horizontalmagnetic field.

Recent studies have revealed the significant influence of finite resistivity on high-energy-density plasmas, contrary to the previous findings of Jukes [J. Fluid Mech. 16, 177 (1963)0022-112010.1017/S0022112063000677]. This paper reexamines Jukes' theory in the context of magneto-Rayleigh-Taylor instability in magnetohydrodynamics with finite resistivity represented by η. The inadequacy of Jukes' approach due to an erroneous boundary condition is demonstrated, and it is shown that although the theory provides some physical insights, it fails to capture crucial features. The dispersion relation proposed in this study highlights that larger growth rates tend to diffuse the magnetic field rapidly, negating its suppressive effect. Moreover, the Atwood number has a significant influence on the growth-rate curves' shape, which differs from those of viscous or elastic flows and ideal magnetohydrodynamics. Additionally, long wavelengths grow proportionally to η^{1/3}, while α indicating growth rates behaves classically when the magnetic field is entirely diffused.

Experimental study on Richtmyer–Meshkov instability at a light–heavy interface over a wide range of Atwood numbers

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).

Inferring the Magnetic Field from the Rayleigh–Taylor Instability

By extending the potential approach of Layzer and Goncharov to the magnetohydrodynamics equations, we find the nonlinear solutions to the single-mode Rayleigh–Taylor instability subjected to uniform magnetic fields at various inclinations. This allows us to derive the analytical prediction of the terminal bubble and spike velocities at arbitrary Atwood numbers, which are assessed against various 2D and 3D direct numerical simulations. Contrary to the linear phase, where the magnetic field inhibits or delays the instability, the growth rate may be enhanced in the nonlinear regime, exhibiting velocities faster than the Alfvén speed. This sheds light on the importance of the nondimensional number expressing the competition between the magnetic and buoyancy effects. Conversely, we show how the orientation and the intensity of the magnetic field can be simply inferred from these solutions.

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux$a_y$is calculated. In addition to mostly negative values of$a_y$, typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of$a_y$due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ($\sim$60 %) is observed.

Large-eddy simulations of compressible Rayleigh–Taylor turbulence with miscible fluids using spatial gradient model

Large-eddy simulations (LES) and implicit LES (ILES) of three-dimensional compressible Rayleigh–Taylor turbulence with miscible fluids are performed and compared with direct numerical simulation (DNS) at the Atwood number At=0.5 and stratification parameters Sr = 1.0 and 4.0. Three sub-grid-scale (SGS) models including constant-coefficient spatial gradient model (CSGM), dynamic Smagorinsky model (DSM), and dynamic mixed model (DMM) are considered. The CSGM model achieves a high accuracy by using the velocity gradients at the neighboring LES grids. The a priori tests show that the CSGM model has significantly higher correlation coefficients and lower relative errors than traditional SGS models. In the a posteriori tests, the probability density functions of SGS terms predicted by the CSGM model are consistent with the filtered DNS results. The CSGM model can accurately predict the small bubble and spike structures, resulting in good predictions of mixing heights and concentration fields. The instantaneous structures, spectra, and statistics of velocity and vorticity fields are also examined, showing the excellent performance of the CSGM model compared to the ILES, DSM, and DMM models. Moreover, the predictions of the temperature and pressure fields by the CSGM model are significantly better than the traditional SGS models and ILES.

Multi-layer Rayleigh–Taylor instability: Consequences for naturally occurring stratified mixing layers

The present study investigates the behavior of multi-layer Rayleigh–Taylor instability (RTI) and enstrophy transport in the flow using a three-dimensional computational framework. The dynamics of RTI are explored in a monotonically unstable stratified fluid system composed of air at different constant temperatures, initially separated by insulating partitions. Our results illustrate the formation of a multi-layer RTI system and the growth of convective mixing layers from interfaces between consecutive layers. The behavior of single-layer and multi-layer RTI is compared by considering the influence of Atwood numbers on growth and characteristics of the mixing layers. We found that the presence of multiple layers affects the onset and development of RTI. The merging of top layers leads to accelerated mixing layer growth, while the bottommost layer experiences early-stage RTI. Furthermore, we utilize the compressible enstrophy transport equation to characterize dominant mechanisms controlling the spatiotemporal evolution of the multi-layer RTI. The results highlight the significance of viscous and compressibility terms, especially in the center. In contrast, the vortex stretching term dominates near side walls for later stages of the RTI. These insights enrich the understanding of multi-layer RTI and its effects on enstrophy transport and mixing in unstably stratified fluid systems, providing valuable data for validating numerical methods and informing future research in related fields.

Reduced mixing in inertial confinement fusion with early-time interface acceleration.

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. This characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

A Level Set-Based Solver for Two-Phase Incompressible Flows: Extension to Magnetic Fluids

Development of a two-phase incompressible solver for magnetic flows in the magnetostatic case is presented. The proposed numerical toolkit couples the Navier–Stokes equations of hydrodynamics with Maxwell's equations of electromagnetism to model the behaviour of magnetic flows in the presence of a magnetic field. To this end, a rigorous implementation of a second-order two-phase solver for incompressible nonmagnetic flows is introduced first. This solver is implemented in the finite-difference framework, where a fifth-order conservative level set method is employed to capture the evolution of the interface, along with an incompressible solver based on the projection scheme to model the fluids. The solver demonstrates excellent performance even with high density ratios across the interface (Atwood number ≈ 1 ), while effectively preserving the mass conservation property. Subsequently, the numerical discretisation of Maxwell's equations under the magnetostatic assumption is described in detail, utilising the vector potential formulation. The primary second-order solver for two-phase flows is extended to the case of magnetic flows, by incorporating the Lorentz force into the momentum equation, accounting for high magnetic permeability ratios across the interface. The implemented solver is then utilised for examining the deformation of ferrofluid droplets in both quiescent and shear flow regimes across various susceptibility values of the droplets. The results suggest that increasing the susceptibility value of the ferrofluid droplet can affect its deformation and rotation in low capillary regimes. In higher capillary flows, increasing the magnetic permeability jump across the interface can further lead to droplet breakup as well. The effect of this property is also investigated for the Rayleigh–Taylor instability growth in magnetic fluids.

Investigation of strong isothermal stratification effects on multi-mode compressible Rayleigh–Taylor instability

Rayleigh–Taylor instability, RTI, occurs at the interface separating two fluids subjected to acceleration when the density gradient and the acceleration are in opposite directions. Previous scientific research primarily considered RTI under the incompressible assumption, which may not be valid in many high-energy-density engineering applications and astrophysical phenomena. In this study, the compressibility effects of the background isothermal stratification strength on multi-mode two-dimensional RTI are explored using fully compressible multi-species direct numerical simulations. Cases under three different isothermal Mach numbers – Ma=0.15, 0.3, and 0.45 – are investigated to explore weakly, moderately, and strongly stratified compressible RTI, respectively, at an Atwood number of 0.04. Unlike incompressible RTI, an increase in the flow compressibility through the strength of the background stratification can suppress the RTI growth and can lead to a termination of the RTI mixing layer growth with a highly molecularly mixed state. Our findings suggest that even at the chosen relatively low Atwood number, the variable-density effects can be significantly enhanced due to an increase in the background stratification for the compressible RTI as different spatial profiles become noticeably asymmetric across the mixing layer for the strongly stratified case. In addition, this study compares the chaotic behavior of the cases by studying the transport of the turbulent kinetic energy as well as the vortex dynamics. The Reynolds number dependence of the results is also examined with three different Reynolds numbers, and the findings for the large-scale mixing and flow quantities of interest are shown to be universal in the range of the Reynolds numbers studied.

Experimental and Numerical Investigation of the Dynamics of Development of Rayleigh–Taylor Instability at Atwood Numbers Close to Unity

Experimental and Numerical Investigation of the Dynamics of Development of Rayleigh–Taylor Instability at Atwood Numbers Close to Unity

Numerical study of the effect of a magnetic field on Rayleigh-Taylor instability with different density disturbances

Rayleigh-Taylor instability (RTI) is a fundamental physical process in fluids and plasmas. RTI is ubiquitous and must be considered in the field of high-energy-density physics, such as in space physics, astrophysics, and inertial confinement fusion. In this work, two-dimensional numerical simulations of laser-driven RTI with different density perturbations are performed using a radiation magnetohydrodynamic simulation program (FLASH). The effect of the applied magnetic field on the evolution of RTI at different Atwood numbers is systematically discussed. The results show that RTI evolves freely without an external magnetic field, and it is accompanied by the generation of secondary Kelvin–Helmholtz instability. Reducing the Atwood number weakens the mixing of fluids and has a strong stabilizing effect on the RTI. Introducing an external magnetic field parallel to the perturbation wave vector further inhibits the development of RTI and Kelvin–Helmholtz instability, with magnetic pressure playing a dominant role. The study results are important to gaining an in-depth understanding of the mixing of magnetic fluids and the magnetic field evolution at the instability interface and provide a reference for subsequent experimental studies on the related magnetization RTI.

Impact dynamics of compound drops of fluids with density contrast

The dynamics of compound drops impacting on a flat substrate is numerically investigated using a ternary-fluid diffuse-interface method, with the aim of assessing the effect of a density difference between the inner and outer droplets (denoted by $\lambda$ ) on the evolution of the interfaces. With the help of numerical simulations, we find that, at the intermediate stage of drop impact, the inner droplet exhibits a self-similar deformation at $\lambda =1$ and relatively high Weber number, and experiences more or less a uniform acceleration for various $\lambda$ . In particular, the acceleration magnitude at $\lambda \ne 1$ can be correlated with the acceleration at $\lambda =1$ and the Atwood number. When the inner droplet is denser than the outer one, a lamella occurs at the spreading front of the inner droplet. We present a scaling analysis of the thickness of the lamella, and the resultant theoretical prediction is in good agreement with numerical results. At the maximal spreading of the compound drop, a bulging structure is formed around the symmetry axis due to the presence of the inner droplet, thereby effectively reducing the liquid supply to the spreading front and leading to a decrease of maximal spreading ratio $\beta _{max}$ as compared with a pure drop. We proposed a corrected Weber number $We^*_\lambda$ by taking account of the combined effects of $\lambda$ , volume fraction of the inner droplet, Weber number and morphology of the compound drop. Integrating $We^*_\lambda$ with the universal model of $\beta _{max}$ for impacting pure drops, we successfully build up a new model for predicting the maximal spreading ratio of impacting compound drops with various $\lambda$ .

Comparison of exchange flow of two miscible fluids in curved and slanted straight tubes

The front velocity in miscible exchange flow along the straight and curved tubes has been investigated. Almost the entire flow lies within the range of inertia dominated regime in the current experimental setup. Available data for straight pipes at different Atwood numbers have been revisited to provide an appropriate ground for better understanding and evaluation of the observations made for the flow in curved tubes. A 19.05 mm diameter and 2.40 m long tube has been formed to follow 55°of an arc in the vertical plane with the radius of curvature of 2.5 meters. The upper and lower parts are separated by a gate valve and filled with heavier and lighter fluids, respectively. The observations and analyses demonstrate how the continuous variations of the flow inclination along the curved path and the change in mixing dynamics associated with it can affect the velocity of downward and upward moving fronts. Most notably it has been detected that the velocity of upward moving front follows a similar trend as in straight pipes due to the continuously increasing mixing. However, the velocity of downward moving heavy front slightly decreases as opposed to the case of straight pipes, signifying the effect of curvature on mixing enhancement.

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

Numerical study of Richtmyer–Meshkov instability of a flat interface driven by perturbed and reflected shock waves

In this paper, the Richtmyer–Meshkov instability of a flat gas interface driven by perturbed and reflected shock waves is numerically investigated. The flat gas interface evolves into a “Λ”-shaped structure with a central N2 cavity and steps on both sides, due to the impaction of the perturbed shock wave. After the secondary collision of the reflected shock wave from the high-density region to the low-density region, the gas interface first undergoes phase inversion, and the “Λ” interface then evolves into a bubble and spike structure. Three cases of different Atwood numbers, N2/SF6, N2/Kr, and N2/CO2, are studied. The collision time and position of the reflected shock wave and the interface, the induced spikes, bubbles and gas mixing, are compared in detail. The formation of the spike and bubble is related to the RM instability developed by the collision of the reflected shock wave and the perturbed interface, in which the effect of baroclinic vorticity is highlighted. With the increase in the Atwood number, the density gradient and the baroclinic vorticity become larger, which induces more vortex along the interface. Kelvin Helmholtz unstable vortices are generated on the “legs” of the spikes due to shearing. The main spike structure is stretched and broken with the effect of the vortex, forming a turbulent mixing zone.

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion [Kuranz et al., Nat. Commun. 9, 1564 (2018)]. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. A terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.