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 thesis defense covers Background/Motivation, High-Temperature SDT Measurements of PBX 9502, Neutron Scattering Measurements of PBX 9502, Reactive Burn Modeling of PBX 9502 Using SURF, Conclusion, and Future Work.
Accurate predictive modeling of high explosives in abnormal conditions, such as a fire, is critical to personnel safety. Modeling depends heavily on precise physical characterization of the high explosive in question. The aim of this thesis was to characterize the detonation response to shock and the microstructure at high-temperature (250°C) of PBX 9502 and use these data to attempt to model these results using a new reactive burn model called SURF (Scaled Uniform Reactive Front). In order to characterize the shock response of PBX 9502 at high-temperature a 1D gas gun experiment was designed and executed that provided shock to detonation transition data. These data were then used to do a preliminary calibration of SURF for PBX 9502 at 250°C. Small-angle neutron scattering (SANS) was employed to characterize the void morphology of PBX 9502 as a function of temperature. This was done in order to feed a modified version of SURF called physically informed - SURF (π-SURF) which uses the void size and quantity distribution data as measured by SANS along with the deflagration and detonation characteristics of a high explosive to inform the SURF model. The π-SURF model was then used to model both ambient and high-temperature PBX 9502 SDT response and the results were compared to the physically measured results. It was found that π-SURF works in ambient cases. At high-temperature, however, it was determined that not enough information is known about hot spot formation to use the π-SURF model for high-temperature scenarios.
The kinetics of TATB decomposition have been investigated by thermal analysis (thermogravimetric analysis, TGA, and differential scanning calorimetry, DSC) at both LLNL and the University of Alabama Birmingham. This report summarizes progress in understanding from those results within the context of identifying the global reaction mechanism, i.e., sequential, autocatalytic, competitive, or some mixture thereof. A particularly important issue is how pressure affects both the reaction rate and enthalpy. The effect of pressure helps identify dominant mechanisms and construct and constrain a global mechanism that will work over a wide range of conditions. This memorandum is simply a snapshot of current progress and not a definitive interpretation. I have written it to organize work of others and my interpretations in a systematic manner with the intention of stimulating further advances in understanding. SDT data comes from Jen Montgomery at LLNL, and high-pressure DSC comes from Victoria Stanford and Sergey Vyazovkin at University of Alabama at Birmingham. Major observations and conclusions, some new and some not, are (1) outgassing below 150 °C prior to decomposition releases adsorbed water vapor, (2) outgassing between 150 and 250 °C may come primarily from intercalated water (3) initial thermal decomposition above 250 °C or so involves dehydration of TATB to create furazans, (4) quickly heating to higher temperatures in open pans with flowing gas leads primarily to TATB evaporation between 270 and 390 °C, (5) using a pinhole crucible at atmospheric pressure largely eliminates TATB sublimation, (6) at typical thermal analysis heat rates, the main stage of decomposition is described by a single mass loss process, and decomposition of residual material is spread over a wide range of thermal severity, (7) at typical thermal analysis heat rates and constant temperatures above 300 °C, heat release occurs in two overlapping peaks, with the first one a little faster than mass loss and the second on a little slower than mass loss, (8) total heat release increases significantly with heating rate and substantially with increased pressure, and (9) isothermal experiments at 340 °C indicate that the first exotherm increases by more with LX-17 than TATB, which suggests a possible chemical role of Kel-F decomposition.
Exploding bridgewire (EBW) detonators have been used in weapon systems since the 1940s but there is huge debate surrounding how energy is transferred throughout the EBW firing system and the mechanism by which the exploding wire leads explosive detonation. This report summarizes the underpinning technologies and physical processes that are currently understood and reviews the various efforts to quantify the mechanism by which the PETN is initiated. The behavior of the firing system is very well understood and predictable. The energy delivered to the wire has been empirically modelled but further investigation is required to understand the role of material heterogeneities and their effect on initiation. The energy delivered by the exploding wire has been quantified in many studies but it is not possible to correlate these to a particular design or firing regime and thus definitive conclusions are impossible. The energy absorbed by the PETN is also not well understood. Therefore, the initiation mechanism within the PETN EBWs has not been determined for EBW detonators. There is strong evidence that a shock-to-detonation (SDT) mechanism is not the sole cause of initiation. There is insufficient understanding of how electrical sparks transmit energy to PETN to make any judgement regarding the role gas ionization plays. Deflagration-to-detonation (DDT) seems like the most likely candidate for initiation but there is still a great lack of evidence for the presence of this mechanism.
Five shock-to-detonation transition (SDT) gas gun experiments for PBX 9012 were performed by Burns and Chiquete [2020]. The data from the experiments is available on the small-scale database. Here we focus on the tracker gauge data for the lead shock trajectory, x(t). In order to determine the shock speed U s (t) and the detonation transition point, a least squared fit to the shock trajectory is used based on a pair of ODEs with a fitting form for the shock acceleration developed by Hill and Gustavsen [2002].
It has previously been shown that the rate parameters of the SURF reactive burn model can be adjusted to fit the gap-stick data for PBX 9501 [Johnson et al., 2018, fig 4]. A recent recalibration of the SURF model for PBX 9501 (lot 730-010 at ρ = 1.837 g/cc 3 ) used different reactants and products EOS and an updated fitting form for the burn rate [Menikoff, 2021]. The new calibration fit gap-stick data [Hill et al., 2021] together with data from SDT experiments and curvature effect data. Here the results of gap-stick simulations using the xRage code are shown. They provide insight into the reactive flow in the gap-stick experiment.
Enhanced heat transfer model has been developed due to the turbulence in the cracks and pores in the burning fractured HE. This turbulence then serves in the initialization of the turbulence in the product gas. Enhanced heat transfer is also modeled in the turbulent product gas. Numerical simulations of the onset of HE burning and transition to detonation have been performed with and without the turbulent heat transfer model. Computational results indicate that onset of chemical reactions and transition to detonation is induced by the shockwave formed by the pressure rise caused by the heating. In other words, processes similar to shock-to-detonation transition (SDT) could become the driving mechanism for HE explosion and detonation in the case of heating rather than the temperature alone.
In this research, we assess the viability of four different, publicly available algorithms for estimating the azimuth and tilt parameters of solar photovoltaic systems using only the associated AC power time series data and site latitude-longitude coordinates. In this work, we curated a benchmarking data set of 44 fixed-tilt systems, comprising 275 measured AC power inverter data streams, with known azimuth and tilt parameters. Additionally, we isolated test cases in the data set with real-world issues, including shading and clipping, to determine how algorithm performance varies based on the presence of these phenomena. Using this data set for benchmarking, we evaluated the estimated vs. actual system characteristics for each algorithm, as well as the associated algorithm execution time using a standardized benchmarking process. The two highest performing algorithms were the Solar Data Tools and the PVWatts 5-based methods, which both achieved a median absolute error of approximately 5 and 1 degrees for azimuth and tilt, respectively. During run time analysis, the SDT method was approximately 5 times faster than the PVWatts 5-based method, with the median execution time for a stream varying between 6 and 8 seconds vs. a median run time of 31 seconds for the PVWatts 5-based method.
Here, this work presents the first application of isotopically labeled LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) to investigate thermal decomposition pathways. Specially synthesized LLM-105 isotopologues were utilized to isolate the influence of labeled 15 NO 2 nitro groups on the formation of lightgas products. Simultaneous differential scanning calorimetry, thermo-gravimetric, and mass spectrometry measurements were employed to track the evolution of product gases, enabling the direct comparison of isotopically shifted species with unlabeled LLM-105. Key findings show that C 2 N 2 production is mainly dependent on nitrogen sources from either the amine groups or the pyrazine ring (i.e., not the nitro groups). The formation of NO, N 2 , and N 2 O all involves the nitro groups to some extent. NO (nitric oxide) was found to be the predominant gas species directly formed from the nitro group of LLM-105. In contrast, mixed nitrogen isotopologues of N 2 and N 2 O (i.e., 14 N 15 N and 15 NNO) formed more readily in comparison to their pure counterparts (i.e., 15 N 2 and 15 N 2 O). This indicates the amine and/or pyrazine groups of LLM-105, in addition to the nitro group, are involved in the decomposition pathways forming N 2 and N 2 O. In addition, our investigation led to the discovery of two previously unreported decomposition products (CHO and HNCO), which were confirmed through hydrogen labelling utilizing deuterium isotopes. These results provide detailed speciation trends of gaseous products during LLM-105 decomposition, offering new insights into reaction pathways. Experimental data reported here will support the development of a detailed chemical kinetics model for LLM-105, essential for the safe handling of high explosives.
Graphite is a classic neutronic material that has been used as both a reactor reflector and moderator in various nuclear reactor systems. The ability to accurately predict the slowing down and thermalization of neutrons in graphite can have significant implications on the safety and operation of such reactor systems. In reactors, the neutron thermalization process is quantified using the thermal scattering law (TSL) and related cross sections for a given moderator. An ideal approach to assess the validity of TSL data is using benchmark measurement based on the pulsed Slowing-Down-Time technique and its comparison with the appropriate graphite nuclear library. In this work, experimental measurement and computational Monte Carlo simulations were performed to benchmark the slowing down characteristics and thermalization of neutrons in nuclear (reactor-grade) graphite. Given the density of graphite, various graphite libraries were selected for the benchmark analysis.