Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPLOSIVE FORMING”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Non-Contact Vapor Detection of Illicit Drugs via Atmospheric Flow Tube-Mass Spectrometry

Real-time, non-contact detection of illicit drugs is a desirable goal for the interdiction of these controlled substances, but the relatively low vapor pressures of such species present a challenge for trace vapor detection technologies. The introduction of atmospheric flow tube-mass spectrometry (AFT-MS), which has previously been demonstrated to detect gas-phase analytes at low parts-per-quadrillion levels for explosives and organophosphorus compounds, also enables the potential for non-contact drug detection. With AFT-MS, direct vapor detection of cocaine and methamphetamine from ~5 µg residues at room temperature is demonstrated herein. Furthermore, thermal desorption of low- to sub-picogram levels of cocaine, methamphetamine, fentanyl, and heroin is observed via AFT-MS using a carrier flow rate of several L/min of air. These low levels can permit non-contact sampling through collection of vapor, effectively preconcentrating the analyte before desorption and analysis. Quantitative evaluation of the thermal desorption approach has yielded limits of detection (LODs) on order of 10 fg for cocaine and fentanyl, 100 fg for methamphetamine, and 1.6 pg for heroin. The LOD for heroin was lowered to 200 fg by using tributyl phosphate as a dopant to form a proton-bound heterodimer with heroin. When used with AFT-MS, the intentional formation of specific drug-dopant adducts has the potential to enhance detection limits and selectivity of additional drug species. Species that are prone to form adducts present a challenge to analysis, but that difficulty can be overcome by the intentional addition of a dopant. Molecules unlikely to form adducts will remain essentially unimpacted, but the adduct-forming species will interact with the dopant to compress the analyte signal into a single peak. This approach would be valuable in the application of non-contact screening for illicit substances via vapor collection followed by thermal desorption for analysis.

methamphetamine, cocaine, heroin, fentanyl, atmosp↗
DOE OSTI2019

Cookoff experiments of a melt cast explosive (Comp-B3)

Validated models of melt cast explosives exposed to accidental fires are essential for safety analysis. In the current work, we provide several experiments that can be used to develop and validate cookoff models of melt cast explosives such as Comp-B3 composed of 60:40 wt% RDX:TNT. We present several vented and sealed experiments from 2.5 mg to 4.2 kg of Comp-B3 in several configurations. We measured pressure, spatial temperature, and ignition time. Some experiments included borescope images obtained during both vented and sealed decomposition. We observed the TNT melt, the suspension of RDX particles in the melt, bubble formation caused by RDX decomposition, and bubble-induced mixing of the suspension. The RDX suspension did not completely dissolve, even as temperatures approached ignition. Our results contrast with published measurements of RDX solubility in hot TNT that suggest RDX would be completely dissolved at these high temperatures. These different observations are attributed to sample purity. We did not observe significant movement of the two-phase mixture until decomposition gases formed bubbles. Bubble generation was inhibited in our sealed experiments and suppressed mixing.

42 ENGINEERING↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mass measurement of 27 P to constrain type-I x-ray burst models and validate the isobaric multiplet mass equation for the $A = 27, T = \frac{3}{2}$ isospin quartet

Background: Light curves are the primary observable of type-I x-ray bursts. Computational x-ray burst models must match simulations to observed light curves. Most of the error in simulated curves comes from uncertainties in rp process reaction rates, which can be reduced via precision mass measurements of neutron-deficient isotopes in the rp process path. Purpose: Perform a precise atomic mass measurement of 27 P. Use this new measurement to calculate rp process reaction rates and input these rates into an x-ray burst model to reduce simulated light curve uncertainty. Use the mass measurement of 27 P to validate the isobaric multiplet mass equation (IMME) for the A = 27 T = 3/2 isospin quartet which 27 P belongs to. Method: High-precision Penning trap mass spectrometry utilizing the time-of-flight ion cyclotron resonance technique was used to determine the atomic mass of 27 P. The MESA code (Modules for Experiments in Stellar Astrophysics) was then used to simulate x-ray bursts using a one-dimensional multizone model to produce updated light curves. Results: The mass excess of 27 P was measured to be –670.7(6) keV, a 14-fold precision increase over the mass reported in the 2020 Atomic Mass Evaluation (AME2020). The 26 Si(p, $γ$) 27 P– 27 P($γ$, p) 26 Si rate equilibrium has been determined to a higher precision based on the precision mass measurement of 27 P. x-ray burst light curves were produced with the MESA code using the new reaction rates. Changes in the mass of 27 P seem to have minimal effect on light curves, even in burster systems tailored to maximize impact. Conclusion: The mass of 27 P does not play a significant role in x-ray burst light curves. It is important to understand that more advanced models do not just provide more precise results, but often qualitatively different ones. This result brings us a step closer to being able to extract stellar parameters from individual x-ray burst observations. In addition, the IMME has been validated for the A = 27, T = 3/2 quartet. The normal quadratic form of the IMME using the latest data yields a reduced $χ$ 2 of 2.9. Finally, the cubic term required to generate an exact fit to the latest data matches theoretical attempts to predict this term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Amorphization and Crystallization of Hexanitroazobenzene (HNAB) Using Conventional DSC and Flash DSC

Abstract The present work presents results from an investigation of the glass transition and crystallization behaviors of HNAB tested over more than five orders of magnitude of cooling rate from 0.005 °C/s to 600 °C/s (0.3 to 36000 °C/min) by a combination of conventional and Flash differential scanning calorimetry (DSC). The work quantifies the influence of the thermal amorphization route on the properties of this high explosive. Cooling rates faster than 100 °C/s (6000 °C/min) result in amorphous HNAB as expected from prior work, but we also find that amorphization of the HNAB occurs at cooling rates slower than 0.008 °C/s (0.5 °C/min). The behavior of the amorphous HNAB made by slow cooling is compared with that of amorphous HNAB made by fast cooling, as well as with that made by solvent casting in terms of glass transition temperature, apparent activation energy of glass transition, and dynamic fragility parameter m . Besides, the non‐isothermal crystallization response as a function of cooling rate is also reported. The thermal stability and decomposition energy of amorphous HNAB are compared with those of the crystalline counterpart, being similar heats of decomposition of 3295 and 3392 J/g, respectively; suggesting that the amorphous HNAB will have similar thermal stability and chemical energy to the crystalline form.

Koh, Yung P.↗
DOE OSTI2022

Observing and modeling the effects of production infrastructure in electromagnetic surveys

Electromagnetic (EM) methods are among the original techniques for subsurface characterization in exploration geophysics because of their particular sensitivity to the earth electrical conductivity, a physical property of rocks distinct yet complementary to density, magnetization, and strength. However, this unique ability also makes them sensitive to metallic artifacts — infrastructure such as pipes, cables, and other forms of cultural clutter — the EM footprint of which often far exceeds their diminutive stature when compared to that of bulk rock itself. In the hunt for buried treasure or unexploded ordnance, this is an advantage; in the long-term monitoring of mature oil fields after decades of production, it is quite troublesome indeed. In this study we consider the latter through the lens of an evolving energy industry landscape in which the traditional methods of EM characterization for the exploration geophysicist are applied toward emergent problems in well-casing integrity, carbon capture and storage, and overall situational awareness in the oil field. We introduce case studies from these exemplars, showing how signals from metallic artifacts can dominate those from the target itself and impose significant burdens on the requisite simulation complexity. Furthermore, we show how recent advances in numerical methods mitigate the computational explosivity of infrastructure modeling, providing feasible and real-time analysis tools for the desktop geophysicist. Lastly, we demonstrate through comparison of field data and simulation results that incorporation of infrastructure into the analysis of such geophysical data is, in a growing number of cases, a requisite but now manageable step.

58 GEOSCIENCES↗
DOE OSTI2023

Analysis of degradation products in thermally treated TATB

Delineating the chemical composition of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) residues produced from the exposure to abnormal thermal environments should lead to a better understanding of the decomposition paths. Identifying and quantifying each compound in thermally produced residues, monitors which compounds are degrading or forming along the decomposition route, as well as providing input for the kinetic models of those pathways. Here, in this paper, we report the methodology of isolating, identifying, and where possible, quantifying soluble compounds present in solid residues of thermally treated TATB (330 °C for tens of minutes). Samples were extracted with DMSO, separated using chromatography, and quantified using their absorption at 354 nm. Identification of unknown compounds was accomplished using high resolution mass spectrometry. TATB, F1 (diamino-dinitro-benzofurazan), HO-TATB (2,4,6-triamino-1-hydroxyl-3,5-dinitrobenzene), and T4A (1-chloro-3,5-dinitro-2,4,6-triaminobenzene) were trace compounds detected in the unreacted TATB. Ten more compounds that formed in the residues were structurally identified including F2 (amino-nitro-difurazan). Several more compounds were observed but not completely identified. We propose possible structures for the unknowns. Of the compounds formed, F1 was the most abundant compound reaching 4.5 % by weight of the degraded solid sample. Other degradation compounds were estimated to sum to trace levels, well below 1 %. Most compounds were new, having not been detected and identified in previous studies of production grade and thermally aged TATB. Many compounds only reached detectable concentrations after several min of thermal exposure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strontium-84 Enrichments in Presolar Grains Provide First Evidence of p -process Nucleosynthesis in Core-collapse Supernovae

This study reports detection of rare p-process isotopes within presolar grains. Presolar grains are relic dust grains from dying stars. These microscopic dust particles are found in primitive solar system materials. Their distinct isotopic compositions record the nucleosynthetic processes in their parent stars and the Galactic chemical environment in which these stars formed. We studied presolar graphite grains of high-density type from the Murchison meteorite and found five grains with subgrains that show enrichments in 84 Sr compared to the solar abundance. 84 Sr is the neutron-deficient isotope of strontium that can be produced in the deep oxygen-rich interior of high-mass stars that end their lives as core-collapse supernovae. The observed 84 Sr excesses cannot be produced in low-mass asymptotic giant branch stars, the source of most high-density presolar graphites found in meteorites. High-density graphites with embedded 84 Sr excesses are, instead, compatible with a core-collapse supernovae origin. The graphite subgrains condensed from carbon-rich materials in the outer layers of core-collapse supernovae, where 84 Sr was destroyed by neutron captures during hydrostatic evolution of the stars and their final explosion. Based on current theoretical stellar models, a few percent of contribution from the inner regions of core-collapse supernovae, which are enriched in p-process nuclides, to the outer carbon-rich regions is the most likely explanation for the observed enrichment of 84 Sr in the subgrains of the high-density graphites. In this study, we present the first observational evidence that core-collapse supernovae produce and eject isotopes made by the p-process.

Astronomy and AstroPhysics↗
DOE OSTI2021

Development of electronically passivating surfaces to enhance battery performance

For many decades, lithium metal batteries (LMBs) have been studied for their high energy density, though they suffer from high reactivity, which safety concerns such as explosions and fires. These events can occur as a result of degradation due to redox reactions and lithium morphology that evolve during cycling. Mostly, the failure is caused by the electrolyte degradation leading to the creation of a solid electrolyte interphase (SEI), which tends to induce electron leakage, lithium dendrite growth, and capacity decay. SEI formation is practically inevitable. This work investigated the minimization of degradation events by modifying different battery configurations. Modifications were made by exchanging electrolytes, electrode surfaces, and even annexing innovative components. Our goal is to achieve a configuration where Li-ion transfer is facilitated, and electron transfer is limited which would enable a longer lithium metal battery life cycle. When performing these modifications and testing their electrochemical performance, a current reduction was observed when the electrodes were spontaneously coated with preformed SEI. Herein several classes of artificial SEI layers are explored for their ability to form passivating layers through use of a redox active molecule.

25 ENERGY STORAGE↗
DOE OSTI2021

Seeding the explosion of a high-current-density conductor in a controlled manner through the addition of micron-scale surface defects

Inhomogeneities in a current-carrying conductor promote non-uniform heating and expansion through the complex feedback between current density, electrical resistivity, Ohmic heating, temperature, and hydrodynamics. Three-dimensional-magnetohydrodynamic (3D-MHD) simulations suggest that μm-scale resistive inclusions or voids seed local overheating and through hydrodynamic explosion generate continuously growing craters which become several times larger than the initial perturbation. The ejected mass is the genesis of an electrothermally driven plasma filament which develops at lower current than plasmas on uniform surfaces adjacent to the defect. This result suggests that 1D or even 2D treatments are largely inadequate for detailed prediction of plasma formation. To test computational predictions, z-pinch experiments driven to 1 MA studied ultra-high-purity aluminum rods which were then machined to include pairs of quasi-hemispherical voids or “engineered defects (ED)” on the current-carrying surface. ED are the dominant current-density perturbation and reproducibly drive local overheating which can be compared with 3D-MHD simulation. Data from high-resolution-gated imagers of visible surface emissions confirm many simulation predictions, including the surface topography of local overheating, and the propensity for neighboring ED to prematurely source plasmas which then connect to form a plasma filament. Results also provide conditional support of theory which suggests heating similarity; that is, heating is independent of ED size for geometrically scaled ED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An Advanced Pretreatment/Anaerobic Digestion (APAD) Technology for Increased Conversion of Sewage Sludge to Bio-natural Gas in Small-scale Wastewater Plants of less than Five tons Sewage Sludge a Day

The problem today of energy production from sewage sludge at small-scale is that conventional Anaerobic Digestion (AD) as used today at Wastewater Treatment Facilities (WWTF) produces too little energy for warrant use of the biogas. It further leaves 50% or more of the waste behind after the treatment. To overcome this problem, we proposed a novel concept based on Advanced Wet Oxidation & Steam Explosion (AWOEx) of the recalcitrant parts of sewage sludge left behind after AD. We further suggest upgrading biogas to renewable natural gas (RNG) using gaseous fermentation of biogas with hydrogen added by a new methanogen. Overall, the DOE funded Advanced Pretreatment & Anaerobic digestion (APAD) project showed significant improvements over current practice. The project demonstrated that AWOEx followed by AD significantly enhanced the carbon conversion efficiency from 37% to 62%, an increase of 68%. This is far higher than the metric for the specific FOA of an increase of 50%. Besides, the project showed high efficiency of our biological conversion of biogas into RNG when using a new isolate of Methanothermobacter wolfeii resulting in a 100% increased production of a refined biogas with maximum 5% CO2. With both AWOEx pretreatment and biogas upgrading, the project showed a CCE of ca. 83%, far higher than any previous work on sewage sludge. Besides over 200% higher amount of energy in the form of RNG, the APAD concept will reduce disposal cost due to significant reduction in the concentration of final sludge product after APAD. The APAD technology can operate as a bolt-on to a conventional AD plant for improving conversion of the residual organics after AD as done in this DOE project. It can further be implemented as a stand-alone process with AWOEx followed by AD for WWTF’s currently operating without AD.

09 BIOMASS FUELS↗
DOE OSTI2025

ALARA Review for Venting of Flanged Tritium Waste Containers (FTWCs) at TA-54

This project will vent headspace hydrogen and oxygen from specialized high-pressure storage vessels called flanged tritium waste containers (FTWCs). There are four of these containers located in Los Alamos National Laboratory’s (LANL’s) Technical Area 54 (TA-54), Building 1028. This building is in the southwest corner of LANL’s Material Disposal Area G. The vented headspace gas is expected to contain tritium in the form of water vapor, elemental hydrogen or a combination of both. The venting operation is expected to take place in 2025. There will be a series of operational readiness reviews prior to venting activities commencing. The FTWCs at TA-54 contain tritium-contaminated metal parts and molecular sieve media, which is a pebble-like material used to absorb water vapor from the air. This molecular sieve media inside the FTWCs is contained in metal canisters, along with some loose media material in bags. Over time, tritiated water vapor that had been adsorbed onto the media can become liberated into the FTWC headspace gas. Radiolysis can cause separation of the water vapor into its hydrogen and oxygen components, resulting in the potentially hazardous gas mixture within the FTWC. LANL has determined that continued tritium storage in these containers can pose an unsafe condition due to possible hydrogen and oxygen gas buildup within the FTWCs with a potential for explosion if sparks are generated within the FTWCs. To mitigate this hazard, the FTWCs will be stabilized by venting them in-place to remove hazardous gases.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

High Solids Loading Aqueous Slurry Formation of Corn Stover Before Pretreatment in a Fed-batch Bioreactor

With the increase on population the world will depend on renewable sources to meet the increasingly energy needs. Use of lignocellulosic biomass as a renewable source has been proven efficient for conversion to cellulosic ethanol and capable of contributing to achieve the threshold on energy demand while reducing greenhouse gases in 90% when compared with fossil fuels (Wang et al., 2007, 2012). However, the processing of biomass encounters limitations in feeding and flow within biorefineries due to the system plugging by compaction and slurry high yield stress, preventing transport of biomass materials and in some cases unexpected plant shutdowns that results in high operational costs (dos Santos et al., 2021; Ximenes et al., 2021). Different solutions for biomass handling and slurry formation from densified materials like pellets have been studied. Cellulosic biomass residues are typically processed in pretreatment reactors to which acid or base is added (Humbird et al., 2002). Other pretreatments such as liquid hot water (LHW) and steam explosion that act without the addition of chemicals use pressures above the saturation vapor of water, to disrupt biomass structure (Ruiz et al., 2021). Although these methods are widely known, the use of pretreatments like acid digestion and LHW present challenges like waste disposal and could be energy inefficient (Mosier et al., 2005). Alternatively other approaches like enzymatic liquefaction have been applied for biomass transformation into slurries with promising results and without the use of pretreatments (dos Santos et al., 2021). Nevertheless, enzymatic liquefaction implementation faces difficulties due to the recalcitrant properties of biomass, its variability, and the release of enzyme inhibitors. While it has been proven that chemical composition properties in biomass have effects on enzyme inhibition (Huang et al., 2022; Kim et al., 2011; Zhai et al., 2018), hindering the efficiency of the liquefaction, not ample research has been conducted in understanding the physical properties as particle size, porosity, water adsorption retention and their effect in the ability to form a slurry at high solids concentration (Yan et al., 2020). In order to improve lignocellulosic biomass handling and formation of biomass slurries, enzyme assisted liquefaction for slurry creation from corn stover at solids loadings up to 30% is reported in this work. Two different kinds of biomass (pelleted corn stover and cobs) were liquefied in a fed-batch process using commercial enzymes Celluclast 1.5L or Ctec-2 at 1FPU or 3 FPU per gram of dry solids in 10 mM sodium citrate buffer solution (pH 4.8). Pellets were fed into a 1 L stirred bioreactor according to a predefined fed-batch protocol over the first 5 hours until reaching 30% of solids loading. After 6, 24 and 96 hours, samples were taken and characterized with respect to their sugar composition, rheology and water absorption. Successful slurry creation with dramatically reduced yield stress was achieved for corn stover for both assessed enzymes. Yield stresses of 178±7 Pa (3 FPU, Celluclast 1.5L) and 79±6 Pa (3 FPU, Ctec-2) were measured for corn stover at 24 hours, compared to 6,000 Pa for samples without enzyme. Yield stress was 155± 29 Pa (3FPU, Ctec-2) and 257 ± 72 Pa (1 FPU, Celluclast 1.5L) for corn cobs at 24 hours. Yield stress decreased when residence time increased with an enhanced fluidity noted for higher enzyme concentrations. A profile for 6, 24 and 96h of yield stress measurements is presented.

Gutierrez, Diana↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗
DOE OSTI2023

Lightning Rings and Gravity Waves: Insights Into the Giant Eruption Plume From Tonga's Hunga Volcano on 15 January 2022

Abstract On 15 January 2022, Hunga Volcano in Tonga produced the most violent eruption in the modern satellite era, sending a water‐rich plume at least 58 km high. Using a combination of satellite‐ and ground‐based sensors, we investigate the astonishing rate of volcanic lightning (>2,600 flashes min −1 ) and what it reveals about the dynamics of the submarine eruption. In map view, lightning locations form radially expanding rings. We show that the initial lightning ring is co‐located with an internal gravity wave traveling >80 m s −1 in the stratospheric umbrella cloud. Buoyant oscillations of the plume's overshooting top generated the gravity waves, which enhanced turbulent particle interactions and triggered high‐current electrical discharges at unusually high altitudes. Our analysis attributes the intense lightning activity to an exceptional mass eruption rate (>5 × 10 9 kg s −1 ), rapidly expanding umbrella cloud, and entrainment of abundant seawater vaporized from magma‐water interaction at the submarine vent.

58 GEOSCIENCES↗