Hostname: page-component-77c89778f8-m8s7h Total loading time: 0 Render date: 2024-07-18T21:12:47.271Z Has data issue: false hasContentIssue false

Variable-density mixing in buoyancy-driven turbulence

Published online by Cambridge University Press:  23 May 2008

D. LIVESCU
Affiliation:
Los Alamos National Laboratory, Los Alamos, NM 87544, USAlivescu@lanl.gov
J. R. RISTORCELLI
Affiliation:
Los Alamos National Laboratory, Los Alamos, NM 87544, USAlivescu@lanl.gov

Abstract

The homogenization of a heterogeneous mixture of two pure fluids with different densities by molecular diffusion and stirring induced by buoyancy-generated motions, as occurs in the Rayleigh–Taylor (RT) instability, is studied using direct numerical simulations. The Schmidt number, Sc, is varied by a factor of 20, 0.1 ≤ Sc ≤ 2.0, and the Atwood number, A, by a factor of 10, 0.05 ≤ A ≤ 0.5. Initial-density intensities are as high as 50% of the mean density. As a consequence of differential accelerations experienced by the two fluids, substantial and important differences between the mixing in a variable-density flow, as compared to the Boussinesq approximation, are observed. In short, the pure heavy fluid mixes more slowly than the pure light fluid: an initially symmetric double delta density probability density function (PDF) is rapidly skewed and, only at long times and low density fluctuations, does it relax to a Gaussian-like PDF. The heavy–light fluid mixing process asymmetry is relevant to the nature of molecular mixing on different sides of a high-Atwood-number RT layer. Diverse mix metrics are used to examine the homogenization of the two fluids. The conventional mix parameter, θ, is mathematically related to the variance of the excess reactant of a hypothetical fast chemical reaction. Bounds relating θ and the normalized product, Ξ, are derived. It is shown that θ underpredicts the mixing, as compared to Ξ, in the central regions of an RT layer; in the edge regions, θ is larger than Ξ. The shape of the density PDF cannot be inferred from the usual mix metrics popular in applications. For example, when θ, Ξ ≥ 0.6, characteristic of the interior of a fully developed RT layer, the PDFs can have vastly different shapes. Bounds on the fluid composition using two low-order moments of the density PDF are derived. The bounds can be used as realizability conditions for low-dimensional models. For the measures studied, the tightest bounds are obtained using Ξ and mean density. The structure of the flow is also examined. It is found that, at early times, the buoyancy production term in the spectral kinetic energy equation is important at all wavenumbers and leads to anisotropy at all scales of motion. At later times, the anisotropy is confined to the largest and smallest scales: the intermediate scales are more isotropic than the small scales. In the viscous range, there is a cancellation between the viscous and nonlinear effects, and the buoyancy production leads to a persistent small-scale anisotropy.

Type
Papers
Copyright
Copyright © Cambridge University Press 2008

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

Batchelor, G. K., Canuto, V. M. & Chasnov, J. R. 1992 Homogeneous buoyancy-generated turbulence. J. Fluid Mech. 235, 349378.Google Scholar
Besnard, D., Harlow, F. H., Rauenzahn, R. M. & Zemach, C. 1992 Turbulence transport equations for variable-density turbulence and their relationships to two field models. Tech. Rep. Los Alamos National Laboratory, lA-12303-MS.CrossRefGoogle Scholar
Cabot, W. H. & Cook, A. W. 2006 Reynolds number effects on Rayleigh–Taylor instability with possible implications for type-ia supernovae. Nature Phys. 2, 562568.Google Scholar
Cook, A. W. & Dimotakis, P. E. 2001 Transition stages of Rayleigh–Taylor instability between miscible fluids. J. Fluid Mech. 443, 6999.Google Scholar
Cook, A. W. & Zhou, Y. 2002 Energy transfer in Rayleigh–Taylor instability. Phys. Rev. E 66, 026312–1026312–12.Google Scholar
Dalziel, S. B., Linden, P. F. & Youngs, D. L. 1999 Self-similarity and internal structure of turbulence induced by Rayleigh–Taylor instability. J. Fluid Mech. 399, 148.CrossRefGoogle Scholar
Dimonte, G. & Schneider, M. 2000 Density ratio dependence of Rayleigh–Taylor mixing for sustained and impulsive acceleration histories. Phys. Fluids 12, 304321.CrossRefGoogle Scholar
George, E. & Glimm, J. 2005 Self-similarity of Rayleigh–Taylor mixing rates. Phys. Fluids 17, 054101–1–13.Google Scholar
Givi, P. 1989 Model-free simulations of turbulent reacting flows. Prog. Energy Combust. Sci. 15, 1107.Google Scholar
Kraichnan, R. H. 1971 Inertial-range transfer in two- and three-dimensional turbulence. J. Fluid Mech. 47, 525535.Google Scholar
Linden, P. F., Redondo, J. M. & Youngs, D. L. 1994 Molecular mixing in Rayleigh–Taylor instability. J. Fluid Mech. 265, 97124.Google Scholar
Liu, X., George, E., Bo, W. & Glimm, J. 2006 Turbulent mixing with physical mass diffusion. Phys. Rev. E 73, 016304–1–5.Google ScholarPubMed
Livescu, D. & Ristorcelli, J. R. 2007 Buoyancy-driven variable-density turbulence. J. Fluid Mech. 591, 4371.Google Scholar
Livescu, D., Jaberi, F. A. & Madnia, C. K. 2000 Passive scalar wake behind a line source in grid turbulence. J. Fluid Mech. 416, 117149.Google Scholar
Livescu, D., Jaberi, F. A. & Madnia, C. K. 2002 The effects of heat release on the energy exchange in reacting turbulent shear flow. J. Fluid Mech. 450, 3566.Google Scholar
Livescu, D., Ristorcelli, J. R., Gore, R. A., Dean, S. H., Cabot, W. H. & Cook, A. W. 2008 High Reynolds number Rayleigh–Taylor turbulence. J. Turbul. (submitted).Google Scholar
Muenshke, N. J., Andrews, M. J. & Schilling, O. 2006 Experimental characterization of spatio-temporal evolution of small Atwood number Rayleigh–Taylor mixing layer. J. Fluid Mech. 567, 2761.Google Scholar
Pope, S. B. 1985 PDF methods for turbulent reacting flows. Prog. Energy Combust. Sci. 11, 119192.Google Scholar
Pope, S. B. 1988 Evolution of surfaces in turbulence. Intl. J. Engng. Sci. 26, 445469.Google Scholar
Ramaprabhu, P. & Andrews, M. J. 2004 Experimental investigation of Rayleigh–Taylor mixing at small Atwood numbers. J. Fluid Mech. 502, 233271.CrossRefGoogle Scholar
Ristorcelli, J. R. 2006 Passive scalar mixing: analytical study of time scale ratio, variance and mix rate. Phys. Fluids 18, 075101–1075101–17.Google Scholar
Ristorcelli, J. R. & Clark, T. T. 2004 Rayleigh–Taylor turbulence: self-similar analysis and direct numerical simulations. J. Fluid Mech. 507, 213253.CrossRefGoogle Scholar
Sandoval, D. L. 1995 The dynamics of variable density turbulence. PhD thesis, University of Washington, LANL Rep. LA-13037-T.CrossRefGoogle Scholar
Sandoval, D. L., Clark, T. T. & Riley, J. J. 1996 Buoyancy-generated variable density turbulence. In Proc. IUTAM Symp. on Variable Density Low Speed Turbulent Flows, pp. 847–864. Kluwer.Google Scholar
Vervisch, L., Bidaux, E., Bray, K. N. C. & Kollman, W. 1995 Surface density function in premixed turbulent combustion modeling, similarities between probability density function and flame surface approaches. Phys. Fluids 7, 24962503.Google Scholar
Warhaft, Z. & Lumley, J. L. 1978 An experimental study of the decay of temperature fluctuations in grid-generated turbulence. J. Fluid Mech. 88, 659684.Google Scholar
Williams, F. A. 1985 Combustion Theory. Perseus Cambridge.Google Scholar
Wilson, P. N. & Andrews, M. J. 2001 Spectral measurements of Rayleigh–Taylor mixing at low Atwood number. Phys. Fluids 14, 938945.Google Scholar
Young, Y. N., Tufo, H., Dubey, A. & Rosner, R. 1999 On the miscible Rayleigh–Taylor instability: two and three dimensions. J. Fluid Mech. 447, 377408.Google Scholar
Youngs, D. L. 1991 Three-dimensional numerical simulation of turbulent mixing by Rayleigh–Taylor instability. Phys. Fluids A 3, 13121320.Google Scholar
Youngs, D. L. 1994 Numerical simulation of mixing in Rayleigh–Taylor and Richtmeyer–Meshkov instabilities. Laser Part. Beams 12, 725750.Google Scholar