Y. Young,H. Tufo,A. Dubey,R. Rosner
Abstract:We investigate the miscible Rayleigh-Taylor (RT) instability in both 2 and 3 dimensions using direct numerical simulations, where the working fluid is assumed incompressible under the Boussinesq approximation. We first consider the case of randomly perturbed interfaces. With a variety of diagnostics, we develop a physical picture for the detailed temporal development of the mixed layer: We identify three distinct evolutionary phases in the development of the mixed layer, which can be related to detailed variations in the growth of the mixing zone. Our analysis provides an explanation for the observed differences between two and three-dimensional RT instability; the analysis also leads us to concentrate on the RT models which (1) work equally well for both laminar and turbulent flows, and (2) do not depend on turbulent scaling within the mixing layer between fluids. These candidate RT models are based on point sources within bubbles (or plumes) and interaction with each other (or the background flow). With this motivation, we examine the evolution of single plumes, and relate our numerical results (of single plumes) to a simple analytical model for plume evolution.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the physical mechanism and evolution characteristics of miscible Rayleigh - Taylor instability, especially the differences under two - dimensional (2D) and three - dimensional (3D) conditions. Specifically, the authors studied the Rayleigh - Taylor instability of incompressible fluids under the Boussinesq approximation through direct numerical simulation and explored the following issues:
1. **Mixing layer development in different dimensions**: The paper compared the development processes of the mixing layer under 2D and 3D conditions, identified three different evolution stages, and analyzed the relationships between these stages and the specific changes in the growth of the mixing zone.
- **Formula representation**: The relationship between the penetration depth \(h_{b, s}\) of the mixing zone, time \(t\), gravitational acceleration \(g\) and Atwood number \(A\) is:
\[
h_{b, s}=\alpha_{b, s} g A t^{2}
\]
where \(\alpha\) is a coefficient for measuring the potential energy release efficiency.
2. **Influence of initial conditions**: It was studied whether the initial conditions have an influence on the long - term evolution of the nonlinear Rayleigh - Taylor instability, especially whether a single scaling law can be expected to be valid throughout the entire process of the nonlinear Rayleigh - Taylor development.
3. **Applicability of different models**: To find models applicable to laminar and turbulent Rayleigh - Taylor instabilities, and these models should not depend on the turbulence scaling within the mixing layer.
4. **Evolution of a single plume**: The evolution of a single plume was studied through numerical simulation and compared with a simple analytical model to understand more complex multi - modal perturbation situations.
### Main contributions of the paper
- **Quantitative comparison**: Through detailed numerical simulations, a quantitative comparison of Rayleigh - Taylor instabilities under 2D and 3D conditions was provided, revealing the differences in the growth rates of the mixing layer widths in different dimensions.
- **Sensitivity to initial conditions**: The importance of initial conditions for the nonlinear evolution was emphasized, and it was pointed out that even over a long period of time, the memory effect of the system still exists and will not be completely erased by the turbulent motion within the mixing interface.
- **Model verification**: A Rayleigh - Taylor model based on point sources was proposed, and its effectiveness in describing the evolution of a single plume was verified through numerical simulation.
In summary, this paper aims to deeply understand the physical mechanism of Rayleigh - Taylor instability, especially the evolution characteristics of the mixing layer under different dimensions and initial conditions, providing an important reference for future experimental and theoretical research.