Transition stages of Rayleigh–Taylor instability between miscible fluids

Direct numerical simulations (DNS) are presented of three-dimensional, Rayleigh–Taylor instability (RTI) between two incompressible, miscible fluids, with a 3:1 density ratio. Periodic boundary conditions are imposed in the horizontal directions of a rectangular domain, with no-slip top and bottom walls. Solutions are obtained for the Navier–Stokes equations, augmented by a species transport-diffusion equation, with various initial perturbations. The DNS achieved outer-scale Reynolds numbers, based on mixing-zone height and its rate of growth, in excess of 3000. Initial growth is diffusive and independent of the initial perturbations. The onset of nonlinear growth is not predicted by available linear-stability theory. Following the diffusive-growth stage, growth rates are found to depend on the initial perturbations, up to the end of the simulations. Mixing is found to be even more sensitive to initial conditions than growth rates. Taylor microscales and Reynolds numbers are anisotropic throughout the simulations. Improved collapse of many statistics is achieved if the height of the mixing zone, rather than time, is used as the scaling or progress variable. Mixing has dynamical consequences for this flow, since it is driven by the action of the imposed acceleration field on local density differences.

[1]  D. Sharp.  An overview of Rayleigh-Taylor instability☆ , 1984 .

[2]  E. Meshkov.  Instability of the interface of two gases accelerated by a shock wave , 1969 .

[3]  J. Lumley,et al.  A First Course in Turbulence , 1972 .

[4]  S. Chandrasekhar.  Hydrodynamic and Hydromagnetic Stability , 1961 .

[5]  C.E. Shannon,et al.  Communication in the Presence of Noise , 1949, Proceedings of the IRE.

[6]  S. Dalziel,et al.  Self-similarity and internal structure of turbulence induced by Rayleigh–Taylor instability , 1999, Journal of Fluid Mechanics.

[7]  S. Chandrasekhar.  The character of the equilibrium of an incompressible heavy viscous fluid of variable density , 1955, Mathematical Proceedings of the Cambridge Philosophical Society.

[8]  G. Taylor.  The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I , 1950, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[9]  D. Youngs,et al.  Numerical simulation of mixing by Rayleigh-Taylor and Richtmyer-Meshkov instabilities , 1994 .

[10]  Dale R. Durran,et al.  The Third-Order Adams-Bashforth Method: An Attractive Alternative to Leapfrog Time Differencing , 1991 .

[11]  L. Shibarshov,et al.  Turbulent mixing at an accelerating interface between liquids of different density , 1978 .

[12]  C. W. Hirt,et al.  Effects of Diffusion on Interface Instability between Gases , 1962 .

[13]  Paul Linden,et al.  Molecular mixing in Rayleigh–Taylor instability , 1994, Journal of Fluid Mechanics.

[14]  R. D. Richtmyer.  Taylor instability in shock acceleration of compressible fluids , 1960 .

[15]  P. Dimotakis,et al.  On the geometry of two-dimensional slices of irregular level sets in turbulent flows , 1999 .

[16]  L. Collatz.  The numerical treatment of differential equations , 1961 .

[17]  P. Dimotakis.  The mixing transition in turbulent flows , 2000, Journal of Fluid Mechanics.

[18]  P. Dimotakis.  Turbulent Free Shear Layer Mixing and Combustion , 1991 .

[19]  F. Beux,et al.  The effect of the numerical scheme on the subgrid scale term in large-eddy simulation , 1998 .

[20]  P. Dimotakis.  Some issues on turbulent mixing and turbulence , 1993 .

[21]  P. Miller.  Mixing in High Schmidt Number Turbulent Jets , 1991 .

[22]  C. L. B. O N D A N.  Turbulent shear-layer mixing at high Reynolds numbers : effects of inflow conditions , 2022 .

[23]  J. Haas,et al.  Experimental investigation into inertial properties of Rayleigh–Taylor turbulence , 1997 .

[24]  R. Rosner,et al.  On the miscible Rayleigh–Taylor instability: two and three dimensions , 2001, Journal of Fluid Mechanics.

[25]  D. Youngs,et al.  Modelling turbulent mixing by Rayleigh-Taylor instability , 1989 .

[26]  J. Bowles,et al.  Fourier Analysis of Numerical Approximations of Hyperbolic Equations , 1987 .

[27]  D. Youngs,et al.  Numerical simulation of turbulent mixing by Rayleigh-Taylor instability , 1984 .

[28]  D. Spalding,et al.  A simple experiment to investigate two-dimensional mixing by Rayleigh-Taylor instability , 1990 .

[29]  T. Matsuno.  Numerical Integrations of the Primitive Equations by a Simulated Backward Difference Method , 1966 .

[30]  S.N.B. Murthy,et al.  Turbulent Free Shear Layer Mixing and Combustion , 1991 .

[31]  Marilyn Schneider,et al.  Density ratio dependence of Rayleigh–Taylor mixing for sustained and impulsive acceleration histories , 2000 .

[32]  D. L. Sandoval,et al.  The dynamics of variable-density turbulence , 1995 .

[33]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[34]  M. Glauser,et al.  Theoretical and Computational Fluid Dynamics the Structure of Inhomogeneous Turbulence in Variable Density Nonpremixed Flames 1 , 2022 .

[35]  K. I. Read,et al.  Experimental investigation of turbulent mixing by Rayleigh-Taylor instability , 1984 .