An official website of the United States government
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.
more »
« less
Interface models for three-dimensional Rayleigh–Taylor instability
We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.
more »
« less
- Award ID(s):
- 2007606
- PAR ID:
- 10451456
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 959
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This work numerically investigates the role of viscosity and resistivity in Rayleigh–Taylor instabilities in magnetized high-energy-density (HED) plasmas for a high Atwood number and high plasma beta regimes surveying across plasma beta and magnetic Prandtl numbers. The numerical simulations are performed using the visco-resistive magnetohydrodynamic equations. Results presented here show that the inclusion of self-consistent viscosity and resistivity in the system drastically changes the growth of the Rayleigh–Taylor instability (RTI) as well as modifies its internal structure at smaller scales. It is seen here that the viscosity has a stabilizing effect on the RTI. Moreover, the viscosity inhibits the development of small-scale structures and also modifies the morphology of the tip of the RTI spikes. On the other hand, the resistivity reduces the magnetic field stabilization, supporting the development of small-scale structures. The morphology of the RTI spikes is seen to be unaffected by the presence of resistivity in the system. An additional novelty of this work is in the disparate viscosity and resistivity profiles that may exist in HED plasmas and their impact on RTI growth, morphology and the resulting turbulence spectra. Furthermore, this work shows that the dynamics of the magnetic field is independent of viscosity and likewise the resistivity does not affect the dissipation of enstrophy and kinetic energy. In addition, power law scalings of enstrophy, kinetic energy and magnetic field energy are provided in both the injection range and inertial sub-range, which could be useful for understanding RTI induced turbulent mixing in HED laboratory and astrophysical plasmas and could aid in the interpretation of observations of RTI-induced turbulence spectra.more » « less
-
The interfacial instability in a two-phase mixing layers between parallel gas and liquid streams is important to two-phase atomization. Depending on the inflow conditions and fluid properties, interfacial instability can be convective or absolute. The goal of the present study is to investigate the impact of gas viscosity on the interfacial instability. Both interface-resolved simulations and linear stability analysis (LSA) have been conducted. In LSA, the Orr–Sommerfeld equation is solved to analyze the spatio-temporal viscous modes. When the gas viscosity decreases, the Reynold number (Re) increases accordingly. The LSA demonstrates that when Re is higher than a critical threshold, the instability transitions from the absolute to the convective (A/C) regimes. Such a Re-induced A/C transition is also observed in the numerical simulations, though the critical Re observed in simulations is significantly lower than that predicted by LSA. The LSA results indicate that the temporal growth rate decreases with Re. When the growth rate reaches zero, the A/C transition will occur. The Re-induced A/C transition is observed in both confined and unconfined mixing layers and also in cases with low and high gas-to-liquid density ratios. In the transition from typical absolute and convective regimes, a weak absolute regime is identified in the simulations, for which the spectrograms show both the absolute and convective modes. The dominant frequency in the weak absolute regime can be influenced by the perturbation introduced at the inlet. The simulation results also show that the wave propagation speed can vary in space. In the absolute instability regime, the wave propagation speed agrees well with the absolute mode celerity near the inlet and increases to the Dimotakis speed further downstream.more » « less
-
Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.more » « less
-
Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.98
