Qu-antitative analysis of diffraction by liquids using a pink-spectrum X-ray source
Explore the source record for details and available documents.
Singh, Saransh↗
SEARCH · Engineering Papers
Results for “shock physics”
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.
Shock-driven amorphization and melting in Fe 2 O 3
We present measurements on Fe 2 O 3 amorphization and melt under laser-driven shock compression up to 209(10) GPa via time-resolved in situ x-ray diffraction. At 122(3) GPa, a diffuse signal is observed indicating the presence of a noncrystalline phase. Structure factors have been extracted up to 182(6) GPa showing the presence of two well-defined peaks. A rapid change in the intensity ratio of the two peaks is identified between 145(12) and 151(12) GPa, indicative of a phase change. The noncrystalline diffuse scattering is consistent with shock amorphization of Fe 2 O 3 between 122(3) and 145(12) GPa, followed by an amorphous-to-liquid transition above 151(12) GPa. Upon release, a noncrystalline phase is observed alongside crystalline α–Fe 2 O 3 . The extracted structure factor and pair distribution function of this release phase resemble those reported for Fe 2 O 3 melt at ambient pressure.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave
Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.
Cambier, Jean-Luc↗
Novel PBX formulations containing thermally-expandable microspheres for on-demand control of explosive behavior
In this work, we present the formulation and analysis of inert Plastic-Bonded eXplosive (PBX)-surrogates loaded with thermally expandable microspheres (TEMs). TEMs consist of a thermoplastic acrylonitrile shell (10-50 µm diameter) encapsulating an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol%. We hypothesize that TEM expansion within a PBX will introduce tunable changes in local density and porosity, and ultimately the shock sensitivity. This paper focuses on microstructural analysis of surrogate PBX-TEM formulations, both before and after thermally- induced TEM expansion. The formulations were analyzed by scanning electron microscopy (SEM), ultra-small angle neutron scattering (USANS), and X-ray Computed Tomography (x-ray CT) under various thermal conditions. Evidence suggests that TEM expansion within a PBX is a viable method to impart microstructural changes that are predicted to have a measurable, on-demand effect on the initiation sensitivity. Future work includes development of PBX-TEM formulations with HMX and shock initiation experiments.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
The synchrotron-maser theory of type II solar radio emission processes - The physical model and generation mechanism
A theory is proposed to explain the generation mechanism of type II solar radio bursts. It is suggested that the shock wave formed at the leading edge of a coronal transient can accelerate electrons. Because of the nature of the acceleration process, the energized electrons can possess a 'hollow-beam' type distribution function. When the electron beam propagates along the ambient magnetic field to lower altitudes and attains larger pitch angles, a synchrotron-maser instability can set in. This instability leads to the amplification of unpolarized or weakly polarized radiation. The present discussion incorporates a model which describes the ambient magnetic field and background plasma by means of MHD simulation. The potential emission regions may be located approximately, according to the time-dependent MHD simulation. Since the average local plasma frequency in the source region can be evaluated from the MHD model, the frequent drift associated with the radiation may be estimated. The result seems to be in good agreement with that derived from observations.
Wu, C. S.↗
Shatter cones formed in large-scale experimental explosion craters
In 1968, a series of 0.5-ton and 100-ton TNT explosion experiments were conducted in granitic rock near Cedar City, Utah, as part of a basic research program on cratering and shock wave propagation. Of special interest was the formation of an important type of shock metamorphic feature, shatter cones. A description is presented of the first reported occurrence of shatter cones in high explosion trials. A background to shatter cone studies is presented and attention is given to the test program, geology and physical properties of the test medium, the observed cratering, and the formational pressures for shatter cones. The high explosion trials conducted demonstrate beyond any doubt, that shatter cones can be formed by shock wave processes during cratering and that average formational pressures in these crystalline rocks are in the 20-60 kb range.
Roddy, D. J.↗
Formation of a spherical plasma liner for plasma-jet-driven magneto-inertial fusion
Plasma-jet-driven magneto-inertial fusion is an alternative approach to controlled nuclear fusion, which aims to utilize a line-replaceable dense plasma liner as a repetitive spherical compression driver. In this experiment, first measurements of the formation of a spherical argon plasma liner formed from 36 discrete pulsed plasma jets are obtained on the Plasma Liner Experiment. Properties including liner uniformity and morphology, plasma density, temperature, and ram pressure are assessed as a function of time throughout the implosion process and indicate an apparent transition from initial kinetic inter-jet interpenetration to collisional regime near stagnation times, in accordance with theoretical expectation. A lack of primary shock structures between adjacent jets during flight implies that arbitrarily smooth liners may be formed by way of corresponding improvements in jet parameters and control. The measurements facilitate the benchmarking of computational models and understanding the scaling of plasma liners toward fusion-relevant energy density.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Delayed Shock-Induced Dust Formation in the Dense Circumstellar Shell Surrounding the Type IIn Supernova SN 2010jl
The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shock waves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as few weeks postexplosion. This emission has been considered by some as evidence for the rapid formation of dust in the cooling postshock CSM. We investigate the physical processes that may inhibit or facilitate the formation of dust in the CSM. When only radiative cooling is considered, the temperature of the dense shocked gas rapidly drops below the dust condensation temperature. However, by accounting for the heating of the postshock gas by the downstream radiation from the shock, we show that dust formation is inhibited until the radiation from the shock weakens as it propagates into the less dense outer regions of the CSM. In SN2010jl, dust formation can therefore only commence after day ∼380. Only the IR emission since that epoch can be attributed to the newly formed CSM dust. Observations on day 460 and later show that the IR luminosity exceeds the UV-optical luminosity. The postshock dust cannot extinct the radiation emitted by the expanding SN shock. Therefore, its IR emission must be powered by an interior source, which we identify as the reverse shock propagating through the SN ejecta. IR emission before day 380 must therefore be an IR echo from preexisting CSM dust.
Sarangi, Arkaprabha↗
Neutrinos from type II supernovae - The first 100 milliseconds
The collapse of a 1.17 solar mass iron core is numerically followed through infall to 100 ms past core bounce, and the emergent neutrino spectra during each phase are highlighted. It is found that, even with fairly optimistic conditions for producing a strong, sustained core-bounce shock wave, the prompt shock stalls within 9 ms of core bounce at a radius of less than 250 km. It appears that a radical change in the character of the progenitor core or in our understanding of the relevant physics of stellar collapse is needed before the direct mechanism for type II supernovae can become viable. Expanding the number of neutrino types from one to six magnifies the debilitating effect of neutrino loss on shock propagation. At shock breakout, prompt bursts of all neutrino types are observed. The luminosities of the nonelectron types show a sudden turn-on in luminosity while that of the electron neutrinos steadily increases throughout infall as a result of accelerating electron capture.
Myra, Eric S.↗
Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression
We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)
Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.
Mathematics and Computing↗
Roadmap for the exposé of radiation flows (Xflows) experiment on $\mathrm{NIF}$
Here, the goal of the Xflows experimental campaign is to study the radiation flow on the National Ignition Facility (NIF) reproducing the sensitivity of the temperature (±8 eV, ±23 μm) and density (±11 mg/cc) measurements of the COAX platform. This new platform will enable future astrophysical experiments involving supernova shock breakout, such as Radishock on OMEGA-60, and stochastic media (such as XFOL on OMEGA). Greater energy and larger physical scale on NIF will enable a greater travel distance of radiation flow, higher density, and more manufacturable foams and enable exploration of a greater range of radiation behavior than achievable in the prior OMEGA experiments. This publication will describe the baseline configuration for the Xflows experimental campaign and the roadmap to achieve its primary objectives.
47 OTHER INSTRUMENTATION↗
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↗
Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms
Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.
0D simulation↗
A cancellation problem in hybrid particle-in-cell schemes due to finite particle size
The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the different spatial discretizations of ions as finite-spatial-size particles and electrons as a grid-based fluid can lead to significant numerical wave dispersion errors in the long wavelength limit (kd i «1, where k is the wavenumber and di is the ion skin-depth). The problem occurs when high-order particle-grid interpolations, or grid-based smoothing, spreads the electric field experienced by the ions across multiple spatial cells and leads to inexact cancellation of electric field terms in the total (ion + electron) momentum equation. Practical requirements on the mesh spacing Δx/d i are suggested to bound these errors from above. The accuracy impact of not respecting these resolution constraints is shown for a non-linear shock problem.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Spatiotemporal evolution of emission and absorption signatures in a laser-produced plasma
Here we report spatiotemporal evolution of emission and absorption signatures of Al species in a nanosecond (ns) laser-produced plasma (LPP). The plasmas were generated from an Inconel target, which contained ~0:4 wt. % Al, using 1064 nm, ≈6 ns full width half maximum pulses from an Nd:YAG laser at an Ar cover gas pressure of ≈34 Torr. The temporal distributions of the Al I (394.4 nm) transition were collected from various spatial points within the plasma employing time-of-flight (TOF) emission and laser absorption spectroscopy, and they provide kinetics of the excited state and ground state population of the selected transition. The emission and absorption signatures showed multiple peaks in their temporal profiles, although they appeared at different spatial locations and times after the plasma onset. The absorption temporal profiles showed an early time signature representing shock wave propagation into the ambient gas. We also used emission and absorption spectral features for measuring various physical properties of the plasma. The absorption spectral profiles are utilized for measuring linewidths, column density, and kinetic temperature, while emission spectra were used to measure excitation temperature. A comparison between excitation and kinetic temperature was made at various spatial points in the plasma. Our results highlight that the TOF measurements provide a resourceful tool for showing the spatiotemporal LPP dynamics with higher spatial and temporal resolution than is possible with spectral measurements but are difficult to interpret without additional information on excitation temperatures and linewidths. The combination of absorption and emission TOF and spectral measurements thus provides a more complete picture of LPP spatiotemporal dynamics than is possible using any one technique alone.
(140.3440) Laser-induced breakdown, (280.4788) Opt↗
Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver
A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗