Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase kinetics”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Decoding Li + /Na + Exchange Route Toward High-Performance Mn-Based Layered Cathodes for Li-Ion Batteries

Li + /Na + exchange has been extensively explored as an effective method to prepare high-performance Mn-based layered cathodes for Li-ion batteries, since the first report in 1996 by P. G. Bruce (Nature, 1996. 381, 499–500). Understanding the detailed structural changes during the ion-exchange process is crucial to implement the synthetic control of high-performance layered Mn-based cathodes, but less studied. Herein, in situ synchrotron X-ray diffraction, density functional theory calculations, and electrochemical tests are combined to conduct the systemic studies into the structural changes during the ion-exchange process of an Mn-only layered cathode O3-type Li 0.67 [Li 0.22 Mn 0.78 ]O 2 (LLMO) from the corresponding counterpart P3-type Na 0.67 [Li 0.22 Mn 0.78 ]O 2 (NLMO). The temperature-resolved observations combined with theoretical calculations reveal that the Li + /Na + exchange is favorable thermodynamically and composited with two tandem topotactic phase transitions: 1) from NLMO to a layered intermediate through ≈60% of Li + /Na + exchange. 2) then to the final layered product LLMO through further Li insertion. Moreover, the intermediate-dominate composite is obtained by slowing down the exchange kinetics below room temperature, showing better electrochemical performance than LLMO obtained by the traditional molten-salt method. Finally, the findings provide guides for the synthetic control of high-performance Mn-based cathodes under mild conditions.

25 ENERGY STORAGE↗
DOE OSTI2023

Prediction of nonlinear specific heat during single crystal HMX phase transition

Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.

42 ENGINEERING↗
DOE OSTI2021

A Kinetic Monte Carlo Approach for Simulating Cascading Transmission Line Failure

In this work, cascading transmission line failures are studied through a dynamical model of the power system operating under fixed conditions. The power grid is modeled as a stochastic dynamical system where first-principles electromechanical dynamics are excited by small Gaussian disturbances in demand and generation around a specified operating point. In this context, a single line failure is interpreted in a large deviation context as a first escape event across a surface in phase space defined by line security constraints. The resulting system of stochastic differential equations admits a transverse decomposition of the drift, which leads to considerable simplification in evaluating the quasipotential (rate function) and, consequently, computation of exit rates. Tractable expressions for the rate of transmission line failure in a restricted network are derived from large deviation theory arguments and validated against numerical simulations. Extensions to realistic settings are considered, and individual line failure models are aggregated into a Markov model of cascading failure inspired by chemical kinetics. Cascades are generated by traversing a graph composed of weighted edges representing transitions to degraded network topologies. Numerical results indicate that the Markov model can produce cascades with qualitative power-law properties similar to those observed in empirical cascades.

Kinetic Monte Carlo, rare events, Cascading Failur↗
DOE OSTI2021

Atomistic Insights of Irreversible Li + Intercalation in MnO 2 Electrode

Tunnel-structured MnO 2 represents open-framed electrode materials for reversible energy storage. Its wide application is limited by its poor cycling stability, whose structural origin is unclear. We tracked the structure evolution of β-MnO 2 upon Li + ion insertion/extraction by combining advanced in situ diagnostic tools at both electrode level (synchrotron X-ray scattering) and single-particle level (transmission electron microscopy). The instability is found to originate from a partially reversible phase transition between β-MnO 2 and orthorhombic LiMnO 2 upon lithiation, causing cycling capacity decay. Moreover, the MnO 2 /LiMnO 2 interface exhibits multiple arrow-headed disordered regions, which severely chop into the host and undermine its structural integrity. Our findings could account for the cycling instability of tunnel-structured materials, based on which future strategies should focus on tuning the charge transport kinetics toward performance enhancement.

25 ENERGY STORAGE↗
DOE OSTI2023

Interpreting inorganic compositional depth profiles to understand the rate-limiting step in vapor phase infiltration processes

Vapor phase infiltration (VPI) is a post-polymerization modification technique that infuses inorganics into polymers to create organic–inorganic hybrid materials with new properties. Much is yet to be understood about the chemical kinetics underlying the VPI process. The aim of this study is to create a greater understanding of the process kinetics that govern the infiltration of trimethyl aluminum (TMA) and TiCl 4 into PMMA to form inorganic-PMMA hybrid materials. To gain insight, this paper initially examines the predicted results for the spatiotemporal concentrations of inorganics computed from a recently posited reaction–diffusion model for VPI. This model provides insight on how the Damköhler number (reaction versus diffusion rates) and non-Fickian diffusional processes (hindering) that result from the material transforming from a polymer to a hybrid can affect the evolution of inorganic concentration depth profiles with time. Subsequently, experimental XPS depth profiles are collected for TMA and TiCl 4 infiltrated PMMA films at 90 °C and 135 °C. The functional behavior of these depth profiles at varying infiltration times are qualitatively compared to various computed predictions and conclusions are drawn about the mechanisms of each of these processes. TMA infiltration into PMMA appears to transition from a diffusion-limited process at low temperatures (90 °C) to a reaction-limited process at high temperatures (135 °C) for the film thicknesses investigated here (200 nm). While TMA appears to fully infiltrate these 200 nm PMMA films within a few hours, TiCl 4 infiltration into PMMA is considerably slower, with full saturation not occurring even after 2 days of precursor exposure. Infiltration at 90 °C is so slow that no clear conclusions about mechanism can be drawn; however, at 135 °C, the TiCl 4 infiltration into PMMA is clearly a reaction-limited process, with TiCl 4 permeating the entire thickness (at low concentrations) within only a few minutes, but inorganic loading continuously increasing in a uniform manner over a course of 2 days. Near-surface deviations from the uniform-loading expected for a reaction-limited process also suggest that diffusional hindering is high for TiCl 4 infiltration into PMMA. In conclusion, these results demonstrate a new, ex situ analysis approach for investigating the rate-limiting process mechanisms for vapor phase infiltration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Charge density waves on a half-filled decorated honeycomb lattice

Tight binding models like the Hubbard Hamiltonian are most often explored in the context of uniform intersite hopping $\textit{t}$. The electron-electron interactions, if sufficiently large compared to this translationally invariant $\textit{t}$, can give rise to ordered magnetic phases and Mott insulator transitions, especially at commensurate filling. The more complex situation of nonuniform $\textit{t}$ has been studied within a number of situations, perhaps most prominently in multiband geometries where there is a natural distinction of hopping between orbitals of different degree of overlap. In this paper we explore related questions arising from the interplay of multiple kinetic energy scales and electron-phonon interactions. Specifically, we use determinant quantum Monte Carlo (DQMC) to solve the half-filled Holstein Hamiltonian on a “decorated honeycomb lattice,” consisting of hexagons with internal hopping $\textit{t}$ coupled together by $\textit{t'}$. This modulation of the hopping introduces a gap in the Dirac spectrum and affects the nature of the topological phases. Here, we determine the range of $\textit{t/t'}$ values which support a charge density wave phase about the Dirac point of uniform hopping $\textit{t = t'}$, as well as the critical transition temperature $T_c$. The QMC simulations are compared with the results of mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Comparing approaches for introducing polycyclic aromatic hydrocarbons to Ge(001) to seed graphene nanoribbon synthesis by CH 4 chemical vapor deposition

Here, we evaluate two approaches for introducing polycyclic aromatic hydrocarbons (PAHs) to a graphene catalyst substrate as part of a two-step chemical vapor deposition (CVD) process for growing graphene nanoribbons (GNRs). In this process, PAHs first form graphene-like seeds on Ge, and then GNRs are subsequently evolved from these PAH-derived seeds via substrate-mediated anisotropic growth kinetics during the CVD of CH4. The first seeding approach sublimes controlled doses of PAH thin films into the CVD chamber, while the second delivers PAHs directly from the vapor phase at set concentrations. Using these two approaches, we measure the dependence of GNR density on PAH dose and compare the experimental results with predictions from a rate model of PAH diffusion, desorption, and clustering. We find that PAHs that more strongly adsorb to the catalyst surface (generally larger PAHs) desirably remain more individualized prior to GNR evolution whereas smaller, more weakly bound PAHs aggregate into larger clusters with various sizes on Ge that are undesirable for synthesizing more monodisperse GNRs. As a result, this work is important because it offers a framework that enables the rational selection and design of seed molecules, advancing anisotropic GNR CVD synthesis.

CVD↗
DOE OSTI2023

Effects of Residual DMSO Adduct on Photonically Cured MAPbI 3 Solar Cells

Defects and impurities in halide perovskite solar cells can negatively impact device performance and long-term stability. In this study, we identify photonically cured methylammonium (MA) lead iodide films contain a residual adduct, MA 2 Pb 3 I 8 (dimethyl sulfoxide) 2 (MA 2 Pb 3 I 8 (DMSO) 2 ), that reduces photocurrent generation in perovskite solar cells (PSCs). This is evidenced by a decrease in the external quantum efficiency (EQE) near 400 nm. Similar EQE reductions were observed in PSCs produced using high-speed processes but have not been thoroughly examined. Through X-ray diffraction patterns and Fourier transform infrared spectroscopy, we establish the photo-inactive MA 2 Pb 3 I 8 (DMSO) 2 as the culprit for the EQE reduction. Combined experimental and simulation results reveal that the MA 2 Pb 3 I 8 (DMSO) 2 phase is located at the hole transport layer/ interface, not on the surface, resulting in lower quantum efficiency and surface photovoltage in the short wavelength region. The residual adduct is kinetically trapped due to the short processing time (20 ms) and crystallization direction but can be removed by an additional photonic pulse. Furthermore, this study highlights the need for careful examination of resulting materials beyond device efficiency when transitioning from laboratory processing to industrial high-speed methods.

14 SOLAR ENERGY↗
DOE OSTI2025

Dissolved Oxygen Redox as the Source of Hydrogen Peroxide and Hydroxyl Radical in Sonicated Emulsive Water Microdroplets

Sonicated emulsive water microdroplets (SEWMs) accelerate and enable a variety of catalyst-free chemical transformations. However, significant unanswered questions remain regarding the chemical intermediates they form and their possible redox origin. In this study, we identified dissolved O 2 as the primary originator of reactive oxygen species (ROS) such as OH• and H 2 O 2. We uncovered the role of dissolved O 2 redox by using a combination of microelectrochemical methods to detect H 2 O 2 , isotopic methods to identify the source of H 2 O 2 , and a combination of electron spin resonance and the DMPO spin trap to detect radicals such as OH•. Notably, we found that H 2 O 2 production is correlated with O 2 content via a reduction pathway enabled by a sufficiently large reducing power that can additionally generate H 2 and even perform Pb electroless deposition on Au and Cu metal substrates. Building on our findings, continuous O 2 bubbling of SEWMs showed accumulation of H 2 O 2 up to ∼88 mM in the aqueous phase within 1 h of sonication, demonstrating the scale-up promise of this method. Distinct to sonochemistry of a single phase, this study advances our understanding of the confluence of redox and chemical reaction mechanisms within SEWMs as a biphasic system. This insight paves the way for improving their reaction kinetics, yield, and selectivity, positioning these attractive redox microreactors as alternatives to traditional electrolyzers.

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

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗
DOE OSTI2019

Effect of Mixed-Solvent Environments on the Selectivity of Acid-Catalyzed Dehydration Reactions

The composition of the liquid phase can alter the rates of individual reaction steps and thus alter the selectivity of acid-catalyzed reactions, but these solvent effects are difficult to anticipate for design purposes. Herein, we report the kinetics and selectivity of Brønsted acid-catalyzed 1,2-propanediol dehydration in pure water and in aqueous mixtures of the polar aprotic cosolvents γ-valerolactone, 1,4-dioxane, tetrahydrofuran, N-methyl-2-pyrrolidone, tetramethylene sulfoxide, and dimethyl sulfoxide at 433 K. We find that the major product of 1,2-propanediol dehydration is propanal in most mixed-solvent environments with selectivities between 1 and 68 mol %. In contrast, 1,2-propanediol dehydration in aqueous mixtures of dimethyl sulfoxide affords acetone as the major product with up to 48% selectivity with minimal propanal formation. We use classical molecular dynamics simulations to probe these solvent effects by computing the difference between the solvation free energies of 1,2-propanediol and propanal in aqueous mixtures of polar aprotic cosolvents and in pure water. We find that the difference in the solvation free energies is correlated with the rates of propanal formation in all mixed-solvent environments, indicating that the solvent-mediated stabilization of the product state relative to the reactant state translates to increased selectivity toward the same product. Similar agreement between simulated solvation free energies and experimental reaction rates/selectivities is established for the acid-catalyzed dehydration of cis- and trans-1,2-cyclohexanediol and 1,3-cyclohexanediol. Finally, analysis of the solvation environment around 1,2-propanediol shows that dimethyl sulfoxide uniquely competes against water to solvate reactive hydroxyl groups, which causes a change in reaction mechanism in this solvent system that leads to the formation of acetone rather than propanal. Here, these results represent a step toward the computationally efficient screening of solvent systems for acid-catalyzed, liquid-phase processes.

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

Intrinsic Kinetics of Polyethylene Terephthalate Pyrolysis via Micropyrolysis and Multivariate Chromatographic Analysis

This study provides an in-depth investigation of the primary decomposition of polyethylene terephthalate (PET) via pyrolysis, employing an experimental-analytic workflow that integrates design of experiments (DoE), micropyrolysis coupled with comprehensive two-dimensional gas chromatography (GC×GC), and multivariate data analysis to verify intrinsic kinetic conditions and elucidate evolving product distributions for mapping key reaction pathways. Peaks that could not be identified using commercial spectral libraries were assigned using Mass Frontier simulations, enabling the identification of divinyl terephthalate, ethyl vinyl terephthalate, and 2-(benzoyloxy)ethyl vinyl terephthalate. A polar×polar (non-orthogonal) column set tailored for the detection of carboxylic acids enhanced the quantification of benzoic acid, 4-vinylbenzoic acid, 4-ethylbenzoic acid, and methylbenzoic acid by up to 6-fold relative to an orthogonal column combination (non-polar×mid-polar). Moreover, pyrolysis variables were systematically evaluated using a Box- Behnken design (BBD), encompassing pyrolysis temperature (500−600 °C), sample weight (50−150 μg), and carrier gas flow rate (100−300 mL min −1 ). Among these, pyrolysis temperature was the only statistically significant factor influencing product yields, ranging from 58.78 to 84.26 wt %. In contrast, neither the sample weight nor the carrier gas flow rate had a significant effect on product yields within the evaluated experimental space. At 600 °C, the major pyrolysis products were benzoic acid (up to 20.20 ± 1.46 wt %) and CO 2 (up to 21.28 ± 1.46 wt %), which can be produced through decarboxylation reactions. These findings underscore the critical importance of selecting appropriate analytical columns for the accurate quantification of heteroatomcontaining products such as carboxylic acids, which may otherwise be underestimated or undetected due to their reactivity with the stationary phase of non-polar and mid-polar columns, as well as other GC components. They also highlight the importance of selecting pyrolysis conditions for investigating the primary decomposition of PET under an isothermal kinetically limited regime.

aromatic compounds↗
DOE OSTI2026

Probing Dynamic Phase Changes in Electrochemical Systems with Cryogenic and in Situ Analytical Electron Microscopy (Final Technical Report)

In this project, we focus on the advancement of our newly developed, unique in situ and cryogenic analytical electron microscopy (AEM) techniques and the development of additional high-end capabilities to quantitatively elucidate the mechanisms of resistive state changes in solid-state electrochemical materials for non-volatile, neuromorphic memory. The objectives of the proposed research are (1) to characterize and elucidate the mechanisms by which electronic state changes occur in various solid-state devices to build guiding principles for electrochemically activated smart materials; (2) to advance our understanding of such mechanisms by coupling cryogenic and in situ techniques to control kinetics of the observed mechanisms while mitigating probe-induced damage; and (3) to utilize additional advanced techniques including electron-beam induced current and quantum mechanical simulations in order to identify the underlying fundamental science of such processes, including roles of defect formation and its effect on phase nucleation, as well as dynamics of these processes in the electrochemical systems at the nanometer scale.

25 ENERGY STORAGE↗
DOE OSTI2022

Mechanistic pathways and role of oxygen in oxidative coupling of methane derived from transient kinetic studies

Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

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

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗