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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis

Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Elliott, Sarah N.↗
DOE OSTI2021

Hamiltonian kinetic-Hall magnetohydrodynamics with fluid and kinetic ions in the current and pressure coupling schemes

We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Control of Photoswitching Kinetics with Strong Light–Matter Coupling in a Cavity

Most photochemistry occurs in the regime of weak light–matter coupling, in which a molecule absorbs a photon and then performs photochemistry from its excited state. In the strong coupling regime, enhanced light–matter interactions between an optical field and multiple molecules lead to collective hybrid light–matter states called polaritons. This strong coupling leads to fundamental changes in the nature of the excited states including multi-molecule delocalized excitations, modified potential energy surfaces, and dramatically altered energy levels relative to non-coupled molecules. Here, the effect of strong light–matter coupling on covalent photochemistry has not been well explored. Photoswitches undergo reversible intramolecular photoreactions that can be readily monitored spectroscopically. In this work, we study the effect of strong light–matter coupling on the kinetics of photoswitching within optical cavities. Reproducing prior experiments, photoswitching of spiropyran/merocyanine photoswitches is decelerated in a cavity. Fulgide photoswitches, however, show the opposite effect, with strong coupling accelerating photoswitching. While modified merocyanine switching can be explained by changes in radiative decay rates or the amount of light in the cavity, modified fulgide switching kinetics suggest direct changes to excited-state reaction kinetics.

Cavities↗
DOE OSTI2022

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Isomeric effects on the reactivity of branched alkenes: An experimental and kinetic modeling study of methylbutenes

Here, a detailed experimental study of the low-to-intermediate temperature combustion of methylbutene isomers, i.e., branched C 5 alkenes, has been undertaken with multiple experimental facilities. Ignition delay times were measured at equivalence ratios 0.5–2.0, 685–1020 K and up to 45 bar condition from two rapid compression machines and showed slight deviation from an Arrhenius behavior for all three isomers, while their reactivity order differs as temperature changes. Sampled intermediates formed during the oxidation process of mixtures at 900–1150 K and 0.82 bar from a flow reactor and at 730 K and 20 bar from a rapid compression machine were analyzed using gas chromatography techniques. Trends in the formation and consumption of sampled intermediates were modeled using a kinetic model developed in this work for all three isomers. Rate of production and sensitivity analyses emphasize the role of double bond-specific reactions governing the global reactivity of these fuels. Additional studies of the addition reactions of HO 2 radicals to the double bond and to allylic radicals may improve the model performance.

2-Methyl-1-butene↗
DOE OSTI2024

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The role of ternary alloying elements in eutectoid transformation of U-10Mo alloy part II. In and Ex-situ Neutron diffraction-based assessment of eutectoid phase transformation kinetics in U-9.8Mo-0.2X alloy (X = Cr, Ni or Co)

HExploring the effects of minor ternary alloying additions, typically impurity elements, on the phase stability of U-10Mo is important for preventing undesirable phase decomposition during processing or during service. This work examines the influence small ternary additions of Cr, Ni, and Co. Both in-situ and ex-situ neutron diffraction measurements made during and after high temperature (450 – 525°C) exposures were used to better define the influence of these elements on the time-temperature-transformation (TTT) behavior of U-10Mo, providing information which is complementary to electron microscopy investigations of the same alloy systems performed in the first part of this work. Minor additions of Ni and Co decrease the ?-phase stability at all temperatures investigated. Signatures of U6X (X = Ni or Co) compounds were shown to be present in amounts of up to 6 wt% in the heat treated alloys, suggesting that the initial precipitation of this phase may catalyze further ?-phase decomposition. On the other hand, the Cr containing alloys were observed to have nearly the same, and in some cases slower, phase transformation kinetics when compared to the binary U-10Mo control samples. The results of the study have enabled preliminary estimates of the TTT curves for the ternary alloys.

phase transformation, kinetics, discontinuous prec↗
DOE OSTI2024

Styrene Thermal Decomposition and Its Reaction with Acetylene under Shock Tube Pyrolysis Conditions: an Experimental and Kinetic Modeling Study

Styrene is an important compound for polymer production and a key intermediate in gas-phase kinetics of polycyclic aromatic hydrocarbons (PAHs). For the first time, the pyrolysis of styrene with and without the presence of acetylene is investigated in a single-pulse shock tube coupled to gas chromatography and mass spectrometry. For each reaction system, quantitative speciation profiles are probed within the temperature range of 1100-1730 K, nominal pressure of 20 bar, and reaction duration of similar to 4 ms. A kinetic model is built to simulate the results. The model explains how styrene is consumed under high-pressure pyrolytic conditions, how the secondary chemistry of intermediate products affect subsequent PAH formation, and how acetylene addition alters the reaction pathways. Throughout the temperature range, styrene breakdown is dominated by the bimolecular interaction between styrene and hydrogen atom, which produces benzene+vinyl or phenyl and ethylene through the stabilization of 2-phenylethyl and its subsequent dissociation. As a result, large amounts of phenyl accumulate, which react with styrene to form C14H12 species while simultaneously releasing H atoms through addition/elimination reactions. The reactivity of fuel consumption is preserved by the regeneration of H atoms as chain carriers. Several C14H10 compounds are formed as a result of the following breakdown of the C14H12 isomers, particularly stilbene, 1,1-diphenyl ethylene, and 9-methyl-9H-fluorene. The presence of acetylene as a co-reactant with styrene allows the Hydrogen-Abstraction-Acetylene-Addition (HACA) pathway to proceed from phenyl radical to enhance the production of phenylacetylene at very low temperatures and acenaphthylene. This hinders the formation of C14H12 isomers, exclusive products from pure styrene dissociation by competing with the styrene+phenyl routes.

Kinetic Modeling↗
DOE OSTI2025

Kinetics for the reaction between the solvated electron and dissolved oxygen in n-dodecane from 2.5 to 40 °C

Temperature-controlled, time-resolved picosecond electron pulse radiolysis was utilized to measure the rate of reaction between the solvated electron (eS–) and dissolved oxygen in n-dodecane solutions from 2.5 to 40 °C for the first time. At 20.0 °C, the reaction rate was determined to be k(eS– + O2) = (4.54 ± 0.21) × 1010 M-1 s-1, with an activation energy of Ea = 14.4 ± 1.3 kJ mol-1. These newly determined kinetic parameters are important for predicting and managing the effects of aerated environments on the degradation of organic solvents used in nuclear fuel reprocessing technologies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Kinetic Probes of the Origin of Activity in MOF-Based C–H Oxidation Catalysis

The development of homogeneous catalysis is enabled by the availability of a rich toolkit of kinetics experiments, such as the Hgdrop test, that differentiate catalytic activity at ligand-supported metal complexes from potential heterogeneous catalysts derived from decomposition of molecular species. Metal-organic frameworks (MOFs) have garnered significant attention as platforms for catalysis at site-isolated, interstitial catalyst sites. Unlike homogeneous catalysis, relatively few strategies have been advanced to evaluate the origin of catalytic activity in MOF-promoted reactions. Many of the MOFs that have been examined as potential catalysts are comprised of molecular constituents that represent viable catalysts in the absence of the extended MOF lattice, and thus interfacial sites and leached homogeneous species represent potential sources of catalyst activity. Here, we demonstrate that analysis of deuterium kinetic isotope effects (KIEs) and olefin epoxidation diastereoselectivity provides probes of the origin of catalytic activity in MOF-promoted oxidation reactions. These analyses support involvement of lattice-based Fe sites in the turnover-limiting step of C–H activation with Fe-MOF-74-based materials (i.e., the MOF functions as a bona fide catalyst) and evidence that Cu 2 -based MOF MIL-125-Cu 2 O 2 functions as a solid-state initiator for solution-phase oxidation chemistry and is not involved in the turnover limiting step (i.e., the MOF does not function as a catalyst for substrate functionalization). furthermore, we anticipate that the simple experiments described here will provide a valuable tool for clarifying the role of MOFs in C–H oxidation reactions.

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

Acidity of surface-infiltrated binary oxides as a sensitive descriptor of oxygen exchange kinetics in mixed conducting oxides

Improving the kinetics of O 2 reduction on oxide surfaces is critical in many energy and fuel conversion technologies. Here we show that the acidity scale for binary oxides is a powerful descriptor for tuning and predicting oxygen surface exchange kinetics on mixed conducting oxides. By infiltrating a selection of binary oxides from strongly basic (Li 2 O) to strongly acidic (SiO 2 ) onto the surface of Pr 0.1 Ce 0.9 O 2-δ samples, it was possible to vary the chemical surface exchange coefficient k chem by 6 orders of magnitude, with basic oxides such as Li 2 O increasing k chem by nearly 1,000 times, with surface concentrations as low as 50 ppm impacting k chem . Strikingly, although the pre-exponential of k chem scales linearly with the acidity of the infiltrated binary oxide, there is nearly no change in the activation energy. Here, the origin of these dramatic changes is proposed to arise from the systematic increase in electron concentration at the Pr 0.1 Ce 0.9 O 2-δ surface with the decreasing acidity of the infiltrated binary oxide.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

A kinetic model of friction in strongly coupled strongly magnetized plasmas

Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge, in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that the “gyro” component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyrofriction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2020

Kinetics and Mechanisms of ZnO to ZIF-8 Transformations in Supercritical CO2 Revealed by in situ X-ray Diffraction

ZIF-8 was synthesized in supercritical carbon dioxide (scCO2). In situ powder X-ray diffraction, ex situ microscopy and simulations, provides a encompassing view of the formation of ZIF-8 and intermediary ZnO@ZIF-8 composites in this nontraditional solvent. Time-resolved imaging exposed divergent physicochemical reaction pathways from previous studies of the growth of anisotropic ZIF-8 core@shell structures in traditional solvents. Synthetically relevant physio-chemical properties of scCO2 were integrated into classical nucleation theory, relating interfacial forces with 3D nu-cleation outcomes. The kinetics of crystallization were examined and displayed a characteristic signature of time- and temperature-dependent mechanisms over the extent of the reaction. Lastly, we show that subtle factors, such as the extent of reaction and the size/shape of sacrificial templates can tailor ZIF-8 composition and size, eliciting control over hierarchical porosity in a nonconventional green solvent.

Sinnwell, Michael A.↗
DOE OSTI2025

Kinetic and structural evidence for specific DMSO interference with reversible binding of uncharged bis-oximes to hAChE and their reactivation kinetics of OP-hAChE

The structural basis of inhibitory effect of organic solvent dimethyl sulfoxide (DMSO) on human acetylcholinesterase (EC 3.1.1.7; hAChE) was inferred from the effect of DMSO on kinetics of reversible inhibition of uncharged, heterocyclic bis-oximes to hAChE, from DMSO effect on rates of reactivation of inactive organophosphate (OP)-hAChE conjugates by bis-oximes and by X-ray structures of bis-oxime and DMSO binding to hAChE. Here, the reversible inhibition constant of DMSO for hAChE in 0.1 M phosphate buffer pH 7.4 at 22 °C, was K i = (0.32 ± 0.04) % (or 45 ± 5 mM). The K i of the bis-oxime LG-703 for hAChE was 3.2-fold larger in 1 % DMSO, consistent with direct competition between LG-703 and DMSO. The X-ray structure of the LG-703∗hAChE complex (PDB ID: 6U3P) shows DMSO and LG-703 bound to individual hAChE monomers, LG-703 in the chain A and DMSO in the chain B. In the co-crystallization both small molecules were present at a similar excess over their corresponding K i values for hAChE (7.8-fold for DMSO and 6.5-fold for LG-703) and formation of two different complexes (DMSO∗hAChE and LG-703∗hAChE), in the same crystal, appears consistent with inhibition kinetics. Furthermore, rates of reactivation of paraoxon-inhibited hAChE (POX-hAChE) and of VX-hAChE by LG-703 and by a novel heterocyclic bis-oxime LG-1922 were reduced 2 – 3-fold in DMSO, consistent with observation of the active-center-bound DMSO molecules in the newly solved structure of the LG-1922∗POX-hAChE complex presented here and in our POX-hAChE structure (PDB ID: 8DT2) showing obstruction of the reactivator access to the conjugated P atom.

Acetylcholinesterase inhibition↗
DOE OSTI2022

Kinetics and Reaction Mechanisms of Acetic Acid Hydrodeoxygenation over Pt and Pt–Mo Catalysts

In this work, kinetic measurements for silica-supported Pt and Pt-Mo catalysts were collected in vapor-phase acetic acid hydrodeoxygenation by varying the hydrogen partial pressure between 18 and 72 kPa and acetic acid partial pressure between 7 and 18 kPa at 423-473 K. At all testing conditions, the Pt-Mo catalyst was more active and selective. In addition, the apparent activation energy for Pt-Mo of 76±3 kJ/mol was lower than that of 84±4 kJ/mol for Pt. The apparent reaction orders were also different. The order in hydrogen of 0.8±0.1 for Pt-Mo changed to zero at higher hydrogen partial pressures while that for Pt remained constant at 0.6±0.1. A near-zero order in acetic acid for Pt-Mo changed to -2.1 at higher acetic acid pressures while that for Pt remained constant at -2.9±0.3. These differences in reaction kinetics as well as in selectivity trends with changes in temperature and feed composition indicated a change in the reaction mechanism for Pt-Mo. The catalysts were characterized with hydrogen temperature programmed desorption, oxygen temperature programmed oxidation and transmission electron microscopy with energy-dispersed X-ray spectroscopy elemental mapping. Mo was present in the form of subnanometer-size clusters on the surface of Pt nanoparticles. Both Pt and Pt-Mo catalysts were stable under the reaction conditions for 10 h, and the size and structure of Pt and Pt-Mo particles remained mostly unchanged, without coke accumulation. Density functional theory calculations show that surface acetate is not a major reaction intermediate on both Pt and Pt-Mo and, instead, C-OH bond splitting with the formation of acetyl is the first reaction step in hydrodeoxygenation. The activation energy for this step is dramatically lower on Pt-Mo. Furthermore, the activity of acetyl on Pt-Mo is different. Unlike on pure Pt, the reaction of C-O bond splitting becomes exothermic with a lower activation energy on Pt-Mo, and the reaction of C-C bond splitting, in contrast, becomes endothermic with a higher activation energy, explaining the experimentally observed higher activities and selectivities of Pt-Mo. In addition, the calculations demonstrate that a pair of Pt-Mo surface atoms acts as a single active site where Mo serves as a preferential binding anchor for O atoms. The presence of Mo atoms changes the structure and energy of adsorbed and reacting surface species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗