Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction rates”

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

ZND Verification Tests for Reactive Burn Models in FLAG

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

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

Kinetics of TATB and LX-17 Decomposition by Thermal Analysis

The kinetics of TATB decomposition have been investigated by thermal analysis (thermogravimetric analysis, TGA, and differential scanning calorimetry, DSC) at both LLNL and the University of Alabama Birmingham. This report summarizes progress in understanding from those results within the context of identifying the global reaction mechanism, i.e., sequential, autocatalytic, competitive, or some mixture thereof. A particularly important issue is how pressure affects both the reaction rate and enthalpy. The effect of pressure helps identify dominant mechanisms and construct and constrain a global mechanism that will work over a wide range of conditions. This memorandum is simply a snapshot of current progress and not a definitive interpretation. I have written it to organize work of others and my interpretations in a systematic manner with the intention of stimulating further advances in understanding. SDT data comes from Jen Montgomery at LLNL, and high-pressure DSC comes from Victoria Stanford and Sergey Vyazovkin at University of Alabama at Birmingham. Major observations and conclusions, some new and some not, are (1) outgassing below 150 °C prior to decomposition releases adsorbed water vapor, (2) outgassing between 150 and 250 °C may come primarily from intercalated water (3) initial thermal decomposition above 250 °C or so involves dehydration of TATB to create furazans, (4) quickly heating to higher temperatures in open pans with flowing gas leads primarily to TATB evaporation between 270 and 390 °C, (5) using a pinhole crucible at atmospheric pressure largely eliminates TATB sublimation, (6) at typical thermal analysis heat rates, the main stage of decomposition is described by a single mass loss process, and decomposition of residual material is spread over a wide range of thermal severity, (7) at typical thermal analysis heat rates and constant temperatures above 300 °C, heat release occurs in two overlapping peaks, with the first one a little faster than mass loss and the second on a little slower than mass loss, (8) total heat release increases significantly with heating rate and substantially with increased pressure, and (9) isothermal experiments at 340 °C indicate that the first exotherm increases by more with LX-17 than TATB, which suggests a possible chemical role of Kel-F decomposition.

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

Determination of Solar System R-Process Abundances using ENDF/B-VIII.0 and TENDL-2015 libraries

Recent multi-messenger detection of the binary neutron star merger (GW170817) energized the astrophysical community and encouraged further research for determination of nuclear physics observables. Comprehensive studies of atomic nuclei in the cosmos provide an opportunity for investigating these astrophysical phenomena and acquiring complementary information on stellar nucleosynthesis processes that can be verified using the latest nuclear data. Evaluated Nuclear Data File (ENDF) libraries contain complete collections of reaction cross sections over the energy range relevant to astrophysics, fission yields and decay data. These data collections have been used worldwide in nuclear science, industry and national security applications. There is great interest in exploring the ENDF/B-VIII.0 and TALYS Evaluated Nuclear Data Library (TENDL-2015) for nuclear astrophysics purposes and comparing findings with the Karlsruhe Astrophysical Database of Nucleosynthesis in Stars (KADoNiS). The Maxwellian-averaged cross sections (MACS) and astrophysical reaction rates have been calculated using the ENDF/B-VIII.0 and TENDL-2015 evaluated data sets. The calculated cross sections were combined with the solar system abundances and fitted using the classical model of stellar nucleosynthesis. Astrophysical rapid- and slow-neutron capture, r- and s-process, respectively, abundances were obtained from present data and compared with available values. Further analysis of MACS reveals potential evaluated libraries data deficiencies and a strong need for new measurements. The current results demonstrate a large nuclear astrophysics potential of evaluated libraries and mutually beneficial relations between nuclear industry and research e orts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

DRACO: An Overview [Slides]

DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.

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

0D Modeling of Reactor Networks Within Zuzax

We report on progress in developing macroscopic balance equations for combustion and electrochemistry systems. A steady state solution capability is described for the macroscopic reactor network, with an associated steady state continuation method and solution storage capability added in. An example is provided of continuation of a hydrogen flame versus the equivalence ratio. The reactor modeling capability is extended to charged fluid systems, with a description of the new ChargedFluidReactor, SubstrateElement, and MetalCurrentElement reactor classes and novel setup of unknowns within these reactors that preserve charge neutrality. Zuzax's setup for electrochemistry is explained including the specification of the electron chemical potential and the adherence to the SHE Reference electrode specification. The description of the different ways to enter electrochemical reaction rates are described, contrasted, and their derivations with respect to one another are derived. An example of using the ChargedFluidReactor within corrosion problems is provided. We present a description of calculations to understand the phenomena of corrosion of copper from a micron sized droplet of NaCl water droplet, where secondary spreading occurs. An analysis of the discrepancies with experiment is carried out, demonstrating that macroscopic balances can be an important tool for understanding what major factors need to be addressed for a better understanding of a physical system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗
DOE OSTI2020

Opportunities to Update the Model of Tank Closure Grout Aging (SRNL-STI-2012-00404) Based on Experimental Results by the Savannah River Ecology Laboratory

The purpose of this report was to review the results of experiments on evolution of tank grout pore fluids performed by the Savannah River Ecology Laboratory (SREL) and reported in Seaman et al. (2020) and highlight opportunities for updating the model of tank grout degradation presented in Denham and Millings (2012). The SREL experimental results validate some assumptions used in the 2012 grout degradation model and provide information useful for updating the model. All experiments on tank grout degradation have limitations because of reaction rates and the enormous disparity between duration of the experiments and the duration of the degradation of the real tank grout system. Given these limitations, the experimental results form the basis for several recommendations for updating the grout degradation model described in this report.

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

Sensitivity/Uncertainty Analysis for Nuclear Criticality Safety Validation [Slides]

Sensitivity coefficients predict the expected change in a response ($\kappa$) due to a change in some input parameter ($\Sigma$). These responses can be k eff (or ratios of reaction rates). These input parameters are typically nuclear data. The coefficients are dimensionless ratios. What would happen to the system k eff if some piece of data were changed by some amount? The coefficient is calculated without making the change.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Criticality and Radiation Shielding Impacts of ENDF/B-VIII.0 [Slides]

Several impacts of ENDF/B-VIII.0 on Radiation Shielding Applications were observed by researchers. Most reaction rates, attenuation ratios, and dose rates are similar for ENDF/B-VII.1 and ENDF/B-VIII.0. A significant reduction in reactor dosimetry results was noted. Investigation continues to reveal underlying data for discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

LLNL FY22 Aging and Lifetimes Exit Criteria (Milestone 8118, GC#5)

This report introduces a novel characterization tool that can rapidly collect in-situ time-dependent oxide thickness and hydride growth data together with spatial data using while-light interferometry. White-light interferometry is a non-destructive and non-contact optical surface profiling technique for materials characterization that offers excellent lateral and vertical resolution. A depleted uranium sample was used to assess the ability of white-light interferometry to provide quantitative corrosion reaction rates at metal surfaces for by hydriding and oxidation and compare this technique to existing characterization tools for capturing early-stage hydride and oxide growth kinetics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Producing ENDF/B-quality Evaluations of 239 Pu(n,f) and 235 U(n,f) Average Prompt Neutron Multiplicities using the CGMF Model

This report documents evaluations of the neutron-induced 239 Pu and 235 U average prompt neutron multiplicities, $\overline{v}_p$, using the CGMF model developed at LANL. These evaluations are not updates of existing ENDF/B-VIII.0 nuclear data, but were both re-done from "scratch". That is, all experimental data were extracted from EXFOR, re-normalized to the newest nuclear data representing monitor observables, and detailed uncertainties were estimated from information in EXFOR as well as from templates of expected measurement uncertainties. We also included experimental data that were not yet available for ENDF/B-VIII.0, for instance, the data of Marini et al. for 239 Pu(n,f) $\overline{v}_p$ or the data of Khoklov et al. for 235 U(n,f) $\overline{v}_p$. Finally, we include the CGMF model; it computes $\overline{v}_p$ based on model parameters that at the same time calculate fission fragments as a function of mass, Y (A), the average total kinetic energy as a function of incident-neutron energy, $\langle TKE \rangle$(E inc ), etc. Such a detailed fission model that ties together many fission quantities has not been used to date for any $\overline{v}_p$ evaluation in ENDF/B. Here, we show that we can get evaluated 235 U(n,f) and 239 Pu $\overline{v}_p$ that not only correspond well to experimental $\overline{v}_p$, but that the associated evaluated model parameters also yield parameterizations of Y (A), $\langle TKE \rangle$(E inc ), etc., that correspond well to their respective experimental data. The new evaluated 239 Pu(n,f) $\overline{v}_p$ was also combined into a 239 Pu test file with the newest nuclear data for the 239 Pu prompt-fission neutron spectrum and fission cross sections. This new 239 Pu file performed reasonably well in predicting PU-MET-FAST assemblies, reaction rates in Jezebel and Flattop and three LLNL pulsed spheres. Due to that, the new 239 Pu(n,f) $\overline{v}_p$ evaluation presented here is being considered for ENDF/B-VIII.1. Hence, we can show here that CGMF is able to produce ENDF/B quality $\overline{v}_p$ nuclear data.

239Pu, 235U, Average Prompt-fission Neutron Multip↗
DOE OSTI2022

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters [Slides]

Sensitivity coefficients of prompt neutron decay constant (α) and effective delayed neutron fraction (β eff ) to Pu-239 nuclear data were calculated using a newly available tool (ACEtk). Uncertainty calculations showed the following trends: (1) The prompt neutron decay constant can be used to reduce nuclear data-induced uncertainty in Pu-239(n,f) and (2) the effective delayed neutron fraction can be used to effectively reduce nuclear data-induced uncertainty in Pu-239 total fission $\tilde{v}$. Sensitivity/uncertainty analysis can be used to determine optimal experiment/detector setup to identify compensating errors in nuclear data. Upcoming work includes investigating other nuclide-reaction pairs and response sensitivities, such as Δ$\rho$, neutron leakage spectra, and reaction rate measurements to constrain nuclear data of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multi-group Examination of Nickel-Reflected HEU System [Slides]

Researchers noted an unusually large bias for 8-in. nickel reflected HEU sphere in HMF-003 between the 252- group library and the CE library in SCALE 6.2.4. The bias was investigated by reviewing reactions that $k_{eff}$ is sensitive to using TSUNAMI and collecting reaction rate data using tools within SCALE. For this system, bias is primarily due to the elastic scattering in nickel. Researchers compared new multigroup structures in SCALE 6.3. For this system, reduction in CE-to-MG bias seen in the 302-group and further improved in the 1597-group structure. When researchers compared ENDF/B-VIII.0 in SCALE 6.3, library showed improved results with new nickel evaluation. However, the MG-to-CE bias can still be large as large biases can occur in any system. This example highlights the importance of validating results with measured systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Reactive Flow Characteristic Equations

The characteristic equations are derived for the reactive flow PDEs in 1-dimension with 1 irreversible reaction. They differ from the non-reacting case only by source terms proportional to the reaction rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

INDEN Fissile Actinide Issues [Slides]

This presentation addresses the outlook for the International Nuclear Data Evaluation Network (INDEN), specifically regarding general issues that have been solved, including Thermal Neutron Constants (TNC) and Prompt Fission Neutron Spectra (PFNS). Other topics discussed include outstanding issues in fissile evaluations, comments on the RPI quasi-differential experiment, and reaction rate testing of trial evaluations in fast assemblies. The presentation concludes by discussing how INDEN interactions strongly helped to find deficiencies in existing evaluations and highlighted potential solutions to existing challenges, that updated trial evaluations for U-235, U-233, Pu-239 are available for testing, that further work is expected on U-235, Pu-239, and U-233, and that performance was tested on the ICSBEP and/or SINBAD benchmarks, during which significant improvement was demonstrated. The INDEN scheme of international collaboration on nuclear data evaluation is working well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

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

Water Mediated Proton Hopping Furazan Mechanism that Accelerates the Decomposition of TATB

The dominant initial decomposition mechanism of TATB is believed to be a dehydration reaction that forms mono- and di-furazans, although other mechanisms have been reported. In this work, a mono-furazan mechanism is discovered that involves a proton hopping process between TATB and water. A first principles nudge elastic band calculation was performed that shows the mechanism has an energy barrier of 44.1 kcal/mol, 25.4 kcal/mol lower than the monofurazan mechanism without water. The reaction rate of the proton hopping mechanism was found to be 10 7 to 10 12 times (at 300-600 K) higher than the mechanism without water. Metadynamics simulations were also performed that indicate the mechanism is qualitatively more probable at lower temperatures (near 300 K) and higher pressures (near 100 bar), although a quantitative confirmation is needed. Future work should focus on modeling the proton hopping mechanism in the condensed state at elevated temperatures (near 600 K).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗