SEARCH · Engineering Papers
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.
This presentation provides an overview of Explosives and Detonators for the NNSA NSTE 270 Course. It includes an overview of how detonators and high explosives are considered as part of Nuclear Explosives Safety. Several accidents involving High Explosives are also discussed. This presentation was adapted from LA-CP-16-20438 and LA-UR-21-20246.
Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.
Response strategies implemented in the first few hours to days after a nuclear detonation on the US homeland may save more than 100,000 lives. Effective planning and pre-event capability development can increase the response efficacy and so the number of lives saved. During such an event, the US Department of Defense, Joint Task Force Civil Support (DoD JTF-CS) provides command and control for the DoD forces supporting civil authority response operations. Among other activities, JTF-CS expects to support FEMA as well as State and local agencies in shelter and evacuation response activities. JTF-CS currently uses a manual method to generate evacuation decision making products. This method compares the dose that would be (a) acquired remaining at a given location to that (b) acquired during evacuation along a small number of routes. In these products, evacuation may be warranted when the projected “remain” dose is greater than the “evacuate” dose. This method only provides a limited consideration of the protection buildings provide their occupants. To inform the JTF-CS nuclear detonation response, the DoD Defense Threat Reduction Agency (DTRA) has tasked Lawrence Livermore National Laboratory (LLNL) to design a set of planning and response products to support shelter-evacuation activities. This report documents these products to facilitate DTRA and JTF-CS planning efforts. To provide context, we also provide appendices that summarize US planning guidance response zones, early (<72 h) response strategies, and other key topics.
Engineering design--good design--is a blend of engineering and science; this is especially true for of explosive detonator design. Detonators are complex components with interrelated mechanical, electrical, and chemical attributes, all of which must be balanced to achieve a robust and reliable, working whole. This design guide presents an integrated, lifecycle approach to the design process, bringing together engineering design with science, emphasizing elegance and subtlety in the design process.
Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.
Soot is a powdery substance composed of solid residues that results from the incomplete combustion of hydrocarbons. Soot holds valuable information because its chemical makeup can provide answers for condensed carbon kinetics in HE detonations by quantitively identifying what carbon species remain. For carbon-rich explosives, this is particularly important because carbon kinetics is a major pathway for energy release. Current efforts have been able to qualitatively identify the carbons and confirm their presence, but we haven’t been able to successfully quantify them. For example, small angle x-ray diffraction has been fielded to qualitatively identify carbon phases as a function of time during a detonation (work by Trevor Willey and Mike Nielson). Other work (ice experiments) has focused on measuring the “effective” static cold curve equation of state (EOS) of the recovered carbon soots.
Experimentally measuring the pressure gain of a rotating detonation combustor through the method of equivalent available pressure is prone to many experimental uncertainties. A detailed analysis of sources of random and systematic uncertainties in the thrust and pressure gain is performed using data from a hydrogen/air-operated rotating detonation combustor. The combustor had an axial air inlet and a 50% exit constriction. The measured (negative) pressure gain values agree with the results in the literature for similar inlet-to-exit area ratios. The base drag correction on the centerbody was found to be the greatest source of uncertainty for the thrust. While thrust is the predominant source of random uncertainty in the pressure gain calculation, significant systematic errors can be committed due to the assumed exit Mach number. The impact of this assumption becomes more detrimental the closer to positive the measured pressure gain is. By experimentally evaluating the exit Mach number, the random and systematic uncertainty in pressure gain both decreased, indicating that this should be an essential step. A definitive demonstration of gain is challenging, given the assessed uncertainties, and recommendations are provided to increase the precision of the equivalent available pressure methodology.
To meet the challenges associated with performing mid-InfraRed (IR) imaging in Rotating Detonation Combustors (RDCs), a novel pulsed mid-IR UpConversion Imaging (UCI) diagnostic has been implemented. UCI is an alternative to direct mid-IR detection that uses nonlinear optical frequency mixing to shift mid-IR wavelengths carrying a target image to shorter wavelengths that can be imaged with high-performance silicon-based CCD/CMOS cameras. This approach offers several favorable properties including high spectral selectivity, high temporal resolution, and superior low-light detectivity. A hydrogen-air research RDC was operated with carbon dioxide addition to allow pulsed UCI imaging of mid-IR luminosity within the combustion channel from spontaneous thermal emissions. The resulting measurements demonstrate high spatiotemporal resolution capable of imaging small structures near the supersonically propagating detonation wave front. The results show how this technique can be used to observe sharp gradients and millimeter-scale structures in the high-temperature, high-pressure zones RDC flow fields.
One of the practical challenges of the studies of rotating detonation engines (RDEs) is the direct estimation of losses from experimental measurements. This study attempts at resolving this limitation by combining a reduced order model (ROM) of the detonation wave characteristics with a Markov chain Monte Carlo parameter estimation framework. The model considers simple deflagration losses and the overall impact of deflagration on RDE performance. To evaluate this model, a Markov Chain Monte Carlo (MCMC) sampling technique was applied to estimate the loss parameters within the model for a set of conditions operated in hydrogen-air over a range of mass flow rates and equivalence ratios. The MCMC parameter estimation framework allowed for the determination of a posterior distribution of the loss parameters for each test condition, an examination of the correlation between the loss parameters and measured performance metrics of the RDE, and an uncertainty propagation of these parameters. The predicted model loss parameters were then compared to indirect experimental measurements of the deflagration combustion fractions to evaluating the validity of the approach, and shed light on the benefits and drawbacks of the model, measurement techniques, and the estimation framework.
Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.
This paper focuses on the development and validation of a combustion model for Computational Fluid Dynamics (CFD) modeling of Rotating Detonation Engines. A zero-dimensional Partially Stirred Reactor (PaSR) with a detailed chemical kinetic mechanism for hydrogen and air is used to model turbulent combustion. The model is computationally efficient and is based on the notion of partial mixing at the sub-grid level with turbulent exchange between mixed and unmixed regions. The ability of the PaSR model to accurately represent both detonative and deflagrative combustion is assessed by validating the results against experimental data. The effects of mesh resolution on the solution are also studied in order to determine if a mesh independent solution is obtainable with the Large Eddy Simulation (LES) approach to modeling turbulence. A comparison is made between the PaSR model and simply ignoring turbulence chemistry interactions which assumes that all species are perfectly mixed at the sub-grid level.
An opto-thermal laser detonator uses resonantly absorptive tuned nano-material associated with secondary explosives for optical absorption and initiation by an integral laser diode. The opto-thermal laser detonator includes main explosive material; resonantly absorptive tuned nano-material; secondary explosive material, wherein the resonantly absorptive tuned nano-material and the secondary explosive material are associated to form associated material made of the resonantly absorptive tuned nano-material and the secondary explosive material; and a laser diode operatively connected to the associated material, wherein the laser diode initiates the associated material which in turn initiates the main explosive material.
This is an extended abstract being submitted for consideration that briefly describes the material that will be included in the final paper. This study will detail the design and testing of a water-cooled Rotating Detonation Combustor (RDC) that permits operation for extended periods of time ensuring that the device has reach a stable operating temperature. Extended run times also provides an opportunity to consider transient behaviors that occur as a result of altering the operating conditions such as equivalence ratio. The experimental setup is also unique in that it consists of a ducted exhaust with a downstream high-temperature valve that can control the pre-combustion pressure in the RDC independent of the combustor annulus and exit geometry. This paper will examine the results from the extended operation of the water-cooled RDC over a range of both transient and steady state equivalence ratios (0.5 - 1.0), pre-combustion operating pressures (0 - 207 kPa), mass flow rates (0.42 – 0.5 kg/sec) and air inlet throat area to combustor channel area ratios (0.09 – 0.32). All tests were performed while operating on hydrogen in air at a nominal air inlet temperature of 340 – 350 K. Limited data obtained at an elevated inlet air temperature of 480 K will also be discussed. Results suggest that under certain conditions the operating mode of the detonation wave can vary as it approaches a stable operating temperature. For the RDC tested, a stable operating temperature is not achieved until approximately 8-10 seconds of run time.
An opto-thermal laser detonator uses resonantly absorptive tuned nano-material associated with secondary explosives for optical absorption and initiation by an integral laser diode. The opto-thermal laser detonator includes main explosive material; resonantly absorptive tuned nano-material; secondary explosive material, wherein the resonantly absorptive tuned nano-material and the secondary explosive material are associated to form associated material made of the resonantly absorptive tuned nano-material and the secondary explosive material; and a laser diode operatively connected to the associated material, wherein the laser diode initiates the associated material which in turn initiates the main explosive material.
This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.
Imaging of Bubble Collapse Effects in Optically Transparent High Explosive as a Method to Study the Detonation Process STL-017-19, Year 2. Presentation to be given a the FY 2020 SDRD Annual Meeting, to be held Sept. 23-24, on WebEx.
A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.
We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.