Role of graphite crystal structure on the shock-induced formation of cubic and hexagonal diamond
Abstract not provided
Volz, Travis J. (ORCID:0000000182249368)↗
SEARCH · Engineering Papers
Results for “shock initiation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB
Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock
The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effect of initial porous density on the shocked equation of state of granular aluminum oxide (Al2O3)
Explore the source record for details and available documents.
Padgiotis, Athena [Sandia National Laboratories (S↗
Shock to detonation transition of pentaerythritol tetranitrate (PETN) initially pressed to 1.65 g/cm 3
A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. In this work, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Revealing the Formation of the <20 MeV Inner Proton Radiation Belt at L ∼ 2 During the 10–11 May 2024 Superstorm
A new proton radiation belt was identified during the geomagnetic superstorm of 10–11 May 2024. To investigate its origin, we use an MHD-test particle simulation to model solar energetic proton (SEP) trapping and the evolution of the initial trapped proton population during the storm. The simulation weights the injected SEP population with interplanetary proton measurements and the initial trapped proton radial profile with energy-resolved Weather System Follow-on-Microwave (WSF-M) proton data, enabling quantitative comparison with observations. Results show that SEP contributions are minimal; the new ∼2–18 MeV belt at L ∼ 2 forms mainly through redistribution and energization of the initial trapped population within hours of the CME shock, driven by subsequent electric field impulses, rather than the initial shock alone. These results provide quantitative evidence based on realistic, plasma-dependent field dynamics and significantly advance understanding of inner radiation belt formation.
79 ASTRONOMY AND ASTROPHYSICS↗
Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations
The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.
42 ENGINEERING↗
Probing the lattice structure of dynamically compressed and released single crystal iron through the alpha to epsilon phase transition
In this work, experiments using broadband Laue x-ray diffraction (XRD) were used to examine the lattice structure of dynamically compressed [100]-oriented single crystal iron samples at the Dynamic Compression Sector at the Advanced Photon Source. These experiments used 1 μm thick iron single crystal samples sandwiched between a polyimide ablator and a polycarbonate window. A 100 J, 10 ns duration laser pulse incident on the polyimide ablator was used to shock compress the iron samples to initial stresses greater than 25 GPa, exceeding the ~13 GPa alpha (body-centered-cubic or bcc structure) to epsilon (hexagonal-close-packed or hcp structure) phase transition stress. XRD measurements were performed at various times relative to the shock wave entering the iron sample: early times, <~150 ps while the initial shock waves propagated through the iron; intermediate times, after the iron equilibrated with the ablator and window reaching a plateau stress state (12 or 17 GPa) lasting several nanoseconds; and late times, during uniaxial strain release. The early time measurements show that in <~150 ps, the high-pressure hcp phase is relaxed with a c/a ratio of 1.61, contrary to previous laser shock experiments where a c/a ratio of 1.7 was inferred. In the plateau stress state and partially released states, XRD measurements showed that the hcp structure retained a c/a ratio of 1.61 with no observable changes in the microstructure. Upon stress release at ~1 GPa/ns release rate, the reverse phase transition (hcp to bcc) to the original single crystal orientation (implying a transformation memory effect) was observed to reach completion somewhere between 13 and 11 GPa, indicating little stress hysteresis under rapid uniaxial strain release. A similar memory effect for the reverse hcp to bcc transformation has been previously observed under hydrostatic compression. However, the bcc/hcp orientation relationships differ somewhat between dynamic and static compression experiments, implying that the transformation pathway under uniaxial dynamic strain differs from the Burgers mechanism.
36 MATERIALS SCIENCE↗
Shock-induced kinetics and cellular structures of liquid nitromethane detonation
Using a combination of high-speed diagnostics, optical pyrometry, velocimetry and video photography, we examine the spatiotemporal reaction kinetics of a prototypical high-energy explosive, liquid nitromethane (NM), as we drive this high explosive into detonation. The detonations were initiated by powerful shock waves whose durations (4 ns) were shorter than the characteristic time associated with the reaction (7 ns). Simple optical spectroscopy alone cannot characterize the kinetics with high time resolution because the reactants and products are flowing at a high velocity of ~6 km/s (6 µm/ns). Additionally, every detection volume behind the shock front contains material initiated over a range of times. High spatial resolution can be obtained by probing interfaces where the shock enters or breaks out of the NM. In addition, we obtain cellular patterns imprinted on the luminous shock front by the two-stage explosion in NM. These spatiotemporal patterns arise naturally as a result of the asymmetry produced by the moving shock front. In this way the shock front serves as a thin moving intrinsic optical gauge that reports and characterizes the two-stage NM explosion behind the front.
42 ENGINEERING↗
Formation of Point Shocks for 3D Compressible Euler
We consider the 3D isentropic compressible Euler equations with the ideal gas law. We provide a constructive proof of the formation of the first point shock from smooth initial datum of finite energy, with no vacuum regions, with nontrivial vorticity present at the shock, and under no symmetry assumptions . We prove that for an open set of Sobolev‐class initial data that are a small L ∞ perturbation of a constant state, there exist smooth solutions to the Euler equations which form a generic stable shock in finite time. The blowup time and location can be explicitly computed, and solutions at the blowup time are smooth except for a single point , where they are of cusp‐type with Hölder C 1/3 regularity. Our proof is based on the use of modulated self‐similar variables that are used to enforce a number of constraints on the blowup profile, necessary to establish global existence and asymptotic stability in self‐similar variables. © 2022 Wiley Periodicals LLC.
Mathematics↗
Compressed Expanded Natural Graphite (CENG) Processing for PCM Composites
The use of phase change materials (PCMs) in thermal energy storage applications has received considerable attention in recent decades. Organic PCMs are popular due to their high latent heat of fusion, noncorrosive properties, and relative stability over many charge and discharge cycles. A primary limitation of these materials is their low thermal conductivity. This has led researchers to develop various methods to increase thermal conductivity by seeding PCM with or impregnating them into more conductive materials. One method is to impregnate PCMs into compressed expanded natural graphite (CENG) matrices, which can improve thermal conductivity by a factor of 100. CENG matrices have received particular interest due to their low cost, high porosity, small (nano/micro) pore size, high pore density, high thermal conductivity, and ability to be compressed into many geometries. PCM/CENG matrix composites have been extensively studied; however, the effect that CENG processing has on PCM saturation and the overall matrix thermal conductivity has not been well investigated. This processing includes four major steps including graphite intercalation, thermal shock, compression, and PCM saturation. Intercalation involves soaking graphite flakes in sulfuric and/or nitric acid to intercalate the acids between the graphene layers. The graphite flakes are then subjected to a high-temperature thermal shock, during which the intercalated acid is gasified rapidly, pushing the graphene layers apart, resulting in accordion-shaped graphite "worms". The "worms" are then compacted to a desired bulk density and then soaked with molten PCM until fully saturated. The properties of the produced CENG matrix, and its ability to allow PCM permeation, are sensitive to the processing parameters, namely, the thermal shock temperature and exposure time, as well as the matrix apparent density or porosity. Here, we study the effect of the thermal shock conditions necessary to expand intercalated graphite flakes on PCM saturation and the expanded graphite's thermal conductivity and morphology. We found that the thermal shock temperature exhibits the greatest influence. At greater shock temperatures, SEM images showed that expanded graphite worms exhibited greater density of pores, thus increasing total surface area within the matrices. Increasing thermal shock temperature yielded greater overall PCM saturation, as well as an increased rate of saturation. Improvements in PCM saturation rate and overall saturation are obtained as the shock temperature is increased. Longer exposure to thermal shock also improves initial saturation rates and is beneficial if a shortened impregnation time is needed. Thermal shock conditions did not impact thermal conductivity; however, conductivity was largely affected by matrix porosity. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to that observed in previous studies.
compressed expanded natural graphite↗
High speed PLIF study of the Richtmyer–Meshkov instability upon re-shock
The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Post-detonation fireball modeling: Validation of freeze out approximations
A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.
97 MATHEMATICS AND COMPUTING↗
Radiation drive designed to extend the pressure ranges measured in Gbar equation of state experiments at the National Ignition Facility
We present the design and demonstration of a Shock-Strengthening hohlraum radiation temperature drive in the Gbar experimental platform at the National Ignition Facility intended to increase the pressure range measured in a single experiment. Previously published experiments by Döppner et al. measured the equation of state in polystyrene from 25 to 60 Mbar. Recent experimental data of the Shock-Strengthening drive initially demonstrated a much larger pressure range from 15 to 110 Mbar using the same peak radiation temperature and experimental platform. The Shock-Strengthening drive starts with a low temperature foot that launches a weak shock into the sample and is followed by a continuous increase in radiation temperature to strengthen the leading shock. The additional strengthening increases the pressure within the sample beyond what is achievable by convergence alone. Design features of the Shock-Strengthening drive and accompanying radiation hydrodynamics simulations are used to illustrate the method by which the pressure range is increased from previous experiments. This method of modifying the radiation temperature drive can be used on the Gbar platform to significantly increase the range for equation of state data collected in a single experiment for many materials.
Physics - Plasma physics↗
INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Initial exploration of a novel transient arrest system involving fuel heating
A preliminary analysis on a novel accident response system to diminish the severity of super- critical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compared to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 to 1.3 dollar, it was found that an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING (Presentation)
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock
The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗