Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics Parameters”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Spectroscopy and Kinetics of Combustion Gases at High Temperatures

This report describes our research program that involves two complementary activities: (1) development and application of cw laser absorption methods for the measurement of concentration time-histories and fundamental spectroscopic parameters for species of interest in combustion; and (2) shock tube studies of reaction kinetics relevant to combustion. This first part of this report covers research during the final three-year support period, i.e. March 2012 – November 2015. The later part of this report summarizes research conducted over multiple-year periods between March 1988 to March 2012. Publications supported by DOE for each period are summarized at the end of that report section.

10 SYNTHETIC FUELS↗
DOE OSTI2024

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1997

Kinetics and Thermochemistry of ClCO Formation from the Cl + CO Association Reaction

Laser flash photolysis of Cl2/CO/M mixtures (M = N2, CO, Ar, CO2) has been employed in conjunction with Cl((sup 2)P(sub J)) detection by time-resolved resonance fluorescence spectroscopy to investigate equilibration kinetics in the reactions Cl((sup 2)P(sub J)) + CO ClCO as a function of temperature (185-260 K) and pressure (14-200 Torr). The association and dissociation reactions are found to be in the low-pressure limit over the range of experimental conditions investigated. In N2 and/or CO buffer gases, the temperature dependences of the ClCO formation and dissociation reaction rate constants are described by the Arrhenius expressions k(sub 1) = (1.05 +/- 0.36) x 10(exp -34) exp[(810 +/- 70)/T] cm(exp 6)/molecules(exp 2).s and k(sub -1) = (4.1 +/- 3.1) x 10(exp -10) exp[(-2960 +/- 60)/T]cu cm/(molecule.s) (errors are 2 sigma). Second- and third-law analyses of the temperature dependence of the equilbrium constant (k/k-1) lead to the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -7.7 +/- 0.6 kcal/mol, Delta-H(sub 0) = -6.9 +/- 0.7 kcal/mol, Delta-S(sub 298) = -23.8 +/- 2.0 cal/mole.K, Delta-H(sub f,298)(ClCO) = 5.2 +/- 0.6 kcal/mol (errors are 2 sigma). The results repported in this study significantly reduce the uncertainties in all reported kinetic and thermodynamic parameters.

Nicovich, J. M.↗
DOE OSTI2019

Surface Wetting Controls Calcium Carbonate Crystallization Kinetics

Because of the widespread presence of foreign substrates in natural settings, mineral precipitation usually occurs via heterogeneous nucleation. This process is controlled by the interplay between the fluid supersaturation and interfacial energies present between the fluid, nucleus, and substrate. Among a number of physicochemical parameters, the surface wetting properties have been shown to be a key parameter controlling heterogeneous nucleation. The present study aims at elucidating the pathway and kinetics of CaCO 3 heterogeneous nucleation on a set of phlogopite micas with and without fluorine/hydroxyl substitutions, yielding substrates with contrasting hydrophilicity. Our results show that, irrespective of surface wetting properties, amorphous calcium carbonate (ACC) is formed during the early stages. The surface wetting properties have a strong effect on the crystallization kinetics: ACC precipitates persist longer on the hydrophilic (hydroxylated) surface than on the less hydrophilic (fluorinated) one. We show that this stabilization could have a thermodynamic origin because of the lower interfacial free energy between the hydrated amorphous precursor and the hydrophilic substrate. Furthermore, these results are highly relevant for biomineralization studies, where differences in wetting properties of organic moieties present in calcifying organisms could be used to accelerate or decelerate the crystallization of the initially formed amorphous precursor phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗
DOE OSTI2016

Effects of Lewis acidity of metal oxide promoters on the activity and selectivity of Co-based Fischer–Tropsch synthesis catalysts

Metal oxides of Ce, Gd, La, Mn, and Zr were investigated as promoters for improving the activity and selectivity of Co-based FTS catalysts. The extent to which these promoters decrease the selectivity toward CH 4 and increase the selectivity toward C 5+ hydrocarbons was found to depend on both the loading and the composition of the oxide promoter. Elemental mapping by STEM–EDS revealed that the propensity for a given metal oxide to associate with Co affects the sensitivity of the product distribution to changes in promoter loading. For all promoters, a sufficiently high loading resulted in the product distributions becoming insensitive to further increases in promoter loading, very likely due to the formation of a half monolayer of promoter oxide over the Co surface. Simulations suggest that the fraction of Co active sites that are adjacent to the promoter moieties approaches unity at this degree of coverage. The oxidation state of the promoter metal cation under reaction conditions, determined by in situ XANES measurements, was used to calculate relative Lewis acidity of the promoter metal cation. A strong positive correlation was found between the C 5+ product selectivity and the Lewis acidity of the promoter metal cations, suggesting that the promotional effects are a consequence of Lewis acid–base interactions between the reaction intermediates and the promoter metal cations. Rate data obtained at different pressures were used to estimate the apparent rate coefficient and the CO adsorption constant appearing in the Langmuir–Hinshelwood expression that describes the CO consumption kinetics for both unpromoted and the metal oxide-promoted catalysts. Both parameters exhibited positive correlations with the promoter Lewis acidity. In conclusion, these results are consistent with the hypothesis that the metal cations of the promoter act as Lewis acids that interact with the O atom of adsorbed CO to facilitate CO adsorption and dissociation.

Johnson, Gregory R.↗
DOE OSTI2019

Cell free biosynthesis of isoprenoids from isopentenol

Cell-free systems are growing in importance for the biosynthesis of complex molecules. These systems combine the precision of traditional chemistry with the versatility of biology in creating superior overall processes. Recently, a new synthetic pathway for the biosynthesis of isoprenoids using the substrate isopentenol, dubbed the isopentenol utilization pathway (IUP), was demonstrated to be a promising alternative to the native 2C-methyl-d-erythritol-4-phosphate (MEP) and mevalonate (MVA) pathways. This simplified pathway, which contains a minimum of four enzymes to produce basic monoterpenes and only depends on ATP and isopentenol as substrates, allows for a highly flexible approach to the commercial synthesis of isoprenoid products. In this work, we use metabolic reconstitution to characterize this new pathway in vitro and demonstrate its use for the cell-free synthesis of mono-, sesquit-, and diterpenoids. Kinetic modeling and sensitivity analysis were also used to identify the most significant parameters for taxadiene productivity, and metabolic control analysis was employed to elucidate protein-level interactions within this pathway, which demonstrated that the IUP enzymatic system is primarily controlled by the concentration and kinetics of choline kinase (CK) and not regulated by any pathway intermediates. We report this is a significant advantage over the natural MEP or MVA pathways as it greatly simplifies future metabolic engineering efforts, both in vitro and in vivo, aiming at improving the kinetics of CK. Finally, we used the insights gathered to demonstrate an in vitro IUP system that can produce 220 mg/L of the diterpene taxadiene, in 9 hr, almost 3-fold faster than any system reported thus far.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Molybdenum Hydride and Dihydride Complexes Bearing Diphosphine Ligands with a Pendant Amine: Formation of Complexes With Bound Amines

CpMo(CO)(PNP)H complexes (PNP = (R2PCH2)2NMe, R = Et or Ph) were synthesized by displacement of two CO ligands of CpMo(CO)3H by the PNP ligand; these complexes were characterized by IR and variable temperature 1H and 31P NMR spectroscopy. CpMo(CO)(PNP)H complexes are formed as mixture of cis and trans-isomers. Both cis-CpMo(CO)(PEtNMePEt)H and trans-CpMo(CO)(PPhNMePPh)H were analyzed by single crystal X-ray diffraction. Electrochemical oxidation of CpMo(CO)(PEtNMePEt)H and CpMo(CO)(PPhNMePPh)H in CH3CN are both irreversible at slow scan rates and quasi-reversible at higher scan rates, with E1/2 = -0.36 V (vs. Cp2Fe+/0) for CpMo(CO)(PEtNMePEt)H and E1/2 = -0.18 V for CpMo(CO)(PPhNMePPh)H. Hydride abstraction from CpMo(CO)(PNP)H with [Ph3C]+[A]- (A = B(C6F5)4 or BArF4; [ArF = 3,5-bis(trifluoromethyl)phenyl]) afforded “tuck-in” [CpMo(CO)(?3-PNP)]+ complexes that feature the amine bound to the metal. Displacement of the ?3 Mo-N bond by CD3CN gives [CpMo(CO)(PNP)(CD3CN)]+. The kinetics of this reaction were studied by NMR spectroscopy, providing the activation parameters ?H‡ = 22.1 kcal/mol, ?S‡ = 1.89 cal/(mol·K), Ea = 22.7 kcal/mol. Protonation of CpMo(CO)(PEtNMePEt)H affords [CpMo(CO)(?2-PEtNMePEt)(H)2]+ as a Mo dihydride complex, which loses H2 to generate [CpMo(CO)(?3-PEtNMePEt)]+ at room temperature. CpMo(CO)(dppp)H (dppp = 1,2-bis(diphenylphosphino)propane) was studied as a Mo diphosphine analogue without a pendant amine, and the product of protonation of this complex gives [CpMo(CO)(dppp)(H)2]+. Our results show that the pendant amine has a strong driving force to form stable “tuck-in” [CpMo(CO)(?3-PNP)]+ complexes, and also promotes hydrogen elimination from [CpMo(CO)(PNP)(H)2]+ complexes by formation of Mo-N dative bond. We thank the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences for support. Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy. We thank Monte Helm, Elliott Hulley and Deanna Miller for help on the crystallography, and Ming Fang for assistance on the electrochemical experiments.

Zhang, Shaoguang↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stochastic Model for High Temperature Oxidation of Cr–Ni Austenitic Steels Assisted by Spallation

Abstract Cr–Ni austenitic steels offer significant high temperature corrosion protection by forming a surface oxide layer. However, above critical service conditions (temperature, atmosphere, thermal cycling), oxidized surface can experience intensive degradation because of scale spallation, which could be detrimental to the in-service life. To predict the effect of scale spallation on oxidation kinetics, a simulation was implemented using a stochastic model. The model considers topological parameters and intensity of spallation which can occur, while delivering a true oxidation constant. The experimental procedure identified the amount of formed spalled scale and topology of spallation based on the use of element mapping of the surface. This information was used to determine a true kinetic constant for a corresponding spallation intensity in oxidized Cr–Ni austenitic steel. To illustrate the capability of the stochastic model, a parametric analysis was performed. The model verified how the spallation parameters could change the oxidation processes from parabolic growth of an adhered oxide layer without spallation to a mixed linear-parabolic, or with a constant thickness of residual scale at high spallation intensity. The spallation model will be used in a separate article to characterize high temperature surface degradation of several Cr–Ni austenitic steels during harsh oxidation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Software Requirements for a Point Kinetics Neutronics Code for Simulating TREAT Reactor Transients

This document provides the software requirements for a computer code to perform point kinetics-based simulations of nuclear reactor transients. These requirements address all aspects of code development and use, including input and output parameters and structure, software capabilities, and supporting documentation. This software is intended to replace the existing point kinetics code TREKIN, which is currently used for the analysis of the Transient Reactor Test Facility (TREAT).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Uniformity of GaAs solar cells grown in a kinetically-limited regime by dynamic hydride vapor phase epitaxy

We demonstrate solar cell devices grown on 50 mm diameter GaAs substrates by dynamic hydride vapor phase epitaxy (D-HVPE). In contrast to our prior D-HVPE devices, grown at 650 °C in a transport-limited regime, these devices were grown at 700 °C in a kinetically-limited growth regime in which the growth rate uniformity is controlled by thermally-activated surface reactions. These devices exhibit open-circuit voltages up to 1.07 V, nearly identical performance to the devices grown at lower temperature in a different growth regime. We evaluate the uniformity of device performance and find only a 1% variation in absolute, device efficiency across the majority of the wafer in devices without anti-reflection coating. We analyze the GaAs and GaInP thickness uniformity and GaInP compositional uniformity, using high-tresolution x-ray diffraction mapping, and find that any non-uniformity in device efficiency is not related to variations in these structural parameters. We combine three-dimensional computational fluid dynamics modeling of our reactor with a kinetic model for GaAs growth, and use the combined model to predict spatial GaAs growth rate over a 50 mm wafer area. We compare these predictions with experimental data from our D-HVPE reactor and find excellent agreement. We use the model to the gain insight into the specific mechanisms that control GaAs spatial uniformity in the kinetically-limitedgrowth regime. Furthermore, these results suggest a large parameter window for the growth of high-performance optoelectronic devices by D-HVPE, possibly with large area uniformity.

14 SOLAR ENERGY↗
DOE OSTI2016

A HELIOSEISMIC SURVEY OF NEAR-SURFACE FLOWS AROUND ACTIVE REGIONS AND THEIR ASSOCIATION WITH FLARES

We use helioseismic holography to study the association of shallow flows with solar flare activity in about 250 large sunspot groups observed between 2010 and 2014 with the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory. Four basic flow parameters: horizontal speed, horizontal component of divergence, vertical component of vorticity, and a vertical kinetic helicity proxy, are mapped for each active region (AR) during its passage across the solar disk. Flow indices are derived representing the mean and standard deviation of these parameters over magnetic masks and compared with contemporary measures of flare X-ray flux. A correlation exists for several of the flow indices, especially those based on the speed and the standard deviation of all flow parameters. However, their correlation with X-ray flux is similar to that observed with the mean unsigned magnetic flux density over the same masks. The temporal variation of the flow indices are studied, and a superposed epoch analysis with respect to the occurrence to 70 M and X-class flares is made. While flows evolve with the passage of the ARs across the disk, no discernible precursors or other temporal changes specifically associated with flares are detected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

An adjoint method for neoclassical stellarator optimization

Stellarators remain an enticing route to steady-state fusion power. However, to achieve the required confinement, the magnetic geometry must be highly optimized. This optimization requires navigating high-dimensional spaces, often necessitating the use of gradient-based methods. The gradient of the neoclassical fluxes is expensive to compute with classical methods, requiring $O(N)$ flux computations, where $N$ is the number of parameters. To reduce the cost of the gradient computation, we present an adjoint method for computing the derivatives of moments of the neoclassical distribution function for stellarator optimization. The linear adjoint method allows derivatives of quantities which depend on solutions of a linear system, such as moments of the distribution function, to be computed with respect to many parameters from the solution of only two linear systems. This reduces the cost of computing the gradient to the point that the finite-collisionality neoclassical fluxes can be used within an optimization loop. With the neoclassical adjoint method, we compute solutions of the drift kinetic equation and an adjoint drift kinetic equation to obtain derivatives of neoclassical quantities with respect to geometric parameters. When the number of parameters in the derivative is large ( $O(10^{2})$ ), this adjoint method provides up to a factor of 200 reduction in cost. We demonstrate adjoint-based optimization of the field strength to obtain minimal bootstrap current on a surface. With adjoint-based derivatives, we also compute the local sensitivity to magnetic perturbations on a flux surface and identify regions where tight tolerances on error fields are required for control of the bootstrap current or radial transport. Moreover, the solve for the ambipolar electric field is accelerated using a Newton method with derivatives obtained from the adjoint method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗