Modeling chemical constraints

Incorporation of Chemical Constraints. To complete the model of the oxide-electrolyte interface, it is necessary to return to the preceding section on models of chemical reactions to find values of ctq, CT-p and 172 calculated from chemical considerations. 

Bond valences can be used in conjunction with other techniques, particularly powder diffraction where, for example, light atoms are difficult to refine in the presence of heavy atoms. Adding the chemical constraints of the bond valence model can stabilize the refinement, particularly in the case of superstructures that have high pseudo-symmetry (Thompson et al. 1999). 

The most recent progress in chemical constraints refers to the implementation of a physicochemical model into the resolution process [64, 66-73], In this manner, the concentration profiles of compounds involved in a kinetic or a thermodynamic process are shaped according to the suitable chemical law (see Figure 11.8). A detailed description of methods for fitting kinetic models to multivariate data is provided in Chapter 7. 

Figure 5.3.9 (A) Simplified geometric model [46, 89] for the preparation of industrial Cu/ZnO catalysts comprising subsequent meso- and nanostructuring of the material from [56], In a first micro structure directing step (mesostructuring), the Cu,Zn coprecipitate crystallizes in the form of thin needles of the zincian malachite precursor, (Cu,Zn)2(0H)C03. In a second step, the individual needles are decomposed and demix into CuO and ZnO. The effectiveness of this nanostructuring step depends critically on a high Zn content in the precursor, which in zincian malachite is limited to Cu Zn ca. 70 30 due to solid-state chemical constraints [75]. Finally, interdispersed CuO/ZnO is reduced to yield active Cu/ZnO. (B) Chemical memory Dependence of catalytic activity in methanol synthesis on the conditions of the coprecipitation and aging steps, from [85].

Many approaches already do this, for example, by incorporating known chemical constraints, densities or hard-sphere repulsions. Many of the emerging methods described below have this flavor, and as time goes on our ability to complex our data and our modeling approaches will only increase. 

The software also evolved as my group caught the interactive modeling bug. I converted the batch program GT to REACT, which was fully interactive. The user entered the chemical constraints for his problem and then typed go to trigger the calculation. Ming-Kuo Lee and I added Pitzer s activity model and a... 

Fig. 6. A section of our electron density map at an early stage of the partial structure Fourier procedure. P is for phosphate and R for ribose. The skeletal model drawing in incomplete for R o and Rvi- Electron densities for Rqq, R70, and C70 are out of the section chosen. Peo, Re9, Ue9, and P70 were not included in the calculation of the phases but came up in the electron density map. After the preliminary fitting as shown here, the refinement procedure optimizes the fitting to the target positions and the chemical constraints.

Constraint definition Here the steric and/or chemical constraints of the design problem are delineated and supplied to the program in an appropriate form. In many cases these constraints will be derived from the active site (or, more generally, receptor site) derived from homology models, receptor models, - pharmacophore models, CoMFA models, or just a single molecule may also be used. 

The optimization formulation (presented in Eqs. I.II-I.I6) consists of the objective function (e.g., minimize TAG Eq. l.II) subjected to process constraints, the process models and constraints (Eqs. I.I-I.IO) of the generic model block mentioned earlier (x is a process variable, the mass flow rate), structural constraints (Eqs. 1.12 and 1.13) representing the superstructure which allows selection of only one process alternative in each step, and cost functions (Eqs. 1.14-1.16) to calculate the operating and capital costs using cost parameters (f l", waste treatment cost utility or chemicals cost , reactor... 

This map has a single quadratic extremum, similar to tliat of tire WR model described in detail earlier. Such maps (togetlier witli tire technical constraint of negative Schwarzian derivative) [23] possess universal properties. In particular, tire universal (U) sequence in which tire periodic orbits appear [24] was observed in tire BZ reaction in accord witli tliis picture of tire chemical dynamics. 

Many additional consistency tests can be derived from phase equiUbrium constraints. From thermodynamics, the activity coefficient is known to be the fundamental basis of many properties and parameters of engineering interest. Therefore, data for such quantities as Henry s constant, octanol—water partition coefficient, aqueous solubiUty, and solubiUty of water in chemicals are related to solution activity coefficients and other properties through fundamental equiUbrium relationships (10,23,24). Accurate, consistent data should be expected to satisfy these and other thermodynamic requirements. Furthermore, equiUbrium models may permit a missing property value to be calculated from those values that are known (2). 

Once the objective and the constraints have been set, a mathematical model of the process can be subjected to a search strategy to find the optimum. Simple calculus is adequate for some problems, or Lagrange multipliers can be used for constrained extrema. When a Rill plant simulation can be made, various alternatives can be put through the computer. Such an operation is called jlowsheeting. A chapter is devoted to this topic by Edgar and Himmelblau Optimization of Chemical Processes, McGraw-HiU, 1988) where they list a number of commercially available software packages for this purpose, one of the first of which was Flowtran. 

No single method or algorithm of optimization exists that can be apphed efficiently to all problems. The method chosen for any particular case will depend primarily on (I) the character of the objective function, (2) the nature of the constraints, and (3) the number of independent and dependent variables. Table 8-6 summarizes the six general steps for the analysis and solution of optimization problems (Edgar and Himmelblau, Optimization of Chemical Processes, McGraw-HiU, New York, 1988). You do not have to follow the cited order exac tly, but vou should cover all of the steps eventually. Shortcuts in the procedure are allowable, and the easy steps can be performed first. Steps I, 2, and 3 deal with the mathematical definition of the problem ideutificatiou of variables and specification of the objective function and statement of the constraints. If the process to be optimized is very complex, it may be necessaiy to reformulate the problem so that it can be solved with reasonable effort. Later in this section, we discuss the development of mathematical models for the process and the objec tive function (the economic model). 

The second classification is the physical model. Examples are the rigorous modiiles found in chemical-process simulators. In sequential modular simulators, distillation and kinetic reactors are two important examples. Compared to relational models, physical models purport to represent the ac tual material, energy, equilibrium, and rate processes present in the unit. They rarely, however, include any equipment constraints as part of the model. Despite their complexity, adjustable parameters oearing some relation to theoiy (e.g., tray efficiency) are required such that the output is properly related to the input and specifications. These modds provide more accurate predictions of output based on input and specifications. However, the interactions between the model parameters and database parameters compromise the relationships between input and output. The nonlinearities of equipment performance are not included and, consequently, significant extrapolations result in large errors. Despite their greater complexity, they should be considered to be approximate as well. 

An example adapted from Verneuil, et al. (Verneuil, V.S., P. Yan, and F. Madron, Banish Bad Plant Data, Chemical Engineeiing Progress, October 1992, 45-51) shows the impact of flow measurement error on misinterpretation of the unit operation. The success in interpreting and ultimately improving unit performance depends upon the uncertainty in the measurements. In Fig. 30-14, the materi balance constraint would indicate that S3 = —7, which is unrealistic. However, accounting for the uncertainties in both Si and S9 shows that the value for S3 is —7 28. Without considering uncertainties in the measurements, analysts might conclude that the flows or model contain bias (systematic) error. 

---