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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Isotopic anomalies from neutron reactions during explosive carbon burning

The possibility that the newly discovered correlated isotopic anomalies for heavy elements in the Allende meteorite were synthesized in the secondary neutron capture episode during the explosive carbon burning, the possible source of the O-16 and Al-26 anomalies, is examined. Explosive carbon burning calculations under typical conditions were first performed to generate time profiles of temperature, density, and free particle concentrations. These quantities were inputted into a general neutron capture code which calculates the resulting isotopic pattern from exposing the preexisting heavy seed nuclei to these free particles during the explosive carbon burning conditions. The interpretation avoids the problem of the Sr isotopic data and may resolve the conflict between the time scales inferred from 1-129, Pu-244, and Al-26.

Lee, T.↗
NASA NTRS1981

Nucleosynthesis of neutron-rich heavy nuclei during explosive helium burning in massive stars

The production of heavy nuclei during explosive helium burning has been calculated using a hydrodynamical model of a 15-solar mass (Type II) supernova and an n-process nuclear reaction network. It is found that the resulting neutron-rich heavy nuclei are not produced in the relative abundances of solar-system r-process material, especially in the vicinity of Pt, nor are any actinides produced. These deficiencies reflect an inadequate supply of neutrons. However, some neutron-rich isotopes, normally associated with the r-process, are produced which may be significant for the production of isotopic anomalies in meteorites.

Blake, J. B.↗
NASA NTRS1981

The elements just beyond iron - Formation during explosive carbon burning

It is pointed out that the details of the synthesis of the elements just beyond iron are unclear. Thus, a reexamination of the role played by neutron-capture reactions during explosive carbon burning is indicated. A description is given of calculations of explosive carbon-burning nucleosynthesis using a complete neutron capture (n-process) computer code with a network extending from Cr through Zr and including not only (n,gamma), (gamma,n) reactions and beta-decay but also (p,n), (n,p), (p,gamma), and (p,alpha) reactions. Initial conditions indicative of the composition of a massive star following core helium burning were selected, and a comparison is made with conditions similar to those used by Howard et al. (1972). It is found that neutron reactions during explosive carbon burning are an important source for the elements just beyond iron.

Wefel, J. P.↗
NASA NTRS1980

Nucleosynthesis of neutron-rich heavy nuclei during explosive helium burning in a 15 solar-mass supernova

The production of heavy nuclei during explosive helium burning has been calculated using the Weaver and Woosley self-consistent model of a complete 15 solar-mass star and the n-process code of Blake and Schramm. It was found that the resulting neutron-rich heavy nuclei are not produced in the relative abundances of solar-system r-process material (such as a Pt peak) nor are any actinides produced. Basically insufficient neutrons are available.

Blake, J. B.↗
NASA NTRS1979

Isotopic anomalies from neutron reactions during explosive carbon burning

The heavy isotopic anomalies observed recently in the fractionation and unknown nuclear inclusions from the Allende meteorite are explained by neutron reactions during the explosive carbon burning (ECB). This model produces heavy anomalies in the same zone where Al-26 and O-16 are produced, thus reducing the number of source zones required for the isotopic anomalies. Unlike the classical r-process, the ECB n-process avoids the problem with the Sr anomaly and may resolve the problem of conflicting time scales between Al-26 and the r-process isotopes I-129 and Pu-244. Experimental studies of Zr and Ce isotopic composition are proposed to test this model.

Lee, T.↗
NASA NTRS1980

Synthesis of Al-26 in explosive hydrogen burning

The possibility of Al-26 synthesis during the explosive processing of hydrogen-rich material in the outer layers of a supernova or in nova envelopes is investigated. It is found that in the peak temperature range of 1-3 x 10 to the 8th deg K and for expansion time scales of the order of 1-1000 s, values of (Al-26)/(Al-27) as high as 0.1-1 can be obtained for values of 0.001-100,000 for the product of the peak density in g/cu cm and the proton mass fraction. Such a level of Al-26 production is considerably higher than that of recent carbon/neon burning nucleosynthesis models, and is sufficient to account for the magnesium anomalies detected in certain inclusions from the Leoville and Allende meteorites. Al abundances resulting from a hydrodynamical calculation performed on the grounds of a 25 solar mass presupernova model are also presented, and the influence of uncertainties in the input physics is discussed; in particular, the rates of the (Mg-25)(p,gamma)(Al-26) and (Al-26)(p,gamma)(Si-27) key reactions are studied.

Arnould, M.↗
DOE OSTI2020

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base

Reactive burn models represent a significant leap in high explosive (HE) modeling capability. The first generation of engineering models of HE detonation are called programmed burn models and they are largely based on the distance between a prescribed detonation point and each zone in a simulation. There have been many advancements to programmed burn models over the years and when the assumptions upon which they are based are met, a properly tuned programmed burn model can be highly accurate but if any of their assumptions is not met, as is the case for corner turning or weakly initiated HE burn, they will give the wrong answer. Reactive burn models represent an entirely new way of modeling HE burn. They use the local conditions of a zone – e.g. temperature, pressure or density – as calculated by a hydrocode to determine if and when the zone is going to detonate and if so, how rapidly. This difference opens up an entirely new set of capabilities for HE modeling. It makes it possible to accurately and predictively model phenomena like the effect of confinement and the formation of dead zones. Reactive burn models have seen sustained development effort at LANL for at least the last decade but several recent developments make it timely to transition reactive burn models from a research topic to a production tool. The main goal of this milestone is to facilitate and accelerate the adoption of reactive burn as a commonly available modeling option, with recommendations on the resolution that will be required and uncertainties associated with their modeling choices. To achieve this, we have performed verification, validation, and uncertainty quantification (UQ) assessments of AWSD and SURF/SURFplus in xRage and FLAG on a variety of different problems.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base (FY2020 L2 Milestone MRT 7129) [Slides]

Advanced reactive burn models enable phenomena to be modeled that aren’t possible with programmed burn. Programmed burn propagates a burn front at a prescribed speed from a prescribed initiation point. Reactive burn uses the hydro variables (pressure, temperature, etc.) to determine when, how rapidly, and to what extent a zone detonates. This enables a completely new set of problems to be simulated. Our goal in this milestone has been to improve the usability of reactive burn modes. Thus, the closure criteria are: 1. Perform code and solution verification of the reactive burn models using appropriate solutions and simplified HE setups. 2. Perform validation assessments using small-scale experiments. 3. Evaluate sensitivities and quantify uncertainties due to model form variations and mesh dependencies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Detonation Wave Profile (update of LA-UR-15-29498-rev1)

The Zel’dovich-von Neumann-Doering (ZND) profile of a detonation wave is derived. Two basic material properties are required: i. An equation of state (EOS) for a partly burned explosive; P(V, e, λ). ii. A burn rate for the reaction progress variable; d/dtλ = R(V, e, λ). For a steady planar detonation wave the reactive flow PDEs can be reduced to ODEs. The detonation wave profile can be determined from an ODE plus algebraic equations for points on the partly burned detonation loci with a specified wave speed. Furthermore, for the CJ detonation speed, the end of the reaction zone is sonic. A solution to the reactive flow equations can be constructed with a rarefaction wave following the detonation wave profile. This corresponds to an underdriven detonation wave, and the rarefaction is know as a Taylor wave.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1994

Convection, nucleosynthesis, and core collapse

We use a piecewise parabolic method hydrodynamics code (PROMETHEUS) to study convective burning in two dimensions in an oxygen shell prior to core collapse. Significant mixing beyond convective boundaries determined by mixing-length theory brings fuel (C-12) into the convective regon, causing hot spots of nuclear burning. Plumes dominate the velocity structure. Finite perturbations arise in a region in which O-16 will be explosively burned to Ni-56 when the star explodes; the resulting instabilities and mixing are likely to distribute Ni-56 throughout the supernova envelope. Inhomogeneities in Y(sub e) may be large enough to affect core collapse and will affect explosive nucleosynthesis. The nature of convective burning is dramatically different from that assumed in one-dimensional simulations; quantitative estimates of nucleosynthetic yields, core masses, and the approach to core collapse will be affected.

Bazan, Grant↗
DOE OSTI2020

Study of the ${{^{25}}}$Mg(d,p)${{^{26}}}$Mg reaction to constrain the ${{^{25}}}$Al(p,${\gamma }$)${{^{26}}}$Si resonant reaction rates in nova burning conditions

The rate of the 25 Al( p,γ) 26 Si reaction is one of the few key remaining nuclear uncertainties required for predicting the production of the cosmic γ-ray emitter 26 Al in explosive burning in novae. This reaction rate is dominated by three key resonances (J π = 0 + , 1 + and 3 + ) in 26 Si. Only the 3 + resonance strength has been directly constrained by experiment. A high resolution measurement of the 25 Mg(d, p) reaction was used to determine spectroscopic factors for analog states in the mirror nucleus, 26 Mg. A first spectroscopic factor value is reported for the 0 + state at 6.256 MeV, and a strict upper limit is set on the value for the 1 + state at 5.691 MeV, that is incompatible with an earlier ( 4 He, 3 He) study. These results are used to estimate proton partial widths, and resonance strengths of analog states in 26 Si contributing to the 25 Al( p,γ) 26 Si reaction rate in nova burning conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
NASA NTRS1994

Evidence for Ni-56 yields Co-56 yields Fe-56 decay in type Ia supernovae

In the prevailing picture of Type Ia supernovae (SN Ia), their explosive burning produces Ni-56, and the radioactive decay chain Ni-56 yields Co-56 yields Fe-56 powers the subsequent emission. We test a central feature of this theory by measuring the relative strengths of a (Co III) emission feature near 5900 A and a (Fe III) emission feature near 4700 A. We measure 38 spectra from 13 SN Ia ranging from 48 to 310 days after maximum light. When we compare the observations with a simple multilevel calculation, we find that the observed Fe/Co flux ratio evolves as expected when the Fe-56/Co-56 abundance ratio follows from Ni-56 yields Co-56 yields Fe-56 decay. From this agreement, we conclude that the cobalt and iron atoms we observe through SN Ia emission lines are produced by the radioactive decay of Ni-56, just as predicted by a wide range of models for SN Ia explosions.

Kuchner, Marc J.↗
DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗
DOE OSTI2021

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↗