General Modelling Procedure

One of the more important features of modelling is the frequent need to reassess both the basic theory (physical model), and the mathematical equations, representing the physical model, (mathematical model), in order to achieve agreement, between the model prediction and actual plant performance (experimental data). 

As shown in Fig. 1.2, the following stages in the modelling procedure can be identified  

Most reactor operations involve many different variables (reactant and product concentrations, temperature, rates of reactant consumption, product formation and heat production) and many vary as a function of time (batch, semi-batch operation). For these reasons the mathematical model will often consist of many differential equations. 

It should be noted that often the model does not have to give an exact fit to data as sometimes it may be sufficient to simply have a qualitative agreement with the process. 

Kapur (1988) has listed thirty-six characteristics or principles of mathematical modelling. Mostly a matter of common sense, it is very important to have them restated, as it is often very easy to lose sight of the principles during the active involvement of modelling. They can be 

The mathematical model can only be an approximation of real-life processes, which are often extremely complex and often only partially understood. Thus models are themselves neither good nor bad but should satisfy a previously well-defined aim. 

Modelling is a process of continuous development, in which it is generally advisable to start off with the simplest conceptual representation of the process and to build in more and more complexities, as the model develops. Starting off with the process in its most complex form often leads to confusion. 

Modelling is an art but also a very important learning process. In addition to a mastery of the relevant theory, considerable insight into the actual functioning of the process is required. One of the most important 

Models must be both realistic and robust. A model predicting effects, which are quite contrary to common sense or to normal experience, is unlikely to be met with confidence. 

---

In this introductory chapter the use of mathematical models in chemical engineering is motivated and examples are given. The general modeling procedure is described, and some important tools that are covered in greater detail later in the book are outlined. 

When more in-depth analysis of environmental fate is required, the analyst must select the modeling procedure that is most appropriate to the circumstances. In general, the more sophisticated models are more data, time, and resource intensive. 

The target level procedure was applied to 16 common air contaminants (Table 6.19). These are common contaminants in the industrial environment, and in many cases are the most critical compounds from the viewpoint of need for control measures. The prevailing concentration data as well as the benchmark levels were taken from Nordic databases, mainly the Finnish sources, and described elsewhere.In addition, a general model for assessing target values for other contaminants is presented in the table. 

Data were subjected to analysis of variance and regression analysis using the general linear model procedure of the Statistical Analysis System (40). Means were compared using Waller-Duncan procedure with a K ratio of 100. Polynomial equations were best fitted to the data based on significance level of the terms of the equations and values. 

The general objective of all radar detection procedures is to get a constant false alarm rate (CFAR) due to the fact that the test cell almost always contains clutter and noise and only in a very few cases contains radar target echo signals. The statistical model and general detection procedure, in which the detector is fixed only with regard to the noise and clutter statistic and independently to the target statistic, has been developed by Neyman and Pearson. 

This example treats a diffusion-reaction process in a spherical biocatalyst bead. The original problem stems from a model of oxygen diffusion and reaction in clumps of animal cells by Keller (1991), but the modelling method also applies to bioflocs and biofilms, which are subject to potential oxygen limitation. Of course, the modelling procedure can also be applied generally to problems in heterogeneous catalysis. 

The sequence of meal consumption was determined by random assignment of diets to subjects. Statistical analysis was performed by a General Linear Models Procedure (20) using split-plot in time analysis with the following non-orthogonal contrasts ... 

Values followed by asterisk are significantly different from the control at the 0.05 level according to the general linear model procedure. 

In this example, 30 calibration factors have to be quantified by the calibration procedure, which requires 30 independent equations and, therefore, 10 calibration samples (gives 3 equations per sample). Martens warned against pitfalls by using MLR. There is, for example, no check whether the general model is adequate to describe the unknown samples. The calculated concentrations may be in error because the sample contains other impurities than present in the calibration samples, or because some... 

A significant step in the numerical solution of packed bed reactor models was taken with the introduction of the method of orthogonal collocation to this class of problems (Finlayson, 1971). Although Finlayson showed the method to be much faster and more accurate than that based on finite differences and to be easily applicable to two-dimensional models with both radial temperature and concentration gradients, the finite difference technique remained the generally accepted procedure for packed bed reactor model solution until about 1977, when the analysis by Jutan et al. (1977) of a complex butane hydrogenolysis reactor demonstrated the real potential of the collocation procedure. 

The present procedure is a specific example of a synthetic method of some generality. The procedure describes an example which is of considerable interest per se but, perhaps more importantly, which also serves as a model for the use of this procedure for the preparation of other a-amino ketones. In the submitters laboratory, this specific procedure is used routinely for the training of persons who will be using this general technique or related techniques.29... 

Analysis of variance appropriate for a crossover design on the pharmacokinetic parameters using the general linear models procedures of SAS or an equivalent program should be performed, with examination of period, sequence and treatment effects. The 90% confidence intervals for the estimates of the difference between the test and reference least squares means for the pharmacokinetic parameters (AUCo-t, AUCo-inf, Cmax should be calculated, using the two one-sided t-test procedure). 

When no Kga or HTU data are available, their values may be estimated by means of a generalized model. A summary of useful models is given in Sec. 5. The values obtained may then be combined by use of Eq. (14-19) to obtain values of H0g and H0L. This simple procedure is not valid when the rate of absorption is limited by chemical reaction. 

Results were analyzed by nested mixed-model ANOVA s using general linear procedures, in the MINITAB 15 statistical program. Nested mixed-model ANOVA was used when multiple leaves per tree and multiple trees per treatment were available. Additional analyses were linear and quadratic regressions (performed in MINITAB 15 and Excel), and when significant differences occurred, means were compared using Student s t-test or nested mixed-model ANOVA. 

In this chapter, we describe the results of our studies we aimed at the development of a general computational procedure to generate automatically and unbiased objective pharmacophore models using the GRID approach and starting with PDB macromolecular complexes. Within the context of structure-based pharmacophore modeling, it represents an approach that is somehow complementary to that described in Chapter 6. We have used logically combined maps... 

Those based on strictly empirical descriptions Mathematical models based on physical and cnemical laws (e.g., mass and energy balances, thermodynamics, chemical reaction kinetics) are frequently employed in optimization applications. These models are conceptually attractive because a general model for any system size can be developed before the system is constructed. On the other hand, an empirical model can be devised that simply correlates input/output data without any physiochemical analysis of the process. For these models, optimization is often used to fit a model to process data, using a procedure called parameter estimation. The well-known least squares curve-fitting procedure is based on optimization theory, assuming that the model parameters are contained linearly in the model. One example is the yield matrix, where the percentage yield of each product in a unit operation is estimated for each feed component... 

The computational procedures for the numerical calculations of the TBI model for specific systems were described in several previous publications (Tan and Chen, 2005, 2008a). Here, we summarize the general computational procedure using a 12-nt DNA duplex as an example. We assume temperature T = 25 °C, dielectric constant e = 78 for solvent and 20 for the nucleic acid interior, and the radii of the (hydrated) ions 4.5 A for Mg2"1" and 3.5 A for Na+. The numerical calculations of the TBI model involve the following five steps ... 

Of course, some general aspects of our treatment could be easily extended to a general form of f b ireJ as in the semi-infinite case [226],but for explicit numerical work a specific form of fs(b ire) ((()) is needed. Equation (10) can be justified for Ising-type lattice models near the critical point [216,220], i.e. when ( ) is near ( >crit=l/2, as well as in the limits f]>—>0 or <()—>1 [11]. The linear term —pj( ) is expected due to the preferential attraction of component B to the walls, and to missing neighbors for the pairwise interactions near the walls while the quadratic term can be attributed to changes in the pairwise interactions near the walls [144,216,227]. We consider Eq. (10) only as a convenient model assumption to illustrate the general theoretical procedures - there is clear evidence that Eq. (10) is not accurate for real polymer mixtures [74,81,82,85]. 

Several vapor pressure osmometers are now commercially available. Although they are mainly used for determining number-average molecular weights in aqueous and organic solvents, they can also be employed to evaluate the total osmolality of biological solutions or dissociation and activity coefficients. Each model has its own technical characteristics. However, all are comparable in terms of general measurement procedure and sensitivity. 

Several authors have addressed the determination of the optical properties of aqueous titanium dioxide suspensions in the context of photoreactor modeling (Brandi et al., 1999 Cabrera et al., 1996 Cured et al., 2002 Salaices et al., 2001, 2002 Satuf et al., 2005 Yokota et al., 1999). Among the determined properties are extinction, scattering, and absorption coefficients, as well as the asymmetry parameter of the scattering phase function. In general the procedures involve fitting of a radiative transfer model to the experimental results for reflectance and transmittance of radiation. 

Regression Analysis. The GLM (General Linear Models) procedure of SAS was used to fit the experimental data to Equation 3. This procedure provides estimates of coefficients and intercept GLM also tests hypotheses and indicates the overall quality of the correlation. Output from the GLM procedure is shown in Tables IV, V, and VI numbers, which are listed to 6 decimal places in the original output, have been rounded off to If- places. 

---