Search NASA⌕ Search

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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 1

We report di-isobutylene has received significant attention as a promising fuel blendstock, as it can be synthesized via biological routes and is a short-listed molecule from the Co-Optima initiative. Di-isobutylene is also popularly used as an alkene representative in multi-component surrogate models for engine studies of gasoline fuels. However, there is limited experimental data available in the literature for neat di-isobutylene under engine-like conditions. Hence, most existing di-isobutylene models have not been extensively validated, particularly at lower temperatures (< 1000 K). Most gasoline surrogate models include the di-isobutylene sub-mechanism published by Metcalfe et al. with little or no modification. The current study is undertaken to develop a detailed kinetic model for di-isobutylene and validate the model using a wide range of relevant experimental data. Part 1 of this study exclusively focuses on the low- to intermediate temperature kinetics of di-isobutylene. An upcoming part 2 discusses the high-temperature model development and validation of the relevant experimental targets. Ignition delay time measurements for the di-isobutylene isomers were performed at pressures ranging from 15 - 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in air and in the temperature range 650 - 900 K using two independent rapid compression machine facilities. In addition, measurements of species identified during the oxidation of these isomers were performed in a jet-stirred reactor and in a rapid compression machine. A detailed kinetic model for the di-isobutylene isomers is developed to capture the wide range of new experimental targets. For the first time, a comprehensive low-temperature chemistry submodel is included. The differences in the important reaction pathways for the accurate prediction of the oxidation of the two DIB isomers are compared using reaction path analysis. In conclusion, the most sensitive reactions controlling the ignition delay times of the DIB isomers under the pressure and temperature conditions necessary for autoignition in engines are identified.

09 BIOMASS FUELS↗
DOE OSTI2020

The Sensitivity of Chemical Kinetic Models on Flame Transfer Functions in Acoustic Fluctuation Environments with Varying Equivalence Ratios

Various factors affect the heat release response of flames – thermodynamics, hydrodynamics, acoustics, etc. The coupling of heat release rate and acoustics, in particular, is essential to the design of combustion devices so as to avoid adverse effects during their operation. Modeling flames in acoustic environments requires an accurate understanding of the fundamental physics that occurs when flames interact with a fluctuating pressure field. Over the past couple of decades, researchers have been quantifying the impact of this fluctuating pressure field on the unsteady heat release rate. The flame transfer functions computed have also been shown to depend on the level of detail in the chemical kinetic model. This paper attempts to bridge the gap between these studies to understand how sensitive the flame transfer function is to various kinetic models. The problem is studied by performing multiple one-dimensional simulations using the low-Mach solver NGA. The detailed GRI-Mech 3.0 will be used for comparison. Also, the pressure amplitude and reactant equivalence ratio are varied to understand how the environment affects this sensitivity.

Girdhar, Anant↗
DOE OSTI2017

An innovative method to build a comprehensive kinetic model for air injection using TGA/DSC experiments

In order to study the AIP (air injection process), a practical method was developed to obtain the kinetic data for an AIP by matching the simulation model results with the TGA (thermogravimetry analysis) experiments in our previous work. However, one shortcoming for that method is the non-uniqueness of the kinetic data obtained. The other shortcoming is that TGA cannot measure the exothermic characteristics of a crude oil. The enthalpy values of the reactions are crucial for developing a comprehensive kinetic model as the enthalpy represents the heat generation of the chemical reaction and dominants the heat effect during the AIP. Here in this study, an innovative method is proposed to build a comprehensive kinetic model by combining both TGA/DSC (differential scanning calorimetry) experiments and simulation. In the low temperature interval, the distillation dominant stage is differentiated from the LTO (low temperature oxidation) stage by comparing the nitrogen purging TGA experiments and the air purging TGA experiments. Then, the LTO stage is further divided into several sub-divided reaction stages with corresponding temperature intervals by analyzing the TGA/DSC data. The kinetic data are estimated by applying the Arrhenius method to sub-divided reaction stages of the TGA experiments. In order to obtain the enthalpy values for oxidation reactions, the negative heat flow of the distillation process is eliminated by subtracting the nitrogen purging DSC data from the air purging DSC data, and the enthalpy values of corresponding reaction stages are obtained from the subtracted heat flow curve whose area of the exothermic peak represents the enthalpy value of a stage. Moreover, the pseudo components are defined based on the nitrogen purging TGA experiments which are associated with the sub-divided reaction stages defined previously to complete the kinetic model. Finally, the kinetic model is validated with the air purging TGA experiments. This study shows a detailed workflow to build a comprehensive kinetic model for AIP by TGA/DSC experiments. An application to a crude oil is also presented to demonstrate its practicability.

42 ENGINEERING↗
DOE OSTI2025

Challenges in integrating dissolved organic matter chemodiversity into kinetic models of soil respiration

The chemodiversity of dissolved organic matter (DOM) in soil has been proposed to influence the microbial metabolism and fate of belowground organic carbon (C). However, integrating DOM chemistry into soil C cycle models to improve predictions of C stocks and fluxes—beyond simply considering DOM pool size—remains a challenge. While recent research suggests that incorporating DOM chemodiversity into models can improve predictions of microbial respiration, there is still a lack of mechanistic understanding describing how DOM chemodiversity affects microbial metabolism and soil respiration. Here, we evaluated whether DOM chemodiversity was a determinant of soil respiration using paired measurements of high-resolution DOM chemistry, obtained from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), and potential soil respiration rates from across the United States (U.S.), all data provided by the Molecular Observation Network. Our objectives were to (1) assess statistical relationships between DOM chemodiversity and microbial respiration, and (2) evaluate the ability of kinetic models to leverage DOM chemistry to explain empirical relationships found in statistical models. Statistical regressions revealed that DOM chemodiversity (alpha diversity) was nonlinearly related to potential soil respiration rates, both independently and through its interactions with DOM and total C concentrations. In soils with relatively high DOM but low total C concentrations, potential soil respiration rates were negatively correlated with DOM alpha diversity, whereas in soils with relatively low DOM and high total C concentrations showed the opposite trend. However, when metabolic transition theory kinetic models were modified to include chemodiversity, their performance was comparable to traditional Monod kinetics approaches, which simulate respiration rates as a function of DOM concentration. The inability to account for nonlinearities in DOM chemodiversity–respiration relationships highlight an opportunity to advance substrate uptake kinetics by establishing causal links between DOM chemodiversity, microbial metabolism trade-offs, and potential interactions under varied environmental conditions.

Bioenergetic model↗
DOE OSTI2020

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

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

Kinetic Modeling of the Reaction Rate for Quartz and Carbon Black Pellet

The kinetic modeling for the carbothermal reduction reaction rate in quartz and carbon black pellets is studied at different temperatures, under varying CO partial pressures in ambient atmosphere, varying carbon contents, different quartz particle sizes, and different crucible opening areas. Carbon black is produced by the cracking of natural gas. The activation energy of the SiC-producing step was determined to be 594 kJ/mol. The averaged pre-exponential factor A obtained from 1898 K, 1923 K, and 1948 K (1625 °C, 1650 °C, and 1675 °C) is 2.62E+16 min -1 . The reaction rate of the gas-solid interface factor, fix-C content (X fix-C ), temperature (T), and CO partial pressure (X CO ) can be expressed as follows:((d pct)/(dt))=(1-0.40×X fix-C -0.86 ×pct)×2.62×10 16 ×exp(-(594000/(RT)))×(2.6-0.015×X co ).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Approach to the Steady State in Kinetic Models with Thermal Reservoirs at Different Temperatures

We continue the investigation of kinetic models of a system in contact via stochastic interactions with several spatially homogeneous thermal reservoirs at different temperatures. Considering models different from those investigated in Carlen et al. (Braz J Probab Stat 29:372–386, 2015), we explicitly compute the unique spatially uniform non-equilibrium steady state (NESS) and prove that it is approached exponentially fast from any uniform initial state. This leaves open the question of whether there exist NESS that are not spatially uniform. Making a further simplification of our models, we then prove non-existence of such NESS and exponential approach to the unique spatially uniform NESS (with a computably boundable rate). The method of proof relies on refined Doeblin estimates and other probabilistic techniques, and is quite different form the analysis in Carlen et al. (Braz J Probab Stat 29:372–386, 2015) that was based on contraction mapping methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Development of an Austenitization Kinetics Model for 22MnB5 Steel

This paper presents a first-order austenitization kinetics model for 22MnB5 steel, commonly used in hot forming die quenching. Model parameters are derived from constant heating rate dilatometry measurements. Vickers hardness measurements made on coupons that were quenched at intermediate stages of the process were used to verify the model, and the Ac{sub 1} and Ac{sub 3} temperatures inferred from dilatometry are consistent with correlations found in the literature. The austenitization model was extended to consider non-constant heating rates typical of industrial furnaces and again showed reasonable agreement between predictions and measurements. Finally, the model is used to predict latent heat evolution during industrial heating and is shown to be consistent with values inferred from thermocouple measurements of furnace-heated 22MnB5 coupons reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Three-dimensional kinetic modeling of streamer propagation in a nitrogen/helium gas mixture

A fully resolved kinetic model (Particle-in-Cell and Direct Simulation Monte Carlo for particle/photon collisions) of a near atmospheric pressure ionization wave is presented here. Fully resolving the required numerical spatial (sub-um) and temporal scales (tens of fs) for atmospheric pressure discharges in three dimensions is still a challenging task on modern super computers. To keep the overall problem tractable, the total number of elements are reduced by only simulating a 10 degree wedge rather than a full 360 degree geometry. The ionization wave is generated in a needle-plane configuration with a gap size of 250 um and a background of nitrogen and helium gas. A voltage of 1500 V is applied to the anode and an initial electron and ion density of 10 9 cm 3 is seeded in a region near the anode electrode tip and extending towards the cathode. As these initial electrons are swept away, photoionization and photoemission create new electrons and allow the ionization front to propagate towards the cathode. Results from the 90% N 2 , 10% He discharge indicate that photoionization has minimal impact on plasma formation processes and cathode photoemission is the dominant mechanism for new electrons. In the 90% He, 10% N 2 discharge case, however, photoionization likely has an impact as the observed locations of photoionization occur far enough away from the ionization front to allow for sufficient avalanche processes that contribute to the propagation of the ionization wave. Additionally, the electron energy distribution functions in the 90% He, 10% N 2 case indicate that there is less energy loss to the low lying molecular N 2 electronic states as well as the vibrational and rotational modes. This leads to higher electron energies and faster plasma development times of ~0.4 ns for the 90% He, 10% N 2 case, and ~1.5 ns for the 90% N 2 , 10% He case. In addition to analysis of the ionization wave results, the overall challenges associated with simulations near atmospheric pressure discharges in three dimensions are discussed, including the limitations of the 10 degree wedge that produces, at least qualitatively, minimal 3D effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Kinetic modeling of ignition in miniature shock tube

The ignition delay times of toluene, n-heptane/iso-octane primary reference fuel blends, and toluene/n-heptane/iso-octane toluene primary reference fuel blends were studied in the miniature high repetition rate shock tube (HRRST). In the HRRST, the temperature and pressure change during the test time, complicating interpretation of ignition delay times. Two methods of chemical kinetic modeling were employed in order to account for the changing thermodynamic states in the HRRST. One was based on direct modeling using the pressure profile and the kinetic mechanism, and the other was based on a Livengood-Wu type approach using modeled constant volume ignition delay times. Two different kinetic mechanisms were used in modeling the data. Both methods showed agreement with experimental data over a range of temperatures, pressures, and blends. While overall agreement was fair, agreement was poorer for one of the mechanisms for the primary reference fuel blends. An evaluation of the ignition delay iso-contours from these mechanisms for representative HRRST experiment thermodynamic trajectories revealed that the HRRST measured IDTs are controlled by the high temperature chemistry and are more sensitive to temperature. Furthermore, the different performance of mechanisms is therefore attributed to their different high temperature reactivity.

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

Kinetic model development for single step ethanol to butene rich olefin process over Cu-Y/Beta catalysts

Here, this study presents the first intrinsic kinetic model for the single-step conversion of ethanol to butene-rich olefins over bifunctional Cu-Y/Beta catalysts, addressing a critical gap in the design and scale-up of Sustainable Aviation Fuel (SAF) processes. The reaction network comprises ten global steps involving dehydrogenation, aldol condensation, hydrogenation, and dehydration reactions, distributed across Cu and rare-earth (Y) active sites. The model incorporates dual-site functionality (Cu and Y site) and explicitly accounts for key intermediates such as crotonaldehyde and butanal. Reaction rates are formulated using Langmuir–Hinshelwood–Hougen–Watson (LHHW) kinetics. Kinetic parameters are extracted by fitting the model to lab-scale packed-bed reactor data across a wide range of temperatures and space velocities, demonstrating strong agreement in ethanol conversion and product selectivity. The reaction kinetics developed in this work provide a foundational basis for constructing reactor models that enable process optimization and scale-up of ethanol-to-jet fuel technologies.

Cu-Y/Beta catalyst↗
DOE OSTI2024

An improved semi-global intrinsic kinetics model for high temperature carbon oxidation

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Experimental and kinetic modeling study of the low- temperature and high-pressure combustion chemistry of straight chain pentanol isomers: 1-, 2- and 3-Pentanol

Pentanols have received significant attention as a potential alternative fuel or fuel additive owing to their high energy densities and low vapor pressure. The development of robust chemical kinetic models for alternative fuels which can provide accurate and efficient predictions of combustion performance across a wide range of engine relevant conditions is important in developing cleaner, more efficient combustors. Although the high temperature oxidation kinetics of pentanol isomers has been researched considerably, their low temperature combustion chemistry needs further investigation. While previously proposed low temperature mechanisms for 1-pentanol based on analogy and rate rules need further refinement, the low temperature oxidation kinetics of 2-pentanol and 3-pentanol has not been studied previously by any means, experimentally or theoretically. A newly developed kinetic mechanism is presented in this work for the three straight chain pentanol isomers: 1-, 2- and 3-pentanol. Low temperature kinetics is based on a recent study by Lockwood et al., 2022 [20] involving theoretical calculations at the CCSD(T)/cc-pV∞Z level of theory for the oxidation pathways involving alcohol peroxy radicals. Rate of production analyses performed in this study highlight the importance of the newly added pressure-dependent reactions of the α-alcohol peroxy radical forming an RO 2 adduct. While the α-alcohol fuel radical reacts with O 2 to directly decompose via a chemically activated pathway at low pressures, the formation of the RO 2 adduct is favored at high pressures. Furthermore, the detailed model is comprehensively validated against new ignition experiments at low temperature and high pressure, together with the wide range of data available in the literature. Both qualitative and quantitative predictions of the experimental data using the proposed kinetic model are satisfactory for all three pentanol isomers studied here.

1-pentanol↗
DOE OSTI2019

Kinetic modeling study of surrogate components for gasoline, jet and diesel fuels: C7-C11 methylated aromatics

Kinetic mechanisms for aromatics are needed to successfully simulate the autoignition of transportation fuels using the surrogate fuel approach. An aromatic detailed kinetic mechanism that describes kinetics of C7-C11 methylated aromatics, including toluene, o-xylene, p-xylene, 1,2,4-trimethylbenzene, 1,3,5, trimethyl-benzene and alpha-methylnaphthalene has been developed in the current study. The kinetic mechanism was built hierarchically using similar set of reaction pathways and reaction rate rules. In the mechanism developed, special emphasis is put on describing the detailed low-temperature ignition chemistry of o-xylene and 1,2,4-trimethylbenzene and, to our knowledge, this is the first attempt to do so in a detailed kinetic mechanism. In addition to kinetic modeling, new experimental data were acquired for toluene, o-xylene, and 1,2,4-trimethylbenzene using a rapid compression machine at low-to-intermediate temperatures and engine relevant pressures. In addition, the mechanism has been compared against data sets from the literature covering ignition delay times, flame speeds, and speciation profiles measured in a jet-stirred reactor and flow reactor. Good agreement is observed between the mechanism predictions and the experimental data. Kinetic analysis demonstrated the importance of including the low temperature chemistry of the benzylperoxy radicals to accurately capture the ignition propensity of o-xylene and 1,2,4-trimethylbenzene at low-to-intermediate temperatures and high pressures. The kinetic mechanism developed in the current study can be used for surrogate modeling of gasoline, jet and diesel fuels. (C) 2018 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

02 PETROLEUM↗
DOE OSTI2019

Validation of a Simple RELAP5-3D Point Kinetics Model of the Full-Slotted MARCH Core in TREAT

A simple RELAP5-3D point kinetics model of the TREAT facility was validated using data from three unshaped transients in the MARCH full-slotted core. These transients were initiated with near step insertions of reactivity of about 1.7%, 2.9%, and 3.9%, respectively. The RELAP5-3D calculations were in reasonable agreement with the measured results. The calculated energy depositions were within the scatter of the data or were slightly conservative for all three transients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectral analysis and kinetic modeling of radioluminescence in air and nitrogen

Here, in this article we present a quantitative analysis of the second positive system of molecular nitrogen and the first negative system of the molecular nitrogen cation excited in the presence of ionizing radiation. Optical emission spectra of atmospheric air and nitrogen surrounding 210 Po sources were measured from 250 to 400 nm. Multi-Boltzmann and non-Boltzmann vibrational distribution spectral models were used to determine the vibrational temperature and vibrational distribution function of the emitting N 2 (C 3 Π u ) and N 2 + (B 2 Σ + u ) states. A zero-dimensional kinetic model, based on the electron energy distribution function (EEDF) and steady-state excitation and de-excitation of N 2 (X 1 Σ + g ), N 2 + (B 2 Σ + u ), N 2 + (X 2 Σ + g ), N 4 + , O 2 + , and N 2 (C 3 Π u , v), was developed for the prediction of the relative spectral intensity of both the N 2 + (B 2 Σ + u → X 2 Σ + g ) emission band and the vibrational bands of N 2 (C 3 Π u → B 3 Π g ) for comparison with the experimental data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2021

Kinetic modeling of anaerobic degradation of plant-derived aromatic mixtures by Rhodopseudomonas palustris

Abstract Rhodopseudomonas palustris is a model microorganism for studying the anaerobic metabolism of aromatic compounds. While it is well documented which aromatics can serve as sole organic carbon sources, co-metabolism of other aromatics is poorly understood. This study used kinetic modeling to analyze the simultaneous degradation of aromatic compounds present in corn stover hydrolysates and model the co-metabolism of aromatics not known to support growth of R. palustris as sole organic substrates. The simulation predicted that p -coumaroyl amide and feruloyl amide were hydrolyzed to p -coumaric acid and ferulic acid, respectively, and further transformed via p -coumaroyl-CoA and feruloyl-CoA. The modeling also suggested that metabolism of p -hydroxyphenyl aromatics was slowed by substrate inhibition, whereas the transformation of guaiacyl aromatics was inhibited by their p -hydroxyphenyl counterparts. It also predicted that substrate channeling may occur during degradation of p -coumaroyl-CoA and feruloyl-CoA, resulting in no detectable accumulation of p -hydroxybenzaldehyde and vanillin, during the transformation of these CoA ligated compounds to p- hydroxybenzoic acid and vanillic acid, respectively. While the simulation correctly represented the known transformation of p -hydroxybenzoic acid via the benzoyl-CoA pathway, it also suggested co-metabolism of vanillic acid and syringic acid, which are known not to serve as photoheterotrophic growth substrate for R. palustris .

59 BASIC BIOLOGICAL SCIENCES↗