Search NASA⌕ Search

SEARCH · Search NASA

Results for “Theoretical 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Property of neoclassical GAMs induced by pellet generated plasma perturbations in the gyrokinetic code XGC

A pellet source model has been implemented into the gyrokinetic code XGC and applied to neoclassical geodesic acoustic mode (GAM) study as the first step to investigate potentially rapid kinetic spread of pellet-born particles following their injection into plasma. Here in this study, GAM oscillations of the radial electric field, as well as of their effect on the radial particle and energy fluxes, are studied and the observed frequencies agree with theoretical expectations. It is verified that (i) GAM oscillations are driven both within the pellet-fuelled region and for locations at a larger radius than this region and not for locations radially inside the pellet-fuelled region, and that, (ii) even though the poloidal spread along the magnetic field lines is fast, the time-averaged radial plasma transport under GAM oscillations is kept to roughly the same level as the neoclassical transport in the absence of plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Photoinduced band renormalization effects in the topological nodal-line semimetal ZrSiS

Out-of-equilibrium effects provide an elegant pathway for probing and understanding the underlying physics of topological materials. In particular, controlling electronic band structure properties using ultrafast optical pulses has shown promise for creating exotic states of matter. Of recent interest is band renormalization in Dirac and Weyl semimetals as it leads to direct physical observables through the enhancement of the effective mass or in the shift of resonant energies. Here, in this work, we provide experimental and theoretical signatures of photoinduced renormalization of the electronic band structure in the topological nodal-line semimetal ZrSiS. Specifically, we show how the change in the transient reflectivity spectra under femtosecond optical excitations is induced by out-of-equilibrium effects that renormalize the kinetic energy of electrons. We associate the observed spectral features with an enhancement of the effective mass and to a redshift of the resonant frequency as a function of pump field strength. Finally, we show that the transient relaxation dynamics of the reflectivity is primarily an electronic effect with a negligible phononic contribution. Our study presents conclusive modifications of electronic properties in ZrSiS using ultrashort pulses and demonstrates the potential of this approach in creating photoinduced phases in topological quantum matter through an all-optical route.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Physics-Informed Neural Network (PINN) Prediction of Mixed Mass-Heat-Crystallization Limited Methane Hydrate Formation and Dissociation in Micro-Confinement

The creation and use of Physics-Informed Neural Networks (PINNs) for simulating the dynamics of methane hydrate formation and dissociation will be presented. The PINN framework's main benefit is its capacity to impose physical consistency with only a partial comprehension of the governing equations. This makes the algorithm especially useful for systems with little experimental evidence or a lack of theoretical knowledge. A strong basis for forecasting methane hydrate behavior over the verified operating ranges of 30.0-80.9 bar pressure and 1.0-4.0 K sub-cooling conditions is provided by the combination of conductive heat transfer equations and mixed mass-transfer–crystallization kinetics. PINNs were more accurate at predicting the mixed mass-heat-crystallization limited kinetics than conventional Artificial Neural Networks (ANNs), demonstrating remarkable predictive accuracy for methane hydrate production over the ANN model. The efficiency of incorporating physical limitations from first principles into machine learning frameworks for methane hydrate crystallizations is reinforced by these findings. For hydrate-related applications in energy generation, carbon sequestration, and climate modelling, our study establishes PINNs as a computational tool that is both scalable and efficient. The proven capacity to close the gap between conventional physics-based simulations and solely data-driven models creates new opportunities for expedited hydrate research and practical applications.

Hartman, Ryan L [NYU Tandon School of Engineering]↗
DOE OSTI2022

Sequential Bond Dissociation Energies of Th + (CO) x , x = 3–6: Guided Ion Beam Collision-Induced Dissociation and Quantum Computational Studies

Collision-induced dissociation (CID) of [Th,xC,xO] + , x = 3–6, with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). Products are formed exclusively by the loss of CO ligands. Analyses of the kinetic energy-dependent CID product cross sections yield bond dissociation energies (BDEs) of (CO) x–1 Th + –CO at 0 K as 1.09 ± 0.05, 0.82 ± 0.07, 0.63 ± 0.05, and 0.70 ± 0.05 eV, respectively. Different structures of [Th,xC,xO] + were explored using various electronic structure methods, and BDEs for CO ligand loss from precursor [Th,xC,xO] + complexes were computed. Both experimental and theoretical results corroborate that the structures of [Th,xC,xO] + , x = 3–6, formed experimentally are homoleptic thorium cation carbonyl complexes, Th + (CO) x . The nonmonotonic trend in experimental BDEs is reproduced theoretically, although ambiguities in the spin states of the x = 4–6 complexes (doublet or quartet) remain. BDEs calculated at the coupled cluster with single, double, and perturbative triple excitations (CCSD(T))/cc-pVXZ//B3LYP/cc-PVXZ (X = T and Q) level and a complete basis set (CBS) extrapolation agree reasonably well with the experimental values for all complexes. Thorium oxide ketenylidene carbonyl cations, OTh + CCO(CO) y , y = 1–4, were calculated to be the most stable structures of [Th,xC,xO] + , x = 3–6, respectively; however, these are not observed in our experiment. Potential energy profiles (PEPs) having either quartet or doublet spin calculated at the B3LYP/cc-pVQZ level suggest that the failure to observe OTh + CCO(CO) y , y = 1–4, is the result of a barrier corresponding to the C–C bond formation, making the formation of OTh + CCO(CO) y inaccessible kinetically under the present experimental conditions.

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

Hydrodechlorination of 1,2-Dichloroethane on Platinum Catalysts: Insights from Reaction Kinetics Experiments, Density Functional Theory, and Microkinetic Modeling

Catalytic hydrodechlorination is a promising strategy for treating industrial 1,2-dichloroethane wastes, for which Pt and Pt-based alloy catalysts are widely used. Here, we performed a detailed mechanistic study for 1,2-dichloroethane hydrodechlorination on Pt using a synergistic approach combining density functional theory (DFT) calculations, reaction kinetics experiments, and microkinetic modeling. Using planewave DFT calculations, we evaluated the reaction energy and activation energy barrier of each elementary step involved in the reaction network on Pt(111). The calculated energetics were then incorporated into a comprehensive mean-field microkinetic model accounting for a total of 65 elementary steps. The model-predicted reaction rates were compared with the results from our reaction kinetics experiments on SiO2-supported Pt catalysts. Our results indicated that the hydrodechlorination of 1,2-dichloroethane on Pt(111) starts with a H-removal step; then, it proceeds through a sequence of alternating dechlorination and dehydrogenation steps until vinylidene (CH 2 C*) is formed; finally, CH 2 C* is hydrogenated to the final product, ethane, sequentially via vinyl (CH 2 CH*), ethylene, and ethyl (CH 3 CH 2 *) intermediates. After model parameter adjustments, we achieved good agreement between our theoretical model and experimental results; the adjustments to the calculated parameters are consistent with the typically anticipated coverage effects. Furthermore, our study offers valuable mechanistic insights, which are useful for improving catalysts for this chemistry.

1,2-dichloroethane↗
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗
DOE OSTI2020

Cool flame chemistry of diesel surrogate compounds: n -Decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane

Elucidating the formation of combustion intermediates is crucial to validate reaction pathways, develop reaction mechanisms and examine kinetic modeling predictions. While high-temperature pyrolysis and oxidation intermediates of alkanes have been thoroughly studied, comprehensive analysis of cool flame intermediates from alkane autoxidation is lacking and challenging due to the complexity of intermediate species produced. Here, jet-stirred reactor autoxidation of four C 10 alkanes: n -decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane, as model compounds of diesel fuel, was investigated from 500 to 630 K using synchrotron vacuum ultraviolet photoionization molecular beam mass spectrometry (SVUV-PIMS). Around 100 intermediates were detected for each fuel. The classes of molecular structures present during the autoxidation of the representative paraffinic functional groups in transport fuels, i.e., n -alkanes, branched alkanes, and cycloalkanes were established and were found to be similar from the oxidation of various alkanes. A theoretical approach was applied to estimate the photoionization cross sections of the intermediates with the same carbon skeleton as the reactants, e.g., alkene, alkenyl keto, cyclic ether, dione, keto-hydroperoxide, diketo-hydroperoxide, and keto-dihydroperoxide. These species are indicators of the first, second, and third O 2 addition reactions for the four C 10 hydrocarbons, as well as bimolecular reactions involving keto-hydroperoxides. Chemical kinetic models for the oxidation of these four fuels were examined by comparison against mole fraction of the reactants and final products obtained in additional experiments using gas chromatography analysis, as well as the detailed species pool and mole fractions of aforementioned seven types of intermediates measured by SVUV-PIMS. This works reveals that the models in the literature need to be improved, not only the prediction of the fuel reactivity and final products, but also the reaction network to predict the formation of many previous undetected intermediates.

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

Neural network interatomic potential-driven analysis of phase stability in Ti–V alloys at the atomistic scale

The evolution of the ω phase in titanium–vanadium (Ti–V) alloys is critical for their mechanical properties, particularly in aerospace and biomedical applications. Here, this study employs a Rapid Artificial Neural Network (RANN) potential to model the ω phase evolution at the atomistic level, demonstrating a high degree of consistency with experimental observations, unlike the Modified Embedded Atom Method (MEAM), which fails to capture this phase transformation accurately. RANN simulations replicate key phenomena such as the nucleation of α precipitates at ω/β interfaces and accurate lattice orientations, enhancing our understanding of phase stability and transformation kinetics. The findings affirm that RANN potentials can significantly improve the prediction accuracy of complex material behaviors, offering a powerful tool for designing advanced materials with tailored properties such as solute effect in various stacking fault energies. This approach not only bridges the gap between theoretical predictions and empirical data but also sets a new direction for future research in materials science, emphasizing the integration of machine learning techniques in the development and optimization of new alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Control of Catalyst Isomers Using an N -Phenyl-Substituted RN(CH 2 CH 2 P i Pr 2 ) 2 Pincer Ligand in CO 2 Hydrogenation and Formic Acid Dehydrogenation

A novel pincer ligand, i Pr PN Ph P [PhN- (CH 2 CH 2 P i Pr 2 ) 2 ], which is an analogue of the versatile MACHO ligand, iPr PN H P [HN(CH 2 CH 2 P i Pr 2 ) 2 ], was synthesized and characterized. The ligand was coordinated to ruthenium, and a series of hydride-containing complexes were isolated and characterized by NMR and IR spectroscopies, as well as X-ray diffraction. Comparisons to previously published analogues ligated by iPr PN H P and iPr PN Me P [CH 3 N(CH 2 CH 2 P i Pr 2 ) 2 ] illustrate that there are large changes in the coordination chemistry that occur when the nitrogen substituent of the pincer ligand is altered. For example, ruthenium hydrides supported by the iPr PN Ph P ligand always form the syn isomer (where syn/anti refer to the relative orientation of the group on nitrogen and the hydride ligand on ruthenium), whereas complexes supported by iPr PN H P form the anti isomer and complexes supported by iPr PN Me P form a mixture of syn and anti isomers. We evaluated the impact of the nitrogen substituent of the pincer ligand in catalysis by comparing a series of iPr PN R P (R = H, Me, Ph)-ligated ruthenium hydride complexes as catalysts for formic acid dehydrogenation and carbon dioxide (CO 2 ) hydrogenation to formate. The iPr PN Ph P-ligated species is the most active for formic acid dehydrogenation, and mechanistic studies suggest that this is likely because there are kinetic advantages for catalysts that operate via the syn isomer. In CO 2 hydrogenation, the iPr PN Ph P-ligated species is again the most active under our optimal conditions, and we report some of the highest turnover frequencies for homogeneous catalysts. Experimental and theoretical insights into the turnover-limiting step of catalysis provide a basis for the observed trends in catalytic activity. Additionally, the stability of our complexes enabled us to detect a previously unobserved autocatalytic effect involving the base that is added to drive the reaction. Overall, by modifying the nitrogen substituent on the MACHO ligand, we have developed highly active catalysts for formic acid dehydrogenation and CO 2 hydrogenation and also provided a framework for future catalyst development.

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

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

Here, we analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call “complexity reduction”. In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature T g of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2–6 nm at T g . Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

cancer↗
DOE OSTI2025

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗
DOE OSTI2021

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗
DOE OSTI2024

Elucidating the Discharge Behavior of Aqueous Zinc Sulfur Batteries in the Presence of Molybdenum(IV) Chalcogenide Catalyst: The Criticality of Interfacial Electrochemistry

The aqueous zinc-sulfur battery holds promise for significant capacity and energy density with low cost and safe operation based on environmentally benign materials. However, it suffers from the sluggish kinetics of the conversion reaction. Here, we highlight the efficacy of molybdenum(IV) sulfide (MoS 2 ) to reduce the overpotential of S-ZnS conversion in aqueous electrolytes and study the discharge products formed at the solid-solid and solid-liquid interfaces using experimental and theoretical approaches. Specifically, the MoS 2 -catalyzed electrochemical conversion reaction is characterized via ex situ X-ray diffraction (XRD), transmission electron microscopy (TEM) with energy dispersive spectroscopy (EDS), Raman spectroscopy, synchrotron-based Mo K-edge X-ray absorption spectroscopy (XAS), and in situ synchrotron-based X-ray computed tomography (XCT). Additionally, operando synchrotron-based S K-edge XAS and X-ray fluorescence (XRF) maps are collected to determine the spatial evolution of sulfur-based species at the electrode-electrolyte interface. Further, coupling the operando S K-edge XAS data with the simulated spectra and fitting the data suggested a possible ZnS 2 intermediate phase.

25 ENERGY STORAGE↗
DOE OSTI2023

Stability of perpendicular magnetohydrodynamic shocks in materials with ideal and nonideal equations of state

Magnetized target fusion approach to inertial confinement fusion involves the formation of strong shocks that travel along a magnetized plasma. Shocks, which play a dominant role in thermalizing the upstream kinetic energy generated in the implosion stage, are seldom free from perturbations, and they wrinkle in response to upstream or downstream disturbances. In Z-pinch experiments, significant plasma instability mitigation was observed with pre-embedded axial magnetic fields. To isolate effects, in this work we theoretically study the impact of perpendicular magnetic fields on the planar shock dynamics for different equations of state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when $\gamma$ > 2 or it weakens them when $\gamma$ < 2. Weak shocks have been found to be stable regardless of the magnetic plasma intensity and gas compressibility; however, for sufficiently strong shocks the magnetic fields can promote a neutral stability/SAE at the shock if the adiabatic index is higher than 1 + $\sqrt{2}$. Finally, results have been validated with numerical simulations performed with the FLASH code.

36 MATERIALS SCIENCE↗