Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Explosive burning”

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.

At least 37 records · Page 2

Physical Sciences Vistas: Issue 4 2021

In this issue of Physical Sciences Vistas we highlight examples of “excellence in nuclear security.” In this issue, highlights of our outstanding work supporting the Laboratory’s nuclear security enterprise include the following: The complete refurbishment of Sigma’s surface finishing lab. Pride in the lab’s revitalized condition on the part of the operations and R&D staff responsible for this complex operation is well justified. Development of advanced materials and processes enabling next-generation weapons designs, much of which has been facilitated by projects supported by the Laboratory Directed Research and Development Program. The use of proton radiography in illuminating elements important to validating high explosive burn models. The experimental series was a collaboration between Los Alamos researchers and staff at Lawrence Livermore and Sandia national laboratories and the United Kingdom’s Atomic Weapons Establishment. Deployment of a sophisticated remote handling unit to safely transfer highly radioactive materials required to support mission-critical projects. The skillful swap of an electrical transformer at the Los Alamos Neutron Science Center (LANSCE). The deliberate operation restored power to the accelerator in less than estimated time. A look at George Goff ’s role in translating fundamental science into applied solutions essential to the Lab’s nuclear security mission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-instability Supernovae

Stars with masses of 80–130 M ⊙ can encounter pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2–3. The light curve peaks at 1.9–2.1 × 10 43 erg s -1 for 50 days and then plateaus at 2–3 × 10 42 erg s -1 for 200 days, depending on the viewing angle. The presence of 12 C and 16 O and the absence of 28 Si and 56 Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe, such as Type IIne, not just PPI SNe.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma Rays as a Signature of r -process Producing Supernovae: Remnants and Future Galactic Explosions

We consider the question of whether core-collapse supernovae (CCSNe) can produce rapid neutron capture process (r-process) elements and how future MeV gamma-ray observations could address this. Rare types of CCSNe characterized by substantial magnetic fields and rotation, known as magnetorotational supernovae (MR-SNe), are theoretically predicted to produce these elements, although direct observational evidence is lacking. We suggest that this critical question be addressed through the study of some of the 11 CCSN remnants located within 10 kpc, as well as through the detection of gamma-ray emission from a future Galactic supernova. We use a two-dimensional MR-SN model to estimate the expected gamma flux stemming from nuclear decays in the range of a few tens of keV to a few MeV. Our results indicate that an observation of 126 Sn ( 126 Sb) in a remnant stands out as a signature of an r-process-producing supernova. Since the neutron-rich conditions that lead to the production of the r-process could also enhance the production of 60 Fe, the detection of substantial 60 Fe ( 60 Co) would be indicative of favorable conditions for the r-process. In the case of a future supernova explosion, when the evolution of the spectrum is studied over 10 days to a few years, a rich picture emerges. At various epochs, the second peak r-process isotopes such as 125 Sn, 131 I, 132 Te, 132 I, and 140 La produce signals that are not obscured by the gamma emission from explosive burning products and electron–positron annihilation. The weak r-process isotopes 95 Nb, 103 Ru, and 106 Rh also have periods of prominence.

Gamma-rays↗

Isotopic compositions of bismuth, lead, thallium, and mercury from mini r-processing

The yields of stable isotopes of Bi, Pb, Tl and Hg as well as yields of Pb-205 are calculated with a parametrized model for 'mini r-processing' in the Ne, O, C-rich zones of explosive burning in massive stars. The Pb isotopic compositions stand out by their comparatively low Pb-207 yields and by the fact that this r-process variant yields Pb-204 quite abundantly. The average Pb-205/Pb-204 yield ratio of 6.1 is the same order of magnitude as yield ratios deduced for s-processing. The Hg from this mini r-process looks like normal solar-system mercury, but with Hg-196 missing and the light s-isotopes A = 198, 199, 200 and 201 depleted (especially the odd-A species).

Heymann, D.↗

Analysis of variability in the burst oscillations of the accreting millisecond pulsar XTE J1814-338

The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the known spin frequency during Type I X-ray bursts. The properties of the burst oscillations reflect the nature of the thermal asymmetry on the stellar surface. We present an analysis of the variability of the burst oscillations of this source, focusing on three characteristics: fractional amplitude, harmonic content and frequency. Fractional amplitude and harmonic content constrain the size, shape and position of the emitting region, whilst variations in frequency indicate motion of the emitting region on the neutron star surface. We examine both long-term variability over the course of the outburst, and short-term variability during the bursts. For most of the bursts, fractional amplitude is consistent with that of the accretion pulsations, implying a low degree of fuel spread. There is however a population of bursts whose fractional amplitudes are substantially lower, implying a higher degree of fuel spread, possibly forced by the explosive burning front of a precursor burst. For the first harmonic, substantial differences between the burst and accretion pulsations suggest that hotspot geometry is not the only mechanism giving rise to harmonic content in the latter. Fractional amplitude variability during the bursts is low; we can only rule out the hypothesis that the fractional amplitude remains constant at the l(sigma) level for bursts that do not exhibit photospheric radius expansion (PRE). There are no significant variations in frequency in any of the bursts except for the one burst that exhibits PRE. This burst exhibits a highly significant but small (= 0.1Hz) drop in frequency in the burst rise. The timescale of the frequency shift is slower than simple burning layer expansion models predict, suggesting that other mechanisms may be at work.

Watts, Anna L.↗

Radioactive nuclei in the early Solar system: analysis of the 15 isotopes produced by core-collapse supernovae

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 and 100 Myr can be used to probe the origin of the Solar system. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: 26 Al, 36 Cl, 41 Ca, 53 Mn, 60 Fe, 92 Nb, 97 Tc, 98 Tc, 107 Pd, 126 Sn, 129 I, 135 Cs, 146 Sm, 182 Hf, and 205 Pb. Here, we probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25 M ⊙ , explosion energies between 3.4 × 10 50 and 1.8 × 10 52 erg and compact remnant masses between 1.5 and 4.89 M ⊙ . We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as 92 Nb and 97 Tc, are the most affected by the remnant mass, 92 Nb varying by five orders of magnitude. Isotopes synthesized primarily in explosive C-burning and explosive He-burning, such as 60 Fe, are most affected by explosion energies. 60 Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15 M ⊙ model. The final yield of each examined SLR is used to compare to literature models.

79 ASTRONOMY AND ASTROPHYSICS↗

A study of shock initiation experiments for the explosive PBX 9502 using three reactive burn models

Shock to detonation transition (SDT) experiments are essential in calibrating and validating reactive burn models for explosives. This work investigates the large collection of SDT test data for the explosive PBX 9502 at ambient temperature that was presented by Gustavsen, Sheffield, and Alcon [Journal of Applied Physics, 99, 114907 (2006)]. We first analyze the experimental data and compare two different methods of determining the shock transition time/distance (namely, the bilinear method and the single-curve method). This reveals some of the uncertainty in estimating shock transition points, which contributes to scatter in Pop-plot data. Next, we compare the WSD, AWSD, and SURFplus reactive burn models for a collection of approximately 20 experimental shots using the FLAG hydrocode. Error estimates are used to quantify how well each reactive burn model (and their respective parameter calibrations) performs at predicting the SDT process for a range of loading conditions. Additionally, the importance of mesh resolution and numerical dissipation in SDT simulations will be assessed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accuracy of High Explosive Data used for Reactive Burn Model Calibration

Reactive burn models for a high explosive (HE) have three components: an equation of state (EOS) for the products, an EOS for the reactants and a reaction rate. For the products, we discuss the uncertainty in data for the planar detonation state and its release isentrope. For the reactants, the uncertainty in the principal shock locus data is discussed. For the reaction rate, we discuss uncertainties in shock initiation data (distance-of-run to detonation as a function of shock pressure), and curvature effect data (detonation speed as a function of front curvature). An important source of uncertainty for calibrating the reaction rate for a plastic-bonded explosive (PBX) results from variations in heterogeneities; that is the mesoscale structure of HE grains, binder and porosity. This shows up in the scatter of both shock initiation data and curvature effect data with small changes in initial density and initial temperature, and systematic changes with PBX lot. As illustrative examples, we use data for the conventional HMX based PBX 9501 and the insensitive TATB based PBX 9502 since these explosives have been extensively studied and high quality data are readily available in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The Deflagration-to-Detonation Transition in Two Dimensions

The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figures 1 and 2 for an overview of the proceeding steps. It has been assumed in the past that the same mechanism which drives DDT in the one dimensional case will apply to both two and three dimensions. However, this has never been experimentally verified and it is possible that an entirely new mechanism drives DDT higher dimensions. In order to test this a series of two dimensional DDT tests were performed to help provide empirical evidence of the mechanism in higher dimensions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗

Supernovae, grains and the formation of the solar system

An investigation is conducted concerning the possibility that observed Mg-26 anomalies in meteorites may be related to a nucleosynthetic event which preceded the formation of the solar system by at most a few million years. The Al-26, which decayed to form the observed excess Mg-26, could have been produced in either explosive carbon burning or in a high temperature carbon burning shell source immediately preceding the explosion. The results of supernova grain condensation calculations are presented and related to the hypothesis that a 'last event' supernova was indeed related to the formation of the solar system and thus might have created the observed isotopic anomalies in magnesium, oxygen, neon, and xenon.

Lattimer, J. M.↗

AWSD Reactive Burn Model for High Explosive LX‐14

ABSTRACT The results of an Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow calibration for the HMX‐based high explosive LX‐14 are presented. The parameters in the AWSD model are calibrated to experimental thermodynamic and gas gun data and to computational results from thermochemical calculations. There is no experimental rate stick data available for LX‐14; therefore, scaled experimental results from other PBX‐based high explosives are used in the calibration to fill this gap in data. Strong agreement is observed between the calibrated AWSD model and experimental data for LX‐14, including validation data that were not used in the calibration procedure. The developed model more accurately describes experimental shock‐to‐detonation results compared to several other reactive flow models for LX‐14 from the literature. The presented results illustrate that the AWSD model is capable of quantitatively describing the reactive burn of LX‐14.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Level Structure of {sup 31}S Via {sup 32}S({ital p,d}) {sup 31}S

Background: Properties of proton-unbound S-31 states determine the P-30(p, gamma) S-31 reaction rate, which has a significant impact on explosive hydrogen burning in classical novae and type-I x-ray bursts. Despite several previous studies, uncertainties still remain with respect to the nuclear structure of S-31 near the proton threshold. Purpose: The level structure of S-31 has been presently investigated via a charged-particle spectroscopy experiment using the S-32(p, d) S-31 reaction. Method: Deuterons corresponding to S-31 excited states with 3.285 <= E-x <= 10.8 MeV were momentum analyzed via an Enge split-pole spectrograph at six laboratory angles between 10 degrees and 62 degrees. Differential cross sections of the S-32(p, d)S-31 reaction were measured at E-p = 34.5 MeV. Distorted-wave Born approximation calculations were performed to constrain the spin-parity assignments of several of the observed levels. Results: We have detected 72 excited states of S-31, out of which 17 are within the astrophysical region of interest corresponding to the temperature range of 0.1-1.5 GK. We have resolved the discrepancy in the spin and parity of an excited state with E-x = 6542 keV, showing that is it not J(pi) = 3/2(-), and therefore the contribution of this state to the P-30(p, gamma) reaction rate is likely much less significant than previously thought owing to the larger angular-momentum transfer required to populate this excited state. Moreover, our measurement results help consolidate the spin-parity assignments for the 6377 and 6636 keV states in S-31. Conclusions: This work presents the most comprehensive spin-parity assignments to date from a single-neutron transfer reaction on S-32 to S-31 excited states in the region between 6 to 7 MeV excitation energy. This region is significant for the determination of the P-30(p, gamma)S-30 reaction rate over the temperatures characteristic of explosive hydrogen burning in novae.

Setoodehnia, K.↗

Sensitivity of 44 Ti and 56 Ni Production in Core-collapse Supernova Shock-driven Nucleosynthesis to Nuclear Reaction Rate Variations

Recent observational advances have enabled high resolution mapping of 44 Ti in core-collapse supernova (CCSN) remnants. Comparisons between observations and models provide stringent constraints on the CCSN mechanism. However, past work has identified several uncertain nuclear reaction rates that influence 44 Ti and 56 Ni production in postprocessing model calculations. We evolved one-dimensional models of 15 M ⊙ , 18 M ⊙ , 22 M ⊙ , and 25 M ⊙ stars from zero age main sequence through CCSN using Modules for Experiments in Stellar Astrophysics and investigated the previously identified reaction rate sensitivities of 44 Ti and 56 Ni production. In this work, we tested the robustness of our results by making various assumptions about the CCSN explosion energy and mass cut. We found a number of reactions that have a significant impact on the nucleosynthesis of 44 Ti and 56 Ni, particularly for lower progenitor masses. Notably, the reaction rates ${}^{13}{\rm{N}}{(\alpha ,p)}^{16}{\rm{O}}$, ${}^{17}{\rm{F}}{(\alpha ,p)}^{20}\mathrm{Ne}$, ${}^{52}\mathrm{Fe}{(\alpha ,p)}^{55}\mathrm{Co}$, ${}^{56}\mathrm{Ni}{(\alpha ,p)}^{59}\mathrm{Cu}$, ${}^{57}\mathrm{Ni}{(n,p)}^{57}\mathrm{Co}$, ${}^{56}\mathrm{Co}{(p,n)}^{56}\mathrm{Ni}$, ${}^{39}{\rm{K}}{(p,\gamma )}^{40}\mathrm{Ca}$, ${}^{47}{\rm{V}}{(p,\gamma )}^{48}\mathrm{Cr}$, ${}^{52}\mathrm{Mn}{(p,\gamma )}^{53}\mathrm{Fe}$, ${}^{57}\mathrm{Co}{(p,\gamma )}^{58}\mathrm{Ni}$, and ${}^{39}{\rm{K}}{(p,\alpha )}^{36}\mathrm{Ar}$ are influential for a large number of model conditions. Furthermore, we found the list of influential reactions identified by previous postprocessing studies of CCSN shock-driven nucleosynthesis is likely incomplete, motivating future larger-scale sensitivity studies.

79 ASTRONOMY AND ASTROPHYSICS↗

The mean mass of the abundant cosmic-ray nuclei from boron to silicon at 1.2 GeV per atomic mass unit

Results are presented for an investigation of cosmic-ray isotopic composition in which a balloon-borne instrument was used to perform an isotopic analysis of the abundant elements from B through Si in the energy range from 1.0 to 1.6 GeV/amu. The technique employed, which depends on the effect of earth's magnetic field on incoming cosmic rays, and reduction of the observational data are described in detail. Mean masses of the elements B, C, N, O, Ne, Mg, and Si are extracted from normalized Cerenkov spectra and evaluated in terms of the 'leaky box' model of galactic propagation and confinement. The implications of the results for cosmic-ray sources are shown to be consistent with the nucleosynthesis of certain isotopes in the cosmic rays in explosive carbon burning under the conditions of peak temperature, density, and neutronization expected in an explosion of a massive highly evolved star. The possibility of alternate interpretations is noted.

Dwyer, R.↗

Low-energy 23 Al β-delayed proton decay and 22 Na destruction in novae

The radionuclide 22 Na is a target of γ-ray astronomy searches, predicted to be produced during thermonuclear runaways driving classical novae. The 22 Na (p, γ) 23 Mg reaction is the main destruction channel of 22 Na during a nova, hence, its rate is needed to accurately predict the 22 Na yield. However, experimental determinations of the resonance strengths have led to inconsistent results. In this Rapid Communication, we report a measurement of the branching ratios of the 23 Al β-delayed protons as a probe of the key 204-keV (center-of-mass) 22 Na (p, γ) 23 Mg resonance strength. We report a factor of 5 lower branching ratio compared to the most recent literature value. The variation in 22 Na yield due to nuclear data inconsistencies was assessed using a series of hydrodynamic nova outburst simulations and has increased to a factor of 3.8, corresponding to a factor of ≈ 2 uncertainty in the maximum detectability distance. Finally, this is the first reported scientific measurement using the Gaseous Detector with Germanium Tagging system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗