skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A numerical study of bubble and spike velocities in shock-driven liquid metals

Abstract

In this study, we use detailed continuum hydrodynamics and molecular dynamics simulations to investigate the dynamics of ejecta that are initialized with large amplitude perturbations and non-sinusoidal shapes. Insights from the simulations are used to suggest a modified expression for the velocity associated with ejected spike structures, whereas a recently suggested model explains the observed bubble velocities. Specifically, we find the asymptotic bubble velocity prediction given by Mikaelian is in excellent agreement with the simulations, when a nonlinear correction for finite amplitudes is used in that model. In contrast, existing models can overpredict observed spike velocities if they do not include the modification of the initial spike growth rates due to nonlinearities. Instead, we find that when potential flow models are corrected with a suitable nonlinear prefactor, this leads to predictions in close agreement with our simulation data. We also propose a simple empirical expression for the nonlinear correction for spike velocities which is able to reproduce results from our simulations and published experimental and simulation data over a wide range of initial conditions and Mach numbers. We discuss extensions of these models to initial interfaces with arbitrary shapes. In particular, for non-sinusoidal shapes, the bubble and spike velocities aremore » still predicted by these models provided we use an effective wavelength λeff which is the wavelength of an equivalent sinusoid that has the same missing area. Lastly, the issues of nonlinearity, non-standard shapes and shock Mach number addressed in this work are relevant to recent experimental campaigns involving twice-shocked targets.« less

Authors:
ORCiD logo [1];  [1]; ORCiD logo [2]; ORCiD logo [2];  [2]
  1. University of North Carolina, Charlotte, NC (United States)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Advanced Research Projects Agency - Energy (ARPA-E)
OSTI Identifier:
1479928
Alternate Identifier(s):
OSTI ID: 1416841
Report Number(s):
LA-UR-17-28876
Journal ID: ISSN 0021-8979
Grant/Contract Number:  
AC52-06NA25396; AC52-06NA2-5396
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Applied Physics
Additional Journal Information:
Journal Volume: 123; Journal Issue: 2; Journal ID: ISSN 0021-8979
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ejecta; Richtmeyer-Meshkov Instability; ejecta models; numerical simulation

Citation Formats

Karkhanis, V., Ramaprabhu, P., Cherne, Frank Joseph, Hammerberg, James Edward, and Andrews, Malcolm John. A numerical study of bubble and spike velocities in shock-driven liquid metals. United States: N. p., 2018. Web. doi:10.1063/1.5008495.
Karkhanis, V., Ramaprabhu, P., Cherne, Frank Joseph, Hammerberg, James Edward, & Andrews, Malcolm John. A numerical study of bubble and spike velocities in shock-driven liquid metals. United States. https://doi.org/10.1063/1.5008495
Karkhanis, V., Ramaprabhu, P., Cherne, Frank Joseph, Hammerberg, James Edward, and Andrews, Malcolm John. 2018. "A numerical study of bubble and spike velocities in shock-driven liquid metals". United States. https://doi.org/10.1063/1.5008495. https://www.osti.gov/servlets/purl/1479928.
@article{osti_1479928,
title = {A numerical study of bubble and spike velocities in shock-driven liquid metals},
author = {Karkhanis, V. and Ramaprabhu, P. and Cherne, Frank Joseph and Hammerberg, James Edward and Andrews, Malcolm John},
abstractNote = {In this study, we use detailed continuum hydrodynamics and molecular dynamics simulations to investigate the dynamics of ejecta that are initialized with large amplitude perturbations and non-sinusoidal shapes. Insights from the simulations are used to suggest a modified expression for the velocity associated with ejected spike structures, whereas a recently suggested model explains the observed bubble velocities. Specifically, we find the asymptotic bubble velocity prediction given by Mikaelian is in excellent agreement with the simulations, when a nonlinear correction for finite amplitudes is used in that model. In contrast, existing models can overpredict observed spike velocities if they do not include the modification of the initial spike growth rates due to nonlinearities. Instead, we find that when potential flow models are corrected with a suitable nonlinear prefactor, this leads to predictions in close agreement with our simulation data. We also propose a simple empirical expression for the nonlinear correction for spike velocities which is able to reproduce results from our simulations and published experimental and simulation data over a wide range of initial conditions and Mach numbers. We discuss extensions of these models to initial interfaces with arbitrary shapes. In particular, for non-sinusoidal shapes, the bubble and spike velocities are still predicted by these models provided we use an effective wavelength λeff which is the wavelength of an equivalent sinusoid that has the same missing area. Lastly, the issues of nonlinearity, non-standard shapes and shock Mach number addressed in this work are relevant to recent experimental campaigns involving twice-shocked targets.},
doi = {10.1063/1.5008495},
url = {https://www.osti.gov/biblio/1479928}, journal = {Journal of Applied Physics},
issn = {0021-8979},
number = 2,
volume = 123,
place = {United States},
year = {Fri Jan 12 00:00:00 EST 2018},
month = {Fri Jan 12 00:00:00 EST 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 18 works
Citation information provided by
Web of Science

Figures / Tables:

Table 1 Table 1: Interface shapes investigated using FLASH (λ= 1 CM) simulations (kh0= 1/8, 0.5, 1 and 2). For Gaussian profiles, we have used the half-height of the perturbation h$$'\atop{0}$$ as discussed in the text.

Save / Share:

Works referenced in this record:

FLASH: An Adaptive Mesh Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes
journal, November 2000


The late time structure of high density contrast, single mode Richtmyer-Meshkov flow
journal, July 2016


On shock driven jetting of liquid from non-sinusoidal surfaces into a vacuum
journal, November 2015


Dynamic Behavior of Materials
book, September 1994


Analytic model of the Grüneisen parameter all densities
journal, August 2004


Unstable Richtmyer–Meshkov growth of solid and liquid metals in vacuum
journal, June 2012


Second shock ejecta measurements with an explosively driven two-shockwave drive
journal, September 2014


Nonlinear Perturbation Theory of the Incompressible Richtmyer-Meshkov Instability
journal, April 1996


Inside the supernova: A powerful convective engine
journal, November 1994


A Source Model for Ejecta
journal, June 2017


A Model of a Source of Shock Wave Metal Ejection Based on Richtmyer–Meshkov Instability Theory
journal, June 2017


Limitations and failures of the Layzer model for hydrodynamic instabilities
journal, July 2008


Analytic Approach to Nonlinear Rayleigh-Taylor and Richtmyer-Meshkov Instabilities
journal, January 1998


Simulations and model of the nonlinear Richtmyer–Meshkov instability
journal, January 2010


Ejecta source model based on the nonlinear Richtmyer-Meshkov instability
journal, January 2013


Ejecta production from shocked Pb surface via molecular dynamics
journal, October 2014


Shock induced polymorphic transition and melting of tin
conference, January 2000


Richtmyer–Meshkov instability growth: experiment, simulation and theory
journal, June 1999


Controlling the Data Glut in Large-Scale Molecular-Dynamics Simulations
journal, January 1997


Semiempirical, Quantum Mechanical Calculation of Hydrogen Embrittlement in Metals
journal, April 1983


Three‐dimensional Numerical Simulations of Rayleigh‐Taylor Unstable Flames in Type Ia Supernovae
journal, October 2005


Atomistic simulations of shock-induced microjet from a grooved aluminium surface
journal, April 2013


Explosively driven two-shockwave tools with applications
journal, May 2014


Ejecta Production from Second Shock: Numerical Simulations and Experiments
journal, February 2017


50 GFlops molecular dynamics on the Connection Machine 5
conference, January 1993


Use of the Richtmyer-Meshkov Instability to Infer Yield Stress at High-Energy Densities
journal, December 2011


The Piecewise Parabolic Method (PPM) for gas-dynamical simulations
journal, April 1984


Instability of the interface of two gases accelerated by a shock wave
journal, January 1972


MOLECULAR DYNAMICS COMES OF AGE: 320 BILLION ATOM SIMULATION ON BlueGene/L
journal, December 2006


Freeze‐out and the effect of compressibility in the Richtmyer–Meshkov instability
journal, January 1994


Estimation of Metal Strength at Very High Rates Using Free-Surface Richtmyer–Meshkov Instabilities
journal, March 2017


Taylor instability in shock acceleration of compressible fluids
journal, May 1960


Foreword to the Special Issue on Ejecta
journal, May 2017


Parallel replica method for dynamics of infrequent events
journal, June 1998


On the Instability of Superposed Fluids in a Gravitational Field.
journal, July 1955


Instabilities and Clumping in Type Ia Supernova Remnants
journal, March 2001


Works referencing / citing this record:

Peculiarities in breakup and transport process of shock-induced ejecta with surrounding gas
journal, May 2019


Ejecta velocities in twice-shocked liquid metals under extreme conditions: A hydrodynamic approach
journal, July 2019