Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

C3MechLite: An integrated component library of compact kinetic mechanisms for low-carbon, carbon neutral and zero-carbon fuels

Based on our latest detailed chemical reaction mechanism, C3MechV4.0, we have developed two reduced reaction mechanisms—C3MechLite and C3MechCore—targeting C 0 –C 3 chemical species including NH 3 . C3MechLite (61 species), contains a number of species comparable to GRI-Mech (53 species), that can accurately predict the combustion characteristics of hydrogen, carbon monoxide, ammonia, methane, natural gas, nitrogen oxides, and their mixtures for a wide range of conditions. C3MechCore (118 species) targets a more comprehensive range of C 0 –C 3 fuels, including ammonia, methanol, ethanol, and dimethyl ether. Both mechanisms demonstrate predictive accuracy comparable to C3MechV4.0 for the combustion characteristics of the target fuels. C3MechLite is designed with a component library structure, enabling further reduction in mechanism size depending on the fuel(s) of interest for 2D/3D numerical simulations. Various combinations of component libraries were validated, and the average prediction error remains within 1 % compared to C3MechLite. Furthermore, the mechanism was applied to 3D LES simulations of H 2 lifted flames and was confirmed to reproduce flame characteristics with high accuracy. C3MechLite and its component library structure enable high-fidelity and computationally efficient chemical kinetic mechanisms, paving the way for application in more complex combustion simulations.

Ammonia↗
DOE OSTI2026

Methyl formate oxidation kinetics up to 100 atm

Methyl formate (MF, CH3OCHO), the simplest ester, is a representative oxygenated fuel with high oxygen content, and low sooting tendency. However, its oxidation behavior under high-pressure and intermediate-temperature conditions remains insufficiently understood, especially where low-temperature peroxy radical chemistry, methanol chemistry, and pressure-dependent reaction pathways play a critical role. In this study, MF oxidation experiments were conducted in the Princeton supercritical-pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm over the temperature range of 400–950 K under both fuel-lean and fuel-rich conditions. Based on the experimental results, an updated HP-Mech was developed by incorporating previous MF sub-mechanisms, expanded low-temperature peroxy pathways, and evaluated pressure-dependent decomposition kinetics. The newly updated HP-Mech shows greatly improved performance in predicting the onset temperature, the key intermediate species fractions, methanol formation, and the progression of MF oxidation across all the experimental conditions. Path flux analysis indicates that MF consumption at the onset stage is dominated by H-abstraction at the methyl site, forming CH2OCHO radicals that lead to the formation and isomerization of O2CH2OCHO, driving low-temperature chain propagation. Moreover, H-abstraction at the formate site forms CH3OCO radicals that preferentially decompose to CH3, initiating the methanol formation pathway linked to CH3O2 and HO2 chemistry. At the same time, HO2 formation is strongly coupled to MF oxidation through multiple MF-derived radical pathways. HCO originates from MF oxidation and acts as a key coupling species linking fuel consumption to HO2 buildup, especially under high-pressure and intermediate-temperature conditions. In addition to this dominant channel, supplementary HO2 formation pathways involving CH3, CH3O, CH2OH, and CH3O2 reacting with O2 further connect methanol chemistry and oxygenated radical chemistry to the HO2 pool, indicating the central role of HO2 in governing MF oxidation. Sensitivity analysis identifies MF with OH/HO2/CH3O2 reactions and the HO2/H2O2/OH sequence as the key factors controlling reactivity in the high-pressure and intermediate-temperature regime. MF directly reacts with OH/HO2/CH3O2 to consume the fuel and produce reactive radicals like CH2OCHO and CH3OCO that undergo subsequent oxidation pathways. Moreover, HO2 recombination suppresses oxidation at lower temperatures, while thermal decomposition of H2O2 accelerates OH production and promotes fuel consumption as temperature increases. The direct formation of active OH from HO2 radicals further completes the mechanism, improving its prediction especially during the oxidation onset stage.

Low-temperature Chemistry↗
DOE OSTI2021

Kinetic, metabolic, and statistical analytics: addressing metabolic transport limitations among organelles and microbial communities

Microbial organisms engage in a variety of metabolic interactions. A crucial part of these interactions is the exchange of molecules between different organelles, cells, and the environment. The main forces mediating this metabolic exchange are transporters. This transport can be difficult to measure experimentally because several transport mechanisms remain opaque. However, theoretical calculations about the inputs and outputs of cells via metabolic exchanges have enabled the successful inference of the workings of intra-organismal and inter-organismal systems. Kinetic, metabolic, and statistical modeling approaches in combination with omics data are enhancing our knowledge and understanding about metabolic exchange and mass resource allocation. Furthermore, this model-driven analytics approach can guide effective experimental design and yield new insights into biological function and control.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Parabolic oxidation kinetics of a plutonium alloy at room temperature

Auger electron spectroscopy (AES) was used to characterize the oxide scale formed on an ~0.5 wt. % Ga plutonium metal alloy. Analysis of changes to the characteristic Auger transitions for Pu and O during sputter depth profiles were used in combination with previously measured sputter yields to determine the thickness of the oxide scale and provide insight into its chemistry. From this data, the parabolic kinetics for the Pu-Ga alloy were determined as a function of relative humidity, RH = 14–92%, and oxygen partial pressure, p(O 2 ) = 31–157 Torr, at 25°C. For each relative humidity investigated, the parabolic rate constant was found to be independent of oxygen partial pressure. Furthermore, power-law fits to the relative humidity effect on the rate constant was found to have a power of n = 1.265, indicating the diffusion species which rate-limits the reaction isn’t molecular water, which subsequent dissociates at the oxide-metal interface, or hydroxyl from dissociation at the gas-surface interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oxidation kinetics and microstructure evolution of high Mn stainless-steel alloy in CO 2 at 700 °C

Oxidation kinetics and scale formation were examined in 21–2N Stainless Steel alloys in a CO 2 environment at 700 °C. Several characterization techniques were used to identify the morphology, crystallographic structure, and chemical composition of the oxide scale formation during exposures up to 1925 h. High manganese content played an important role in homogenous oxide scale formation of Mn 3 O 4 and MnCr 2 O 4 at the scale/gas and scale/alloy interfaces and as an FCC stabilizer in the alloy. Depletion of Mn resulted in an FCC to BCC phase transformation at the scale/alloy interface. No carburization was observed, indicating that the oxide layer is preventing carbon uptake.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revisiting a minimally destructive analytic approach for determining electrochemical kinetic parameters: Measuring aluminum corrosion across a wide pH range based on the Butler-Volmer equation

Here, this study revisits the three-point sampling of the simplified Butler-Volmer equation to address the limitations of strong potentiodynamic polarization, which can introduce irreversible damage and uncertainty in corrosion analysis. The method extracts electrochemical kinetic parameters while minimizing polarization effects, evaluates noise sensitivity relative to overpotential, and accounts for errors from signal noise, OCP drift, ohmic resistance, and mass-transfer constraints. Verified against the Tafel extrapolation method for aluminum corrosion across a wide pH range, this low-polarization approach enables accurate evaluations with specific error estimates, offering a robust alternative to linear polarization resistance methods that assume constant Tafel slopes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗
DOE OSTI2023

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗
DOE OSTI2020

Boosting energy efficiency of Li-rich layered oxide cathodes by tuning oxygen redox kinetics and reversibility

In developing electrode materials for next-generation Li-ion batteries, significant efforts have been given to the energy, power density and cycling stability, with much less (if any) attention paid to the energy efficiency – arguably, the most important practical measure for large-scale applications. This is particularly true for the oxygen-redox active electrodes, such as Li 1.2 Ni 0.13 Co 0.13 Mn 0.54 O 2 , the notorious energy-inefficient cathode that has an extremely high capacity but comes with large voltage hysteresis and voltage decay. Herein, we report the rational design of an energy-efficient Li-rich layered cathode along with high energy, power density and cycling stability enabled by tuning oxygen redox activity. Specifically, the target material Li 1.12 Ni 0.22 Co 0.13 Mn 0.52 O 2 exhibits an ultrahigh energy efficiency at 1 C (90.6%), high capacity (> 200 mAh g -1 ) with 98.9% retention and less than 150 mV decay at the extended 200 cycles. Overall, through direct comparison between the material and Li 1.2 Ni 0.13 Co 0.13 Mn 0.54 O 2 , we show that the compositional change, although slightly, greatly improves the oxygen redox kinetics and reversibility, thereby boosts energy efficiency. The findings offer a strategy to narrow the gap between scientific interest and practical application of oxygen-redox chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Observations of multi-ion physics and kinetic effects in a surrogate to the solar CNO reactions

The ‘CNO process’ occurs in heavier stars with finite metallicity in which hydrogen burning is catalyzed in the presence of 12 C. These reactions are more strongly dependent on temperature than the pp cycle reactions, and thus the CNO cycle dominates only in massive stars. For these types of reactions to be studied at ICF facilities such as OMEGA, an implosion platform using heavier nuclei in the fuel and capable of creating ion temperatures on the order of at least 20 keV is required. A potential route to reach these conditions is to take advantage of kinetic effects in low-convergence shock-driven ‘exploding pusher’ implosions. In this experiment, shots were conducted at the OMEGA laser facility using the surrogate reaction 13 C + D. Its cross section is substantially higher than the actual astrophysical CNO reactions. The yield of this reaction in these implosions was much lower than expected. Finally, physical explanations are discussed, with significant species stratification the likely explanation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effect of particle characteristics, kinetics and transport phenomena on the prediction of particle mass loss and products yields during biomass fast pyrolysis

Here a low computational cost model for the pyrolysis of biomass anisotropic particles was used to study the effect of thermal conductivity of the particle, CRECK kinetic scheme implementation, intraparticle secondary reactions, heats of reaction, and advection on particle conversion prediction. The model considers a shrinking anisotropic cylindrical particle with biomass pseudo-components, the appearance of liquid intermediates, intra-particle mass and heat transport and secondary reactions of volatiles. The model was validated against oak and birchwood single particle experiments. More accurate predictions for mass loss evolution are found when considering the change of the particle thermal conductivity anisotropy with conversion and the temperature dependency of the thermal conductivities of the biomass and the gas phase, compared to ignoring these phenomena. The CRECK biomass pyrolysis reaction scheme provides insights into the effect of biomass composition on the evolution of the liquid intermediate phase allowing future studies about the presence of heavy compounds in the bio-oil and other phenomena such as aerosol ejection. Secondary intraparticle vapor phase reactions showed a negligible effect in volatiles conversion while heterogeneous reactions slightly overestimate its conversion. Heats of reaction and product advection are important for accurate prediction of particle mass loss.

09 BIOMASS FUELS↗
DOE OSTI2020

Kinetics of low-temperature methane activation on IrO 2 (1 1 0): Role of local surface hydroxide species

The ability of the IrO 2 (1 1 0) surface to promote CH 4 activation at low temperatures (~150 K) suggests possibilities for developing IrO 2 -based catalysts to selectively convert light alkanes to high-value chemicals. In this study, we present experimental results and microkinetic simulations, based on density functional theory (DFT) calculations, of temperature programmed reaction spectra (TPRS) obtained for CH 4 σ-complexes adsorbed on IrO 2 (1 1 0), and focus on clarifying how surface OH groups modify the branching between CH 4 desorption, activation and subsequent reactions during TPRS. Our DFT results predict that surface OH groups strongly destabilize CH 4 σ-complexes on IrO 2 (1 1 0) in addition to deactivating surface O-atoms that are needed to achieve CH 4 activation at low temperature. We demonstrate that a microkinetic model that incorporates the influence of surface OH groups on the CH 4 binding and reactivity reproduces experimental TPRS results which show that adsorbed CH 4 σ-complexes preferentially dissociate at low CH 4 coverage, but that an increasing fraction of the adsorbed CH 4 desorbs at low temperature (~120 K) with increasing initial CH 4 or H coverage. The simulations reveal that low-temperature CH4 desorption during TPRS arises primarily from CH 4 σ-complexes that are kinetically-trapped between adjacent surface OH groups, and are thus destabilized and unable to access reactive O-atoms. In addition, when the effect of adjacent OH groups is incorporated we find that the PBE-D3 functional incorporating dispersion provides CH 4 binding energies that agree more closely with that obtained from the TPRS experiments. Overall, our results demonstrate that the local effect of adjacent OH groups must be incorporated into any microkinetic models to properly capture the selectivity between extensive and partial oxidation in alkane conversion on IrO 2 (1 1 0) under reaction conditions.

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

Promoting the water dissociation of nickel sulfide electrocatalyst through introducing cationic vacancies for accelerated hydrogen evolution kinetics in alkaline media

Developing cost-effective non-precious electrocatalysts with excellent HER performance in alkaline media is of significance for industrial hydrogen production. Herein, we successfully endow Ni 3 S 2 with unprecedentedly enhanced alkaline HER activities by introducing nickel vacancies. The Ni 3 S 2 catalyst with nickel vacancies exhibits an ultralow overpotential of 35 mV at 10 mA cm –2 in 1 M KOH electrolyte, outperforming commercial Pt/C (56 mV) and most reported non-precious electrocatalysts. The electron paramagnetic resonance and X-ray absorption fine structure etc. validates the formation of Ni vacancies and the change of the geometric/electronic structure of Ni 3 S 2 . Additionally, the theoretical calculations reveal that the Ni sites with higher valence state neighbouring the Ni vacancies can promote the water dissociation, playing a vital role in accelerating HER kinetics in alkaline media. Finally, this work may provide a universal strategy to enhance the HER activity of transition-metal compound catalysts in alkaline media by means of cationic defect engineering.

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

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

In this paper, we consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗