Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 OSTI2021

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Koopman wavefunctions and Clebsch variables in Vlasov–Maxwell kinetic theory

Motivated by recent discussions on the possible role of quantum computation in plasma simulations, here, we present different approaches to Koopman's Hilbert-space formulation of classical mechanics in the context of Vlasov–Maxwell kinetic theory. The celebrated Koopman–von Neumann construction is provided with two different Hamiltonian structures: one is canonical and recovers the usual Clebsch representation of the Vlasov density, the other is non-canonical and appears to overcome certain issues emerging in the canonical formalism. Furthermore, the canonical structure is restored for a variant of the Koopman–von Neumann construction that carries a different phase dynamics. Going back to van Hove's prequantum theory, the corresponding Koopman–van Hove equation provides an alternative Clebsch representation which is then coupled to the electromagnetic fields. Finally, the role of gauge transformations in the new context is discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A kinetic model for pure electron plasma compression via the rotating wall technique

An approximate model for pure electron plasma compression is developed for the case where the rotating wall (RW) electric field couples to the ${E}\times{B}$rotation and axial bounce motion of the electrons. The key assumption in the model is that, throughout the compression, the plasma remains in a slowly evolving thermal equilibrium defined by the plasma temperature and angular momentum. Linearized drift kinetic theory is employed to derive an expression for torque exerted by the RW field on the plasma through coupling to the resonant plasma particles, and averaging is used to find the torque that both compresses and heats the plasma. The evolution equations for the angular velocity and temperature of the plasma include the compression and heating from the torque and cooling from cyclotron radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2025

Isotopic Fractionation and Kinetic Isotope Effects of a Purified Bacterial Nitric Oxide Reductase (NOR)

Nitrous oxide (N 2 O) is a serious concern due to its role in global warming and ozone destruction. Agricultural practices account for ~80% of all anthropogenic N 2 O produced in the US, due in large part to the stimulation of microbial denitrification. Stable isotopes are uniquely suited to examine both microbial N 2 O sources and the mechanism of N 2 O biosynthesis through the use of 2 Site Preference (δ 15 N SP ; the difference in δ 15 N between the central and outer N atoms in N 2 O) and kinetic isotope effects (KIEs), respectively. Using trace gas isotope ratio mass spectrometry (TG-IRMS), we determined the δ 15 N, δ 15 N α , δ 15 N β , and δ 18 O of N 2 O produced by a purified cytochrome c nitric oxide reductase (cNOR) from Paracoccus denitrificans. We also calculated δ 15 N SP , the KIEs, and associated isotopic enrichment factors (ε) for N bulk , N α , and N β . A normal isotope effect was observed for bulk 15 N, with a KIE value of 1.0086 ± 0.0009 (ε = -8.6 ± 0.9‰). The isotope effects for both 15 N α and 15 N β were also normal, with position-specific KIEs of 1.0072 ± 0.0010 (ε = -7.2 ± 1.0‰) and 1.0100 ± 0.0010 (ε = -9.9 ± 1.0‰), respectively, and δ 15 N SP values ranged from 0.5 to 8.7‰ with no significant trend as the reaction proceeded. Values of δ 18 O increased with N 2 O production (slope of δ 18 O against [-flnf/(1-f)] = -19.9 ± 1.9‰). We present implications for the mechanism of N 2 O production from cNOR based on our data.

Rivett, Elise D. [Michigan State Univ., East Lansi↗
DOE OSTI2020

Controlling Ion Coordination Structure and Diffusion Kinetics for Optimized Electrode-Electrolyte Interphases and High Performance Si Anodes

The cycling performance of batteries is largely determined by electrode-electrolyte interphase associated with the chemical and electrochemical properties of electrolyte salts and solvents. In this work, we used a localized high concentration electrolyte (LHCE) that enables high performance of Si anodes as a model system and investigated the scientific mechanism behind the optimization of electrolyte. The correlation between electrode (cathode/anode)-electrolyte interphase and ion coordination structure and diffusion kinetics were elucidated systematically. Among all the LHCEs tested, the one with 1.8M LiFSI concentration and 1:2 molar ratio of carbonate solvents and fluorinated diluents has the highest proportion of solvent-separated ion pairs and contact ion pairs and fastest ion diffusion. It enables the most stable electrode-electrolyte interphase and hence the best Si anode performance. These results reveal the fundamental mechanism behind the stability of LHCEs and provide guidance on the new electrolyte design for Si anodes and other battery systems.

Jia, Haiping↗
DOE OSTI2023

Thermodynamic and Kinetic Barriers Limiting Solid-State Reactions Resolved through In Situ Synchrotron Studies of Lithium Halide Salts

Although halide salts such as LiCl and LiBr are routinely used as a source of Li ions during ion exchange reactions, a detailed understanding of the processes controlling the rates of these reactions is presently lacking. Recently, it has been discovered that the rate limiting barriers for ion exchange are commonly associated with these salts rather than the ceramic target of ion exchange, making it important to quantitatively understand of salt processes. Here it is demonstrated that in situ synchrotron studies of ion exchange reactions can be used to precisely quantify the thermodynamic activation energies associated with these solid state reactions in a manner that can be directly compared with predictions from density functional theory (DFT). While the temperature dependence of the LiCl reaction rate is found to be set by a barrier associated with ion hopping, it was discovered that for LiBr it is also affected by the defect formation energy – an energy found to be substantially lower than predicted by DFT. Furthermore, it is shown that when reaction rates for different relative amounts of reactants are varied, it is possible to identify the rate-limiting reagent and to elucidate an overall scaling relationship that controls the concentration-dependence of the reaction rate. Also, it is demonstrated that global fits across doped and undoped salts can be used to probe both intrinsic and extrinsic vacancy concentrations. This improved understanding of ion exchange mechanisms can be used to accelerate ion exchange reaction rates by orders of magnitude. In conclusion, the techniques demonstrated here can be broadly applied to probe the kinetics and thermodynamics of solid state reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE↗
DOE OSTI2023

TGA Study of the Reaction Kinetics of a CuFeMnAlO 4+δ Oxygen Carrier for Methane Chemical Looping Combustion under High Pressure and Atmospheric Pressure Conditions

Reduction and oxidation kinetics of a CuFeMnAlO 4+δ oxygen carrier for methane were studied using an atmospheric thermogravimetric analyzer (TGA) and a pressurized TGA. Important parameters such as temperature (500–900 °C), particle size distribution (150–650 μm), and partial pressures of reactants (methane and oxygen) were evaluated. Results showed that the reduction of the oxygen carrier was not affected by particle size distribution studied in this work, and that above 750 °C, the effect of the temperature on oxygen carrier reduction started to diminish. The reaction time approaches an asymptotic limit as methane partial pressure increases. For oxygen carrier oxidation, it was found that the rate of oxidation increased with increasing oxygen concentration. The effect of the oxygen concentration was more noticeable when the oxygen concentration was below 7.5%. Because of the inherent feature of the oxygen carrier, its chemical looping with oxygen uncoupling property was also assessed. A phase boundary-controlled shrinking sphere model was used to model the test data from the TGAs. Modeling results showed that the oxygen carrier reduction had a reaction order of 1 with an activation energy of ∼46 kJ/mol. Oxygen carrier oxidation demonstrated a reaction order of 1.1 with an activation energy of ∼55 kJ/mol. Additionally, the differences of the experimental results between the atmospheric pressure TGA and the pressurized TGA were discussed in detail and the limitations of the pressurized TGA were highlighted.

Hydrocarbons↗
DOE OSTI2020

Calcium-Uranyl-Carbonato Species Kinetically Limit U(VI) Reduction by Fe(II) and Lead to U(V)-Bearing Ferrihydrite

Reaction conditions and mechanisms promoting or inhibiting U reduction exert a central control on U solubility and, therefore, U transport and its associated risks. Here, we vary and track common aqueous uranium species to show that a kinetic restriction inhibits homogeneous reduction of the calcium-uranyl-carbonato species (CaUO$_2$(CO$_3$)$_3$$^{2–}$ and Ca$_2$UO$_2$(CO$_3$)$_3$) by Fe(II)$_{(aq)}$, while ferrihydrite surface-catalyzed reduction of all aqueous uranyl by Fe(II) proceeds. Using U L$_3$ high energy resolution fluorescence detection (HERFD) X-ray absorption near edge structure (XANES) spectroscopy, U L$_3$ extended X-ray absorption fine structure (EXAFS) spectroscopy, and transmission electron microscopy (TEM), we also show that U(V) is generated and incorporated into ferrihydrite formed from homogeneous U(VI) reduction by Fe(II)(aq). Through elucidation of the mechanisms that inhibit reduction of the calcium-uranyl-carbonato species and promote stabilization of U(V), we advance our understanding of the controls on U solubility and thus improve prediction of U transport in surface and subsurface systems.

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

Methane Catalytic Pyrolysis by Microwave and Thermal Heating over Carbon Nanotube-Supported Catalysts: Productivity, Kinetics, and Energy Efficiency

Methane catalytic pyrolysis, which is the reaction to produce hydrogen and carbon without emitting CO 2 , represents an approach for decarbonization using natural gas as an energy resource. In this work, the endothermic pyrolysis reaction was carried out under two heating scenarios: convective thermal heating and microwave-driven irradiative heating. The pyrolysis reaction was conducted at 550-600 °C over carbon nanotube-supported Ni-Pd and Ni-Cu catalysts. On both catalysts, an enhanced methane conversion rate was observed under microwave irradiation. The enhanced catalytic activity was hypothetically caused by the presence of free electrons in the carbon atoms within CNT that enabled the CNT support to absorb microwave energy effectively and to be heated efficiently by microwave. The microwave catalytic pyrolysis has shown improvement in kinetics, where the apparent activation energy dropped from 45.5 kJ/mol under conventional convective heating to 24.8 kJ/mol under microwave irradiation. When the methane conversion rate is increased by 37 %, the microwave power consumption only changed by 10.8 %. The research demonstrated the potential of transforming natural gas to clean hydrogen and value-added carbon in a more energy-efficient way. Process simulation and techno-economic analysis showed that potentially hydrogen minimum selling price of about $1 /kg H 2 could be achieved.

03 NATURAL GAS↗
DOE OSTI2026

A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability

This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.

petrochemicals↗
DOE OSTI2023

Changing Directions: Influence of Ligand Electronics on the Directionality and Kinetics of Photoinduced Charge Transfer in Cu(I)Diimine Complexes

A key challenge to the effective utilization of solar energy is to promote efficient photoinduced charge transfer, specifically avoiding unproductive, circuitous electron-transfer pathways and optimizing the kinetics of charge separation and recombination. We hypothesize that one way to address this challenge is to develop a fundamental understanding of how to initiate and control directional photoinduced charge transfer, particularly for earth-abundant first-row transition-metal coordination complexes, which typically suffer from relatively short excited-state lifetimes. Here, in this paper, we report a series of functionalized heteroleptic copper(I)bis(phenanthroline) complexes, which have allowed us to investigate the directionality of intramolecular photoinduced metal-to-ligand charge transfer (MLCT) as a function of the substituent Hammett parameter. Ultrafast transient absorption suggests a complicated interplay of MLCT localization and solvent interaction with the Cu(II) center of the MLCT state. This work provides a set of design principles for directional charge transfer in earth-abundant complexes and can be used to efficiently design pathways for connecting the molecular modules to catalysts or electrodes and integration into systems for light-driven catalysis.

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

Pynta-An Automated Workflow for Calculation of Surface and Gas-Surface Kinetics

Many important industrial processes rely on heterogeneous catalytic systems. However, given all possible catalysts and conditions of interest, it is impractical to optimize most systems experimentally. Automatically generated microkinetic models can be used to efficiently consider many catalysts and conditions. However, these microkinetic models require accurate estimation of many thermochemical and kinetic parameters. Manually calculating these parameters is tedious and error prone, involving many interconnected computations. Here, we present Pynta, a workflow software for automating the calculation of surface and gas–surface reactions. Pynta takes the reactants, products, and atom maps for the reactions of interest, generates sets of initial guesses for all species and saddle points, runs all optimizations, frequency, and IRC calculations, and computes the associated thermochemistry and rate coefficients. It is able to consider all unique adsorption configurations for both adsorbates and saddle points, allowing it to handle high index surfaces and bidentate species. Pynta implements a new saddle point guess generation method called harmonically forced saddle point searching (HFSP). HFSP defines harmonic potentials based on the optimized adsorbate geometries and which bonds are breaking and forming that allow initial placements to be optimized using the GFN1-xTB semiempirical method to create reliable saddle point guesses. This method is reaction class agnostic and fast, allowing Pynta to consider all possible adsorbate site placements efficiently. We demonstrate Pynta on 11 diverse reactions involving monodenate, bidentate, and gas-phase species, many distinct reaction classes, and both a low and a high index facet of Cu. Our results suggest that it is very important to consider reactions between adsorbates adsorbed in all unique configurations for interadsorbate group transfers and reactions on high index surfaces.

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

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2020

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗