Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical explosions”

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 127 records · Page 7

Liquid rocket actuators and operators

All the types of actuators and associated operators used in booster, upper stage, and spacecraft propulsion and reaction-control systems except for chemical-explosive actuators and turbine actuators are discussed. Discussion of static and dynamic seals, mechanical transmission of motion, and instrumentation is included to the extent that actuator or operator design is affected. Selection of the optimum actuator configuration is discussed for specific application which require a tradeoff study that considers all the relevant factors: available energy sources, load capacity, stroke, speed of response, leakage limitations, environmental conditions, chemical compatibility, storage life and conditions, size, weight, and cost. These factors are interrelated with overall control-system design evaluations that are beyond the scope of this monograph; however, literature references are cited for a detailed review of the general considerations. Perinent advanced-state-of-the-art design concepts are surveyed briefly.

Source record↗

Rugged nanoparticle tracers for mass tracking in explosive events

Tracing the flow of solid matter during an explosion requires a rugged tag that can be measured by a unique identifiable signature. Silica-covered semiconductor quantum dots (QDs) provide a unique and tunable photoluminescent signature that emits from within a sacrificial outer layer. Five types of silica-covered zinc sulfide QDs were synthesized and covalently bound to commercial luminescent powders. The combination of five dots and five powders enables a matrix of 25 unique tags. The tracers are shown to be tolerant of environments associated with chemical explosives and provides a unique tag to evaluate debris fields.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The Low-Yield Nuclear Monitoring (LYNM) Experimental Science Plan

The Low-Yield Nuclear Monitoring (LYNM) Program is a long-term NNSA research and development effort designed to improve the United States’ explosion monitoring capabilities, particularly with respect to low-yield and potentially evasive underground nuclear testing. The LYNM Program focuses on researching, discovering, and exploiting unique and useful signatures, from all available technologies and sensors (e.g., seismic, acoustic, electromagnetic, gases, and particulates (both stable and radioactive)). Four Department of Energy laboratories participate in LYNM and together are referred to as the ‘quad-lab’. The R&D program execution is performed under NNSA defined structures known as “ventures”. Four science ventures were established: 1) Explosion Source Functions; 2) Containment of Low-Yield Underground Tests; 3) Local Signatures; 4) Dynamic Monitoring Networks. To organize and execute the large field scale experiment a fifth venture was established: 5) Physics Experiment One (PE-1). As part of the LYNM Program, a series of experiments are planned at various scales and levels of venture involvement. These vary from those that involve a subset of labs and/or LYNM ventures (e.g., small experiments), to full quad-lab LYNM Program field-scale integrated experiments. Such experiments may involve chemical explosions with tracer materials or other means of simulating the expected signals from a nuclear explosion. The LYNM Program does not conduct actual nuclear explosions. Since 1992, the U.S. has observed a moratorium on underground nuclear explosions. This document is intended to provide the underlying scientific basis for the LYNM planned experimental work. Each specific LYNM experiment will develop a goals, objectives, and requirements (GOR) plan following the guidance in this document. The LYNM technical staff will define the numbers and types of experiments required over the course of the Program based on technical needs and within funding constraints. As with any scientific experiment series, the number and types of experiments may change based upon the experimental results obtained. An experiment that agrees with models/codes/software signature predictions may need fewer repetitions/variations, depending upon the level of statistical rigor desired, as compared to one in which the predictions and experimental data do not match. The large LYNM field-scale integrated experiments require the longest lead-time for planning, and these are discussed in more detail near the end of this document.

58 GEOSCIENCES↗

Finite-Difference Simulation for Infrasound Generated by Finite-Extent Ground Motions

Underground explosions can produce infrasound in the atmosphere, and the wavefield characteristics are often governed by the ground surface motions. Finite-difference methods are popular for infrasound simulation as their generality and robustness allow for complex atmospheric structures and surface topography. A simple point-source approximation is often used because infrasound wavelengths tend to be large relative to the source dimensions. However, this assumption may not be able to capture the complexity of explosion-induced ground motions if the surface area is not compact, and appropriate source models must be incorporated into the finite-difference simulations for accurate infrasound prediction. In this study, we develop a point source representation of the complex ground motions for infrasound sources. Instead of a single point source, we use a series of point sources distributed over the source area. These distributed point sources can be equivalent to air volume changes produced by the ground motions in the atmosphere. Further, we apply the distributed point-source method to a series of buried chemical explosions conducted during the Source Physics Experiment Phase I. Epicentral ground-motion measurements during the experiments provide a way to calculate accurate distributed point sources. We validate and evaluate the accuracy of distributed point source approach for infrasound simulations by direct comparison with acoustic observations in the field experiment.

58 GEOSCIENCES↗

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational diagnostics and characterization of combustion recession in diesel sprays

While low-temperature combustion (LTC) strategies have been found to mitigate nitrogen oxides and particulate matter emissions in diesel engines, studies have also associated LTC with an increase in unburnt hydrocarbons. With more recent studies on diesel after end-of-injection (AEOI), combustion recession is identified as a phenomenon where at near nozzle region, high-temperature ignition (HTI) combustion can propagate back to the nozzle tip consuming the unburnt hydrocarbons AEOI. Current literature has suggested that combustion recession is controlled by auto-ignition. However, high-fidelity simulations and detailed analysis of such a mechanism are missing. In this study, comprehensive Large Eddy Simulations of a reacting spray at “Spray A” conditions are performed, where detailed analysis of combustion recession concerning flame morphology and propagation modes are included. In particular, this study demonstrated for the first time that while combustion recession is mainly auto-ignition dominated (consistent with the literature), a cool flame was found to deflagrate towards the richer regions of the mixture, promoting mixing and increasing the mixture temperature. This leads to HTI kernels, which then grow and develop as deflagrative waves, therefore sustaining the combustion recession process. The study also detailed the extinction mechanism of combustion: the entrainment wave will overlean the near-nozzle mixtures, rendering it unable to support HTI, which leads to the extinction of the upstream flame AEOI in lower reactivity mixtures. Combustion recession is also observed to be contingent on the chemical and diffusion processes, even at low scalar dissipation rates. Finally, a new criterion for combustion recession based on chemical explosive mode is proposed and validated with previous combustion recession index to quantify the extent of HTI in near-nozzle mixtures AEOI. The newly developed metric combined with a previous experimentally-based metric can provide simple but valuable measurements of the degree and propensity of the upstream flame AEOI.

CEMA↗

Analysis of 127 Xe tracer measurements using a net counts method

A suite of measurement systems were deployed as part of the Physical Experiment 1 series of experiments, which involved detonating chemical explosives along with radionuclide tracers in an underground cavity, at the Nevada National Security Site (NNSS) in the United States. One of the radionuclide tracers, 127 Xe was released from the containment following the explosion and detected on a SAUNA Q B sampler situated approximately 3.5 km away. The system uses a beta-gamma coincidence detector system to measure fission product radioisotopes of xenon relevant to nuclear explosion monitoring. In this work we use the coincidence measurement data to analyse and interpret the results from the SAUNA Q B system, to calculate the measured 127 Xe activity concentration(s).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DNS Of the ignition process of n-heptane/air premixed combustion with low-temperature chemistry in turbulent boundary layer

In the present work, three-dimensional direct numerical simulation (DNS) of n-heptane/air premixed combustion in turbulent boundary layer was performed to explore the near-wall ignition process with low -temperature chemistry. A reduced chemical mechanism with 58 species and 387 elementary reactions for n-heptane combustion was used in the DNS. The general characteristics of the ignition process near the wall were examined. Here, it was found that low-temperature ignition (LTI) dominates the upstream region, and high -temperature ignition (HTI) appears in the downstream region. The ignition process and the low-temperature chemistry pathways of the DNS are compared with those of a corresponding laminar case. It was found that the ignition process was affected by turbulence, which results in thickened reaction zones. However, the carbon flow analysis of low-temperature chemistry showed that turbulence rarely affects the low-temperature chemistry pathway. The combustion modes of various regions were scrutinized based on the budget terms of species transport equations and the chemical explosion mode analysis (CEMA). It was shown that the reaction term of RO 2 is significant during the LTI process of the upstream region, and the reaction terms of CH 2 O and CO 2 are evident in the downstream region, indicating the occurrence of HTI. It was also shown that auto-ignition is dominant in the upstream region. With increasing streamwise distance, the contribution of flame propagation increases, which takes over that of auto-ignition in the near-wall region.

33 ADVANCED PROPULSION SYSTEMS↗

Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase

Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topographically Scattered Infrasound Waves Observed on Microbarometer Arrays in the Lower Stratosphere

When an acoustic wave strikes a topographic feature, some of its energy is scattered. Sensors on the ground cannot capture these scattered signals when they propagate at high angles. We report observations of upwardly-scattered acoustic waves prior to refraction back to the ground, intercepting them with a set of balloon-borne infrasound microbarometers in the lower stratosphere over northern Sweden. We show that these scattered waves generate a coda whose presence can be related to topography beneath balloons and low-altitude acoustic ducts. The inclination of the coda signals changes systematically with time, as expected from waves arriving from scatterers successively closer to receivers. The codas are present when a temperature inversion channels infrasound from a set of ground chemical explosions along the ground, but are absent following the inversion's dissipation. Since scattering partitions energy away from the main arrival, these observations imply a mechanism of amplitude loss that had previously been inaccessible to measurement. As such, these data and results allow for a better comprehension of interactions between atmospheric infrasound propagation and the solid earth

58 GEOSCIENCES↗

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

Infrasound direction of arrival determination using a balloon-borne aeroseismometer

Free-floating balloons are an emerging platform for infrasound recording, but they cannot host arrays sufficiently wide for multi-sensor acoustic direction finding techniques. Because infrasound waves are longitudinal, the balloon motion in response to acoustic loading can be used to determine the signal azimuth. This technique, called “aeroseismometry,” permits sparse balloon-borne networks to geolocate acoustic sources. This is demonstrated by using an aeroseismometer on a stratospheric balloon to measure the direction of arrival of acoustic waves from successive ground chemical explosions. A geolocation algorithm adapted from hydroacoustics is then used to calculate the location of the explosions.

47 OTHER INSTRUMENTATION↗

Luminescent silica microagglomerates, synthesis, and environmental testing

Abstract Tracking mass through harsh environments requires surrogate particles that withstand the event and endure until sampling. Silica-covered quantum dots have been shown to withstand a range of environmental pHs from months to years; in this work they are shown to endure in anticipated local environments. Two methods of particle synthesis were employed to produce luminescent silica with particle diameters 0.1–4 μm. These tracer particles scale for mass production, tolerate harsh environments, and endure in debris. They could be deployed in places such as chemical explosions, industrial processes, geologic test beds, oil and gas fields, nuclear reactors, and geothermal plants to track mass under harsh conditions. Graphical abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effects of Atmospheric Models on the Estimation of Infrasonic Source Functions at the Source Physics Experiment

We invert infrasound signals for an equivalent seismoacoustic source function using different atmospheric models to produce the necessary Green’s functions. The infrasound signals were produced by a series of underground chemical explosions as part of the Source Physics Experiment (SPE). In a previous study, we inverted the infrasound data using so-called predictive atmospheric models, which were based on historic, regional-scaled, publicly available weather observations interpolated onto a 3D grid. For the work presented here, we invert the same infrasound data, but using atmospheric models based on weather data collected in a time window that includes the approximate time of the explosion experiments, which we term postdictive models. We build two versions of the postdictive models for each SPE event: one that is based solely on the regional scaled observations, and one that is based on both regional scaled observations combined with on-site observations obtained by a weather sonde released at the time of the SPE. We then invert the observed data set three times, once for each atmospheric model type. We find that the estimated seismoacoustic source functions are relatively similar in waveform shape regardless of which atmospheric model that we used to construct the Green’s functions. However, we find that the amplitude of the estimated source functions is systematically dependent on the atmospheric model type: using the predictive atmospheric models to invert the data generally yields estimated source functions that are larger in amplitude than those estimated using the postdictive models.

54 ENVIRONMENTAL SCIENCES↗

How Good Is Your Location? Comparing and Understanding the Uncertainties in Location for the 1993 Rock Valley Sequence

Accurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.

58 GEOSCIENCES↗

Inverting infrasound data for the seismoacoustic source time functions and surface spall at the Source Physics Experiments Phase II: Dry Alluvium Geology

This report presents the infrasound data recorded as part of the Source Physics Experiment - Phase 2, Dry Alluvium Geology. This experiment, also known colloquially as DAG, consisted of four underground chemical explosions at the Nevada National Security Site. We focus our analysis on only the fourth explosion (DAG-4) as we determined that this was the only event that produced clear source-generated infrasound energy as recorded by the DAG sensors. We analyze the data using two inversion methods. The first method is designed to estimate the point-source seismoacoustic source time functions, and the second inversion method is designed to estimate the first-order characteristics (e.g. horizontal dimensions and maximum amplitude) of the actual spall surface. For both analysis methods, we are able to fit the data reasonably well, with various assumptions of the source model. The estimated seismoacoustic source appears to be a combination of a buried, isotropic explosion with a maximum amplitude of ~2 x 10 9 Nm and a vertically oriented force, applied to the Earth's surface with a maximum amplitude of 4 x 10 7 N. We use the vertically oriented force to simulate surface spall. The estimated spall surface has an approximate radius of ~40 m with a maximum acceleration magnitude in the range of 0.8 to 1.5 m/s/s. These estimates are approximately similar to the measured surface acceleration at the site.

58 GEOSCIENCES↗