Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explosive flow”

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.

410 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Recurrence Rate and Magma Effusion Rate for the Latest Volcanism on Arsia Mons, Mars

Magmatism and volcanism have evolved the Martian lithosphere, surface, and climate throughout the history of Mars. Constraining the rates of magma generation and timing of volcanism on the surface clarifies the ways in which magma and volcanic activity have shaped these Martian systems. The ages of lava flows on other planets are often estimated using impact crater counts, assuming that the number and size-distribution of impact craters per unit area reflect the time the lava flow has been on the surface and exposed to potential impacts. Here we show that impact crater age model uncertainty is reduced by adding stratigraphic information observed at locations where neighboring lavas abut each other, and demonstrate the significance of this reduction in age uncertainty for understanding the history of a volcanic field comprising 29 vents in the 110-kilometer-diameter caldera of Arsia Mons, Mars. Each vent within this caldera produced lava flows several to tens of kilometers in length; these vents are likely among the youngest on Mars, since no impact craters in their lava flows are larger than 1 kilometer in diameter. First, we modeled the age of each vent with impact crater counts performed on their corresponding lava flows and found very large age uncertainties for the ages of individual vents, often spanning the estimated age for the entire volcanic field. The age model derived from impact crater counts alone is broad and unimodal, with estimated peak activity in the field around 130Ma (megaannum, 1 million years). Next we applied our volcano event age model (VEAM), which uses a directed graph of stratigraphic relationships and random sampling of the impact crater age determinations to create alternative age models. Monte Carlo simulation was used to create 10,000 possible vent age sets. The recurrence rate of volcanism is calculated for each possible age set, and these rates are combined to calculate the median recurrence rate of all simulations. Applying this approach to the 29 volcanic vents, volcanism likely began around 200-300Ma then first peaked around 150Ma, with an average production rate of 0.4 vents per Myr (million years). The recurrence rate estimated including stratigraphic data is distinctly bimodal, with a second, lower peak in activity around 100Ma. Volcanism then waned until the final vents were produced 10-90Ma. Based on this model, volume flux is also bimodal, reached a peak rate of 1-8 cubic kilometers per million years by 150Ma and remained above half this rate until about 90Ma, after which the volume flux diminished greatly. The onset of effusive volcanism from 200-150Ma might be due to a transition of volcanic style away from explosive volcanism that emplaced tephra on the western flank of Arsia Mons, while the waning of volcanism after the 150Ma peak might represent a larger-scale diminishing of volcanic activity at Arsia Mons related to the emplacement of flank apron lavas.

Richardson, Jacob A.↗
NASA NTRS1997

Cratering on Titan: A Pre-Cassini Perspective

The NASA-ESA Cassini mission, comprising a formidably instrumented orbiter and parachute-borne probe to be launched this October, promises to reveal a crater population on Titan that has been heretofore hidden by atmospheric haze. This population on the largest remaining unexplored surface in the solar system will be invaluable in comparative planetological studies, since it introduces evidence of the atmospheric effects of cratering on an icy satellite. Here, I highlight some impact features we may hope to find and could devote some modeling effort toward. Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s Atmosphere -94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface.Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s; Atmosphere about 94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface. Titan is comparable to Callisto and Ganymede for strength/gravity, Mars/Earth/Venus for atmospheric interaction, and Hyperion, Rhea, and Iapetus for impactor distribution. The leading/trailing asymmetry of crater density from heliocentric impactors is expected to be about 5-6, in the absence of resurfacing. Any Saturnocentric impactor population is likely to alter this. In particular the impact disruption of Hyperion is noted; because of the 3:4 orbital resonance with Titan, fragments from the proto-Hyperion breakup would have rapidly accreted onto Titan. Titan's resurfacing history is of course unknown. The disruption of impactors into fragments that individually create small craters is expected to occur. A crude estimate suggests a maximum separation of about 2 km (compared with 4 km on Venus, or 0.5 km on Earth). Crater chains are unlikely on Titan, since impactors must pass close enough to Saturn to be tidally disrupted; as a result, they would suffer aerodynamic disruption. Crater counting on adjacent satellites gives densities of about 200 per 10 (exp 6) square km for 20-km-diameter craters. However, the presence of a thick atmosphere leads to atmospheric shielding, depleting the relative abundance of small craters. This has been evaluated by models, and the relative abundance of small craters may be due to a diagnostic atmospheric collapse. A number of radar-dark "splotches" have been detected on Venus; these have been attributed to the interaction of the surface with the atmospheric shockwave produced by the Tunguska-like explosion of a bolide in the atmosphere. Simple analogy suggests that similar features might occur on Titan, but the shocked mass density (which controls the momentum coupling between the surface and the shockwave) of Titan's cold N2 atmosphere is about 20x smaller than that of Venus's hot CO2 atmosphere. Unless ice is much more easily turned to rubble than is rock, such features seem less probable on Titan. When the energy deposited by an impact forms a fireball with an equilibrate greater than one scale height, the fireball expands upward and can distribute ejecta. on ballistic exoatmospheric trajectories. On Venus this process is believed to be responsible for the parabolic features; the interaction of various-sized particles falling through the atmosphere with the zonal wind field winnows the particles to form a parabolic deposit. Although such a process is possible on Titan, the large scale height at higher altitudes would make it more difficult. Comparison with craters on other icy satellites suggests that craters on Titan will be fairly shallow (depth/diameter about 0.1) and craters greater than 10 km in diameter will have central peaks or domed bases, perhaps with central pits. The formation of ejecta. blankets may involve the atmosphere in a significant way, both by restraining the expansion of the ejecta cloud and by influencing the thermal history of the ejecta. Compared with Venus, Titan's atmosphere will chill an impact melt somewhat quickly, so the long ejecta flows seen on Venus seem less likely, detailed modeling needs to be performed to determine the impact melt production. Crater topography on Titan may be highlighted by the influence of liquids forming crater lakes. Craters with central peaks will typically form ring-shaped lakes, although horseshoe-shaped takes may be common; domed craters with central pits may even form bullseye lakes with islands with central ponds. If liquids have covered a substantial part of Titan's surface for a substantial period, hydroblemes and tsunami deposits may be common.

Lorenz, R. D.↗
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

LX-17 Thermal Decomposition-Characterization of Solid Residues from Cook-Off in a Small-Scale Vessel Under Confinement

Concerns surround whether insensitive (or any) energetic materials are more dangerous to handle when exposed to abnormal thermal environments. This study characterizes the residual material remaining after LX-17 (92.5 % 1,3,5-triamino 2,4,6-trinitro benzene (TATB) and 7.5 % Kel-F) is exposed to various thermal environments in a sealed small-scale vessel cook-off test reactor (heated at 0.1 to 100 °C/min until the reactor opened at 3000 psi (20.7 MPa)). Previous work has shown no additional sensitivity of these residues as evaluated by small-scale safety analysis, but characterization on the molecular scale indicates the TATB is transformed to more reactive compounds as well as the residue could be precursors to toxic gases. The solids and chars were characterized by various analytical methods. Heat-flow measurements indicated exothermic release is due to a mixture of residual TATB and related decomposition products (which may be more energetic). The N/C and O/N ratios indicated a material much more degraded than TATB. Primarily, the solids were a network of amorphous C inter-dispersed with N and O. Types of bonding include C–C, C–N, N–H, N–C, N=C, N≡C, C–O=, and –OH. Solvent extracts of the solids showed TATB decomposition intermediates benzo-furazans and benzo-furoxans, substituted TATB (mono-nitroso, hydroxyl, and chlorinated) along with several unidentified smaller molecules. These results indicate thermal treatment produces an amorphous carbon residue with heteroatoms incorporated through differing functionality, varying depending upon the thermal severity of exposure. These structures also could further decompose producing toxic light gases (such as cyanide).

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

Independent Review of Groundwater Remediation Strategy for Hexavalent Chromium and RDX Groundwater Plumes at Los Alamos National Laboratory (Rev. 1)

Site operations at the Los Alamos National Laboratory (LANL) resulted in the release of oxidized chromium, Cr(VI), into Sandia Canyon from cooling tower effluent from 1956 until 1972. The chromium traveled with the surface water approximately 3 miles downstream before migrating below ground surface. Chromium concentrations exceed 50 μg/L in the upper portion of the aquifer. Another LANL groundwater plume of concern is associated with RDX (Royal Demolition Explosives, 1,3,5-trinitro-1,3,5-triazine). Between 1951 and 1996, RDX was released to the mesa-top facilities' process water outfall, adjacent and underlying soils, and alluvial sediments, along with surface water in Cañon de Valle. Between 2000 and 2010, two remedial actions were deployed, removing much of the near-surface RDX, however, recharge due to precipitation has transported RDX into the perched-intermediate zone and into the regional aquifer. The report documents an independent technical review by scientists from the Department of Energy (DOE) Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) to provide recommendations for potential near term actions to address and optimize remediation for both the Cr(VI) and RDX plumes. The proposed near-term remedial actions include design of pump and treat systems for Cr(VI) and monitoring and study for natural attenuation for RDX. The review assesses existing data, conceptual and numerical modeling, and it recommends a technical integration process to support identifying and implementing strategic, effective and efficient remedies. The DOE Environmental Management Los Alamos Field Office (EM-LA) and their cleanup contractor Newport News Nuclear-BWTX, LLC Los Alamos (N3B) provided the information required for the review. Interviews were also conducted with regulators to obtain the full spectrum of technical, regulatory and scientific perspectives. The independent review team was impressed by the capabilities, experiences, innovativeness, and insightfulness of the technical representatives from both the regulator, the New Mexico Environment Department (NMED) and N3B. Incorporation of vadose zone flow pathways in the conceptual site model (CSM) and configuring the numerical modeling for the site was generally state-of-the-practice (or better). This could be considered state-of-the-art by addressing uncertainties related to spatial extent of hydraulic windows. The reviews from the regulators were thorough and often provided useful concepts for consideration and future/study resolutions. The overarching consensus recommendation of independent review team is that the LANL groundwater plumes should be addressed in context of the emerging "management of complex sites" paradigm.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2003

Fragmentation under the Scaling Symmetry and Turbulent Cascade with Intermittency

Fragmentation plays an important role in a variety of physical, chemical, and geological processes. Examples include atomization in sprays, crushing of rocks, explosion and impact of solids, polymer degradation, etc. Although each individual action of fragmentation is a complex process, the number of these elementary actions is large. It is natural to abstract a simple 'effective' scenario of fragmentation and to represent its essential features. One of the models is the fragmentation under the scaling symmetry: each breakup action reduces the typical length of fragments, r (right arrow) alpha r, by an independent random multiplier alpha (0 < alpha < 1), which is governed by the fragmentation intensity spectrum q(alpha), integral(sup 1)(sub 0) q(alpha)d alpha = 1. This scenario has been proposed by Kolmogorov (1941), when he considered the breakup of solid carbon particle. Describing the breakup as a random discrete process, Kolmogorov stated that at latest times, such a process leads to the log-normal distribution. In Gorokhovski & Saveliev, the fragmentation under the scaling symmetry has been reviewed as a continuous evolution process with new features established. The objective of this paper is twofold. First, the paper synthesizes and completes theoretical part of Gorokhovski & Saveliev. Second, the paper shows a new application of the fragmentation theory under the scale invariance. This application concerns the turbulent cascade with intermittency. We formulate here a model describing the evolution of the velocity increment distribution along the progressively decreasing length scale. The model shows that when the turbulent length scale gets smaller, the velocity increment distribution has central growing peak and develops stretched tails. The intermittency in turbulence is manifested in the same way: large fluctuations of velocity provoke highest strain in narrow (dissipative) regions of flow.

Gorokhovski, M.↗
NASA NTRS1998

The Coronae of Moderate-Mass Giants in the Hertzsprung Gap and the Clump

We have used the Roentgensatellit (ROSAT), the Extreme Ultraviolet Explorer (EUVE), and the Hubble Space Telescope (HST) to measure X-ray and ultraviolet emissions of moderate-mass (APPROX. 2-3 solar mass) giants in the Hertzsprung gap (spectral types early F to mid-G) and the post-helium flash "clump" (approx. G8-K0). Our motivation was to document the evolution of hot coronae (T greater than 10(exp 6)K) along the post-main-sequence trajectories traveled by such stars in order to gain insight concerning the "X-ray deficiency" of the F-GO giants and the strong braking of stellar rotation at the red edge of the Hertzsprung gap. With few exceptions, Hertzsprung gap and clump giants observed by ROSAT PSPC show hot (T approx. 10(exp 7)K) coronal energy distributions, regardless of any X-ray deficiency, EUVE spectra of gap star 31 Com (G0 111) indicate a broad coronal emission measure hump at approx. 10(exp 7.2)K, while the active clump giant beta Ceti (K0 III) displays a sharp peak at approx. 10(exp 6.8)K, as seen previously in the mixed clump/gap binary Capella (alpha Aur: G8 III + G0 III). The gap giants upsilon Peg (F8 III) and 24 UMa (G4 III) have EUV emissions of intermediate temperature (approx. 10(exp 7.0)K). The stars 31 Com, psi(sup 3) Psc (G0 III), and beta Cet exhibit redshifted transition zone (TZ: approx. 10(exp 5)K) lines in HST GHRS spectra, as reported earlier in Procyon (alpha CMi: F5 IV-V) and Capella G0. Such redshifts on the Sun are thought to signify flows in magnetic loops. beta Cas (F2 III)-a rare soft coronal source among the gap stars-displays blueshifts of C iv and 0 iv], although emissions at cooler and hotter temperatures are near the photospheric velocity. The remarkably broad line profiles of the fastest rotating gap giants suggest that the 10(exp 5) K "subcoronal " emission zones extend to h approx. R(sub *) above the photosphere, about 50 scale heights. In contrast to the TZ line redshifts, the upper chromospheric emissions (e.g., Mg II and Si III) of 31 Com and upsilon(sup 3) Psc have blueshifted cores. Blue-asymmetric peaks in the solar Mg a lines are thought to indicate dynamical heating in the chromosphere. Observations of the H(sub I) Ly(alpha) feature of 31 Com taken 9 months apart reveal striking profile changes, reminiscent of those noted previously in the Ly(alpha) blue peak of the Capella G0 star. We used the far-ultraviolet diagnostics, in combination with ROSAT X-ray photometry and EUVE high-excitation line strengths, to constrain physical models of the stellar outer atmospheres. Quasi-static magnetic loops can simulate the empirical coronal emission measures of the giant stars, but the inferred pressures for sensible loop lengths conflict with direct measurements of subcoronal densities. Furthermore, the high rate of emission at approx. 10(exp 5) K cannot be explained by thermal conduction down the legs of hot quasi-static loops. On the other hand, the possible existence of elongated (l - R(sub *) emission structures on the gap giants leads to a speculative scenario to explain the X-ray deficiency. It is based on the increased importance of the dynamical filling phase ("explosive evaporation") of the loop life cycle; conductive cooling, yielding TZ emissions at the footpoints, when the heating is interrupted; and the possibility for transitions between " hot " and " cool " energy balance solutions owing to dynamical suspension and centrifugal trapping of the cooling gas. The long loops might represent a vestigial global " magnetosphere " inherited from the main-sequence phase, which ultimately is disrupted near approx. G0 by the deepening convective envelope and growth of a more solar-like dynamo. Coronal emissions might be boosted temporarily as the X-ray deficiency is removed but soon would be quenched by wind braking previously inhibited by the magnetospheric "dead zone."

Ayres, Thomas R.↗
DOE OSTI2020

Investigating the minimum post-burst energy required to function an exploding bridgewire detonator

The minimum energy required after bridge-burst to function a commercially available RP-1 exploding bridgewire detonator with a gold bridgewire was investigated using two capacitor discharge units (CDUs) making use of state-of-the-art on/off high-voltage high-current switches. Both switches are capable of shutting off the current flow in ≈100ns and the conduction time was varied systematically to limit the duration of electrical energy delivered from the CDU until a detonation go/no-go threshold was reached. It was found that the minimal extra electrical energy required for a soft-fire (a statistically marginal go) was approximately half that required for a hard-fire (a definite go). Furthermore, the measured energy required to function the RP-1 detonator from the 11.5 kV CDU was approximately 20 mJ greater than for the 600 V CDU, and it is postulated that this was due to arcing across the detonator terminals and an associated current shunt that was not efficient at depositing energy conducive to detonator functioning. Measurements of the bridgewire temperatures bursting in air showed that the maximum temperature reached was consistent (≈18kK) when sufficient energy was delivered that the detonator functioned (either soft- or hard-fire), but that the temperature for a sub-threshold bridgewire burst was significantly cooler at approximately 8 kK. Furthermore, it was demonstrated that the intensity of the light emitted at burst scaled with deposited energy. A related finding is that the light emission from the bridgewire peaks some tens of nanoseconds after both peak voltage and peak electrical power.

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

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Hydrogen Leak Modeling for Development of Smart Distributed Monitoring Under Unintended Releases

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower. Hydrogen has the potential to play a crucial role in decarbonizing industrial processes that are currently reliant on fossil fuels and provide long-duration and/or seasonal energy storage to enable electricity decarbonization. Hydrogen can also be used as a fuel for fuel cell vehicles, providing a zero-emission alternative to traditional internal combustion engines. DOE launched the Hydrogen Energy Earthshot (Hydrogen Shot) in June 2021 to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). While promising, Hydrogen is highly-flammable, and in the presence of oxygen, it can form explosive mixtures. . Therefore, understanding leak scenarios is essential to evaluate and mitigate the safety risks associated with potential hydrogen leaks. An increased understanding of leak behavior, and having tools to model leaks, can help assess how hydrogen would disperse in different environments, influencing emergency response plans and safety measures, and identify potential issues with materials and design systems that can withstand the challenges posed by hydrogen. Recently, researchers have attempted to study hydrogen leaks for development of risk management strategies. However, the focus has been on closed or semi-closed spaces like storage rooms, vehicles, garages, and fueling stations - all promising locations for future hydrogen infrastructure. In this presentation, the modeling environment extends the span of research further by modeling hydrogen leak in an outdoor, open space. We will present the key challenges with modeling hydrogen leaks in an uncontrollable environment, how they were handled, and how modeling results informed sensor selection and placement. A Hydrogen research facility at the National Renewable Energy Laboratory (NREL) was used as a case study to model hydrogen leaks. In the future, Hydrogen wide area detection methodologies will be developed and tested at this site to monitor for unintended and operational hydrogen releases. The data generated from modeling will be used to develop a predictive model to detect hydrogen leak location based on concentration measured by sensors in this open space. Furthermore, the facility was also chosen because controlled hydrogen releases can be performed. A computational fluid dynamics (CFD) based modeling approach was taken to model hydrogen leak. The full-scale hydrogen facility was modeled with a large ambient domain. The electrolyzer at the facility can produce a controlled release rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind direction, wind speed at various altitudes, and temperature were used as inputs to the model. To capture the variability of weather conditions, a subset of the weather conditions experienced during daytime hours without precipitation over the course of three months was generated; using established data clustering techniques, a total of 100 condition sets were chosen. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. These distributions are compared to experimental data from a constant mass flow, controlled hydrogen release at the facility. The stochastic wind conditions of the release make direct validation difficult, therefore, statistical comparison approaches were used. Wind conditions are found to significantly impact the release behavior, including direction and concentration. Sensor selection and placement is proposed for the facility and is now based on release behavior predicted for the facility given its weather patterns; this is much more informed than without the modeling results. The methodology and analysis procedure can be translated to other facilities using modified geometries and site-specific weather conditions. Hydrogen holds great promise as a renewable energy fuel, but ensuring safety in its production, storage, and use is paramount. Studying potential leak scenarios in an open space will help develop sensors to detect hydrogen on a large spectrum of concentration and eventually build a smart distributed monitoring system.

CFD↗
DOE OSTI2023

Review of particle deposition to and removal from clothing, skin, and hair after a radioactive airborne dispersal event

Explosive Radiological Dispersal Devices (RDD) – aka dirty bombs – are seen as a credible method to carry out a radiological terror attack. After exploding a radioactive source, the radionuclide-laden plume will be blown downwind of ground zero, with particles falling out and potentially depositing on people caught in and under the cloud. Some of these people may not show any sign of radiation sickness and therefore not realize they have been contaminated and may take the radioactive particulate with them on their daily activities, thus spreading the radioactive particulate outside the initially contaminated area. This paper reviews the scientific literature to better understand the rate at which particulate deposits on and is removed from the different “surfaces” of a person, i.e., hair, skin, and clothing. Prior research indicates that: 1) particle deposition is usually higher on skin than on hair and clothing; 2) particle deposition is greater for a person with higher skin moisture, 3) stronger wind increases the deposition flux onto a person, and 4) the fraction of particulate deposited on the hair, skin, and clothing respectively depends on the length of the hair, assuming all the hair surface is available for deposition. The studies taken into consideration show that the largest uncertainty in particulate deposition onto a person is due to clothing type because of the different possible weave arrangements and tightness which translate into differences in actual surface area and surface roughness. A factor of 2-to-20 variation in deposition rate was found. Removal of the particulate from the contaminated person may be due to wind, a person's movement, and/or contact transfer, i.e., by touching a different clean surface. Experiments show that the majority of the particulate is resuspended within 2–6 h mostly depending on the intensity of physical activity. The largest uncertainty in particulate removal from skin depends on the skin moisture, transfer rate of single-contact, and how many objects/people a person touches per hour. No data for hair were found for particle removal and resuspension. The studies considered did not utilize radionuclides directly; however, data on adhesion of radioactive vs. their non-radioactive counterpart have shown that the uncertainty due to the radioactivity of the particles is lower than that due to other factors. In conclusion, an idealized scenario involving a single building in the path of the cloud showed the impact of building-influenced flow on the cloud transport path and mixing, which affects the radiological dose the downwind population is exposed to and consequently the health effects.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

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

Effect of Sample Mass, Confinement, and Preheating Time on the Thermal Response of LLM‐105: Experiments and Kinetic Analysis

Various small-scale experiments were performed to provide data for developing a model to predict the thermal response of LLM-105 over a wide range of conditions. The thermal decomposition of LLM-105 was studied as a function of sample mass, confinement of volatile products, and preheating time in both isothermal and ramped heating experiments. The thermal decomposition of LLM-105 is a two-step process, as shown by the two exothermic peaks in the heat flow profiles, which were fitted to two nth-order autocatalytic reaction models with a similar activation energy of ∼289 kJ/mol. The magnitude and shape of these peaks varied with sample mass and confinement. Increasing sample mass enhanced the second exotherm with respect to the first one, while increasing the level of confinement promoted a transition from a sublimation-dominated regime towards thermal decomposition. The effect of LLM-105 particle size on the rate of weight loss was evident for open-pan experiments, where bigger particles sublimed at lower temperatures than smaller particles. Thermal response and solid residue composition of LLM-105 samples were analyzed following preheating for different durations. Longer preheating times caused a shift of the second exotherm to lower temperatures and a decrease in the reaction enthalpy, confirming that LLM-105 decay is a consecutive reaction mechanism, probably autocatalytic. In conclusion, the kinetic model derived from ramped experiments was validated against the measured LLM-105 fraction remaining and the enthalpy remaining of the solid residue as a function of preheating times and showed good agreement.

Chemistry - Chemical explosives↗
DOE OSTI2025

Vulcan Test Platform: Demonstrating the Data Center as a Flexible Grid Asset

Explosive data center demand is outpacing grid infrastructure development. AI workloads and hyperscale cloud growth are creating unprecedented power requirements, while traditional grid expansion faces multi-year development timelines, regulatory hurdles, and decarbonization challenges. Sidewalk Infrastructure Partners recognized this impending crisis years ago and founded Verrus to develop an innovative solution: data centers that function as grid assets rather than passive loads. The Verrus approach integrates proprietary grid-aware controls with battery energy storage systems (BESS) in a medium-voltage architecture that delivers three critical capabilities: Fast-responding demand flexibility that can service requests from the utility within 10 seconds, Uninterrupted transition to islanded operation during grid disturbances, Continuous uptime assurance while maintaining all customer service level commitments Through Verrus' strategic partnership with the National Renewable Energy Laboratory (NREL), we validated these capabilities on Vulcan, a 70 MW utility-scale test platform powered by NREL's ARIES Virtual Emulation Environment. This deployment-ready technology has successfully demonstrated that Verrus data centers can deliver meaningful grid services while maintaining mission-critical reliability. This technical report outlines the design, methodology, and results of this emulated deployment, demonstrating that data centers can play a pivotal role in enhancing grid flexibility and reliability, without sacrificing service level guarantees.

24 POWER TRANSMISSION AND DISTRIBUTION↗