Chemistry Reference
In-Depth Information
with
aN .N 1 C S /
2.N 1 C ˇ S /
S .2.1 ˇ S / N/
2.N 1 C ˇ S /
A 11 D 1
and A 12 D
:
The characteristic polynomial of this matrix is the quadratic
2
.A 11 C ˇ S
C NA 12 / C ˇ S A 11 ;
so the conditions for asymptotic stability are (see Appendix F)
.1 ˇ S /.1 A 11 / NA 12 >0;
(4.48)
.1 C ˇ S /.1 C A 11 / C NA 12 >0;
(4.49)
ˇ S A 11 <1;
(4.50)
which reduce to the conditions,
.1 ˇ S / aN .N 1 C S /
aNˇ S .2.1 ˇ S / N/
2.N 1 C ˇ S /
2.N 1 C ˇ S /
>0; (4.51)
.1 C ˇ S / 2
aN .N 1 C S /
2.N 1 C ˇ S /
aNˇ S .2.1 ˇ S / N/
2.N 1 C ˇ S /
C
>0; (4.52)
aNˇ S .N 1 C S /
2.N 1 C ˇ S /
S
1:
(4.53)
It is obvious that due to the algebraic complexity of these stability conditions
further analytical calculations will become quite involved. Therefore, instead we
give a brief numerical study that will give us a flavor of the results one might expect
to hold in general.
First of all we investigate the effect of the bifurcation that marks the exchange of
the equilibrium E 3 with the equilibrium E 1 occurring at
A.1 ˇ S /.N 1 C ˇ S /
cN 2
L bif
D
.1/ . This is not a usual
transcritical bifurcation because E 1 and E 3 are fixed points of two different maps,
and the merging occurs along a line of non-differentiability. Indeed, this is a typical
border collision bifurcation, the effect of which is quite difficult to forecast. This is
shown by three different bifurcation diagrams obtained for increasing values of the
capacity limit L across the bifurcation value. The first bifurcation diagram, shown in
Fig. 4.13 is obtained with the set of parameters N D 4, A D 2, c D 0:15, ˇ S D 0:6,
and a D 0:5, with L in the range Œ0:5;2. For this set of parameters the bifurcation
value is L D 1:2 and, as it can be seen in Fig. 4.13, when L crosses the bifurcation
.3/
along the boundary that separates the regions
D
and
D
Search WWH ::




Custom Search