Engineering topics

Hobbs, Michael L.

Publications and source records attributed to Hobbs, Michael L..

5 records

DOE OSTI2024

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. In conclusion, we present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

cohesive powders↗

DOE OSTI2024

Modeling delayed thermal runaway in nitric acid-soaked cat litter mixed with radioactive waste

Thermal ignition of radioactive waste within a 55-gallon drum was simulated by using a pressure-dependent waste decomposition model (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) calibrated with data from full-scale drum experiments (Parker et al. in The thermolytic response of a surrogate RNS waste mixture at the drum scale. Los Alamos National Laboratory Report LA-UR-16-21760, 2016) and validated with experiments from multiple laboratories (Hobbs et al. in Thermal analysis of aged nitric acid-soaked kitty litter in TRU waste drums-23370.WM2023 Conference, Phoenix, AZ, 2023). The acceleration of nitric acid chemistry reacting with an organic cat litter leading to thermal ignition was likely triggered by a restricted vent in the drum. Here, we address whether the form of the rate equation in (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) is sufficient to extrapolate thermal ignition within aged drums of similar content that have been stored in Texas for over nine years by investigating four different reaction rate forms for waste decomposition. A critical reaction rate reduction analysis is performed on each of these models to determine if delayed thermal runaway within vented aged waste is possible after nine years. We found that a pressure-dependent first-order rate expression not only predicted the accidental ignition of the waste drum, but the form also matches multiple experiments from different laboratories. Even though the waste composition decreases over time, the model predicts that acceleration leading to thermal runaway is possible if the waste is confined, even after 9 years. In conclusion, waste containing oxidizers such as nitric acid should not be mixed with organic adsorbents, especially if the waste is confined.

cookoff↗

DOE OSTI2023

Predicting EBW detonator failure using DSC data

Exploding bridgewire detonators (EBWs) containing pentaerythritol tetranitrate (PETN) exposed to high temperatures may not function following discharge of the design electrical firing signal from a charged capacitor. Knowing functionality of these arbitrarily facing EBWs is crucial when making safety assessments of detonators in accidental fires. Orientation effects are only significant when the PETN is partially melted. Here, the melting temperature can be measured with a differential scanning calorimeter. Nonmelting EBWs will be fully functional provided the detonator never exceeds 406 K (133 °C) for at least 1 h. Conversely, EBWs will not be functional once the average input pellet temperature exceeds 414 K (141 °C) for a least 1 min which is long enough to cause the PETN input pellet to completely melt. Functionality of the EBWs at temperatures between 406 and 414 K will depend on orientation and can be predicted using a stratification model for downward facing detonators but is more complex for arbitrary orientations. A conservative rule of thumb would be to assume that the EBWs are fully functional unless the PETN input pellet has completely melted.

42 ENGINEERING↗

DOE OSTI2023

Swelling and permeability effects during propellant cookoff

Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. Furthermore, we found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS2000

Uncertainty Analysis of Decomposing Polyurethane Foam

Sensitivity/uncertainty analyses are necessary to determine where to allocate resources for improved predictions in support of our nation's nuclear safety mission. Yet, sensitivity/uncertainty analyses are not commonly performed on complex combustion models because the calculations are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, a variety of sensitivity/uncertainty analyses were used to determine the uncertainty associated with thermal decomposition of polyurethane foam exposed to high radiative flux boundary conditions. The polyurethane used in this study is a rigid closed-cell foam used as an encapsulant. Related polyurethane binders such as Estane are used in many energetic materials of interest to the JANNAF community. The complex, finite element foam decomposition model used in this study has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state decomposition front velocity calculated as the derivative of the decomposition front location versus time. An analytical mean value sensitivity/uncertainty (MV) analysis was used to determine the standard deviation by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation was essentially determined from a second derivative that was extremely sensitive to numerical noise. To minimize the numerical noise, 50-micrometer element dimensions and approximately 1-msec time steps were required to obtain stable uncertainty results. As an alternative method to determine the uncertainty and sensitivity in the decomposition front velocity, surrogate response surfaces were generated for use with a constrained Latin Hypercube Sampling (LHS) technique. Two surrogate response surfaces were investigated: 1) a linear surrogate response surface (LIN) and 2) a quadratic response surface (QUAD). The LHS techniques do not require derivatives of the response variable and are subsequently relatively insensitive to numerical noise. To compare the LIN and QUAD methods to the MV method, a direct LHS analysis (DLHS) was performed using the full grid and timestep resolved finite element model. The surrogate response models (LIN and QUAD) are shown to give acceptable values of the mean and standard deviation when compared to the fully converged DLHS model.

Hobbs, Michael L.↗