Comparison of two- and three-dimensional simulations of miscible Richtmyer-Meshkov instability with multimode initial conditions

Olson, Britton J.; Greenough, Jeffrey A.

doi:10.1063/1.4898157

Cited by 38 publications

(12 citation statements)

References 17 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…At this Reynolds number, one expects three-dimensional effects to become important, and the flow to become turbulent. Nevertheless, consistent with previous studies (see for example [28,29]), we neglect the three-dimensional effects, leaving that for further study in the future. The shock-tube numerical experiments were performed in a two-dimensional domain of size [15λ 0 × 2.5λ 0 ] discretized with a base mesh of resolution [480 × 80].…”

Section: Physical Domain and Simulation Setupsupporting

confidence: 57%

Effects of solid circular cylinders on the dynamics of the Richtmyer–Meshkov instability

Al-Marouf

Samtaney

2020

Shock Waves

View full text Add to dashboard Cite

We discuss the effect on the transient dynamics of the Richtmyer-Meshkov instability due to small-scale perturbations introduced by a set of solid circular cylinders located along the material interface. The arrangement of the cylinders mimics (in a twodimensional domain) the presence of the solid supporting grid wires used in the formation of the material interface in an experimental setup. The numerical experiments are performed by utilizing an embedded boundary-Adaptive mesh refinement method. We quantify the effect of introducing the cylinders on the mixing layer growth and the mixedness level, and qualitatively demonstrate the effect of the perturbation introduced by the cylinders on the mixing layer structure. We empirically model the effect of the solid cylinders based on their diameter, and find that the effect of the cylinders on the mixing layer growth is a linear function of the cylinder diameter. The linear model allows us to predict the change in the mixing layer growth within a ±10% range, except for very small ratios of the initial perturbation amplitude of the density interface to the cylinder diameter.

show abstract

Section: Physical Domain and Simulation Setupsupporting

confidence: 57%

Effects of solid circular cylinders on the dynamics of the Richtmyer–Meshkov instability

Al-Marouf

Samtaney

2020

Shock Waves

View full text Add to dashboard Cite

show abstract

“…The decrease of the mixing width around 3.3 ms corresponds to compression from the reflected shock, while the inflection around 4.5 ms corresponds to interaction of the mixing layer with the rarefaction wave. Such features are observed in both LES [46,62,63] as well as RANS simulations [29,57,64,65] of reshocked RM mixing layers. Morgan and Wickett previously demonstrated a similar level of agreement between simulation and experiment for the k-L-a model [56].…”

Section: B Richtmyer-meshkov Mixingmentioning

confidence: 79%

Two-length-scale turbulence model for self-similar buoyancy-, shock-, and shear-driven mixing

2018

View full text Add to dashboard Cite

The three-equation k-L-a turbulence model [B. Morgan and M. Wickett, Three-equation model for the self-similar growth of Rayleigh-Taylor and Richtmyer-Meshkov instabilities, Phys. Rev. E 91, 043002 (2015)PLEEE81539-375510.1103/PhysRevE.91.043002] is extended by the addition of a second length scale equation. It is shown that the separation of turbulence transport and turbulence destruction length scales is necessary for simultaneous prediction of the growth parameter and turbulence intensity of a Kelvin-Helmholtz shear layer when model coefficients are constrained by similarity analysis. Constraints on model coefficients are derived that satisfy an ansatz of self-similarity in the low-Atwood-number limit and allow the determination of model coefficients necessary to recover expected experimental behavior. The model is then applied in one-dimensional simulations of Rayleigh-Taylor, reshocked Richtmyer-Meshkov, Kelvin-Helmholtz, and combined Rayleigh-Taylor-Kelvin-Helmholtz instability mixing layers to demonstrate that the expected growth rates are recovered numerically. Finally, it is shown in the case of combined instability that the model predicts a mixing width that is a linear combination of Rayleigh-Taylor and Kelvin-Helmholtz mixing processes.

show abstract

“…The Miranda code solves the hydrodynamics equations presented in section II A with a tenth-order compact differencing scheme for spatial discretization and a fourthorder explicit Runge-Kutta scheme for temporal integration. It has been utilized extensively in previous studies of RT and RM mixing [40][41][42][43][44][45][46][47]. To model the subgrid scale (SGS) transfer of energy, Miranda utilizes an artificial fluid LES (AFLES) approach in which artificial transport terms are added to the fluid viscosity, the bulk viscosity, the thermal conductivity, and the molecular diffusivity [48,49].…”

Section: B the Miranda Codementioning

confidence: 99%

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of a reacting Rayleigh-Taylor mixing layer in a spherical geometry

et al. 2018

Self Cite

View full text Add to dashboard Cite

Tenth-order compact difference code Miranda is used to perform large-eddy simulation (LES) of a hydrogen gas/plastic mixing layer in a spherical geometry. Once the mixing layer has achieved selfsimilar growth, it is heated to 1 keV, and the second-order arbitrary Lagrangian/Eulerian (ALE) code Ares is used to simulate mixing layer evolution as it undergoes thermonuclear (TN) burn. Both premixed (in which deuterium and tritium are initially present in the gas) and non-premixed (in which deuterium is initially present only in the plastic) variants are considered at Atwood numbers 0.05 and 0.50. The impact of turbulent mixing on mean TN reaction rate is examined, and a four-equation k-L-a-V Reynolds-averaged Navier Stokes (RANS) model is presented. The k-L-a-V model, which represents an extension of the k-L-a model [B. Morgan and M. Wickett, "Three-equation model for the self-similar growth of Rayleigh-Taylor and Richtmyer-Meshkov instabilities," Phys. Rev. E 91, 043002 (2015)] by the addition of a transport equation for the scalar mass fraction variance, is then applied in one-dimensional simulations of the reacting mixing layer under consideration. Excellent agreement is obtained between LES and RANS in total TN neutron production when fluctuations in reaction cross-section can be neglected.

show abstract

Comparison of two- and three-dimensional simulations of miscible Richtmyer-Meshkov instability with multimode initial conditions

Cited by 38 publications

References 17 publications

Effects of solid circular cylinders on the dynamics of the Richtmyer–Meshkov instability

Effects of solid circular cylinders on the dynamics of the Richtmyer–Meshkov instability

Two-length-scale turbulence model for self-similar buoyancy-, shock-, and shear-driven mixing

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of a reacting Rayleigh-Taylor mixing layer in a spherical geometry

Contact Info

Product

Resources

About