Complex Fluid Mixing Flows: Simulation vs. Theory vs. Experiment

J. Glimm,D. Sharp

Abstract:It is an old saw that science has as its three pillars theory, experiment, and simulation. That we can believe simulation when the three do not agree is a testament to a newfound credence for simulations. Our story concerns widely different simulation and experimental results in a problem of long standing: determination of the rates of chaotic mixing resulting from acceleration across a fluid density discontinuity, namely the classic Rayleigh–Taylor (RT) instability. RT and other interface instabilities [1,9] have long attracted the attention of leading applied mathematicians and physicists, not only Rayleigh and Taylor but also Kelvin, Chandrasekhar, Fermi, Garrett Birkhoff, and John Wheeler, among others. These instabilities occur in flows important both to basic science (supernova explosions) and to technology (inertial confinement fusion). The simplest and most important variable used to describe the mixing rate of fluids of different densities is the dimensionless parameter α, which characterizes the penetration distance h of the lighter fluid (in structures called bubbles) into the heavier fluid: h = α Agt. Here, g is the acceleration force, such as gravity; A is the Atwood number, a dimensionless measure of the contrast in the densities (ρ1 and ρ2) of the two fluids: A = (ρ2 – ρ1)/(ρ2 + ρ1); and t is time. In gravitational free fall, α = 1/2; this value of α describes the penetration of the heavier fluid (in structures called spikes) in the case of an extreme density contrast, A = 1. The formula predicts constant acceleration for the mixing front, as explained below (Theory).