Empirical rate law models

V.Gaspar and K.Showalter, The Oscillatory Reaction. Empirical Rate Law Model and Detailed Mechanism, Journal of the American Chemical Society,... 

One alternative to working with a full elementary step mechanism that has proved useful in the study of a number of complex reactions involves the use of an empirical rate law model. In this approach, we describe the overall reaction not in terms of a complete set of elementary steps, but in terms of a set of processes that add up to the full reaction and for each of which a rate law is known experimentally. The set of rate equations that describes the system is then the set of rate laws for all of these processes. Typically, there will be far fewer rate equations than for a full mechanism, but each equation will have a more complex form, since we are... 

Intermediate in complexity and accuracy between true mechanisms and abstract models are empirical rate law models. Here, the modeler eschews the identification of elementary steps and, instead, works with experimentally established rate laws for the component overall stoichiometric processes that make up a particular reaction. Each process may consist of several elementary steps and involve many reaction intermediates, but it enters the model only as a single empirical rate equation, and only those species that appear in the rate equation need be included in the model. Assuming that the empirical rate laws have been accurately determined, this approach will give results for the species contained in the rate laws that are identical to the results from the full mechanism, so long as no intermediate builds up to a significant concentration and so long as the component processes are independent of one another. This last requirement implies that no intermediate that is omitted from the model is involved in more than one process, and that there are no cross-reactions between intermediates involved in different processes. 

Caspar and Showalter identify five overall stoichiometric processes that make up the total reaction. These are summarized in Table 5.1. The rate laws used for processes A, B, and C have been simplified from the full empirical multiterm rate laws (Reynolds, 1958 von Bunau and Eigen, 1962 Liebhafsky and Roe, 1979) by taking into account the conditions of the oscillatory EOE reaction. When the stoichiometries of the model equations for processes A and D are simplified by dividing by 3, we obtain the following set of p.seudo-elementary reactions as our final empirical rate law model. By chance, the model obeys mass action kinetics, though the rate constants are products of rate constants for the nonelementary component processes and the constant concentrations. 

We can model this behaviour with a set of three reactions and their differential equations, (a) In the first reaction the sheep are breeding. Note, that there is a constant supply of grass and this reaction could go on forever. As it is written, this reaction violates the law of conservation of mass, it is only an empirical rate law. In a second reaction (b), wolves eat sheep and breed themselves. The third reaction (c) completes the system, wolves have to die a natural death. 

Mox represents the metal ion catalyst in its oxidised form (Ceexperimentally determined empirical rate law and does clearly not comprise stoichiometrically correct elementary processes. The five reactions in the model provide the means to kinetically describe the four essential stages of the BZ reaction ... 

Function (2.2) can be considered as an empirical model used to best fit the experimental concentration-time data. In practice, laws different from (2.2) are also encountered, especially when dependence on the concentration is considered however, a simple theory based on the kinetic theory of gases can only explain the simplest of these empirical rate laws. The general idea of this theory is that reaction occurs as a consequence of a collision between adequately energized molecules of reactants. The frequency of collision of two molecules can explain simple reaction... 

The fact that some kinetic profiles are fitted by sums of exponentials, and others are fitted by power functions, suggests that different types of basic mechanisms are at work. In fact, as concluded in Chapter 7, while kinetics from homogeneous media can be fitted by sums of exponentials, heterogeneity shapes kinetic profiles best represented by empirical power-law models. Conversely, when power laws fit the observed data, they suggest that the rate at which a material leaves the site of a process is itself a function of time in the process, i.e., age of material in the process. 

In a similar spirit, Lengyel et al. (1990) have proposed and analyzed a particularly elegant model of another oscillating reaction, the chlorine dioxide-iodine-malonic acid (CIO2 -I2 -MA) reaction. Their experiments show that the following three reactions and empirical rate laws capture the behavior of the system ... 

As already noted, these reactions are characterized by hyperbolic rate forms. Empirical power law models can also be used, but their applicability is restricted to the ranges of parameter values used in their formulation. Note in the table that mechanically agitated contactors (MACs) and bubble-column reactors (BCRs) are the most commonly used reactors. The design of such reactors is considered in Chapter 16. 

An elegant example of the empirical rate law approach is the model developed by Gaspar and Showalter (1990) for the mixed Landolt or Edblom-Orbdn-Epstein (EOE) reaction (Edblom et ah, 1986). The reaction of iodate, sulfite, and ferrocyanide ions in a CSTR gives rise to sustained oscillation. In a flow reactor, the concentrations of two of these input reactants -iodate and ferrocyanide—remain nearly constant, while the pH and the concentrations of sulfite. 

Empirical Stoichiometry and Rate Law Model Stoichiometry and Rate Law... 

In solution kinetics an important distinction between solvent and reactants must be maintained. The solvent is continually interacting with the reactants if such interactions were incorporated in the mechanistic model all solution mechanisms would perforce be multimolecular. By considering the solvent as a medium and not as a participant in the reaction (unless, of course, solvent actually takes part in a reaction step), the problem of mechanism is greatly simplified. In this sense isomerizations, rearrangements, and conformational changes, like the chairi chair2 interconversion in cyclohexane, are first-order reactions for which the empirical rate law is a direct indication of the only important elementary step. Most solution reactions proceed via bimolecular steps. There are countless examples for which only one such step is needed and for which the rate law reflects that process. 

Strong evidence for this model is provided by an elegant substrate-catalyzed single-turnover experiment presented by Sadow and coworkers [81]. The four-coordinate Cs-symmetric magnesium tris(oxazoline)borate alkyl complex [To MgMe] (To = tris(4,4-dimethyl-2-oxazolinyl)phenylborate) excludes additional ligands from the coordination sphere and is a competent precatalyst for the hydroamination of l-aminopent-4-enes. In line with the expectations furnished by the aforementioned systems, the reaction of l-amino-2,2-diphenylpent-4-ene follows the empirical rate law v = / obs[aminoalkene] [catalyst], with a KIE of... 

Empirical Models. In the case of an empirical equation, the model is a power law rate equation that expresses the rate as a product of a rate constant and the reactant concentrations raised to a power (17), such as... 

Hydrogen adsorption was also described as irreversible in our previous mechanism,10 and an empirical kinetic law was used to describe the rate of this step. However, a deeper analysis of literature data revealed that this step is likely in equilibrium, too. On the basis of this evidence, the previously developed model has been modified in this work in order to improve the physical consistency of the proposed mechanism. 

These solutions to the one-dimensional advection-diffusion model can be used to estimate reaction rate constants Ck) from the pore-water concentrations of S, if and s are known. More sophisticated approaches have been used to define the reaction rate term as the sum of multiple removals and additions whose functionalities are not necessarily first-order. Information on the reaction kinetics is empirically obtained by determining which algorithmic representation of the rate law best fits the vertical depth concentration data. The best-fit rate law can then be used to provide some insight into potential... 

When a simple, fast and robust model with global kinetics is the aim, the reaction kinetics able to predict correctly the rate of CO, H2 and hydrocarbons oxidation under most conditions met in the DOC consist of semi-empirical, pseudo-steady state kinetic expressions based on Langmuir-Hinshelwood surface reaction mechanism (cf., e.g., Froment and Bischoff, 1990). Such rate laws were proposed for CO and C3H6 oxidation in Pt/y-Al203 catalytic mufflers in the presence of NO already by Voltz et al. (1973) and since then this type of kinetics has been successfully employed in many models of oxidation and three-way catalytic monolith converters... 

The works of Whitlow and Roth (1988) as well as of Beltran et al. (1990) employed an empirical approach for modeling. This procedure uses a global rate law of nh order for the observed disappearance of all target contaminants... 

Process-scale models represent the behavior of reaction, separation and mass, heat, and momentum transfer at the process flowsheet level, or for a network of process flowsheets. Whether based on first-principles or empirical relations, the model equations for these systems typically consist of conservation laws (based on mass, heat, and momentum), physical and chemical equilibrium among species and phases, and additional constitutive equations that describe the rates of chemical transformation or transport of mass and energy. These process models are often represented by a collection of individual unit models (the so-called unit operations) that usually correspond to major pieces of process equipment, which, in turn, are captured by device-level models. These unit models are assembled within a process flowsheet that describes the interaction of equipment either for steady state or dynamic behavior. As a result, models can be described by algebraic or differential equations. As illustrated in Figure 3 for a PEFC-base power plant, steady-state process flowsheets are usually described by lumped parameter models described by algebraic equations. Similarly, dynamic process flowsheets are described by lumped parameter models comprising differential-algebraic equations. Models that deal with spatially distributed models are frequently considered at the device... 

The Power Law model (excluding temperature dependence) is a two-parameter empirical model proposed by Ostwald and de Waele (39). It is based on the experimental observation that by plotting In tj(y)vs. In(7), a straight line is obtained in the high shear rate region for... 

In chemical kinetics, semi-empirical non-linear models for reaction rates are commonly used. For example, Boudart [3] summarizes the laws of reaction rates in the formulae... 

The first-order rate coefficient, k, of this pseudo-elementary process is assumed to vary with temperature according to an Arrhenius law. Model parameters are the stoichiometric coefficients v/ and the Arrhenius parameters of the rate coefficient, k. The estimation of the decomposition rate coefficient, k, requires a knowledge of the feed conversion, which is not directly measurable due to the complexity of analyzing both reactants and reaction products. Thus, a supplementary empirical relationship is needed to relate the feed conversion (conversion of A) to some experimentally accessible variable (Ross and Shu have chosen the yield of C3 and lighter hydrocarbons). It is observed that the rate coefficient, k, is not constant and decreases with increasing conversion. Furthermore, the zero-conversion rate coefficient depends on feed specifications (such as average carbon number, hydrogen content, isoparaffin/normal-paraffin ratio). Stoichiometric coefficients are also correlated with conversion. Of course, it is necessary to write supplementary empirical relationships to account for these effects. 

The proportionality factor k is called the experimental rate constant with catalytic reactions, this constant is frequently a complex quantity which may be a product of rate constants of several steps or may include equilibrium constants of the fast steps. The exponents m,n,... in the power-law equations may be any fraction or small integer (positive, negative or zero). The constants K, in the denominator of the equations of type 7 are often but not always related to adsorption equilibrium constants. While they have to be evaluated from the experimental data, the values of exponents a, b, a, / and d arc derived from the assumed mechanism in the case of model rate equations. In empirical rate equations these constants can attain any value (fractional or small integer, usually positive) and... 

It should be noted that the effects of fillers may be incorporated into the cure and shear-rate effects. The main forms of combined-effects model consist of WLF, power-law or Carreau shear effects, Arrhenius or WLF thermal effects and molecular, conversion or empirical cure effects. Nguyen (1993) and Peters et al. (1993) used a modified Cox-Merz relationship to propose a modified power-law model for highly filled epoxy-resin systems. Nguyen (1993) also questions the validity of the separability of thermal and cure effects in the derivation of combined models. 

Surprisii y, although empirical, these rate laws give satisfoctoty results in simulations. This si psts that sinqtle rate eiqiressions are probably. sufficient for simulation purposes. However, introducing new reactions (new hydrocarbons, steam-reformii, etc.) in the model would require time-consumii experiments. The same problem arises when acUqiting the fiequemy foctor, koi, and the apparent activation enei, Eai, to the catalyst under use. There is thus a need for a simple method allowing the estimation of... 

---