Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Hydrogen oxidation mechanism with applications to (1) the chaperon efficiency of carbon dioxide and (2) vitiated air testing

Ignition delay times for the hydrogen/oxygen/carbon dioxide/argon system were obtained behind reflected shock waves. A detailed kinetic mechanism modeled the experimental hydrogen/oxygen data, Skinner and Ringrose's high-pressure data, and Slack and Grillo's hydrogen/air data. A carbon dioxide chaperon efficiency of 7.0 +/- 0.2 was determined. The reaction pathway H2O yields H2O2 yields OH yields H was required to model the high-pressure data. It is suggested that some of the lowest temperature data points (1.0 and 0.5 atm) for Slack and Grillo's hydrogen/air experiments are in error. It was found that the technique of simplifying a detailed kinetic mechanism for a limited range of experimental data may render the model useless for other test conditions.

Brabbs, Theodore A.↗
DOE OSTI2022

High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions

Abstract Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12// ω B97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol −1 ) relative to those calculated at ω B97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients $${k}_{\infty }(T)$$ k ∞ ( T ) between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.

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

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2024

Comparison of second-order spectral line widths formulae

Frequently used second-order spectral line width formulae from the projector operator and kinetic theory methods have not been formally compared previously. Here it is shown that a systematic second-order expansion of the projection operator expression including initial correlations agrees with the second-order kinetic theory result. The agreement assumes a common approximation in the projector operator method that introduces a screened radiator–perturber interaction to account for neglected electron–electron correlations. On the other hand, it is shown that the usual width expression from the projection operator approach neglecting initial correlations differs from kinetic theory. The differences, however, are at least third order in the radiator–perturber interaction. The comparisons suggest using the more compact width expression from kinetic theory, which includes initial correlations and a systematic screening of the radiator–perturber interactions, as the starting point for second-order width calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spontaneous Reorganization of DNA-Based Polymers in Higher Ordered Structures Fueled by RNA

We demonstrate a strategy that allows for the spontaneous reconfiguration of self-assembled DNA polymers exploiting RNA as chemical fuel. To do this, we have rationally designed orthogonally addressable DNA building blocks that can be transiently deactivated by RNA fuels and subtracted temporarily from participation in the self-assembly process. Through a fine modulation of the rate at which the building blocks are reactivated we can carefully control the final composition of the polymer and convert a disordered polymer in a higher order polymer, which is disfavored from a thermodynamic point of view. We measure the dynamic reconfiguration via fluorescent signals and confocal microscopy, and we derive a kinetic model that captures the experimental results. Our approach suggests a novel route toward the development of biomolecular materials in which engineered chemical reactions support the autonomous spatial reorganization of multiple components.

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

Probing the proton exchange kinetics of BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3− δ ceramic electrolyte by operando diffuse reflectance infrared Fourier transform spectroscopy

Proton exchange kinetics plays an important role in governing the performance of intermediate-temperature protonic ceramic electrolysis cells (PCECs) for hydrogen production. Our understanding of the nature of the surface hydration reaction at the single-cell level, however, remains very limited, hampering further efficiency improvements. Here, in this study, we developed a custom operando diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) platform that operates under high temperature and steam conditions with applied bias. Quantitative investigations of surface H 2 O/D 2 O isotope exchange in a BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3−δ (BZCYYb1711) protonic electrolyte-based single cell were conducted under different applied voltages using this DRIFTS platform, to gain molecular-level insight into hydration kinetics. The findings show that the application of an external voltage significantly enhances the surface proton exchange rate, decreasing the apparent activation energy from 29.1 kJ mol −1 at open-circuit voltage (OCV) to 6.8 kJ mol −1 at 1.3 V. In addition, distinct voltage-induced spectral shifts in O–D vibrations point to dynamic changes in surface hydration. These findings demonstrate a sensitive spectroscopic platform for probing interfacial proton processes and reveal strong electrochemical control over surface proton kinetics, offering new opportunities for probing electrolyte hydration behavior in PCECs.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

The Contribution of Proton-Donor pKa on Reactivity Profiles of [FeFe]-hydrogenases

The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H 2 coupled to the reduction–oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H 2 oxidation or H 2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H 2 oxidation or H 2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H 2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H 2 evolution. Overall, the changes in reactivity profiles were profound, spanning 10 5 in ratio of the H 2 oxidation-to-H 2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1979

Some aspects of condensation in clouds

Several general features of nucleation characteristics of low-density cosmic clouds are discussed. These are: (1) tendency for metastable condensates to form; (2) nonoccurrence of nominal refractory molecules in the gas; (3) a strong temperature dependence of condensation at relatively low temperatures; and (4) significant vibrational disequilibrium in cosmic clouds. These support previous analyses which indicate that equilibrium calculations have restricted applicability. A kinetic treatment of condensation is required for cosmic grains, and the possibility of formulating such an analysis is pointed out.

Donn, B.↗
NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
NASA NTRS1990

Precision tethered satellite attitude control

Tethered spacecraft possess unique dynamic characteristics which make them advantageous for certain classes of experiments. One use for which tethers are particularly well suited is to provide an isolated platform for spaceborne observatories. The advantages of tethering a pointing platform 1 or 2 km from a space shuttle or space station are that, compared to placing the observatory on the parent spacecraft, vibrational disturbances are attenuated and contamination is eliminated. In practice, all satellites have some requirement on the attitude control of the spacecraft, and tethered satellites are no exception. It has previously been shown that conventional means of performing attitude control for tethered satellites are insufficient for any mission with pointing requirements more stringent than about 1 deg. This is due mainly to the relatively large force applied by the tether to the spacecraft. A particularly effective method of implementing attitude control for tethered satellites is to use this tether tension force to generate control torques by moving the tether attach point relative to the subsatellite center of mass. A demonstration of this attitude control technique on an astrophysical pointing platform has been proposed for a space shuttle flight test project and is referred to as the Kinetic Isolation Tether Experiment (KITE).

Kline-Schoder, Robert J.↗
NASA NTRS1995

Rotational And Rovibrational Energy Transfer In Electron Collisions With Molecules

Air flows around a hypervelocity reentry vehicle undergo dissociation, rovibrational excitation and ionization. More specifically the air, initially 80% N2 and 20% O2, in the shock layer consists of species such as N, O, N2, O2, NO, N+, O+, N+, O+, NO+ and 2 free electrons. It was pointed out in multi temperature models'' that the temperature of the rotational energy modes and the gas-kinetic translational temperature are quickly equilibrated by a few collisions and rise rapidly to high temperatures as 50000K before falling off to equilibrium value of 10000K. Contrary, the electronic and vibrational temperatures state energy distributions remain low (less than 15000K) because of the slow equilibration. Electron vibrational energy transfer is thought to play a crucial role in such a ionizing flow regime since chemical reaction rates and dissociation depend strongly on the vibrational temperatures. Modeling of these flowfields in principle require the rovibrational excitation and de-excitation cross section data for average electron energies from threshold up to several eV (leV=11605.4 K). In this lecture we focus on theoretical description of rotational effects i.e. energy transfer of electrons to molecules such that the molecular rotational (vojo goes to voj) or vibrational and rotational (v(sub 0)j(sub 0) goes to vj) states are changed. Excitation and de-excitation of electronic states was discussed in a previous talk at this conference.

Thuemmel, Helmar T.↗
NASA NTRS2016

Transient Growth Analysis of Compressible Boundary Layers with Parabolized Stability Equations

The linear form of parabolized linear stability equations (PSE) is used in a variational approach to extend the previous body of results for the optimal, non-modal disturbance growth in boundary layer flows. This methodology includes the non-parallel effects associated with the spatial development of boundary layer flows. As noted in literature, the optimal initial disturbances correspond to steady counter-rotating stream-wise vortices, which subsequently lead to the formation of stream-wise-elongated structures, i.e., streaks, via a lift-up effect. The parameter space for optimal growth is extended to the hypersonic Mach number regime without any high enthalpy effects, and the effect of wall cooling is studied with particular emphasis on the role of the initial disturbance location and the value of the span-wise wavenumber that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary layer equations, mean flow solutions based on the full Navier-Stokes (NS) equations are used in select cases to help account for the viscous-inviscid interaction near the leading edge of the plate and also for the weak shock wave emanating from that region. These differences in the base flow lead to an increasing reduction with Mach number in the magnitude of optimal growth relative to the predictions based on self-similar mean-flow approximation. Finally, the maximum optimal energy gain for the favorable pressure gradient boundary layer near a planar stagnation point is found to be substantially weaker than that in a zero pressure gradient Blasius boundary layer.

Compressible boundary layer↗
NASA NTRS2004

Computation Molecular Kinetics Model of HZE Induced Cell Cycle Arrest

Cell culture models play an important role in understanding the biological effectiveness of space radiation. High energy and charge (HZE) ions produce prolonged cell cycle arrests at the G1/S and G2/M transition points in the cell cycle. A detailed description of these phenomena is needed to integrate knowledge of the expression of DNA damage in surviving cells, including the determination of relative effectiveness factors between different types of radiation that produce differential types of DNA damage and arrest durations. We have developed a hierarchical kinetics model that tracks the distribution of cells in various cell phase compartments (early G1, late G1, S, G2, and M), however with transition rates that are controlled by rate-limiting steps in the kinetics of cyclin-cdk's interactions with their families of transcription factors and inhibitor molecules. The coupling of damaged DNA molecules to the downstream cyclin-cdk inhibitors is achieved through a description of the DNA-PK and ATM signaling pathways. For HZE irradiations we describe preliminary results, which introduce simulation of the stochastic nature of the number of direct particle traversals per cell in the modulation of cyclin-cdk and cell cycle population kinetics. Comparison of the model to data for fibroblast cells irradiated photons or HZE ions are described.

Cucinotta, Francis A.↗
NASA NTRS1987

The carbon dioxide chaperon efficiency for the reaction H + O2 + M yields HO2 + M from ignition delay times behind reflected shock waves

Ignition delay times for stoichiometric hydrogen-oxygen in argon with and without carbon dioxide were measured behind reflected shock waves. A 20-reaction kinetic mechanism models the measured hydrogen-oxygen delay times over the temperature range 950 to 1300 K. The chaperon efficiency for carbon dioxide determined for the hydrogen-oxygen carbon dioxide mixture was 7.0. This value is in agreement with literature values but much less than a recent value obtained from flow tube experiments. Delay times measured behind a reflected shock wave were about 20% longer than those measured behind incident shock waves. The kinetic mechanism successfully modeled the high-pressure data of Skinner and the hydrogen-air data of Stack. It is suggested that the lowest temperature points for the hydrogen-air data of Slack are unreliable and that the 0.27-atm data may illustrate a case where vibrational relaxation of nitrogen is important. The reaction pathway HO2 yields H2O2 yields OH yields H was required to model the high-pressure data of Skinner. The successful modeling of the stoichiometric hydrogen-air data demonstrates the appropriateness of deriving kinetic models from data for gas mixtures highly diluted with argon. The technique of reducing a detailed kinetic mechanism to only the important reactions for a limited range of experimental data may render the mechanism useless for other test conditions.

Brabbs, Theodore A.↗
DOE OSTI2020

Global flow structure and exact formal transseries of the Gubser flow in kinetic theory

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗