Global reaction models

Global reactions models are those where the reactants and products are defined by available analytical separation schemes. They generally represent the interconversion of lumps or pseudocomponents, i.e., aggregates of many molecules with common attributes. In general, the characteristics which assemble given sets of molecules into a lump will not be reactivity. For example, perhaps the two most commonly found globally lumped models are based on boiling point or solubility characteristics. 

Global reactions models have historical significance. Their formulation is largely due to the global nature of analytical output of hydrocarbon mixtures in the 1940s and 50s. Additional advantages of these models include their easy formulation and solution. 

The most detailed modelling approach summarized in Figure 1 is found at the mechanistic level. These models are explicit accounts of the chemistry of elementary steps. Thus the hierarchy of the levels, i.e., reaction models in Figure 1, now becomes quite clear. Mechanistic models, which provide the temporal and many times spatial variation of the composition of each component and reaction intermediate, are based at the lowest modelling level. Their output, however, is typically phrased in terms of ensembles of stable molecular constituents which is more characteristic of the intermediate level molecular models. The molecular models, in turn, require subsequent organization in order to connect to the global reaction models and relevant product fractions at the top or global level. 

In summary, the prospects of mechanistic-based modelling appear to be several years away. However, it is clear that global reaction models are already losing favor because of their inability to answer current and future questions. We may look back and see the Clean Air Act as the beginning of the need for molecule-based models, where the input and output has to be commensurate with the nature of the questions being asked. 

Nitrogen oxides are mainly formed through three paths. In the fuel-N path nitrogen containing species in the fuel can form NO or N2. In the two other paths the fixation of N2 in air is involved. One of these is the well known thermal-NO path, where radicals react at high temperatures with N2 to form NO. The other one is the often called prompt-NOx path, where hydrocarbon radicals react with N2. For most of these paths global reaction models are available (Mitchell and Tarbell, 1982 De Soete, 1974 Bowman, 1979) which can be also found in most current CFD codes. If a certain path is... 

The aim of this study is to develop model reaction for the characterization of the acidity and basicity of various transition aluminas, the experimental conditions being close to that for catalysis use. Among various model reactions, the transformation of cyclopentanol and cyclohexanone mixture was chosen for this work. Indeed, this reaction was well known for estimating simultaneously the acid-base properties of oxide catalysts [1], Two reactions take place the hydrogen transfer (HT) on basic sites and the alcohol dehydration (DEH) on acid sites. The global reaction scheme is shown in Figure 1. 

The time-to-rimaway can be calculated using dT/dt and Tmax values. This calculated time is a measure of the possible global reaction rate. ARC experimental results may also be used to develop required mathematical models for process design. 

In the washcoat, reaction rates are modeled via global reaction mechanisms. In such a global or macrokinetic reaction mechanism, several microkinetic adsorption, reaction and desorption steps are lumped together, reducing the overall number of kinetic parameters considerably. For some catalysts,... 

The employed physical and chemically based modeling approach enables a relatively large variation of geometrical catalyst parameters, e.g. catalyst length, diameter, etc., with the global reaction kinetics and therefore the overall... 

The aim of the present section is to illustrate the procedures employed for the derivation of dynamic kinetic models appropriate for simulation of exhaust aftertreatment devices according to the converter models illustrated in the previous section. In particular, it will be shown how to derive global reaction kinetics which are based on a fundamental study aimed at the elucidation of the reaction mechanism. In principle, this approach enables a greater model adherence to the real behavior of the reacting system, which should eventually afford better results when validating the model in a wide range of operating conditions, as typically required for automotive applications. 

In this section we shall discuss the development of a global kinetic model for DOC. The basic model reactions considered in the DOC model are summarized in Table II. Here the real HC mixture is modeled by two characteristic hydrocarbons—propene and decane. Propene represents more reactive, light hydrocarbons, which practically do not adsorb during cold start, while decane is a representative of heavier hydrocarbons with significant adsorption on... 

For a comprehensive kinetic description of the NH3 + N0/N02 reacting system in a wide range of temperatures and N02/N0X feed ratios a global kinetic model was developed, based on the whole set of reactions in Table V. 

Fig. 1.15 Schematic of the energy curves in the Marcus-Hush model with a single, global reaction coordinate q such that the potential energy hypersurface reduces to two parabolas and the activation energy can he calculated from the intersection point between them. The electronic coupling (Sect. 1.7.2.2) and the continuum of electronic levels in the metal electrode (Sect. 1.7.2.1) are not shown...

To have a quantitative idea of the higher unit surface area activity of the Monolith catalyst, rate constants based on surface area were considered essential to know. To accomplish this, the global reaction kinetics of desulfurization and denitrogen-ation were determined. For the desulfurization the following three kinetic models, as suggested in the literature, were tested to determine which best represented the data of this study. 

In modeling studies of the decomposition it is usually assumed that RDX enters the vapor phase before dissociating [27-31]. These models usually assume a global reaction for the decomposition and thus are not dependent upon the details of the dissociation mechanism. Presumably the availability of such a mechanism could lead to refinement in the modeling. 

This CPU issue along with historical limitations in analytical chemistry have often forced practice-oriented reaction models to be at a global level. These models traditionally involve the relevant boiling point or solubility defined reactants, and not the controlling molecules or intermediates. Thus they will often appear to suffer from the point of view of chemical significance however they are generally much easier to solve. Most importantly, these models are relevant the reactants and products in these models are bought and sold ... 

Asphaltene and resid pyrolysis provide two relevant examples of global pyrolysis models. The pyrolysis of an isolated asphaltene feedstock typically yields the type of data summarized in Figure 2, a plot of the temporal variation of weight based product fractions as a function of time (7). This figure illustrates the exponential disappearance of asphaltene accompanied by the formation of coke, maltene and gas product fractions. Consideration of the initial slopes for the formation of coke, maltene and gas fractions led to the type of reaction network shown in Figure 3. Since resid and its reaction products can likewise be defined in terms of the solubility and volatility-based product groups asphaltene,... 

Hydrotreating reaction models are nicely exemplified by the Amoco Easy-Hard lump model of Figure 5. The global paths and the kinetics reported in this figure were taken from the review by Beaton and Bertolacini (3), which summarizes experimental results, the modelling, and commercialization of Amoco s hydrotreating unit In short, the subdivisions of resid into hard and easy fractions and gas oil into reactant and product fractions were required to account for kinetic nonlinearities. The seven lumps and 14 rate constants provided significant flexibility and were able to describe the kinetics of several relevant product fractions. 

It is worth reconsidering the reaction modelling hierarchy of Figure 1 in the light of the foregoing discussion. Clearly, the relevant global models are easily formulated and nicely adaptable to process engineering models. However, they do so at the expense of chemical fundamentals and therefore can only be used as interpretative correlations. At the other extreme, mechanistic models provide detailed chemical analyses of reaction systems. However, they bear the burden of extensive computational requirements and face the reality of the paucity of experimental data. Molecular models provide a convenient compromise between the two. The CPU requirements are reasonable, while at the same time chemical fundamentals are retained because of the implicit connection to reaction mechanisms. 

Experimentation begins with the definition of the kinetic problem to be studied and of the experimental methods used to study it. There are experimental methods for studying either quasi-isolated elementary processes or a global reaction under experimental conditions similar to those encountered in industry. A description of the most important experimental methods can be found in refs. 1—13. A discussion of some problems connected with reactor modelling is given in Sect. 3. 

TABLE II. QUASIGLOBAL TOLUENE AND ISO-OCTANE MODEL (continued). RATE PARAMETERS FOR GLOBAL REACTIONS... 

Quasiglobal kinetics models, which have previously been shown to represent lean and stoichiometric combustion of a variety of hydrocarbon fuels, have been extended to represent lean and rich combustion of toluene and iso-octane. The model predicts the thermal state of the flow and emissions of CO, soot, and N0X. The thermal state of the flow and the stable species were shown to be accurately predicted for jet-stirred combustor experiments. For rich combustion, hydrocarbon intermediates and soot are additional combustion products. The global reactions and rates were developed to represent near-adiabatic jet-stirred combustor data and were then verified by comparison to the near iso-thermal jet-stirred combustor data. N0X emissions behavior was investigated with the quasiglobal kinetics model to represent rich combustion... 

A combustion model for gaseous fuels based on the ratio of components was implemented assuming that all fluids which are mixed are also combusted [10]. The hypothesis is that mixing of reaction partners by diflusion is much slower than reaction kinetics. To calculate increasing of temperature during oxidation the global reaction... 

About 20 years ago, the theory of the be- and r-matrices, a global algebraic model of the logical structure of constitutional chemistry was formulated. This theory is the first direct mathematical approach to chemistry which also accentuates its dynamic aspect. The representation by mathematics comprises the individual objects of chemistry and also their relations, including their interconvertibility by chemical reactions. A decade later, the theory of the chemical identity groups was published in a monograph. It is a unified theory of stereochemistry that is primarily devoted to relations between molecular systems. 

The pseudo-homogeneous assumption means that both the solid and fluid phases are are considered a single phase. Therefore, we avoid considering mass and heat transfer from and to the catalytic pellets. This model assumes that the conqionent concentrations and the temperature in the pellets are the same as those in the fluid phase. This assumption is approximated when the catalyst pellet is small and mass and heat transfer between the pellets and the fluid phase are rapid. The reaction rate for this model, called the global reaction rate, includes heat and mass transfer. If heat and mass transfer are made insignificant, then the reaction rate is called the intrinsic reaction rate. 

---