Engineering topics
Publications and source records attributed to Hammerberg, James Edward.
In this work, we measure continuous thermal radiance from evolving clouds of liquid metal fragments ejected into vacuum, nonreactive, and reactive gas. We implement a model for the thermalization of the ejecta and gas and use this to constrain the absolute temperature of the ejecta cloud. This model enables further analyses of ejecta thermal behavior under a variety of conditions.
Here, we present a density functional theory study of cerium deuterides, varying deuterium content and temperature. Results are discussed in relation to ejecta experiments recently described in [J. Dynam. Mat. Behav. 1, 12 (2017)]. Using quantum molecular dynamics, we calculate diffusion coefficients, radial distribution functions, and heat capacities of these materials. Properties are calculated over temperatures ranging from 500-4000K at near ambient pressure and at 2.37 fold compression. We find that deuterium diffusion occurs rapidly accross this temperature and compression region and that heat capacity calculations are in agreement with approximations made in experiments. For cerium metal, a comparison of self-diffusion coefficient and viscosity with experimental data and analytical models shows excellent agreement. Lastly, we discuss the stucture of cerium deuterides with respect to radial distribution functions and provide fitted equations for diffusion coefficients over temperature.
We present the experimental validation of ejecta source, mass and velocity models using proton radiography. Ejecta mass and velocity measurements were performed under 7 atm of He gas on Sn targets with 4 concentric sinusoidal features with radial wavelength of λ=2 mm and amplitude η 0 = 0.159 mm. 42 images at 275 and 400 ns intervals were obtained between 0 and 14 µs from the time of shock breakout. Volume density profiles for the spikes were determined through Abel inversion and spike, bubble and surface velocities were determined using proton radiography and also verified through Photon Doppler Velocimetry (PDV). Ejecta sheets evolve to a self-similar density distribution that depends on a scaling variable z/η s t where η˙s is the spike tip velocity, z is the distance from the free surface and t is the time after shock breakout. Both the density profiles and the time dependence of the mass per unit area in the evolving spikes are in good agreement with a Richtmyer-Meshkov instability based model for ejecta production and evolution.
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.
Summary and Future Direction: Diffusion coefficients - Calculated and compared with experimental viscosity; More work on solid state diffusion coefficients. Heat capacity - Calculated and compared with Ce heat capacity; Need to calculate the solid heat capacity and find more experimental data. Heat of formation - Calculated for one and two Ce atoms; Slab calculations are under way for a surface. Reaction rate - Calculation of a single Ce atom reaction is in progress; Simulations of larger clusters for comparison with ionic reaction experiments. Emissivity and Thermal Conductivity - Single snapshot calculated; Averaging over ionic configurations at different T is in progress.
The experiments of Juanicotena and Szarynski, namely T101, T102, and T105 are modeled for purposes of gaining a better understanding of the FLAG friction model frictmultiscale2. This exercise has been conducted as a first step toward model validation. It is shown that with inclusion of the friction model in the numerical analysis, the results of Juanicotena and Szarynski are predicted reasonably well. Without the friction model, simulation results do not match the experimental data nearly as well. Suggestions for follow-on work are included.
Here, we report experimental results from an initial study of reactive and nonreactive metal fragments—ejecta—transporting in vacuum, and in reactive and nonreactive gases. We postulate that reactive metal fragments ejected into a reactive gas, such as H 2 , will break up into smaller fragments in situations where they are otherwise hydrodynamically stable in a nonreactive gas such as He. To evaluate the hypothesis we machined periodic perturbations onto thin Ce and Zn coupons and then explosively shocked them to eject hot, micron-scale fragments from the perturbations. The ejecta masses were diagnosed with piezoelectric pressure transducers, and their transport in H 2 and He was imaged with visible and infrared (IR) cameras. Because Ce + H 2 → CeH 2 + ΔH, where ΔH is the enthalpy of formation, an observed increase of the relative IR (radiance) temperature T R between the Ce–H 2 and Ce–He gas systems can be used to estimate the amount of Ce that converts to CeH 2 . As a result, the experiments sought to determine whether dynamic chemical effects should be included in ejecta-transport models.