I'm a passionate theoretical scientist who finds joy in developing theories and algorithms that try to capture the complex physical nature of chemistry.
The Brusselator is a model developed by I. Prigogine and R. Lefever at the Universite Libre de Bruxelles in Belgium. [16] developed for the description of nonlinear systems. The model includes some assumptions [12, 11], such as irreversible reactions ( for and constant reactant and product species concentrations A, B, C and D. Realistically, this happens by constant addition of A and B and removal of C and D.
With the assumptions of chapter 2.1, the dierential kinetic time laws of the species X and Y can be derived. The concentration of a species is noted with square brackets. The concentration notation is assumed to be dimensionless in the following, since there is a mathematical consideration of the issue. With and for the differential time laws follow.
Spatially Homogeneous Brusselator
Dynamic Analysis of the Spatially Homogeneous Brusselator
This part is rather long and thus shortened. The full version is shown in detail in my thesis. The fixpoint can be calculated.
It is not possibile to analyse non-linear differential equations. Because of this, the System needs to be linearised by using an approximation. The simplified Idea behind my procedure was using a disturbance and investigating its behaviour.
Next both functions and are substituted and a Taylor series approximation is used.
Same for .
By neglecting any higher therms () it is obvious, that the matrix of partial derivatives is the Jacobian. This is only possible as long as the position and type of the fixed point of the linearized system is also exactly the one of the non-linear system.
Eigenvalues can easily be calculated and simplified using trace and determinant.
Now, the stability and the type of the fixpoint can be calculated.
For the spatially homogeneous Brusselator there are four domains depending on [A] and [B]. The domains are shown in Tab. 2.1 and Fig. 2.6.
For the spatially homogeneous Brusselator there are four domains depending on [A] and [B]. The domains are shown in Tab. 2.1 and Fig. 2.6.
The system behavior is independent of starting conditions and as long as and holds. For and the system would already be in equilibrium.
Figure 2.6: [A]-[B] phase space plot with marked regions providing stable fixed points and unstable fixed points as well as with the marked region resulting in complex eigenvalues.
Limit Cycle
With the use of Poincaré–Bendixson theorem, it can be shown, that a limit cycle is possibile. This depends on the domain present. The proof for the Brusselator model is carried out in [25, p. 42] or more generally in [26, p. 196].
Figure 2.7: Poincare-Bendixon theorem schematically applied to the phase space plot Chap. 5.1 with constraints [A] = 1 and [B] = 3.
Reaction-Diffusion System of Brusselator
Both differential equations can now be advanced by using the second Fick’s law.
There is no analytical solution for the RDS of Brusselator.
Dynamic Analysis of the Reaction-Diffusion System
It is found, that the Fixpoint equals the one of the spartially homogenous system.
The linearisation works the same way as with the spartially homogenous system, just a little more complicated.
A disturbance as follows is beeing used.
With the use of , and
and
will be the result.
With the use of , and
can be formed. Because all reaction terms are spratially independent and without any spartial operators the vector functin can be aproximated as the Jacobian.
This can be shortened.
The Stability of can now be aproximated with the use of
The domains of the spartially homogenous Brusselator are now extended. Four different functions can be found, wich sepreate the domains. This time, two of them refer to the –-phase space.
Abbildung 2.9: Dx-Dy-Plot mit markierten Bereichen für und , welche stabile Fixpunkte und instabile Fixpunkte liefern sowie die Turing-StabilitätFigure 2.10: A-B phase space plot with marked regions for , , and , which provide stable fixed points and unstable fixed points, and Turing stability.