Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal network model”

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

Accelerated computation of lattice thermal conductivity using neural network interatomic potentials

We report with the development of the density functional theory (DFT) and ever-increasing computational capacity, an accurate prediction of lattice thermal conductivity based on the Boltzmann transport theory becomes computationally feasible, contributing to a fundamental understanding of thermal conductivity as well as a choice of the optimal materials for specific applications. However, steep costs in evaluating third-order force constants limit the theoretical investigation to crystals with high symmetry and few atoms in the unit cell. Currently, machine learning potentials are garnering attention as a computationally efficient high-fidelity model of DFT, and several studies have demonstrated that the lattice thermal conductivity could be computed accurately via the machine learning potentials. However, test materials were mostly crystals with high symmetries, and the applicability of machine learning potentials to a wide range of materials has yet to be demonstrated. Furthermore, establishing a standard training set that provides consistent accuracy and computational efficiencies across a wide range of materials would be useful. To address these issues, herein we compute lattice thermal conductivities at 300 K using neural network interatomic potentials. As test materials, we select 25 materials with diverse symmetries and a wide range of lattice thermal conductivities between 10 -1 and 10 3 Wm -1 K -1 . Among various choices of training sets, we find that molecular dynamics trajectories at 50–700 K consistently provide results at par with DFT for the test materials. In contrast to pure DFT approaches, the computational cost in the present approach is uniform over the test materials, yielding a speed gain of 2–10 folds. When a smaller reduced training set is used, the relative efficiency increases by up to ~50 folds without sacrificing accuracy significantly. The current work will broaden the application scope of machine learning potentials by establishing a robust framework for accurately computing lattice thermal conductivity with machine learning potentials.

36 MATERIALS SCIENCE↗

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural and Optical Response of Polymer-Stabilized Blue Phase Liquid Crystal Films to Volatile Organic Compounds

Engineering useful mechanical properties into stimuli-responsive soft materials without compromising their responsiveness is, in many cases, an unresolved challenge. For example, polymer networks formed within blue-phase liquid crystals (BPs) have been shown to form mechanically robust films, but the impact of polymer networks on the response of these soft materials to chemical stimuli has not been explored. In this work, we report on the response of polymer-stabilized BPs (PSBPs) to volatile organic compounds (VOCs, using toluene as a model compound) and compare the response to BPs without polymer stabilization and to polymerized nematic and cholesteric phases. We find that PSBPs generate an optical response to toluene vapor (change in reflection intensity under crossed polars) that is sixfold greater in sensitivity than the polymerized nematic or cholesteric phases and with a limit of detection (140 ± 10 ppm at 25 °C) that is relevant to the measurement of permissible exposure limits for humans. Additionally, when compared to BPs that have not been polymerized, PSBPs respond to a broader range of toluene vapor concentrations (5000 vs <1000 ppm) over a wider temperature interval (25–45 vs 45–53 °C). We place these experimental observations into the context of a simple thermodynamic model to explore how the PSBP response reflects the effect of toluene on competing contributions of double-twisted LC cylinders, disclinations, and polymer network to the free energy that controls the PSBP lattice spacing. Overall, we conclude that the mechanical and thermal stability of PSBPs, when combined with their optical responsiveness to toluene, make this class of self-supporting LCs a promising one as the basis of passive and compact (e.g., wearable) sensors for VOCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of Analytical and Numerical Performance Predictions for an International Space Station Node 3 Internal Active Thermal Control System Regenerative Heat Exchanger

The complexity of International Space Station (ISS) systems modeling often necessitates the concurrence of various dissimilar, parallel analysis techniques to validate modeling. This was the case with a feasibility and performance study of the ISS Node 3 Regenerative Heat Exchanger (RHX). A thermo-hydraulic network model was created and analyzed in SINDA/FLUINT. A less complex, closed form solution of the systems dynamics was created using an Excel Spreadsheet. The purpose of this paper is to provide a brief description of the modeling processes utilized, the results and benefits of each to the ISS Node 3 RHX study.

Wise, Stephen A.↗

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗

Primordial black hole dark matter from catastrogenesis with unstable pseudo-Goldstone bosons

We propose a new scenario for the formation of asteroid-mass primordial black holes (PBHs). Our mechanism is based on the annihilation of the string-wall network associated with the breaking of a U(1) global symmetry into a discrete Z N symmetry. If the potential has multiple local minima ( N > 1) the network is stable, and the annihilation is guaranteed by a bias among the different vacua. The collapse of the string-wall network is accompanied by catastrogenesis , a large production of pseudo-Goldstone bosons (pGBs) — e.g. axions, ALPs, or majorons — gravitational waves, and PBHs. If pGBs rapidly decay into products that thermalize, as predicted e.g. in the high-quality QCD axion and heavy majoron models, they do not contribute to the dark matter population, but we show that PBHs can constitute 100% of the dark matter. The gravitational wave background produced by catastrogenesis with heavy unstable axions, ALPs, or majorons could be visible in future interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗

Isomerization kinetics of azobenzene crosslinked to a liquid crystal polymer network

Azobenzene is an important chromophore having numerous applications that are a direct implication of its ability to change shape upon irradiation with a suitable wavelength. Here, in this study, photoisomerization kinetics of a diacrylate azobenzene molecule (A3MA) cross-linked to a liquid crystal polymer network (LCN) is investigated. Isomerization experiments were conducted on planar LCN thin films with varying cross-linking densities for a fixed azobenzene concentration. The experiments were conducted at approximately 26 °C, which is significantly below each film’s glass transition temperature. The forward trans-cis and reverse cis-trans isomerization slows down in crosslinked azobenzene liquid crystal polymer networks (ALCNs) as compared to a solution and an uncrosslinked polymer mixture. The observed slowdown is due to the reduction in the free volume and increased steric hindrances with an increase in crosslinking. A single-term first-order kinetics rate equation is known to describe the isomerization kinetics in solutions or uncrosslinked mixtures while the same form fails to represent the reaction kinetics in moderate to highly cross-linked ALCNs. Our results for the first time indicate that a two-term first-order reaction rate is an accurate way to describe the forward trans-cis and reverse thermal cis-trans isomerization in crosslinked azobenzene liquid crystal polymer networks. In addition, the effect of film thickness and incident intensities on the reaction kinetics is investigated. Finally, the two-term isomerization kinetics is incorporated into a finite element photo-chemo-mechanical model to reveal the consequent slowdown in the actuation response of ALCN films. The current findings will be of particular interest to theoreticians and experimentalists alike, for it will help predict and understand accurately the responses of azobenzene-based coatings, and actuators and may have potential applications in underwater soft robotics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

1979 January 13 - An intense gamma-ray burst with a possible associated optical transient

The time history, energy spectra, and precise localization of the intense January 13, 1979 gamma-ray burst, observed by five spacecraft in the interplanetary network, are presented. The time history displays a pulse-afterpulse structure suggestive of periodicity. Fine time resolution spectral analysis shows that a wide variety of models can be used to fit individual spectra, while spectra integrated over longer periods are well fitted by a thermal synchrotron law, making it difficult to identify any specific emission mechanism. The precise localization may be consistent with an archival optical transient having no quiescent counterpart down to 22d mag on recent plates.

Barat, C.↗

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

The Effect of Grain Boundary Facet Junctions on Segregation and Embrittlement

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

36 MATERIALS SCIENCE↗

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Photometric Studies of Heat Flow at the Photosphere

Continuum photometry is carried out and the results of comparing these observations with models of photospheric heat flow are described. The main results are: (1) a possible detection of weak bright rings around some spot penumbrae; (2) no evidence is found for large scale photospheric brightness inhomogeneities exceeding 2-3 K which places tighter constraints on models of global scale convection; (3) supergranular scale continuum structures observed across the photosphere appear mainly due to random clumping of granules; (4) the one case observed of a sunspot emergence shows no thermal shadow exceeding 1.5 K rms one day prior to umbra appearance; (5) network and faculae are found to show a small excess brightness even at mu = 1, so detection of faculae at mu = 1 by differential photometry indicates a gentler temperature gradient near tau = 1 in the facular (relative to cell) atmosphere; (6) the limb darkening study shows no significant global variations to within 0.1% rms.

Foukal, P.↗

Investigation of mesoscale pre-convective forcing mechanisms

The installation and testing of a numerical mesoscale boundary layer model is discussed. A data analysis is first made using the rawinsonde network which shows a strong baroclinic zone developing in the region where convection was initiated. The observed wind field shows an easterly component developing in opposition to the synoptic flow producing convergence. The numerical boundary layer model is used to diagnose the physical development of the thermal and flow fields. Two physical mechanisms are examined; topography and cloud shading. The results shows that while heated sloping topography alone does create a baroclinic zone and convergence values, cloud shading is a stronger mechanism for producing a thermal field representative of observed conditions

Mcnider, R. T.↗

NASA Tech Briefs, April 2002

The contents include: 1) Application Briefs; 2) Sneak Preview of Sensors Expo; 3) The Complexity of the Diagnosis Problem; 4) Design Concepts for the ISS TransHab Module; 5) Characteristics of Supercritical Transitional Mixing Layers; 6) Electrometer for Triboelectric Evaluation of Materials; 7) Infrared CO2 Sensor With Built-In Calibration Chambers; 8) Solid-State Potentiometric CO Sensor; 9) Planetary Rover Absolute Heading Detection Using a Sun Sensor; 10) Concept for Utilizing Full Areas of STJ Photodetector Arrays; 11) Development of Cognitive Sensors; 12) Enabling Higher-Voltage Operation of SOl CMOS Transistors; 13) Estimating Antenna-Pointing Errors From Beam Squints; 14) Advanced-Fatigue-Crack-Growth and Fracture- Mechanics Program; 15) Software for Sequencing Spacecraft Actions; 16) Program Distributes and Tracks Organizational Memoranda; 16) Flat Membrane Device for Dehumidification of Air; 17) Inverted Hindle Mount Reduces Sag of a Large, Precise Mirror; 18) Heart-Pump-Outlet/Cannula Coupling; 19) Externally Triggered Microcapsules Release Drugs In Situ; 20) Combinatorial Drug Design Augmented by Information Theory; 21) Multiple-Path-Length Optical Absorbance Cell; 22) Model of a Fluidized Bed Containing a Mixture of Particles; 23) Refractive Secondary Concentrators for Solar Thermal Systems; 24) Cold Flow Calorimeter; 25) Methodology for Tracking Hazards and Predicting Failures; 26) Estimating Heterodyne-Interferometer Polarization Leakage; 27) An Efficient Algorithm for Propagation of Temporal- Constraint Networks; 28) Software for Continuous Replanning During Execution; 29) Surface-Launched Explorers for Reconnaissance/Scouting; 30) Firmware for a Small Motion-Control Processor; 31) Gear Bearings and Gear-Bearing Transmissions; and 32) Linear Dynamometer With Variable Stroke and Frequency.

Source record↗

Benchmarking materials property prediction methods: the Matbench test set and Automatminer reference algorithm

Abstract We present a benchmark test suite and an automated machine learning procedure for evaluating supervised machine learning (ML) models for predicting properties of inorganic bulk materials. The test suite, Matbench, is a set of 13 ML tasks that range in size from 312 to 132k samples and contain data from 10 density functional theory-derived and experimental sources. Tasks include predicting optical, thermal, electronic, thermodynamic, tensile, and elastic properties given a material’s composition and/or crystal structure. The reference algorithm, Automatminer, is a highly-extensible, fully automated ML pipeline for predicting materials properties from materials primitives (such as composition and crystal structure) without user intervention or hyperparameter tuning. We test Automatminer on the Matbench test suite and compare its predictive power with state-of-the-art crystal graph neural networks and a traditional descriptor-based Random Forest model. We find Automatminer achieves the best performance on 8 of 13 tasks in the benchmark. We also show our test suite is capable of exposing predictive advantages of each algorithm—namely, that crystal graph methods appear to outperform traditional machine learning methods given ~10 4 or greater data points. We encourage evaluating materials ML algorithms on the Matbench benchmark and comparing them against the latest version of Automatminer.

36 MATERIALS SCIENCE↗

Recent Developments at the Raft River Geothermal Field

The Raft River geothermal field, located in Cassia County in southwestern Idaho, is the site of a Department of Energy Enhanced Geothermal System project. U.S. Geothermal, Inc. currently produces about 11 MWe from Precambrian metamorphic rocks. These lie beneath ~5,000 ft of Quaternary and Tertiary volcaniclastic and volcanic deposits. Maximum temperatures range from 271°F to 300°F. Well RRG-9 ST1, the well targeted for stimulation is located approximately 1 mile south of the main bore field. The open hole section of the well, from 5,551 to 5,900 ft MD, consists of Precambrian Elba Quartzite, the stimulation target, granite and minor diabase. Prior to setting the casing acoustic, gamma ray, and density logs were run. After completing the well, a step rate/step down test was conducted. The maximum injection rate achieved was 18 bpm at a wellhead pressure of 1,150 psig. A borehole televiewer run in the open hole section showed evidence of more than eighty fractures. The majority of these fractures trend from N20⁰W to N20°E and dip from 40° to 60°W. Permeable fractures were encountered in the Elba Quartzite at 5,640-5,660 ft MD. Analysis of the injection test indicates that the minimum in-situ principal stress in this zone is 3,050-3,200 psi, corresponding to a fracture gradient of 0.59-0.62 psi/ft. A discrete fracture network model was developed using measured and inferred fracture orientations, distributions and dimensions. A three-phase stimulation program is proposed for RRG-9 ST-1. During the first two stages, water at 140°F, and later 40°F, will be injected to pre-condition and thermally fracture the reservoir. Here, the third stage will consist of a high rate, large volume conventional hydraulic stimulation.

Bradford, Jacob↗

Features of a SINDA/FLUINT model of a liquid oxygen supply line

The modeling features used in a steady-state heat transfer problem using SINDA/FLUINT are described. The problem modeled is a 125 feet long, 3 inch diameter pipe, filled with liquid oxygen flow driven by a given pressure gradient. The pipe is fully insulated in five sections. Three sections of 1 inch thick spray-on foam and two sections of vacuum jacket. The model evaluates friction, turns losses and convection heat transfer between the fluid and the pipe wall. There is conduction through the foam insulation with temperature dependent thermal conductivity. The vacuum space is modeled with radiation and gas molecular conduction, if present, in the annular gap. Heat is transferred between the outer surface and surrounding ambient by natural convection and radiation; and, by axial conduction along the pipe and through the vacuum jacket spacers and welded seal flanges. The model makes extensive use of SINDA/FLUINT basic capabilities such as the GEN option for nodes and conductors (to generate groups of nodes or conductors), the SIV option (to generate single, temperature varying conductors), the SIM option (for multiple, temperature varying conductors) and the M HX macros for fluids (to generate strings of lumps, paths, and ties representing a diabatic duct). It calls subroutine CONTRN (returns the relative location in the G-array of a network conductor, given an actual conductor number) enabling an extensive manipulation of conductor (calculates an assignment of their values) with DO loops. Models like this illustrate to the new and even to the old SINDA/FLUINT user, features of the program that are not so obvious or known, and that are extremely handy when trying to take advantage of both, the automation of the DATA headers and make surgical modifications to specific parameters of the thermal or fluid elements in the OPERATIONS portion of the model.

Simmonds, Boris G.↗

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗