Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The mosaic structure of plasma bulk flows in the Earth's magnetotail

Moments of plasma distributions observed in the magnetotail vary with different time scales. In this paper we attempt to explain the observed variability on intermediate timescales of approximately 10-20 min that result from the simultaneous energization and spatial structuring of solar wind plasma in the distant magnetotail. These processes stimulate the formation of a system of spatially disjointed. highly accelerated filaments (beamlets) in the tail. We use the results from large-scale kinetic modeling of magnetotail formation from a plasma mantle source to calculate moments of ion distribution functions throughout the tail. Statistical restrictions related to the limited number of particles in our system naturally reduce the spatial resolution of our results, but we show that our model is valid on intermediate spatial scales Delta(x) x Delta(z) equal to approximately 1 R(sub E) x 1000 km. For these spatial scales the resulting pattern, which resembles a mosaic, appears to be quite variable. The complexity of the pattern is related to the spatial interference between beamlets accelerated at various locations within the distant tail which mirror in the strong near-Earth magnetic field. Global motion of the magnetotail results in the displacement of spacecraft with respect to this mosaic pattern and can produce variations in all of the moments (especially the x-component of the bulk velocity) on intermediate timescales. The results obtained enable us to view the magnetotail plasma as consisting of two different populations: a tailward-Earthward system of highly accelerated beamlets interfering with each other, and an energized quasithermal population which gradually builds as the Earth is approached. In the near-Earth tail, these populations merge into a hot quasi-isotropic ion population typical of the near-Earth plasma sheet. The transformation of plasma sheet boundary layer (PSBL) beam energy into central plasma sheet (CPS) quasi-thermal energy occurs in the absence of collisions or noise. This paper also clarifies the relationship between the global scale where an MHD description might be appropriate and the lower intermediate scales where MHD fails and large-scale kinetic theory should be used.

Ashour-Abdalla, M.↗
DOE OSTI2021

Efficient and reversible electron bifurcation with either normal or inverted potentials at the bifurcating cofactor

A longstanding mystery surrounding electron bifurcation is the significance of inverted (or “crossed”) reduction potentials of the two-electron bifurcating cofactor. Using a many-electron open-system kinetic model, we show that reversible and efficient electron bifurcation is possible without inverted reduction potentials at the bifurcating site if the absolute value of the difference between first and second reduction potentials of the bifurcating species is sufficiently large (on the scale of the redox-potential span of the high- and low-potential branches). Surprisingly, the case with strong, normally ordered potentials at the bifurcating cofactor can produce electron bifurcation that is just as effective as the case with strongly inverted potentials. Lastly, this finding amplifies the puzzle surrounding the recruitment of inverted potentials in the few well-characterized bifurcating systems of nature and suggests that electron bifurcating cofactors without strongly inverted potentials may yet be discovered.

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

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
DOE OSTI2025

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2016

Kinetics of Levoglucosenone Isomerization

In this paper, we studied the acid-catalyzed isomerization of levoglucosenone (LGO) to 5-hydroxymethylfurfural (HMF) and developed a reaction kinetics model that describes the experimental data across a range of conditions (100–150 °C, 50–100 mm H 2 SO 4 , 50–150 mm LGO). LGO and its hydrated derivative exist in equilibrium under these reaction conditions. Thermal and catalytic degradation of HMF are the major sources of carbon loss. Within the range of conditions studied, higher temperatures and shorter reaction times favor the production of HMF. The yields of HMF and levulinic acid decrease monotonically as tetrahydrofuran is added to the aqueous solvent system, indicating that water plays a role in the LGO isomerization reaction. Initial-rate analyses show that HMF is produced solely from LGO rather than from the hydrated derivative of LGO. Finally, the results of this study are consistent with a mechanism for LGO isomerization that proceeds through hydration of the anhydro bridge, followed by ring rearrangement analogous to the isomerization of glucose to fructose.

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

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗
DOE OSTI2023

A projection-based analytical Jacobian framework for chemical kinetics applications

A major challenge in simulating complex combustion systems with detailed chemical kinetic models is the cost of integrating the chemical source terms, often done using stiff ODE solvers that require frequent Jacobian evaluations. Using analytically derived Jacobian matrices instead of divided-difference-based numerical Jacobian approximations can significantly reduce the associated computational cost. However, ambiguities arise in the formulation of analytical Jacobians because the chemical state of the system, or state vector, can be expressed in multiple ways, involving variables that are typically not independent from one another. Here, in this work, the consequences of those ambiguities on practical calculations are characterized in detail, and a generalized, projection-based framework is proposed as a mitigation strategy. Performances are assessed in a series of test cases involving a variety of configurations and numerical solution approaches. Results show that with proper treatment, commonly used analytical Jacobian formulations can be considered as equivalent for practical purposes, thereby alleviating concerns that the state vector chosen to express the governing equations and corresponding analytical Jacobian may significantly impact the accuracy of the simulations.

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

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

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

Prediction and Measurement of the Kinetics of Isothermal Crystallization and Devitrification of Mold Powder Slags

The rate and extent of formation of crystalline phases in the slag film that develops between the mold and the strand in steel continuous casting determine heat transfer and lubrication between them and therefore play crucial roles in shell integrity. In this work, a commercial mold powder for casting medium carbon (MC) thin steel slabs is selected for carrying out detailed comparisons between predicted and experimental results of the rates of isothermal crystallization and devitrification of liquid and glassy slag, respectively. Predictions are obtained with a kinetic model that combines the induction period (of Šimon and Kolman) and the Kissinger methods to evaluate time–temperature–transformation (TTT) diagrams for the devitrification and crystallization reactions associated with the peaks appearing in non-isothermal thermo-analysis traces. To obtain experimental validation of the diagram, a high-temperature confocal laser scanning microscope (HT-CLSM) is used to rapidly heat or cool a sample, ~ 140 mg, of glassy powder or liquid slag, respectively, to selected isothermal conditions, where it is held for a prescribed length of time before quenching. Identification and quantification of the phases present in the quenched samples are done by a quantitative X-ray powder diffraction analysis (QXRPDA) method that uses a PONKCS (partial or not known crystal structure) procedure based on Le Bail and Rietveld strategies for refinement of diffraction patterns. The many experimental results indicate that the predicted TTT diagrams for devitrification and crystallization of the mold slag give a good estimation of the rate of the transformations that take place at different temperatures and reveal the differences in the reactions that occur along both paths.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Methanol Formation via Oxygen Insertion Chemistry in Ices

We present experimental constraints on the insertion of oxygen atoms into methane to form methanol in astrophysical ice analogs. In gas-phase and theoretical studies this process has previously been demonstrated to have a very low or nonexistent energy barrier, but the energetics and mechanisms have not yet been characterized in the solid state. We use a deuterium UV lamp filtered by a sapphire window to selectively dissociate O{sub 2} within a mixture of O{sub 2}:CH{sub 4} and observe efficient production of CH{sub 3}OH via O({sup 1}D) insertion. CH{sub 3}OH growth curves are fit with a kinetic model, and we observe no temperature dependence of the reaction rate constant at temperatures below the oxygen desorption temperature of 25 K. Through an analysis of side products we determine the branching ratio of ice-phase oxygen insertion into CH{sub 4}: ∼65% of insertions lead to CH{sub 3}OH, with the remainder leading instead to H{sub 2}CO formation. There is no evidence for CH{sub 3} or OH radical formation, indicating that the fragmentation is not an important channel and that insertions typically lead to increased chemical complexity. CH{sub 3}OH formation from O{sub 2} and CH{sub 4} diluted in a CO-dominated ice similarly shows no temperature dependence, consistent with expectations that insertion proceeds with a small or nonexistent barrier. Oxygen insertion chemistry in ices should therefore be efficient under low-temperature ISM-like conditions and could provide an important channel to complex organic molecule formation on grain surfaces in cold interstellar regions such as cloud cores and protoplanetary disk midplanes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Influence of gasoline fuel formulation on lean autoignition in a mixed-mode-combustion (deflagration/autoignition) engine

We report stoichiometric spark-ignition engines suffer efficiency penalties due to throttling losses at low loads, a low specific-heat ratio of the stoichiometric working fluid, and limits on compression ratio due to end-gas autoignition leading to undesirable knocking. Mixed-Mode Combustion (MMC) mitigates these shortcomings by using a lean working fluid where a spark-initiated pilot-stabilized deflagrative flame front is followed by controlled end-gas autoignition. This MMC study investigates the effects of initial conditions (intake air temperature, intake pressure, equivalence ratio, and intake oxygen fraction) on autoignition tendency of four gasoline-range fuels with varying properties and composition. The use of fuels with varying octane sensitivity (S) allowed exploring the importance of low-temperature heat release in triggering autoignition. Fuels with high S were less reactive for conditions that promote low-temperature chemistry (operation at high intake air pressure or without N2 dilution). Conversely, an Alkylate fuel with low S showed a greater autoignition resistance at operating conditions that were unfavorable for low-temperature chemistry. Next, the effect of residual gas composition on autoignition tendency of fuels was examined with a chemical-kinetics model. Among the various molecules in the residual gas, nitric oxide (NO) enhanced the low-temperature chemistry and increased the autoignition tendency most significantly. The fuels’ autoignition response to increasing NO amount corroborates the experimental observations. Next, the sequential autoignition of the end-gas was assessed to be less impacted by thermal stratification because of lean mixtures showing relatively less low-temperature chemistry, when compared to stoichiometric mixtures. Next, the effect of changing equivalence ratio on the autoignition was found to be similar for all fuels, regardless of their S. With changing intake air temperature, the response of fuels’ autoignition tendency depended on the dilution level used. At high dilution (i.e. low intake [O2]), fuels’ reactivity increased with increasing intake air temperature. In contrast, for operation without dilution, the autoignition tendency of the low-S Alkylate fuel decreased with increasing intake air temperature, while that of high-S High Cycloalkane fuel still increased with increasing intake air temperature. In conclusion, conventional octane metrics (RON and MON) have utility in assessing the autoignition tendency under lean MMC operation. Moreover, the fuel requirements for MMC align with that of stoichiometric operation: i.e., high RON and high S fuels are desirable for stable non-knocking operation.

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

Narrowing transient testing pulse widths to enhance LWR RIA experiment design in the TREAT facility

The Transient Reactor Test (TREAT) Facility provided thousands of transient irradiations and plays a crucial role in nuclear-heated safety research. The facility’s flexible design and multi-mission nature saw historic experiments for numerous reactor fuels and transient types but was never specifically adapted to address very-brief pulse transients akin to postulated Light Water Reactor (LWR) Reactivity Initiated Accidents (RIA). Since the behaviors of fuel under these conditions depends strongly on energy input duration and resulting cladding time/temperature response under pellet cladding mechanical interactions, this three-year project was conceived to investigate new pulse-narrowing capabilities. Kinetic models showed incremental improvements for minor facility enhancements including increased reactivity step insertions (to initiate the power pulse) and slightly-increased transient rod drive speed for pulse termination (“clipping”). Replacing peak fuel assemblies, repositioning non-transient control rods to hold down the “hot side” of the core, and balancing against required excess reactivity needs, the limiting fuel assembly power can be reduced by ~20 %. Furthermore, this is a remarkable discovery of latent capability in TREAT, not only in enabling the subject capabilities for reduced pulse widths, but also significantly “uprating” the core’s transient energy capacity. Incremental improvements can likely enhance TREAT’s capability into BWR-relevant missions; briefer PWR pulse shapes can only be achieved with an advanced clipping system. Helium-3 ( 3 He) was found to offer the greatest benefits for clipping design and overall performance. A unique concept showed great promise for enabling a 3 He-based system with a reasonable cost; this design concept will likely become the focal point of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

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

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

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

Switching kinetics in ferroelectric zinc magnesium oxide thin films

The switching kinetics of RF-magnetron reactively sputtered ~200 nm thick Zn 1-x Mg x O (ZMO) ferroelectric thin films with x = 0.41 and x = 0.27 prepared on Pt/Ti/SiO 2 /Si substrates were studied at applied fields near the coercive fields, ranging from 3.4 to 5.1 MV cm -1 , and at temperatures ranging from room temperature to 100 °C. Polarization reversal in ZMO followss the Kolmogorov-Avrami-Ishibashi kinetics model for nucleation and growth at applied fields near the coercive field, and obeys an individual column switching (ICS) model at higher applied fields required to switch most of the spontaneous polarization. Switching in the high-applied field regime can be described using either the Gaussian or inverse gamma distribution functions depending on the mole fraction of magnesium in the film. The switching current transients of the high-Mg content ZMO film (x = 0.41) are always bi-modal, whereas the low Mg composition film (x = 0.27) is described by the Gaussian distribution initially following wake-up, and becomes bimodal with continued cycling. Rayleigh-like behavior of the dielectric constant revealed a simultaneous increase of the irreversible and decrease of the reversible contributions to the dielectric constant, which was ascribed to an increase in the density of mobile domain walls with cycling and resulted in faster switching.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1984

Numerical modeling of D-mappings with applications to chemical kinetics

Numerical modeling of D-mappings was studied and applied to solving nonlinear stiff systems. These mappings were locally linearized for convergence analysis, and some applications were made to chemical kinetics. The technique avoids using multistep implicit codes that require inversion of Jacobian matrices, but depends on the Jacobians for its convergence analysis.

Dey, S. K.↗
DOE OSTI2024

Precision Calcination Mechanism of CaCO 3 to High-Porosity Nanoscale CaO CO 2 Sorbent Revealed by Direct In Situ Observations

Deploying energy storage and carbon capture at scale is hindered by the substantial endothermic penalty of decomposing CaCO 3 to CaO and CO 2 , and the rapid loss of CO 2 absorption capacity by CaO sorbent particles due to sintering at the high requisite decomposition temperatures. The decomposition reaction mechanism underlying sorbent deactivation remains unclear at the atomic level and nanoscale due to past reliance on postmortem characterization methods with insufficient spatial and temporal resolution. Thus, elucidating the important CaCO 3 decomposition reaction pathway requires direct observation by time-resolved (sub-)nanoscale methods. Here, chemical and structural dynamics during the decomposition of CaCO 3 nanoparticles to nanoporous CaO particles comprising high-surface-area CaO nanocrystallites are examined. Comparing in situ transmission electron microscopy (TEM) and synchrotron X-ray diffraction experiments gives key insights into the dynamics of nanoparticle calcination, involving anisotropic CaCO 3 thermal distortion before conversion to thermally dilated energetically stable CaO crystallites. Time-resolved TEM uncovered a novel CaO formation mechanism involving heterogeneous nucleation at extended CaCO 3 defects followed by sweeping reaction front motion across the initial CaCO 3 particle. These observations clarify longstanding, yet incomplete, reaction mechanisms and kinetic models lacking accurate information about (sub-)nanoscale dynamics, while also demonstrating calcination of CaCO 3 without sintering through rapid heating and precise temperature control.

36 MATERIALS SCIENCE↗