“Development of a Python program for investigation of reaction-diffusion coupling in oscillating reactions”

### Github

### Introduction

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.

### 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].## 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.

## Leave a Reply