Kinetic models and their analysis

Fields and Noyes considered the case of the Belousov-Zhabotinskii reaction proceeding in well stirred solution, at a constant temperature and pressure, in a closed system. They further assumed that the changes in [A] = [Br03-] may be disregarded in a first approximation, that the product P = HOBr has no effect on the reaction kinetics, and that the reversibility of the above reaction can be neglected. 

---

A reverse kinetic problem consists in identifying the type of kinetic models and their parameters according to experimental (steady-state and unsteady-state) data. So far no universal method to solve reverse problems has been suggested. The solutions are most often obtained by selecting a series of direct problems. Mathematical treatment is preceded by a qualitative analysis of experimental data whose purpose is to reduce drastically the number of hypotheses under consideration [31]. 

T. Ogiso, M. Iwaki, T. Tanino, A. Yono, A. Ito, In vitro Skin Penetration and Degradation of Peptides and Their Analysis Using a Kinetic Model , Biol. Pharm. Bull. 2000, 23, 1346-1351. 

The first two sections of Chapter 5 give a practical introduction to dynamic models and their numerical solution. In addition to some classical methods, an efficient procedure is presented for solving systems of stiff differential equations frequently encountered in chemistry and biology. Sensitivity analysis of dynamic models and their reduction based on quasy-steady-state approximation are discussed. The second central problem of this chapter is estimating parameters in ordinary differential equations. An efficient short-cut method designed specifically for PC s is presented and applied to parameter estimation, numerical deconvolution and input determination. Application examples concern enzyme kinetics and pharmacokinetic compartmental modelling. 

Metabolic flux analysis Cellular metabolites and metabolic fluxes can be combined into a series of balance equations, not unlike a series of (bio)chemical reactions in a kinetic model. Metabolic flux analysis is the description of the components and their connections in a metabolic network. 

This chapter gave an overview of how to simplify complex processes sufficiently to allow the use of analytical models for their analysis and optimization. These models are based on mass, momentum, energy and kinetic balance equations, with simplified constitutive models. At one point, as the complexity and the depth of these models increases by introducing more realistic geometries and conditions, the problems will no longer have an analytical solution, and in many cases become non-linear. This requires the use of numerical techniques which will be covered in the third part of this book, and for the student of polymer processing, perhaps in a more advanced course. 

Biosimulation has a dominant role to play in systems biology. In this chapter, we briefly outline two approaches to systems biology and the role that mathematical models has to play in them. Our focus is on kinetic models, and silicon cell models in particular. Silicon cell models are kinetic models that are firmly based on experiment. They allow for a test of our knowledge and identify gaps and the discovery of unanticipated behavior of molecular mechanisms. These models are very complicated to analyze because of the high level of molecular-mechanistic detail included in them. To facilitate their analysis and understanding of their behavior, model reduction is an important tool for the analysis of silicon cell models. We present balanced truncation as one method to perform model reduction and apply it to a silicon cell model of glycolysis in Saccharomyces cerevisiae. 

Gunn and coworkers (1,2) were the first to propose a steady-state model, and their predictions agreed very well with the Hanna UCG test results. In an analysis of the different versions of permeation models that have appeared in the literature, Haynes (3) judged the steady-state model superior for most applications since reaction kinetics are taken into account and only a modest computational effort is required. Despite these desirable features, applications of the steady-state model have not been as widespread as one might anticipate. 

At the same time calculations on the modified MEIS are possible without additional kinetic models and do not require extra experimental data for calculations, which makes it possible to use less initial information and obviously reduces the time and labor spent for computing experiment. Furthermore, there arise principally new possibilities for the analysis of methods to mitigate emissions from pulverized-coal boilers, since at separate modeling of different mechanisms of NO formation the measures taken can result in different consequences for each in terms of efficiency. Consideration of kinetic constraints in MEIS will substantially expand the sphere of their application to study other methods of coal combustion (fluidized bed, fixed bed, etc.) and to model processes of forming other pollutants such as polyaromatic hydrocarbons, CO, soot, etc. 

The analysis of kinetic data and their interpretation by models have been discussed in books dealing with pure chemical kinetics [1—15],... 

We suppose that the value method of analyzing the reaction kinetic models may enhance considerably the efficiency of a similar procedure, as it improves the capabilities to retrieve information from the kinetic models of reactions. On the other hand, this method enables to rank the steps, and accordingly, the rate constants of the kinetic model, by their sensitivity in describing specific experimental results. This will permit to determine acceptability of applying the rate constants with a prescribed accuracy. One can look through the analysis of imcertainty of the kinetic model due to the uncertainty of rate constant of steps in references [58,65,66]. 

The value analysis of the kinetic models, and as a consequence, the ranking of individual chemical transformations by their contribution into the resultant indicator of a chemical multisteps reaction, provides a basis which enables to retrieve information on the reactivity of the reacting particles. For example, it is possible to identify what chemical reactions involving an initial reagent and its intermediates play an important role and then to forecast how the change in the molecular structure of a reaction species will influence on the result of a complex chemical process. 

A third example of the role of kinetic modelling in the analysis of practical combustion systems at high temperatures is in simulations of pulse combustion. Pulse combustors have been used for many years in propulsion, industrial processes such as drying, and in home furnaces, where they are valued for their thermal efficiency and low NOx production rates. However, until kinetic modelling was used to analyze the role of thermal ignition in these systems, it had been impossible to understand the principles of pulse combustion and the real reasons for their good performance. 

Lutz and Matyjaszewski [18] have followed the evolution of the bromine end functionality during the bulk ATRP of styrene, in the presence of the CuBr/4,4 -di-(5-nonyl)-2,2 -bipyridine catalyst. The retention of the bromide chain-end functionality was monitored through the withdrawal of aliquots at given times from the polymerization mixture and their analysis by H NMR (600 MHz). A decrease in the functionality with conversion was observed, significant at high monomer conversion (90%). The experimental data allowed, by comparison with a kinetic model of styrene ATRP, better understanding of the side reactions that led to the loss in catalyst functionality and helped in the design of the most suitable reaction conditions in order to optimize the reaction kinetics and end-product properties. 

In comparative studies, Mydlarz and Jones (1991, 1994) point out that use of different size-dependent crystal growth rate models and methods for the analysis of curved log-linear MSMPR CSD data can lead to gross variations in the magnitudes of crystallization kinetics so determined. The common occurrence of such anomalous MSMPR CSD data and their analysis by differing methods may in large part explain variations in reported crystallization kinetics. The use of exponential size-dependent equations provides improved data fitting and CSD prediction compared with all other models tested, but given their purely empirical basis they should nevertheless still be used with caution. 

---