Multi-fluid mixing in a system with three stratified layers is explored using two-dimensional compressible direct numerical simulations (DNS) by solving fully compressible multi-species Navier-Stokes equations. All configuration cases under investigation consist of at least one acceleration-driven Rayleigh-Taylor unstable interface where the direction of the acceleration is from the heavier to the lighter fluids to initiate the chaotic mixing. The DNS are initialized with isopycnic background stratification, where the species densities of the three fluids are initially constant over the domain, and also they are benefited from the adaptive mesh refinement (AMR) to reduce the computational cost. We investigate four cases with the global Atwood number of 0.04 where the stratified layers along the acceleration direction (from bottom to top) have heavy-intermediate-light, heavy-light-light, heavy-light-intermediate, and heavy-light-heavy densities, and each case is studied with two different distances between the first and third fluid layers.
It is found that the globally unstable case with heavy-intermediate-light densities, which has two unstable interfaces, exhibits symmetric mixing during the flow evolution and eventually behaves similarly to the classical two-layer Rayleigh-Taylor instability (RTI) flows. This finding is consistent with the observation in previous three-layer RTI experiments. The case with heavy-light-light densities is performed for comparison purposes, and it practically only has a single-stratified interface similar to the two-layer RTI problem. For the remaining two cases, the top interface of the three-layer RTI is stably stratified. It is shown that for the cases where the top interface is neutrally stratified (e.g., heavy-light-light case) or weakly stratified (e.g., heavy-light-intermediate case), upward pure fluids penetration is larger compared to the heavy-light-heavy case, whose top interface is initially strongly stable. In addition to the study on large-scale RTI entrainment, we also present mixedness and vortical dynamics of the flows.
more »
« less
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.
more » « less- Award ID(s):
- 2234415
- NSF-PAR ID:
- 10477451
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 35
- Issue:
- 8
- ISSN:
- 1070-6631
- Page Range / eLocation ID:
- 084116
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Direct Numerical Simulation (DNS) of spatially-developing turbulent boundary layers (SDTBL) is performed over isothermal/adiabatic flat plates for incompressible and compressible-subsonic (M∞ = 0.5 and 0.8) flow regimes. Similar low Reynolds numbers are considered in all cases with the purpose of assessing modest flow compressibility on low/high order flow statistics of Zero Pressure Gradient (ZPG) flows. The considered molecular Prandtl number is 0.72. Additionally, temperature is regarded as a passive scalar in the incompressible SDTBL with the purpose to examine differences in the thermal transport phenomena of subsonic flows, i.e., passive vs. active scalar. It was found that the Van Driest transform and Morkovin scaling are able to collapse incompressible and subsonic quantities very well.more » « less
-
The current work evaluates the effectiveness of two leading-edge dynamic stall criteria in mild to moderately compressible regimes using numerical simulations. The two criteria under consideration, namely, the maximum magnitudes of the leading edge suction parameter (max(𝐿𝐸 𝑆𝑃)) and boundary enstrophy flux (max(|𝐵𝐸𝐹|)), have previously been found to be effective at signaling dynamic stall in the incompressible regime. Based on unsteady Reynolds-averaged Navier-Stokes simulations at a Reynolds number of 2 × 105 and freestream Mach numbers between 0.1 - 0.5, we observe that these criteria are directly applicable in the mild to moderately compressible regimes, since they are reached shortly after suction collapse at the leading edge and well in advance of dynamic stall vortex formation for all the cases. This is attributed to compressibility effects promoting adverse-pressure-gradient(APG)-induced stall for the flow conditions considered. For the highest Mach number of 0.5, shock wave interactions with the separated shear layer are observed. It is noted that although compressibility leads to separation at a lower APG, the maximum APG scaled by the local flow density remains in the same range for all the cases.more » « less
-
The asymptotic behaviour of Reynolds stresses close to walls is well established in incompressible flows owing to the constraint imposed by the solenoidal nature of the velocity field. For compressible flows, thus, one may expect a different asymptotic behaviour, which has indeed been noted in the literature. However, the transition from incompressible to compressible scaling, as well as the limiting behaviour for the latter, is largely unknown. Thus, we investigate the effects of compressibility on the near-wall, asymptotic behaviour of turbulent fluxes using a large direct numerical simulation (DNS) database of turbulent channel flow at higher than usual wall-normal resolutions. We vary the Mach number at a constant friction Reynolds number to directly assess compressibility effects. We observe that the near-wall asymptotic behaviour for compressible turbulent flow is different from the corresponding incompressible flow even if the mean density variations are taken into account and semi-local scalings are used. For Mach numbers near the incompressible regimes, the near-wall asymptotic behaviour follows the well-known theoretical behaviour. When the Mach number is increased, turbulent fluxes containing wall-normal components show a decrease in the slope owing to increased dilatation effects. We observe that $R_{vv}$ approaches its high-Mach-number asymptote at a lower Mach number than that required for the other fluxes. We also introduce a transition distance from the wall at which turbulent fluxes exhibit a change in scaling exponents. Implications for wall models are briefly presented.more » « less
-
Direct Numerical Simulation (DNS) of compressible spatially-developing turbulent boundary layers (SDTBL) is performed at a Mach number of 2.5 and low/high Reynolds numbers over isothermal Zero-Pressure Gradient (ZPG) flat plates. Turbulent inflow information is generated via a dynamic rescaling-recycling approach (J. Fluid Mech., 670, pp. 581-605, 2011), which avoids the use of empirical correlations in the computation of inlet turbulent scales. The range of the low Reynolds number case is approximately 400-800, based on the momentum thickness, freestream velocity and wall viscosity. DNS at higher Reynolds numbers (~3,000, about four-fold larger) is also carried out with the purpose of analyzing the effect of Reynolds number on the transport phenomena in the supersonic regime. Additionally, low/high order flow statistics are compared with DNS of an incompressible isothermal ZPG boundary layer at similar low Reynolds numbers and the temperature regarded as a passive scalar. Peaks of turbulence intensities move closer to the wall as the Reynolds number increases in the supersonic flat plate. Furthermore, Reynolds shear stresses depict a much larger "plateau" (constant shear layer) at the highest Reynolds number considered in present study.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0164504
