Model formulation, reaction-rate equations

It is usually assumed in the derivation of isothermal rate equations based on geometric reaction models, that interface advance proceeds at constant rate (Chap. 3 Sects. 2 and 3). Much of the early experimental support for this important and widely accepted premise derives from measurements for dehydration reactions in which easily recognizable, large and well-defined nuclei permitted accurate measurement. This simple representation of constant rate of interface advance is, however, not universally applicable and may require modifications for use in the formulation of rate equations for quantitative kinetic analyses. Such modifications include due allowance for the following factors, (i) The rate of initial growth of small nuclei is often less than that ultimately achieved, (ii) Rates of interface advance may vary with crystallographic direction and reactant surface, (iii) The impedance to water vapour escape offered by... 

The principles underlying the formulation of rate equations applicable to the decompositions of solids are presented in Section 5.4. In summary, these result in the replacement of the concentration terms generally applicable in homogeneous rate processes by geometric or diffusion parameters. It is possible, in principle, to formulate a further set of kinetic models that describe concurrent reactions proceeding... 

This representation summarizes the conceptual foundation for the theory of reaction kinetics in solids. Essential models used for the formulation of rate equations include those described in the following text more detailed accounts are given in the references cited earlier. 

Vanag and Epstein have formulated a four-variable model to understand pattern formation in the BZ-AOT system [449,453]. Their model builds on the Oregonator, see Sect. 1.4.8. It assumes that the chemistry within the water core of the droplets is well described by the two-variable Oregonator rate equations (1.131). It further assumes that the species in the oil phase are inert, since they lack reaction partners, the key reactants all being confined to the aqueous core of the droplets. Consequently, only transfer reactions occur for the activator B1O2 and inhibitor Br2 in the oil phase. The rate terms for the two transfer reactions are added to the rate terms of the two-variable Oregonator model. The reaction-diffusion equations of the four-variable model of the BZ-AOT system are given in nondimensionalized form by... 

Minimizing the cycle time in filament wound composites can be critical to the economic success of the process. The process parameters that influence the cycle time are winding speed, molding temperature and polymer formulation. To optimize the process, a finite element analysis (FEA) was used to characterize the effect of each process parameter on the cycle time. The FEA simultaneously solved equations of mass and energy which were coupled through the temperature and conversion dependent reaction rate. The rate expression accounting for polymer cure rate was derived from a mechanistic kinetic model. 

Erdey-Gruz and Volmer (2) derived the current-potential relationship in 1930 using the Arrhenius equation (1889) for the reaction rate constant and introduced the transfer coefficient. They also formulated the nucleation model of electrochemical crystal growth. 

In practice, a gray-box model is developed in steps. One early step is to decide which variables and interactions to include. This is often done by the sketching of an interaction-graph. It must then be decided if a variable should be a state or a dependent variable, and how the interactions should be formulated. In the case of metabolic reactions, the expression forms for the reactions have often been characterized in in-vitro experiments. If this has been done, there are also often in-vitro estimates of the kinetic parameters. For enzymatic networks, however, such in-vitro studies are much more rare, and it is hence typically less known which expression to choose for the reaction rates, and what a good estimate for the kinetic parameters is. In any case, the standard method of combining reaction rates, r,-, and an interaction graph into a set of differential equations is to use the stoichiometric coefficients, Sij... 

Formulating appropriate rate laws for CO adsorption, OH adsorption and the reaction between these two surface species, a set of four coupled ordinary differential equations is obtained, whereby the dependent variables are the average coverages of CO and OH, the concentration of CO in the reaction plane and the electrode potential. In accordance with the experiments, the model describes the S-shaped I/U curve and thus also bistability under potentiostatic control. However, neither oscillatory behavior is found for realistic parameter values (see the discussion above) nor can the nearly current-independent, fluctuating potential be reproduced, which is observed for slow galvanodynamic sweeps (c.f. Fig. 30b). As we shall discuss in Section 4.2.2, this feature might again be the result of a spatial instability. 

Most often the problem encountered in describing the particle formation quantitatively involves formulating a mathematical model that is close enough to the real process but not so con Iex as to make a solution impossible. Common to the chain-reaction models is the fact that a simple solution for the particle number is difficult to obtain analytically without considerable simplification. Below are given the rate equations for the chain model based on a previous treatment (Hansen and Ugelstad, 1978). In addition the present model includes the desorption of monomer radicals and the different absorption rates, as outlined in Table I. 

The demonstration that yield-time data for a particular rate process are satisfactorily represented by a rate equation from the set in Table 3.3. is a usefiil initial result (as in all kinetic analyses), but formulation of a reaction model requires more support than a statistical comparison with the fit of the same data to other rate equations. Some important aspects of kinetic analysis are summarized below. 

The objective of Chaps. 10 and 11 is to combine intrinsic rate equations with intrapellet and fluid-to-pellet transport rates in order to obtain global rate equations useful for design. It is at this point that models of porous catalyst pellets and effectiveness factors are introduced. Slurry reactors offer an excellent example of the interrelation between chemical and physical processes, and such systems are used to illustrate the formulation of global rates of reaction. 

Thermogravimetry is an attractive experimental technique for investigations of the thermal reactions of a wide range of initially solid or liquid substances, under controlled conditions of temperature and atmosphere. TG measurements probably provide more accurate kinetic (m, t, T) values than most other alternative laboratory methods available for the wide range of rate processes that involve a mass loss. The popularity of the method is due to the versatility and reliability of the apparatus, which provides results rapidly and is capable of automation. However, there have been relatively few critical studies of the accuracy, reproducibility, reliability, etc. of TG data based on quantitative comparisons with measurements made for the same reaction by alternative techniques, such as DTA, DSC, and EGA. One such comparison is by Brown et al. (69,70). This study of kinetic results obtained by different experimental methods contrasts with the often-reported use of multiple mathematical methods to calculate, from the same data, the kinetic model, rate equation g(a) = kt (29), the Arrhenius parameters, etc. In practice, the use of complementary kinetic observations, based on different measurable parameters of the chemical change occurring, provides a more secure foundation for kinetic data interpretation and formulation of a mechanism than multiple kinetic analyses based on a single set of experimental data. 

Eq. (15) is the rate equation of the reaction (also called the kinetic model ). The formulation of such a differential equation for all reacting substances is the basic step in describing the kinetics of chemical/biochemical reactions. These rate equations include concentration values of the relevant reaction partners and kinetic parameters such as the rate constant k. An investigation of enzyme kinetics includes the measurement of reaction rates, the choosing of an appropriate kinetic model and the identification of the kinetic parameters. 

Avoidance of the appearance of the cold-boundary difficulty can be achieved only by revision of the physical model on which the mathematical formulation is based. There are many ways in which this can be done. In one approach, von Karman and Millan [17] replace t = 0 by t = r, (where 0 < T < 1) as the position at which to apply the cold-boundary condition e = 0. The introduction of this artifice is equivalent to employing the physical concept of an ignition temperature I], below which the chemical reaction rate vanishes [see equations (19) and (33)]. An alternative procedure for determining a finite, nonzero value for A involves the assumption that a flame holder serves as a weak heat sink [18], removing an amount of heat per unit area per second equal to (XdT/dx)i = mCp(T - TQXdx/d )i. Since e = 0 at the flame holder (that is, there is no reaction upstream), equation (19) implies that dx/d i = rj, whence the heat-sink concept is seen to be mathematically equivalent to the ignition-temperature concept at the cold boundary. 

Heterogeneous reactions. Components of water or air pollution are usually in the fluid phase. Hence we may write equations such as equations 6.2, 6.6, 6.8, and 6.12 for the fluid. The fluid may have non-permeable boundaries (the reactor walls) and permeable boundaries (entrances and exits of the system as well as catalytic surfaces where mass fluxes must be equal to the superficial reaction rates). Usually, these reaction rates are modeled as pseudo-homogeneous and, moreover, concentration measurements are almost always made in the fluid phase. Heterogeneous reactions are the result of a process that occurs at phase interfaces. This means that for the differential equation written for the fluid phase, heterogeneous reactions (surface reactions, for example) are just boundary conditions. The problem is very simple to formulate at steady state and at the boundary of an active surface, the normal mass or molar fluxes must be made equal to the heterogeneous, superficial reaction rate. Then,... 

---