-
Local Wave Number Model for Inhomogeneous Two-Fluid Mixing
Authors:
Nairita Pal,
Ismael Boureima,
Noah Braun,
Susan Kurien,
Praveen Ramaprabhu,
Andrew Lawrie
Abstract:
We present a study of a two-point spectral turbulence model (Local Wave-Number model or LWN model) for the Rayleigh-Taylor (RT) instability. The model outcomes are compared with statistical quantities extracted from three-dimensional simulation of the RT problem. These simulations are initialized with high wavenumber perturbations at the interface of a heavy fluid placed on top of a light fluid so…
▽ More
We present a study of a two-point spectral turbulence model (Local Wave-Number model or LWN model) for the Rayleigh-Taylor (RT) instability. The model outcomes are compared with statistical quantities extracted from three-dimensional simulation of the RT problem. These simulations are initialized with high wavenumber perturbations at the interface of a heavy fluid placed on top of a light fluid so that the density gradient is in the direction opposite to acceleration due to gravity. We consider flows of low to medium density contrast and compare the LWN model against simulation data using the mix-width evolution as the primary metric. The original model specified physically reasonable but largely \emph{ad hoc} terms to account for the inhomogeneous mechanisms involved in growing the mixing layer. We systematically assess the role of each of the terms in the LWN model equations by comparison with simulation. Two of these, the kinematic source term, introduced to maintain a finite covariance between density and specific volume, and a spectral distortion term, introduced as spectral modifications of the density-specific-volume covariance, both result in severely over-predicting the mix layer growth. A simplified model eliminating those two terms is shown to improve the capture of both mix-width evolution as well as the turbulent mass flux velocity profiles across the mix layer at different times. However, this simplification reveals that fidelity to other metrics such as the density-specific-volume covariance, and the turbulent kinetic energy are somewhat compromised. The implications of this outcome are discussed with respect to the physics of the RT problem, and we provide this study as a guide for the practical use of such a model.
△ Less
Submitted 16 August, 2021; v1 submitted 18 March, 2020;
originally announced March 2020.
-
Turbulent transport and mixing in the multimode narrowband Richtmyer-Meshkov instability
Authors:
B. Thornber,
J. Griffond,
P. Bigdelou,
I. Boureima,
P. Ramaprabhu,
O. Schilling,
R. J. R. Williams
Abstract:
The mean momentum and heavy mass fraction, turbulent kinetic energy, and heavy mass fraction variance fields, as well as the budgets of their transport equations, are examined at several times during the evolution of a narrowband Richtmyer-Meshkov instability initiated by a Mach 1.84 shock traversing a perturbed interface separating gases with a density ratio of 3. The results are computed using t…
▽ More
The mean momentum and heavy mass fraction, turbulent kinetic energy, and heavy mass fraction variance fields, as well as the budgets of their transport equations, are examined at several times during the evolution of a narrowband Richtmyer-Meshkov instability initiated by a Mach 1.84 shock traversing a perturbed interface separating gases with a density ratio of 3. The results are computed using the `quarter scale' data from four algorithms presented in the θ-group study of Thornber et al. [Phys. Fluids 29, 105107 (2017)]. The present study is inspired by a previous similar study of Rayleigh-Taylor instability and mixing using direct numerical simulation data by Schilling and Mueschke [Phys. Fluids 22, 105102 (2010)]. In addition to comparing the predictions of the data from four implicit large-eddy simulation codes, the budgets are used to quantify the relative importance of the terms in the transport equations, and the balance of the terms is employed to infer the numerical dissipation. Terms arising from the compressibility of the flow are examined in particular, i.e., the pressure-dilatation. The results are useful for validation of large-eddy simulation and Reynolds-averaged modeling of Richtmyer-Meshkov instability.
△ Less
Submitted 31 October, 2019;
originally announced October 2019.
-
Late-time growth rate, mixing and anisotropy in the multimode narrowband Richtmyer--Meshkov Instability: the $θ$-Group Collaboration
Authors:
B. Thornber,
J. Griffond,
O. Poujade,
N. Attal,
H. Varshochi,
P. Bigdelou,
P. Ramaprabhu,
B. Olson,
J. Greenough,
Y. Zhou,
O. Schilling,
K. A. Garside,
R. J. R. Williams,
C. A. Batha,
P. A. Kuchugov,
M. E. Ladonkina,
V. F. Tishkin,
N. V. Zmitrenko,
V. B. Rozanov,
D. L. Youngs
Abstract:
Turbulent Richtmyer--Meshkov instability (RMI) is investigated through a series of high resolution three dimensional smulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a p…
▽ More
Turbulent Richtmyer--Meshkov instability (RMI) is investigated through a series of high resolution three dimensional smulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm, and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regime combined. Overall, the growth rate exponent is determined as $θ=0.292 \pm 0.009$, in good agreement with prior studies; however, the exponent is decaying slowly in time. $θ$ is shown to be relatively insensitive to the choice of mixing layer width measurement. The asymptotic integral molecular mixing measures $Θ=0.792\pm 0.014$, $Ξ=0.800 \pm 0.014$ and $Ψ=0.782\pm 0.013$ which are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49\% and 66\% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single $θ$. Quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.
△ Less
Submitted 29 June, 2017;
originally announced June 2017.
-
Scrambled and Unscrambled Turbulence
Authors:
P. Ramaprabhu,
V. Karkhanis,
A. G. W. Lawrie
Abstract:
The linked fluid dynamics videos depict Rayleigh-Taylor turbulence when driven by a complex acceleration profile involving two stages of acceleration interspersed with a stage of stabilizing deceleration. Rayleigh-Taylor (RT) instability occurs at the interface separating two fluids of different densities, when the lighter fluid is accelerated in to the heavier fluid. The turbulent mixing arising…
▽ More
The linked fluid dynamics videos depict Rayleigh-Taylor turbulence when driven by a complex acceleration profile involving two stages of acceleration interspersed with a stage of stabilizing deceleration. Rayleigh-Taylor (RT) instability occurs at the interface separating two fluids of different densities, when the lighter fluid is accelerated in to the heavier fluid. The turbulent mixing arising from the development of the miscible RT instability is of key importance in the design of Inertial Confinement Fusion capsules, and to the understanding of astrophysical events, such as Type Ia supernovae. By driving this flow with an accel-decel-accel profile, we have investigated how structures in RT turbulence are affected by a sudden change in the direction of the acceleration first from destabilizing acceleration to deceleration, and followed by a restoration of the unstable acceleration. By studying turbulence under such highly non-equilibrium conditions, we hope to develop an understanding of the response and recovery of self-similar turbulence to sudden changes in the driving acceleration.
△ Less
Submitted 11 October, 2013;
originally announced October 2013.