Reaction networks procedure

The simultaneous determination of a great number of constants is a serious disadvantage of this procedure, since it considerably reduces the reliability of the solution. Experimental results can in some, not too complex cases be described well by means of several different sets of equations or of constants. An example would be the study of Wajc et al. (14) who worked up the data of Germain and Blanchard (15) on the isomerization of cyclohexene to methylcyclopentenes under the assumption of a very simple mechanism, or the simulation of the course of the simplest consecutive catalytic reaction A — B —> C, performed by Thomas et al. (16) (Fig. 1). If one studies the kinetics of the coupled system as a whole, one cannot, as a rule, follow and express quantitatively mutually influencing single reactions. Furthermore, a reaction path which at first sight is less probable and has not been therefore considered in the original reaction network can be easily overlooked. 

For one catalytic cycle, relaxation in the subspace = 0 is approximated by relaxation of a chain that is produced from the cycle by cutting the limiting step (Section 2). For reaction networks under consideration (with one cyclic attractor in auxiliary discrete dynamical system) the direct generalization works for approximation of relaxation in the subspace = 0 it is sufficient to perform the following procedures ... 

The procedure illustrated here may be applied directly to any system of n components which are interconverted via a single, central intermediate. A similar approach may be used to integrate the differential equations for any network of unimolecular reactions, plus a central intermediate. Their procedure involves using the Laplace-the kinetic analysis of reaction networks, including the three components... 

Multiple steady states as discussed in the previous subsection are related to the nonisothermicity of the CSTR. However, even in the isothermal case, a CSTR is known to be able to exhibit multiple steady states, periodic orbits, and chaotic behavior for sufficiently complex reaction network structures (see, e.g.. Gray and Scott, 1990). When the number of reactions is very large, the problem becomes a formidable one. In a series of papers (Feinberg 1987, 1988, and the literature quoted therein), Feinberg and his coworkers have developed a procedure for CSTRs that can be applied to systems with arbitrarily large numbers of reactants and reactions. The procedure is based on the deficiency concept discussed in Appendix C. 

I. Names, T. Vidoczy, L. Botar and D. G l, A Possible Construction of a Complex Chemical Reaction Network I. Definitions and Procedure for Construction, Theor. Chim. Acta 45 (1977) 215-223. 

A similar combinatorial approach can be also applied for complex catalytic reactions (Shalgunov et al., 1999 Temkin, 2000 Temkin et al., 1996 Zeigarnik and Valdes-Perez, 1998). It is based on the notion that the reaction network involves the formation of certain products (final or intermediate) from one precursor via different routes. Only a complete accounting of all junctions and connecting reactions in such a network treated together with the appropriate kinetic parameters can give a realistic representation of the overall process. Network compilation is a subject for a formal computerized procedure. 

Even "only" 120 rate coefficients and their activation energies are far more than were known or could be measured directly when this work was initiated. The majority had to be determined with statistical optimization routines from observed kinetics, a procedure that relies on the postulated reaction networks to be correct and the calculation not yielding a false or physically incorrect optimum. Not surprisingly, here and there a coefficient was later found to be seriously off the mark as better experimental techniques were developed and more information became available. This is an ongoing endeavor being continuously perfected. 

There are several advantages in adopting this new viewpoint of a reaction network for its analysis, visualization and reduction. The main advantage is that within this representation a reaction network becomes analogous to a general linear electrical circuit network, so that the procedures of electric circuit analysis are applicable. 

When the rate expression includes additive terms, each of which contains the above scaling problems, the situation becomes much more difficult. The fitting of data by a sum of exponential terms is a notoriously unsatisfactory procedure. Unfortunately this has to be done in many catalytic rate expressions as well as in many reaction networks. Although the problem of fitting this type of function has not been solved satisfactorily there are two established methods of dealing with it in a reasonably effective way ... 

With the generality of the procedures discussed in this chapter, it seems reasonable to suggest that GA algorithms are useful and promising for the determination of reaction mechanisms and rate coefficients of complex reaction networks. 

K. Yamasaki, A. Watanabe, T. Kakuda, and I. Tokue [ItU. J. Chem. Kinet., 30, 47-54 (1988)] developed a new procedure for analyzing the rates of very rapid reactions. In their analysis, these investigators considered a variety of reaction networks. A particularly interesting case is that represented by the following sequence of reactions ... 

Nemes, I., Vidoczy, T., Botar, L. Gal, D. (1977a). A possible construction of a complex chemical reaction network, I. Definitions and procedure for construction, Theoret. Chim. Acta Berl), 45, 215-23. 

This preparation procedure also creates solid-state phases that are key to the performance of the Mo-V-Nb-Te-0 catalyst for propane ammoxidation. High activity and selectivity result when the x-ray powder diffraction pattern shows the presence of specific diffraction lines attributed to two separate phases denoted as Ml and M2 by Mitsubishi Chemical Corp. The diffraction lines assigned to these two phases are given in Table 7 (146). The coexistence of these two phases is viewed as key to the successful functioning of the catalyst. Specifically, the Ml phase is purportedly responsible for the oxidative dehydrogenation of propane to propylene, the key intermediate in the reaction network. This reaction sequence, in which the first step is the formation of a propylene intermediate, is the same as noted previously with other propane ammoxidation catalysts, most notably with the V-Sb-0 catalyst (see above). The M2 phase of the Mo-V-Nb-Te-0 catalyst is reportedly the center for the selective ammoxidation of the propylene intermediate to acrylonitrile. As the first-formed intermediate, propylene is apparently the source of all the observed reaction products. Although a detailed kinetic analysis has not been presented, a cursory report, published in Japan, summarized the kinetic experiments for the conversion of propane and propylene over a... 

The situation is even more complex in the case of reaction related parameters. Independent of the appropriate reactor model, for complex reactions not only the reaction scheme but also the adequate type of rate equation for each reaction step has to be chosen before the parameters can be estimated. As a rule, in reaction engineering only the analytically measureable reactants (and reactions) should form the basis of the reaction network in contrast to physico-chemical research where the true" reaction mechanism (involving radical intermediates or active complexes) is sought. Certainly, the stoichiometry of the experimental product spectrum is important, but also the concentration/reactiontime dependencies like those given in Figure 6 are helpful. In contrast to parameter estimation and model discrimination, there exists no unique and straight forward analytic procedure for the built-up of even a simplified reaction scheme. The intuition of the chemical reaction engineer is therefore heavily relied upon. 

---