Search NASA⌕ Search

SEARCH · Search NASA

Results for “explosive”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Methods for Color Center Preserving Hydrogen‐Termination of Diamond

Abstract Chemical functionalization of diamond surfaces by hydrogen is an important method for controlling the charge state of near‐surface fluorescent color centers, an essential process in fabricating devices such as diamond field‐effect transistors and chemical sensors, and a required first step for realizing families of more complex terminations through subsequent chemical processing. In all these cases, termination is typically achieved using hydrogen plasma sources that can etch or damage the diamond, as well as deposited materials or embedded color centers. This work explores alternative methods for lower‐damage hydrogenation of diamond surfaces, specifically the annealing of diamond samples in high‐purity, non‐explosive mixtures of nitrogen and hydrogen gas, and the exposure of samples to microwave hydrogen plasmas in the absence of intentional stage heating. The effectiveness of these methods are characterized by x‐ray photoelectron spectroscopy (XPS), and comparison of the results to density‐functional modelling of the surface hydrogenation energetics implicates surface oxygen ligands as the primary factor limiting the termination quality of annealed samples. Finally, photoluminescence (PL) spectroscopy is used to verify that both the annealing and reduced sample temperature plasma methods are non‐destructive to near‐surface ensembles of nitrogen‐vacancy (NV) centers, in stark contrast to plasma treatments that use heated sample stages.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Divergent Responses of Carbon Nitride Dot‐Based Amorphous Species and Small Molecule Hybrids to Trace Level Analytes

Bottom-up synthesis of carbon nitride dots (CNDs) offers a versatile platform for the creation of diverse nanomaterials with tunable properties. Here, we report a facile hydrothermal approach using citric acid (CA) and urea (U) as precursors to synthesize CNDs with varying degrees of condensation and crystallinity. By carefully controlling reaction conditions and post-synthetic treatments, we obtained two distinct fractions: a polycrystalline fraction composed of small-molecule hybrids and an amorphous fraction containing CNDs. We then sought to understand how the dominant species in these fractions impact sensing abilities using trace-level explosive exemplars. In conclusion, the results have important implications for sensing and related applications where understanding the complex interplay between synthetic conditions and post-synthetic processing play vital roles in determining the final properties of CND materials.

carbon nitride dots↗
DOE OSTI2026

Energetic Copolymers From LLM-105 and Aliphatic Isocyanates

Energetic polymers typically feature fuel-rich backbones with pendant oxidizing explosophores, which are used to modify pressure and temperature characteristics in energetic formulations. However, these pendant explosophores generally increase sensitivity and decrease the thermal stability of the polymer in the condensed phase. Direct polymerization of insensitive high explosives (IHEs) bearing polymerizable functionalities, such as amines, provides a synthetic platform for deriving tunable energetic materials. Here, this work demonstrates the first direct polymerization of the IHE 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) with aliphatic isocyanate comonomers to yield an energetic copolyurea (PUa1) and an energetic copolyurea network (PUa2). 1 H and 13 C nuclear magnetic resonance (NMR) and Fourier transform infrared spectroscopy (FTIR) confirm copolyurea synthesis. PUa1 exhibited branching reactions commonly found in polyurea syntheses, and these side reactions were suppressed in PUa2 through use of a polyfunctional isocyanate. Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) demonstrate the energetic IHE based polyurea preserves the energetic decomposition of the parent monomer LLM-105 in PUa1 (1203 J g −1 ) and PUa2 (687 J g −1 ), albeit with reduced thermal stability with peak decomposition temperatures of 235°C and 233°C, respectively. This work presents a new approach for generating energetic polymers from IHE cores with tunable energetic, thermal, and structural characteristics.

Chemistry↗
DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗
DOE OSTI2024

Nitrogen‐Nitrogen Bond Breaking in Irradiation Products of Hexanitrohexaazaisowurtzitane (CL‐20)

While the primary result of an interaction of ionizing radiation with an organic material is the ejection of electrons from molecular orbitals producing cations, the free electrons can further interact, yielding excited states; and, potentially, anions. In this work, we computed reaction barriers and energies for N−N bond breaking in the CL-20 cation, anion, and excited state to investigate if CL-20 degradation is accelerated after radiation exposure. We focused on N−N bond dissociation because it is the rate-determining step in the thermal degradation of CL-20. We found that N−N cleavage rapidly takes place in the CL-20 cation and anion with greatly reduced reaction energies as compared to CL-20. We also outlined a potential path for photodissociation of the N−N bond. While CL-20 is kinetically stable, initial degradation occurs readily in its irradiation products. In conclusion, this is of importance for the performance of the explosive after high-dose exposure and influences aging if CL-20 is irradiated at a low dose over extended periods of time.

CL-20 anion↗
DOE OSTI2025

Automated Scoring of Morphological Changes in Images of Pentaerythritol Tetranitrate

Recent advances in characterization techniques that generate large datasets of material microstructure images require robust, automated image-processing. We applied an unsupervised anomaly detection method called feature anomaly detection system (FADS) to automatically detect and quantify microstructure changes in images of the explosive pentaerythritol tetranitrate (PETN) aged at various temperatures. We demonstrated the FADS approach on two-dimensional images extracted from computed tomography scans, but the same technique can be readily applied to other imaging modalities. FADS calculates anomaly scores on the basis of differences in filter activations of nominal and test data in pretrained convolutional neural networks. The FADS scores successfully differentiated between pristine PETN and PETN aged at a temperature where material coarsening occurred. Morphological metric analysis of segmented images verified observed trends in FADS scores as a function of aging temperature and aging time, specifically by calculating volume fractions, specific boundary lengths, two-point correlation functions, and local thicknesses. Here, the FADS technique has two important advantages compared to traditional morphological analysis: First, it uses grayscale images as input, rather than images that are segmented to separate the appropriate phases; and second, FADS scores capture any type of changes among image sets, rather than requiring prior knowledge or selection of a relevant set of metrics.

Accelerated aging↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gamma-Ray and Cosmic Ray Muon Modalities for Cargo Inspection

Screening and inspection of cargo containers are two essential methods to nondestructively examine the contents of shipment. These methods enable the detection of illicit transportation of unauthorized materials such as nuclear and radioactive materials, explosives, drugs, and so on, typically at borders or secure facilities. Although high-energy X-ray transmission is a standard system and is widely used for cargo inspection, the inherent challenges of high false-positive rates and high attenuation factors necessitate the development of complementary techniques that can increase the detection efficiency and accuracy in large and dense materials. Gamma-rays, which possess higher penetration characteristics because of their high energy, offer an alternative nonintrusive modality for cargo scanning. They represent a promising inspection method when compared to X-rays for three reasons: (1) improved ability to detect nuclear and radioactive materials, (2) higher inspection throughput rates, and (3) lower false-positive rates. Currently, there are two main gamma-ray inspection techniques, active and passive interrogation. Active interrogation can be further grouped into (1) gamma-ray transmission imaging and (2) neutron-induced gamma-ray emission detection. Gamma-ray transmission imaging utilizes differences in material densities for mapping the shipment contents and detecting anomalies. It is analogous to the X-ray transmission method; however, the high-energy photons make it more difficult to shield against, which enables more efficient performance in large and dense material inspection. Neutron-induced gamma-ray emission inspection is designed for the detection of nuclear and radioactive material because those materials emit characteristic gamma-rays when they are activated by neutron absorption. On the other hand, passive interrogation techniques rely on high-efficiency detectors to detect radiation emitted from hidden special nuclear or other radioactive materials. Similar to passive interrogation, cosmic ray muon monitoring and imaging are relatively new techniques that do not require external radioactive sources. These techniques have received attention as a potential next-generation radiographic probe to identify illicit transportation of nuclear and radioactive materials in cargo containers. Cosmic ray muons have unique features, (1) much higher energies than X-rays or gamma-rays (on the order of 10−1—104 GeV), (2) enhanced penetration capability, and (3) natural occurrence, thereby eliminating the need for induced radiation sources. These features enable cosmic ray muons to be utilized for detection of special nuclear materials in high-background-noise environments. By analyzing incoming and outgoing muon trajectories, scattering angles, and energies, it has been shown that it would be possible to locate hidden and well-shielded materials in cargo containers via three-dimensional muon tomography images or signal analysis. Gamma-rays, cosmic ray muons, and other nonintrusive cargo inspection modalities are complementary to each other, allowing them to address various cargo inspection conditions (i.e., scanning time, cost, radiation exposure level, and types of target materials). This chapter presents a detailed review of the theoretical fundamentals and technical principles behind the current gamma-ray and cosmic ray muon modalities for cargo inspection. Additionally, critical assessments and suggestions for the future directions to advance the use of gamma and muon modalities are discussed.

Bae, Junghyun↗
DOE OSTI2024

Bosenovae with quadratically-coupled scalars in quantum sensing experiments

Abstract Ultralight dark matter (ULDM) particles of massm ϕ ≲ 1 eV can form boson stars in DM halos. Collapse of boson stars leads to explosive bosenova emission of copious relativistic ULDM particles. In this work, we analyze the sensitivity of terrestrial and space-based experiments to detect such relativistic scalar ULDM particles interacting through quadratic couplings with Standard Model constituents, including electrons, photons, and gluons. We highlight key differences with searches for linear ULDM couplings. Screening of ULDM with quadratic couplings near the surface of the Earth can significantly impact observations in terrestrial experiments, motivating future space-based experiments. We demonstrate excellent ULDM discovery prospects, especially for quantum sensors, which can probe quadratic couplings orders below existing constraints by detecting bosenova events in the ULDM mass range 10 −23 eV ≲m ϕ ≲ 10 −5 eV. We also report updated constraints on quadratic couplings of ULDM in case it comprises cold DM.

Physics↗
DOE OSTI2025

Effect of a Fine-Scale Layered Structure of the Atmosphere on Infrasound Signals from Fragmenting Meteoroids

We investigate the influence of a fine-scale (FS) layered structure in the atmosphere on the propagation of infrasound signals generated by fragmenting meteoroids. Using a pseudo-differential parabolic equation (PPE) approach, we model broadband acoustic signals from point sources at altitudes of 35–100 km. The presence of FS fluctuations in the stratosphere (37–45 km) and the lower thermosphere (100–120 km) modifies ray trajectories, causing multiple arrivals and prolonged signal durations at ground stations. In particular, meteoroids fragmenting at 80–100 km can produce two distinct thermospheric arrivals beyond 150 km range, while meteoroids descending to 50 km or below yield weak, long-lived arrivals within the acoustic shadow zone via antiguiding propagation and diffraction. Comparison with observed infrasound data confirms that FS-layered inhomogeneities can account for multi-arrival “N-waves,” broadening potential interpretations of meteoroid signals. The results also apply to other atmospheric-entry objects, such as sample return capsules, emphasizing how FS structure impacts shock wave propagation. In conclusion, our findings advance understanding of wavefield evolution in a layered atmosphere and have broad relevance for global infrasound monitoring of diverse phenomena (e.g., re-entry capsules, rocket launches, and large-scale explosions).

Aeroacoustics↗
DOE OSTI2024

Mechanical and poroelastic behavior of porous tuff under drained and undrained conditions

A series of drained and undrained water-saturated constant mean-stress tests were performed to investigate the strength, elasticity, and poroelastic response of a water-saturated high porosity nonwelded tuff. Drained strengths are found to increase with increasing effective confining pressures. Elastic moduli increase with increasing mean stress. Undrained strengths are small due to development of high pore pressures that generate low effective confining pressures. Skempton’s values are pressure dependent and appear to reflect the onset of inelastic deformation. Permeabilities decrease after deformation from ~ 10 –14 to ~ 10 –16 m 2 and are a function of the applied confining pressure. Deformation is dominated by pore collapse, compaction, and intense microfracturing, with the undrained tests favoring microfracture-dominant deformation and the drained tests favoring compaction-dominant deformation. These property determinations and observations are used to develop/parameterize physics-based models for underground explosives testing.

58 GEOSCIENCES↗
DOE OSTI2026

Mixed-precision numerics in scientific applications: survey and perspectives

The explosive demand for artificial intelligence (AI) workloads has led to a significant increase in silicon area dedicated to lower-precision computations on recent high-performance computing hardware designs. However, mixed-precision capabilities, which can achieve performance improvements of up to 8x compared to double-precision in extreme compute-intensive workloads, remain largely untapped in most scientific applications. A growing number of efforts have shown that mixed-precision algorithmic innovations can deliver superior performance without sacrificing accuracy. These developments should prompt computational scientists to seriously consider whether their scientific modeling and simulation applications could benefit from the acceleration offered by new hardware and mixed-precision algorithms. In this survey, we (1) review progress across diverse scientific domains—fluid dynamics, weather and climate, quantum chemistry, and computational genomics—that have begun adopting mixed-precision strategies; (2) examine state-of-the-art algorithmic techniques such as iterative refinement, splitting and emulation schemes, and adaptive precision solvers; (3) assess their implications for accuracy, performance, and resource utilization; and (4) survey the emerging software ecosystem that enables mixed-precision methods at scale. We conclude with perspectives and recommendations on cross-cutting opportunities, domain-specific challenges, and the role of co-design between application scientists, numerical analysts, and computer scientists. Collectively, this survey underscores that mixed-precision numerics can reshape computational science by aligning algorithms with the evolving landscape of hardware capabilities.

Graphics processing units↗
DOE OSTI2026

NASA analysis of space mission options for the 2025 planetary defense conference hypothetical asteroid impact threat scenario

The 2025 Planetary Defense Conference (PDC) hypothetical asteroid impact threat exercise is being conducted in coordination with the United Nations-endorsed Space Mission Planning Advisory Group (SMPAG), exercising SMPAG's processes for assessing space mission options and communicating to decision makers. Here, in this paper, we describe the work performed by our NASA-led team and present the results we contributed to the SMPAG effort for the exercise. Our NASA-led team assessed mission options for asteroid reconnaissance (flyby and rendezvous), deflection, and robust disruption using several viable techniques: kinetic impactors, ion beam deflection, and nuclear explosive devices. Our simulations and analyses considered how much change-in-velocity an asteroid can tolerate before fragmentation onset, and we worked towards establishing requirements for robustly disrupting an asteroid. Heuristics informed by the simulation results are incorporated into optimization of deflection and disruption mission campaign options. Finally, we make some observations about useful generalizations from these results with potential applicability to any planetary defense scenario.

Asteroid deflection↗
DOE OSTI2025

127 Xe quantification method development and intercomparison exercise

Monitoring of the atmosphere for fission products ( 131m Xe, 133m Xe, 133 Xe, and 135 Xe) is performed by various laboratories to detect nuclear explosions. Quantification of 127 Xe is not routinely performed by laboratories measuring atmospheric radioxenon because it is not a fission product. 127 Xe was recently detected by a ground-based beta-gamma air monitoring system. When measured using beta-gamma coincidence detector systems, such as those in use on the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty (CTBT), 127 Xe can interfere with the quantification of fission product radioxenon due to overlap of the 127 Xe beta-gamma coincidence signatures with those of fission product radioxenon. Here, this work demonstrates quantification of 127 Xe at different laboratories with different measurement techniques. Production and purification of 127 Xe was performed by neutron activation of enriched 126 Xe. The purified 127 Xe was then split between laboratories, and detection and quantification methods were developed. At Idaho National Laboratory, a quantification method involving high purity germanium detectors was devised that included self-attenuation correction. At AWE, a beta-gamma coincidence counting method, as used in support of the IMS, was modified to enable the measurement and analysis of the 127 Xe samples. Corrections were made for self-attenuation, which showed a strong xenon volume dependency, for some coincidence signatures. The gas sample activity concentration was used as the comparison metric and it showed excellent agreement between the methods.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
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 OSTI2026

Improving anaerobic digestion of sewage sludge to renewable natural gas by the Advanced Pretreatment & Anaerobic Digestion technology (APAD): Pilot testing

Conventional anaerobic digestion (AD) of sewage sludge in wastewater treatment facilities suffers from low carbon conversion efficiency (CCE = 40%) and requires costly CO2 removal for injection of the produced CH4 into the natural gas grid. To address these limitations, we developed the Advanced Pretreatment and Anaerobic Digestion (APAD) process. This integrates Advanced Wet Oxidation & Steam Explosion (AWOEx) pretreatment of residual sludge after conventional AD, followed by biogas upgradation using a novel methanogenic strain, Methanothermobacter wolfeii BSEL, converting CO2 with H2 into CH4 or RNG (renewable natural gas). Pilot-scale results demonstrated that AWOEx pretreatment achieved a CCE of 62% for the residual sludge, 68% higher than the conventional AD process. The CH4 production was further increased by 79%. Subsequent biogas upgrading in a trickling bed reactor with H2 further enhanced total methane output by 100% and resulted in a final CO2 concentration of =3%. The integrated APAD process achieved a remarkable overall CCE of 83%, resulting in a 200% increase in RNG output when compared to conventional AD. Techno-economic analysis revealed that AWOEx pretreatment alone reduced sludge treatment costs from $494 to $253 per ton of dry solids. The complete APAD process incurred a higher cost of treatment of $530 per ton, driven by prices for bottled H2. The process did, however, show gains in energy recovery and decarbonization. Renewable H2, which may reduce in price in the near future, can positively improve the economics of biogas upgrading for the APAD process.

Life Cycle Assessment (LCA)↗
DOE OSTI2025

Computational diagnostics and characterization of combustion recession in diesel sprays

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

CEMA↗
DOE OSTI2025

In situ multi-tier auto-ignition detection applied to dual-fuel combustion simulations

Here we use an anomaly detection methodology that is centered on analyzing fourth-order joint moments (co-kurtosis), particularly focusing on its application in auto-ignition of combustion problems with large numbers of species. Unsupervised anomaly detection is challenging to generalize across problem types and domains. A recent technique, centered on analyzing information in the fourth-order joint moment co-kurtosis, has shown promise, especially for high-dimensional scientific data. In this work we present developments to the co-kurtosis based anomaly detection method needed to make it effective and scalable for large-scale distributed scientific data, such as those generated by massively parallel simulations. An in situ co-kurtosis algorithm is employed as the anomaly detection method for identifying ignition kernels in simulations of turbulent combustion. Here, we extend an existing methodology which identifies regions of the domain where anomalies are present, and add another tier of anomaly detection where the individual samples contributing to the anomaly are identified. We apply this algorithm on-the-fly to a variety of turbulent reacting flow problems and compare it to the widely used (but significantly more expensive) chemical explosive mode analysis (CEMA). We demonstrate the ability of the method to detect and identify the onset of low and high temperature ignition which can be used for computational steering, as chemical and combustion anomalies occur intermittently at spatio-temporal locations unknown a priori. Finally, we apply our lightweight in situ algorithm to an exascale high-fidelity simulation with a total of 2.4 Trillion degrees of freedom, performed using an adaptive mesh refinement solver. Furthermore, through a scalability analysis, we show that the relative computational cost of this in-situ anomaly detection algorithm compared to an iteration of the reacting flow solver is negligible.

97 MATHEMATICS AND COMPUTING↗