SEARCH · Engineering Papers

Results for “Kinetic modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Turbulent hydrocarbon combustions kinetics - Stochastic modeling and verification

Idealized reactors, that are designed to ensure perfect mixing and are used to generate the combustion kinetics for complex hydrocarbon fuels, may depart from the ideal and influence the kinetics model performance. A complex hydrocarbon kinetics model that was established by modeling a jet-stirred combustor (JSC) as a perfectly stirred reactor (PSR), is reevaluated with a simple stochastic process in order to introduce the unmixedness effect quantitatively into the reactor system. It is shown that the comparisons of the predictions and experimental data have improved dramatically with the inclusion of the unmixedness effect in the rich combustion region. The complex hydrocarbon kinetics is therefore verified to be mixing effect free and be applicable to general reacting flow calculations.

Wang, T. S.↗

DOE OSTI2021

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.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

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.

36 MATERIALS SCIENCE↗

DOE OSTI2021

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.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2025

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.

autoencoder neural networks↗

DOE OSTI2025

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.

Chemical kinetics↗

DOE OSTI2024

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.

Rizzo, Thomas G. [SLAC National Accelerator Labora↗

DOE OSTI2021

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Isotherm, Kinetic, Process Modeling, and Techno-Economic Analysis of a Diamine-Appended Metal-Organic Framework for CO 2 Capture Using Fixed Bed Contactors

In this report, diamine-appended metal-organic frameworks exhibiting step-shaped CO 2 adsorption are exceptional candidates for energy-efficient carbon capture. However, there are few studies examining their performance in real-world capture scenarios, in part due to the challenge inherent in modeling their CO 2 uptake behavior. Here, we develop a dual-site Sips model to fit experimental CO 2 adsorption data for dmpn-Mg 2 (dobpdc) (dmpn = 2,2-dimethyl-1,3-diaminopropane; dobpdc 4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) and develop a linear driving force model for the adsorption kinetics based on available experimental data. These models are used to develop a dynamic, fixed bed, nonisothermal contactor model using shaped particles of the material, which is validated with experimental breakthrough data. We also examine the effects of the high heat of adsorption of the material on CO 2 uptake performance and find that heat removal is essential to maximize capture performance. We finally investigate "basic"(no bed cooling during adsorption) and "modified"(bed cooling during adsorption) temperature swing adsorption (TSA) processes using dmpn-Mg 2 (dobpdc), and their process economics are compared to a state-of-the-art monoethanolamine (MEA) capture system with and without heat recovery. In the absence of heat recovery, the adsorbent systems are more costly than established technology. However, with 85% heat recovery, both adsorbent-based TSA processes are projected to cost less than the MEA system. This work highlights that thermal management is vital for implementation of dmpn-Mg 2 (dobpdc) as a viable CO 2 capture technology. Investigation of other contactor technologies that can provide unique ways to manage system heat represent promising future areas of study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Toward UV models of kinetic mixing and portal matter. VII. A light dark photon in the 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 model

The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new 𝑈⁢(1) 𝐷 gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter 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, also exist. The presence of such types of states and the strong suggestion of a need to embed 𝑈⁢(1) 𝐷 into a non-Abelian gauge structure not too far above the TeV scale based on the renormalization-group equation running of the 𝑈⁢(1)𝐷 gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group 𝐺 = 𝑆⁢𝑈⁢(3) 𝑐 × 𝑆⁢𝑈⁢(3) 𝐿 × 𝑈⁢(1) 𝐴 × 𝑈⁢(1) 𝐵 = 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 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 𝐺 = 𝐺 SM × 𝐺 𝐷 encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability

This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.

petrochemicals↗

NASA NTRS1981

Effects of turbulence on a kinetic auroral arc model

A plasma kinetic model of an inverted-V auroral arc structure which includes the effects of electrostatic turbulence is proposed. In the absence of turbulence, a parallel potential drop is supported by magnetic mirror forces and charge quasi neutrality, with energetic auroral ions penetrating to low altitudes; relative to the electrons, the ions' pitch angle distribution is skewed toward smaller pitch angles. The electrons energized by the potential drop form a current which excites electrostatic turbulence. In equilibrium the plasma is marginally stable. The conventional anomalous resistivity contribution to the potential drop is very small. Anomalous resistivity processes are far too dissipative to be powered by auroral particles. It is concluded that under certain circumstances equilibrium may be impossible and relaxation oscillations set in.

Cornwall, J. M.↗

DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Thermodynamics and kinetics of precipitation and austenite reversion during aging of Ti-free grade 300 maraging steel manufactured by laser powder bed fusion (LPBF)

Austenite reversion and intermetallic precipitation are two key phase transformations in maraging steels which determine the strength and ductility. Usually, an optimum aging temperature/time is determined using trial and error approach with inputs from thermodynamic databases and a-posteriori microscopy/hardness evaluations. In this study, a thermo-kinetic model is developed to predict the precipitation and austenite reversion kinetics during various post-fabrication heat treatments in a Ti-free variant of grade 300 maraging steel manufactured using laser powder bed fusion. The thermo-kinetic model was based on the principles of Scheil solidification and the classical nucleation theory. The predictions from the thermo-kinetic model were verified using suitable microscopy techniques, and the results from the thermo-kinetic model agree well with the experimental measurements. Overall, results highlight the fact that an Integrated Computational Materials Engineering approach can be used to predict the precipitation and austenite reversion kinetics, and be used as an input for an efficient design of heat treatments for achieving optimum mechanical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Reaction Pathways and Energy Consumption in NH 3 Decomposition for H 2 Production by Low Temperature, Atmospheric Pressure Plasma

Pathways for NH 3 decomposition to N 2 and N 2 H 4 by atmospheric pressure nonthermal plasma are analyzed using a combination of molecular beam mass spectrometry measurements and zero-dimensional kinetic modeling. Experimental measurements show that NH 3 conversion and selectivity towards N 2 formation scale monotonically with the specific energy input into the plasma with ~ 100% selectivity to N 2 formation achieved at specific energy inputs above 0.12 J cm −3 (3.1 eV (molecule NH 3 ) −1 ). The kinetic model recovers these trends, although it underpredicts N 2 selectivity at low specific energy input. These discrepancies can be explained by the underestimation of reaction rate coefficients for reactions that consume N 2 H x species in collisions with H radicals and/or radial nonuniformities in power deposition, gas temperature, and species concentrations that are not represented by the plug flow approximation used in the model. The kinetic model shows that N 2 formation proceeds through N 2 H x decomposition pathways rather than NH x decomposition pathways in low temperature, atmospheric pressure plasma. Higher selectivity toward N 2 production can be achieved by operating at higher NH 3 conversion and with a higher gas temperature. Furthermore, the high energy cost of NH 3 decomposition by atmospheric pressure nonthermal plasma found in this work (25–50 eV (molecule NH 3 converted) −1 ; 17–33 eV (molecule H 2 formed) −1 ) is a result of the energy requirement for electron-impact dissociation of NH 3 and the significant re-formation of NH 3 by three-body recombination reactions between NH 2 and H.

Nonthermal plasma↗

DOE OSTI2025

Real-fluid behavior in rapid compression machines: Does it matter?

Rapid compression machines (RCMs) have been extensively used to quantify fuel autoignition chemistry and validate chemical kinetic models at high-pressure conditions. Historically, the analyses of experimental and modeling RCM autoignition data have been conducted based on the adiabatic core hypothesis with ideal gas assumption, where real-fluid behavior has been completely overlooked, though this might be significant at common RCM test conditions. Here, this work presents a first-of-its-kind study that addresses two significant but overlooked questions for autoignition studies within RCMs in the fundamental combustion community: (i) experiment-wise, can unaccounted-for real-fluid behavior in RCMs affect the interpretation and analysis of RCM experimental data? and (ii) simulation-wise, can unaccounted-for real-fluid behavior in RCMs affect RCM autoignition modeling and the validation of chemical kinetic models? To this end, theories for real-fluid isentropic change are newly proposed and derived based on high-order Virial EoS, and are further incorporated into an effective-volume real-fluid autoignition modeling framework newly developed for RCMs. With detailed analyses, the strong real-fluid behavior in representative RCM tests is confirmed, which can greatly influence the interpretation of RCM autoignition experiments, particularly the determination of end-of-compression temperature and evolution of the adiabatic core in the reaction chamber. Furthermore, real-fluid RCM modeling results reveal that considerable error can be introduced into simulating RCM autoignition experiments when following the community-wide accepted effective-volume approach by assuming ideal-gas behavior, which can be as high as 64% in the simulated ignition delay time at compressed pressure of 125 bar and lead to contradictory validation results of chemical kinetic models. Therefore, we recommend the community to adopt frameworks with real-fluid behavior fully accounted for (e.g., the one developed in this study) to analyze and simulate past and future RCM experiments, so as to avoid misinterpretation of RCM autoignition experiments and eliminate the potential errors that can be introduced into the simulation results with the existing RCM modeling frameworks.

High-order Virial equation of state↗

DOE OSTI2024

Forward and inverse predictive transient models of TREAT using surrogate reactivity models

In this work, we present two novel approaches for predicting the power evolution and control rod height of the Transient Reactor Test Facility (TREAT) to support experiment modeling; specifically, transient analysis of the NASA-sponsored Sirius series of experiments. These approaches utilize steady-state Monte Carlo model, point kinetics model, and surrogate models to predict power evolution and control rod axial position during transient experiments. Both approaches were tested and validated against several Sirius experiments that were performed in the TREAT facility at different power levels. The validation test results show very good agreement with the experimental data, and the models were able to accurately predict the power evolution and the axial control rod position with an average error within 3.0%. Finally, this indicates that these approaches will help the reactor engineering team of the TREAT facility in preparing and predicting the power and temperature of the experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

maranasgroup/yeastsMFA

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

Dinh, Hoang↗