Selectivity kinetic parameter estimation

During the selectivity kinetic parameter estimation, the relationship for x in terms of C5 - is determined from Eq. (12). For an assumed set of rate constants K, x is calculated for each composition data point such that the experimentally measured C5- equals that estimated from Eq. (12). Selectivity composition profiles as a function of C5- are generated in this manner. The proper selectivity matrix K will be that which minimizes the deviation between experimental and predicted profiles for the hydrocarbons other than C5-, as illustrated in Fig. 10. 

Obtaining Kinetic Samples for Reactive Extrusion. To develop and test kinetic models, homogeneous samples with a well defined temperature-time history are required. Temperature history does not necessarily need to be isothermal. In fact, well defined nonisothermal histories can provide very good test data for models. However, isothermal data is very desirable at the initial stages of model building to simplify both model selection and parameter estimation problems. 

The desire to formulate reaction schemes in terms of molecular processes taking place on a catalyst surface must be balanced with the need to express the reaction scheme in terms of kinetic parameters that are accessible to experimental measurement or theoretical prediction. This compromise between mechanistic detail and kinetic parameter estimation plays an important role in the use of reaction kinetics analysis to describe the reaction chemistry for a catalytic process. Here, we discuss four case studies in which different compromises are made to develop an adequate kinetic model that describes the available observations determined experimentally and/or theoretically. For convenience, we selected these examples from our work in this field however, this selection is arbitrary, and many other examples could have been chosen from the literature. 

The kinetic parameters associated with the synthesis of norbomene are determined by using the experimental data obtained at elevated temperatures and pressures. The reaction orders with respect to cyclopentadiene and ethylene are estimated to be 0.96 and 0.94, respectively. According to the simulation results, the conversion increases with both temperature and pressure but the selectivity to norbomene decreases due to the formation of DMON. Therefore, the optimal reaction conditions must be selected by considering these features. When a CSTR is used, the appropriate reaction conditions are found to be around 320°C and 1200 psig with 4 1 mole ratio of ethylene to DCPD in the feed stream. Also, it is desirable to have a Pe larger than 50 for a dispersed PFR and keep the residence time low for a PFR with recycle stream. 

Estimation of parameters. Model parameters in the selected model are then estimated. If available, some model parameters (e.g. thermodynamic properties, heat- and mass-transfer coefficient, etc.) are taken from literature. This is usually not possible for kinetic parameters. These should be estimated based on data obtained from laboratory expieriments, if possible carried out isothermal ly and not falsified by heat- and mass-transport phenomena. The methods for parameter estimation, also the kinetic parameters in complex organic systems, and for discrimination between models are discussed in more detail in Section 5.4.4. More information on parameter estimation the reader will find in review papers by Kittrell (1970), or Froment and Hosten (1981) or in the book by Froment and Bischoff (1990). 

These four procedures are all recommended to be performed in the order shown to achieve optimal parameter estimation followed by a final validation of the gravity sewer process model (Figure 7.7). In the case of design of a new sewer system, procedure number 4 is, of course, not relevant and kinetic parameters for the sewer biofilm must be evaluated and selected based on information from comparative systems. 

Figure 1. Kinetic parameters for the selection of antibody-catalyzed aldol and retro-aldol reactions, reflecting the biocatalyst s ability to accept substrates that differ clearly with respect to their molecular geometry. No background reaction was observed for the self-condensation of cyclopentanone. The indicated value for cyclopentanone addition to pentanal was estimated using the published kuncat value of 2.28 X 10 M s for the aldol addition of acetone to an aldehyde. Reproduced with permission of the authors and the American Association for the Advancement of Science.

Selected entries from Methods in Enzymology [vol, page(s)] Computer programs, 240, 312 infrared S-H stretch bands for hemoglobin A, 232, 159-160 determination of enzyme kinetic parameter, 240, 314-319 kinetics program, in finite element analysis of hemoglobin-CO reaction, 232, 523-524, 538-558 nonlinear least-squares method, 240, 3-5, 10 to oxygen equilibrium curve, 232, 559, 563 parameter estimation with Jacobians, 240, 187-191. 

The preceding equations can be used with appropriate experimental data to determine the kinetic parameters k, E, n, k , j , n i, Kw, and Kh. The pseudo-monomolecular form of Eq. (5) can be used to simplify the parameter estimation into a selectivity problem (determine K) and an activity problem (determine 0). This reduces the original highly nonlinear problem into two... 

The deactivation kinetics were determined through a series of seven separate parameter estimation problems. As with the start-of-cycle case, separate estimating problems resulted from uncoupling the reactions of each carbon number by properly selecting the charge stock. This allowed the independent determination of submatrices in the rate constant matrix Dp [Eq. (37)]. 

The microkinetic analysis is certainly a scientifically interesting approach which will contribute to the identification and selection of catalytic compounds even in more complex situations as described above. One problem still to be solved is the experimental procurement and/or estimation of the parameters used in microkinetic simulations, which limits the wide applicability of the method. Providing kinetic parameters for a complex reaction network from kinetic experiments for an analogous catalyst is a time-consuming process. Despite the availability of modem experimental equipment and efficient computers, a complex reaction demands at least one man year of work [51]. The estimation of parameters by ab initio or semiempirical methods has to be considered with caution because ideal surfaces are usually assumed. 

In this chapter the aspects of model selection/discrimination and parameter estimation and the experimental acquisition of kinetic data are not dealt with, since they fall fell outside its scope. Moreover, in interpreting the observed temperature dependency of the rate coefficients in this chapter it was assumed that we are dealing with intrinsic kinetic data. As will be shown in a Chapter 7, other, parasitic, phenomena of mass and heat transfer may interfere, disguising the intrinsic kinetics. Criteria will be presented there, however, to avoid this experimental problem. 

Figure 30.6 shows a prediction of the plasma concentration of ARA-C and total radioactivity (ARA-C plus ARA-U) following administration of two separate bolus intravenous injections of 1.2 mg/kg to a 70-kg woman. All compartment sizes and blood flow rates were estimated a -priori, and all enzyme kinetic parameters were determined from published in vitro studies. None of the parameters was selected specifically for this patient only the dose per body weight was used in the simulation. The prediction has the correct general shape and magnitude. It can be made quantitative by relatively minor changes in model parameters with no requirement to adjust the parameters describing metabolism. 

Pyrolysis kinetics of all the selected lignocellulosic wastes is properly described over the wide thermal degradation range 25°C 900°C by a model that considers an increasing dependence of the activation energy on the temperature and waste conversion with the process course. Appreciable differences in the estimated kinetic parameters are found. 

An integral part of the parameter estimation methodology is mechanism discrimination, i.e., selection of the best mechanism that would result in the best kinetic model. Nonlinear parameter estimation is an extensive topic and will not be further discussed here. For more details see Froment and Hosten, "Catalytic Kinetics —Modeling, in Catalysis—Science and Technology, Springer-Ver-lag, New York, 1981. 

For homogeneous systems, the kinetic parameters k and tq are commonly determined by lab-scale batch experiments in which all reagents are combined at the start of the reaction. Following the concentration of all components (reactants and products) over the course of the reaction then allows for the estimation of the kinetic parameters. Since water has limited solubility in the reaction mixture at the start, conventional kinetic batch experiments could result in erroneous calculation of kj and tq if the limits for homogeneity are crossed. To ensure reaction homogeneity and reliable kinetic measurements, the gradual and continuous addition of water was selected as a suitable method for experimentation (semi-batch mode). The kinetic parameters were then recovered using an appropriate mathematical model with parameter estimation module. 

Next, it was declared that the comparison of calculations with selected experiments cannot be used as a universal criterion for model development and selection of rate constants, or should be restricted and at least performed very carefully. The values of the kinetic parameters were taken from independent measurements, theoretical calculations, or semi-empirical estimations. So, the invariability of rate constants for the description of any particular experiments was accepted as the one of the main principle. 

Microsome preparations from the livers of rats gavaged with coal tar creosote were used to assay the activities of two glucuronosyltransferases, 1-hydroxypyrene UGT and />nitrophcnol UGT, and to estimate the kinetic parameters of the two enzymes (Luukanen et al. 1997). Pretreatment with creosote increased the ratio of V /K by 18-fold for 1-hydroxypyrene UGT and by 2-3-fold for/ -nitrophenol, suggesting that a highly efficient form of glucuronosyltransferase was selectively induced by creosote. 

In practice individual molecules can rarely be identified and they are instead collected into molecular classes. For reactions that occur in a particular process step, the methodology includes reaction rules that identify the structures which will undergo reaction and converts the reactant vector into that of the corresponding product(s). In addition, the appropriate kinetic parameters are selected for each reactant vector for that reactant. The physical properties of each product vector can be calculated so that product properties can be arrived at and the quality estimated as well as yield. 

Based on the obtained results, immobilized enzyme prepared at pH 4 (CALB-4) and 5 (CALB-5) were selected for thermal stability studies at 60 °C. Table 2 shows the kinetic parameters of thermal deactivation at 60 °C estimated by fitting deactivation model, Eq. 2, to experimental data. It can be observed that the deactivation profiles of CALB-4 and CALB-5 were similar, as ki and 2 are almost the same. However, by comparing the half-lives of CALB-4 and CALB-5 to CALB-7A, it can be observed that the biocatalyst prepared at pH 7 is 3- and 3.8-fold more stable than the ones prepared at pH 4 and 5, respectively. Considering that CALB has an optimum pH between 7 and 8 [35], when the enzyme is adsorbed on coconut fiber at pH 7, a favorable molecule conformation is preserved, being, therefore, more thermal stable than CALB-4 and CALB-5. Therefore, lipase immobilized at pH 7 was selected for fiirther kinetic, operational, and thermal stability studies. 

The SSF parameters are listed in Table 5 and the measured and calculated data can be seen in Fig. 8. In modeling the fermentation reaction kinetics, it was found that for the parameter estimation simulations to match Matlab Simulink simulations of the same model, the data had to be weighted. Giving higher weights to a data point means that the selected data point has more influence over the parameter estimates. If certain data points are more precisely... 

A simultaneous estimation of the kinetic parameters selected in the previous step is performed using the QN, whereas all other parameters remained fixed in the global optimum region, calculated by the RGA. The FORTRAN IMSL routine DBCONF was... 

---