Linearised analyses

Guha [5] pointed out some limitations in the linearised analyses developed by Horlock and Woods to determine the changes in optimum conditions with the three parameters n (and n ),/ and Not only is the accurate determination of (Cpg)i3 (and hence n ) important but also the fuel-air ratio although small, it cannot be assumed to be a constant as r is varied. Guha presented more accurate analyses of how the optimum conditions are changed with the introduction of specific heat variations with temperature and with the fuel-air ratio. 

In some cases, however, it is possible, by analysing the equations of motion, to determine the criteria by which one flow pattern becomes unstable in favor of another. The mathematical technique used most often is linearised stabiHty analysis, which starts from a known solution to the equations and then determines whether a small perturbation superimposed on this solution grows or decays as time passes. 

The types of system behaviour predicted, by the above analysis are depicted in Figs. 3.16 and 3.17. The phase-plane plots of Fig. 3.17 give the relation of the dependant variables C and T. Detained explanation of phase-plane plots is given in control textbooks (e.g., Stephanopoulos, 1984). Linearisation of the reactor model equations is used in the simulation example, HOMPOLY. 

This analysis is limited, since it is based on a steady-state criterion. The linearisation approach, outlined above, also fails in that its analysis is restricted to variations, which are very close to the steady state. While this provides excellent information on the dynamic stability, it cannot predict the actual trajectory of the reaction, once this departs from the near steady state. A full dynamic analysis is, therefore, best considered in terms of the full dynamic model equations and this is easily effected, using digital simulation. The above case of the single CSTR, with a single exothermic reaction, is covered by the simulation examples, THERMPLOT and THERM. Other simulation examples, covering aspects of stirred-tank reactor stability are COOL, OSCIL, REFRIG and STABIL. 

Most laboratories now have access to powerful computers and an extensive array of commercially available data analysis software (e.g., Prism (GraphPad, San Diego, CA), Sigma Plot (San Rafael, CA)). This provides ready access to the use of nonlinear regression techniques for the direct analysis of binding data, together with appropriate statistical analyses. However, there remains a valuable place for the manual methods, which involve linearisation, particularly in the undergraduate arena and these have been rehearsed in the above text. 

Prom the results presented in this chapter, it has been shown that the first step in the control problem of a CSTR should be the use of an appropriate mathematical model of the reactor. The analysis of the stability condition at the steady states is a previous consideration to obtain a linearised model for control purposes. The analysis of a CSTR linear model is carried out trough a scaling up reactor s volume in order to investigate the difference between the reactor and jacket equilibrium temperatures as the volume of the reactor changes from small to high value. 

Since the velocity of flow has a parabolic function, the velocity profile near the walls is nonlinear. However, in many works, this is approximated by a linearised form, as the gradient right at the walls. This makes the mathematical analysis of diffusion near one of the walls easier. Differentiating (13.1) and setting y = —h (that is, considering the bottom surface), we obtain... 

In Chap. 2, the concept of the diffusion layer was established. It is a thickness, within which a large fraction of diffusional changes take place, and at a distance of several times this thickness, practically no more diffusional changes are observed. This layer will here be given the symbol 8b (D for diffusion). In fluid dynamics, there is a similar layer, within which most of the velocity changes occur. This is the hydrodynamic layer 8. It turns out that for diffusive mass transfer, 8b is usually much smaller than <5/,. This is fortunate, because it justifies to some extent the linearised velocity profiles often assumed near walls, making analysis easier. These relations are very lucidly discussed in a classic paper by Vielstich [560]. 

Calibration is necessary for in-situ spectrometry in TLC. Either the peak height or the peak area data are measured, and used for calculation. Although the nonlinear calibration curve with an external standard method is used, however, it shows only a small deviation from linearity at small concentrations [94.95 and fulfils the requirement of routine pharmaceutical analysis 96,97J. One problem may be the saturation function of the calibration curve. Several linearisation equations have been constructed, which serve to calculate the point of determination on the basis of the calibration line and these linearisation equations are used in the software of some scanners. A more general problem is the saturation function of the calibration curve. It is a characteristic of a wide variety of adsorption-type phenomena, such as the Langmuir and the Michaelis-Menten law for enzyme kinetics as detailed in the literature [98. Saturation is also evident for the hyperbolic shape of the Kubelka-Munk equation that has to be taken into consideration when a large load is applied and has to be determined. 

As shown in Table 19, calculations [88] using the full-potential linearised augmented plane wave (FLAPW) method for the Al(lOO)—c(2 x 2)—Li phase formed by adsorption of 0.5 ML Li at room temperature con rm quantitatively the results of the LEED analysis described in Sec. 4.3 where, as shown in Fig. 15, Li was found to adsorb in a four-fold, substitutional site. 

Fig, 1. Thallium measured in the musculature of the Bream (Abramis brama) from the Saarland conurbation. Closed circles Rehlingen barrage weir open circles Gudingen barrage weir. The solid line shows the linear component of the trend as estimated with the Theil estimator and the dotted line depicts the trend estimator for data linearised prior to the analysis... 

For small disturbances of the adsorption layer from equilibrium, Lucassen (1976) derived an analytical solution (cf Section 6.1.1). An analysis of the effect of a micellar kinetics mechanism of stepwise aggregation-disintegration and the role of polydispersity of micelles was made by Dushkin Ivanov (1991) and Dushkin et al. (1991). Although it results in analytical expressions, it is based on some restricting linearisations, for example with respect to adsorption isotherm, and therefore, it is valid only for states close to equilibrium. 

An analysis of the complex problem of liquid flow and adsorption kinetics was given by Ziller et al. (1985) (cf. Section 4.8). The system of linearised hydrodynamic and transport equations was solved numerically under the assumption of different boundary conditions. The results show that a quantitative interpretation of inclined plate experiments requires a time consuming analysis. However, a rather good approximation is possible by neglecting the Marangoni flow induced by surface tension gradients against the flow direction. 

Because the mathematical theory of nonlinear dynamic systems is complex, such models are used more for dynamic simulation purposes than for the analysis and design of control systems. Simpler models can be obtained by linearisation around a nominal operating point. 

A scaled linearised state-space description around the nominal operating point of the dynamic model has been generated. The matrices A, B, C, and D have been exported in MATLAB , where a controllability analysis as fimction of fi equency has been performed. Alternatively, transfer fimctions have been generated. The results are similar. 

In subsequent treatments, Teramoto et al (25, 26) developed an analysis of both semi-batch and continuous reactor performance which incorporated a quantitative discrimination of the role of film and bulk reaction. The intractability of the non-linear product terms in the diffusion/reaction equations was ultimately avoided by a linearisation method, identical to that proposed by Hikita and Asai (22). In this approach the profiles Co(x) and Cg(x) are replaced by their interfacial values, so t t the diffusion/reaction equations become... 

Sharma et al (28) point out that the simple mass and chemical series resistance model previously used by Schmitz and Amundsen (29) is completely inadequate for describing such a system, since it is precisely the complex nature of the interactions of diffusion and reaction which give rise to the complex reactor multiplicity phenomena. Sharma et al s analysis is somewhat incomplete since they do not evaluate the film/bulk interactions, nor do they make an appraisal of the validity of the Teramoto linearisation approximation. Again, this particular form of approximation could be highly inappropriate if the concentration of reagent B falls to zero in the mass transfer film. This probably occurs at their highest temperature steady states at a /M(1) value of around 100. 

Using the data from Sect. 3.4.3, perform a more detailed analysis of the effect of linearisation on the model. Consider how linearisation changes the distribution of the original data points and how this could impact the parameter estimates obtained. What kind of transformations will cause this behaviour ... 

The numerical analysis method used for this work Is based on that described in reference (2). Full details of the assumptions made are given in reference (2), these Include Incompressible, isovisoous lubricant of negligible inertia, rigid circular journal and bearing, etc. This analysis method has been successfully applied to steadily loaded bearings, and to small journal displacement and velocity perturbations required for the computation of linearised stiffness and damping coefficients. 

It is important to note that the linearised displacement and velocity coefficients, commonly used in lateral vibration analysis, only facilitate the estimation of change of oil film force components from an equilbrium condition. In contrast with this, the oil film force components given by equations [5] and [6] are the total values. The estimation of the oil film force components at any location of the journal within the bearing... 

In the next section extensive simulations are carried out in order to determine an operating range and conditions for the use of a linearised model in control system analysis and design. These simulations cover all aspects of vehicle motion under realistic operating conditions considering driver s inputs as well as external disturbances. For that purpose transient and steady-state analysis is performed. However only a few results are present for the reason of space. Full details are given in [2]. 

---