A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

Thomas J. Rehagen; Jeffrey A. Greenough; Britton J. Olson

doi:10.1115/1.4035944

Journal of Fluids Engineering | Volume 139 | Issue 6 | research-article

< Previous Article Next Article >

Research Papers: Fundamental Issues and Canonical Flows

A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

Thomas J. Rehagen, Jeffrey A. Greenough and Britton J. Olson

[+] Author and Article Information

Thomas J. Rehagen

Lawrence Livermore National Laboratory,
Livermore, CA 94550
e-mail: rehagen1@llnl.gov

Jeffrey A. Greenough

Lawrence Livermore National Laboratory,
Livermore, CA 94550,

Britton J. Olson

Lawrence Livermore National Laboratory,
Livermore, CA 94550

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received October 20, 2016; final manuscript received January 13, 2017; published online April 20, 2017. Assoc. Editor: Praveen Ramaprabhu.The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States Government purposes.

J. Fluids Eng 139(6), 061204 (Apr 20, 2017) (9 pages) Paper No: FE-16-1689; doi: 10.1115/1.4035944 History: Received October 20, 2016; Revised January 13, 2017

Article

References

Figures

Tables

Abstract

The compressible Rayleigh–Taylor (RT) instability is studied by performing a suite of large eddy simulations (LES) using the Miranda and Ares codes. A grid convergence study is carried out for each of these computational methods, and the convergence properties of integral mixing diagnostics and late-time spectra are established. A comparison between the methods is made using the data from the highest resolution simulations in order to validate the Ares hydro scheme. We find that the integral mixing measures, which capture the global properties of the RT instability, show good agreement between the two codes at this resolution. The late-time turbulent kinetic energy and mass fraction spectra roughly follow a Kolmogorov spectrum, and drop off as k approaches the Nyquist wave number of each simulation. The spectra from the highest resolution Miranda simulation follow a Kolmogorov spectrum for longer than the corresponding spectra from the Ares simulation, and have a more abrupt drop off at high wave numbers. The growth rate is determined to be between around 0.03 and 0.05 at late times; however, it has not fully converged by the end of the simulation. Finally, we study the transition from direct numerical simulation (DNS) to LES. The highest resolution simulations become LES at around t/τ ≃ 1.5. To have a fully resolved DNS through the end of our simulations, the grid spacing must be 3.6 (3.1) times finer than our highest resolution mesh when using Miranda (Ares).

FIGURES IN THIS ARTICLE

Purchase this Content

$25.00

Purchase

Learn about subscription and purchase options

Your Session has timed out. Please sign back in to continue.

References

Rayleigh, L. , 1883, “ Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density,” Proc. R. Math. Soc., 14(1), pp. 170–177.

Taylor, G. , 1950, “ The Instability of Liquid Surfaces When Accelerated in a Direction Perpendicular to Their Planes—I,” Proc. R. Soc. London Ser. A, 201(1065), pp. 192–196. [CrossRef]

Chandrasekhar, S. , 1955, “ The Character of the Equilibrium of an Incompressible Heavy Viscous Fluid of Variable Density,” Math. Proc. Cambridge Philos. Soc., 51(1), pp. 162–178. [CrossRef]

Youngs, D. L. , 1984, “ Numerical Simulation of Turbulent Mixing by Rayleigh–Taylor Instability,” Phys. D, 12(1–3), pp. 32–44. [CrossRef]

Youngs, D. L. , 1989, “ Modelling Turbulent Mixing by Rayleigh–Taylor Instability,” Physica D, 37(1), pp. 270–287. [CrossRef]

Sharp, D. H. , 1984, “ An Overview of Rayleigh–Taylor Instability,” Physica D, 12(1–3), p. 3. [CrossRef]

Haan, S. W. , 1989, “ Onset of Nonlinear Saturation for Rayleigh–Taylor Growth in the Presence of a Full Spectrum of Modes,” Phys. Rev. A, 39(11), pp. 5812–5825. [CrossRef]

Haan, S. W. , 1991, “ Weakly Nonlinear Hydrodynamic Instabilities in Inertial Fusion,” Phys. Fluids B, 3(8), pp. 2349–2355. [CrossRef]

Ofer, D. , Alon, U. , Shvarts, D. , McCrory, R. L. , and Verdon, C. P. , 1996, “ Modal Model for the Nonlinear Multimode Rayleigh–Taylor Instability,” Phys. Plasmas, 3(8), pp. 3073–3090. [CrossRef]

Shvarts, D. , Alon, U. , Ofer, D. , McCrory, R. L. , and Verdon, C. P. , 1995, “ Nonlinear Evolution of Multimode Rayleigh–Taylor Instability in Two and Three Dimensions,” Phys. Plasmas, 2(6), pp. 2465–2472. [CrossRef]

Kartoon, D. , Oron, D. , Arazi, L. , and Shvarts, D. , 2003, “ Three-Dimensional Multimode Rayleigh–Taylor and Richtmyer–Meshkov Instabilities at All Density Ratios,” Laser Part. Beams, 21(3), pp. 327–334. [CrossRef]

Abarzhi, S. I. , 2010, “ On Fundamentals of Rayleigh–Taylor Turbulent Mixing,” Europhys. Lett., 91(3), p. 35001. [CrossRef]

Anisimov, S. I. , Drake, R. P. , Gauthier, S. , Meshkov, E. E. , and Abarzhi, S. I. , 2013, “ What is Certain and What is not so Certain in Our Knowledge of Rayleigh–Taylor Mixing?,” Philos. Trans. R. Soc. London Ser. A, 371(2003), p. 20130266. [CrossRef]

Layzer, D. , 1955, “ On the Instability of Superposed Fluids in a Gravitational Field,” Astrophys. J., 122, p. 1. [CrossRef]

Hecht, J. , Alon, U. , and Shvarts, D. , 1994, “ Potential Flow Models of Rayleigh–Taylor and Richtmyer–Meshkov Bubble Fronts,” Phys. Fluids, 6(12), pp. 4019–4030. [CrossRef]

Oron, D. , Arazi, L. , Kartoon, D. , Rikanati, A. , Alon, U. , and Shvarts, D. , 2001, “ Dimensionality Dependence of the Rayleigh–Taylor and Richtmyer–Meshkov Instability Late-Time Scaling Laws,” Phys. Plasmas, 8(6), pp. 2883–2889. [CrossRef]

Goncharov, V. N. , 2002, “ Analytical Model of Nonlinear, Single-Mode, Classical Rayleigh–Taylor Instability at Arbitrary Atwood Numbers,” Phys. Rev. Lett., 88(13), p. 134502. [CrossRef] [PubMed]

Mikaelian, K. O. , 2003, “ Explicit Expressions for the Evolution of Single-Mode Rayleigh–Taylor and Richtmyer–Meshkov Instabilities at Arbitrary Atwood Numbers,” Phys. Rev. E, 67(2), p. 026319. [CrossRef]

Srebro, Y. , Elbaz, Y. , Sadot, O. , Arazi, L. , and Shvarts, D. , 2003, “ A General Buoyancy-Drag Model for the Evolution of the Rayleigh–Taylor and Richtmyer–Meshkov Instabilities,” Laser Part. Beams, 21(3), pp. 347–353.

Jin, H. , Liu, X. F. , Lu, T. , Cheng, B. , Glimm, J. , and Sharp, D. H. , 2005, “ Rayleigh–Taylor Mixing Rates for Compressible Flow,” Phys. Fluids, 17(2), p. 024104. [CrossRef]

George, E. , and Glimm, J. , 2005, “ Self-Similarity of Rayleigh–Taylor Mixing Rates,” Phys. Fluids, 17(5), p. 054101. [CrossRef]

Lafay, M.-A. , LeCreurer, B. , and Gauthier, S. , 2007, “ Compressibility Effects on the Rayleigh–Taylor Instability Between Miscible Fluids,” Europhys. Lett., 79(6), p. 64002. [CrossRef]

Livescu, D. , and Ristorcelli, J. R. , 2008, “ Variable-Density Mixing in Buoyancy-Driven Turbulence,” J. Fluid Mech., 605, pp. 145–180. [CrossRef]

Livescu, D. , Ristorcelli, J. R. , Petersen, M. R. , and Gore, R. A. , 2010, “ New Phenomena in Variable-Density Rayleigh–Taylor Turbulence,” Phys. Scr. Vol. T, 142(1), p. 014015. [CrossRef]

Gauthier, S. , 2013, “ Compressibility Effects in Rayleigh–Taylor Flows: Influence of the Stratification,” Phys. Scr. Vol. T, 155(1), p. 014012. [CrossRef]

Reckinger, S. J. , Livescu, D. , and Vasilyev, O. V. , 2016, “ Comprehensive Numerical Methodology for Direct Numerical Simulations of Compressible Rayleigh–Taylor Instability,” J. Comput. Phys., 313, pp. 181–208. [CrossRef]

Mellado, J. P. , Sarkar, S. , and Zhou, Y. , 2005, “ Large-Eddy Simulation of Rayleigh–Taylor Turbulence With Compressible Miscible Fluids,” Phys. Fluids, 17(7), p. 076101. [CrossRef]

Olson, B. J. , and Cook, A. W. , 2007, “ Rayleigh–Taylor Shock Waves,” Phys. Fluids, 19(12), p. 128108. [CrossRef]

Cook, A. W. , Cabot, W. , and Miller, P. L. , 2004, “ The Mixing Transition in Rayleigh–Taylor Instability,” J. Fluid Mech., 511, pp. 333–362. [CrossRef]

Cabot, W. H. , and Cook, A. W. , 2006, “ Reynolds Number Effects on Rayleigh–Taylor Instability With Possible Implications for Type IA Supernovae,” Nat. Phys., 2(8), pp. 562–568. [CrossRef]

Olson, B. J. , Larsson, J. , Lele, S. K. , and Cook, A. W. , 2011, “ Nonlinear Effects in the Combined Rayleigh–Taylor/Kelvin–Helmholtz Instability,” Phys. Fluids, 23(11), p. 114107. [CrossRef]

Cook, A. W. , 2007, “ Artificial Fluid Properties for Large-Eddy Simulation of Compressible Turbulent Mixing,” Phys. Fluids, 19(5), p. 055103. [CrossRef]

Cook, A. W. , 2009, “ Enthalpy Diffusion in Multicomponent Flows,” Phys. Fluids, 21(5), p. 055109. [CrossRef]

Mani, A. , Larsson, J. , and Moin, P. , 2009, “ Suitability of Artificial Bulk Viscosity for Large-Eddy Simulation of Turbulent Flows With Shocks,” J. Comput. Phys., 228(19), pp. 7368–7374. [CrossRef]

Johnsen, E. , Larsson, J. , Bhagatwala, A. V. , Cabot, W. H. , Moin, P. , Olson, B. J. , Rawat, P. S. , Shankar, S. K. , Sjögreen, B. , Yee, H. C. , Zhong, X. , and Lele, S. K. , 2010, “ Assessment of High-Resolution Methods for Numerical Simulations of Compressible Turbulence With Shock Waves,” J. Comput. Phys., 229(4), pp. 1213–1237. [CrossRef]

Olson, B. J. , and Lele, S. K. , 2013, “ Directional Artificial Fluid Properties for Compressible Large-Eddy Simulation,” J. Comput. Phys., 246, pp. 207–220. [CrossRef]

Darlington, R. M. , McAbee, T. L. , and Rodrigue, G. , 2001, “ A Study of ALE Simulations of Rayleigh–Taylor Instability,” Comput. Phys. Commun., 135(1), pp. 58–73. [CrossRef]

Wilkins, M. , 1963, “ Calculation of Elastic-Plastic Flow,” Lawrence Radiation Laboratory, Report No. UCRL-7322.

Kolev, T. V. , and Rieben, R. N. , 2009, “ A Tensor Artificial Viscosity Using a Finite Element Approach,” J. Comput. Phys., 228(22), pp. 8336–8366. [CrossRef]

Sharp, R. , and Barton, R. , 1981, “ Hemp Advection Model,” Lawrence Livermore Laboratory, Livermore, CA, Report No. UCID 17809.

Berger, M. J. , and Oliger, J. , 1984, “ Adaptive Mesh Refinement for Hyperbolic Partial Differential Equations,” J. Comput. Phys., 53(3), pp. 484–512. [CrossRef]

Berger, M. J. , and Colella, P. , 1989, “ Local Adaptive Mesh Refinement for Shock Hydrodynamics,” J. Comput. Phys., 82(1), pp. 64–84. [CrossRef]

Dimonte, G. , 2004, “ Dependence of Turbulent Rayleigh–Taylor Instability on Initial Perturbations,” Phys. Rev. E, 69(5), p. 056305. [CrossRef]

Cook, A. W. , and Dimotakis, P. E. , 2001., “ Transition Stages of Rayleigh–Taylor Instability Between Miscible Fluids,” J. Fluid Mech., 443(9), pp. 69–99.

Ristorcelli, J. R. , and Clark, T. T. , 2004, “ Rayleigh–Taylor Turbulence: Self-Similar Analysis and Direct Numerical Simulations,” J. Fluid Mech., 507(5), pp. 213–253. [CrossRef]

Dimonte, G. , Youngs, D. L. , Dimits, A. , Weber, S. , Marinak, M. , Wunsch, S. , Garasi, C. , Robinson, A. , Andrews, M. J. , Ramaprabhu, P. , Calder, A. C. , Fryxell, B. , Biello, J. , Dursi, L. , MacNeice, P. , Olson, K. , Ricker, P. , Rosner, R. , Timmes, F. , Tufo, H. , Young, Y.-N. , and Zingale, M. , 2004, “ A Comparative Study of the Turbulent Rayleigh–Taylor Instability Using High-Resolution Three-Dimensional Numerical Simulations: The Alpha-Group Collaboration,” Phys. Fluids, 16(5), pp. 1668–1693. [CrossRef]

Ramaprabhu, P. , Dimonte, G. , and Andrews, M. J. , 2005, “ A Numerical Study of the Influence of Initial Perturbations on the Turbulent Rayleigh–Taylor Instability,” J. Fluid Mech., 536(8), pp. 285–319. [CrossRef]

Olson, D. H. , and Jacobs, J. W. , 2009, “ Experimental Study of Rayleigh–Taylor Instability With a Complex Initial Perturbation,” Phys. Fluids, 21(3), p. 034103. [CrossRef]

Thornber, B. , 2016, “ Impact of Domain Size and Statistical Errors in Simulations of Homogeneous Decaying Turbulence and the Richtmyer–Meshkov Instability,” Phys. Fluids, 28(4), p. 045106. [CrossRef]

Youngs, D. L. , 1991, “ Three-Dimensional Numerical Simulation of Turbulent Mixing by Rayleigh–Taylor Instability,” Phys. Fluids A, 3(5), pp. 1312–1320. [CrossRef]

Tritschler, V. K. , Olson, B. J. , Lele, S. K. , Hickel, S. , Hu, X. Y. , and Adams, N. A. , 2014, “ On the Richtmyer–Meshkov Instability Evolving From a Deterministic Multimode Planar Interface,” J. Fluid Mech., 755, pp. 429–462. [CrossRef]

Olson, B. J. , and Greenough, J. A. , 2014, “ Comparison of Two- and Three-Dimensional Simulations of Miscible Richtmyer–Meshkov Instability With Multimode Initial Conditions,” Phys. Fluids, 26(10), p. 101702. [CrossRef]

Morgan, B. E. , and Greenough, J. A. , 2015, “ Large-Eddy and Unsteady RANS Simulations of a Shock-Accelerated Heavy Gas Cylinder,” Shock Waves, 26(4), pp. 355–383. [CrossRef]

Dittrich, T. R. , Hurricane, O. A. , Callahan, D. A. , Dewald, E. L. , Döppner, T. , Hinkel, D. E. , Berzak Hopkins, L. F. , Le Pape, S. , Ma, T. , Milovich, J. L. , Moreno, J. C. , Patel, P. K. , Park, H.-S. , Remington, B. A. , Salmonson, J. D. , and Kline, J. L. , 2014, “ Design of a High-Foot High-Adiabat ICF Capsule for the National Ignition Facility,” Phys. Rev. Lett., 112(5), p. 055002. [CrossRef] [PubMed]

Ma, T. , Patel, P. K. , Izumi, N. , Springer, P. T. , Key, M. H. , Atherton, L. J. , Benedetti, L. R. , Bradley, D. K. , Callahan, D. A. , Celliers, P. M. , Cerjan, C. J. , Clark, D. S. , Dewald, E. L. , Dixit, S. N. , Döppner, T. , Edgell, D. H. , Epstein, R. , Glenn, S. , Grim, G. , Haan, S. W. , Hammel, B. A. , Hicks, D. , Hsing, W. W. , Jones, O. S. , Khan, S. F. , Kilkenny, J. D. , Kline, J. L. , Kyrala, G. A. , Landen, O. L. , Le Pape, S. , MacGowan, B. J. , Mackinnon, A. J. , MacPhee, A. G. , Meezan, N. B. , Moody, J. D. , Pak, A. , Parham, T. , Park, H.-S. , Ralph, J. E. , Regan, S. P. , Remington, B. A. , Robey, H. F. , Ross, J. S. , Spears, B. K. , Smalyuk, V. , Suter, L. J. , Tommasini, R. , Town, R. P. , Weber, S. V. , Lindl, J. D. , Edwards, M. J. , Glenzer, S. H. , and Moses, E. I. , 2013, “ Onset of Hydrodynamic Mix in High-Velocity, Highly Compressed Inertial Confinement Fusion Implosions,” Phys. Rev. Lett., 111(8), p. 085004. [CrossRef] [PubMed]

Olson, B. J. , and Greenough, J. , 2014, “ Large Eddy Simulation Requirements for the Richtmyer–Meshkov Instability,” Phys. Fluids, 26(4), p. 044103. [CrossRef]

Figures

Fig. 1

The initial perturbation spectrum given in Eq. (13)

View Large | Save Figure | Download Slide (.ppt)

Fig. 2

The mixing width as a function of time computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c): (a) Miranda, (b) Ares, and (c) mesh D for Miranda and Ares

View Large | Save Figure | Download Slide (.ppt)

Fig. 3

The mixedness as a function of time computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c): (a) Miranda, (b) Ares, and (c) mesh D for Miranda and Ares

View Large | Save Figure | Download Slide (.ppt)

Fig. 4

The mixed mass as a function of time computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c): (a) Miranda, (b) Ares, and (c) mesh D for Miranda and Ares

View Large | Save Figure | Download Slide (.ppt)

Fig. 5

Spectra of the turbulence kinetic energy at t/τ = 13 computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c). The thin black line above the spectra is thek^–5∕3 fiducial. (a) Miranda (b) Ares, and (c) mesh D for Miranda and Ares

View Large | Save Figure | Download Slide (.ppt)

Fig. 6

Spectra of the mass fraction of the heavier fluid at t/τ = 13 computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c). The thin black line above the spectra is the k^−5∕3 fiducial.

$Spectra of the mass fraction of the heavier fluid at t/τ = 13 computed from the Miranda simulations (a), the Ares simulations (b), and a comparison of the results from the highest resolution mesh, mesh D (c). The thin black line above the spectra is the k−5∕3 fiducial.$

View Large | Save Figure | Download Slide (.ppt)

Fig. 7

The growth rate, α_b as a function of time, computed from the highest resolution Miranda (solid lines) and Ares (dotted lines) simulations, using Eq. (30) (thick lines) and Eq. (31) (thin lines)

View Large | Save Figure | Download Slide (.ppt)

Fig. 8

The effective numerical viscosity as a function of the grid Reynolds number expression, similar to Fig. 15(b) from Ref. [56]. Blue symbols are data from Miranda and red are from Ares. The circle, square, triangle, and diamond values are calculated at t/τ = 0.5, 1, 1.5, and 3 respectively. See text for details.

View Large | Save Figure | Download Slide (.ppt)

Fig. 9

The time dependence of the slope of the data points in Fig. 8. The black points are calculated from the highest resolution (mesh D) data, while the red, blue, and purple points are from mesh C, mesh B, and mesh A, respectively. The dashed gray line is at a slope of −1.4.

View Large | Save Figure | Download Slide (.ppt)

Fig. 10

The equilibrium density profile given in Eq. (A5) as a function of z (solid line), with the values of g, M_H, M_L, T₀, and δ given in Sec. 2. The dashed line shows the density profile if the two fluids were not mixed across the interface initially, i.e., if there is a sharp transition between the heavy and light fluid.

View Large | Save Figure | Download Slide (.ppt)

Tables

Errata

Discussions

Web of Science® Times Cited: 0

Some tools below are only available to our subscribers or users with an online account.

PDF
Email

You must be logged in as an individual user to share content.

Return to: A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

Copyright in the material you requested is held by the American Society of Mechanical Engineers (unless otherwise noted). This email ability is provided as a courtesy, and by using it you agree that you are requesting the material solely for personal, non-commercial use, and that it is subject to the American Society of Mechanical Engineers' Terms of Use. The information provided in order to email this topic will not be used to send unsolicited email, nor will it be furnished to third parties. Please refer to the American Society of Mechanical Engineers' Privacy Policy for further information.
Share
Get Citation

Rehagen TJ, Greenough JA, Olson BJ. A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes. ASME. J. Fluids Eng. 2017;139(6):061204-061204-9. doi:10.1115/1.4035944.

Download citation file:

RIS (Zotero)

EndNote

BibTex

Medlars

ProCite

RefWorks

Reference Manager

© Copyright
Get Alerts

Processing your request... Please Wait.

A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

Abstract

FIGURES IN THIS ARTICLE

References

Figures

Tables

Errata

Discussions

Related Content

Journals

Conference Proceedings

eBooks