Science.gov

Sample records for chemical kinetic modeling

  1. LLNL Chemical Kinetics Modeling Group

    SciTech Connect

    Pitz, W J; Westbrook, C K; Mehl, M; Herbinet, O; Curran, H J; Silke, E J

    2008-09-24

    The LLNL chemical kinetics modeling group has been responsible for much progress in the development of chemical kinetic models for practical fuels. The group began its work in the early 1970s, developing chemical kinetic models for methane, ethane, ethanol and halogenated inhibitors. Most recently, it has been developing chemical kinetic models for large n-alkanes, cycloalkanes, hexenes, and large methyl esters. These component models are needed to represent gasoline, diesel, jet, and oil-sand-derived fuels.

  2. Chemical kinetics modeling

    SciTech Connect

    Westbrook, C.K.; Pitz, W.J.

    1993-12-01

    This project emphasizes numerical modeling of chemical kinetics of combustion, including applications in both practical combustion systems and in controlled laboratory experiments. Elementary reaction rate parameters are combined into mechanisms which then describe the overall reaction of the fuels being studied. Detailed sensitivity analyses are used to identify those reaction rates and product species distributions to which the results are most sensitive and therefore warrant the greatest attention from other experimental and theoretical research programs. Experimental data from a variety of environments are combined together to validate the reaction mechanisms, including results from laminar flames, shock tubes, flow systems, detonations, and even internal combustion engines.

  3. Chemical kinetics and combustion modeling

    SciTech Connect

    Miller, J.A.

    1993-12-01

    The goal of this program is to gain qualitative insight into how pollutants are formed in combustion systems and to develop quantitative mathematical models to predict their formation rates. The approach is an integrated one, combining low-pressure flame experiments, chemical kinetics modeling, theory, and kinetics experiments to gain as clear a picture as possible of the process in question. These efforts are focused on problems involved with the nitrogen chemistry of combustion systems and on the formation of soot and PAH in flames.

  4. Chemical Kinetic Modeling of Advanced Transportation Fuels

    SciTech Connect

    PItz, W J; Westbrook, C K; Herbinet, O

    2009-01-20

    Development of detailed chemical kinetic models for advanced petroleum-based and nonpetroleum based fuels is a difficult challenge because of the hundreds to thousands of different components in these fuels and because some of these fuels contain components that have not been considered in the past. It is important to develop detailed chemical kinetic models for these fuels since the models can be put into engine simulation codes used for optimizing engine design for maximum efficiency and minimal pollutant emissions. For example, these chemistry-enabled engine codes can be used to optimize combustion chamber shape and fuel injection timing. They also allow insight into how the composition of advanced petroleum-based and non-petroleum based fuels affect engine performance characteristics. Additionally, chemical kinetic models can be used separately to interpret important in-cylinder experimental data and gain insight into advanced engine combustion processes such as HCCI and lean burn engines. The objectives are: (1) Develop detailed chemical kinetic reaction models for components of advanced petroleum-based and non-petroleum based fuels. These fuels models include components from vegetable-oil-derived biodiesel, oil-sand derived fuel, alcohol fuels and other advanced bio-based and alternative fuels. (2) Develop detailed chemical kinetic reaction models for mixtures of non-petroleum and petroleum-based components to represent real fuels and lead to efficient reduced combustion models needed for engine modeling codes. (3) Characterize the role of fuel composition on efficiency and pollutant emissions from practical automotive engines.

  5. Chemical kinetics models for semiconductor processing

    SciTech Connect

    Coltrin, M.E.; Creighton, J.R.; Meeks, E.; Grcar, J.F.; Houf, W.G.; Kee, R.J.

    1997-12-31

    Chemical reactions in the gas-phase and on surfaces are important in the deposition and etching of materials for microelectronic applications. A general software framework for describing homogeneous and heterogeneous reaction kinetics utilizing the Chemkin suite of codes is presented. Experimental, theoretical and modeling approaches to developing chemical reaction mechanisms are discussed. A number of TCAD application modules for simulating the chemically reacting flow in deposition and etching reactors have been developed and are also described.

  6. Chemical Kinetic Modeling of Hydrogen Combustion Limits

    SciTech Connect

    Pitz, W J; Westbrook, C K

    2008-04-02

    A detailed chemical kinetic model is used to explore the flammability and detonability of hydrogen mixtures. In the case of flammability, a detailed chemical kinetic mechanism for hydrogen is coupled to the CHEMKIN Premix code to compute premixed, laminar flame speeds. The detailed chemical kinetic model reproduces flame speeds in the literature over a range of equivalence ratios, pressures and reactant temperatures. A series of calculation were performed to assess the key parameters determining the flammability of hydrogen mixtures. Increased reactant temperature was found to greatly increase the flame speed and the flammability of the mixture. The effect of added diluents was assessed. Addition of water and carbon dioxide were found to reduce the flame speed and thus the flammability of a hydrogen mixture approximately equally well and much more than the addition of nitrogen. The detailed chemical kinetic model was used to explore the detonability of hydrogen mixtures. A Zeldovich-von Neumann-Doring (ZND) detonation model coupled with detailed chemical kinetics was used to model the detonation. The effectiveness on different diluents was assessed in reducing the detonability of a hydrogen mixture. Carbon dioxide was found to be most effective in reducing the detonability followed by water and nitrogen. The chemical action of chemical inhibitors on reducing the flammability of hydrogen mixtures is discussed. Bromine and organophosphorus inhibitors act through catalytic cycles that recombine H and OH radicals in the flame. The reduction in H and OH radicals reduces chain branching in the flame through the H + O{sub 2} = OH + O chain branching reaction. The reduction in chain branching and radical production reduces the flame speed and thus the flammability of the hydrogen mixture.

  7. Chemical kinetics and modeling of planetary atmospheres

    NASA Technical Reports Server (NTRS)

    Yung, Yuk L.

    1990-01-01

    A unified overview is presented for chemical kinetics and chemical modeling in planetary atmospheres. The recent major advances in the understanding of the chemistry of the terrestrial atmosphere make the study of planets more interesting and relevant. A deeper understanding suggests that the important chemical cycles have a universal character that connects the different planets and ultimately link together the origin and evolution of the solar system. The completeness (or incompleteness) of the data base for chemical kinetics in planetary atmospheres will always be judged by comparison with that for the terrestrial atmosphere. In the latter case, the chemistry of H, O, N, and Cl species is well understood. S chemistry is poorly understood. In the atmospheres of Jovian planets and Titan, the C-H chemistry of simple species (containing 2 or less C atoms) is fairly well understood. The chemistry of higher hydrocarbons and the C-N, P-N chemistry is much less understood. In the atmosphere of Venus, the dominant chemistry is that of chlorine and sulfur, and very little is known about C1-S coupled chemistry. A new frontier for chemical kinetics both in the Earth and planetary atmospheres is the study of heterogeneous reactions. The formation of the ozone hole on Earth, the ubiquitous photochemical haze on Venus and in the Jovian planets and Titan all testify to the importance of heterogeneous reactions. It remains a challenge to connect the gas phase chemistry to the production of aerosols.

  8. A kinetic model for chemical neurotransmission

    NASA Astrophysics Data System (ADS)

    Ramirez-Santiago, Guillermo; Martinez-Valencia, Alejandro; Fernandez de Miguel, Francisco

    Recent experimental observations in presynaptic terminals at the neuromuscular junction indicate that there are stereotyped patterns of cooperativeness in the fusion of adjacent vesicles. That is, a vesicle in hemifusion process appears on the side of a fused vesicle and which is followed by another vesicle in a priming state while the next one is in a docking state. In this talk we present a kinetic model for this morphological pattern in which each vesicle state previous to the exocytosis is represented by a kinetic state. This chain states kinetic model can be analyzed by means of a Master equation whose solution is simulated with the stochastic Gillespie algorithm. With this approach we have reproduced the responses to the basal release in the absence of stimulation evoked by the electrical activity and the phenomena of facilitation and depression of neuromuscular synapses. This model offers new perspectives to understand the underlying phenomena in chemical neurotransmission based on molecular interactions that result in the cooperativity between vesicles during neurotransmitter release. DGAPA Grants IN118410 and IN200914 and Conacyt Grant 130031.

  9. Chemical Kinetic Modeling of Biofuel Combustion

    NASA Astrophysics Data System (ADS)

    Sarathy, Subram Maniam

    Bioalcohols, such as bioethanol and biobutanol, are suitable replacements for gasoline, while biodiesel can replace petroleum diesel. Improving biofuel engine performance requires understanding its fundamental combustion properties and the pathways of combustion. This study's contribution is experimentally validated chemical kinetic combustion mechanisms for biobutanol and biodiesel. Fundamental combustion data and chemical kinetic mechanisms are presented and discussed to improve our understanding of biofuel combustion. The net environmental impact of biobutanol (i.e., n-butanol) has not been studied extensively, so this study first assesses the sustainability of n-butanol derived from corn. The results indicate that technical advances in fuel production are required before commercializing biobutanol. The primary contribution of this research is new experimental data and a novel chemical kinetic mechanism for n-butanol combustion. The results indicate that under the given experimental conditions, n-butanol is consumed primarily via abstraction of hydrogen atoms to produce fuel radical molecules, which subsequently decompose to smaller hydrocarbon and oxygenated species. The hydroxyl moiety in n-butanol results in the direct production of the oxygenated species such as butanal, acetaldehyde, and formaldehyde. The formation of these compounds sequesters carbon from forming soot precursors, but they may introduce other adverse environmental and health effects. Biodiesel is a mixture of long chain fatty acid methyl esters derived from fats and oils. This research study presents high quality experimental data for one large fatty acid methyl ester, methyl decanoate, and models its combustion using an improved skeletal mechanism. The results indicate that methyl decanoate is consumed via abstraction of hydrogen atoms to produce fuel radicals, which ultimately lead to the production of alkenes. The ester moiety in methyl decanoate leads to the formation of low molecular

  10. Chemical kinetic modeling of exhaust hydrocarbon oxidation

    SciTech Connect

    Wu, K.C.; Hochgreb, S.; Norris, M.B. . Dept. of Mechanical Engineering)

    1995-01-01

    Numerical simulations of the oxidation of unburned hydrocarbons from spark ignition engines were made based on full-chemistry, zero-dimensional models and compared with experiments for engine-out hydrocarbons and exhaust port oxidation. Simple correlations can be drawn between calculated results for hydrocarbon oxidation half-lives in plug or stirred reactors and measured hydrocarbon emissions. the extent of reaction through the exhaust port was simulated using calculated temperature histories for each burned gas mass element leaving the cylinder, coupled to detailed chemical kinetic rate equations. The results show that, for the fuels considered, the extent of oxidation of the remaining unburned fuel measured through the exhaust can be bracketed by the calculated results for the well-mixed (average) and core (adiabatically expanded) temperatures in the exhaust. Most of the oxidation is shown to occur at the very early exhaust times. For the paraffins considered, comparisons of simulations and experiments suggest that fuel oxidation is partially controlled by the mixing of cold gases at the initial stages of exhaust, where temperatures are high and the cold unburned mixture emerges from the wall layers into the exhaust jet. These conclusions are supported by the relatively small measured dependence on fuel type of the extent of oxidation in the exhaust, and by the resulting ratio of fuel to nonfuel hydrocarbons in the exhaust port exit.

  11. Detailed Chemical Kinetic Modeling of Cyclohexane Oxidation

    SciTech Connect

    Silke, E J; Pitz, W J; Westbrook, C K; Ribaucour, M

    2006-11-10

    A detailed chemical kinetic mechanism has been developed and used to study the oxidation of cyclohexane at both low and high temperatures. Reaction rate constant rules are developed for the low temperature combustion of cyclohexane. These rules can be used for in chemical kinetic mechanisms for other cycloalkanes. Since cyclohexane produces only one type of cyclohexyl radical, much of the low temperature chemistry of cyclohexane is described in terms of one potential energy diagram showing the reaction of cyclohexyl radical + O{sub 2} through five, six and seven membered ring transition states. The direct elimination of cyclohexene and HO{sub 2} from RO{sub 2} is included in the treatment using a modified rate constant of Cavallotti et al. Published and unpublished data from the Lille rapid compression machine, as well as jet-stirred reactor data are used to validate the mechanism. The effect of heat loss is included in the simulations, an improvement on previous studies on cyclohexane. Calculations indicated that the production of 1,2-epoxycyclohexane observed in the experiments can not be simulated based on the current understanding of low temperature chemistry. Possible 'alternative' H-atom isomerizations leading to different products from the parent O{sub 2}QOOH radical were included in the low temperature chemical kinetic mechanism and were found to play a significant role.

  12. Detailed Chemical Kinetic Modeling of Hydrazine Decomposition

    NASA Technical Reports Server (NTRS)

    Meagher, Nancy E.; Bates, Kami R.

    2000-01-01

    The purpose of this research project is to develop and validate a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. Hydrazine is used extensively in aerospace propulsion, and although liquid hydrazine is not considered detonable, many fuel handling systems create multiphase mixtures of fuels and fuel vapors during their operation. Therefore, a thorough knowledge of the decomposition chemistry of hydrazine under a variety of conditions can be of value in assessing potential operational hazards in hydrazine fuel systems. To gain such knowledge, a reasonable starting point is the development and validation of a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. A reasonably complete mechanism was published in 1996, however, many of the elementary steps included had outdated rate expressions and a thorough investigation of the behavior of the mechanism under a variety of conditions was not presented. The current work has included substantial revision of the previously published mechanism, along with a more extensive examination of the decomposition behavior of hydrazine. An attempt to validate the mechanism against the limited experimental data available has been made and was moderately successful. Further computational and experimental research into the chemistry of this fuel needs to be completed.

  13. Elimination kinetic model for organic chemicals in earthworms.

    PubMed

    Dimitrova, N; Dimitrov, S; Georgieva, D; Van Gestel, C A M; Hankard, P; Spurgeon, D; Li, H; Mekenyan, O

    2010-08-15

    Mechanistic understanding of bioaccumulation in different organisms and environments should take into account the influence of organism and chemical depending factors on the uptake and elimination kinetics of chemicals. Lipophilicity, metabolism, sorption (bioavailability) and biodegradation of chemicals are among the important factors that may significantly affect the bioaccumulation process in soil organisms. This study attempts to model elimination kinetics of organic chemicals in earthworms by accounting for the effects of both chemical and biological properties, including metabolism. The modeling approach that has been developed is based on the concept for simulating metabolism used in the BCF base-line model developed for predicting bioaccumulation in fish. Metabolism was explicitly accounted for by making use of the TIMES engine for simulation of metabolism and a set of principal transformations. Kinetic characteristics of transformations were estimated on the basis of observed kinetics data for the elimination of organic chemicals from earthworms. PMID:20185163

  14. A Detailed Chemical Kinetic Model for TNT

    SciTech Connect

    Pitz, W J; Westbrook, C K

    2005-01-13

    A detailed chemical kinetic mechanism for 2,4,6-tri-nitrotoluene (TNT) has been developed to explore problems of explosive performance and soot formation during the destruction of munitions. The TNT mechanism treats only gas-phase reactions. Reactions for the decomposition of TNT and for the consumption of intermediate products formed from TNT are assembled based on information from the literature and on current understanding of aromatic chemistry. Thermodynamic properties of intermediate and radical species are estimated by group additivity. Reaction paths are developed based on similar paths for aromatic hydrocarbons. Reaction-rate constant expressions are estimated from the literature and from analogous reactions where the rate constants are available. The detailed reaction mechanism for TNT is added to existing reaction mechanisms for RDX and for hydrocarbons. Computed results show the effect of oxygen concentration on the amount of soot precursors that are formed in the combustion of RDX and TNT mixtures in N{sub 2}/O{sub 2} mixtures.

  15. Chemical Kinetic Models for HCCI and Diesel Combustion

    SciTech Connect

    Pitz, W J; Westbrook, C K; Mehl, M; Sarathy, S M

    2010-11-15

    Predictive engine simulation models are needed to make rapid progress towards DOE's goals of increasing combustion engine efficiency and reducing pollutant emissions. These engine simulation models require chemical kinetic submodels to allow the prediction of the effect of fuel composition on engine performance and emissions. Chemical kinetic models for conventional and next-generation transportation fuels need to be developed so that engine simulation tools can predict fuel effects. The objectives are to: (1) Develop detailed chemical kinetic models for fuel components used in surrogate fuels for diesel and HCCI engines; (2) Develop surrogate fuel models to represent real fuels and model low temperature combustion strategies in HCCI and diesel engines that lead to low emissions and high efficiency; and (3) Characterize the role of fuel composition on low temperature combustion modes of advanced combustion engines.

  16. Kinetic model of continuous-wave flow chemical lasers

    NASA Astrophysics Data System (ADS)

    Gao, Z.; X., E.

    1982-02-01

    A kinetic approach to modeling the gain in a chemical wave continuous laser when the lasing frequency is coincident with the center of the line shape is presented. Governing equations are defined for the relaxing behavior of an initially nonequilibrium distribution toward the local equilibrium Boltzmann-Maxwellian distribution. A new gain is introduced which is related to the thermal motion of the molecules and cold-reaction and premixed CW models are discussed. Coincidence of the lasing frequency with the line shape is demonstrated to result in a radiative intensity within the homogeneous broadening limit. The rate model predictions are compared with those of the kinetic model. It is found that when the broadening parameter is less than 0.2 the kinetic model more accurately describes the behavior of the CW chemical laser.

  17. An efficient chemical kinetics solver using high dimensional model representation

    SciTech Connect

    Shorter, J.A.; Ip, P.C.; Rabitz, H.A.

    1999-09-09

    A high dimensional model representation (HDMR) technique is introduced to capture the input-output behavior of chemical kinetic models. The HDMR expresses the output chemical species concentrations as a rapidly convergent hierarchical correlated function expansion in the input variables. In this paper, the input variables are taken as the species concentrations at time t{sub i} and the output is the concentrations at time t{sub i} + {delta}, where {delta} can be much larger than conventional integration time steps. A specially designed set of model runs is performed to determine the correlated functions making up the HDMR. The resultant HDMR can be used to (1) identify the key input variables acting independently or cooperatively on the output, and (2) create a high speed fully equivalent operational model (FEOM) serving to replace the original kinetic model and its differential equation solver. A demonstration of the HDMR technique is presented for stratospheric chemical kinetics. The FEOM proved to give accurate and stable chemical concentrations out to long times of many years. In addition, the FEOM was found to be orders of magnitude faster than a conventional stiff equation solver. This computational acceleration should have significance in many chemical kinetic applications.

  18. Model reduction for chemical kinetics: An optimization approach

    SciTech Connect

    Petzold, L.; Zhu, W.

    1999-04-01

    The kinetics of a detailed chemically reacting system can potentially be very complex. Although the chemist may be interested in only a few species, the reaction model almost always involves a much larger number of species. Some of those species are radicals, which are very reactive species and can be important intermediaries in the reaction scheme. A large number of elementary reactions can occur among the species; some of these reactions are fast and some are slow. The aim of simplified kinetics modeling is to derive the simplest reaction system which retains the essential features of the full system. An optimization-based method for reduction of the number of species and reactions in chemical kinetics model is described. Numerical results for several reaction mechanisms illustrate the potential of this approach.

  19. Computer-Aided Construction of Chemical Kinetic Models

    SciTech Connect

    Green, William H.

    2014-12-31

    The combustion chemistry of even simple fuels can be extremely complex, involving hundreds or thousands of kinetically significant species. The most reasonable way to deal with this complexity is to use a computer not only to numerically solve the kinetic model, but also to construct the kinetic model in the first place. Because these large models contain so many numerical parameters (e.g. rate coefficients, thermochemistry) one never has sufficient data to uniquely determine them all experimentally. Instead one must work in “predictive” mode, using theoretical rather than experimental values for many of the numbers in the model, and as appropriate refining the most sensitive numbers through experiments. Predictive chemical kinetics is exactly what is needed for computer-aided design of combustion systems based on proposed alternative fuels, particularly for early assessment of the value and viability of proposed new fuels before those fuels are commercially available. This project was aimed at making accurate predictive chemical kinetics practical; this is a challenging goal which requires a range of science advances. The project spanned a wide range from quantum chemical calculations on individual molecules and elementary-step reactions, through the development of improved rate/thermo calculation procedures, the creation of algorithms and software for constructing and solving kinetic simulations, the invention of methods for model-reduction while maintaining error control, and finally comparisons with experiment. Many of the parameters in the models were derived from quantum chemistry calculations, and the models were compared with experimental data measured in our lab or in collaboration with others.

  20. Chemical kinetic modeling of H{sub 2} applications

    SciTech Connect

    Marinov, N.M.; Westbrook, C.K.; Cloutman, L.D.

    1995-09-01

    Work being carried out at LLNL has concentrated on studies of the role of chemical kinetics in a variety of problems related to hydrogen combustion in practical combustion systems, with an emphasis on vehicle propulsion. Use of hydrogen offers significant advantages over fossil fuels, and computer modeling provides advantages when used in concert with experimental studies. Many numerical {open_quotes}experiments{close_quotes} can be carried out quickly and efficiently, reducing the cost and time of system development, and many new and speculative concepts can be screened to identify those with sufficient promise to pursue experimentally. This project uses chemical kinetic and fluid dynamic computational modeling to examine the combustion characteristics of systems burning hydrogen, either as the only fuel or mixed with natural gas. Oxidation kinetics are combined with pollutant formation kinetics, including formation of oxides of nitrogen but also including air toxics in natural gas combustion. We have refined many of the elementary kinetic reaction steps in the detailed reaction mechanism for hydrogen oxidation. To extend the model to pressures characteristic of internal combustion engines, it was necessary to apply theoretical pressure falloff formalisms for several key steps in the reaction mechanism. We have continued development of simplified reaction mechanisms for hydrogen oxidation, we have implemented those mechanisms into multidimensional computational fluid dynamics models, and we have used models of chemistry and fluid dynamics to address selected application problems. At the present time, we are using computed high pressure flame, and auto-ignition data to further refine the simplified kinetics models that are then to be used in multidimensional fluid mechanics models. Detailed kinetics studies have investigated hydrogen flames and ignition of hydrogen behind shock waves, intended to refine the detailed reactions mechanisms.

  1. Chemical Kinetic Modeling of Combustion of Automotive Fuels

    SciTech Connect

    Pitz, W J; Westbrook, C K; Silke, E J

    2006-11-10

    The objectives of this report are to: (1) Develop detailed chemical kinetic reaction models for components of fuels, including olefins and cycloalkanes used in diesel, spark-ignition and HCCI engines; (2) Develop surrogate mixtures of hydrocarbon components to represent real fuels and lead to efficient reduced combustion models; and (3) Characterize the role of fuel composition on production of emissions from practical automotive engines.

  2. A chemical kinetic modeling study of chlorinated hydrocarbon combustion

    SciTech Connect

    Pitz, W.J.; Westbrook, C.K.

    1990-09-05

    The combustion of chloroethane is modeled as a stirred reactor so that we can study critical emission characteristics of the reactor as a function of residence time. We examine important operating conditions such as pressure, temperature, and equivalence ratio and their influence on destructive efficiency of chloroethane. The model uses a detailed chemical kinetic mechanism that we have developed previously for C{sub 3} hydrocarbons. We have added to this mechanism the chemical kinetic mechanism for C{sub 2} chlorinated hydrocarbons developed by Senkan and coworkers. In the modeling calculations, sensitivity coefficients are determined to find which reaction-rate constants have the largest effect on destructive efficiency. 24 refs., 6 figs., 1 tab.

  3. Chemical Kinetic Models for HCCI and Diesel Combustion

    SciTech Connect

    Pitz, W J; Westbook, C K; Mehl, M

    2008-10-30

    Hydrocarbon fuels for advanced combustion engines consist of complex mixtures of hundreds or even thousands of different components. These components can be grouped into a number of chemically distinct classes, consisting of n-paraffins, branched paraffins, cyclic paraffins, olefins, oxygenates, and aromatics. Biodiesel contains its own unique chemical class called methyl esters. The fractional amounts of these chemical classes are quite different in gasoline, diesel fuel, oil-sand derived fuels and bio-derived fuels, which contributes to the very different combustion characteristics of each of these types of combustion systems. The objectives of this project are: (1) Develop detailed chemical kinetic models for fuel components used in surrogate fuels for diesel and HCCI engines; (2) Develop surrogate fuel models to represent real fuels and model low temperature combustion strategies in HCCI and diesel engines that lead to low emissions and high efficiency; and (3) Characterize the role of fuel composition on low temperature combustion modes of advanced combustion engines.

  4. Progress in Chemical Kinetic Modeling for Surrogate Fuels

    SciTech Connect

    Pitz, W J; Westbrook, C K; Herbinet, O; Silke, E J

    2008-06-06

    Gasoline, diesel, and other alternative transportation fuels contain hundreds to thousands of compounds. It is currently not possible to represent all these compounds in detailed chemical kinetic models. Instead, these fuels are represented by surrogate fuel models which contain a limited number of representative compounds. We have been extending the list of compounds for detailed chemical models that are available for use in fuel surrogate models. Detailed models for components with larger and more complicated fuel molecular structures are now available. These advancements are allowing a more accurate representation of practical and alternative fuels. We have developed detailed chemical kinetic models for fuels with higher molecular weight fuel molecules such as n-hexadecane (C16). Also, we can consider more complicated fuel molecular structures like cyclic alkanes and aromatics that are found in practical fuels. For alternative fuels, the capability to model large biodiesel fuels that have ester structures is becoming available. These newly addressed cyclic and ester structures in fuels profoundly affect the reaction rate of the fuel predicted by the model. Finally, these surrogate fuel models contain large numbers of species and reactions and must be reduced for use in multi-dimensional models for spark-ignition, HCCI and diesel engines.

  5. Towards cleaner combustion engines through groundbreaking detailed chemical kinetic models

    PubMed Central

    Battin-Leclerc, Frédérique; Blurock, Edward; Bounaceur, Roda; Fournet, René; Glaude, Pierre-Alexandre; Herbinet, Olivier; Sirjean, Baptiste; Warth, V.

    2013-01-01

    In the context of limiting the environmental impact of transportation, this paper reviews new directions which are being followed in the development of more predictive and more accurate detailed chemical kinetic models for the combustion of fuels. In the first part, the performance of current models, especially in terms of the prediction of pollutant formation, is evaluated. In the next parts, recent methods and ways to improve these models are described. An emphasis is given on the development of detailed models based on elementary reactions, on the production of the related thermochemical and kinetic parameters, and on the experimental techniques available to produce the data necessary to evaluate model predictions under well defined conditions. PMID:21597604

  6. Bio-butanol: Combustion properties and detailed chemical kinetic model

    SciTech Connect

    Black, G.; Curran, H.J.; Pichon, S.; Simmie, J.M.; Zhukov, V.

    2010-02-15

    Autoignition delay time measurements were performed at equivalence ratios of 0.5, 1 and 2 for butan-1-ol at reflected shock pressures of 1, 2.6 and 8 atm at temperatures from 1100 to 1800 K. High-level ab initio calculations were used to determine enthalpies of formation and consequently bond dissociation energies for each bond in the alcohol. A detailed chemical kinetic model consisting of 1399 reactions involving 234 species was constructed and tested against the delay times and also against recent jet-stirred reactor speciation data with encouraging results. The importance of enol chemistry is highlighted. (author)

  7. Incorporation of chemical kinetic models into process control

    SciTech Connect

    Herget, C.J.; Frazer, J.W.

    1981-07-08

    An important consideration in chemical process control is to determine the precise rationing of reactant streams, particularly when a large time delay exists between the mixing of the reactants and the measurement of the product. In this paper, a method is described for incorporating chemical kinetic models into the control strategy in order to achieve optimum operating conditions. The system is first characterized by determining a reaction rate surface as a function of all input reactant concentrations over a feasible range. A nonlinear constrained optimization program is then used to determine the combination of reactants which produces the specified yield at minimum cost. This operating condition is then used to establish the nominal concentrations of the reactants. The actual operation is determined through a feedback control system employing a Smith predictor. The method is demonstrated on a laboratory bench scale enzyme reactor.

  8. High Temperature Chemical Kinetic Combustion Modeling of Lightly Methylated Alkanes

    SciTech Connect

    Sarathy, S M; Westbrook, C K; Pitz, W J; Mehl, M

    2011-03-01

    Conventional petroleum jet and diesel fuels, as well as alternative Fischer-Tropsch (FT) fuels and hydrotreated renewable jet (HRJ) fuels, contain high molecular weight lightly branched alkanes (i.e., methylalkanes) and straight chain alkanes (n-alkanes). Improving the combustion of these fuels in practical applications requires a fundamental understanding of large hydrocarbon combustion chemistry. This research project presents a detailed high temperature chemical kinetic mechanism for n-octane and three lightly branched isomers octane (i.e., 2-methylheptane, 3-methylheptane, and 2,5-dimethylhexane). The model is validated against experimental data from a variety of fundamental combustion devices. This new model is used to show how the location and number of methyl branches affects fuel reactivity including laminar flame speed and species formation.

  9. Upper D region chemical kinetic modeling of LORE relaxation times

    NASA Astrophysics Data System (ADS)

    Gordillo-Vázquez, F. J.; Luque, A.; Haldoupis, C.

    2016-04-01

    The recovery times of upper D region electron density elevations, caused by lightning-induced electromagnetic pulses (EMP), are modeled. The work was motivated from the need to understand a recently identified narrowband VLF perturbation named LOREs, an acronym for LOng Recovery Early VLF events. LOREs associate with long-living electron density perturbations in the upper D region ionosphere; they are generated by strong EMP radiated from large peak current intensities of ±CG (cloud to ground) lightning discharges, known also to be capable of producing elves. Relaxation model scenarios are considered first for a weak enhancement in electron density and then for a much stronger one caused by an intense lightning EMP acting as an impulsive ionization source. The full nonequilibrium kinetic modeling of the perturbed mesosphere in the 76 to 92 km range during LORE-occurring conditions predicts that the electron density relaxation time is controlled by electron attachment at lower altitudes, whereas above 79 km attachment is balanced totally by associative electron detachment so that electron loss at these higher altitudes is controlled mainly by electron recombination with hydrated positive clusters H+(H2O)n and secondarily by dissociative recombination with NO+ ions, a process which gradually dominates at altitudes >88 km. The calculated recovery times agree fairly well with LORE observations. In addition, a simplified (quasi-analytic) model build for the key charged species and chemical reactions is applied, which arrives at similar results with those of the full kinetic model. Finally, the modeled recovery estimates for lower altitudes, that is <79 km, are in good agreement with the observed short recovery times of typical early VLF events, which are known to be associated with sprites.

  10. Integration Strategies for Efficient Multizone Chemical Kinetics Models

    SciTech Connect

    McNenly, M J; Havstad, M A; Aceves, S M; Pitz, W J

    2009-10-15

    Three integration strategies are developed and tested for the stiff, ordinary differential equation (ODE) integrators used to solve the fully coupled multizone chemical kinetics model. Two of the strategies tested are found to provide more than an order of magnitude of improvement over the original, basic level of usage for the stiff ODE solver. One of the faster strategies uses a decoupled, or segregated, multizone model to generate an approximate Jacobian. This approach yields a 35-fold reduction in the computational cost for a 20 zone model. Using the same approximate Jacobian as a preconditioner for an iterative Krylov-type linear system solver, the second improved strategy achieves a 75-fold reduction in the computational cost for a 20 zone model. The faster strategies achieve their cost savings with no significant loss of accuracy. The pressure, temperature and major species mass fractions agree with the solution from the original integration approach to within six significant digits; and the radical mass fractions agree with the original solution to within four significant digits. The faster strategies effectively change the cost scaling of the multizone model from cubic to quadratic, with respect to the number of zones. As a consequence of the improved scaling, the 40 zone model offers more than a 250-fold cost savings over the basic calculation.

  11. Chemical Kinetic Modeling of HMX and TATB Laser Ignition Tests

    SciTech Connect

    Tarver, C M

    2004-03-02

    Recent laser ignition experiments on octahydro-1,3,5,7-tetranitro-1,3,5,7-terrazocine (HMX) and 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) subjected to laser fluxes ranging from 10 to 800 W/cm{sup 2} produced ignition times from seconds to milliseconds. Global chemical kinetic thermal decomposition models for HMX and TATB have been developed to calculate times to thermal explosion for experiments in the seconds to days time frame. These models are applied to the laser ignition experimental data in this paper. Excellent agreement was obtained for TATB, while the calculated ignition times were longer than experiment for HMX at lower laser fluxes. At the temperatures produced in the laser experiments, HMX melts. Melting generally increases condensed phase reaction rates so faster rates were used for three of the HMX reaction rates. This improved agreement with experiments at the lower laser fluxes but yielded very fast ignition at high fluxes. The calculated times to ignition are in reasonable agreement with the laser ignition experiments, and this justifies the use of these models for estimating reaction times at impact and shock ''hot spot'' temperatures.

  12. Chemical kinetic modeling of H{sub 2} applications

    SciTech Connect

    Westbrook, C.K.; Marinov, N.; Pitz, W.J.; Curran, H.

    1996-10-01

    This project is intended to develop detailed and simplified kinetic reaction mechanisms for the combustion of practical systems fueled by hydrogen, and then to use those mechanisms to examine the performance, efficiency, pollutant emissions, and other characteristics of those systems. During the last year, a H2/NOx mechanism has been developed that gives much improved predictions of NOx emissions compared to experimental data. Preliminary chemical kinetic and equilibrium calculations have been performed in support of Br2-H2O experiments to be conducted at NREL. Hydrogen, hydrogen/methane and hydrogen/natural gas mixtures have been investigated in a knock-rating engine to assess their automotive knock characteristics. The authors are currently developing the simplified analog reaction mechanisms that are computationally simple, yet still reproduce many of the macroscopic features of flame propagation.

  13. Parameter Estimates in Differential Equation Models for Chemical Kinetics

    ERIC Educational Resources Information Center

    Winkel, Brian

    2011-01-01

    We discuss the need for devoting time in differential equations courses to modelling and the completion of the modelling process with efforts to estimate the parameters in the models using data. We estimate the parameters present in several differential equation models of chemical reactions of order n, where n = 0, 1, 2, and apply more general…

  14. Detailed Chemical Kinetic Modeling of Diesel Combustion with Oxygenated Fuels

    SciTech Connect

    Curran, H J; Fisher, E M; Glaude, P-A; Marinov, N M; Pitz, W J; Westbrook, C K; Flynn, P F; Durrett, R P; zur Loye, A O; Akinyemi, O C; Dryer, F L

    2000-01-11

    Emission standards for diesel engines in vehicles have been steadily reduced in recent years, and a great deal of research and development effort has been focused on reducing particulate and nitrogen oxide emissions. One promising approach to reducing emissions involves the addition of oxygen to the fuel, generally by adding an oxygenated compound to the normal diesel fuel. Miyamoto et al. [1] showed experimentally that particulate levels can be significantly reduced by adding oxygenated species to the fuel. They found the Bosch smoke number (a measure of the particulate or soot levels in diesel exhaust) falls from about 55% for conventional diesel fuel to less than 1% when the oxygen content of the fuel is above about 25% by mass, as shown in Figure 1. It has been well established that addition of oxygenates to automotive fuel, including both diesel fuel as well as gasoline, reduces NOx and CO emissions by reducing flame temperatures. This is the basis for addition of oxygenates to produce reformulated gasoline in selected portions of the country. Of course, this is also accompanied by a slight reduction in fuel economy. A new overall picture of diesel combustion has been developed by Dec [2], in which laser diagnostic studies identified stages in diesel combustion that had not previously been recognized. These stages are summarized in Figure 2. The evolution of the diesel spray is shown, starting as a liquid jet that vaporizes and entrains hot air from the combustion chamber. This relatively steady process continues as long as fuel is being injected. In particular, Dec showed that the fuel spray vaporizes and mixes with air and products of earlier combustion to provide a region in which a gas phase, premixed fuel-rich ignition and burn occurs. The products of this ignition are then observed experimentally to lead rapidly to formation of soot particles, which subsequently are consumed in a diffusion flame. Recently, Flynn et al. [3] used a chemical kinetic and

  15. Recent Advances in Detailed Chemical Kinetic Models for Large Hydrocarbon and Biodiesel Transportation Fuels

    SciTech Connect

    Westbrook, C K; Pitz, W J; Curran, H J; Herbinet, O; Mehl, M

    2009-03-30

    n-Hexadecane and 2,2,4,4,6,8,8-heptamethylnonane represent the primary reference fuels for diesel that are used to determine cetane number, a measure of the ignition property of diesel fuel. With the development of chemical kinetics models for these two primary reference fuels for diesel, a new capability is now available to model diesel fuel ignition. Also, we have developed chemical kinetic models for a whole series of large n-alkanes and a large iso-alkane to represent these chemical classes in fuel surrogates for conventional and future fuels. Methyl decanoate and methyl stearate are large methyl esters that are closely related to biodiesel fuels, and kinetic models for these molecules have also been developed. These chemical kinetic models are used to predict the effect of the fuel molecule size and structure on ignition characteristics under conditions found in internal combustion engines.

  16. Chemical kinetic modeling of component mixtures relevant to gasoline

    SciTech Connect

    Mehl, M; Curran, H J; Pitz, W J; Westbrook, C K

    2009-02-13

    Real fuels are complex mixtures of thousands of hydrocarbon compounds including linear and branched paraffins, naphthenes, olefins and aromatics. It is generally agreed that their behavior can be effectively reproduced by simpler fuel surrogates containing a limited number of components. In this work, a recently revised version of the kinetic model by the authors is used to analyze the combustion behavior of several components relevant to gasoline surrogate formulation. Particular attention is devoted to linear and branched saturated hydrocarbons (PRF mixtures), olefins (1-hexene) and aromatics (toluene). Model predictions for pure components, binary mixtures and multi-component gasoline surrogates are compared with recent experimental information collected in rapid compression machine, shock tube and jet stirred reactors covering a wide range of conditions pertinent to internal combustion engines. Simulation results are discussed focusing attention on the mixing effects of the fuel components.

  17. Detailed and global chemical kinetics model for hydrogen

    SciTech Connect

    Marinov, N.M.; Westbrook, C.K.; Pitz, W.J.

    1995-03-01

    Detailed and global chemical kinetic computations for hydrogen-air mixtures have been performed to describe flame propagation, flame structure and ignition phenomena. Simulations of laminar flame speeds, flame compositions and shock tube ignition delay times have been successfully performed. Sensitivity analysis was applied to determine the governing rate-controlling reactions for the experimental data sets examined. In the flame propagation and structure studies, the reactions, OH + H{sub 2} = H{sub 2}0 + H, 0 + H{sub 2} = OH + H and 0 + OH = 0{sub 2} + H were the most important in flames. The shock tube ignition delay time study indicated the H + 0{sub 2} + M = H0{sub 2} + M (M = N{sub 2}, H{sub 2}) and 0 + OH = 0{sub 2} + H reactions controlled ignition. A global rate expression for a one-step overall reaction was developed and validated against experimental hydrogen-air laminar flame speed data. The global reaction expression was determined to be 1.8 {times} 10{sup 13} exp({minus}17614K/T)[H{sub 2}]{sup 1.0}[O{sub 2}]{sup 0.5} for the single step reaction H{sub 2} + 1/2O{sub 2} = H{sub 2}O.

  18. Chemical kinetic modeling of component mixtures relevant to gasoline

    SciTech Connect

    Mehl, M; Curran, H J; Pitz, W J; Dooley, S; Westbrook, C K

    2008-05-29

    Detailed kinetic models of pyrolysis and combustion of hydrocarbon fuels are nowadays widely used in the design of internal combustion engines and these models are effectively applied to help meet the increasingly stringent environmental and energetic standards. In previous studies by the combustion community, such models not only contributed to the understanding of pure component combustion, but also provided a deeper insight into the combustion behavior of complex mixtures. One of the major challenges in this field is now the definition and the development of appropriate surrogate models able to mimic the actual features of real fuels. Real fuels are complex mixtures of thousands of hydrocarbon compounds including linear and branched paraffins, naphthenes, olefins and aromatics. Their behavior can be effectively reproduced by simpler fuel surrogates containing a limited number of components. Aside the most commonly used surrogates containing iso-octane and n-heptane only, the so called Primary Reference Fuels (PRF), new mixtures have recently been suggested to extend the reference components in surrogate mixtures to also include alkenes and aromatics. It is generally agreed that, including representative species for all the main classes of hydrocarbons which can be found in real fuels, it is possible to reproduce very effectively in a wide range of operating conditions not just the auto-ignition propensity of gasoline or Diesel fuels, but also their physical properties and their combustion residuals [1]. In this work, the combustion behavior of several components relevant to gasoline surrogate formulation is computationally examined. The attention is focused on the autoignition of iso-octane, hexene and their mixtures. Some important issues relevant to the experimental and modeling investigation of such fuels are discussed with the help of rapid compression machine data and calculations. Following the model validation, the behavior of mixtures is discussed on the

  19. Chemical kinetic modeling of propane oxidation behind shock waves

    NASA Technical Reports Server (NTRS)

    Mclain, A. G.; Jachimowski, C. J.

    1977-01-01

    The stoichiometric combustion of propane behind incident shock waves was studied experimentally and analytically over a temperature range from 1700 K to 2600 K and a pressure range from 1.2 to 1.9 atm. Measurements of the concentrations of carbon monoxide (CO) and carbon dioxide (CO2) and the product of the oxygen atom and carbon dioxide concentrations (O)(CO) were made after passage of the incident shock wave. A kinetic mechanism was developed which, when used in a computer program for a flowing, reacting gas behind an incident shock wave predicted experimentally measured results quite well. Ignition delay times from the literature were also predicted quite well. The kinetic mechanism consisted of 59 individual kinetic steps.

  20. Computer Simulation in Chemical Kinetics

    ERIC Educational Resources Information Center

    Anderson, Jay Martin

    1976-01-01

    Discusses the use of the System Dynamics technique in simulating a chemical reaction for kinetic analysis. Also discusses the use of simulation modelling in biology, ecology, and the social sciences, where experimentation may be impractical or impossible. (MLH)

  1. Chemical Kinetics for Modeling Silicon Epitaxy from Chlorosilanes

    SciTech Connect

    Balakrishna, A.; Chacin, J.M.; Comita, P.B.; Haas, B.; Ho, P.; Thilderkvist, A.

    1998-11-24

    A reaction mechanism has been developed that describes the gas-phas 0971 and surface reactions involved in the chemical vapor deposition of Si from chlorosilanes. Good agreement with deposition rate data from a single wafer reactor with no wafer rotation has been attained over a range of gas mixtures, total flow rates, and reactor temperatures.

  2. KEMOD: A mixed chemical kinetic and equilibrium model of aqueous and solid phase geochemical reactions

    SciTech Connect

    Yeh, G.T.; Iskra, G.A.; Szecsody, J.E.; Zachara, J.M.; Streile, G.P.

    1995-01-01

    This report presents the development of a mixed chemical Kinetic and Equilibrium MODel in which every chemical species can be treated either as a equilibrium-controlled or as a kinetically controlled reaction. The reaction processes include aqueous complexation, adsorption/desorption, ion exchange, precipitation/dissolution, oxidation/reduction, and acid/base reactions. Further development and modification of KEMOD can be made in: (1) inclusion of species switching solution algorithms, (2) incorporation of the effect of temperature and pressure on equilibrium and rate constants, and (3) extension to high ionic strength.

  3. Chemical Kinetics Database

    National Institute of Standards and Technology Data Gateway

    SRD 17 NIST Chemical Kinetics Database (Web, free access)   The NIST Chemical Kinetics Database includes essentially all reported kinetics results for thermal gas-phase chemical reactions. The database is designed to be searched for kinetics data based on the specific reactants involved, for reactions resulting in specified products, for all the reactions of a particular species, or for various combinations of these. In addition, the bibliography can be searched by author name or combination of names. The database contains in excess of 38,000 separate reaction records for over 11,700 distinct reactant pairs. These data have been abstracted from over 12,000 papers with literature coverage through early 2000.

  4. Kinetic Monte Carlo models for the study of chemical reactions in the Earth's upper atmosphere

    NASA Astrophysics Data System (ADS)

    Turchak, L. I.; Shematovich, V. I.

    2016-06-01

    A stochastic approach to study the non-equilibrium chemistry in the Earth's upper atmosphere is presented, which has been developed over a number of years. Kinetic Monte Carlo models based on this approach are an effective tool for investigating the role of suprathermal particles both in local variations of the atmospheric chemical composition and in the formation of the hot planetary corona.

  5. Chemical kinetics and photochemical data for use in stratospheric modeling

    NASA Technical Reports Server (NTRS)

    Demore, W. B.; Margitan, J. J.; Molina, M. J.; Watson, R. T.; Golden, D. M.; Hampson, R. F.; Kurylo, M. J.; Howard, C. J.; Ravishankara, A. R.

    1985-01-01

    Rate constants and photochemical cross sections are presented. The primary application of the data is for modeling of the stratospheric processes, with particular emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena.

  6. Chemical kinetic and photochemical data for use in stratospheric modelling

    NASA Technical Reports Server (NTRS)

    Demore, W. B.; Stief, L. J.; Kaufman, F.; Golden, D. M.; Hampton, R. F.; Kurylo, M. J.; Margitan, J. J.; Molina, M. J.; Watson, R. T.

    1979-01-01

    An evaluated set of rate constants and photochemical cross sections were compiled for use in modelling stratospheric processes. The data are primarily relevant to the ozone layer, and its possible perturbation by anthropogenic activities. The evaluation is current to, approximately, January, 1979.

  7. Chemical kinetics and photochemical data for use in stratospheric modeling

    NASA Technical Reports Server (NTRS)

    Demore, W. B.; Sander, S. P.; Golden, D. M.; Hampson, R. F.; Kurylo, M. J.; Howard, C. J.; Ravishankara, A. R.; Kolb, C. E.; Molina, M. J.

    1992-01-01

    As part of a series of evaluated sets, rate constants and photochemical cross sections compiled by the NASA Panel for Data Evaluation are provided. The primary application of the data is in the modeling of stratospheric processes, with particular emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena. Copies of this evaluation are available from the Jet Propulsion Laboratory.

  8. Transport Properties of a Kinetic Model for Chemical Reactions without Barriers

    SciTech Connect

    Alves, Giselle M.; Kremer, Gilberto M.; Soares, Ana Jacinta

    2011-05-20

    A kinetic model of the Boltzmann equation for chemical reactions without energy barrier is considered here with the aim of evaluating the reaction rate and characterizing the transport coefficient of shear viscosity for the reactive system. The Chapman-Enskog solution of the Boltzmann equation is used to compute the chemical reaction effects, in a flow regime for which the reaction process is close to the final equilibrium state. Some numerical results are provided illustrating that the considered chemical reaction without energy barrier can induce an appreciable influence on the reaction rate and on the transport coefficient of shear viscosity.

  9. Chemical kinetic modeling of chlorinated hydrocarbons under stirred-reactor conditions

    SciTech Connect

    Pitz, W.J.; Westbrook, C.K.

    1990-10-04

    The combustin of chloroethane is modeled as a stirred reactor so that we can study critical emission characteristics of the reactor as a function of residence time. We examine important operating conditions such as pressure, temperature, and equivalence ratio and their influence on destructive efficiency of chloroethane and production of other chlorinated products. The model uses a detailed chemical kinetic mechanism that we have developed previously for C{sub 3} hydrocarbons. We have added to this mechanism the chemical kinetic mechanism for C{sub 2} chlorinated hydrocarbons developed by Senkan and coworkers. Some reactions have been added to Senkan's mechanism and some of the reaction-rate expressions have been updated to reflect recent developments in the literature. In the modeling calculations, sensitivity coefficients are determined to find which reaction-rate constants have the largest effect on destructive efficiency. 25 refs., 6 figs., 1 tab.

  10. History and Philosophy of Science through Models: The Case of Chemical Kinetics

    NASA Astrophysics Data System (ADS)

    Justi, Rosária; Gilbert, John K.

    The case for a greater role for the history and philosophy of science in science education is reviewed. It is argued that such a role can only be realised if it is based on both a credible analytical approach to the history and philosophy of science and if the evolution of a sufficient number of major themes in science is known in suitable detail. Adopting Lakatos' Theory of Scientific Research Programmes as the analytical approach, it is proposed that the development, use, and replacement, of specific models forms the core of such programmes.Chemical kinetics was selected as an exemplar major topic in chemistry. Eight models which have played a central role in the evolution of the study of chemical kinetics were identified by an analysis of the literature. The implications that these models have for the teaching and learning of chemistry today are discussed.

  11. Chemical Looping Combustion Kinetics

    SciTech Connect

    Edward Eyring; Gabor Konya

    2009-03-31

    One of the most promising methods of capturing CO{sub 2} emitted by coal-fired power plants for subsequent sequestration is chemical looping combustion (CLC). A powdered metal oxide such as NiO transfers oxygen directly to a fuel in a fuel reactor at high temperatures with no air present. Heat, water, and CO{sub 2} are released, and after H{sub 2}O condensation the CO{sub 2} (undiluted by N{sub 2}) is ready for sequestration, whereas the nickel metal is ready for reoxidation in the air reactor. In principle, these processes can be repeated endlessly with the original nickel metal/nickel oxide participating in a loop that admits fuel and rejects ash, heat, and water. Our project accumulated kinetic rate data at high temperatures and elevated pressures for the metal oxide reduction step and for the metal reoxidation step. These data will be used in computational modeling of CLC on the laboratory scale and presumably later on the plant scale. The oxygen carrier on which the research at Utah is focused is CuO/Cu{sub 2}O rather than nickel oxide because the copper system lends itself to use with solid fuels in an alternative to CLC called 'chemical looping with oxygen uncoupling' (CLOU).

  12. Chemical Kinetics Laboratory Discussion Worksheet

    PubMed Central

    Demoin, Dustin Wayne; Jurisson, Silvia S.

    2013-01-01

    A laboratory discussion worksheet and its answer key provide instructors and students a discussion model to further the students’ understanding of chemical kinetics. This discussion worksheet includes a section for students to augment their previous knowledge about chemical kinetics measurements, an initial check on students’ understanding of basic concepts, a group participation model where students work on solving complex-conceptual problems, and a conclusion to help students connect this discussion to their laboratory or lecture class. Additionally, the worksheet has a detailed solution to a more advanced problem to help students understand how the concepts they have put together relate to problems they will encounter during later formal assessments. PMID:24092948

  13. Calibration of Chemical Kinetic Models Using Simulations of Small-Scale Cookoff Experiments

    SciTech Connect

    Wemhoff, A P; Becker, R C; Burnham, A K

    2008-02-26

    Establishing safe handling limits for explosives in elevated temperature environments is a difficult problem that often requires extensive simulation. The largest influence on predicting thermal cookoff safety lies in the chemical kinetic model used in these simulations, and these kinetic model reaction sequences often contain multiple steps. Several small-scale cookoff experiments, notably Differential Scanning Calorimetry (DSC), Thermogravimetric Analysis (TGA), One-Dimensional Time-to-Explosion (ODTX), and the Scaled Thermal Explosion (STEX) have been performed on various explosives to aid in cookoff behavior determination. Past work has used a single test from this group to create a cookoff model, which does not guarantee agreement with the other experiments. In this study, we update the kinetic parameters of an existing model for the common explosive 2,4,6-Trinitrotoluene (TNT) using DSC and ODTX experimental data at the same time by minimizing a global Figure of Merit based on hydrodynamic simulated data. We then show that the new kinetic model maintains STEX agreement, reduces DSC agreement, and improves ODTX and TGA agreement when compared to the original model. In addition, we describe a means to use implicit hydrodynamic simulations of DSC experiments to develop a reaction model for TNT melting.

  14. Measurement and chemical kinetic model predictions of detonation cell size in methanol-oxygen mixtures

    NASA Astrophysics Data System (ADS)

    Eaton, R.; Zhang, B.; Bergthorson, J. M.; Ng, H. D.

    2012-03-01

    In this study, detonation cell sizes of methanol-oxygen mixtures are experimentally measured at different initial pressures and compositions. Good agreement is found between the experiment data and predictions based on the chemical length scales obtained from a detailed chemical kinetic model. To assess the detonation sensitivity in methanol-oxygen mixtures, the results are compared with those of hydrogen-oxygen and methane-oxygen mixtures. Based on the cell size comparison, it is shown that methanol-oxygen is more detonation sensitive than methane-oxygen but less sensitive than hydrogen-oxygen.

  15. Development and validation of a model for the chemical kinetics of graphite oxidation

    NASA Astrophysics Data System (ADS)

    El-Genk, Mohamed S.; Tournier, Jean-Michel P.

    2011-04-01

    A model of the chemical kinetics and primary reactions of graphite oxidation is developed and successfully validated for pyrolytic carbon thin films. The model uses Gaussian distributions of the activation energies for adsorption and desorption and the measured active surface area (ASA) as a function of burn-off. The activation energies distributions and the pre-exponential rate coefficients for the four elementary oxidation kinetics reactions in the model are obtained from the reported measurements of the gases yields and adsorbed oxygen using a multi-parameter optimization algorithm. The model calculates the production rates of CO and CO 2 and the gasification rate as functions of temperature and oxygen partial pressure, and its predictions are in excellent agreement with reported experimental measurements. Results for pyrolytic carbon thin films show that when the oxygen pressure is kept constant, the gasification rate depends on both temperature and ASA until a full burn-off is reached. By contrast, in a depleting oxygen environment, only partial burn-off is possible; gasification ceases following the consumption of the free oxygen in the enclosure. This model represents the first phase in an ongoing effort to develop a model for predicting the oxidation kinetics of nuclear graphite following a massive air ingress in high temperature reactors.

  16. Modeling of the HiPco process for carbon nanotube production. I. Chemical kinetics.

    PubMed

    Dateo, Christopher E; Gökçen, Tahir; Meyyappan, M

    2002-10-01

    A chemical kinetic model is developed to help understand and optimize the production of single-walled carbon nanotubes via the high-pressure carbon monoxide (HiPco) process, which employs iron pentacarbonyl as the catalyst precursor and carbon monoxide as the carbon feedstock. The model separates the HiPco process into three steps, precursor decomposition, catalyst growth and evaporation, and carbon nanotube production resulting from the catalyst-enhanced disproportionation of carbon monoxide, known as the Boudouard reaction: 2 CO(g)-->C(s) + CO2(g). The resulting detailed model contains 971 species and 1948 chemical reactions. A second model with a reduced reaction set containing 14 species and 22 chemical reactions is developed on the basis of the detailed model and reproduces the chemistry of the major species. Results showing the parametric dependence of temperature, total pressure, and initial precursor partial pressures are presented, with comparison between the two models. The reduced model is more amenable to coupled reacting flow-field simulations, presented in the following article. PMID:12908291

  17. Modeling of the HiPco process for carbon nanotube production. I. Chemical kinetics

    NASA Technical Reports Server (NTRS)

    Dateo, Christopher E.; Gokcen, Tahir; Meyyappan, M.

    2002-01-01

    A chemical kinetic model is developed to help understand and optimize the production of single-walled carbon nanotubes via the high-pressure carbon monoxide (HiPco) process, which employs iron pentacarbonyl as the catalyst precursor and carbon monoxide as the carbon feedstock. The model separates the HiPco process into three steps, precursor decomposition, catalyst growth and evaporation, and carbon nanotube production resulting from the catalyst-enhanced disproportionation of carbon monoxide, known as the Boudouard reaction: 2 CO(g)-->C(s) + CO2(g). The resulting detailed model contains 971 species and 1948 chemical reactions. A second model with a reduced reaction set containing 14 species and 22 chemical reactions is developed on the basis of the detailed model and reproduces the chemistry of the major species. Results showing the parametric dependence of temperature, total pressure, and initial precursor partial pressures are presented, with comparison between the two models. The reduced model is more amenable to coupled reacting flow-field simulations, presented in the following article.

  18. A Shock Tube and Chemical Kinetic Modeling Study of the Oxidation of 2,5-Dimethylfuran

    PubMed Central

    Sirjean, Baptiste; Fournet, René; Glaude, Pierre-Alexandre; Battin-Leclerc, Frédérique; Wang, Weijing; Oehlschlaeger, Matthew A.

    2013-01-01

    A detailed kinetic model describing the oxidation of 2,5-dimethylfuran (DMF), a potential second-generation biofuel, is proposed. The kinetic model is based upon quantum chemical calculations for the initial DMF consumption reactions and important reactions of intermediates. The model is validated by comparison to new DMF shock tube ignition delay time measurements (over the temperature range 1300 – 1831 K and at nominal pressures of 1 and 4 bar) and the DMF pyrolysis speciation measurements of Lifshitz et al. [J. Phys. Chem. A 102 (52) (1998) 10655-10670] Globally, modeling predictions are in good agreement with the considered experimental targets. In particular, ignition delay times are predicted well by the new model, with model experiment deviations of at most a factor of two, and DMF pyrolysis conversion is predicted well, to within experimental scatter of the Lifshitz et al. data. Additionally, comparisons of measured and model predicted pyrolysis speciation provides validation of theoretically calculated channels for the oxidation of DMF. Sensitivity and reaction flux analyses highlight important reactions as well as the primary reaction pathways responsible for the decomposition of DMF and formation and destruction of key intermediate and product species. PMID:23327724

  19. Learning Chemical Kinetics with Spreadsheets.

    ERIC Educational Resources Information Center

    Blickensderfer, Roger

    1990-01-01

    Presented are several simple kinetic systems together with the spreadsheets used to solve them. A set of exercises in chemical kinetics appropriate for an introductory course in physical chemistry is given. Error propagation calculations with experimental data are illustrated. (CW)

  20. A comparative study of chemical kinetics models for HMX in mesoscale simulations of shock initiation due to void collapse

    NASA Astrophysics Data System (ADS)

    Rai, Nirmal; Schweigert, Igor; Udaykumar, H. S.

    2015-06-01

    The development of chemical kinetics schemes for use in modeling the reactive mechanics of energetic materials such as HMX has been an active area of research. Decomposition, deflagration and detonation models need to predict time to ignition and locations of onset of chemical reaction in energetic materials when used in meso- and macro-scale simulations. Modeling the chemical processes and development of appropriate kinetic law is challenging work because of lack of experimental data. However, significant work has been done in this area. Multistep kinetic models by Tarver and Tran, Henson and Smilowitz have provided plausible chemical kinetic rate laws for HMX. These models vary in the way they model the details of the decomposition process. Hence, a comparative study of different models will provide an understanding of the uncertainties involved in predicting ignition in HMX. In the current work, hot-spot ignition due to void collapse in shock compressed HMX has been analyzed using several reaction rate models, including the Tarver-Tran 4-equation model, the Henson-Smilowitz 7-equation model, and a new rate model that combines the condensed-phase decomposition rates measured by Brill et al and the detailed mechanism of nitramine flame chemistry due to Yetter et al. The chemical models have been incorporated in a massively parallel Eulerian code SCIMITAR3D. The variations in the predicted thresholds due to differences in the rate models will be discussed.

  1. HCCI experiments with gasoline surrogate fuels modeled by a semidetailed chemical kinetic model

    SciTech Connect

    Andrae, J.C.G.; Head, R.A.

    2009-04-15

    Experiments in a homogeneous charge compression ignition (HCCI) engine have been conducted with four gasoline surrogate fuel blends. The pure components in the surrogate fuels consisted of n-heptane, isooctane, toluene, ethanol and diisobutylene and fuel sensitivities (RON-MON) in the fuel blends ranged from two to nine. The operating conditions for the engine were p{sub in}=0.1 and 0.2 MPa, T{sub in}=80 and 250 C, {phi}=0.25 in air and engine speed 1200 rpm. A semidetailed chemical kinetic model (142 species and 672 reactions) for gasoline surrogate fuels, validated against ignition data from experiments conducted in shock tubes for gasoline surrogate fuel blends at 1.0{<=} p{<=}5.0MPa, 700{<=} T{<=}1200 K and {phi}=1.0, was successfully used to qualitatively predict the HCCI experiments using a single zone modeling approach. The fuel blends that had higher fuel sensitivity were more resistant to autoignition for low intake temperature and high intake pressure and less resistant to autoignition for high intake temperature and low intake pressure. A sensitivity analysis shows that at high intake temperature the chemistry of the fuels ethanol, toluene and diisobutylene helps to advance ignition. This is consistent with the trend that fuels with the least Negative Temperature Coefficient (NTC) behavior show the highest octane sensitivity, and become less resistant to autoignition at high intake temperatures. For high intake pressure the sensitivity analysis shows that fuels in the fuel blend with no NTC behavior consume OH radicals and acts as a radical scavenger for the fuels with NTC behavior. This is consistent with the observed trend of an increase in RON and fuel sensitivity. With data from shock tube experiments in the literature and HCCI modeling in this work, a correlation between the reciprocal pressure exponent on the ignition delay to the fuel sensitivity and volume percentage of single-stage ignition fuel in the fuel blend was found. Higher fuel

  2. Modeling of simulated photochemical smog with kinetic mechanisms. Volume 2. Chemk: a computer modeling scheme for chemical kinetics. Final report, July 1978-September 1979

    SciTech Connect

    Whitten, G.Z.; Hogo, H.

    1980-02-01

    Mechanisms that describe the formation of photochemical smog are developed using a computer modeling technique directed toward the simulation of data collected in two smog chambers: an indoor chamber and a dual outdoor chamber. Individual compounds for which specific experiments were simulated and mechanisms developed include the following: formaldehyde, acetaldehyde, ethylene, propylene, butane, and toluene. Volume 2 contains the user's manual and coding for a chemical kinetics computer program, CHEMK.

  3. Chemical kinetic model of hydrocarbon generation, expulsion, and destruction applied to the Maracaibo basin, Venezuela

    SciTech Connect

    Sweeney, J.J.; Braun, R.L.; Burnham, A.K.

    1995-10-01

    This paper describes the development and application of a compositional chemical model of hydrocarbon generation, expulsion,a nd destruction for the Cretaceous La Luna Formation source rock of the Maraciabo basin, Venezuela. Applications include both laboratory and geological settings. Laboratory pyrolysis experiments were used to study bulk oil generation, expulsion, and associated changes in composition of the kerogen, extractable organic matter, and expelled and unexpelled hydrocarbons. The laboratory experiments were also used to determine kinetic parameters to quantitatively describe organic reactions, via a computer model that also includes simulation of pressure-driven primary expulsion, over widely varying conditions. We show that the chemical model accuratley simulates the experimental results. Thermal history models for wells in the Maraciabo basin were used to simulate hydrocarbon generation and pore pressure development in the La Luna Formation and expulsion into nearby Cretaceous reservoirs. Results of the modeling indicate that both compaction disequilibrium and organic maturation play important roles in the development of excess pore pressure in the La Luna Formation. The model simulation of the variation of indicators such as Rock-Eval parameters and extract and oil compositions shows generally good agreement with measurements from remaining kerogen, oils, and extracts recovered from the La Luna Formation and from nearby Cretaceous reservoirs.

  4. Insitu Measurements and Modeling of Carbon Nanotube Array Growth Kinetics during Chemical Vapor Deposition

    SciTech Connect

    Puretzky, Alexander A; Geohegan, David B; Jesse, Stephen; Ivanov, Ilia N; Eres, Gyula

    2005-01-01

    Direct measurements of carbon nanotube growth kinetics are described based upon time-resolved reflectivity (TRR) of a HeNe laser beam from vertically aligned nanotube arrays (VANTAs) as they grow during chemical vapor deposition (CVD). Growth rates and terminal lengths were measured in situ for VANTAs growing during CVD between 535 C and 900 C on Si substrates with evaporated Al/Fe/Mo multi-layered catalysts and acetylene feedstock at different feedstock partial pressures. Methods of analysis of the TRR signals are presented to interpret catalyst particle formation and oxidation, as well as the porosity of the VANTAs. A rate-equation model is developed to describe the measured kinetics in terms of activation energies and rate constants for surface carbon formation and diffusion on the catalyst nanoparticle, nanotube growth, and catalyst over-coating. Taken together with the TRR data, this model enables basic understanding and optimization of growth conditions for any catalyst/feedstock combination. The model lends insight into the main processes responsible for the growth of VANTAs, the measured number of walls in the nanotubes at different temperatures, conditions for growth of single-wall carbon nanotube arrays, and likely catalyst poisoning mechanisms responsible for the sharp decline in growth rates observed at high temperatures.

  5. A Sequential Fluid-mechanic Chemical-kinetic Model of Propane HCCI Combustion

    SciTech Connect

    Aceves, S M; Flowers, D L; Martinez-Frias, J; Smith, J R; Westbrook, C; Pitz, W; Dibble, R; Wright, J F; Akinyemi, W C; Hessel, R P

    2000-11-29

    We have developed a methodology for predicting combustion and emissions in a Homogeneous Charge Compression Ignition (HCCI) Engine. This methodology combines a detailed fluid mechanics code with a detailed chemical kinetics code. Instead of directly linking the two codes, which would require an extremely long computational time, the methodology consists of first running the fluid mechanics code to obtain temperature profiles as a function of time. These temperature profiles are then used as input to a multi-zone chemical kinetics code. The advantage of this procedure is that a small number of zones (10) is enough to obtain accurate results. This procedure achieves the benefits of linking the fluid mechanics and the chemical kinetics codes with a great reduction in the computational effort, to a level that can be handled with current computers. The success of this procedure is in large part a consequence of the fact that for much of the compression stroke the chemistry is inactive and thus has little influence on fluid mechanics and heat transfer. Then, when chemistry is active, combustion is rather sudden, leaving little time for interaction between chemistry and fluid mixing and heat transfer. This sequential methodology has been capable of explaining the main characteristics of HCCI combustion that have been observed in experiments. In this paper, we use our model to explore an HCCI engine running on propane. The paper compares experimental and numerical pressure traces, heat release rates, and hydrocarbon and carbon monoxide emissions. The results show an excellent agreement, even in parameters that are difficult to predict, such as chemical heat release rates. Carbon monoxide emissions are reasonably well predicted, even though it is intrinsically difficult to make good predictions of CO emissions in HCCI engines. The paper includes a sensitivity study on the effect of the heat transfer correlation on the results of the analysis. Importantly, the paper also

  6. Modeling Multiphase Chemical Kinetics of OH Radical Reacting with Biomass Burning Organic Aerosol

    NASA Astrophysics Data System (ADS)

    Arangio, Andrea; Slade, Jonathan H.; Berkemeier, Thomas; Knopf, Daniel A.; Shiraiwa, Manabu

    2014-05-01

    Levoglucosan, abietic acid and nitroguaiacol are commonly used as molecular tracers of biomass burning in source apportionment. Recent studies have demonstrated the decay of levoglucosan when the particles were exposed to atmospherically relevant concentration of OH radicals [1-3]. However, multiphase chemical kinetics of OH radical reacting with such compounds has not fully understood. Here we apply the kinetic multi-layer model for gas-particle interactions (KM-GAP) [4] to experimental data of OH exposure to levoglucosan, abietic acid and nitroguaiacol [1]. KM-GAP resolves the following mass transport and chemical reactions explicitly: gas-phase diffusion, reversible surface adsorption, surface reaction, surface-bulk transport, bulk diffusion and reaction. The particle shrink due to the evaporation of volatile reaction products is also considered. The time- and concentration-dependence of reactive uptake coefficient of OH radicals were simulated by KM-GAP. The measured OH uptake coefficients were fitted by a Monte Carlo (MC) filtering coupled with a genetic algorithm (GA) to derive physicochemical parameters such as bulk diffusion coefficient, Henry's law coefficient and desorption lifetime of OH radicals. We assessed the relative contribution of surface and bulk reactions to the overall uptake of OH radicals. Chemical half-life and the evaporation time scale of these compounds are estimated in different scenarios (dry, humid and cloud processing conditions) and at different OH concentrations. REFERENCES [1] J. H. Slade, D. A. Knopf, Phys. Chem. Chem. Phys., 2013, 15, 5898. [2] S. H. Kessler, J. D. Smith, D.L. Che, D.R. Worsnop, K. R. Wilson, J. H. Kroll, Environ. Sci. Technol., 2010, 44, 7005. [3] C. J. Hennigan, A. P. Sullivan, J. L. Collett Jr, A. L. Robinson, Geophys. Res. Lett., 2010, 37, L09806. [4] M. Shiraiwa, C. Pfrang, T. Koop, U. Pöschl, Atmos. Chem. Phys, 2012, 12, 2777.

  7. Time-resolved simplified chemical kinetics modelling using computational singular perturbation

    NASA Technical Reports Server (NTRS)

    Lam, S. H.; Goussis, D. A.; Konopka, D.

    1989-01-01

    A CO-CH4-air reaction system is used to demonstrate the computational singular perturbation (CSP) method for deriving time-resolved simplified chemical kinetics models. CSP provides a programmable algorithm to group the given collection of elementary reactions into reaction groups which are ordered according to their speed. The concept of Importance Index k(m)exp s is introduced: k(m)exp s is defined to be a number between 0 and 1 which measures the importance of the m-th reaction group to the s-th reactant and can readily be computed from data generated by CSP. It is suggested that the robustness of the solutions of the reaction system can be qualitatively assessed by inspecting the Importance Index data.

  8. Chemical kinetics parameters and model validation for the gasification of PCEA nuclear graphite

    SciTech Connect

    El-Genk, Mohamed S; Tournier, Jean-Michel; Contescu, Cristian I

    2014-01-01

    A series of gasification experiments, using two right cylinder specimens (~ 12.7 x 25.4 mm and 25.4 x 25.4 mm) of PCEA nuclear graphite in ambient airflow, measured the total gasification flux at weight losses up to 41.5% and temperatures (893-1015 K) characteristics of those for in-pores gasification Mode (a) and in-pores diffusion-limited Mode (b). The chemical kinetics parameters for the gasification of PCEA graphite are determined using a multi-parameters optimization algorithm from the measurements of the total gasification rate and transient weight loss in experiments. These parameters are: (i) the pre-exponential rate coefficients and the Gaussian distributions and values of specific activation energies for adsorption of oxygen and desorption of CO gas; (ii) the specific activation energy and pre-exponential rate coefficient for the breakup of stable un-dissociated C(O2) oxygen radicals to form stable (CO) complexes; (iii) the specific activation energy and pre-exponential coefficient for desorption of CO2 gas and; (iv) the initial surface area of reactive free sites per unit mass. This area is consistently 13.5% higher than that for nuclear graphite grades of NBG-25 and IG-110 and decreases inversely proportional with the square root of the initial mass of the graphite specimens in the experiments. Experimental measurements successfully validate the chemical-reactions kinetics model that calculates continuous Arrhenius curves of the total gasification flux and the production rates of CO and CO2 gases. The model results at different total weight losses agree well with measurements and expand beyond the temperatures in the experiments to the diffusion-limited mode of gasification. Also calculated are the production rates of CO and CO2 gases and their relative contributions to the total gasification rate in the experiments as functions of temperature, for total weight losses of 5% and 10%.

  9. Asymptotic expansions of slow invariant manifolds and reduction of chemical kinetics models

    NASA Astrophysics Data System (ADS)

    Sobolev, V. A.; Tropkina, E. A.

    2012-01-01

    Methods of the geometric theory of singular perturbations are used to reduce the dimensions of problems in chemical kinetics. The methods are based on using slow invariant manifolds. As a result, the original system is replaced by one on an invariant manifold, whose dimension coincides with that of the slow subsystem. Explicit and implicit representations of slow invariant manifolds are applied. The mathematical apparatus described is used to develop N.N. Semenov's fundamental ideas related to the method of quasi-stationary concentrations and is used to study particular problems in chemical kinetics.

  10. Dominant particles and reactions in a two-temperature chemical kinetic model of a decaying SF6 arc

    NASA Astrophysics Data System (ADS)

    Wang, Xiaohua; Gao, Qingqing; Fu, Yuwei; Yang, Aijun; Rong, Mingzhe; Wu, Yi; Niu, Chunping; Murphy, Anthony B.

    2016-03-01

    This paper is devoted to the computation of the non-equilibrium composition of an SF6 plasma, and determination of the dominant particles and reactions, at conditions relevant to high-voltage circuit breakers after current zero (temperatures from 12 000 K to 1000 K and a pressure of 4 atm). The non-equilibrium composition is characterized by departures from both thermal and chemical equilibrium. In thermal non-equilibrium process, the electron temperature (T e) is not equal to the heavy-particle temperature (T h), while for chemical non-equilibrium, a chemical kinetic model is adopted. In order to evaluate the reasonableness and reliability of the non-equilibrium composition, calculation methods for equilibrium composition based on Gibbs free energy minimization and kinetic composition in a one-temperature kinetic model are first considered. Based on the one-temperature kinetic model, a two-temperature kinetic model with the ratio T e/T h varying as a function of the logarithm of electron density ratio (n e/n\\text{e}\\max ) was established. In this model, T* is introduced to allow a smooth transition between T h and T e and to determine the temperatures for the rate constants. The initial composition in the kinetic models is obtained from the asymptotic composition as infinite time is approached at 12 000 K. The molar fractions of neutral particles and ions in the two-temperature kinetic model are consistent with the equilibrium composition and the composition obtained from the one-temperature kinetic model above 10 000 K, while significant differences appear below 10 000 K. Based on the dependence of the particle distributions on temperature in the two-temperature kinetic model, three temperature ranges, and the dominant particles and reactions in the respective ranges, are determined. The full model is then simplified into three models and the accuracy of the simplified models is assessed. The simplified models reduce the number of species and

  11. Acceleration of the KINETICS Integrated Dynamical/Chemical Computational Model Using MPI

    NASA Technical Reports Server (NTRS)

    Grossman, Max; Willacy, Karen; Allen, Mark

    2011-01-01

    Understanding the evolution of a planet's atmosphere not only provides a better theoretical understanding of planetary physics and the formation of planets, but also grants useful insight into Earth's own atmosphere. One of the tools used at JPL for the modeling of planetary atmospheres and protostellar disks is KINETICS. KINETICS can simulate years of complex dynamics and chemistry.

  12. Chemical kinetics on extrasolar planets.

    PubMed

    Moses, Julianne I

    2014-04-28

    Chemical kinetics plays an important role in controlling the atmospheric composition of all planetary atmospheres, including those of extrasolar planets. For the hottest exoplanets, the composition can closely follow thermochemical-equilibrium predictions, at least in the visible and infrared photosphere at dayside (eclipse) conditions. However, for atmospheric temperatures approximately <2000K, and in the uppermost atmosphere at any temperature, chemical kinetics matters. The two key mechanisms by which kinetic processes drive an exoplanet atmosphere out of equilibrium are photochemistry and transport-induced quenching. I review these disequilibrium processes in detail, discuss observational consequences and examine some of the current evidence for kinetic processes on extrasolar planets. PMID:24664912

  13. History and Philosophy of Science through Models: The Case of Chemical Kinetics.

    ERIC Educational Resources Information Center

    Justi, Rosaria; Gilbert, John K.

    1999-01-01

    A greater role for the history and philosophy of science in science education can only be realized if it is based on both a credible analytical approach--such as that of Lakatos--and if the evolution of a sufficient number of major themes in science is known in suitable detail. Considers chemical kinetics as an example topic. Contains 62…

  14. A new methodology to determine kinetic parameters for one- and two-step chemical models

    NASA Technical Reports Server (NTRS)

    Mantel, T.; Egolfopoulos, F. N.; Bowman, C. T.

    1996-01-01

    In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.

  15. CFD analysis of municipal solid waste combustion using detailed chemical kinetic modelling.

    PubMed

    Frank, Alex; Castaldi, Marco J

    2014-08-01

    Nitrogen oxides (NO x ) emissions from the combustion of municipal solid waste (MSW) in waste-to-energy (WtE) facilities are receiving renewed attention to reduce their output further. While NO x emissions are currently 60% below allowed limits, further reductions will decrease the air pollution control (APC) system burden and reduce consumption of NH3. This work combines the incorporation of the GRI 3.0 mechanism as a detailed chemical kinetic model (DCKM) into a custom three-dimensional (3D) computational fluid dynamics (CFD) model fully to understand the NO x chemistry in the above-bed burnout zones. Specifically, thermal, prompt and fuel NO formation mechanisms were evaluated for the system and a parametric study was utilized to determine the effect of varying fuel nitrogen conversion intermediates between HCN, NH3 and NO directly. Simulation results indicate that the fuel nitrogen mechanism accounts for 92% of the total NO produced in the system with thermal and prompt mechanisms accounting for the remaining 8%. Results also show a 5% variation in final NO concentration between HCN and NH3 inlet conditions, demonstrating that the fuel nitrogen intermediate assumed is not significant. Furthermore, the conversion ratio of fuel nitrogen to NO was 0.33, revealing that the majority of fuel nitrogen forms N2. PMID:25005043

  16. A reduced chemical kinetic model for IC engine combustion simulations with primary reference fuels

    SciTech Connect

    Ra, Youngchul; Reitz, Rolf D.

    2008-12-15

    A reduced chemical kinetic mechanism for the oxidation of primary reference fuel (PRF) has been developed and applied to model internal combustion engines. Starting from an existing reduced reaction mechanism for n-heptane oxidation, a new reduced n-heptane mechanism was generated by including an additional five species and their relevant reactions, by updating the reaction rate constants of several reactions pertaining to oxidation of carbon monoxide and hydrogen, and by optimizing reaction rate constants of selected reactions. Using a similar approach, a reduced mechanism for iso-octane oxidation was built and combined with the n-heptane mechanism to form a PRF mechanism. The final version of the PRF mechanism consists of 41 species and 130 reactions. Validation of the present PRF mechanism was performed with measurements from shock tube tests, and HCCI and direct injection engine experiments available in the literature. The results show that the present PRF mechanism gives reliable performance for combustion predictions, as well as computational efficiency improvements for multidimensional CFD simulations. (author)

  17. A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions

    SciTech Connect

    Zhang, Fan; Yeh, Gour-Tsyh; Parker, Jack C; Brooks, Scott C; Pace, Molly; Kim, Young Jin; Jardine, Philip M; Watson, David B

    2007-01-01

    This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing NE equilibrium reactions and a set of reactive transport equations of M-NE kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.

  18. Premixed hydrocarbon stagnation flames : experiments and simulations to validate combustion chemical-kinetic models

    NASA Astrophysics Data System (ADS)

    Benezech, Laurent Jean-Michel

    A methodology based on the comparison of flame simulations relying on reacting flow models with experiment is applied to C1-C3 stagnation flames. The work reported targets the assessment and validation of the modeled reactions and reaction rates relevant to (C1-C3)-flame propagation in several detailed combustion kinetic models. A concensus does not, as yet, exist on the modeling of the reasonably well-understood oxidation of C1-C2 flames, and a better knowledge of C3 hydrocarbon combustion chemistry is required before attempting to bridge the gap between the oxidation of C1-C2 hydrocarbons and the more complex chemistry of heavier hydrocarbons in a single kinetic model. Simultaneous measurements of velocity and CH-radical profiles were performed in atmospheric propane(C3H8)- and propylene(C3H6)-air laminar premixed stagnation flames stabilized in a jet-wall configuration. These nearly-flat flames can be modeled by one-dimensional simulations, providing a means to validate kinetic models. Experimental data for these C3 flames and similar experimental data for atmospheric methane(CH4)-, ethane(C2H6)-, and ethylene(C2H4)-air flames are compared to numerical simulations performed with a one-dimensional hydrodynamic model, a multi-component transport formulation including thermal diffusion, and different detailed-chemistry models, in order to assess the adequacy of the models employed. A novel continuation technique between kinetic models was developed and applied successfully to obtain solutions with the less-robust models. The 2005/12 and 2005/10 releases of the San Diego mechanism are found to have the best overall performance in C3H8 and C3H6 flames, and in CH4, C2H6, and C2H4 flames, respectively. Flame position provides a good surrogate for flame speed in stagnation-flow stabilized flames. The logarithmic sensitivities of the simulated flame locations to variations in the kinetic rates are calculated via the "brute-force" method for fifteen representative flames

  19. A Chemical Kinetic Modeling Study of the Effects of Oxygenated Hydrocarbons on Soot Emissions from Diesel Engines

    SciTech Connect

    Westbrook, C K; Pitz, W J; Curran, H J

    2005-11-14

    A detailed chemical kinetic modeling approach is used to examine the phenomenon of suppression of sooting in diesel engines by addition of oxygenated hydrocarbon species to the fuel. This suppression, which has been observed experimentally for a few years, is explained kinetically as a reduction in concentrations of soot precursors present in the hot products of a fuel-rich diesel ignition zone when oxygenates are included. Oxygenates decrease the overall equivalence ratio of the igniting mixture, producing higher ignition temperatures and more radical species to consume more soot precursor species, leading to lower soot production. The kinetic model is also used to show how different oxygenates, ester structures in particular, can have different soot-suppression efficiencies due to differences in molecular structure of the oxygenated species.

  20. Effects of correlated parameters and uncertainty in electronic-structure-based chemical kinetic modelling.

    PubMed

    Sutton, Jonathan E; Guo, Wei; Katsoulakis, Markos A; Vlachos, Dionisios G

    2016-04-01

    Kinetic models based on first principles are becoming common place in heterogeneous catalysis because of their ability to interpret experimental data, identify the rate-controlling step, guide experiments and predict novel materials. To overcome the tremendous computational cost of estimating parameters of complex networks on metal catalysts, approximate quantum mechanical calculations are employed that render models potentially inaccurate. Here, by introducing correlative global sensitivity analysis and uncertainty quantification, we show that neglecting correlations in the energies of species and reactions can lead to an incorrect identification of influential parameters and key reaction intermediates and reactions. We rationalize why models often underpredict reaction rates and show that, despite the uncertainty being large, the method can, in conjunction with experimental data, identify influential missing reaction pathways and provide insights into the catalyst active site and the kinetic reliability of a model. The method is demonstrated in ethanol steam reforming for hydrogen production for fuel cells. PMID:27001728

  1. Effects of correlated parameters and uncertainty in electronic-structure-based chemical kinetic modelling

    NASA Astrophysics Data System (ADS)

    Sutton, Jonathan E.; Guo, Wei; Katsoulakis, Markos A.; Vlachos, Dionisios G.

    2016-04-01

    Kinetic models based on first principles are becoming common place in heterogeneous catalysis because of their ability to interpret experimental data, identify the rate-controlling step, guide experiments and predict novel materials. To overcome the tremendous computational cost of estimating parameters of complex networks on metal catalysts, approximate quantum mechanical calculations are employed that render models potentially inaccurate. Here, by introducing correlative global sensitivity analysis and uncertainty quantification, we show that neglecting correlations in the energies of species and reactions can lead to an incorrect identification of influential parameters and key reaction intermediates and reactions. We rationalize why models often underpredict reaction rates and show that, despite the uncertainty being large, the method can, in conjunction with experimental data, identify influential missing reaction pathways and provide insights into the catalyst active site and the kinetic reliability of a model. The method is demonstrated in ethanol steam reforming for hydrogen production for fuel cells.

  2. Thermochemical analysis and kinetics aspects for a chemical model for camphene ozonolysis

    NASA Astrophysics Data System (ADS)

    Oliveira, R. C. de M.; Bauerfeldt, G. F.

    2012-10-01

    In this work, a chemical model for the camphene ozonolysis, leading to carbonyl final products, is proposed and discussed on the basis of the thermochemical properties and kinetic data obtained at density functional theory levels of calculation. The mechanism is initiated by the electrophilic attack of ozone to the double bond in camphene leading to a 1,2,3-trioxolane intermediate, which decomposes to peroxy radicals and carbonyl compounds in a total of 10 elementary reactions. The thermodynamic properties (enthalpy and entropies differences) are calculated at 298 K. For the thermochemical evaluation, theoretical calculations are performed with the B3LYP, MPW1PW91, and mPW1K density functionals and the basis sets 6-31G(d), 6-31G(2d,2p), 6-31+G(d,p), and 6-31+G(2d,2p). Eventually, single point calculations adopting the 6-311++G(2d,2p) basis set are performed in order to improve the electronic energies. The enthalpy profiles suggest highly exothermic reactions for the individual steps, with a global enthalpy difference of -179.18 kcal mol-1, determined at the B3LYP/6-31+G(2d,2p) level. The Gibbs free energy differences for each step, at 298 K, calculated at the B3LYP/6-311++G(2d,2p)//B3LYP/6-31+G(2d,2p) level, are used to estimate the composition of a final product mixture under equilibrium conditions as 58% of camphenilone and 42% of 6,6-dimethyl-ɛ-caprolactone-2,5-methylene. For the reaction kinetics, the bimolecular O3 + camphene step is assumed to be rate determining in the global mechanism. A saddle point for the ozone addition to the double bond is located and rate constants are determined on the basis of the transition state theory. This saddle point is well represented by a loosely bound structure and corrections for the basis set superposition error (BSSE) are calculated, either by considering the effect over the geometry optimization procedure (here referred as CP1 procedure), or the effect of the BSSE over the electronic energy of a previously optimized

  3. Kinetic models of hydrocarbon generation

    SciTech Connect

    Burnham, A.K.; Sweeney, J.J.

    1990-10-25

    We are carrying out an integrated program of laboratory experiments, kinetics modeling, and basin thermal history modeling in order to better understand the natural breakdown of organic matter into oil and gas. Our kinetic models of organic maturation are being used to better understand the coupling of generation, cracking, expulsion, and overpressuring in both the laboratory and geologic setting. Currently we are carrying out chemical experiments and developing more efficient chemical kinetic modeling schemes to obtain a better understanding of expulsion and cracking from lean source rocks and from hydrogen-poor (terrestrial) organic source material. We verify the chemical kinetic models by integrating them with thermal history models of hydrocarbon-producing sediments and comparing predicted and observed characteristics of the hydrocarbon occurrence in a variety of settings. We intend to apply this approach to evaluate the potential for deep gas resources in the Pacific Northwest and in the Louisiana Gulf Coast. 11 refs., 4 figs.

  4. Modeling the kinetics of microbial degradation of deicing chemicals in porous media under flow conditions.

    PubMed

    Wehrer, Markus; Jaesche, Philipp; Totsche, Kai Uwe

    2012-09-01

    A quantitative knowledge of the fate of deicing chemicals in the subsurface can be provided by joint analysis of lab experiments with numerical simulation models. In the present study, published experimental data of microbial degradation of the deicing chemical propylene glycol (PG) under flow conditions in soil columns were simulated inversely to receive the parameters of degradation. We evaluated different scenarios of an advection-dispersion model including different terms for degradation, such as zero order, first order and inclusion of a growing and decaying biomass for their ability to explain the data. The general break-through behavior of propylene glycol in soil columns can be simulated well using a coupled model of solute transport and degradation with growth and decay of biomass. The susceptibility of the model to non-unique solutions was investigated using systematical forward and inverse simulations. We found that the model tends to equifinal solutions under certain conditions. PMID:22609860

  5. HCCI experiments with toluene reference fuels modeled by a semidetailed chemical kinetic model

    SciTech Connect

    Andrae, J.C.G.; Brinck, T.; Kalghatgi, G.T.

    2008-12-15

    A semidetailed mechanism (137 species and 633 reactions) and new experiments in a homogeneous charge compression ignition (HCCI) engine on the autoignition of toluene reference fuels are presented. Skeletal mechanisms for isooctane and n-heptane were added to a detailed toluene submechanism. The model shows generally good agreement with ignition delay times measured in a shock tube and a rapid compression machine and is sensitive to changes in temperature, pressure, and mixture strength. The addition of reactions involving the formation and destruction of benzylperoxide radical was crucial to modeling toluene shock tube data. Laminar burning velocities for benzene and toluene were well predicted by the model after some revision of the high-temperature chemistry. Moreover, laminar burning velocities of a real gasoline at 353 and 500 K could be predicted by the model using a toluene reference fuel as a surrogate. The model also captures the experimentally observed differences in combustion phasing of toluene/n-heptane mixtures, compared to a primary reference fuel of the same research octane number, in HCCI engines as the intake pressure and temperature are changed. For high intake pressures and low intake temperatures, a sensitivity analysis at the moment of maximum heat release rate shows that the consumption of phenoxy radicals is rate-limiting when a toluene/n-heptane fuel is used, which makes this fuel more resistant to autoignition than the primary reference fuel. Typical CPU times encountered in zero-dimensional calculations were on the order of seconds and minutes in laminar flame speed calculations. Cross reactions between benzylperoxy radicals and n-heptane improved the model predictions of shock tube experiments for {phi}=1.0 and temperatures lower than 800 K for an n-heptane/toluene fuel mixture, but cross reactions had no influence on HCCI simulations. (author)

  6. The mathematical properties of the quasi-chemical model for microorganism growth-death kinetics in foods.

    PubMed

    Ross, E W; Taub, I A; Doona, C J; Feeherry, F E; Kustin, K

    2005-03-15

    Knowledge of the mathematical properties of the quasi-chemical model [Taub, Feeherry, Ross, Kustin, Doona, 2003. A quasi-chemical kinetics model for the growth and death of Staphylococcus aureus in intermediate moisture bread. J. Food Sci. 68 (8), 2530-2537], which is used to characterize and predict microbial growth-death kinetics in foods, is important for its applications in predictive microbiology. The model consists of a system of four ordinary differential equations (ODEs), which govern the temporal dependence of the bacterial life cycle (the lag, exponential growth, stationary, and death phases, respectively). The ODE system derives from a hypothetical four-step reaction scheme that postulates the activity of a critical intermediate as an antagonist to growth (perhaps through a quorum sensing biomechanism). The general behavior of the solutions to the ODEs is illustrated by several examples. In instances when explicit mathematical solutions to these ODEs are not obtainable, mathematical approximations are used to find solutions that are helpful in evaluating growth in the early stages and again near the end of the process. Useful solutions for the ODE system are also obtained in the case where the rate of antagonist formation is small. The examples and the approximate solutions provide guidance in the parameter estimation that must be done when fitting the model to data. The general behavior of the solutions is illustrated by examples, and the MATLAB programs with worked examples are included in the appendices for use by predictive microbiologists for data collected independently. PMID:15734564

  7. Incorporation of a Chemical Kinetics Model for Composition B in a Parallel Finite-Element Algorithm

    NASA Astrophysics Data System (ADS)

    Kallman, Elizabeth; Pauler, Denise

    2009-06-01

    A thermal degradation model for Composition B (Comp B) explosive is being evaluated for incorporation into a finite-element algorithm [1]. The RDX component of Comp B dominates the thermal degradation since its decomposition process occurs at lower temperatures than TNT. The model assumes that solid and liquid RDX decompose by the same mechanisms, but along different reaction pathways [2, 3]. A steady-state approximation is applied to the gaseous intermediates and is compared to the full transient analysis for the entire reaction scheme. The parallel finite-element algorithm is used to predict the pressure increase on the interior of the metal casing of confined Comp B due to the production of gases during thermal decomposition. =0pt References [1] E. M. Kallman, ``Scalable Cluster-Based Galerkin Analysis for Kinetics Models of Energetic Materials,'' SIAM CSE, March 2-6, 2009. [2] D. K. Zerkle, ``Composition B Decomposition and Ignition Model,'' 13th International Detonation Symposium, July 23-28, 2006. [3] J. M. Zucker, A. J. Barra, D. K. Zerkle, M. J. Kaneshige and P. M. Dickson, ``Thermal Decomposition Models for High Explosive Compositions,'' 14th APS Topical Conference on Shock Compression of Condensed Matter, July 31-August 5, 2005.

  8. Quantum logics and chemical kinetics

    NASA Astrophysics Data System (ADS)

    Ivanov, C. I.

    1981-06-01

    A statistical theory of chemical kinetics is presented based on the quantum logical concept of chemical observables. The apparatus of Boolean algebra B is applied for the construction of appropriate composition polynomials referring to any stipulated arrangement of the atomic constituents. A physically motivated probability measure μ( F) is introduced on the field B of chemical observables, which considers the occurrence of the yes response of a given F ɛ B. The equations for the time evolution of the species density operators and the master equations for the corresponding number densities are derived. The general treatment is applied to a superposition of elementary substitution reactions (AB) α + C ⇄ (AC) β + B. The expressions for the reaction rate coefficients are established.

  9. Chemical Kinetic Modeling of Dimethyl Carbonate in an Opposed-Flow Diffusion Flame

    SciTech Connect

    Glaude, P A; Pitz, W J; Thomson, M J

    2003-12-08

    Dimethyl carbonate (DMC) has been of interest as an oxygenate additive to diesel fuel because of its high oxygen content. In this study, a chemical kinetic mechanism for DMC was developed for the first time and used to understand its combustion under conditions in an opposed flow diffusion flame. Computed results were compared to experimental results from an opposed flow diffusion flame. It was found that the decomposition rate DMC {yields} H{sub 3}COC(=O)O. + CH{sub 3} in the flame was much slower than originally thought because resonance stabilization in the H{sub 3}COC(=O)O. radical was less than expected. Also, a new molecular elimination path for DMC is proposed and its rate calculated by quantum chemical methods. In the simulations of DMC in the flame, it was determined that much of the oxygen in dimethyl carbonate goes directly to CO{sub 2}. This characteristic indicates that DMC would not be an effective oxygenate additive for reducing soot emissions from diesel engines. In an ideal oxygenate additive for diesel fuel, each oxygen atom stays bonded to one carbon atom in the products thereby preventing the formation of carbon-carbon bonds that can lead to soot. When CO2 is formed directly, two oxygen atoms are bonded to one carbon atom thereby wasting one oxygen atom in the oxygenate additive. To determine how much CO{sub 2} is formed directly, the branching ratio of the key reaction, CH{sub 3}OC.=O going to the products CH{sub 3} + CO{sub 2} or CH{sub 3}O + CO was determined by ab initio methods. The A-factors of the rate constant of this reaction were found to be about 20 times higher than previous factors estimates. The new reaction rate constants obtained can be used as reaction rate rules for all oxygenates that contain the ester moiety including biodiesel.

  10. Kinetic-quantum chemical model for catalytic cycles: the Haber-Bosch process and the effect of reagent concentration.

    PubMed

    Kozuch, Sebastian; Shaik, Sason

    2008-07-01

    A combined kinetic-quantum chemical model is developed with the goal of estimating in a straightforward way the turnover frequency (TOF) of catalytic cycles, based on the state energies obtained by quantum chemical calculations. We describe how the apparent activation energy of the whole cycle, so-called energetic span (delta E), is influenced by the energy levels of two species: the TOF determining transition state (TDTS) and the TOF determining intermediate (TDI). Because these key species need not be adjoining states, we conclude that for catalysis there are no rate-determining steps, only rate determining states. In addition, we add here the influence of reactants concentrations. And, finally, the model is applied to the Haber-Bosch process of ammonia synthesis, for which we show how to calculate which catalyst will be the most effective under specific reagents conditions. PMID:18537227

  11. Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

    NASA Technical Reports Server (NTRS)

    Bittker, D. A.

    1993-01-01

    A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

  12. Mathematical description of complex chemical kinetics and application to CFD modeling codes

    NASA Technical Reports Server (NTRS)

    Bittker, D. A.

    1993-01-01

    A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

  13. Kinetics of microbial degradation of deicing chemicals in percolated porous media - the modeling perspective

    NASA Astrophysics Data System (ADS)

    Wehrer, Markus; Lissner, Heidi; Totsche, Kai

    2013-04-01

    A quantitative knowledge of the fate of deicing chemicals in the subsurface can be provided by analysis of laboratory and field experiments with numerical simulation models. In the present study, experimental data of microbial degradation of the deicing chemical propylene glycol (PG) under flow conditions in soil columns and field lysimeters were simulated to analyze the process conditions of degradation and to obtain the according parameters. Results from the column experiment were evaluated applying different scenarios of an advection-dispersion model using HYDRUS-1D. To reconstruct the data, different competing degradation models were included, i.e., zero order, first order and inclusion of a growing and decaying biomass. The general breakthrough behavior of propylene glycol in soil columns can be simulated well using a coupled model of solute transport and degradation with growth and decay of biomass. The susceptibility of the model to non-unique solutions was investigated using systematical forward and inverse simulations. We found that the model tends to equifinal solutions under certain conditions. Complex experimental boundary conditions can help to avoid this. Under field conditions, the situation is far more complex than in the laboratory. Studying the fate of PG with undisturbed lysimeters we found that aerobic and anaerobic degradation occurs simultaneously. We attribute this to the physical structure and the aggregated nature of the undisturbed soil material . This results in the presence of spatially disjoint oxidative and reductive regions of microbial activity and requires, but is not fully reflected by a dual porosity model. Currently, the numerical simulation of this system is in progress, considering several flow and transport models. A stochastic global search algorithm (DREAM-ZS) is used in conjuction with HYDRUS-1D to avoid local minima in the inverse simulations. The study shows the current limitations and potentials of modeling degradation

  14. Kinetics study of solid ammonia borane hydrogen release--modeling and experimental validation for chemical hydrogen storage.

    PubMed

    Choi, Young Joon; Rönnebro, Ewa C E; Rassat, Scot; Karkamkar, Abhi; Maupin, Gary; Holladay, Jamie; Simmons, Kevin; Brooks, Kriston

    2014-05-01

    Ammonia borane (AB), NH3BH3, is a promising material for chemical hydrogen storage with 19.6 wt% gravimetric hydrogen capacity of which maximum 16.2 wt% hydrogen can be released via an exothermic thermal decomposition below 200 °C. We have investigated the kinetics of hydrogen release from AB and from an AB-methyl cellulose (AB/MC) composite at temperatures of 160-300 °C using both experiments and modeling. The hydrogen release rate at 300 °C is twice as fast as at 160 °C. The purpose of our study was to show safe hydrogen release without thermal runaway effects and to validate system model kinetics. AB/MC released hydrogen at ∼20 °C lower than neat AB and at a faster release rate in that temperature range. Based on the experimental results, the kinetics equations were revised to better represent the growth and nucleation process during decomposition of AB. We explored two different reactor concepts; auger and fixed bed. The current auger reactor concept turned out to not be appropriate, however, we demonstrated safe self-propagation of the hydrogen release reaction of solid AB/MC in a fixed bed reactor. PMID:24647627

  15. Kinetics Study of Solid Ammonia Borane Hydrogen Release – Modeling and Experimental Validation for Chemical Hydrogen Storage

    SciTech Connect

    Choi, Yong-Joon; Ronnebro, Ewa; Rassat, Scot D.; Karkamkar, Abhijeet J.; Maupin, Gary D.; Holladay, Jamelyn D.; Simmons, Kevin L.; Brooks, Kriston P.

    2014-02-24

    Ammonia borane (AB), NH3BH3, is a promising material for chemical hydrogen storage with 19.6 wt% gravimetric hydrogen capacity of which 16.2 wt% hydrogen can be utilized below 200°C. We have investigated the kinetics of hydrogen release from AB and from an AB-methyl cellulose (AB/MC) composite at temperatures of 160-300°C using both experiments and modeling. The purpose of our study was to show safe hydrogen release without thermal runaway effects and to validate system model kinetics. AB/MC released hydrogen at ~20°C lower than neat AB and at a rate that is two times faster. Based on the experimental results, the kinetics equations were revised to better represent the growth and nucleation process during decomposition of AB. We explored two different reactor concepts; Auger and fixed bed. The current Auger reactor concept turned out to not be appropriate, however, we demonstrated safe self-propagation of the hydrogen release reaction of solid AB/MC in a fixed bed reactor.

  16. Detailed chemical kinetic models for large n-alkanes and iso-alkanes found in conventional and F-T diesel fuels

    SciTech Connect

    Westbrook, C K; Pitz, W J; Mehl, M; Curran, H J

    2009-03-09

    n-Hexadecane and 2,2,4,4,6,8,8-heptamethylnonane represent the primary reference fuels for diesel that are used to determine cetane number, a measure of the ignition property of diesel fuel. With the development of chemical kinetics models for both primary reference fuels, a new capability is now available to model diesel fuel ignition. Additionally, we have developed chemical kinetic models for a whole series of large n-alkanes and a large iso-alkane to represent these chemical classes in fuel surrogates for conventional and future fuels. These chemical kinetic models are used to predict the effect of the aforementioned fuel components on ignition characteristics under conditions found in internal combustion engines.

  17. Fast Prediction of HCCI and PCCI Combustion with an Artificial Neural Network-Based Chemical Kinetic Model

    SciTech Connect

    Piggott, W T; Aceves, S M; Flowers, D L; Chen, J Y

    2007-09-26

    We have added the capability to look at in-cylinder fuel distributions using a previously developed ignition model within a fluid mechanics code (KIVA3V) that uses an artificial neural network (ANN) to predict ignition (The combined code: KIVA3V-ANN). KIVA3V-ANN was originally developed and validated for analysis of Homogeneous Charge Compression Ignition (HCCI) combustion, but it is also applicable to the more difficult problem of Premixed Charge Compression Ignition (PCCI) combustion. PCCI combustion refers to cases where combustion occurs as a nonmixing controlled, chemical kinetics dominated, autoignition process, where the fuel, air, and residual gas mixtures are not necessarily as homogeneous as in HCCI combustion. This paper analyzes the effects of introducing charge non-uniformity into a KIVA3V-ANN simulation. The results are compared to experimental results, as well as simulation results using a more physically representative and computationally intensive code (KIVA3V-MPI-MZ), which links a fluid mechanics code to a multi-zone detailed chemical kinetics solver. The results indicate that KIVA3V-ANN produces reasonable approximations to the more accurate KIVA3V-MPI-MZ at a much reduced computational cost.

  18. TACK—a program coupling chemical kinetics with a two-dimensional transport model in geochemical systems

    NASA Astrophysics Data System (ADS)

    Källvenius, Göran; Ekberg, Christian

    2003-05-01

    The Transport And Chemical Kinetics (TACK) program has been designed to make predictions of the chemistry in the vicinity of a planned repository for nuclear waste, i.e. SFL 3-5, where SFL is the Swedish abbreviation for "Swedish repository for long-lived waste". This implies modelling transport and chemistry in fractured rock. The system concerned in the modelling of SFL is leaching water from decommissioning waste in concrete. The concrete will raise the pH in the water to between 12 and 13.5. So far, only a few calculations have been made on such systems. Coupled transport and chemical reaction programs should be used, since the system is important for safety assessments of the repository. At least two of programs can be used for this kind of problem, for example OS3D/GIMRT and PHAST. As it is also important to consider the uncertainty of the model, the TACK program fills an important purpose here. A slightly different approach to the problem may give significantly different results. Because validation is generally not possible, using several programs is the only key to identifying conceptual uncertainties. To illustrate this point, comparative calculations have been made between TACK and the PHAST program. The calculations gave qualitatively similar result but quantitatively somewhat differing results. The TACK program couples the well known PHREEQC geochemical program with a two-dimensional transport model. The PHREEQC calculations include speciation of solutions and mineral reactions involving kinetics. The reasons for choosing this program are that it is quite a general one and is relatively stable at the high pH values present in the systems used. The transport phenomena taken into account in the model are advection, diffusion and dispersion in two dimensions.

  19. Modelling cycle to cycle variations in an SI engine with detailed chemical kinetics

    SciTech Connect

    Etheridge, Jonathan; Mosbach, Sebastian; Kraft, Markus; Wu, Hao; Collings, Nick

    2011-01-15

    This paper presents experimental results and a new computational model that investigate cycle to cycle variations (CCV) in a spark ignition (SI) engine. An established stochastic reactor model (SRM) previously used to examine homogeneous charge compression ignition (HCCI) combustion has been extended by spark initiation, flame propagation and flame termination sub-models in order to simulate combustion in SI engines. The model contains a detailed chemical mechanism but relatively short computation times are achieved. The flame front is assumed to be spherical and centred at the spark location, and a pent roof and piston bowl geometry are accounted for. The model is validated by simulating the pressure profile and emissions from an iso-octane fuelled single cylinder research engine that showed low CCV. The effects of key parameters are investigated. Experimental results that show cycle to cycle fluctuations in a four-cylinder naturally aspirated gasoline fuelled SI engine are presented. The model is then coupled with GT-Power, a one-dimensional engine simulation tool, which is used to simulate the breathing events during a multi-cycle simulation. This allows an investigation of the cyclic fluctuations in peak pressure. The source and magnitude of nitric oxide (NO) emissions produced by different cycles are then investigated. It was found that faster burning cycles result in increased NO emissions compared with cycles that have a slower rate of combustion and that more is produced in the early stages of combustion compared with later in the cycle. The majority of NO was produced via the thermal mechanism just after combustion begins. (author)

  20. Inflation Rates, Car Devaluation, and Chemical Kinetics.

    ERIC Educational Resources Information Center

    Pogliani, Lionello; Berberan-Santos, Mario N.

    1996-01-01

    Describes the inflation rate problem and offers an interesting analogy with chemical kinetics. Presents and solves the car devaluation problem as a normal chemical kinetic problem where the order of the rate law and the value of the rate constant are derived. (JRH)

  1. Enhancing Thai Students' Learning of Chemical Kinetics

    ERIC Educational Resources Information Center

    Chairam, Sanoe; Somsook, Ekasith; Coll, Richard K.

    2009-01-01

    Chemical kinetics is an extremely important concept for introductory chemistry courses. The literature suggests that instruction in chemical kinetics is often teacher-dominated at both the secondary school and tertiary levels, and this is the case in Thailand--the educational context for this inquiry. The work reported here seeks to shift students…

  2. Updated Chemical Kinetics and Sensitivity Analysis Code

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, Krishnan

    2005-01-01

    An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

  3. Methyl Formate Oxidation: Speciation Data, Laminar Burning Velocities, Ignition Delay Times and a Validated Chemical Kinetic Model

    SciTech Connect

    Dooley, S.; Burke, M. P.; Chaos, M.; Stein, Y.; Dryer, F. L.; Zhukov, V. P.; Finch, O.; Simmie, J. M.; Curran, H. J.

    2010-07-16

    The oxidation of methyl formate (CH{sub 3}OCHO) has been studied in three experimental environments over a range of applied combustion relevant conditions: 1. A variable-pressure flow reactor has been used to quantify reactant, major intermediate and product species as a function of residence time at 3 atm and 0.5% fuel concentration for oxygen/fuel stoichiometries of 0.5, 1.0, and 1.5 at 900 K, and for pyrolysis at 975 K. 2. Shock tube ignition delays have been determined for CH{sub 3}OCHO/O{sub 2}/Ar mixtures at pressures of ≈ 2.7, 5.4, and 9.2 atm and temperatures of 1275–1935 K for mixture compositions of 0.5% fuel (at equivalence ratios of 1.0, 2.0, and 0.5) and 2.5% fuel (at an equivalence ratio of 1.0). 3. Laminar burning velocities of outwardly propagating spherical CH{sub 3}OCHO/air flames have been determined for stoichiometries ranging from 0.8–1.6, at atmospheric pressure using a pressure-release-type high-pressure chamber. A detailed chemical kinetic model has been constructed, validated against, and used to interpret these experimental data. The kinetic model shows that methyl formate oxidation proceeds through concerted elimination reactions, principally forming methanol and carbon monoxide as well as through bimolecular hydrogen abstraction reactions. The relative importance of elimination versus abstraction was found to depend on the particular environment. In general, methyl formate is consumed exclusively through molecular decomposition in shock tube environments, while at flow reactor and freely propagating premixed flame conditions, there is significant competition between hydrogen abstraction and concerted elimination channels. It is suspected that in diffusion flame configurations the elimination channels contribute more significantly than in premixed environments.

  4. Comprehensive chemical kinetic modeling of the oxidation of C8 and larger n-alkanes and 2-methylalkanes

    SciTech Connect

    Sarathy, S M; Westbrook, C K; Pitz, W J; Mehl, M; Togbe, C; Dagaut, P; Wang, H; Oehlschlaeger, M; NIemann, U; Seshadri, K; Veloo, P S; Ji, C; Egolfopoulos, F; Lu, T

    2011-03-16

    Conventional petroleum jet and diesel fuels, as well as alternative Fischer-Tropsch (FT) fuels and hydrotreated renewable jet (HRJ) fuels, contain high molecular weight lightly branched alkanes (i.e., methylalkanes) and straight chain alkanes (n-alkanes). Improving the combustion of these fuels in practical applications requires a fundamental understanding of large hydrocarbon combustion chemistry. This research project presents a detailed and reduced chemical kinetic mechanism for singly methylated iso-alkanes (i.e., 2-methylalkanes) ranging from C{sub 8} to C{sub 20}. The mechanism also includes an updated version of our previously published C{sub 8} to C{sub 16} n-alkanes model. The complete detailed mechanism contains approximately 7,200 species 31,400 reactions. The proposed model is validated against new experimental data from a variety of fundamental combustion devices including premixed and nonpremixed flames, perfectly stirred reactors and shock tubes. This new model is used to show how the presence of a methyl branch affects important combustion properties such as laminar flame propagation, ignition, and species formation.

  5. Experimental and chemical kinetic modeling study of small methyl esters oxidation: Methyl (E)-2-butenoate and methyl butanoate

    SciTech Connect

    Gail, S.; Sarathy, S.M.; Thomson, M.J.; Dievart, P.; Dagaut, P.

    2008-12-15

    This study examines the effect of unsaturation on the combustion of fatty acid methyl esters (FAME). New experimental results were obtained for the oxidation of methyl (E)-2-butenoate (MC, unsaturated C{sub 4} FAME) and methyl butanoate (MB, saturated C{sub 4} FAME) in a jet-stirred reactor (JSR) at atmospheric pressure under dilute conditions over the temperature range 850-1400 K, and two equivalence ratios ({phi}=0.375,0.75) with a residence time of 0.07 s. The results consist of concentration profiles of the reactants, stable intermediates, and final products, measured by probe sampling followed by on-line and off-line gas chromatography analyses. The oxidation of MC and MB in the JSR and under counterflow diffusion flame conditions was modeled using a new detailed chemical kinetic reaction mechanism (301 species and 1516 reactions) derived from previous schemes proposed in the literature. The laminar counterflow flame and JSR (for {phi}=1.13) experimental results used were from a previous study on the comparison of the combustion of both compounds. Sensitivity analyses and reaction path analyses, based on rates of reaction, were used to interpret the results. The data and the model show that MC has reaction pathways analogous to that of MB under the present conditions. The model of MC oxidation provides a better understanding of the effect of the ester function on combustion, and the effect of unsaturation on the combustion of fatty acid methyl ester compounds typically found in biodiesel. (author)

  6. A comprehensive experimental and detailed chemical kinetic modelling study of 2,5-dimethylfuran pyrolysis and oxidation

    PubMed Central

    Somers, Kieran P.; Simmie, John M.; Gillespie, Fiona; Conroy, Christine; Black, Gráinne; Metcalfe, Wayne K.; Battin-Leclerc, Frédérique; Dirrenberger, Patricia; Herbinet, Olivier; Glaude, Pierre-Alexandre; Dagaut, Philippe; Togbé, Casimir; Yasunaga, Kenji; Fernandes, Ravi X.; Lee, Changyoul; Tripathi, Rupali; Curran, Henry J.

    2013-01-01

    The pyrolytic and oxidative behaviour of the biofuel 2,5-dimethylfuran (25DMF) has been studied in a range of experimental facilities in order to investigate the relatively unexplored combustion chemistry of the title species and to provide combustor relevant experimental data. The pyrolysis of 25DMF has been re-investigated in a shock tube using the single-pulse method for mixtures of 3% 25DMF in argon, at temperatures from 1200–1350 K, pressures from 2–2.5 atm and residence times of approximately 2 ms. Ignition delay times for mixtures of 0.75% 25DMF in argon have been measured at atmospheric pressure, temperatures of 1350–1800 K at equivalence ratios (ϕ) of 0.5, 1.0 and 2.0 along with auto-ignition measurements for stoichiometric fuel in air mixtures of 25DMF at 20 and 80 bar, from 820–1210 K. This is supplemented with an oxidative speciation study of 25DMF in a jet-stirred reactor (JSR) from 770–1220 K, at 10.0 atm, residence times of 0.7 s and at ϕ = 0.5, 1.0 and 2.0. Laminar burning velocities for 25DMF-air mixtures have been measured using the heat-flux method at unburnt gas temperatures of 298 and 358 K, at atmospheric pressure from ϕ = 0.6–1.6. These laminar burning velocity measurements highlight inconsistencies in the current literature data and provide a validation target for kinetic mechanisms. A detailed chemical kinetic mechanism containing 2768 reactions and 545 species has been simultaneously developed to describe the combustion of 25DMF under the experimental conditions described above. Numerical modelling results based on the mechanism can accurately reproduce the majority of experimental data. At high temperatures, a hydrogen atom transfer reaction is found to be the dominant unimolecular decomposition pathway of 25DMF. The reactions of hydrogen atom with the fuel are also found to be important in predicting pyrolysis and ignition delay time experiments. Numerous proposals are made on the mechanism and kinetics of the previously

  7. A comprehensive experimental and detailed chemical kinetic modelling study of 2,5-dimethylfuran pyrolysis and oxidation.

    PubMed

    Somers, Kieran P; Simmie, John M; Gillespie, Fiona; Conroy, Christine; Black, Gráinne; Metcalfe, Wayne K; Battin-Leclerc, Frédérique; Dirrenberger, Patricia; Herbinet, Olivier; Glaude, Pierre-Alexandre; Dagaut, Philippe; Togbé, Casimir; Yasunaga, Kenji; Fernandes, Ravi X; Lee, Changyoul; Tripathi, Rupali; Curran, Henry J

    2013-11-01

    The pyrolytic and oxidative behaviour of the biofuel 2,5-dimethylfuran (25DMF) has been studied in a range of experimental facilities in order to investigate the relatively unexplored combustion chemistry of the title species and to provide combustor relevant experimental data. The pyrolysis of 25DMF has been re-investigated in a shock tube using the single-pulse method for mixtures of 3% 25DMF in argon, at temperatures from 1200-1350 K, pressures from 2-2.5 atm and residence times of approximately 2 ms. Ignition delay times for mixtures of 0.75% 25DMF in argon have been measured at atmospheric pressure, temperatures of 1350-1800 K at equivalence ratios (ϕ) of 0.5, 1.0 and 2.0 along with auto-ignition measurements for stoichiometric fuel in air mixtures of 25DMF at 20 and 80 bar, from 820-1210 K. This is supplemented with an oxidative speciation study of 25DMF in a jet-stirred reactor (JSR) from 770-1220 K, at 10.0 atm, residence times of 0.7 s and at ϕ = 0.5, 1.0 and 2.0. Laminar burning velocities for 25DMF-air mixtures have been measured using the heat-flux method at unburnt gas temperatures of 298 and 358 K, at atmospheric pressure from ϕ = 0.6-1.6. These laminar burning velocity measurements highlight inconsistencies in the current literature data and provide a validation target for kinetic mechanisms. A detailed chemical kinetic mechanism containing 2768 reactions and 545 species has been simultaneously developed to describe the combustion of 25DMF under the experimental conditions described above. Numerical modelling results based on the mechanism can accurately reproduce the majority of experimental data. At high temperatures, a hydrogen atom transfer reaction is found to be the dominant unimolecular decomposition pathway of 25DMF. The reactions of hydrogen atom with the fuel are also found to be important in predicting pyrolysis and ignition delay time experiments. Numerous proposals are made on the mechanism and kinetics of the previously unexplored

  8. Detailed chemical kinetic models for large n-alkanes and iso-alkanes found in conventional and F-T diesel fuels

    SciTech Connect

    Westbrook, C K; Pitz, W J; Curran, H J; Mehl, M

    2008-12-15

    Detailed chemical kinetic models are needed to simulate the combustion of current and future transportation fuels. These models should represent the various chemical classes in these fuels. Conventional diesel fuels are composed of n-alkanes, iso-alkanes, cycloalkanes and aromatics (Farrell et al. 2007). For future fuels, there is a renewed interest in Fischer-Tropsch (F-T) processes which can be used to synthesize diesel and other transportation fuels from biomass, coal and natural gas. F-T diesel fuels are expected to be similar to F-T jet fuels which are commonly comprised of iso-alkanes with some n-alkanes (Smith and Bruno, 2008). Thus, n-alkanes and iso-alkanes are common chemical classes in these conventional and future fuels. This paper reports on the development of chemical kinetic models of large n-alkanes and iso-alkanes to represent these chemical classes in conventional and future fuels. Two large iso-alkanes are 2,2,4,4,6,8,8-heptamethylnonane, which is a primary reference fuel for diesel, and isooctane, a primary reference fuel for gasoline. Other iso-alkanes are branched alkanes with a single methyl side chain, typical of most F-T fuels. The chemical kinetic models are then used to predict the effect of these fuel components on ignition characteristics under conditions found in internal combustion engines.

  9. The Monte Carlo Integration Computer as an Instructional Model for the Simulation of Equilibrium and Kinetic Chemical Processes: The Development and Evaluation of a Teaching Aid.

    ERIC Educational Resources Information Center

    Wood, Dean Arthur

    A special purpose digital computer which utilizes the Monte Carlo integration method of obtaining simulations of chemical processes was developed and constructed. The computer, designated as the Monte Carlo Integration Computer (MCIC), was designed as an instructional model for the illustration of kinetic and equilibrium processes, and was…

  10. Kinetic-fluid dynamics modeling of I{sub 2} dissociation in supersonic chemical oxygen-iodine lasers

    SciTech Connect

    Waichman, K.; Barmashenko, B. D.; Rosenwaks, S.

    2009-09-15

    The mechanism of I{sub 2} dissociation in supersonic chemical oxygen-iodine lasers (COILs) is studied applying kinetic-fluid dynamics modeling, where pathways involving the excited species I{sub 2}(X {sup 1}SIGMA{sub g}{sup +},10<=v<25), I{sub 2}(X {sup 1}SIGMA{sub g}{sup +},25<=v<=47), I{sub 2}(A{sup '} {sup 3}PI{sub 2u}), I{sub 2}(A {sup 3}PI{sub 1u}), O{sub 2}(X {sup 3}SIGMA{sub g}{sup -},v), O{sub 2}(a {sup 1}DELTA{sub g},v), O{sub 2}(b {sup 1}SIGMA{sub g}{sup +},v), and I({sup 2}P{sub 1/2}) as intermediate reactants are included. The gist of the model is adding the first reactant and reducing the contribution of the second as compared to previous models. These changes, recently suggested by Azyazov, et al. [J. Chem. Phys. 130, 104306 (2009)], significantly improve the agreement with the measurements of the gain in a low pressure supersonic COIL for all I{sub 2} flow rates that have been tested in the experiments. In particular, the lack of agreement for high I{sub 2} flow rates, which was encountered in previous models, has been eliminated in the present model. It is suggested that future modeling of the COIL operation should take into account the proposed contribution of the above mentioned reactants.

  11. Diesel Surrogate Fuels for Engine Testing and Chemical-Kinetic Modeling: Compositions and Properties

    PubMed Central

    Mueller, Charles J.; Cannella, William J.; Bays, J. Timothy; Bruno, Thomas J.; DeFabio, Kathy; Dettman, Heather D.; Gieleciak, Rafal M.; Huber, Marcia L.; Kweon, Chol-Bum; McConnell, Steven S.; Pitz, William J.; Ratcliff, Matthew A.

    2016-01-01

    The primary objectives of this work were to formulate, blend, and characterize a set of four ultralow-sulfur diesel surrogate fuels in quantities sufficient to enable their study in single-cylinder-engine and combustion-vessel experiments. The surrogate fuels feature increasing levels of compositional accuracy (i.e., increasing exactness in matching hydrocarbon structural characteristics) relative to the single target diesel fuel upon which the surrogate fuels are based. This approach was taken to assist in determining the minimum level of surrogate-fuel compositional accuracy that is required to adequately emulate the performance characteristics of the target fuel under different combustion modes. For each of the four surrogate fuels, an approximately 30 L batch was blended, and a number of the physical and chemical properties were measured. This work documents the surrogate-fuel creation process and the results of the property measurements. PMID:27330248

  12. Chemical kinetics and photochemical data for use in stratospheric modeling. Evaluation number 6

    NASA Technical Reports Server (NTRS)

    Demore, W. B.; Molina, M. J.; Watson, R. T.; Golden, D. M.; Hampson, R. F.; Kurylo, M. J.; Howard, C. J.; Ravishankara, A. R.

    1983-01-01

    Evaluated sets of rate constants and photochemical cross sections are presented. The primary application of the data is in the modeling of stratospheric processes, with particular emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena.

  13. Chemical kinetics and photochemical data for use in stratospheric modeling: Evaluation number 5

    NASA Technical Reports Server (NTRS)

    Demore, W. B.

    1982-01-01

    Sets of rate constants and photochemical cross sections compiled which were evaluated. The primary application of the data is in the modeling of stratospheric processes on the ozone layer and its possible perturbation by anthropogenic and natural phenomena are emphasized.

  14. Determination of Chemical Kinetic Rate Constants of a Model for Carbothermal Processing of Lunar Regolith Simulant Using Methane

    NASA Technical Reports Server (NTRS)

    Balasubramaniam, R; Gokoglu, S.; Hegde, U.

    2009-01-01

    We have previously developed a chemical conversion model of the carbothermal processing of lunar regolith using methane to predict the rate of production of carbon monoxide. In this carbothermal process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Hydrogen is carried away by the exiting stream of gases and carbon is deposited on the melt surface. The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. In our model, we assume that the flux of carbon deposited is equal to the product of the surface reaction rate constant gamma and the concentration of methane adjacent to the melt surface. Similarly, the rate of consumption of carbon per unit volume in the melt is equal to the product of the melt reaction rate constant k and the concentrations of carbon and metal oxide in the melt. In this paper, we describe our effort to determine gamma and k by comparison of the predictions from our model with test data obtained by ORBITEC (Orbital Technologies Corporation). The concentration of methane adjacent to the melt surface is a necessary input to the model. It is inferred from the test data by a mass balance of methane, adopting the usual assumptions of the continuously-stirred-tank-reactor model, whereby the average concentration of a given gaseous species equals its exit concentration. The reaction rates gamma and k have been determined by a non-linear least-squares fit to the test data for the production of carbon monoxide and the fraction of the incoming methane that is converted. The comparison of test data with our model predictions using the determined chemical kinetic rate constants provides a consistent interpretation of the process over the full range of temperatures, pressures, and methane flow rates used in the tests, thereby increasing our confidence to use the model for scale-up purposes.

  15. Chemical kinetics of geminal recombination

    SciTech Connect

    Levin, P.P.; Khudyakov, I.V.; Brin, E.F.; Kuz'min, V.A.

    1988-09-01

    The kinetics of geminal recombination of triplet radical pairs formed in photoreduction of benzophenone by p-cresol in glycerin solution was studied by pulsed laser photolysis. The experiments were conducted at several temperatures and in a constant magnetic field of H = 0.34 T. The parameters in six kinetic equations describing geminal recombination were determined with a computer. The values of the sums of the squares of the residual deviations of the approximation were obtained. It was found that the kinetics are best described by the functions proposed by Noyes and Shushin. It was shown that it is necessary to use the mutual diffusion coefficient of the radicals, which is significantly smaller than the sum of the estimations of the experimental values of the radical diffusion coefficients, for describing the kinetics due to the correlations of the molecular motions of the radicals in the cage.

  16. DETAILED CHEMICAL KINETIC MODELING OF ISO-OCTANE SI-HCCI TRANSITION

    SciTech Connect

    Havstad, M A; Aceves, S M; McNenly, M J; Piggott, W T; Edwards, K D; Wagner, R M; Daw, C S; Finney, C A

    2009-10-12

    The authors describe a CHEMKIN-based multi-zone model that simulates the expected combustion variations in a single-cylinder engine fueled with iso-octane as the engine transitions from spark-ignited (ST) combustion to homogeneous charge compression ignition (HCCI) combustion. The model includes a 63-species reaction mechanism and mass and energy balances for the cylinder and the exhaust flow. For this study they assumed that the SI-to-HCCI transition is implemented by means of increasing the internal exhaust gas recirculation (EGR) at constant engine speed. This transition scneario is consistent with that implemented in previously reported experimental measurements on an experimental engine equipped with variable valve actuation. They find that the model captures many of the important experimental trends, including stable SI combustion at low EGR ({approx} 0.10), a transition to highly unstable combustion at intermediate EGR, and finally stable HCCI combustion at very high EGR ({approx} 0.75). Remaining differences between the predicted and experimental instability patterns indicate that there is further room for model improvement.

  17. Chemical Kinetics and Photochemical Data for Use in Stratospheric Modeling. Evaluation No. 12

    NASA Technical Reports Server (NTRS)

    DeMore, W. B.; Sander, S. P.; Golden, D. M.; Hampson, R. F.; Kurylo, M. J.; Howard, C. J.; Ravishankara, A. R.; Kolb, C. E.; Molina, M. J.

    1997-01-01

    This is the twelfth in a series of evaluated sets of rate constants and photochemical cross sections compiled by the NASA Panel for Data Evaluation. The primary application of the data is in the modeling of stratospheric processes, with special emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena.

  18. Detailed Chemical Kinetic Modeling of Iso-octane SI-HCCI Transition

    SciTech Connect

    Havstad, Mark A; Aceves, Salvador M; McNenly, Matthew J; Piggott, William T; Edwards, Kevin Dean; Wagner, Robert M; Daw, C Stuart; FINNEY, Charles E A

    2010-01-01

    We describe a CHEMKIN-based multi-zone model that simulates the expected combustion variations in a single-cylinder engine fueled with iso-octane as the engine transitions from spark-ignited (SI) combustion to homogenous charge compression ignition (HCCI) combustion. The model includes a 63-species reaction mechanism and mass and energy balances for the cylinder and the exhaust flow. For this study we assumed that the SI-to-HCCI transition is implemented by means of increasing the internal exhaust gas recirculation (EGR) at constant engine speed. This transition scenario is consistent with that implemented in previously reported experimental measurements on an experimental engine equipped with variable valve actuation. We find that the model captures many of the important experimental trends, including stable SI combustion at low EGR (-0.10), a transition to highly unstable combustion at intermediate EGR, and finally stable HCCI combustion at very high EGR (-0.75). Remaining differences between the predicted and experimental instability patterns indicate that there is further room for model improvement.

  19. Chemical kinetics and photochemical data for use in stratospheric modeling evaluation Number 8

    NASA Technical Reports Server (NTRS)

    Demore, W. B.; Molina, M. J.; Sander, S. P.; Golden, D. M.; Hampson, R. F.; Kurylo, M. J.; Howard, C. J.; Ravishankara, A. R.

    1987-01-01

    This is the eighth in a series of evaluated sets of rate constants and photochemical cross sections compiled by the NASA Panel for Data Evaluation. The primary application of the data is in the modeling of stratospheric processes, with particular emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena. Copies of this evaluation are available from the Jet Propulsion Laboratory, Documentation Section, 111-116B, California Institute of Technology, Pasadena, California, 91109.

  20. Detailed Chemical Kinetic Modeling of Surrogate Fuels for Gasoline and Application to an HCCI Engine

    SciTech Connect

    Naik, C V; Pitz, W J; Sj?berg, M; Dec, J E; Orme, J; Curran, H J; Simmie, J M; Westbrook, C K

    2005-01-07

    Gasoline consists of many different classes of hydrocarbons, such as paraffins, olefins, aromatics, and cycloalkanes. In this study, a surrogate gasoline reaction mechanism is developed, and it has one representative fuel constituent from each of these classes. These selected constituents are iso-octane, n-heptane, 1-pentene, toluene, and methyl-cyclohexane. The mechanism was developed in a step-wise fashion, adding submechanisms to treat each fuel component. Reactions important for low temperature oxidation (<1000K) and cross-reactions among different fuels are incorporated into the mechanism. The mechanism consists of 1214 species and 5401 reactions. A single-zone engine model is used to evaluate how well the mechanism captures autoignition behavior for conditions corresponding to homogeneous charge compression ignition (HCCI) engine operation. Experimental data are available for both how the combustion phasing changes with fueling at a constant intake temperature, and also how the intake temperature has to be changed with pressure in order to maintain combustion phasing for a fixed equivalence ratio. Three different surrogate fuel mixtures are used for the modeling. Predictions are in reasonably good agreement with the engine data. In addition, the heat release rate is calculated and compared to the data from experiments. The model predicts less low-temperature heat release than that measured. It is found that the low temperature heat-release rate depends strongly on engine speed, reactions of RO{sub 2}+HO{sub 2}, fuel composition, and pressure boost.

  1. Chemical kinetics and oil shale process design

    SciTech Connect

    Burnham, A.K.

    1993-07-01

    Oil shale processes are reviewed with the goal of showing how chemical kinetics influences the design and operation of different processes for different types of oil shale. Reaction kinetics are presented for organic pyrolysis, carbon combustion, carbonate decomposition, and sulfur and nitrogen reactions.

  2. Chemical kinetic mechanistic models to investigate cancer biology and impact cancer medicine

    NASA Astrophysics Data System (ADS)

    Stites, Edward C.

    2013-04-01

    Traditional experimental biology has provided a mechanistic understanding of cancer in which the malignancy develops through the acquisition of mutations that disrupt cellular processes. Several drugs developed to target such mutations have now demonstrated clinical value. These advances are unequivocal testaments to the value of traditional cellular and molecular biology. However, several features of cancer may limit the pace of progress that can be made with established experimental approaches alone. The mutated genes (and resultant mutant proteins) function within large biochemical networks. Biochemical networks typically have a large number of component molecules and are characterized by a large number of quantitative properties. Responses to a stimulus or perturbation are typically nonlinear and can display qualitative changes that depend upon the specific values of variable system properties. Features such as these can complicate the interpretation of experimental data and the formulation of logical hypotheses that drive further research. Mathematical models based upon the molecular reactions that define these networks combined with computational studies have the potential to deal with these obstacles and to enable currently available information to be more completely utilized. Many of the pressing problems in cancer biology and cancer medicine may benefit from a mathematical treatment. As work in this area advances, one can envision a future where such models may meaningfully contribute to the clinical management of cancer patients.

  3. Modeling the chemical kinetics of atmospheric plasma for cell treatment in a liquid solution

    NASA Astrophysics Data System (ADS)

    Kim, H. Y.; Lee, H. W.; Kang, S. K.; Wk. Lee, H.; Kim, G. C.; Lee, J. K.

    2012-07-01

    Low temperature atmospheric pressure plasmas have been known to be effective for living cell inactivation in a liquid solution but it is not clear yet which species are key factors for the cell treatment. Using a global model, we elucidate the processes through which pH level in the solution is changed from neutral to acidic after plasma exposure and key components with pH and air variation. First, pH level in a liquid solution is changed by He+ and He(21S) radicals. Second, O3 density decreases as pH level in the solution decreases and air concentration decreases. It can be a method of removing O3 that causes chest pain and damages lung tissue when the density is very high. H2O2, HO2, and NO radicals are found to be key factors for cell inactivation in the solution with pH and air variation.

  4. Fluid flow and chemical reaction kinetics in metamorphic systems

    SciTech Connect

    Lasaga, A.C.; Rye, D.M. )

    1993-05-01

    The treatment and effects of chemical reaction kinetics during metamorphism are developed along with the incorporation of fluid flow, diffusion, and thermal evolution. The interplay of fluid flow and surface reaction rates, the distinction between steady state and equilibrium, and the possible overstepping of metamorphic reactions are discussed using a simple analytic model. This model serves as an introduction to the second part of the paper, which develops a reaction model that solves the coupled temperature-fluid flow-chemical composition differential equations relevant to metamorphic processes. Consideration of stable isotopic evidence requires that such a kinetic model be considered for the chemical evolution of a metamorphic aureole. A general numerical scheme is discussed to handle the solution of the model. The results of this kinetic model allow us to reach several important conclusions regarding the factors controlling the chemical evolution of mineral assemblages during a metamorphic event. 41 refs., 19 figs., 5 tabs.

  5. Modeling the chemical kinetics of atmospheric plasma for cell treatment in a liquid solution

    SciTech Connect

    Kim, H. Y.; Kang, S. K.; Lee, H. Wk.; Lee, H. W.; Kim, G. C.; Lee, J. K.

    2012-07-15

    Low temperature atmospheric pressure plasmas have been known to be effective for living cell inactivation in a liquid solution but it is not clear yet which species are key factors for the cell treatment. Using a global model, we elucidate the processes through which pH level in the solution is changed from neutral to acidic after plasma exposure and key components with pH and air variation. First, pH level in a liquid solution is changed by He{sup +} and He(2{sup 1}S) radicals. Second, O{sub 3} density decreases as pH level in the solution decreases and air concentration decreases. It can be a method of removing O{sub 3} that causes chest pain and damages lung tissue when the density is very high. H{sub 2}O{sub 2}, HO{sub 2}, and NO radicals are found to be key factors for cell inactivation in the solution with pH and air variation.

  6. Mass Conservation and Chemical Kinetics.

    ERIC Educational Resources Information Center

    Barbara, Thomas M.; Corio, P. L.

    1980-01-01

    Presents a method for obtaining all mass conservation conditions implied by a given mechanism in which the conditions are used to simplify integration of the rate equations and to derive stoichiometric relations. Discusses possibilities of faulty inference of kinetic information from a given stoichiometry. (CS)

  7. Application of a Genetic Algorithm to the Optimization of Rate Constants in Chemical Kinetic Models for Combustion Simulation of HCCI Engines

    NASA Astrophysics Data System (ADS)

    Kim, Sang-Kyu; Ito, Kazuma; Yoshihara, Daisuke; Wakisaka, Tomoyuki

    For numerically predicting the combustion processes in homogeneous charge compression ignition (HCCI) engines, practical chemical kinetic models have been explored. A genetic algorithm (GA) has been applied to the optimization of the rate constants in detailed chemical kinetic models, and a detailed kinetic model (592 reactions) for gasoline reference fuels with arbitrary octane number between 60 and 100 has been obtained from the detailed reaction schemes for iso-octane and n-heptane proposed by Golovitchev. The ignition timing in a gasoline HCCI engine has been predicted reasonably well by zero-dimensional simulation using the CHEMKIN code with this detailed kinetic model. An original reduced reaction scheme (45 reactions) for dimethyl ether (DME) has been derived from Curran’s detailed scheme, and the combustion process in a DME HCCI engine has been predicted reasonably well in a practical computation time by three-dimensional simulation using the authors’ GTT code, which has been linked to the CHEMKIN subroutines with the proposed reaction scheme and also has adopted a modified eddy dissipation combustion model.

  8. Chemical kinetics of Estane aging in PBX

    SciTech Connect

    Pack, R.T.; Hanson, D.E.; Redondo, A.

    1997-12-01

    The Plastic-Bonded Explosive PBX 9501 is about 95% HMX, 2.5% Estane 5703, 2.5% nitroplasticizer (NP), and 0.1% stabilizer by weight. The NP, BDNPA/F, is a eutectic mixture of bis(2,2-dinitropropyl) acetal and bis(2,2-dinitropropyl)formal. The stabilizer is diphenylamine (DPA) or Irganox 1010. The Estane, a polyester-polyurethane, slowly degrades with time. Knowledge of the effect of the Estane aging on the mechanical properties of the PBX 9501 is required to predict with confidence the useful lifetime of the explosive with respect to safety and reliability. A detailed master equation model of the chemical mechanisms and kinetics of the aging of Estane 5703 in PBX 9501 is being developed. Its output will be used as input into other models being developed to calculate the changes in the mechanical properties of the PBX.

  9. Multidimensional reactor kinetics modeling

    SciTech Connect

    Diamond, D.J.

    1996-11-01

    There is general agreement that for many light water reactor transient calculations, it is-necessary to use a multidimensional neutron kinetics model coupled to a thermal-hydraulics model for satisfactory results. These calculations are needed for a variety of applications for licensing safety analysis, probabilistic risk assessment (PRA), operational support, and training. The latter three applications have always required best-estimate models, but in the past applications for licensing could be satisfied with relatively simple models. By using more sophisticated best-estimate models, the consequences of these calculations are better understood, and the potential for gaining relief from restrictive operating limits increases. Hence, for all of the aforementioned applications, it is important to have the ability to do best-estimate calculations with multidimensional neutron kinetics models. coupled to sophisticated thermal-hydraulic models. Specifically, this paper reviews the status of multidimensional neutron kinetics modeling which would be used in conjunction with thermal-hydraulic models to do core dynamics calculations, either coupled to a complete NSSS representation or in isolation. In addition, the paper makes recommendations as to what should be the state-of-the-art for the next ten years. The review is an update to a previous review of the status as of ten years ago. The general requirements for a core dynamics code and the modeling available for such a code, discussed in that review, are still applicable. The emphasis in the current review is on the neutron kinetics assuming that the necessary thermal-hydraulic capability exists. In addition to discussing the basic neutron kinetics, discussion is given of related modeling (other than thermal- hydraulics). The capabilities and limitations of current computer codes are presented to understand the state-of-the-art and to help clarify the future direction of model development in this area.

  10. S3 and S4 abundances and improved chemical kinetic model for the lower atmosphere of Venus

    NASA Astrophysics Data System (ADS)

    Krasnopolsky, Vladimir A.

    2013-07-01

    Mixing ratios of S3 and S4 are obtained from reanalysis of the spectra of true absorption in the visible range retrieved by Maiorov et al. (Maiorov, B.S. et al. [2005]. Solar Syst. Res. 39, 267-282) from the Venera 11 observations. These mixing ratios are fS3 = 11 ± 3 ppt at 3-10 km and 18 ± 3 ppt at 10-19 km, fS4 = 4 ± 4 ppt at 3-10 km and 6 ± 2 ppt at 10-19 km, and show a steep decrease in both S3 and S4 above 19 km. Photolysis rates of S3 and S4 at various altitudes are calculated using the Venera 11 spectra and constant photolysis yields as free parameters. The chemical kinetic model for the Venus lower atmosphere (Krasnopolsky, V.A. [2007]. Icarus 191, 25-37) has been improved by inclusion of the S4 cycle from Yung et al. (Yung, Y.L. et al. [2009]. J. Geophys. Res. 114, E00B34), reduction of the H2SO4 and CO fluxes at the upper boundary of 47 km by a factor of 4 in accord with the recent photochemical models for the middle atmosphere, by using a closed lower boundary for OCS instead of a free parameter for this species at the surface, and some minor updates. Our model with the S4 cycle but without the SO3 + 2 OCS reaction suggested by Krasnopolsky and Pollack (Krasnopolsky, V.A., Pollack, J.B. [1994]. Icarus 109, 58-78) disagrees with the observations of OCS, CO, S3, and S4. However, inclusion of the S4 cycle improves the model fit to all observational constraints. The best-fit activation energy of 7800 K for thermolysis of S4 supports the S4 enthalpy from Mills (Mills, K.C. [1974]. Thermodynamic Data for Inorganic Sulfides, Selenides and Tellurides. Butterworths, London). Chemistry of the Venus lower atmosphere is initiated by disequilibrium products H2SO4 and CO from the middle atmosphere, photolysis of S3 and S4, and thermochemistry in the lowest scale height. The chemistry is mostly driven by sulfur that is formed in a slow reaction SO + SO, produces OCS, and results in dramatic changes in abundances of OCS, CO, and free sulfur allotropes. The SX + OCS

  11. Trajectory Calculations in Chemical Kinetics

    ERIC Educational Resources Information Center

    Hemphill, Gregory L.; White, John M.

    1972-01-01

    This exercise, suitable for an advanced undergraduate physical chemistry lab, examines the detailed theoretical description of a chemical reaction. Mathematical techniques of moderate complexity serve to introduce some aspects of theoretical chemistry. (Author/TS)

  12. Chemical, physical, and theoretical kinetics of an ultrafast folding protein

    PubMed Central

    Kubelka, Jan; Henry, Eric R.; Cellmer, Troy; Hofrichter, James; Eaton, William A.

    2008-01-01

    An extensive set of equilibrium and kinetic data is presented and analyzed for an ultrafast folding protein—the villin subdomain. The equilibrium data consist of the excess heat capacity, tryptophan fluorescence quantum yield, and natural circular-dichroism spectrum as a function of temperature, and the kinetic data consist of time courses of the quantum yield from nanosecond-laser temperature-jump experiments. The data are well fit with three kinds of models—a three-state chemical-kinetics model, a physical-kinetics model, and an Ising-like theoretical model that considers 105 possible conformations (microstates). In both the physical-kinetics and theoretical models, folding is described as diffusion on a one-dimensional free-energy surface. In the physical-kinetics model the reaction coordinate is unspecified, whereas in the theoretical model, order parameters, either the fraction of native contacts or the number of native residues, are used as reaction coordinates. The validity of these two reaction coordinates is demonstrated from calculation of the splitting probability from the rate matrix of the master equation for all 105 microstates. The analysis of the data on site-directed mutants using the chemical-kinetics model provides information on the structure of the transition-state ensemble; the physical-kinetics model allows an estimate of the height of the free-energy barrier separating the folded and unfolded states; and the theoretical model provides a detailed picture of the free-energy surface and a residue-by-residue description of the evolution of the folded structure, yet contains many fewer adjustable parameters than either the chemical- or physical-kinetics models. PMID:19033473

  13. Modelling of the physico-chemical behaviour of clay minerals with a thermo-kinetic model taking into account particles morphology in compacted material.

    NASA Astrophysics Data System (ADS)

    Sali, D.; Fritz, B.; Clément, C.; Michau, N.

    2003-04-01

    Modelling of fluid-mineral interactions is largely used in Earth Sciences studies to better understand the involved physicochemical processes and their long-term effect on the materials behaviour. Numerical models simplify the processes but try to preserve their main characteristics. Therefore the modelling results strongly depend on the data quality describing initial physicochemical conditions for rock materials, fluids and gases, and on the realistic way of processes representations. The current geo-chemical models do not well take into account rock porosity and permeability and the particle morphology of clay minerals. In compacted materials like those considered as barriers in waste repositories, low permeability rocks like mudstones or compacted powders will be used : they contain mainly fine particles and the geochemical models used for predicting their interactions with fluids tend to misjudge their surface areas, which are fundamental parameters in kinetic modelling. The purpose of this study was to improve how to take into account the particles morphology in the thermo-kinetic code KINDIS and the reactive transport code KIRMAT. A new function was integrated in these codes, considering the reaction surface area as a volume depending parameter and the calculated evolution of the mass balance in the system was coupled with the evolution of reactive surface areas. We made application exercises for numerical validation of these new versions of the codes and the results were compared with those of the pre-existing thermo-kinetic code KINDIS. Several points are highlighted. Taking into account reactive surface area evolution during simulation modifies the predicted mass transfers related to fluid-minerals interactions. Different secondary mineral phases are also observed during modelling. The evolution of the reactive surface parameter helps to solve the competition effects between different phases present in the system which are all able to fix the chemical

  14. Kinetic Modeling of Biological Systems

    PubMed Central

    Petzold, Linda; Pettigrew, Michel F.

    2010-01-01

    The dynamics of how the constituent components of a natural system interact defines the spatio-temporal response of the system to stimuli. Modeling the kinetics of the processes that represent a biophysical system has long been pursued with the aim of improving our understanding of the studied system. Due to the unique properties of biological systems, in addition to the usual difficulties faced in modeling the dynamics of physical or chemical systems, biological simulations encounter difficulties that result from intrinsic multiscale and stochastic nature of the biological processes. This chapter discusses the implications for simulation of models involving interacting species with very low copy numbers, which often occur in biological systems and give rise to significant relative fluctuations. The conditions necessitating the use of stochastic kinetic simulation methods and the mathematical foundations of the stochastic simulation algorithms are presented. How the well-organized structural hierarchies often seen in biological systems can lead to multiscale problems, and possible ways to address the encountered computational difficulties are discussed. We present the details of the existing kinetic simulation methods, and discuss their strengths and shortcomings. A list of the publicly available kinetic simulation tools and our reflections for future prospects are also provided. PMID:19381542

  15. Inflation Rates, Car Devaluation, and Chemical Kinetics

    NASA Astrophysics Data System (ADS)

    Pogliani, Lionello; Berberan-Santos, Màrio N.

    1996-10-01

    The inflation rate problem of a modern economy shows quite interesting similarities with chemical kinetics and especially with first-order chemical reactions. In fact, capital devaluation during periods of rather low inflation rates or inflation measured over short periods shows a dynamics formally similar to that followed by first-order chemical reactions and they can thus be treated by the aid of the same mathematical formalism. Deviations from this similarity occurs for higher inflation rates. The dynamics of price devaluation for two different types of car, a compact car and a luxury car, has been followed for seven years long and it has been established that car devaluation is a process that is formally similar to a zeroth-order chemical kinetic process disregarding the type of car, if car devaluation is much faster than money devaluation. In fact, expensive cars devaluate with a faster rate than inexpensive cars.

  16. The applications of chemical thermodynamics and chemical kinetics to planetary atmospheres research

    NASA Technical Reports Server (NTRS)

    Fegley, Bruce, Jr.

    1990-01-01

    A review of the applications of chemical thermodynamics and chemical kinetics to planetary atmospheres research during the past four decades is presented with an emphasis on chemical equilibrium models and thermochemical kinetics. Several current problems in planetary atmospheres research such as the origin of the atmospheres of the terrestrial planets, atmosphere-surface interactions on Venus and Mars, deep mixing in the atmospheres of the gas giant planets, and the origin of the atmospheres of outer planet satellites all require laboratory data on the kinetics of thermochemical reactions for their solution.

  17. Chemical Dosing and First-Order Kinetics

    ERIC Educational Resources Information Center

    Hladky, Paul W.

    2011-01-01

    College students encounter a variety of first-order phenomena in their mathematics and science courses. Introductory chemistry textbooks that discuss first-order processes, usually in conjunction with chemical kinetics or radioactive decay, stop at single, discrete dose events. Although single-dose situations are important, multiple-dose events,…

  18. Kinetic studies of elementary chemical reactions

    SciTech Connect

    Durant, J.L. Jr.

    1993-12-01

    This program concerning kinetic studies of elementary chemical reactions is presently focussed on understanding reactions of NH{sub x} species. To reach this goal, the author is pursuing experimental studies of reaction rate coefficients and product branching fractions as well as using electronic structure calculations to calculate transition state properties and reaction rate calculations to relate these properties to predicted kinetic behavior. The synergy existing between the experimental and theoretical studies allow one to gain a deeper insight into more complex elementary reactions.

  19. A high temperature and atmospheric pressure experimental and detailed chemical kinetic modelling study of 2-methyl furan oxidation

    PubMed Central

    Somers, Kieran P.; Simmie, John M.; Gillespie, Fiona; Burke, Ultan; Connolly, Jessica; Metcalfe, Wayne K.; Battin-Leclerc, Frédérique; Dirrenberger, Patricia; Herbinet, Olivier; Glaude, Pierre-Alexandre; Curran, Henry J.

    2013-01-01

    An experimental ignition delay time study for the promising biofuel 2-methyl furan (2MF) was performed at equivalence ratios of 0.5, 1.0 and 2.0 for mixtures of 1% fuel in argon in the temperature range 1200–1800 K at atmospheric pressure. Laminar burning velocities were determined using the heat-flux method for mixtures of 2MF in air at equivalence ratios of 0.55–1.65, initial temperatures of 298–398 K and atmospheric pressure. A detailed chemical kinetic mechanism consisting of 2059 reactions and 391 species has been constructed to describe the oxidation of 2MF and is used to simulate experiment. Accurate reproduction of the experimental data has been obtained over all conditions with the developed mechanism. Rate of production and sensitivity analyses have been carried out to identify important consumption pathways of the fuel and key kinetic parameters under these conditions. The reactions of hydrogen atom with the fuel are highlighted as important under all experimental conditions studied, with abstraction by the hydrogen atom promoting reactivity and hydrogen atom addition to the furan ring inhibiting reactivity. This work, to the authors knowledge, is the first to combine theoretical and experimental work to describe the oxidation of any of the alkylated furans. The mechanism developed herein to describe 2MF combustion should also function as a sub-mechanism to describe the oxidation of 2,5-dimethyl furan whilst also providing key insights into the oxidation of this similar biofuel candidate. PMID:23814505

  20. Mechanism of chemical degradation and determination of solubility by kinetic modeling of the highly unstable sesquiterpene lactone nobilin in different media.

    PubMed

    Thormann, Ursula; De Mieri, Maria; Neuburger, Markus; Verjee, Sheela; Altmann, Peter; Hamburger, Matthias; Imanidis, Georgios

    2014-10-01

    The objective of this work was first to investigate the chemical degradation of the sesquiterpene lactone nobilin and determine its solubility under conditions of concurrent degradation for partially amorphous starting material; second, to determine the effect of biorelevant media used in the in vitro measurement of intestinal absorption on degradation and solubility of nobilin. Purely aqueous medium (aq-TMCaco ), fasted and fed state simulated intestinal fluid (FaSSIF-TMCaco and FeSSIF-TMCaco ), and two liposomal formulations (LiposomesFaSSIF and LiposomesFeSSIF ) with the same lipid concentration as FaSSIF-TMCaco and FeSSIF-TMCaco were used. Degradation products were identified by nuclear magnetic resonance and X-ray crystallography and the order of reaction kinetics was determined. Solubility was deduced with a mathematical model encompassing dissolution, degradation, and reprecipitation kinetics that took into account particle size distribution of the solid material. Degradation mechanism of nobilin involved water-catalyzed opening of the lactone ring and transannular cyclization resulting in five degradation products. Degradation followed first-order kinetics in aq-TMCaco and FaSSIF-TMCaco , and higher-order kinetics in FeSSIF-TMCaco and the two liposomal formulations, whereas degradation in the latter media was diminished. Solubility of nobilin increased in the order: aq-TMCaco < FaSSIF-TMCaco , < LiposomesFaSSIF < FeSSIF-TMCaco < LiposomesFeSSIF . Improvement of stability and solubility was consistent with the incorporation of the nobilin molecule into colloidal lipid particles. The developed kinetic model is proposed to be a useful tool for deducing solubility under dynamic conditions. PMID:25176605

  1. Chemical Kinetics of Hydrocarbon Ignition in Practical Combustion Systems

    SciTech Connect

    Westbrook, C.K.

    2000-07-07

    Chemical kinetic factors of hydrocarbon oxidation are examined in a variety of ignition problems. Ignition is related to the presence of a dominant chain branching reaction mechanism that can drive a chemical system to completion in a very short period of time. Ignition in laboratory environments is studied for problems including shock tubes and rapid compression machines. Modeling of the laboratory systems are used to develop kinetic models that can be used to analyze ignition in practical systems. Two major chain branching regimes are identified, one consisting of high temperature ignition with a chain branching reaction mechanism based on the reaction between atomic hydrogen with molecular oxygen, and the second based on an intermediate temperature thermal decomposition of hydrogen peroxide. Kinetic models are then used to describe ignition in practical combustion environments, including detonations and pulse combustors for high temperature ignition, and engine knock and diesel ignition for intermediate temperature ignition. The final example of ignition in a practical environment is homogeneous charge, compression ignition (HCCI) which is shown to be a problem dominated by the kinetics intermediate temperature hydrocarbon ignition. Model results show why high hydrocarbon and CO emissions are inevitable in HCCI combustion. The conclusion of this study is that the kinetics of hydrocarbon ignition are actually quite simple, since only one or two elementary reactions are dominant. However, there are many combustion factors that can influence these two major reactions, and these are the features that vary from one practical system to another.

  2. Modeling the Human Kinetic Adjustment Factor for Inhaled Volatile Organic Chemicals: Whole Population Approach versus Distinct Subpopulation Approach

    PubMed Central

    Valcke, M.; Nong, A.; Krishnan, K.

    2012-01-01

    The objective of this study was to evaluate the impact of whole- and sub-population-related variabilities on the determination of the human kinetic adjustment factor (HKAF) used in risk assessment of inhaled volatile organic chemicals (VOCs). Monte Carlo simulations were applied to a steady-state algorithm to generate population distributions for blood concentrations (CAss) and rates of metabolism (RAMs) for inhalation exposures to benzene (BZ) and 1,4-dioxane (1,4-D). The simulated population consisted of various proportions of adults, elderly, children, neonates and pregnant women as per the Canadian demography. Subgroup-specific input parameters were obtained from the literature and P3M software. Under the “whole population” approach, the HKAF was computed as the ratio of the entire population's upper percentile value (99th, 95th) of dose metrics to the median value in either the entire population or the adult population. Under the “distinct subpopulation” approach, the upper percentile values in each subpopulation were considered, and the greatest resulting HKAF was retained. CAss-based HKAFs that considered the Canadian demography varied between 1.2 (BZ) and 2.8 (1,4-D). The “distinct subpopulation” CAss-based HKAF varied between 1.6 (BZ) and 8.5 (1,4-D). RAM-based HKAFs always remained below 1.6. Overall, this study evaluated for the first time the impact of underlying assumptions with respect to the interindividual variability considered (whole population or each subpopulation taken separately) when determining the HKAF. PMID:22523487

  3. Simulating flame lift-off characteristics of diesel and biodiesel fuels using detailed chemical-kinetic mechanisms and LES turbulence model.

    SciTech Connect

    Som, S; Longman, D. E.; Luo, Z; Plomer, M; Lu, T; Senecal, P.K.; Pomraning, E

    2012-01-01

    Combustion in direct-injection diesel engines occurs in a lifted, turbulent diffusion flame mode. Numerous studies indicate that the combustion and emissions in such engines are strongly influenced by the lifted flame characteristics, which are in turn determined by fuel and air mixing in the upstream region of the lifted flame, and consequently by the liquid breakup and spray development processes. From a numerical standpoint, these spray combustion processes depend heavily on the choice of underlying spray, combustion, and turbulence models. The present numerical study investigates the influence of different chemical kinetic mechanisms for diesel and biodiesel fuels, as well as Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) turbulence models on predicting flame lift-off lengths (LOLs) and ignition delays. Specifically, two chemical kinetic mechanisms for n-heptane (NHPT) and three for biodiesel surrogates are investigated. In addition, the RNG k-{epsilon} (RANS) model is compared to the Smagorinsky based LES turbulence model. Using adaptive grid resolution, minimum grid sizes of 250 {micro}m and 125 {micro}m were obtained for the RANS and LES cases respectively. Validations of these models were performed against experimental data from Sandia National Laboratories in a constant volume combustion chamber. Ignition delay and flame lift-off validations were performed at different ambient temperature conditions. The LES model predicts lower ignition delays and qualitatively better flame structures compared to the RNG k-{epsilon} model. The use of realistic chemistry and a ternary surrogate mixture, which consists of methyl decanoate, methyl 9-decenoate, and NHPT, results in better predicted LOLs and ignition delays. For diesel fuel though, only marginal improvements are observed by using larger size mechanisms. However, these improved predictions come at a significant increase in computational cost.

  4. Modelling Heart Rate Kinetics

    PubMed Central

    Zakynthinaki, Maria S.

    2015-01-01

    The objective of the present study was to formulate a simple and at the same time effective mathematical model of heart rate kinetics in response to movement (exercise). Based on an existing model, a system of two coupled differential equations which give the rate of change of heart rate and the rate of change of exercise intensity is used. The modifications introduced to the existing model are justified and discussed in detail, while models of blood lactate accumulation in respect to time and exercise intensity are also presented. The main modification is that the proposed model has now only one parameter which reflects the overall cardiovascular condition of the individual. The time elapsed after the beginning of the exercise, the intensity of the exercise, as well as blood lactate are also taken into account. Application of the model provides information regarding the individual’s cardiovascular condition and is able to detect possible changes in it, across the data recording periods. To demonstrate examples of successful numerical fit of the model, constant intensity experimental heart rate data sets of two individuals have been selected and numerical optimization was implemented. In addition, numerical simulations provided predictions for various exercise intensities and various cardiovascular condition levels. The proposed model can serve as a powerful tool for a complete means of heart rate analysis, not only in exercise physiology (for efficiently designing training sessions for healthy subjects) but also in the areas of cardiovascular health and rehabilitation (including application in population groups for which direct heart rate recordings at intense exercises are not possible or not allowed, such as elderly or pregnant women). PMID:25876164

  5. Modelling heart rate kinetics.

    PubMed

    Zakynthinaki, Maria S

    2015-01-01

    The objective of the present study was to formulate a simple and at the same time effective mathematical model of heart rate kinetics in response to movement (exercise). Based on an existing model, a system of two coupled differential equations which give the rate of change of heart rate and the rate of change of exercise intensity is used. The modifications introduced to the existing model are justified and discussed in detail, while models of blood lactate accumulation in respect to time and exercise intensity are also presented. The main modification is that the proposed model has now only one parameter which reflects the overall cardiovascular condition of the individual. The time elapsed after the beginning of the exercise, the intensity of the exercise, as well as blood lactate are also taken into account. Application of the model provides information regarding the individual's cardiovascular condition and is able to detect possible changes in it, across the data recording periods. To demonstrate examples of successful numerical fit of the model, constant intensity experimental heart rate data sets of two individuals have been selected and numerical optimization was implemented. In addition, numerical simulations provided predictions for various exercise intensities and various cardiovascular condition levels. The proposed model can serve as a powerful tool for a complete means of heart rate analysis, not only in exercise physiology (for efficiently designing training sessions for healthy subjects) but also in the areas of cardiovascular health and rehabilitation (including application in population groups for which direct heart rate recordings at intense exercises are not possible or not allowed, such as elderly or pregnant women). PMID:25876164

  6. Perspective: Stochastic algorithms for chemical kinetics

    PubMed Central

    Gillespie, Daniel T.; Hellander, Andreas; Petzold, Linda R.

    2013-01-01

    We outline our perspective on stochastic chemical kinetics, paying particular attention to numerical simulation algorithms. We first focus on dilute, well-mixed systems, whose description using ordinary differential equations has served as the basis for traditional chemical kinetics for the past 150 years. For such systems, we review the physical and mathematical rationale for a discrete-stochastic approach, and for the approximations that need to be made in order to regain the traditional continuous-deterministic description. We next take note of some of the more promising strategies for dealing stochastically with stiff systems, rare events, and sensitivity analysis. Finally, we review some recent efforts to adapt and extend the discrete-stochastic approach to systems that are not well-mixed. In that currently developing area, we focus mainly on the strategy of subdividing the system into well-mixed subvolumes, and then simulating diffusional transfers of reactant molecules between adjacent subvolumes together with chemical reactions inside the subvolumes. PMID:23656106

  7. Mechanistic quantitative structure-activity relationship model for the photoinduced toxicity of polycyclic aromatic hydrocarbons. 1: Physical model based on chemical kinetics in a two-compartment system

    SciTech Connect

    Krylov, S.N.; Huang, X.D.; Zeiler, L.F.; Dixon, D.G.; Greenberg, B.M.

    1997-11-01

    A quantitative structure-activity relationship model for the photoinduced toxicity of 16 polycyclic aromatic hydrocarbons (PAHs) to duckweed (Lemna gibba) in simulated solar radiation (SSR) was developed. Lemna gibba was chosen for this study because toxicity could be considered in two compartments: water column and leaf tissue. Modeling of photoinduced toxicity was described by photochemical reactions between PAHs and a hypothetical group of endogenous biomolecules (G) required for normal growth, with damage to G by PAHs and/or photomodified PAHs in SSR resulting in impaired growth. The reaction scheme includes photomodification of PAHs, uptake of PAHs into leaves, triplet-state formation of intact PAHs, photosensitization reactions that damage G, and reactions between photomodified PAHs and G. The assumptions used were: the PAH photomodification rate is slower than uptake of chemicals into leaves, the PAH concentration in aqueous solution is nearly constant during a toxicity test, the fluence rate of actinic radiation is lower within leaves than in the aqueous phase, and the toxicity of intact PAHs in the dark is negligible. A series of differential equations describing the reaction kinetics of intact and photomodifed PAHs with G was derived. The resulting equation for PAH toxicity was a function of treatment period, initial PAH concentration, relative absorbance of SSR by each PAH, quantum yield for formation of triplet-state PAH, and rate of PAH photomodification. Data for growth in the presence of intact and photomodified PAHs were used to empirically solve for a photosensitization constant (PSC) and a photomodification constant (PMC) for each of the 16 PAHs tested. For 9 PAHs the PMC dominates and for 7 PAHs the PSC dominates.

  8. High Fidelity Modeling of Turbulent Mixing and Chemical Kinetics Interactions in a Post-Detonation Flow Field

    NASA Astrophysics Data System (ADS)

    Sinha, Neeraj; Zambon, Andrea; Ott, James; Demagistris, Michael

    2015-06-01

    Driven by the continuing rapid advances in high-performance computing, multi-dimensional high-fidelity modeling is an increasingly reliable predictive tool capable of providing valuable physical insight into complex post-detonation reacting flow fields. Utilizing a series of test cases featuring blast waves interacting with combustible dispersed clouds in a small-scale test setup under well-controlled conditions, the predictive capabilities of a state-of-the-art code are demonstrated and validated. Leveraging physics-based, first principle models and solving large system of equations on highly-resolved grids, the combined effects of finite-rate/multi-phase chemical processes (including thermal ignition), turbulent mixing and shock interactions are captured across the spectrum of relevant time-scales and length scales. Since many scales of motion are generated in a post-detonation environment, even if the initial ambient conditions are quiescent, turbulent mixing plays a major role in the fireball afterburning as well as in dispersion, mixing, ignition and burn-out of combustible clouds in its vicinity. Validating these capabilities at the small scale is critical to establish a reliable predictive tool applicable to more complex and large-scale geometries of practical interest.

  9. Three-stage autoignition of gasoline in an HCCI engine: An experimental and chemical kinetic modeling investigation

    SciTech Connect

    Machrafi, Hatim; Cavadias, Simeon

    2008-12-15

    The alternative HCCI combustion mode presents a possible means for decreasing the pollution with respect to conventional gasoline or diesel engines, while maintaining the efficiency of a diesel engine or even increasing it. This paper investigates the possibility of using gasoline in an HCCI engine and analyzes the autoignition of gasoline in such an engine. The compression ratio that has been used is 13.5, keeping the inlet temperature at 70 C, varying the equivalence ratio from 0.3 to 0.54, and the EGR (represented by N{sub 2}) ratio from 0 to 37 vol%. For comparison, a PRF95 and a surrogate containing 11 vol% n-heptane, 59 vol% iso-octane, and 30 vol% toluene are used. A previously validated kinetic surrogate mechanism is used to analyze the experiments and to yield possible explanations to kinetic phenomena. From this work, it seems quite possible to use the high octane-rated gasoline for autoignition purposes, even under lean inlet conditions. Furthermore, it appeared that gasoline and its surrogate, unlike PRF95, show a three-stage autoignition. Since the PRF95 does not contain toluene, it is suggested by the kinetic mechanism that the benzyl radical, issued from toluene, causes this so-defined ''obstructed preignition'' and delaying thereby the final ignition for gasoline and its surrogate. The results of the kinetic mechanism supporting this explanation are shown in this paper. (author)

  10. Stochastic kinetic mean field model

    NASA Astrophysics Data System (ADS)

    Erdélyi, Zoltán; Pasichnyy, Mykola; Bezpalchuk, Volodymyr; Tomán, János J.; Gajdics, Bence; Gusak, Andriy M.

    2016-07-01

    This paper introduces a new model for calculating the change in time of three-dimensional atomic configurations. The model is based on the kinetic mean field (KMF) approach, however we have transformed that model into a stochastic approach by introducing dynamic Langevin noise. The result is a stochastic kinetic mean field model (SKMF) which produces results similar to the lattice kinetic Monte Carlo (KMC). SKMF is, however, far more cost-effective and easier to implement the algorithm (open source program code is provided on

  11. Evaluating temperature and fuel stratification for heat-release rate control in a reactivity-controlled compression-ignition engine using optical diagnostics and chemical kinetics modeling

    SciTech Connect

    Musculus, Mark P. B.; Kokjohn, Sage L.; Reitz, Rolf D.

    2015-04-23

    We investigated the combustion process in a dual-fuel, reactivity-controlled compression-ignition (RCCI) engine using a combination of optical diagnostics and chemical kinetics modeling to explain the role of equivalence ratio, temperature, and fuel reactivity stratification for heat-release rate control. An optically accessible engine is operated in the RCCI combustion mode using gasoline primary reference fuels (PRF). A well-mixed charge of iso-octane (PRF = 100) is created by injecting fuel into the engine cylinder during the intake stroke using a gasoline-type direct injector. Later in the cycle, n-heptane (PRF = 0) is delivered through a centrally mounted diesel-type common-rail injector. This injection strategy generates stratification in equivalence ratio, fuel blend, and temperature. The first part of this study uses a high-speed camera to image the injection events and record high-temperature combustion chemiluminescence. Moreover, the chemiluminescence imaging showed that, at the operating condition studied in the present work, mixtures in the squish region ignite first, and the reaction zone proceeds inward toward the center of the combustion chamber. The second part of this study investigates the charge preparation of the RCCI strategy using planar laser-induced fluorescence (PLIF) of a fuel tracer under non-reacting conditions to quantify fuel concentration distributions prior to ignition. The fuel-tracer PLIF data show that the combustion event proceeds down gradients in the n-heptane distribution. The third part of the study uses chemical kinetics modeling over a range of mixtures spanning the distributions observed from the fuel-tracer fluorescence imaging to isolate the roles of temperature, equivalence ratio, and PRF number stratification. The simulations predict that PRF number stratification is the dominant factor controlling the ignition location and growth rate of the reaction zone. Equivalence ratio has a smaller, but still significant

  12. Optimization of KINETICS Chemical Computation Code

    NASA Technical Reports Server (NTRS)

    Donastorg, Cristina

    2012-01-01

    NASA JPL has been creating a code in FORTRAN called KINETICS to model the chemistry of planetary atmospheres. Recently there has been an effort to introduce Message Passing Interface (MPI) into the code so as to cut down the run time of the program. There has been some implementation of MPI into KINETICS; however, the code could still be more efficient than it currently is. One way to increase efficiency is to send only certain variables to all the processes when an MPI subroutine is called and to gather only certain variables when the subroutine is finished. Therefore, all the variables that are used in three of the main subroutines needed to be investigated. Because of the sheer amount of code that there is to comb through this task was given as a ten-week project. I have been able to create flowcharts outlining the subroutines, common blocks, and functions used within the three main subroutines. From these flowcharts I created tables outlining the variables used in each block and important information about each. All this information will be used to determine how to run MPI in KINETICS in the most efficient way possible.

  13. Chemical Weathering Kinetics of Basalt on Venus

    NASA Technical Reports Server (NTRS)

    Fegley, Bruce, Jr.

    1997-01-01

    The purpose of this project was to experimentally measure the kinetics for chemical weathering reactions involving basalt on Venus. The thermochemical reactions being studied are important for the CO2 atmosphere-lithosphere cycle on Venus and for the atmosphere-surface reactions controlling the oxidation state of the surface of Venus. These reactions include the formation of carbonate and scapolite minerals, and the oxidation of Fe-bearing minerals. These experiments and calculations are important for interpreting results from the Pioneer Venus, Magellan, Galileo flyby, Venera, and Vega missions to Venus, for interpreting results from Earth-based telescopic observations, and for the design of new Discovery class (e.g., VESAT) and New Millennium missions to Venus such as geochemical landers making in situ elemental and mineralogical analyses, and orbiters, probes and balloons making spectroscopic observations of the sub-cloud atmosphere of Venus.

  14. Chemical Conversion Pathways and Kinetic Modeling for the OH-Initiated Reaction of Triclosan in Gas-Phase

    PubMed Central

    Zhang, Xue; Zhang, Chenxi; Sun, Xiaomin; Kang, Lingyan; Zhao, Yan

    2015-01-01

    As a widely used antimicrobial additive in daily consumption, attention has been paid to the degradation and conversion of triclosan for a long time. The quantum chemistry calculation and the canonical variational transition state theory are employed to investigate the mechanism and kinetic property. Besides addition and abstraction, oxidation pathways and further conversion pathways are also considered. The OH radicals could degrade triclosan to phenols, aldehydes, and other easily degradable substances. The conversion mechanisms of triclosan to the polychlorinated dibenzopdioxin and furan (PCDD/Fs) and polychlorinated biphenyls (PCBs) are clearly illustrated and the toxicity would be strengthened in such pathways. Single radical and diradical pathways are compared to study the conversion mechanism of dichlorodibenzo dioxin (DCDD). Furthermore, thermochemistry is discussed in detail. Kinetic property is calculated and the consequent ratio of kadd/ktotal and kabs/ktotal at 298.15 K are 0.955 and 0.045, respectively. Thus, the OH radical addition reactions are predominant, the substitute position of OH radical on triclosan is very important to generate PCDD and furan, and biradical is also a vital intermediate to produce dioxin. PMID:25867482

  15. Chemical conversion pathways and kinetic modeling for the OH-initiated reaction of triclosan in gas-phase.

    PubMed

    Zhang, Xue; Zhang, Chenxi; Sun, Xiaomin; Kang, Lingyan; Zhao, Yan

    2015-01-01

    As a widely used antimicrobial additive in daily consumption, attention has been paid to the degradation and conversion of triclosan for a long time. The quantum chemistry calculation and the canonical variational transition state theory are employed to investigate the mechanism and kinetic property. Besides addition and abstraction, oxidation pathways and further conversion pathways are also considered. The OH radicals could degrade triclosan to phenols, aldehydes, and other easily degradable substances. The conversion mechanisms of triclosan to the polychlorinated dibenzopdioxin and furan (PCDD/Fs) and polychlorinated biphenyls (PCBs) are clearly illustrated and the toxicity would be strengthened in such pathways. Single radical and diradical pathways are compared to study the conversion mechanism of dichlorodibenzo dioxin (DCDD). Furthermore, thermochemistry is discussed in detail. Kinetic property is calculated and the consequent ratio of k add/k total and k abs/k total at 298.15 K are 0.955 and 0.045, respectively. Thus, the OH radical addition reactions are predominant, the substitute position of OH radical on triclosan is very important to generate PCDD and furan, and biradical is also a vital intermediate to produce dioxin. PMID:25867482

  16. Detailed chemical kinetic oxidation mechanism for a biodiesel surrogate

    SciTech Connect

    Herbinet, O; Pitz, W J; Westbrook, C K

    2007-09-17

    A detailed chemical kinetic mechanism has been developed and used to study the oxidation of methyl decanoate, a surrogate for biodiesel fuels. This model has been built by following the rules established by Curran et al. for the oxidation of n-heptane and it includes all the reactions known to be pertinent to both low and high temperatures. Computed results have been compared with methyl decanoate experiments in an engine and oxidation of rapeseed oil methyl esters in a jet stirred reactor. An important feature of this mechanism is its ability to reproduce the early formation of carbon dioxide that is unique to biofuels and due to the presence of the ester group in the reactant. The model also predicts ignition delay times and OH profiles very close to observed values in shock tube experiments fueled by n-decane. These model capabilities indicate that large n-alkanes can be good surrogates for large methyl esters and biodiesel fuels to predict overall reactivity, but some kinetic details, including early CO2 production from biodiesel fuels, can be predicted only by a detailed kinetic mechanism for a true methyl ester fuel. The present methyl decanoate mechanism provides a realistic kinetic tool for simulation of biodiesel fuels.

  17. Detailed chemical kinetic oxidation mechanism for a biodiesel surrogate

    SciTech Connect

    Herbinet, O; Pitz, W J; Westbrook, C K

    2007-09-20

    A detailed chemical kinetic mechanism has been developed and used to study the oxidation of methyl decanoate, a surrogate for biodiesel fuels. This model has been built by following the rules established by Curran et al. for the oxidation of n-heptane and it includes all the reactions known to be pertinent to both low and high temperatures. Computed results have been compared with methyl decanoate experiments in an engine and oxidation of rapeseed oil methyl esters in a jet stirred reactor. An important feature of this mechanism is its ability to reproduce the early formation of carbon dioxide that is unique to biofuels and due to the presence of the ester group in the reactant. The model also predicts ignition delay times and OH profiles very close to observed values in shock tube experiments fueled by n-decane. These model capabilities indicate that large n-alkanes can be good surrogates for large methyl esters and biodiesel fuels to predict overall reactivity, but some kinetic details, including early CO{sub 2} production from biodiesel fuels, can be predicted only by a detailed kinetic mechanism for a true methyl ester fuel. The present methyl decanoate mechanism provides a realistic kinetic tool for simulation of biodiesel fuels.

  18. Degradation of di(2-ethyl hexyl) phthalate by Fusarium culmorum: Kinetics, enzymatic activities and biodegradation pathway based on quantum chemical modelingpathway based on quantum chemical modeling.

    PubMed

    Ahuactzin-Pérez, Miriam; Tlecuitl-Beristain, Saúl; García-Dávila, Jorge; González-Pérez, Manuel; Gutiérrez-Ruíz, María Concepción; Sánchez, Carmen

    2016-10-01

    Di(2-ethylhexyl) phthalate (DEHP) is a plasticizer widely used in the manufacture of plastics, and it is an environmental contaminant. The specific growth rate (μ), maximum biomass (Xmax), biodegradation constant of DEHP (k), half-life (t1/2) of DEHP biodegradation and removal efficiency of DEHP, esterase and laccase specific activities, and enzymatic yield parameters were evaluated for Fusarium culmorum grown on media containing glucose and different concentrations of DEHP (0, 500 and 1000mg/L). The greatest μ and the largest Xmax occurred in media supplemented with 1000mg of DEHP/L. F. culmorum degraded 95% of the highest amount of DEHP tested (1000mg/L) within 60h of growth. The k and t1/2 were 0.024h(-1) and 28h, respectively, for both DEHP concentrations. The removal efficiency of DEHP was 99.8% and 99.9% for 1000 and 500mg/L, respectively. Much higher specific esterase activity than specific laccase activity was observed in all media tested. The compounds of biodegradation of DEHP were identified by GC-MS. A DEHP biodegradation pathway by F. culmorum was proposed on the basis of the intermolecular flow of electrons of the identified intermediate compounds using quantum chemical modeling. DEHP was fully metabolized by F. culmorum with butanediol as the final product. This fungus offers great potential in bioremediation of environments polluted with DEHP. PMID:27277206

  19. Maximum Probability Reaction Sequences in Stochastic Chemical Kinetic Systems

    PubMed Central

    Salehi, Maryam; Perkins, Theodore J.

    2010-01-01

    The detailed behavior of many molecular processes in the cell, such as protein folding, protein complex assembly, and gene regulation, transcription and translation, can often be accurately captured by stochastic chemical kinetic models. We investigate a novel computational problem involving these models – that of finding the most-probable sequence of reactions that connects two or more states of the system observed at different times. We describe an efficient method for computing the probability of a given reaction sequence, but argue that computing most-probable reaction sequences is EXPSPACE-hard. We develop exact (exhaustive) and approximate algorithms for finding most-probable reaction sequences. We evaluate these methods on test problems relating to a recently-proposed stochastic model of folding of the Trp-cage peptide. Our results provide new computational tools for analyzing stochastic chemical models, and demonstrate their utility in illuminating the behavior of real-world systems. PMID:21629860

  20. Chemical Kinetic Reaction Mechanisms for Combustion of Hydrocarbon and Other Types of Chemical Fuels

    DOE Data Explorer

    The central feature of the Combustion Chemistry project at LLNL is the development, validation, and application of detailed chemical kinetic reaction mechanisms for the combustion of hydrocarbon and other types of chemical fuels. For the past 30 years, LLNL's Chemical Sciences Division has built hydrocarbon mechanisms for fuels from hydrogen and methane through much larger fuels including heptanes and octanes. Other classes of fuels for which models have been developed include flame suppressants such as halons and organophosphates, and air pollutants such as soot and oxides of nitrogen and sulfur. Reaction mechanisms have been tested and validated extensively through comparisons between computed results and measured data from laboratory experiments (e.g., shock tubes, laminar flames, rapid compression machines, flow reactors, stirred reactors) and from practical systems (e.g., diesel engines, spark-ignition engines, homogeneous charge, compression ignition (HCCI) engines). These kinetic models are used to examine a wide range of combustion systems.

  1. Role of substrate inhibition kinetics in enzymatic chemical oscillations.

    PubMed Central

    Shen, P; Larter, R

    1994-01-01

    Two chemical kinetic models are investigated using standard nonlinear dynamics techniques to determine the conditions under which substrate inhibition kinetics can lead to oscillations. The first model is a classical substrate inhibition scheme based on Michaelis-Menten kinetics and involves a single substrate. Only when this reaction takes place in a flow reactor (i.e., both substrate and product are taken to follow reversible flow terms) are oscillations observed; however, the range of parameter values over which such oscillations occur is so narrow it is experimentally unobservable. A second model based on a general mechanism applied to the kinetics of many pH-dependent enzymes is also studied. This second model includes both substrate inhibition kinetics as well as autocatalysis through the activation of the enzyme by hydrogen ion. We find that it is the autocatalysis that is always responsible for oscillatory behavior in this scheme. The substrate inhibition terms affect the steady-state behavior but do not lead to oscillations unless product inhibition or multiple substrates are present; this is a general conclusion we can draw from our studies of both the classical substrate inhibition scheme and the pH-dependent enzyme mechanism. Finally, an analysis of the nullclines for these two models allows us to prove that the nullcline slopes must have a negative value for oscillatory behavior to exist; this proof can explain our results. From our analysis, we conclude with a brief discussion of other enzymes that might be expected to produce oscillatory behavior based on a pH-dependent substrate inhibition mechanism. Images FIGURE 8 FIGURE 10 PMID:7819481

  2. Spectral Quasi-Equilibrium Manifold for Chemical Kinetics.

    PubMed

    Kooshkbaghi, Mahdi; Frouzakis, Christos E; Boulouchos, Konstantinos; Karlin, Iliya V

    2016-05-26

    The Spectral Quasi-Equilibrium Manifold (SQEM) method is a model reduction technique for chemical kinetics based on entropy maximization under constraints built by the slowest eigenvectors at equilibrium. The method is revisited here and discussed and validated through the Michaelis-Menten kinetic scheme, and the quality of the reduction is related to the temporal evolution and the gap between eigenvalues. SQEM is then applied to detailed reaction mechanisms for the homogeneous combustion of hydrogen, syngas, and methane mixtures with air in adiabatic constant pressure reactors. The system states computed using SQEM are compared with those obtained by direct integration of the detailed mechanism, and good agreement between the reduced and the detailed descriptions is demonstrated. The SQEM reduced model of hydrogen/air combustion is also compared with another similar technique, the Rate-Controlled Constrained-Equilibrium (RCCE). For the same number of representative variables, SQEM is found to provide a more accurate description. PMID:27116566

  3. Reduced chemical kinetic mechanisms for hydrocarbon fuels

    SciTech Connect

    Montgomery, C J; Cremer, M A; Heap, M P; Chen, J -Y; Westbrook, C K; Maurice, L Q

    1999-12-10

    Using CARM (Computer Aided Reduction Method), a computer program that automates the mechanism reduction process, a variety of different reduced chemical kinetic mechanisms for ethylene and n-heptane have been generated. The reduced mechanisms have been compared to detailed chemistry calculations in simple homogeneous reactors and experiments. Reduced mechanisms for combustion of ethylene having as few as 10 species were found to give reasonable agreement with detailed chemistry over a range of stoichiometries and showed significant improvement over currently used global mechanisms. The performance of reduced mechanisms derived from a large detailed mechanism for n-heptane was compared to results from a reduced mechanism derived from a smaller semi-empirical mechanism. The semi-empirical mechanism was advantageous as a starting point for reduction for ignition delay, but not for PSR calculations. Reduced mechanisms with as few as 12 species gave excellent results for n-heptane/air PSR calculations but 16-25 or more species are needed to simulate n-heptane ignition delay.

  4. Kinetics investigations of atmospheric chemical reactions

    SciTech Connect

    Hills, A.J.

    1987-01-01

    Two separate gas-phase kinetics investigations were performed using a low-pressure fast-flow system with mass spectrometer detection. The first part of this research was a study of the atmospheric reactivity of diatomic sulfur, S/sub 2/. Rates of the reactions of sulfur with O, O/sub 2/, O/sub 3/, N/sub 2/O, NO, and NO/sub 2/ were investigated at 409 K and low pressure (0.89-3.0 Torr) in a discharge-flow system with mass spectrometric detection. The second investigation involves a study of the synergistic coupling of atmospheric bromine and chlorine chemistry. Recent measurements of ozone in the stratosphere over Antarctica have shown that the springtime ozone column decreased by 40% from 1960 to 1985. Both dynamical and chemical theories have been advanced to explain the formation of the Antarctic ozone hole. Prominent among these theories is that a synergistic interaction between gas-phase BrO and ClO radicals may be responsible for springtime ozone loss. The overall rate constant for the reaction, BrO + ClO ..-->.. Br + OClO ..-->.. Br + Cl + O/sub 2/ ..-->.. BrCl + O/sub 2/, has been measured over the temperature range 241-408 K. The rate constant for the overall reaction equals (8.2 +/- 1.0) 10/sup -12/ cm/sup 3//molecule s, independent of temperature.

  5. Detailed Chemical Kinetic Mechanisms for Combustion of Oxygenated Fuels

    SciTech Connect

    Fisher, E.M.; Pitz, W.J.; Curran, H.J.; Westbrook, C.K.

    2000-01-11

    Thermodynamic properties and detailed chemical kinetic models have been developed for the combustion of two oxygenates: methyl butanoate, a model compound for biodiesel fuels, and methyl formate, a related simpler molecule. Bond additivity methods and rules for estimating kinetic parameters were adopted from hydrocarbon combustion and extended. The resulting mechanisms have been tested against the limited combustion data available in the literature, which was obtained at low temperature, subatmospheric conditions in closed vessels, using pressure measurements as the main diagnostic. Some qualitative agreement was obtained, but the experimental data consistently indicated lower overall reactivities than the model, differing by factors of 10 to 50. This discrepancy, which occurs for species with well-established kinetic mechanisms as well as for methyl esters, is tentatively ascribed to the presence of wall reactions in the experiments. The model predicts a region of weak or negative dependence of overall reaction rate on temperature for each methyl ester. Examination of the reaction fluxes provides an explanation of this behavior, involving a temperature-dependent competition between chain-propagating unimolecular decomposition processes and chain-branching processes, similar to that accepted for hydrocarbons. There is an urgent need to obtain more complete experimental data under well-characterized conditions for thorough testing of the model.

  6. Reaction Mechanism Generator: Automatic construction of chemical kinetic mechanisms

    NASA Astrophysics Data System (ADS)

    Gao, Connie W.; Allen, Joshua W.; Green, William H.; West, Richard H.

    2016-06-01

    Reaction Mechanism Generator (RMG) constructs kinetic models composed of elementary chemical reaction steps using a general understanding of how molecules react. Species thermochemistry is estimated through Benson group additivity and reaction rate coefficients are estimated using a database of known rate rules and reaction templates. At its core, RMG relies on two fundamental data structures: graphs and trees. Graphs are used to represent chemical structures, and trees are used to represent thermodynamic and kinetic data. Models are generated using a rate-based algorithm which excludes species from the model based on reaction fluxes. RMG can generate reaction mechanisms for species involving carbon, hydrogen, oxygen, sulfur, and nitrogen. It also has capabilities for estimating transport and solvation properties, and it automatically computes pressure-dependent rate coefficients and identifies chemically-activated reaction paths. RMG is an object-oriented program written in Python, which provides a stable, robust programming architecture for developing an extensible and modular code base with a large suite of unit tests. Computationally intensive functions are cythonized for speed improvements.

  7. The know unknowns: Detailed simulations and low-order modeling to characterize facility-induced non-idealities in chemical-kinetics experiments

    NASA Astrophysics Data System (ADS)

    Ihme, Matthias

    2013-11-01

    Experimental investigations to study chemical-kinetics processes, reaction-rates or ignition properties are frequently accompanied by facility-induced non-idealities. Examples are turbulence and thermo-viscous boundary layers in rapid compression machines, temperature fluctuations and mixture inhomogeneities in flow-reactors, or shock-bifurcations and pressure drifts in shock-tubes. Although experimental investigations are carefully conducted to mitigate these effects, they are difficult to quantify experimentally. Simulations can assist in identifying these non-idealities and in guiding experimental instrumentation to improve measurement accuracies. This presentation discusses three different modeling approaches to characterize facility-effects in rapid compression machines, flow reactors, and shock-tubes. After providing an overview about these facilities and describing the underlying models, examples are presented to illustrate effects of turbulence, mixture-inhomogeneities, heat-losses, and thermal stratification on the ignition dynamics in these facilities. Diagnostics is developed to assess the sensitivity of the induction chemistry and to quantify reliable operating regimes that are not contaminated by these non-ideal processes.

  8. Evaluating temperature and fuel stratification for heat-release rate control in a reactivity-controlled compression-ignition engine using optical diagnostics and chemical kinetics modeling

    DOE PAGESBeta

    Musculus, Mark P. B.; Kokjohn, Sage L.; Reitz, Rolf D.

    2015-04-23

    We investigated the combustion process in a dual-fuel, reactivity-controlled compression-ignition (RCCI) engine using a combination of optical diagnostics and chemical kinetics modeling to explain the role of equivalence ratio, temperature, and fuel reactivity stratification for heat-release rate control. An optically accessible engine is operated in the RCCI combustion mode using gasoline primary reference fuels (PRF). A well-mixed charge of iso-octane (PRF = 100) is created by injecting fuel into the engine cylinder during the intake stroke using a gasoline-type direct injector. Later in the cycle, n-heptane (PRF = 0) is delivered through a centrally mounted diesel-type common-rail injector. This injectionmore » strategy generates stratification in equivalence ratio, fuel blend, and temperature. The first part of this study uses a high-speed camera to image the injection events and record high-temperature combustion chemiluminescence. Moreover, the chemiluminescence imaging showed that, at the operating condition studied in the present work, mixtures in the squish region ignite first, and the reaction zone proceeds inward toward the center of the combustion chamber. The second part of this study investigates the charge preparation of the RCCI strategy using planar laser-induced fluorescence (PLIF) of a fuel tracer under non-reacting conditions to quantify fuel concentration distributions prior to ignition. The fuel-tracer PLIF data show that the combustion event proceeds down gradients in the n-heptane distribution. The third part of the study uses chemical kinetics modeling over a range of mixtures spanning the distributions observed from the fuel-tracer fluorescence imaging to isolate the roles of temperature, equivalence ratio, and PRF number stratification. The simulations predict that PRF number stratification is the dominant factor controlling the ignition location and growth rate of the reaction zone. Equivalence ratio has a smaller, but still

  9. LSENS - GENERAL CHEMICAL KINETICS AND SENSITIVITY ANALYSIS CODE

    NASA Technical Reports Server (NTRS)

    Bittker, D. A.

    1994-01-01

    which provides the relationships between the predictions of a kinetics model and the input parameters of the problem. LSENS provides for efficient and accurate chemical kinetics computations and includes sensitivity analysis for a variety of problems, including nonisothermal conditions. LSENS replaces the previous NASA general chemical kinetics codes GCKP and GCKP84. LSENS is designed for flexibility, convenience and computational efficiency. A variety of chemical reaction models can be considered. The models include static system, steady one-dimensional inviscid flow, reaction behind an incident shock wave including boundary layer correction, and the perfectly stirred (highly backmixed) reactor. In addition, computations of equilibrium properties can be performed for the following assigned states, enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static problems LSENS computes sensitivity coefficients with respect to the initial values of the dependent variables and/or the three rates coefficient parameters of each chemical reaction. To integrate the ODEs describing chemical kinetics problems, LSENS uses the packaged code LSODE, the Livermore Solver for Ordinary Differential Equations, because it has been shown to be the most efficient and accurate code for solving such problems. The sensitivity analysis computations use the decoupled direct method, as implemented by Dunker and modified by Radhakrishnan. This method has shown greater efficiency and stability with equal or better accuracy than other methods of sensitivity analysis. LSENS is written in FORTRAN 77 with the exception of the NAMELIST extensions used for input. While this makes the code fairly machine independent, execution times on IBM PC compatibles would be unacceptable to most users. LSENS has been successfully implemented on a Sun4 running SunOS and a DEC VAX running VMS. With minor modifications, it should also be easily implemented on other

  10. Chemical kinetic and photochemical data for use in stratospheric modeling evaluation number 4: NASA panel for data evaluation

    NASA Technical Reports Server (NTRS)

    1981-01-01

    Evaluated sets of rate constants and photochemical cross sections compiled by the Panel are presented. The primary application of the data is in the modelling of stratospheric processes, with particular emphasis on the ozone layer and its possible perturbation by anthropogenic and natural phenomena.

  11. Theory of Chemical Modeling

    NASA Astrophysics Data System (ADS)

    Kühn, Michael

    In order to deal with the complexity of natural systems simplified models are employed to illustrate the principal and regulatory factors controlling a chemical system. Following the aphorism of Albert Einstein: Everything should be made as simple as possible, but not simpler, models need not to be completely realistic to be useful (Stumm and Morgan 1996), but need to meet a successful balance between realism and practicality. Properly constructed, a model is neither too simplified that it is unrealistic nor too detailed that it cannot be readily evaluated and applied to the problem of interest (Bethke 1996). The results of a model have to be at least partially observable or experimentally verifiable (Zhu and Anderson 2002). Geochemical modeling theories are presented here in a sequence of increasing complexity from geochemical equilibrium models to kinetic, reaction path, and finally coupled transport and reaction models. The description is far from complete but provides the needs for the set up of reactive transport models of hydrothermal systems as done within subsequent chapters. Extensive reviews of geochemical models in general can be found in the literature (Appelo and Postma 1999, Bethke 1996, Melchior and Bassett 1990, Nordstrom and Ball 1984, Paschke and van der Heijde 1996).

  12. Analytical Derivation of Moment Equations in Stochastic Chemical Kinetics

    PubMed Central

    Sotiropoulos, Vassilios; Kaznessis, Yiannis N.

    2011-01-01

    The master probability equation captures the dynamic behavior of a variety of stochastic phenomena that can be modeled as Markov processes. Analytical solutions to the master equation are hard to come by though because they require the enumeration of all possible states and the determination of the transition probabilities between any two states. These two tasks quickly become intractable for all but the simplest of systems. Instead of determining how the probability distribution changes in time, we can express the master probability distribution as a function of its moments, and, we can then write transient equations for the probability distribution moments. In 1949, Moyal defined the derivative, or jump, moments of the master probability distribution. These are measures of the rate of change in the probability distribution moment values, i.e. what the impact is of any given transition between states on the moment values. In this paper we present a general scheme for deriving analytical moment equations for any N-dimensional Markov process as a function of the jump moments. Importantly, we propose a scheme to derive analytical expressions for the jump moments for any N-dimensional Markov process. To better illustrate the concepts, we focus on stochastic chemical kinetics models for which we derive analytical relations for jump moments of arbitrary order. Chemical kinetics models are widely used to capture the dynamic behavior of biological systems. The elements in the jump moment expressions are a function of the stoichiometric matrix and the reaction propensities, i.e the probabilistic reaction rates. We use two toy examples, a linear and a non-linear set of reactions, to demonstrate the applicability and limitations of the scheme. Finally, we provide an estimate on the minimum number of moments necessary to obtain statistical significant data that would uniquely determine the dynamics of the underlying stochastic chemical kinetic system. The first two moments

  13. A Kinetic-fluid Model

    SciTech Connect

    First Author = C.Z. Cheng; Jay R. Johnson

    1998-07-10

    A nonlinear kinetic-fluid model for high-beta plasmas with multiple ion species which can be applied to multiscale phenomena is presented. The model embeds important kinetic effects due to finite ion Larmor radius (FLR), wave-particle resonances, magnetic particle trapping, etc. in the framework of simple fluid descriptions. When further restricting to low frequency phenomena with frequencies less than the ion cyclotron frequency the kinetic-fluid model takes a simpler form in which the fluid equations of multiple ion species collapse into single-fluid density and momentum equations and a low frequency generalized Ohm's law. The kinetic effects are introduced via plasma pressure tensors for ions and electrons which are computed from particle distribution functions that are governed by the Vlasov equation or simplified plasma dynamics equations such as the gyrokinetic equation. The ion FLR effects provide a finite parallel electric field, a perpendicular velocity that modifies the ExB drift, and a gyroviscosity tensor, all of which are neglected in the usual one-fluid MHD description. Eigenmode equations are derived which include magnetosphere-ionosphere coupling effects for low frequency waves (e.g., kinetic/inertial Alfven waves and ballooning-mirror instabilities).

  14. Kinetic Modeling of Divertor Plasma

    NASA Astrophysics Data System (ADS)

    Ishiguro, Seiji; Hasegawa, Hiroki; Pianpanit, Theerasarn

    2015-11-01

    Particle-in-Cell (PIC) simulation with the Monte Carlo collisions and the cumulative scattering angle coulomb collision can present kinetic dynamics of divertor plasmas. We are developing two types of PIC codes. The first one is the three dimensional bounded PIC code where three dimensional kinetic dynamics of blob is studied and current flow structures related to sheath formation are unveiled. The second one is the one spatial three velocity space dimensional (1D3V) PIC code with the Monte Carlo collisions where formation of detach plasma is studied. First target of our research is to construct self-consistent full kinetic simulation modeling of the linear divertor simulation experiments. This work is performed with the support and under the auspices of NIFS Collaboration Research program (NIFS15KNSS059, NIFS14KNXN279, and NIFS13KNSS038) and the Research Cooperation Program on Hierarchy and Holism in Natural Science at NINS.

  15. Comparing microbial and chemical kinetics for modelling soil organic carbon decomposition using the DecoChem v1.0 and DecoBio v1.0 models

    NASA Astrophysics Data System (ADS)

    Xenakis, G.; Williams, M.

    2014-07-01

    Soil organic matter is a vast store of carbon, with a critical role in the global carbon cycle. Despite its importance, the dynamics of soil organic carbon decomposition, under the impact of climate change or changing litter inputs, are poorly understood. Current biogeochemical models usually lack microbial processes and thus miss an important feedback when considering the fate of carbon. Here we use a series of modelling experiments to evaluate two different model structures: one with a standard first-order kinetic representation of soil decomposition (DecoChem v1.0, hereafter chemical model) and one with control of soil decomposition through microbial activity (DecoBio v1.0, hereafter biological model). The biological model includes cycling of organic matter into and out of microbial biomass, and simulates the decay rate as a functional of microbial activity. We tested two hypotheses. First, we hypothesized different responses in the two models to increased litter inputs and glucose additions. In the microbial model we hypothesized that this perturbation would prime microbial activity and reduce soil carbon stocks; in the chemical model we expected this perturbation to increase C stocks. In the biological model, responses to changed litter quantity were more rapid, but with the residence time of soil C altering such that soil C stocks were buffered. However, in the biological model there was a strong response to increased glucose additions (i.e. changes in litter quality), with significant losses to soil C stocks over time, driven by priming. Secondly, we hypothesized that warming will stimulate decomposition in the chemical model and loss of C, but in the biological model soil C will be less sensitive to warming, due to complex microbial feedbacks. The numerical experiments supported this hypothesis, with the chemical model soil C residence times and steady-state C stocks adjusting strongly with temperature changes, extending over decades. On the other hand, the

  16. Understanding Chemical Reaction Kinetics and Equilibrium with Interlocking Building Blocks

    ERIC Educational Resources Information Center

    Cloonan, Carrie A.; Nichol, Carolyn A.; Hutchinson, John S.

    2011-01-01

    Chemical reaction kinetics and equilibrium are essential core concepts of chemistry but are challenging topics for many students, both at the high school and undergraduate university level. Visualization at the molecular level is valuable to aid understanding of reaction kinetics and equilibrium. This activity provides a discovery-based method to…

  17. Symmetry Relations in Chemical Kinetics Arising from Microscopic Reversibility

    NASA Astrophysics Data System (ADS)

    Adib, Artur B.

    2006-01-01

    It is shown that the kinetics of time-reversible chemical reactions having the same equilibrium constant but different initial conditions are closely related to one another by a directly measurable symmetry relation analogous to chemical detailed balance. In contrast to detailed balance, however, this relation does not require knowledge of the elementary steps that underlie the reaction, and remains valid in regimes where the concept of rate constants is ill defined, such as at very short times and in the presence of low activation barriers. Numerical simulations of a model of isomerization in solution are provided to illustrate the symmetry under such conditions, and potential applications in protein folding or unfolding are pointed out.

  18. Kinetic Modeling of Microbiological Processes

    SciTech Connect

    Liu, Chongxuan; Fang, Yilin

    2012-08-26

    Kinetic description of microbiological processes is vital for the design and control of microbe-based biotechnologies such as waste water treatment, petroleum oil recovery, and contaminant attenuation and remediation. Various models have been proposed to describe microbiological processes. This editorial article discusses the advantages and limiation of these modeling approaches in cluding tranditional, Monod-type models and derivatives, and recently developed constraint-based approaches. The article also offers the future direction of modeling researches that best suit for petroleum and environmental biotechnologies.

  19. The Multiplexed Chemical Kinetic Photoionization Mass Spectrometer: A New Approach To Isomer-resolved Chemical Kinetics

    SciTech Connect

    Osborne, David L.; Zou, Peng; Johnsen, Howard; Hayden, Carl C.; Taatjes, Craig A.; Knyazev, Vadim D.; North, Simon W.; Peterka, Darcy S.; Ahmed, Musahid; Leone, Stephen R.

    2008-08-28

    We have developed a multiplexed time- and photon-energy?resolved photoionizationmass spectrometer for the study of the kinetics and isomeric product branching of gasphase, neutral chemical reactions. The instrument utilizes a side-sampled flow tubereactor, continuously tunable synchrotron radiation for photoionization, a multi-massdouble-focusing mass spectrometer with 100percent duty cycle, and a time- and positionsensitive detector for single ion counting. This approach enables multiplexed, universal detection of molecules with high sensitivity and selectivity. In addition to measurement of rate coefficients as a function of temperature and pressure, different structural isomers can be distinguished based on their photoionization efficiency curves, providing a more detailed probe of reaction mechanisms. The multiplexed 3-dimensional data structure (intensity as a function of molecular mass, reaction time, and photoionization energy) provides insights that might not be available in serial acquisition, as well as additional constraints on data interpretation.

  20. Kinetic model of HIV infection

    SciTech Connect

    Zhdanov, V. P.

    2007-10-15

    Recent experiments clarifying the details of exhaustion of CD8 T cells specific to various strains of human immunodeficiency virus (HIV) are indicative of slow irreversible (on a one-year time scale) deterioration of the immune system. The conventional models of HIV kinetics do not take this effect into account. Removing this shortcoming, we show the likely influence of such changes on the escape of HIV from control of the immune system.

  1. Hungarian University Students' Misunderstandings in Thermodynamics and Chemical Kinetics

    ERIC Educational Resources Information Center

    Turanyi, Tamas; Toth, Zoltan

    2013-01-01

    The misunderstandings related to thermodynamics (including chemical equilibrium) and chemical kinetics of first and second year Hungarian students of chemistry, environmental science, biology and pharmacy were investigated. We demonstrated that Hungarian university students have similar misunderstandings in physical chemistry to those reported in…

  2. Chemical kinetics computer program for static and flow reactions

    NASA Technical Reports Server (NTRS)

    Bittker, D. A.; Scullin, V. J.

    1972-01-01

    General chemical kinetics computer program for complex gas mixtures has been developed. Program can be used for any homogeneous reaction in either one dimensional flow or static system. It is flexible, accurate, and easy to use. It can be used for any chemical system for which species thermodynamic data and reaction rate constant data are known.

  3. Detailed Chemical Kinetic Reaction Mechanisms for Incineration of Organophosphorus and Fluoro-Organophosphorus Compounds

    SciTech Connect

    Glaude, P A; Melius, C; Pitz, W J; Westbrook, C K

    2001-12-13

    A detailed chemical kinetic reaction mechanism is developed to describe incineration of the chemical warfare nerve agent sarin (GB), based on commonly used principles of bond additivity and hierarchical reaction mechanisms. The mechanism is based on previous kinetic models of organophosphorus compounds such as TMP, DMMP and DIMP that are often used as surrogates to predict incineration of GB. Kinetic models of the three surrogates and GB are then used to predict their consumption in a perfectly stirred reactor fueled by natural gas to simulate incineration of these chemicals. Computed results indicate that DIMP is the only one of these surrogates that adequately describes combustion of GB under comparable conditions. The kinetic pathways responsible for these differences in reactivity are identified and discussed. The most important reaction in GB and DIMP that makes them more reactive than TMP or DMMP is found to be a six-center molecular elimination reaction producing propene.

  4. A kinetic-theory approach to turbulent chemically reacting flows

    NASA Technical Reports Server (NTRS)

    Chung, P. M.

    1976-01-01

    The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

  5. Kinetic models of conjugated metabolic cycles

    NASA Astrophysics Data System (ADS)

    Ershov, Yu. A.

    2016-01-01

    A general method is developed for the quantitative kinetic analysis of conjugated metabolic cycles in the human organism. This method is used as a basis for constructing a kinetic graph and model of the conjugated citric acid and ureapoiesis cycles. The results from a kinetic analysis of the model for these cycles are given.

  6. Modelling reaction kinetics inside cells

    PubMed Central

    Grima, Ramon; Schnell, Santiago

    2009-01-01

    In the past decade, advances in molecular biology such as the development of non-invasive single molecule imaging techniques have given us a window into the intricate biochemical activities that occur inside cells. In this article we review four distinct theoretical and simulation frameworks: (1) non-spatial and deterministic, (2) spatial and deterministic, (3) non-spatial and stochastic and (4) spatial and stochastic. Each framework can be suited to modelling and interpreting intracellular reaction kinetics. By estimating the fundamental length scales, one can roughly determine which models are best suited for the particular reaction pathway under study. We discuss differences in prediction between the four modelling methodologies. In particular we show that taking into account noise and space does not simply add quantitative predictive accuracy but may also lead to qualitatively different physiological predictions, unaccounted for by classical deterministic models. PMID:18793122

  7. Ernest Rutherford, Avogadro's Number, and Chemical Kinetics

    NASA Astrophysics Data System (ADS)

    Leenson, I. A.

    1998-08-01

    The paper presents a way for students to use data from Rutherford's works (1908 - 1911) in order to determine one of the most precise values of Avogadro Constant available at the beginning of the century. A brief discussion of earlier and modern methods for the determination of this fundamental constant is followed by vast quotations from the works of Rutherford, Boltwood and Geiger. Then there are given a dozen of problems and questions for students about these classical experiments; they vary in complexity from rather simple to quite challenging. Additional information and hints are provided to help the students in solving the problems. The last part contains detailed answers and solutions to all problems. The article will be useful for students of general chemistry, radiochemistry and physical chemistry (kinetics).

  8. HCCI in a CFR engine: experiments and detailed kinetic modeling

    SciTech Connect

    Flowers, D; Aceves, S; Smith, R; Torres, J; Girard, J; Dibble, R

    1999-11-05

    Single cylinder engine experiments and chemical kinetic modeling have been performed to study the effect of variations in fuel, equivalence ratio, and intake charge temperature on the start of combustion and the heat release rate. Neat propane and a fuel blend of 15% dimethyl-ether in methane have been studied. The results demonstrate the role of these parameters on the start of combustion, efficiency, imep, and emissions. Single zone kinetic modeling results show the trends consistent with the experimental results.

  9. A coke oven model including thermal decomposition kinetics of tar

    SciTech Connect

    Munekane, Fuminori; Yamaguchi, Yukio; Tanioka, Seiichi

    1997-12-31

    A new one-dimensional coke oven model has been developed for simulating the amount and the characteristics of by-products such as tar and gas as well as coke. This model consists of both heat transfer and chemical kinetics including thermal decomposition of coal and tar. The chemical kinetics constants are obtained by estimation based on the results of experiments conducted to investigate the thermal decomposition of both coal and tar. The calculation results using the new model are in good agreement with experimental ones.

  10. Spectral method for a kinetic swarming model

    NASA Astrophysics Data System (ADS)

    Gamba, Irene M.; Haack, Jeffrey R.; Motsch, Sebastien

    2015-09-01

    In this paper 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. We observe that the kinetic model captures key features such as vortex formation and traveling waves.

  11. Spectral method for a kinetic swarming model

    SciTech Connect

    Gamba, Irene M.; Haack, Jeffrey R.; Motsch, Sebastien

    2015-04-28

    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.

  12. Prospective Chemistry Teachers' Conceptions of Chemical Thermodynamics and Kinetics

    ERIC Educational Resources Information Center

    Sozbilir, Mustafa; Pinarbasi, Tacettin; Canpolat, Nurtac

    2010-01-01

    This study aimed at identifying specifically prospective chemistry teachers' difficulties in determining the differences between the concepts of chemical thermodynamics and kinetics. Data were collected from 67 prospective chemistry teachers at Kazim Karabekir Education Faculty of Ataturk University in Turkey during 2005-2006 academic year. Data…

  13. Chemical Kinetics at the Single-Molecule Level

    ERIC Educational Resources Information Center

    Levitus, Marcia

    2011-01-01

    For over a century, chemists have investigated the rates of chemical reactions using experimental conditions involving huge numbers of molecules. As a consequence, the description of the kinetics of the reaction in terms of average values was good enough for all practical purposes. From the pedagogical point of view, such a description misses the…

  14. Chemical TOPAZ: Modifications to the heat transfer code TOPAZ: The addition of chemical reaction kinetics and chemical mixtures

    SciTech Connect

    Nichols, A.L. III.

    1990-06-07

    This is a report describing the modifications which have been made to the heat flow code TOPAZ to allow the inclusion of thermally controlled chemical kinetics. This report is broken into parts. The first part is an introduction to the general assumptions and theoretical underpinning that were used to develop the model. The second section describes the changes that have been implemented into the code. The third section is the users manual for the input for the code. The fourth section is a compilation of hints, common errors, and things to be aware of while you are getting started. The fifth section gives a sample problem using the new code. This manual addenda is written with the presumption that most readers are not fluent with chemical concepts. Therefore, we shall in this section endeavor to describe the requirements that must be met before chemistry can occur and how we have modeled the chemistry in the code.

  15. Chemical kinetic performance losses for a hydrogen laser thermal thruster

    NASA Technical Reports Server (NTRS)

    Mccay, T. D.; Dexter, C. E.

    1985-01-01

    Projected requirements for efficient, economical, orbit-raising propulsion systems have generated investigations into several potentially high specific impulse, moderate thrust, advanced systems. One of these systems, laser thermal propulsion, utilizes a high temperature plasma as the enthalpy source. The plasma is sustained by a focused laser beam which maintains the plasma temperature at levels near 20,000 K. Since such temperature levels lead to total dissociation and high ionization, the plasma thruster system potentially has a high specific impulse decrement due to recombination losses. The nozzle flow is expected to be sufficiently nonequilibrium to warrant concern over the achievable specific impluse. This investigation was an attempt at evaluation of those losses. The One-Dimensional Kinetics (ODK) option of the Two-Dimensional Kinetics (TDK) Computer Program was used with a chemical kinetics rate set obtained from available literature to determine the chemical kinetic energy losses for typical plasma thruster conditions. The rates were varied about the nominal accepted values to band the possible losses. Kinetic losses were shown to be highly significant for a laser thermal thruster using hydrogen. A 30 percent reduction in specific impulse is possible simply due to the inability to completely extract the molecular recombination energy.

  16. Chemical kinetic reaction mechanism for the combustion of propane

    NASA Technical Reports Server (NTRS)

    Jachimowski, C. J.

    1984-01-01

    A detailed chemical kinetic reaction mechanism for the combustion of propane is presented and discussed. The mechanism consists of 27 chemical species and 83 elementary chemical reactions. Ignition and combustion data as determined in shock tube studies were used to evaluate the mechanism. Numerical simulation of the shock tube experiments showed that the kinetic behavior predicted by the mechanism for stoichiometric mixtures is in good agrement with the experimental results over the entire temperature range examined (1150-2600K). Sensitivity and theoretical studies carried out using the mechanism revealed that hydrocarbon reactions which are involved in the formation of the HO2 radical and the H2O2 molecule are very important in the mechanism and that the observed nonlinear behavior of ignition delay time with decreasing temperature can be interpreted in terms of the increased importance of the HO2 and H2O2 reactions at the lower temperatures.

  17. Infrared absorption spectroscopy and chemical kinetics of free radicals

    SciTech Connect

    Curl, R.F.; Glass, G.P.

    1993-12-01

    This research is directed at the detection, monitoring, and study of chemical kinetic behavior by infrared absorption spectroscopy of small free radical species thought to be important intermediates in combustion. During the last year, infrared kinetic spectroscopy using excimer laser flash photolysis and color-center laser probing has been employed to study the high resolution spectrum of HCCN, the rate constant of the reaction between ethynyl (C{sub 2}H) radical and H{sub 2} in the temperature region between 295 and 875 K, and the recombination rate of propargyl (CH{sub 2}CCH) at room temperature.

  18. Fourth-Order Vibrational Transition State Theory and Chemical Kinetics

    NASA Astrophysics Data System (ADS)

    Stanton, John F.; Matthews, Devin A.; Gong, Justin Z.

    2015-06-01

    Second-order vibrational perturbation theory (VPT2) is an enormously successful and well-established theory for treating anharmonic effects on the vibrational levels of semi-rigid molecules. Partially as a consequence of the fact that the theory is exact for the Morse potential (which provides an appropriate qualitative model for stretching anharmonicity), VPT2 calculations for such systems with appropriate ab initio potential functions tend to give fundamental and overtone levels that fall within a handful of wavenumbers of experimentally measured positions. As a consequence, the next non-vanishing level of perturbation theory -- VPT4 -- offers only slight improvements over VPT2 and is not practical for most calculations since it requires information about force constants up through sextic. However, VPT4 (as well as VPT2) can be used for other applications such as the next vibrational correction to rotational constants (the ``gammas'') and other spectroscopic parameters. In addition, the marriage of VPT with the semi-classical transition state theory of Miller (SCTST) has recently proven to be a powerful and accurate treatment for chemical kinetics. In this talk, VPT4-based SCTST tunneling probabilities and cumulative reaction probabilities are give for the first time for selected low-dimensional model systems. The prospects for VPT4, both practical and intrinsic, will also be discussed.

  19. A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reaction

    NASA Technical Reports Server (NTRS)

    Lian, Yongsheng; Xu, Kun

    1999-01-01

    This paper concerns the extension of the multicomponent gas-kinetic BGK-type scheme to multidimensional chemical reactive flow calculations. In the kinetic model, each component satisfies its individual gas-kinetic BGK equation and the equilibrium states of both components are coupled in space and time due to the momentum and energy exchange in the course of particle collisions. At the same time, according to the chemical reaction rule one component can be changed into another component with the release of energy, where the reactant and product could have different gamma. Many numerical test cases are included in this paper, which show the robustness and accuracy of kinetic approach in the description of multicomponent reactive flows.

  20. Kinetic modelling of mitochondrial translation.

    PubMed

    Korla, Kalyani; Mitra, Chanchal K

    2014-01-01

    Mitochondrial genome contains 13 protein coding genes, all being part of the oxidative phosphorylation complexes. The process of translation of these protein coding mRNAs in mitochondrial matrix is a good miniature model of translation in cytoplasm. In this work, we have simulated three phases of mitochondrial translation viz. initiation, elongation and termination (including ribosome recycling). The kinetic equations for these phases have been deduced based on the information available in literature. Various factors involved in the process have been included explicitly. Kinetic simulation was done using Octave, open source software. Scripts were written individually for each phase. Initiation begins with mitoribosome, mRNA, fMet-tRNA and initiation factors. The final product of the initiation script, the initiation complex, was introduced as the start point in the successive step, i.e. elongation. Elongation is a particular extensive process where the various aminoacyl-tRNAs already present in the matrix check for matching with the triplet codon in A-site of mitoribosome. This script consists of two parts: one with the time behaviour of the factors involved in the molecular process (using ordinary differential equation solver) and the other including the reading of triplet codon on the mRNA and incorporating the corresponding aminoacyl-tRNA, and then at each step elongating the peptide chain (using loops and conditions). The peptide chain thus formed in the elongation step (in the loops and conditions segment) was released in the termination step. This was followed by mitoribosome recycling where the mitoribosome reached the native state and was ready for the next cycle of translation. PMID:24028553

  1. Chemical kinetics of cetane number improving agents

    SciTech Connect

    Hashimoto, K.; Akutsu, Y.; Arai, M.; Tamura, M.

    1996-12-31

    The increasing demand for diesel fuels has resulted in the use of greater percentage of cracked distillates having poor ignition properties. The ignition properties of diesel fuels can be rated in terms of their cetane number and diesel fuels having low cetane number may have poor ignition properties such as diesel knock, difficulties to start engines in the cold weather and so on. Such diesel fuels need cetane number improving agents. In the 1940s and 1950s alkyl nitrates, alkyl nitrites and organic peroxides were found to be effective cetane number improving additives. Our recent study suggests that free radicals produced from thermal decomposition just before ignition should have an important role to improve their ignition properties. However no studies on the reaction mechanism for improving effect of these additives have been attempted because of complex nature of spontaneous ignition reaction of hydrocarbons. In order to clarify the reaction mechanism for improving effects of cetane number improving agents. We here have attempted to simulate the spontaneous ignition of n-butane as a model compound in the presence of alkyl nitrites as cetane number improving agents.

  2. Chemical Dynamics, Molecular Energetics, and Kinetics at the Synchrotron

    SciTech Connect

    Leone, Stephen R.; Ahmed, Musahid; Wilson, Kevin R.

    2010-03-14

    Scientists at the Chemical Dynamics Beamline of the Advanced Light Source in Berkeley are continuously reinventing synchrotron investigations of physical chemistry and chemical physics with vacuum ultraviolet light. One of the unique aspects of a synchrotron for chemical physics research is the widely tunable vacuum ultraviolet light that permits threshold ionization of large molecules with minimal fragmentation. This provides novel opportunities to assess molecular energetics and reaction mechanisms, even beyond simple gas phase molecules. In this perspective, significant new directions utilizing the capabilities at the Chemical Dynamics Beamline are presented, along with an outlook for future synchrotron and free electron laser science in chemical dynamics. Among the established and emerging fields of investigations are cluster and biological molecule spectroscopy and structure, combustion flame chemistry mechanisms, radical kinetics and product isomer dynamics, aerosol heterogeneous chemistry, planetary and interstellar chemistry, and secondary neutral ion-beam desorption imaging of biological matter and materials chemistry.

  3. Biomass torrefaction: modeling of volatile and solid product evolution kinetics.

    PubMed

    Bates, Richard B; Ghoniem, Ahmed F

    2012-11-01

    The aim of this work is the development of a kinetics model for the evolution of the volatile and solid product composition during torrefaction conditions between 200 and 300°C. Coupled to an existing two step solid mass loss kinetics mechanism, this model describes the volatile release kinetics in terms of a set of identifiable chemical components, permitting the solid product composition to be estimated by mass conservation. Results show that most of the volatiles released during the first stage include highly oxygenated species such as water, acetic acid, and carbon dioxide, while volatiles released during the second step are composed primarily of lactic acid, methanol, and acetic acid. This kinetics model will be used in the development of a model to describe reaction energy balance and heat release dynamics. PMID:23026268

  4. CERENA: ChEmical REaction Network Analyzer—A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics

    PubMed Central

    Kazeroonian, Atefeh; Fröhlich, Fabian; Raue, Andreas; Theis, Fabian J.; Hasenauer, Jan

    2016-01-01

    Gene expression, signal transduction and many other cellular processes are subject to stochastic fluctuations. The analysis of these stochastic chemical kinetics is important for understanding cell-to-cell variability and its functional implications, but it is also challenging. A multitude of exact and approximate descriptions of stochastic chemical kinetics have been developed, however, tools to automatically generate the descriptions and compare their accuracy and computational efficiency are missing. In this manuscript we introduced CERENA, a toolbox for the analysis of stochastic chemical kinetics using Approximations of the Chemical Master Equation solution statistics. CERENA implements stochastic simulation algorithms and the finite state projection for microscopic descriptions of processes, the system size expansion and moment equations for meso- and macroscopic descriptions, as well as the novel conditional moment equations for a hybrid description. This unique collection of descriptions in a single toolbox facilitates the selection of appropriate modeling approaches. Unlike other software packages, the implementation of CERENA is completely general and allows, e.g., for time-dependent propensities and non-mass action kinetics. By providing SBML import, symbolic model generation and simulation using MEX-files, CERENA is user-friendly and computationally efficient. The availability of forward and adjoint sensitivity analyses allows for further studies such as parameter estimation and uncertainty analysis. The MATLAB code implementing CERENA is freely available from http://cerenadevelopers.github.io/CERENA/. PMID:26807911

  5. CERENA: ChEmical REaction Network Analyzer--A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics.

    PubMed

    Kazeroonian, Atefeh; Fröhlich, Fabian; Raue, Andreas; Theis, Fabian J; Hasenauer, Jan

    2016-01-01

    Gene expression, signal transduction and many other cellular processes are subject to stochastic fluctuations. The analysis of these stochastic chemical kinetics is important for understanding cell-to-cell variability and its functional implications, but it is also challenging. A multitude of exact and approximate descriptions of stochastic chemical kinetics have been developed, however, tools to automatically generate the descriptions and compare their accuracy and computational efficiency are missing. In this manuscript we introduced CERENA, a toolbox for the analysis of stochastic chemical kinetics using Approximations of the Chemical Master Equation solution statistics. CERENA implements stochastic simulation algorithms and the finite state projection for microscopic descriptions of processes, the system size expansion and moment equations for meso- and macroscopic descriptions, as well as the novel conditional moment equations for a hybrid description. This unique collection of descriptions in a single toolbox facilitates the selection of appropriate modeling approaches. Unlike other software packages, the implementation of CERENA is completely general and allows, e.g., for time-dependent propensities and non-mass action kinetics. By providing SBML import, symbolic model generation and simulation using MEX-files, CERENA is user-friendly and computationally efficient. The availability of forward and adjoint sensitivity analyses allows for further studies such as parameter estimation and uncertainty analysis. The MATLAB code implementing CERENA is freely available from http://cerenadevelopers.github.io/CERENA/. PMID:26807911

  6. Accelerating the Computation of Detailed Chemical Reaction Kinetics for Simulating Combustion of Complex Fuels

    SciTech Connect

    Grout, Ray W

    2012-01-01

    Combustion of hydrocarbon fuels has been a very challenging scientific and engineering problem due to the complexity of turbulent flows and hydrocarbon reaction kinetics. There is an urgent need to develop an efficient modeling capability to accurately predict the combustion of complex fuels. Detailed chemical kinetic models for the surrogates of fuels such as gasoline, diesel and JP-8 consist of thousands of chemical species and Arrhenius reaction steps. Oxygenated fuels such as bio-fuels and heavier hydrocarbons, such as from newer fossil fuel sources, are expected to have a much more complex chemistry requiring increasingly larger chemical kinetic models. Such models are beyond current computational capability, except for homogeneous or partially stirred reactor type calculations. The advent of highly parallel multi-core processors and graphical processing units (GPUs) promises a steep increase in computational performance in the coming years. This paper will present a software framework that translates the detailed chemical kinetic models to high- performance code targeted for GPU accelerators.

  7. Accelerating the Computation of Detailed Chemical Reaction Kinetics for Simulating Combustion of Complex Fuels

    SciTech Connect

    Sankaran, R.; Grout, R.

    2012-01-01

    Combustion of hydrocarbon fuels has been a very challenging scientific and engineering problem due to the complexity of turbulent flows and hydrocarbon reaction kinetics. There is an urgent need to develop an efficient modeling capability to accurately predict the combustion of complex fuels. Detailed chemical kinetic models for the surrogates of fuels such as gasoline, diesel and JP-8 consist of thousands of chemical species and Arrhenius reaction steps. Oxygenated fuels such as bio-fuels and heavier hydrocarbons, such as from newer fossil fuel sources, are expected to have a much more complex chemistry requiring increasingly larger chemical kinetic models. Such models are beyond current computational capability, except for homogeneous or partially stirred reactor type calculations. The advent of highly parallel multi-core processors and graphical processing units (GPUs) promises a steep increase in computational performance in the coming years. This paper will present a software framework that translates the detailed chemical kinetic models to high-performance code targeted for GPU accelerators.

  8. Fast algorithm for calculating chemical kinetics in turbulent reacting flow

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.; Pratt, D. T.

    1986-01-01

    This paper addresses the need for a fast batch chemistry solver to perform the kinetics part of a split operator formulation of turbulent reacting flows, with special attention focused on the solution of the ordinary differential equations governing a homogeneous gas-phase chemical reaction. For this purpose, a two-part predictor-corrector algorithm which incorporates an exponentially fitted trapezoidal method was developed. The algorithm performs filtering of ill-posed initial conditions, automatic step-size selection, and automatic selection of Jacobi-Newton or Newton-Raphson iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm, termed CREK1D (combustion reaction kinetics, one-dimensional), compared favorably with the code LSODE when tested on two representative problems drawn from combustion kinetics, and is faster than LSODE.

  9. Chemical kinetic mechanism for the oxidation of paraffinic hydrocarbons needed for primary reference fuels

    SciTech Connect

    Westbrook, C.K.; Pitz, W.J.

    1993-03-01

    A detailed chemical kinetic reaction mechanism is described which simulates the oxidation of the primary reference fuels n-heptane and iso-octane. The high temperature subset of these mechanisms is identified, and the extensions to deal with low temperature conditions are also explained. The algorithms used to assign reaction rates to elementary steps in the reaction mechanism are described, and the means of identifying the different chemical species and the relevant reactions are outlined. Finally, we show how interested kinetic modeling researchers can obtain copies of this reaction mechanism.

  10. [Study on Chemical Kinetic Effect of Dielectric Barrier Discharge Plasma].

    PubMed

    Zrang, Peng; Hong, Yan-ji; Shen, Shuang-yan; Ding, Xiao-yu; Ma, Di

    2015-03-01

    To reveal the mechanism of plasma (assisted the ignition process of methane/air further, schematic of dielectric barrier discharge plasma system with atmospheric air was designed and set up, the emission spectrum of dielectric barrier discharge plasma with atmospheric air was measured, and the active particles produced by the interaction of dielectric barrier discharge plasma with atmospheric air were analyzed with the spectrum technology, the ignition model and calculation methods of sensitivity analysis and reaction path analysis were given, effects of NO and O3 on the ignition delay time were simulated, and the chemical kinetics mechanism of NO and O3 assisted ignition was revealed via sensitivity analysis and reaction path analysis. The results show that main excited particles of N2 and O3 are generated via effect of plasma on the atmospheric air, which are converted into active particles of NO(ξ) and O3 in the end, the life of which are longer than any other active particles, effects of plasma on the ignition is simplified as effects of NO(ξ) and O3 on the ignition; NO and O3 could reduce the ignition delay time significantly, but the amplitude decrease with increase of the initial temperature, this is because the rate of ignition is decided by the oxidation rate of CH3, the oxidized pathway of CH3 is R155 and R156 for auto-ignition and their rates are slower when temperature is low, so the ignition delay time of methane/air is longer; NO could reduce the ignition delay time significantly because of the oxidized pathway of CH3 is changed to R327 CH3O2 + NO = CH3O + NO2, R328 CH3 + NO2 = CH3O + NO for NO(ξ) (assisted ignition process from R155 and R156 for auto-ignition; and the chemical kinetic effect is the dominating factor of O3 on the ignition and which change the reaction path. PMID:26117883

  11. A multipurpose reduced chemical-kinetic mechanism for methanol combustion

    NASA Astrophysics Data System (ADS)

    Fernández-Tarrazo, Eduardo; Sánchez-Sanz, Mario; Sánchez, Antonio L.; Williams, Forman A.

    2016-07-01

    A multipurpose reduced chemical-kinetic mechanism for methanol combustion comprising 8 overall reactions and 11 reacting chemical species is presented. The development starts by investigating the minimum set of elementary reactions needed to describe methanol combustion with reasonable accuracy over a range of conditions of temperature, pressure, and composition of interest in combustion. Starting from a 27-step mechanism that has been previously tested and found to give accurate predictions of ignition processes for these conditions, it is determined that the addition of 11 elementary reactions taken from its basis (San Diego) mechanism extends the validity of the description to premixed-flame propagation, strain-induced extinction of non-premixed flames, and equilibrium composition and temperatures, giving results that compare favourably with experimental measurements and also with computations using the 247-step detailed San Diego mechanism involving 50 reactive species. Specifically, premixed-flame propagation velocities and extinction strain rates for non-premixed counterflow flames calculated with the 38-step mechanism show departures from experimental measurements and detailed-chemistry computations that are roughly on the order of 10%, comparable with expected experimental uncertainties. Similar accuracy is found in comparisons of autoignition times over the range considered, except at very high temperatures, under which conditions the computations tend to overpredict induction times for all of the chemistry descriptions tested. From this 38-step mechanism, the simplification is continued by introducing steady-state approximations for the intermediate species CH3, CH4, HCO, CH3O, CH2OH, and O, leading to an 8-step reduced mechanism that provides satisfactory accuracy for all conditions tested. The flame computations indicate that thermal diffusion has a negligible influence on methanol combustion in all cases considered and that a mixture-average species

  12. Kinetic model for erythrocyte aggregation.

    PubMed

    Bertoluzzo, S M; Bollini, A; Rasia, M; Raynal, A

    1999-01-01

    It is well known that light transmission through blood is the most widely utilized method for the study of erythrocyte aggregation. The curves obtained had been considered empirically as exponential functions. In consequence, the process becomes characterized by an only parameter that varies with all the process factors without discrimination. In the present paper a mathematical model for RBC aggregation process is deduced in accordance with von Smoluchowski's theory about the kinetics of colloidal particles agglomeration. The equation fitted the experimental pattern of the RBC suspension optical transmittance closely and contained two parameters that estimate the most important characteristics of the aggregation process separately, i.e., (1) average size of rouleaux at equilibrium and (2) aggregation rate. The evaluation of the method was assessed by some factors affecting erythrocyte aggregation, such as temperature, plasma dilutions, Dextran 500, Dextran 70 and PVP 360, at different media concentrations, cellular membrane alteration by the alkylating agent TCEA, and decrease of medium osmolarity. Results were interpreted considering the process characteristics estimated by the parameters, and there were also compared with similar studies carried out by other authors with other methods. This analysis allowed us to conclude that the equation proposed is reliable and useful to study erythrocyte aggregation. PMID:10660481

  13. Prediction of Kinetic Isotope Effects for Various Hydride Transfer Reactions Using a New Kinetic Model.

    PubMed

    Shen, Guang-Bin; Xia, Ke; Li, Xiu-Tao; Li, Jun-Ling; Fu, Yan-Hua; Yuan, Lin; Zhu, Xiao-Qing

    2016-03-24

    In this work, kinetic isotope effect (KIEself) values of 68 hydride self-exchange reactions, XH(D) + X(+) → X(+) + XH(D), in acetonitrile at 298 K were determined using a new experimental method. KIE values of 4556 hydride cross transfer reactions, XH(D) + Y(+) → X(+) + YH(D), in acetonitrile were estimated from the 68 determined KIEself values of hydride self-exchange reactions using a new KIE relation formula derived from Zhu's kinetic equation and the reliability of the estimations was verified using different experimental methods. A new KIE kinetic model to explain and predict KIE values was developed according to Zhu's kinetic model using two different Morse free energy curves instead of one Morse free energy curve in the traditional KIE theories to describe the free energy changes of X-H bond and X-D bond dissociation in chemical reactions. The most significant contribution of this paper to KIE theory is to build a new KIE kinetic model, which can be used to not only uniformly explain the various (normal, enormous and inverse) KIE values but also safely prodict KIE values of various chemical reactions. PMID:26938149

  14. The Role of Comprehensive Detailed Chemical Kinetic Reaction Mechanisms in Combustion Research

    SciTech Connect

    Westbrook, C K; Pitz, W J; Curran, H J; Mehl, M

    2008-07-16

    Recent developments by the authors in the field of comprehensive detailed chemical kinetic reaction mechanisms for hydrocarbon fuels are reviewed. Examples are given of how these mechanisms provide fundamental chemical insights into a range of combustion applications. Practical combustion consists primarily of chemical heat release from reactions between a fuel and an oxidizer, and computer simulations of practical combustion systems have become an essential tool of combustion research (Westbrook et al., 2005). At the heart of most combustion simulations, the chemical kinetic submodel frequently is the most detailed, complex and computationally costly part of a system model. Historically, the chemical submodel equations are solved using time-implicit numerical algorithms, due to the extreme stiffness of the coupled rate equations, with a computational cost that varies roughly with the cube of the number of chemical species in the model. While early mechanisms (c. 1980) for apparently simple fuels such as methane (Warnatz, 1980) or methanol (Westbrook and Dryer, 1979) included perhaps 25 species, current detailed mechanisms for much larger, more complex fuels such as hexadecane (Fournet et al., 2001; Ristori et al., 2001; Westbrook et al., 2008) or methyl ester methyl decanoate (Herbinet et al., 2008) have as many as 2000 or even 3000 species. Rapid growth in capabilities of modern computers has been an essential feature in this rapid growth in the size and complexity of chemical kinetic reaction mechanisms.

  15. The combustion chemistry of a fuel tracer: Measured flame speeds and ignition delays and a detailed chemical kinetic model for the oxidation of acetone

    SciTech Connect

    Pichon, S.; Black, G.; Simmie, J.M.; Curran, H.J.; Chaumeix, N.; Yahyaoui, M.; Donohue, R.

    2009-02-15

    Acetone ignition delay and stretch-free laminar flame speed measurements have been carried out and a kinetic model has been developed to simulate these and literature data for acetone and for ketene, which was found to be an important intermediate in its oxidation. The mechanism has been based on one originally devised for dimethyl ether and modified through validation of the hydrogen, carbon monoxide and methane sub-mechanisms. Acetone oxidation in argon was studied behind reflected shock waves in the temperature range 1340-1930 K, at 1 atm and at equivalence ratios of 0.5, 1 and 2; it is also shown that the addition of up to 15% acetone to a stoichiometric n-heptane mixture has no effect on the measured ignition delay times. Flame speeds at 298 K and 1 atm of pure acetone in air were measured in a spherical bomb; a maximum flame speed of {proportional_to}35 cm s{sup -1} at {phi}=1.15 is indicated. (author)

  16. Seizure modeling of Pb(II) and Cd(II) from aqueous solution by chemically modified sugarcane bagasse fly ash: isotherms, kinetics, and column study.

    PubMed

    Shah, Bhavna; Mistry, Chirag; Shah, Ajay

    2013-04-01

    Heavy metal pollution is a common environmental problem all over the world. The purpose of the research is to examine the applicability of bagasse fly ash (BFA)-an agricultural waste of sugar industry used for the synthesis of zeolitic material. The zeolitic material are used for the uptake of Pb(II) and Cd(II) heavy metal. Bagasse fly ash is used as a native material for the synthesis of zeolitic materials by conventional hydrothermal treatment without (conventional zeolitic bagasse fly ash (CZBFA)) and with electrolyte (conventional zeolitic bagasse fly ash in electrolyte media (ECZBFA)) media. Heavy metal ions Pb(II) and Cd(II) were successfully seized from aqueous media using these synthesized zeolitic materials. In this study, the zeolitic materials were well characterized by different instrumental methods such as Brunauer-Emmett-Teller, XRF, Fourier transform infrared spectroscopy, powder X-ray diffraction, and scanning electron microscopic microphotographs. The presence of analcime, phillipsite, and zeolite P in adsorbents confirms successful conversion of native BFA into zeolitic materials. Seizure modeling of Pb(II) and Cd(II) was achieved by batch sorption experiments, isotherms, and kinetic studies. These data were used to compare and evaluate the zeolitic materials as potential sorbents for the uptake of heavy metal ions from an aqueous media. The Langmuir isotherm correlation coefficient parameters best fit the equilibrium data which indicate the physical sorption. Pseudo-second-order and intra-particle diffusion model matches best which indicates that the rate of sorption was controlled by film diffusion. The column studies were performed for the practical function of sorbents, and breakthrough curves were obtained, which revealed higher sorption capacity as compared to batch method. Synthesized zeolitic material (CZBFA and ECZBFA), a low-cost sorbent, was proven as potential sorbent for the uptake of Pb(II) and Cd(II) heavy metal ions. PMID:22739768

  17. Theoretical validation of chemical kinetic mechanisms : combustion of methanol.

    SciTech Connect

    Skodje, R. T.; Tomlin, A. S.; Klippenstein, S. J.; Harding, L. B.; Davis, M. J.; Chemical Sciences and Engineering Division; Univ. of Colorado; Univ. of Leeds

    2010-08-19

    A new technique is proposed that uses theoretical methods to systematically improve the performance of chemical kinetic mechanisms. Using a screening method, the chemical reaction steps that most strongly influence a given kinetic observable are identified. The associated rate coefficients are then improved by high-level quantum chemistry and transition-state-theory calculations, which leads to new values for the coefficients and smaller uncertainty ranges. This updating process is continued as new reactions emerge as the most important steps in the target observable. The screening process employed is a global sensitivity analysis that involves Monte Carlo sampling of the full N-dimensional uncertainty space of rate coefficients, where N is the number of reaction steps. The method is applied to the methanol combustion mechanism of Li et al. (Int. J. Chem. Kinet. 2007, 39, 109.). It was found that the CH{sub 3}OH + HO{sub 2} and CH{sub 3}OH + O{sub 2} reactions were the most important steps in setting the ignition delay time, and the rate coefficients for these reactions were updated. The ignition time is significantly changed for a broad range of high-concentration methanol/oxygen mixtures in the updated mechanism.

  18. MODELING CHEMICAL INTERACTIONS IN ANAEROBIC BIOFILM SYSTEMS

    EPA Science Inventory

    Rigorous steady-state models of acetate- and methanol-utilizing methanogenic biofilms are developed taking into account the mass transfer of neutral and ionic species, pH changes within the biofilm, pH-dependent Monod kinetics, chemical equilibrium,, electronueutrality, gas produ...

  19. Site-specific functionalization for chemical speciation of Cr(III) and Cr(VI) using polyaniline impregnated nanocellulose composite: equilibrium, kinetic, and thermodynamic modeling

    NASA Astrophysics Data System (ADS)

    Jain, Priyanka; Varshney, Shilpa; Srivastava, Shalini

    2015-10-01

    Site-specific functionalizations are the emergent attention for the enhancement of sorption latent of heavy metals. Limited chemistry has been applied for the fabrication of diafunctionalized materials having potential to tether both environmentally stable oxidation states of chromium (Cr(III) and Cr(VI). Polyaniline impregnated nanocellulose composite (PANI-NCC) has been fabricated using click chemistry and explored for the removal of Cr(III) and Cr(VI) from hydrological environment. The structure, stability, morphology, particle size, surface area, hydrophilicity, and porosity of fabricated PANI-NCC were characterized comprehensively using analytical techniques and mathematical tools. The maximum sorption performance of PANI-NCC was procured for (Cr(III): 47.06 mg g-1; 94.12 %) and (Cr(VI): 48.92 mg g-1; 97.84 %) by equilibrating 0.5 g sorbent dose with 1000 mL of 25 mg L-1 chromium conc. at pH 6.5 and 2.5 for Cr(III) and Cr(VI), respectively. The sorption data showed a best fit to the Langmuir isotherm and pseudo-second-order kinetic model. The negative value of ∆ G° (-8.59 and -11.16 kJ mol-1) and ∆ H° (66.46 × 10-1 and 17.84 × 10-1 kJ mol-1), and positive value of ∆ S° (26.66 and 31.46 J mol-1K-1) for Cr(III) and Cr(VI), respectively, reflect the spontaneous, feasibility, and exothermic nature of the sorption process. The application of fabricated PANI-NCC for removing both the forms of chromium in the presence of other heavy metals was also tested at laboratory and industrial waste water regime. These findings open up new avenues in the row of high performance, scalable, and economic nanobiomaterial for the remediation of both forms of chromium from water streams.

  20. LSENS, a general chemical kinetics and sensitivity analysis code for gas-phase reactions: User's guide

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, Krishnan; Bittker, David A.

    1993-01-01

    A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS, are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include static system, steady, one-dimensional, inviscid flow, shock initiated reaction, and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method, which works efficiently for the extremes of very fast and very slow reaction, is used for solving the 'stiff' differential equation systems that arise in chemical kinetics. For static reactions, sensitivity coefficients of all dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters can be computed. This paper presents descriptions of the code and its usage, and includes several illustrative example problems.

  1. In vitro-in vivo extrapolation of quantitative hepatic biotransformation data for fish - I. A review of methods, and strategies for incorporating intrinsic clearance estimates into chemical kinetic models

    SciTech Connect

    Nichols, John W.; Schultz, Irv R.; Fitzsimmons, Patrick N..

    2006-06-10

    Mammalian researchers have developed a stepwise approach to predict in vivo hepatic clearance from measurements of in vitro hepatic metabolism. The resulting clearance estimates have been used to screen drug candidates, identify potential drug-drug interactions, investigate idiosyncratic drug responses, and support toxicology risk assessments. In this report we review these methods, discuss their potential application to studies with fish, and describe how extrapolated values could be incorporated into well-known compartmental kinetic models. Empirical equations that relate extrapolation factors to chemical log Kow are given to facilitate the incorporation of metabolism data into bioconcentration and bioaccumulation models. Because they explicitly incorporate the concept of clearance, compartmental clearance volume models are particularly well suited for incorporating hepatic clearance estimates. The manner in which these clearance values are incorporated into a given model depends, however, on the measurement frame of reference. Procedures for the incorporation of in vitro metabolism data into physiologically based toxicokinetic (PBTK) models are also described. Unlike most compartmental models, PBTK models are developed to describe the effects of metabolism in the tissue where it occurs. In addition, PBTK models are well suited to modeling metabolism in more than one tissue.

  2. Summary of the LLNL one-dimensional transport-kinetics model of the troposphere and stratosphere: 1981

    SciTech Connect

    Wuebbles, D.J.

    1981-09-01

    Since the LLNL one-dimensional coupled transport and chemical kinetics model of the troposphere and stratosphere was originally developed in 1972 (Chang et al., 1974), there have been many changes to the model's representation of atmospheric physical and chemical processes. A brief description is given of the current LLNL one-dimensional coupled transport and chemical kinetics model of the troposphere and stratosphere.

  3. Viral kinetic modeling: state of the art

    SciTech Connect

    Canini, Laetitia; Perelson, Alan S.

    2014-06-25

    Viral kinetic modeling has led to increased understanding of the within host dynamics of viral infections and the effects of therapy. Here we review recent developments in the modeling of viral infection kinetics with emphasis on two infectious diseases: hepatitis C and influenza. We review how viral kinetic modeling has evolved from simple models of viral infections treated with a drug or drug cocktail with an assumed constant effectiveness to models that incorporate drug pharmacokinetics and pharmacodynamics, as well as phenomenological models that simply assume drugs have time varying-effectiveness. We also discuss multiscale models that include intracellular events in viral replication, models of drug-resistance, models that include innate and adaptive immune responses and models that incorporate cell-to-cell spread of infection. Overall, viral kinetic modeling has provided new insights into the understanding of the disease progression and the modes of action of several drugs. In conclusion, we expect that viral kinetic modeling will be increasingly used in the coming years to optimize drug regimens in order to improve therapeutic outcomes and treatment tolerability for infectious diseases.

  4. Viral kinetic modeling: state of the art

    DOE PAGESBeta

    Canini, Laetitia; Perelson, Alan S.

    2014-06-25

    Viral kinetic modeling has led to increased understanding of the within host dynamics of viral infections and the effects of therapy. Here we review recent developments in the modeling of viral infection kinetics with emphasis on two infectious diseases: hepatitis C and influenza. We review how viral kinetic modeling has evolved from simple models of viral infections treated with a drug or drug cocktail with an assumed constant effectiveness to models that incorporate drug pharmacokinetics and pharmacodynamics, as well as phenomenological models that simply assume drugs have time varying-effectiveness. We also discuss multiscale models that include intracellular events in viralmore » replication, models of drug-resistance, models that include innate and adaptive immune responses and models that incorporate cell-to-cell spread of infection. Overall, viral kinetic modeling has provided new insights into the understanding of the disease progression and the modes of action of several drugs. In conclusion, we expect that viral kinetic modeling will be increasingly used in the coming years to optimize drug regimens in order to improve therapeutic outcomes and treatment tolerability for infectious diseases.« less

  5. An Experimental and Chemical Kinetics Study of the Combustion of Syngas and High Hydrogen Content Fuels

    SciTech Connect

    Santoro, Robers; Dryer, Frederick; Ju, Yiguang

    2013-09-30

    An integrated and collaborative effort involving experiments and complementary chemical kinetic modeling investigated the effects of significant concentrations of water and CO2 and minor contaminant species (methane [CH4], ethane [C2H6], NOX, etc.) on the ignition and combustion of HHC fuels. The research effort specifically addressed broadening the experimental data base for ignition delay, burning rate, and oxidation kinetics at high pressures, and further refinement of chemical kinetic models so as to develop compositional specifications related to the above major and minor species. The foundation for the chemical kinetic modeling was the well validated mechanism for hydrogen and carbon monoxide developed over the last 25 years by Professor Frederick Dryer and his co-workers at Princeton University. This research furthered advance the understanding needed to develop practical guidelines for realistic composition limits and operating characteristics for HHC fuels. A suite of experiments was utilized that that involved a high-pressure laminar flow reactor, a pressure-release type high-pressure combustion chamber and a high-pressure turbulent flow reactor.

  6. Theory of homogeneous nucleation - A chemical kinetic view

    NASA Technical Reports Server (NTRS)

    Yang, C. H.; Qiu, H.

    1986-01-01

    A simple function with two undetermined parameters has been used in place of the Thomson-Gibbs relation to relate the activation energy of the vaporization reaction to cluster size. The parameters are iterated to assume optimum values in numerical computation so experimental data may be correlated. Calculations show this approach closely predicts and correlates available data for water, benzene, and ethanol. The nucleation formulism is redeveloped with an emphasis on the chemical kinetic view. Surface tension of the liquid and free energy of droplet formation are not used in its derivation.

  7. Pre-equilibrium approximation in chemical and photophysical kinetics

    NASA Astrophysics Data System (ADS)

    Rae, Margaret; Berberan-Santos, Mário N.

    2002-07-01

    For most mechanisms of chemical reactions and molecular photophysical processes the time evolution of the concentration of the intervening species cannot be obtained analytically. The pre-equilibrium approximation is one of several useful approximation methods that allow the derivation of explicit solutions and simplify numerical solutions. In this work, a general view of the pre-equilibrium approximation is presented, along with the respective analytical solution. It is also shown that the kinetic behavior of systems subject to pre-equilibration can be obtained by the application of perturbation theory. Several photophysical systems are discussed, including excimer formation, thermally activated delayed fluorescence, and external-heavy atom quenching of luminescence.

  8. Hard-sphere kinetic models for inert and reactive mixtures.

    PubMed

    Polewczak, Jacek

    2016-10-19

    I consider stochastic variants of a simple reacting sphere (SRS) kinetic model (Xystris and Dahler 1978 J. Chem. Phys. 68 387-401, Qin and Dahler 1995 J. Chem. Phys. 103 725-50, Dahler and Qin 2003 J. Chem. Phys. 118 8396-404) for dense reacting mixtures. In contrast to the line-of-center models of chemical reactive models, in the SRS kinetic model, the microscopic reversibility (detailed balance) can be easily shown to be satisfied, and thus all mathematical aspects of the model can be fully justified. In the SRS model, the molecules behave as if they were single mass points with two internal states. Collisions may alter the internal states of the molecules, and this occurs when the kinetic energy associated with the reactive motion exceeds the activation energy. Reactive and non-reactive collision events are considered to be hard sphere-like. I consider a four component mixture A, B, A (*), B (*), in which the chemical reactions are of the type [Formula: see text], with A (*) and B (*) being distinct species from A and B. This work extends the joined works with George Stell to the kinetic models of dense inert and reactive mixtures. The idea of introducing smearing-type effect in the collisional process results in a new class of stochastic kinetic models for both inert and reactive mixtures. In this paper the important new mathematical properties of such systems of kinetic equations are proven. The new results for stochastic revised Enskog system for inert mixtures are also provided. PMID:27545341

  9. Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

    NASA Technical Reports Server (NTRS)

    Hilst, G. R.

    1973-01-01

    Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

  10. A green approach towards adoption of chemical reaction model on 2,5-dimethyl-2,5-di-(tert-butylperoxy)hexane decomposition by differential isoconversional kinetic analysis.

    PubMed

    Das, Mitali; Shu, Chi-Min

    2016-01-15

    This study investigated the thermal degradation products of 2,5-dimethyl-2,5-di-(tert-butylperoxy) hexane (DBPH), by TG/GC/MS to identify runaway reaction and thermal safety parameters. It also included the determination of time to maximum rate under adiabatic conditions (TMR(ad)) and self-accelerating decomposition temperature obtained through Advanced Kinetics and Technology Solutions. The apparent activation energy (Ea) was calculated from differential isoconversional kinetic analysis method using differential scanning calorimetry experiments. The Ea value obtained by Friedman analysis is in the range of 118.0-149.0 kJ mol(-1). The TMR(ad) was 24.0 h with an apparent onset temperature of 82.4°C. This study has also established an efficient benchmark for a thermal hazard assessment of DBPH that can be applied to assure safer storage conditions. PMID:26368796

  11. Algebraic operator approach to gas kinetic models

    NASA Astrophysics Data System (ADS)

    Il'ichov, L. V.

    1997-02-01

    Some general properties of the linear Boltzmann kinetic equation are used to present it in the form ∂ tϕ = - †Âϕ with the operators Âand† possessing some nontrivial algebraic properties. When applied to the Keilson-Storer kinetic model, this method gives an example of quantum ( q-deformed) Lie algebra. This approach provides also a natural generalization of the “kangaroo model”.

  12. Chemical Kinetic Study of Toluene Oxidation Under Premixed and Nonpremixed Conditions

    SciTech Connect

    Costa, I D; Bozzelli, J W; Seiser, R; Pitz, W J; Westbrook, C K; Chen, C -; Fournet, R; Seshadri, K; Battin-Leclerc, F; Billaud, F

    2003-12-10

    A study was performed to elucidate the chemical-kinetic mechanism of combustion of toluene. A detailed chemical-kinetic mechanism for toluene was improved by adding a more accurate description of the phenyl + O{sub 2} reaction channels, toluene decomposition reactions and the benzyl + O reaction. Results of the chemical kinetic mechanism are compared with experimental data obtained from premixed and non-premixed systems. Under premixed conditions, predicted ignition delay times are compared with new experimental data obtained in shock tube. Also, calculated species concentration histories are compared to experimental flow reactor data from the literature. Under non-premixed conditions, critical conditions of extinction and autoignition were measured in strained laminar flows in the counterflow configuration. Numerical calculations are performed using the chemical-kinetic mechanism at conditions corresponding to those in the experiments. Critical conditions of extinction and autoignition are predicted and compared with the experimental data. Comparisons between the model predictions and experimental results of ignition delay times in shock tube, and extinction and autoignition in non-premixed systems show that the chemical-kinetic mechanism predicts that toluene/air is overall less reactive than observed in the experiments. For both premixed and non-premixed systems, sensitivity analysis was used to identify the reaction rate constants that control the overall rate of oxidation in each of the systems considered. Under shock tube conditions, the reactions that influence ignition delay time are H + O{sub 2} chain branching, the toluene decomposition reaction to give an H atom, and the toluene + H abstraction reaction. The reactions that influence autoignition in non-premixed systems involve the benzyl + HO{sub 2} reaction and the phenyl + O{sub 2} reaction.

  13. CHEMKIN-III: A FORTRAN chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics

    SciTech Connect

    Kee, R.J.; Rupley, F.M.; Meeks, E.; Miller, J.A.

    1996-05-01

    This document is the user`s manual for the third-generation CHEMKIN package. CHEMKIN is a software package whose purpose is to facilitate the formation, solution, and interpretation of problems involving elementary gas-phase chemical kinetics. It provides a flexible and powerful tool for incorporating complex chemical kinetics into simulations of fluid dynamics. The package consists of two major software components: an Interpreter and a Gas-Phase Subroutine Library. The Interpreter is a program that reads a symbolic description of an elementary, user-specified chemical reaction mechanism. One output from the Interpreter is a data file that forms a link to the Gas-Phase Subroutine Library. This library is a collection of about 100 highly modular FORTRAN subroutines that may be called to return information on equations of state, thermodynamic properties, and chemical production rates. CHEMKIN-III includes capabilities for treating multi-fluid plasma systems, that are not in thermal equilibrium. These new capabilities allow researchers to describe chemistry systems that are characterized by more than one temperature, in which reactions may depend on temperatures associated with different species; i.e. reactions may be driven by collisions with electrons, ions, or charge-neutral species. These new features have been implemented in such a way as to require little or no changes to CHEMKIN implementation for systems in thermal equilibrium, where all species share the same gas temperature. CHEMKIN-III now has the capability to handle weakly ionized plasma chemistry, especially for application related to advanced semiconductor processing.

  14. Thermoreversible associating polymer networks. I. Interplay of thermodynamics, chemical kinetics, and polymer physics

    NASA Astrophysics Data System (ADS)

    Hoy, Robert S.; Fredrickson, Glenn H.

    2009-12-01

    Hybrid molecular dynamics/Monte Carlo simulations are used to study melts of unentangled, thermoreversibly associating supramolecular polymers. In this first of a series of papers, we describe and validate a model that is effective in separating the effects of thermodynamics and chemical kinetics on the dynamics and mechanics of these systems, and is extensible to arbitrarily nonequilibrium situations and nonlinear mechanical properties. We examine the model's quiescent (and heterogeneous) dynamics, nonequilibrium chemical dynamics, and mechanical properties. Many of our results may be understood in terms of the crossover from diffusion-limited to kinetically limited sticky bond recombination, which both influences and is influenced by polymer physics, i.e., the connectivity of the parent chains.

  15. Development of Detailed Kinetic Models for Fischer-Tropsch Fuels

    SciTech Connect

    Westbrook, C K; Pitz, W J; Carstensen, H; Dean, A M

    2008-10-28

    Fischer-Tropsch (FT) fuels can be synthesized from a syngas stream generated by the gasification of biomass. As such they have the potential to be a renewable hydrocarbon fuel with many desirable properties. However, both the chemical and physical properties are somewhat different from the petroleum-based hydrocarbons that they might replace, and it is important to account for such differences when considering using them as replacements for conventional fuels in devices such as diesel engines and gas turbines. FT fuels generally contain iso-alkanes with one or two substituted methyl groups to meet the pour-point specifications. Although models have been developed for smaller branched alkanes such as isooctane, additional efforts are required to properly capture the kinetics of the larger branched alkanes. Recently, Westbrook et al. developed a chemical kinetic model that can be used to represent the entire series of n-alkanes from C{sub 1} to C{sub 16} (Figure 1). In the current work, the model is extended to treat 2,2,4,4,6,8,8-heptamethylnonane (HMN), a large iso-alkane. The same reaction rate rules used in the iso-octane mechanism were incorporated in the HMN mechanism. Both high and low temperature chemistry was included so that the chemical kinetic model would be applicable to advanced internal combustion engines using low temperature combustion strategies. The chemical kinetic model consists of 1114 species and 4468 reactions. Concurrently with this effort, work is underway to improve the details of specific reaction classes in the mechanism, guided by high-level electronic structure calculations. Attention is focused upon development of accurate rate rules for abstraction of the tertiary hydrogens present in branched alkanes and properly accounting for the pressure dependence of the ?-scission, isomerization, and R + O{sub 2} reactions.

  16. From prelife to life: how chemical kinetics become evolutionary dynamics

    PubMed Central

    Chen, Irene A.

    2015-01-01

    Conspectus Life is that which evolves. Living systems are the products of evolutionary processes and are capable of undergoing further evolution. A crucial question for the origin of life is the following: when do chemical kinetics become evolutionary dynamics? In this paper we review properties of ‘prelife’ and discuss the transition from prelife to life. We describe prelife as a chemical system where activated monomers can co-polymerize into macromolecules (such as RNA). These macromolecules are information carriers. Their physical and chemical properties depend to a certain extent on their particular sequence of monomers. We consider prelife as a logical precursor of life, where macromolecules are formed by copolymerization, but they are not capable of replication. Prelife can undergo ‘prevolutionary dynamics’. There can be mutation, selection and cooperation. Prelife selection, however, is blunt: small differences in rate constants lead to small differences in abundance. Life emerges with the ability of replication. In the resulting evolutionary dynamics selection is sharp: small differences in rate constants can lead to large differences in abundance. We also study the competition of different ‘prelives’ and find that there can be selection for those systems which ultimately give rise to replication. The transition from prelife to life can occur over an extended period of time. There may not have been a single moment which marks the origin of life. Instead prelife seeds many attempts for the origin of life. Eventually life takes over and destroys prelife. PMID:22335792

  17. Significance of Xenobiotic Metabolism for Bioaccumulation Kinetics of Organic Chemicals in Gammarus pulex

    PubMed Central

    2012-01-01

    Bioaccumulation and biotransformation are key toxicokinetic processes that modify toxicity of chemicals and sensitivity of organisms. Bioaccumulation kinetics vary greatly among organisms and chemicals; thus, we investigated the influence of biotransformation kinetics on bioaccumulation in a model aquatic invertebrate using fifteen 14C-labeled organic xenobiotics from diverse chemical classes and physicochemical properties (1,2,3-trichlorobenzene, imidacloprid, 4,6-dinitro-o-cresol, ethylacrylate, malathion, chlorpyrifos, aldicarb, carbofuran, carbaryl, 2,4-dichlorophenol, 2,4,5-trichlorophenol, pentachlorophenol, 4-nitrobenzyl-chloride, 2,4-dichloroaniline, and sea-nine (4,5-dichloro-2-octyl-3-isothiazolone)). We detected and identified metabolites using HPLC with UV and radio-detection as well as high resolution mass spectrometry (LTQ-Orbitrap). Kinetics of uptake, biotransformation, and elimination of parent compounds and metabolites were modeled with a first-order one-compartment model. Bioaccumulation factors were calculated for parent compounds and metabolite enrichment factors for metabolites. Out of 19 detected metabolites, we identified seven by standards or accurate mass measurements and two via pathway analysis and analogies to other compounds. 1,2,3-Trichlorobenzene, imidacloprid, and 4,6-dinitro-o-cresol were not biotransformed. Dietary uptake contributed little to overall uptake. Differentiation between parent and metabolites increased accuracy of bioaccumulation parameters compared to total 14C measurements. Biotransformation dominated toxicokinetics and strongly affected internal concentrations of parent compounds and metabolites. Many metabolites reached higher internal concentrations than their parents, characterized by large metabolite enrichment factors. PMID:22321051

  18. Oxidation Kinetics of Chemically Vapor-Deposited Silicon Carbide in Wet Oxygen

    NASA Technical Reports Server (NTRS)

    Opila, Elizabeth J.

    1994-01-01

    The oxidation kinetics of chemically vapor-deposited SiC in dry oxygen and wet oxygen (P(sub H2O) = 0.1 atm) at temperatures between 1200 C and 1400 C were monitored using thermogravimetric analysis. It was found that in a clean environment, 10% water vapor enhanced the oxidation kinetics of SiC only very slightly compared to rates found in dry oxygen. Oxidation kinetics were examined in terms of the Deal and Grove model for oxidation of silicon. It was found that in an environment containing even small amounts of impurities, such as high-purity Al2O3 reaction tubes containing 200 ppm Na, water vapor enhanced the transport of these impurities to the oxidation sample. Oxidation rates increased under these conditions presumably because of the formation of less protective sodium alumino-silicate scales.

  19. Polar Coordinate Lattice Boltzmann Kinetic Modeling of Detonation Phenomena

    NASA Astrophysics Data System (ADS)

    Lin, Chuan-Dong; Xu, Ai-Guo; Zhang, Guang-Cai; Li, Ying-Jun

    2014-11-01

    A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier—Stokes equations, including contribution of chemical reaction, via the Chapman—Enskog expansion. For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.

  20. Uncovering Oscillations, Complexity, and Chaos in Chemical Kinetics Using Mathematica

    NASA Astrophysics Data System (ADS)

    Ferreira, M. M. C.; Ferreira, W. C., Jr.; Lino, A. C. S.; Porto, M. E. G.

    1999-06-01

    Unlike reactions with no peculiar temporal behavior, in oscillatory reactions concentrations can rise and fall spontaneously in a cyclic or disorganized fashion. In this article, the software Mathematica is used for a theoretical study of kinetic mechanisms of oscillating and chaotic reactions. A first simple example is introduced through a three-step reaction, called the Lotka model, which exhibits a temporal behavior characterized by damped oscillations. The phase plane method of dynamic systems theory is introduced for a geometric interpretation of the reaction kinetics without solving the differential rate equations. The equations are later numerically solved using the built-in routine NDSolve and the results are plotted. The next example, still with a very simple mechanism, is the Lotka-Volterra model reaction, which oscillates indefinitely. The kinetic process and rate equations are also represented by a three-step reaction mechanism. The most important difference between this and the former reaction is that the undamped oscillation has two autocatalytic steps instead of one. The periods of oscillations are obtained by using the discrete Fourier transform (DFT)-a well-known tool in spectroscopy, although not so common in this context. In the last section, it is shown how a simple model of biochemical interactions can be useful to understand the complex behavior of important biological systems. The model consists of two allosteric enzymes coupled in series and activated by its own products. This reaction scheme is important for explaining many metabolic mechanisms, such as the glycolytic oscillations in muscles, yeast glycolysis, and the periodic synthesis of cyclic AMP. A few of many possible dynamic behaviors are exemplified through a prototype glycolytic enzymatic reaction proposed by Decroly and Goldbeter. By simply modifying the initial concentrations, limit cycles, chaos, and birhythmicity are computationally obtained and visualized.

  1. Infrared Absorption Spectroscopy and Chemical Kinetics of Free Radicals

    SciTech Connect

    Curl, Robert F; Glass, Graham

    2004-11-01

    This research was directed at the detection, monitoring, and study of the chemical kinetic behavior by infrared absorption spectroscopy of small free radical species thought to be important intermediates in combustion. Work on the reaction of OH with acetaldehyde has been completed and published and work on the reaction of O({sup 1}D) with CH{sub 4} has been completed and submitted for publication. In the course of our investigation of branching ratios of the reactions of O({sup 1}D) with acetaldehyde and methane, we discovered that hot atom chemistry effects are not negligible at the gas pressures (13 Torr) initially used. Branching ratios of the reaction of O({sup 1}D) with CH{sub 4} have been measured at a tenfold higher He flow and fivefold higher pressure.

  2. Model-Independent Bounds on Kinetic Mixing

    DOE PAGESBeta

    Hook, Anson; Izaguirre, Eder; Wacker, Jay G.

    2011-01-01

    New Abelimore » an vector bosons can kinetically mix with the hypercharge gauge boson of the Standard Model. This letter computes the model-independent limits on vector bosons with masses from 1 GeV to 1 TeV. The limits arise from the numerous e + e − experiments that have been performed in this energy range and bound the kinetic mixing by ϵ ≲ 0.03 for most of the mass range studied, regardless of any additional interactions that the new vector boson may have.« less

  3. Model Independent Bounds on Kinetic Mixing

    SciTech Connect

    Hook, Anson; Izaguirre, Eder; Wacker, Jay G.; /SLAC

    2011-08-22

    New Abelian vector bosons can kinetically mix with the hypercharge gauge boson of the Standard Model. This letter computes the model independent limits on vector bosons with masses from 1 GeV to 1 TeV. The limits arise from the numerous e{sup +}e{sup -} experiments that have been performed in this energy range and bound the kinetic mixing by {epsilon} {approx}< 0.03 for most of the mass range studied, regardless of any additional interactions that the new vector boson may have.

  4. Some aspects of mathematical and chemical modeling of complex chemical processes

    NASA Technical Reports Server (NTRS)

    Nemes, I.; Botar, L.; Danoczy, E.; Vidoczy, T.; Gal, D.

    1983-01-01

    Some theoretical questions involved in the mathematical modeling of the kinetics of complex chemical process are discussed. The analysis is carried out for the homogeneous oxidation of ethylbenzene in the liquid phase. Particular attention is given to the determination of the general characteristics of chemical systems from an analysis of mathematical models developed on the basis of linear algebra.

  5. Kinetic Analysis of Protein Folding Lattice Models

    NASA Astrophysics Data System (ADS)

    Chen, Hu; Zhou, Xin; Liaw, Chih Young; Koh, Chan Ghee

    Based on two-dimensional square lattice models of proteins, the relation between folding time and temperature is studied by Monte Carlo simulation. The results can be represented by a kinetic model with three states — random coil, molten globule, and native state. The folding process is composed of nonspecific collapse and final searching for the native state. At high temperature, it is easy to escape from local traps in the folding process. With decreasing temperature, because of the trapping in local traps, the final searching speed decreases. Then the folding shows chevron rollover. Through the analysis of the fitted parameters of the kinetic model, it is found that the main difference between the energy landscapes of the HP model and the Go model is that the number of local minima of the Go model is less than that of the HP model.

  6. Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

    NASA Technical Reports Server (NTRS)

    Goussis, D. A.; Lam, S. H.; Gnoffo, P. A.

    1990-01-01

    The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

  7. APOLLO: A computer program for the calculation of chemical equilibrium and reaction kinetics of chemical systems

    SciTech Connect

    Nguyen, H.D.

    1991-11-01

    Several of the technologies being evaluated for the treatment of waste material involve chemical reactions. Our example is the in situ vitrification (ISV) process where electrical energy is used to melt soil and waste into a ``glass like`` material that immobilizes and encapsulates any residual waste. During the ISV process, various chemical reactions may occur that produce significant amounts of products which must be contained and treated. The APOLLO program was developed to assist in predicting the composition of the gases that are formed. Although the development of this program was directed toward ISV applications, it should be applicable to other technologies where chemical reactions are of interest. This document presents the mathematical methodology of the APOLLO computer code. APOLLO is a computer code that calculates the products of both equilibrium and kinetic chemical reactions. The current version, written in FORTRAN, is readily adaptable to existing transport programs designed for the analysis of chemically reacting flow systems. Separate subroutines EQREACT and KIREACT for equilibrium ad kinetic chemistry respectively have been developed. A full detailed description of the numerical techniques used, which include both Lagrange multiplies and a third-order integrating scheme is presented. Sample test problems are presented and the results are in excellent agreement with those reported in the literature.

  8. APOLLO: A computer program for the calculation of chemical equilibrium and reaction kinetics of chemical systems

    SciTech Connect

    Nguyen, H.D.

    1991-11-01

    Several of the technologies being evaluated for the treatment of waste material involve chemical reactions. Our example is the in situ vitrification (ISV) process where electrical energy is used to melt soil and waste into a glass like'' material that immobilizes and encapsulates any residual waste. During the ISV process, various chemical reactions may occur that produce significant amounts of products which must be contained and treated. The APOLLO program was developed to assist in predicting the composition of the gases that are formed. Although the development of this program was directed toward ISV applications, it should be applicable to other technologies where chemical reactions are of interest. This document presents the mathematical methodology of the APOLLO computer code. APOLLO is a computer code that calculates the products of both equilibrium and kinetic chemical reactions. The current version, written in FORTRAN, is readily adaptable to existing transport programs designed for the analysis of chemically reacting flow systems. Separate subroutines EQREACT and KIREACT for equilibrium ad kinetic chemistry respectively have been developed. A full detailed description of the numerical techniques used, which include both Lagrange multiplies and a third-order integrating scheme is presented. Sample test problems are presented and the results are in excellent agreement with those reported in the literature.

  9. An open-source chemical kinetics network: VULCAN

    NASA Astrophysics Data System (ADS)

    Tsai, Shang-Min; Lyons, James; Heng, Kevin

    2015-12-01

    I will present VULCAN, an open-source 1D chemical kinetics code suited for the temperature and pressure range relevant to observable exoplanet atmospheres. The chemical network is based on a set of reduced rate coefficients for C-H-O systems. Most of the rate coefficients are based on the NIST online database, and validated by comparing withthermodynamic equilibrium codes (TEA, STANJAN). The difference between the experimental rates and those from the thermodynamical data is carefully examined and discussed. For the numerical method, a simple, quick, semi-implicit Euler integrator is adopted to solve the stiff chemical reactions, within an operator-splitting scheme for computational efficiency.Several test runs of VULCAN are shown in a hierarchical way: pure H, H+O, H+O+C, including controlled experiments performed with a simple analytical temperature-pressure profiles, so that different parameters, such as the stellar irradiation, atmospheric opacities and albedo can be individually explored to understand how these properties affect the temperaturestructure and hence the chemical abundances. I will also revisit the "transport-induced-quenching” effects, and discuss the limitation of this approximation and its impact on observations. Finally, I will discuss the effects of C/O ratio and compare with published work in the literature.VULCAN is written in Python and is part of the publicly-available set of community tools we call the Exoclimes Simulation Platform (ESP; www.exoclime.org). I am a Ph.D student of Kevin Heng at the University of Bern, Switzerland.

  10. A Steady-State Approximation to the Two-Dimensional Master Equation for Chemical Kinetics Calculations.

    PubMed

    Nguyen, Thanh Lam; Stanton, John F

    2015-07-16

    In the field of chemical kinetics, the solution of a two-dimensional master equation that depends explicitly on both total internal energy (E) and total angular momentum (J) is a challenging problem. In this work, a weak-E/fixed-J collisional model (i.e., weak-collisional internal energy relaxation/free-collisional angular momentum relaxation) is used along with the steady-state approach to solve the resulting (simplified) two-dimensional (E,J)-grained master equation. The corresponding solutions give thermal rate constants and product branching ratios as functions of both temperature and pressure. We also have developed a program that can be used to predict and analyze experimental chemical kinetics results. This expedient technique, when combined with highly accurate potential energy surfaces, is cable of providing results that may be meaningfully compared to experiments. The reaction of singlet oxygen with methane proceeding through vibrationally excited methanol is used as an illustrative example. PMID:25815602