Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic model data”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Homomolecular Triplet–Triplet Annihilation in Metalloporphyrin Photosensitizers

Metalloporphyrins are ubiquitous in their applications as triplet photosensitizers, particularly for promoting sensitized photochemical upconversion processes. In this study, bimolecular excited state triplet–triplet quenching kinetics, termed homomolecular triplet–triplet annihilation (HTTA), exhibited by the traditional triplet photosensitizers–zinc(II) tetraphenylporphyrin (ZnTPP), palladium(II) octaethylporphyrin (PdOEP), platinum(II) octaethylporphyrin (PtOEP), and platinum(II) tetraphenyltetrabenzoporphyrin (PtTPBP)–were revealed using conventional transient absorption spectroscopy. Nickel(II) tetraphenylporphyrin was used as a control sample as it is known to be rapidly quenched intramolecularly through ligand-field state deactivation and, therefore, cannot result in triplet–triplet annihilation (TTA). The single wavelength transients associated with the metalloporphyrin triplet excited state decay–measured as a function of incident laser pulse energy in toluene–were well modeled using parallel first- and second-order kinetics, consistent with HTTA being operable. The combined transient kinetic data enabled the determination of the first-order rate constants (k T ) for excited triplet decay in ZnTPP (4.0 × 10 3 s –1 ), PdOEP (3.6 × 10 3 s –1 ), PtOEP (1.2 × 10 4 s –1 ), and PtTPBP (2.1 × 10 4 s –1 ) as well as the second-order rate constant (k TT ) for HTTA in ZnTPP (5.5 × 10 9 M –1 s –1 ), PdOEP (1.1 × 10 10 M –1 s –1 ), PtOEP (7.1 × 10 9 M –1 s –1 ), and PtTPBP (1.6 × 10 10 M –1 s –1 ). In most instances, triplet excited state extinction coefficients are either reported for the first time or have been revised using ultrafast transient absorption spectroscopy and singlet depletion: ZnTPP (78,000 M –1 cm –1 ) at 470 nm, PdOEP (67,000 M –1 cm –1 ) at 430 nm, PtOEP (51,000 M –1 cm –1 ) at 418 nm, and PtTPBP (100,000 M –1 cm –1 ) at 460 nm. Furthermore, the combined experimental results establish competitive time scales for homo- and heteromolecular TTA rate constants, implying the significance of considering HTTA processes in future research endeavors harnessing TTA photochemistry using common metalloporphyrin photosensitizers.

14 SOLAR ENERGY↗
DOE OSTI2021

Long-Term 13 C Uptake by 12 C-Enriched Calcite

Knowledge of the exchange of carbon isotopes between dissolved inorganic carbon and calcite minerals is of long-standing importance for the interpretation of sedimentary paleoclimate records and 14 C transport in the geosphere. To assess the mechanism and rates of carbon isotope exchange, we equilibrated 12 C-pure synthetic calcite particles in water, first in a glovebox and then in contact with atmospheric P CO 2 and 13 C/ 12 C ratio, at two different temperatures. Cavity ring-down infrared spectroscopy δ 13 C measurements of the solid revealed sustained 13C incorporation for over a period of 500 days (21 °C) and 125 days (50 °C). We developed a quantitative model for recrystallization and isotope exchange, assuming that the interfacial free energy provides a thermodynamic driving force for the growth of larger particles at the expense of smaller ones. Furthermore, this Ostwald ripening model did not reproduce the kinetics of 13C uptake and required greater coarsening than observed. Rather, the data were best explained by a mechanism involving surface exchange and solid-phase diffusion of 13C into the particles with an inferred effective diffusion constant at 21 °C of about 10 –25 m 2 /s. Although this work cannot rule out the possibility that structural or chemical aspects of the synthetic particles enabled faster 13 C uptake than could be observed in natural systems, this study adds to the body of the recent work, suggesting that fast exchange processes are possible, likely through grain boundaries and other defects.

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

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

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

Investigation of the Potential of Repurposing Medium-Density Fiberboard Waste as an Adsorbent for Heavy Metal Ion Removal

Medium-density fiberboard (MDF) waste generation has increased steadily over the past decades, and therefore, the investigation of novel methods to recycle this waste is very important. The potential of repurposing MDF waste as an adsorbent for the treatment of Cd(II), Cu(II), Pb(II), and Zn(II) ions in water was investigated using MDF offcuts. The highest adsorption potential in single-metal ion solution systems was observed for Pb(II) ions. The experimental data of Pb(II) ions fit well with the Freundlich isotherm and pseudo-second-order kinetic models. Complexation and electrostatic interactions were identified as the adsorption mechanisms. The adsorption behavior of multi-metal ion adsorption systems was investigated by introducing Cd(II) ions as a competitive metal ion. The presence of the Cd(II) ions reduced the adsorption potential of Pb(II) ions, yet the preference for the Pb(II) ions remained. Regeneration studies were performed by using 0.1 M HCl as a regeneration agent for both systems. Even though a significant amount of adsorbed metal ions were recovered, the adsorption potential of the MDF was reduced in the subsequent adsorption cycles. Based on these results, MDF fines have the potential to be used as an economical adsorbent for remediation of wastewater containing heavy metal ions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LigninGraphs: lignin structure determination with multiscale graph modeling

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

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

Knock Mitigation Effectiveness of EGR across the Pressure-Temperature Domain

Exhaust gas recirculation (EGR) has been shown to enable efficiency improvements in SI engines through multiple different mechanisms, including decreasing the knock propensity at high load, which allows higher compression ratio. While many of the benefits of EGR are applicable to both low and high power density engines, including reductions in pumping work and improved specific heat ratio, the knock benefits and corresponding compression ratio increases have been limited to low power density naturally aspirated engines primarily intended for hybrid vehicle architectures. An earlier study [1] indicated that there may be a kinetic limitation for the ability of EGR to mitigate knock under these conditions, but that study only considered a small number of conditions. In this investigation, we expand on that study while also providing data for model validation for the new light-duty combustion consortium from the U.S. Department of Energy: Partnership for Advancing Combustion Engines (PACE). In this investigation, the effectiveness of EGR to mitigate knock is studied with regards to the effect of engine speed (1,500 and 3,000 rpm), changing trajectory in the pressure-temperature domain by varying the intake manifold temperature (35, 60, and 90 deg C), and by considering the effect of minor species by studying the effect of untreated EGR vs. EGR that has been treated by an automotive three-way catalyst. Additionally, to increase the relevance of these data for future modeling studies, the performance of the full boiling range gasoline was compared relative to a surrogate formulation. The study found that the fuel surrogate performs well, confirmed the kinetic limitations of EGR to mitigate knock under boost, and showed improvements in EGR performance with catalyzed EGR.

42 ENGINEERING↗
DOE OSTI2022

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

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

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2024

Structure Sensitive Reaction Kinetics of Chiral Molecules on Intrinsically Chiral Surfaces

Enantiospecific heterogeneous catalysis utilizes chiral surfaces to resolve enantiomers via structure sensitive surface chemistry. The catalyst design challenge is the identification of chiral surface structures that maximize enantiospecificity. Herein, we develop data driven models for the enantiospecificity of tartaric acid reactions on chiral Cu(hkl) R&S surfaces. Measurements of enantiospecific rate constants were obtained by using curved Cu(hkl) R&S surfaces that enable kinetic measurements on hundreds of chiral surface orientations. One model uses feature vectors derived from generalized coordination numbers to capture the local structure around Cu atoms exposed by the Cu(hkl) R&S surfaces. The second model introduces the use of chiral cubic harmonic functions to capture the symmetry constraints of the face-centered cubic Cu structure. The model using 58 generalized coordination numbers has a fitting error similar to that of the model using only 5 cubic harmonic functions. The two models predict maxima in the enantiospecificity on surfaces with very similar surface orientations. The models developed in this work are applicable for any enantiospecific reaction happening on any chiral material with a cubic lattice structure, opening the way to understanding the surface structure sensitivity of the enantiospecific reaction kinetics.

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

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

ATOMIC Simulations and Experimental Data for Sodium and Copper Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1]. Our overall suite of simulations contains 12 sets of simulations: training and validation sets of simulations for each of two resolutions of each of sodium, copper, and mixtures of the two. The training data were produced using a 500-run design that varies input parameters temperature (T), electron density (Ne), proportion of sodium (pNa), and proportion of copper (pCu). The latter two variables sum to one and are unused in the single-element simulations. The validation data was produced with a 25-run design. The coarse simulations produce spectra over a range of 240nm - 880nm that roughly mimics the range collected by the ChemCam instrument on the Mars rover Curiosity. The fine resolution simulations produce spectra over a range of 550nm to 600nm and were produced to mimic the included experimental results. The experimental results cover four mixtures of sodium and copper. All data sets are kept in zip files whose names indicate the elemental composition (NaCu, Na, or Cu), resolution (coarse or fine), and purpose (training or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. See attached abstract for figures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Leaching Kinetics of Rare Earth Elements from Fire Clay Seam Coal

Recovery of rare earth elements (REEs) from coal samples collected from the Fire Clay coal seam using diluted mineral acid solutions was investigated. The initial processing step was coal recovery using conventional froth flotation which concentrated the REEs in tailing material resulting in an upgrade to values around 700 ppm on a dry whole mass basis. Leaching experiments were performed on the flotation tailings material using a 1.2 M sulfuric acid solution adjusted to a temperature of 75 °C to study the extractability of REEs from coal material. The effect of particle size, leaching time, leaching temperature, and solid concentration on REE leaching recovery were evaluated. The kinetic data obtained from leaching over a range of temperatures suggested that the leaching process follows the shrinking core model with possibly a mixed control mechanism that may be a result of several heterogenous materials leaching simultaneously. Leaching recovery increased rapidly at the beginning of the reaction then slowed as the system reached equilibrium. The apparent activation energy determined from test data obtained over a range of temperatures using 1 M sulfuric acid was 36 kJ/mol for the first 20 min of reaction time and 27 kJ/mol for the leaching period between 20 and 120 min. The leaching of light REEs during the initial stage was determined to be driven by a chemical reaction, followed by the formation of a product layer, which required lower activation energy in the later stage of leaching. In regards to the heavy REEs, the major mechanism for leaching is desorption and the product layer formation does not affect the heavy REEs significantly.

58 GEOSCIENCES↗
DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

High-pressure neutron powder diffraction study of an arsenolite deuterium inclusion compound: structure and formation kinetics

High-pressure variable-temperature neutron powder diffraction was employed to investigate the synthesis and structural properties of arsenolite (As 4 O 6 ), a cubic arsenic( III ) oxide polymorph, and its inclusion compound with deuterium. The compound forms with the stoichiometry As 4 O 6 ·2D 2 , consistent with earlier findings. At 295 K and 1.96(2) GPa, D 2 molecules occupy the 16c site (0, 0, 0) and are aligned along the <111> direction. The diffraction data suggest that the D 2 molecules at this site become rotationally disordered upon decompression to 0.78(2) GPa at 120 K, while no occupancy is observed at the 8a position (⅛, ⅛, ⅛) under any of the studied conditions. The kinetic data of the high-pressure inclusion compound synthesis at 1.56(2) and 1.96(2) GPa were analyzed using the Avrami model. The Avrami exponent, 0.135(16), remains invariant with pressure, demonstrating that the reaction topochemistry is unaffected by compression. Similarly, the reaction rates, 0.29(3) and 0.23(3) min −0.135 at 1.56(2) and 1.96(2) GPa, respectively, exhibit negligible pressure dependence and suggest a surprisingly low activation volume of 1.6 cm 3 mol −1 for the reaction. These results provide detailed insight into the structural dynamics and kinetics of deuterium inclusion in arsenolite under high pressure.

Guńka, Piotr A. [Warsaw University of Technology (↗
DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

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

Metabolic Pathway Engineering

Modern microbial and enzyme engineering and their advancement are increasingly dependent on the marriage of a wide range of sophisticated technologies. For students entering the field of biotechnology, the outlook is indeed daunting. Expertise at levels beyond that of simple familiarity will be needed to conduct competitive research. It goes without saying that to be competitive, all new researchers in this field will need basic preparation in molecular biology, biochemistry, and genetics. Further, experience working with concepts and experimental tools in enzyme biochemistry and kinetics, gene editing, computational metabolic pathway modeling, experimental pathway flux analysis, and computational clustering tools to process complex data sets will be vital for success. We speculate that most biotechnology researchers in early career at this time will build teams of collaborators to address these disparate science fields rather than attempt to become experts in one lab.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗