We study the Rayleigh–Taylor (RT) mixing layer, presenting simulations in agreement with experimental data. This problem is an idealized subproblem of important scientific and engineering problems, such as gravitationally induced mixing in oceanography and performance assessment for inertial confinement fusion. Engineering codes commonly achieve correct simulations through the calibration of adjustable parameters. In this sense, they are interpolative and not predictive. As computational science moves from the interpolative to the predictive and reduces the reliance on experiment, the quality of decision making improves. The diagnosis of errors in a multi-parameter, multi-physics setting is daunting, so we address this issue in the proposed idealized setting. The validation tests presented are thus a test for engineering codes, when used for complex problems containing RT features. The RT growth rate, characterized by a dimensionless but non-universal parameter α , describes the outer edge of the mixing zone. Increasingly accurate front tracking/large eddy simulations reveal the non-universality of the growth rate and agreement with experimental data. Increased mesh resolution allows reduction in the role of key subgrid models. We study the effect of long-wavelength perturbations on the mixing growth rate. A self-similar power law for the initial perturbation amplitudes is here inferred from experimental data. We show a maximum ±5% effect on the growth rate. Large (factors of 2) effects, as predicted in some models and many simulations, are inconsistent with the experimental data of Youngs and co-authors. The inconsistency of the model lies in the treatment of the dynamics of bubbles, which are the shortest-wavelength modes for this problem. An alternative theory for this shortest wavelength, based on the bubble merger model, was previously shown to be consistent with experimental data.
Rayleigh-Taylor mixing is a classical hydrodynamic instability that occurs when a light fluid pushes against a heavy fluid. The two main sources of nonideal behavior in Rayleigh-Taylor (RT) mixing are regularizations (physical and numerical), which produce deviations from a pure Euler equation, scale invariant formulation, and nonideal (i.e., experimental) initial conditions. The Kolmogorov theory of turbulence predicts stirring at all length scales for the Euler fluid equations without regularization. We interpret mathematical theories of existence and nonuniqueness in this context, and we provide numerical evidence for dependence of the RT mixing rate on nonideal regularizations; in other words, indeterminacy when modeled by Euler equations. Operationally, indeterminacy shows up as nonunique solutions for RT mixing, parametrized by Schmidt and Prandtl numbers, in the large Reynolds number (Euler equation) limit. Verification and validation evidence is presented for the large eddy simulation algorithm used here. Mesh convergence depends on breaking the nonuniqueness with explicit use of the laminar Schmidt and Prandtl numbers and their turbulent counterparts, defined in terms of subgrid scale models. The dependence of the mixing rate on the Schmidt and Prandtl numbers and other physical parameters will be illustrated. We demonstrate numerically the influence of initial conditions on the mixing rate. Both the dominant short wavelength initial conditions and long wavelength perturbations are observed to play a role. By examination of two classes of experiments, we observe the absence of a single universal explanation, with long and short wavelength initial conditions, and the various physical and numerical regularizations contributing in different proportions in these two different contexts.large eddy simulations | subgrid scale models | turbulence R ayleigh-Taylor (RT) instability is a classical hydrodynamic instability (1, 2) that occurs when a light fluid pushes against a heavy fluid. The density contrast is measured dimensionlessly by the Atwood number A ¼ ðρ 2 − ρ 1 Þ∕ðρ 2 þ ρ 1 Þ, with ρ i the density of fluid i. The mixing rate can be characterized by the dimensionless parameter α, defined through the equationwhere h is the penetration distance of the light fluid into the heavy and g is the acceleration. The mixing rapidly becomes turbulent, with Reynolds numbers Re ≈ 50;000 observed in typical experiments. For this reason, simulations based on compressible codes, with hyperbolic time step (CFL) restrictions, are unable to resolve the viscous length scales, and such codes are run in an under-resolved manner. Such simulations are called large eddy simulations (LES). LES require subgrid scale (SGS) models to describe the effect of the small (subgrid) scales on the large (grid-resolved) ones. Models for the Reynolds stress based on eddy viscosity are the most familiar of this type.Much of the work on RT mixing has been influenced by ideas of universality and self similarity. If α is to be a universal physical...
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.