Local stability analysis

An important outcome of these simulations is the location of HB points (largely ignored in previous work), which is important for the development of extinction theory. In particular, the turning point E lies on a locally unstable stationary solution branch and does not coincide with the actual extinction, as previously thought. The actual extinction point is the termination point of oscillations. Thus, local stability analysis is essential to properly analyze flame stability and develop extinction theory. 

Of the various methods of weighted residuals, the collocation method and, in particular, the orthogonal collocation technique have proved to be quite effective in the solution of complex, nonlinear problems of the type typically encountered in chemical reactors. The basic procedure was used by Stewart and Villadsen (1969) for the prediction of multiple steady states in catalyst particles, by Ferguson and Finlayson (1970) for the study of the transient heat and mass transfer in a catalyst pellet, and by McGowin and Perlmutter (1971) for local stability analysis of a nonadiabatic tubular reactor with axial mixing. Finlayson (1971, 1972, 1974) showed the importance of the orthogonal collocation technique for packed bed reactors. 

Equations (3.20) and (3.21) with their stationary-state solutions (3.24) and (3.25) are simple enough to provide a good introduction to some of the mathematical techniques which can serve us so well in analysing these sorts of chemical models. In the next sections we will explain the ideas of local stability analysis ( 3.2) and then apply them to our specific model ( 3.3). After that we introduce the basic aspects of a technique known as the Hopf bifurcation analysis ( 3.4) which enables us to locate the conditions under which oscillatory states are likely to appear. We set out only those aspects that are required within this book, without any pretence at a complete... 

We have now seen how local stability analysis can give us useful information about any given state in terms of the experimental conditions (i.e. in terms of the parameters p and ku for the present isothermal autocatalytic model). The methods are powerful and for low-dimensional systems their application is not difficult. In particular we can recognize the range of conditions over which damped oscillatory behaviour or even sustained oscillations might be observed. The Hopf bifurcation condition, in terms of the eigenvalues k2 and k2, enabled us to locate the onset or death of oscillatory behaviour. Some comments have been made about the stability and growth of the oscillations, but the details of this part of the analysis will have to wait until the next chapter. 

Many of the chemistry references appropriate to this chapter have been given in chapter 2. Local stability analysis is covered in most advanced mathematical texts on non-linear ordinary differential equations, for example ... 

Each component of the perturbations has been separated into two terms a time-dependent amplitude An and Tm, and a time-dependent spatial term cos (nnx). If the uniform state is stable, all the time-dependent coefficients will tend in time to zero. If the uniform state is temporally unstable even in the well-stirred case, but stable to spatial patterning, then the coefficients A0 and T0 will grow but the other amplitudes Ax-Ax and 7 1-7 0O will again tend to zero. If the uniform state becomes unstable to pattern formation, at least some of the higher coefficients will grow. This may all sound rather technical but is really only a generalization of the local stability analysis of chapter 3. 

Local stability analysis and Hopf bifurcation for n variables... 

The local stability analysis for a two-dimensional system has been given in some detail in chapter 3 and used frequently in later sections. The same general principles apply to any number of variables. Provided the perturbations imposed upon the stationary state are infinitesimal, the growth or decay for an n variable system will be governed by the sum of n exponential terms of the form e 1 1 with i = 1, n. 

Villadsen, J. and Michelsen, M. L. (1972) Diffusion and reaction in spherical catalyst pellets steady state and local stability analysis. Chem. Engng Sci. 27, 751. 

Local stability analysis of steady states Figure 3.29... 

Local stability analysis of the steady state is performed by linearizing equations (1) to (4) around the steady state. If we define the deviation variables as follows... 

As in the model of Section 2, the problem can be studied on its omega limit set with three rest points Eq,Ei,E2. A local stability analysis and, for some special cases, the asymptotic behavior of solutions were given in [E]. However, the populations cannot invade each other simultaneously El and E2 cannot be simultaneously unstable), so the persistence theory does not hold [E]. Moreover, for Michaelis-Menten dynamics, when one of the boundary rest points is locally stable and the other unstable, the locally stable one is globally stable [HWE]. In particular, the oscillation observed in the case of system (3.2) does not occur with (3.4). Indeed, the delayed system seems to behave much like the simple chemostat. 

THERMOKINETIC FEEDBACK OSCILLATIONS AND LOCAL STABILITY ANALYSIS... 

As we have seen, there may be more than one solution to Eq. [12]. We will consider each steady state solution in turn, that is, we will carry out a local stability analysis for each that satisfies Eq. [12]. Because the local analysis will involve a linearization of the full equations, this type of analysis is also called a linear stability analysis. 

Global StaMlity in the CSTR.— The failure of linear stability analysis to cover the macroscopic behaviour of the CSTR is well illustrated by the oscillatory states computed by Aris and Amundson for such a reactor operating with feedback control. Local stability analysis indicates an unstable equilibrium state but in the large this is surrounded by a stable limit cycle and the resultant behaviour is one of temperatures and concentrations oscillating about an unstable state, rather than approaching a stable one. 

It is then possible to perfonn a local stability analysis to characterize the convergence to equilibrium of 0 and f to the equilibrium (NakBT, 0). Although the analysis makes certain stringent assumptions they are verified in [225] for a Lennard-Jones molecular system. 

A local stability analysis can be performed on the transformation defined in Eq. (31). In the neighborhood of the fixed point the linear stability matrix is... 

Local Stability Analysis of Lumped Process Models... 

This section contains the basic notions and techniques which are used for local stability analysis of lumped process models. 

Open loop local stability analysis of DAE models... 

A homoclinic bifurcation is a composite construction. Its first stage is based on the local stability analysis for determining whether the equilibrirun state is a saddle or a saddle-focus, as well as what the first and second saddle values are, and so on. On top of that, one deals with the evolution of a -limit sets of separatrices as parameters of the system change. A special consideration should also be given to the dimension of the invariant manifolds of saddle periodic trajectories bifurcating from a homoclinic loop. It directly correlates with the ratio of the local expansion versus contraction near the saddle point, i.e. it depends on the signs of the saddle values. 

---