Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical reaction dynamics”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Single Atom Catalysts: A Review of Characterization Methods

Abstract Single atom catalysts (SACs) harbor a potential to exceed nanoparticle catalysts in terms of activity, stability and selectivity in a growing number of chemical reactions. Although their investigation is attracting significant attention, important fundamental questions focusing on key physicochemical properties of SACs (e. g., structure – property relationships, structural dynamics, reaction‐driven restructuring) remain unanswered. A main challenge for research in the field is how to reliably characterize the environments of single atoms in the presence of complicating factors such as low weight loadings, strong metal‐support interactions, and atomic and multiscale heterogeneity of bonding in the single atom sites. This review addresses this challenge – identifying catalytically relevant features of physicochemical properties of single atoms (charge state, electronic structure, atomic configuration, bonding interactions with a support) and surveying advanced tools/methods for characterizing them. The review places a strong emphasis on multimodal methods exploiting X‐ray absorption, emission and photoelectron spectroscopies, and provides several examples from the authors’ research that demonstrate their use as powerful tools for SAC characterization.

Kottwitz, Matthew↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS2003

Quantum Dynamics Scattering Study of AB+CDE Reactions: A Seven Dimensional Treatment for the H2+C2H Reaction

A time-dependent wave-packet approach is presented for the quantum dynamics study of the AB+CDE reaction system for zero total angular momentum. A seven-degree-of-freedom calculation is employed to study the chemical reaction of H2+C2H yields H + C2H2 by treating C2H as a linear molecule. Initial state selected reaction probabilities are presented for various initial ro-vibrational states. This study shows that vibrational excitation of H2 enhances the reaction probability, whereas the excitation of C2H has only a small effect on the reactivity. An integral cross section is also reported for the initial ground states of H2 and C2H. The theoretical and experimental results agree with each other very well when the calculated seven dimensional results are adjusted to account for the lower transition state barrier heights found in recent ab initio calculations.

Wang, Dunyou↗
DOE OSTI2022

Vibrational control of the reaction pathway in the H + CHD 3 → H 2 + CD 3 reaction

An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant’s CD 3 umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state–specific level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS2004

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Credible Criticality Safety Margin in the 30B Package with LEU+ UF 6 and Hypothetical Water Ingress

The commercial nuclear industry is pursuing advancements in fuel and reactor design that increase the uranium enrichment above 5 wt. % 235 U. These advancements will necessitate the ability to transport bulk quantities of UF 6 at increased enrichments. Currently, the 30B cylinder is the primary container used by the industry for UF 6 storage and transportation and has a long history of successful shipments. This container can support up to 2,277 kg of UF 6 at a maximum enrichment of 5 wt. % 235 U. Previous evaluations have assessed the potential impact of criticality safety for 30B transport at higher enrichments but assumed moderator intrusion would not require evaluation. Although current regulations allow for the exception of moderator intrusion for UF 6 packages through the design and quality control of the package content, this exception is limited to enrichments up to 5 wt. % 235 U. Thus, an investigation of moderator intrusion into a 30B cylinder should be performed. Moderator intrusion into a 30B cylinder is a unique condition for criticality safety evaluation in transportation because of the violent chemical reactions that occur between UF 6 and H 2 O. The resulting intrusion is strongly dependent on the breach size, breach location, breach interface solid/ullage of the UF 6 content, the UF 6 distribution (which is temperature dependent), and temperature/pressure conditions which are dynamic in accident conditions. Additionally, the complexity of the HF-UO 2 F 2 -H 2 O interface that occurs during the event can influence the potential solubility of uranium in the system, as well as influence the amount of UF 6 reaction with H 2 O. With the operating experience from Orano Federal Services LLC, the UF 6 chemical expertise from Oak Ridge National Laboratory (ORNL), and the criticality safety expertise from ORNL, this paper evaluates the neutronic conditions (i.e., k eff ) that apply the understood chemistry and experimental conditions that occur during moderator intrusion of a breached 30B cylinder under postulated accidents. This report examines the historical evaluations of UF 6 transport and expands these evaluations for the enrichments expected for nuclear industry advancement. These simulations primarily feature a homogeneous mixture of UF 6 and H 2 O as an infinite media system and as a sphere with water reflection given an impurity limit of 0.5 wt. % UF 6 . This report demonstrates that this H/U limit is valid up to 8 wt. % enrichment for both H 2 O and HF as the moderating mixture. The water-reflected homogeneous mixed spheres evaluated demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. Additionally, this report evaluates moderator ingress scenarios of a 30B cylinder. The 30B cylinder simulations vary the cylinder orientation, the mass of UF 6 in the system, the mass of H 2 O in the system, and how much H 2 O has reacted with UF 6 . The moderator ingress of a 30B cylinder was evaluated for the water-reflected homogeneous mixed spheres to demonstrate the amount of safety margin applied by restricting the 30B cylinder impurity limit to 0.5 wt. %. These simulations are standard practice and are independent of the UF 6 -H 2 O reaction. This work also explored more complex simulations that layer the H 2 O over the UF 6 , which incorporates the amount of reacted UF 6 as a separate layer and considers solubility limits of UO 2 F 2 in H 2 O as it is produced. This approach is intended to simulate the postulated event of a large hairline crack occurring on a submerged cylinder. Finally, this work considers the potential of H 2 O mixing into the UF 6 solid heterogeneously using a sponge-like model. The water ingress into the UF 6 is treated as random size spheres, and the reaction products form into a layer over the spheres.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Plasmon Dynamics Driven by Aggregation of Tris(2,2′-bipyridine)ruthenium(II)-Functionalized Gold Nanoparticles Probed by XANES and Transient Absorption Spectroscopy

Energy conversion dynamics is critical for advancing next-generation photovoltaics, optoelectronics, and light-harvesting technologies. Noble metal plasmonic nanoparticles play a pivotal role as nanoscale electromagnetic confinement structures, driving photon-induced chemical reactions. In this study, we explore the effects of [Ru(bpy) 3 ] 2+ functionalization and aggregation on citrate-capped gold nanoparticles (AuNPs) of 40 and 100 nm diameters, focusing on molecule-plasmon interactions and their influence on electronic and energy dissipation properties. X-ray absorption near-edge spectroscopy (XANES) revealed that [Ru(bpy) 3 ] 2+ functionalization induces controlled aggregation without altering the oxidation state of gold. A more pronounced white-line intensity is observed in 40 nm AuNPs, consistent with greater s–p–d hybridization and a higher density of surface states, likely influenced by both nanoparticle size and aggregation. Transient absorption (TA) spectroscopy highlights faster electron–phonon relaxation dynamics in aggregated 40 nm nanoparticles, which is attributed to increased electron delocalization and more efficient coupling to the phonon bath. In contrast, 100 nm nanoparticles exhibit minimal changes due to a lower degree of aggregation. Interestingly, we observe that enhanced electron–phonon coupling in aggregated nanoparticles coincides with a slowing of electron–electron scattering. These observations suggest a competitive interplay between the two relaxation pathways, where enhanced energy transfer to the lattice in aggregated systems can suppress electronic thermalization. In conclusion, these findings underscore the critical role of nanoparticle size, aggregation, and molecule–surface interactions in modulating plasmonic dynamics and excited-state lifetimes and further provide valuable insights into designing tailored plasmonic systems with transformative potential for sensing, catalysis, and energy conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimizing the Solvent Reorganization Free Energy by Metal Substitution for Nanocage Catalysis

The supramolecular capsule Ga 4 L 6 12– has been found to catalyze a number of important chemical reactions, but the close proximity of a poorly organized external solvent environment can serve as an impediment to the reaction chemistry. Here, using ab initio molecular dynamic calculations and electric field analyses, we find that metal substitution of indium for gallium lowers the total and electrostatic activation barriers by ~3–4 kcal/mol in water solvent. Using an energy decomposition analysis of the interaction of water solvent with the metal vertices, we find that Pauli repulsion between the metal and water decreases, allowing the water coordination with the metal to increase, upon substitution of In for Ga. We therefore determine that the stabilization of the transition state is due to a better arrangement of water molecules around the metal vertices of In 4 L 6 12– and provides a proof of concept of how to redesign the nanocage scaffold in order to reduce the solvent reorganization energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Liquid rocket combustor computer code development

The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.

Liang, P. Y.↗
NASA NTRS2002

Multi-Zone Liquid Thrust Chamber Performance Code with Domain Decomposition for Parallel Processing

Computational Fluid Dynamics (CFD) has considerably evolved in the last decade. There are many computer programs that can perform computations on viscous internal or external flows with chemical reactions. CFD has become a commonly used tool in the design and analysis of gas turbines, ramjet combustors, turbo-machinery, inlet ducts, rocket engines, jet interaction, missile, and ramjet nozzles. One of the problems of interest to NASA has always been the performance prediction for rocket and air-breathing engines. Due to the complexity of flow in these engines it is necessary to resolve the flowfield into a fine mesh to capture quantities like turbulence and heat transfer. However, calculation on a high-resolution grid is associated with a prohibitively increasing computational time that can downgrade the value of the CFD for practical engineering calculations. The Liquid Thrust Chamber Performance (LTCP) code was developed for NASA/MSFC (Marshall Space Flight Center) to perform liquid rocket engine performance calculations. This code is a 2D/axisymmetric full Navier-Stokes (NS) solver with fully coupled finite rate chemistry and Eulerian treatment of liquid fuel and/or oxidizer droplets. One of the advantages of this code has been the resemblance of its input file to the JANNAF (Joint Army Navy NASA Air Force Interagency Propulsion Committee) standard TDK code, and its automatic grid generation for JANNAF defined combustion chamber wall geometry. These options minimize the learning effort for TDK users, and make the code a good candidate for performing engineering calculations. Although the LTCP code was developed for liquid rocket engines, it is a general-purpose code and has been used for solving many engineering problems. However, the single zone formulation of the LTCP has limited the code to be applicable to problems with complex geometry. Furthermore, the computational time becomes prohibitively large for high-resolution problems with chemistry, two-equation turbulence model, and two-phase flow. To overcome these limitations, the LTCP code is rewritten to include the multi-zone capability with domain decomposition that makes it suitable for parallel processing, i.e., enabling the code to run every zone or sub-domain on a separate processor. This can reduce the run time by a factor of 6 to 8, depending on the problem.

Homayun K. Navaz↗
DOE OSTI2025

Time-resolved chemically-selective spectroscopic investigation of the redox reaction between hematite and aluminium

Thermite reactions –highly energetic redox processes between a metal and an oxide—are used in welding, propulsion, and the fabrication of advanced materials. When reduced to the nanoscale, these reactions exhibit enhanced energetic performance, but their ultrafast dynamics remain poorly understood. Gaining insight into charge transfer during these processes is essential for advancing applications in energy conversion and materials design. Here we show that the reaction between aluminium and hematite, a common iron oxide, can be tracked with femtosecond resolution using extreme ultraviolet (EUV) time-resolved absorption spectroscopy at the Fe M 2,3 and Al L 2,3 edges. By exciting the system with an ultrashort optical pulse and probing element-specific absorption changes, we observe an early spectral shift that reveals the formation of localized charge carriers (polarons). Comparing samples with different supporting substrates highlights ultrafast electron transfer from aluminium to hematite. These results demonstrate an approach to investigating charge flow in energetic materials and provide a basis for studying fast chemical reactions with chemical specificity.

Electron transfer↗
DOE OSTI2022

Mechanistic study of chemical looping reactions between solid carbon fuels and CuO

Copper (Cu) based chemical-looping combustion (CLC) is a promising process that utilizes solid carbon fuel such as coal and biomass. Understanding the reaction kinetics in this process can facilitate the industrial design of CLC units. Here, in this work, molecular dynamics (MD) simulations were performed to investigate the reaction kinetics of n-butane and two simplified solid carbon fuels (lignite and anthracite) with/without copper oxide (CuO) nanoparticle. In addition, experiments were conducted on the thermal characteristics, flammability, and flame speeds of CuO and solid carbon mixtures. For n-butane, oxidation on the CuO surface is the primary reaction since the activation energy of the surface reaction is much lower than that of oxygen (O 2 ) combustion (9.2 vs. 53.3 kcal/mol). On the other hand, for lignite, there is a smaller difference in the activation energy of O 2 combustion and surface reaction (23.04 vs. 6.34 kcal/mol). We hypothesize that some solid carbon fuels can have different reaction kinetics dependent upon temperature. This is proven by the case of a simplified anthracite coal: an Arrhenius plot shows that this solid carbon fuel has two different reaction kinetic regimes and the critical temperature for the change in kinetics is related to the oxygen uncoupling of the CuO nanoparticle. Like the modeling simulations, a change in activation energy is observed in the experimental results, where desorption of molecular oxygen from CuO becomes important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

The rate of the reaction between CN and C2H2 at interstellar temperatures

The rate coefficient for the important interstellar reaction between CN and C2H2 has been calculated as a function of temperature between 10 and 300 K. The potential surface for this reaction has been determined through ab initio quantum chemical techniques; the potential exhibits no barrier in the entrance channel but does show a small exit channel barrier, which lies below the energy of reactants. Phase-space calculations for the reaction dynamics, which take the exit channel barrier into account, show the same unusual temperature dependence as determined by experiment, in which the rate coefficient at first increases as the temperature is reduced below room temperature and then starts to decrease as the temperature drops below 50-100 K. The agreement between theory and experiment provides strong confirmation that the reaction occurs appreciably at cool interstellar temperatures.

NASA Discipline Exobiology↗
DOE OSTI2024

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗