SEARCH · Search NASA
Results for “Kinetic Modelling”
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.
A standardized workflow for kinetic metabolic model curation and dissemination
Kinetic metabolic models provide invaluable insights into cellular metabolism, supporting applications in synthetic biology, metabolic engineering, and systems biology. However, reproducibility and utility of these models hinge on clear and rigorous documentation, standardized annotation, and accessible visualization. This paper presents a workflow for building, annotating, visualizing, and sharing kinetic metabolic models. Our method integrates community standards and open-source tools to ensure reproducibility, interoperability, and user accessibility. This procedure enables researchers to produce reusable and well-documented kinetic models, advancing their role as powerful tools in metabolic research.
Low temperature oxidation of methylpropyl ether
Ethers are potential high-performance fuels that can be produced from renewable carbon sources such as biomass, but their combustion properties are not clearly understood. The cetane numbers (CN) for ethers, a measure of autoignition tendency, can range from less than 10 to greater than 100, and predictive tools based upon the molecular structures are desirable in order to focus fuel development research. In this study, we used molecular modeling, chemical kinetic modeling, and experimental measurements to understand the important reactions that lead to the autoignition of methylpropyl ether (MPE), a model for ethers containing alkyl chains long enough for intramolecular hydrogen transfer. These reactions are known to be important in the autoignition of alkanes and ethers and we propose to investigate the relationship between alkyl and ether reaction pathways. A reaction mechanism for MPE was created using MIT’s Reaction Mechanism Generator (RMG) and G4 quantum calculation results and models were tested using the ChemKin Pro suite of software. Results from atmospheric pressure, flow tube reactor experiments and rapid compression machine (RCM) experiments were used to test the kinetic model and suggest refinements.
A further experimental and modeling study of acetaldehyde combustion kinetics
Acetaldehyde is an important intermediate and a toxic emission in the combustion of fuels, especially for biofuels. To better understand its combustion characteristics, a detailed chemical kinetic model describing the oxidation of acetaldehyde has been developed and comprehensively validated against various types of literature data including laminar flame speeds, oxidation and pyrolysis in shock tubes, chemical structure of premixed flames, and low-temperature oxidation in jet-stirred reactors. To extend the validation range, the chemical structure of a counterflow flame fueled by acetaldehyde at 600 Torr has been measured using vacuum ultra-violet photoionization molecular-beam mass spectrometry. In addition, ignition delay times at 10 atm and 700-1100 K were measured in a rapid compression machine, and a negative temperature coefficient (NTC) behavior was observed. Here, the present kinetic model well reproduces the results of various acetaldehyde combustion experiments covering wide ranges of temperatures (300–2300 K) and pressures (0.02–10 atm), and explains well the observed NTC behavior based on the competition between multiple oxidation pathways for the methyl radicals and their self-recombination forming ethane, a relatively stable species at temperatures below 1000 K.
Spectral method for a kinetic swarming model
Here we present the first numerical method for a kinetic description of the Vicsek swarming model. The kinetic model poses a unique challenge, as there is a distribution dependent collision invariant to satisfy when computing the interaction term. We use a spectral representation linked with a discrete constrained optimization to compute these interactions. To test the numerical scheme we investigate the kinetic model at different scales and compare the solution with the microscopic and macroscopic descriptions of the Vicsek model. Lastly, we observe that the kinetic model captures key features such as vortex formation and traveling waves.
Dynamic Metabolic Model Building Based on the Ensemble Modeling Approach
Ensemble modeling of kinetic systems addresses the challenges of kinetic model construction, with respect to parameter value selection, and still allows for the rich insights possible from kinetic models. This project aimed to show that constructing, implementing, and analyzing such models is a useful tool for the metabolic engineering toolkit, and that they can result in actionable insights from models. Key concepts are developed and deliverable publications and results are presented.
Low Temperature Oxidation of Methylpropyl Ether
Ethers are potential high-performance fuels that can be produced from renewable carbon sources such as biomass, but their combustion properties are not clearly understood. The cetane numbers (CN) for ethers, a measure of autoignition tendency, can range from less than 10 to greater than 100, and predictive tools based upon the molecular structures are desirable in order to focus fuel development research. In this study, we used molecular modeling, chemical kinetic modeling, and experimental measurements to understand the important reactions that lead to the autoignition of methylpropyl ether (MPE), a model for ethers containing alkyl chains long enough for intramolecular hydrogen transfer. These reactions are known to be important in the autoignition of alkanes and ethers and we propose to investigate the relationship between alkyl and ether reaction pathways. A reaction mechanism for MPE was created using MIT's Reaction Mechanism Generator (RMG) and G4 quantum calculation results and models were tested using the ChemKin Pro suite of software. Results from atmospheric pressure, flow tube reactor experiments and rapid compression machine (RCM) experiments were used to test the kinetic model and suggest refinements.
A Reduced Resistive Wall Mode Kinetic Stability Model for Disruption Forecasting
Kinetic modification of ideal stability theory from stabilizing resonances of mode-particle interaction has had success in explaining resistive wall mode (RWM) stability limits in tokamaks. With the goal of real-time stability forecasting, a reduced kinetic stability model has been implemented in the new Disruption Event Characterization and Forecasting (DECAF) code, which has been written to analyze disruptions in tokamaks. The reduced model incorporates parameterized models for ideal limits on beta, a ratio of plasma pressure to magnetic pressure, which are shown to be in good agreement with DCON code calculations. Increased beta between these ideal limits causes a shift in the unstable region of delta W_K space, where delta W_K is the change in potential energy due to kinetic effects that is solved for by the reduced model, such that it is possible for plasmas to be unstable at intermediate beta but stable at higher beta. Gaussian functions for delta W_K are defined as functions of E cross B frequency and collisionality, with parameters reflecting the experience of the National Spherical Torus Experiment (NSTX). The reduced model was tested on a database of discharges from NSTX and experimentally stable and unstable discharges were separated noticeably on a stability map in E cross B frequency, collisionality space. The reduced model only failed to predict an unstable RWM in 15.6% of cases with an experimentally unstable RWM and performed well on predicting stability for experimentally stable discharges as well.
A genome-scale Escherichia coli kinetic metabolic model k-ecoli457 satisfying flux data for multiple mutant strains
Kinetic models of metabolism at a genome scale that faithfully recapitulate the effect of multiple genetic interventions would be transformative in our ability to reliably design novel overproducing microbial strains. Here, we introduce k-ecoli457, a genome-scale kinetic model of Escherichia coli metabolism that satisfies fluxomic data for wild-type and 25 mutant strains under different substrates and growth conditions. The k-ecoli457 model contains 457 model reactions, 337 metabolites and 295 substrate-level regulatory interactions. Parameterization is carried out using a genetic algorithm by simultaneously imposing all available fluxomic data (about 30 measured fluxes per mutant). Furthermore, the Pearson correlation coefficient between experimental data and predicted product yields for 320 engineered strains spanning 24 product metabolites is 0.84. This is substantially higher than that using flux balance analysis, minimization of metabolic adjustment or maximization of product yield exhibiting systematic errors with correlation coefficients of, respectively, 0.18, 0.37 and 0.47.
Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model
A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.
Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth
Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.
Uranium hydride corrosion. II Modeling spatially random growth using reaction and nucleation kinetics
A uranium hydride (UH 3 ) corrosion kinetics model has been developed, based on dimensionally restricted nucleation kinetics, which combines heretofore separate reaction and nucleation kinetics models into a single predictive construct. This theoretical framework accounts for the superficial and random pitting behaviour of UH 3 corrosion. The model can be used to generate the net corrosion surface area and volume fractions for a uranium specimen, generally accounting for its geometric form factor through the surface-to-volume ratio. Furthermore, the combined corrosion kinetics model enables the generation of quantities which may be directly compared to those obtained by non-destructive, non-contact evaluation methods.
Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling
In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.
A comprehensive iso-octane combustion model with improved thermochemistry and chemical kinetics
Iso-Octane (2,2,4-trimethylpentane) is a primary reference fuel and an important component of gasoline fuels. Furthermore, it is a key component used in surrogates to study the ignition and burning characteristics of gasoline fuels. This paper presents an updated chemical kinetic model for iso-octane combustion. Specifically, the thermodynamic data and reaction kinetics of iso-octane have been re-assessed based on new thermodynamic group values and recently evaluated rate coefficients from the literature. The adopted rate coefficients were either experimentally measured or determined by analogy to theoretically calculated values. New alternative isomerization pathways for peroxy-alkyl hydroperoxide ($\dot{O}$OQOOH) radicals were added to the reaction mechanism. The updated kinetic model was compared against new ignition delay data measured in rapid compression machines (RCM) and a high-pressure shock tube. Our experiments were conducted at pressures of 20 and 40 atm, at equivalence ratios of 0.4 and 1.0, and at temperatures in the range of 632–1060 K. The updated model was further compared against shock tube ignition delay times, jet-stirred reactor oxidation speciation data, premixed laminar flame speeds, counterflow diffusion flame ignition, and shock tube pyrolysis speciation data available in the literature. Finally, the updated model was used to investigate the importance of alternative isomerization pathways in the low temperature oxidation of highly branched alkanes. When compared to available models in the literature, the present model represents the current state-of-the-art in fundamental thermochemistry and reaction kinetics of iso-octane; and thus provides the best prediction of wide ranging experimental data and fundamental insights into iso-octane combustion chemistry.
Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks
Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.
Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite
This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.
Towards UV-Models of Kinetic Mixing and Portal Matter VII: A Light Dark Photon in the 3c3L1A1B Model
The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new U(1)DU(1)D gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter (DM) with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter (PM), also exist. The presence of such types of states and the strong suggestion of a need to embed U(1)DU(1)D into a non-abelian gauge structure not too far above the TeV scale based on the RGE running of the U(1)DU(1)D gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group G=SU(3)c×SU(3)L×U(1)A×U(1)B=3c3L1A1BG=SU(3)c×SU(3)L×U(1)A×U(1)B=3c3L1A1B is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-abelian fashion that is not a simple product group of the form G=GSM×GDG=GSM×GD encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.