# How can I convert a Brusselator type model of reaction-diffusion to include vortical motion?

A system called Model G, which is similar to the Brusselator, describes reaction-diffusion processes taking place among hypothetical sub-quantum units. We are investigating whether it would be possible to add vortical motion to this system?

The proposed sub-quantum units differ from quarks in that they exhibit only standard classical properties such as reaction and diffusion. As a result of their processes, quantum structures are found to emerge which have properties and behaviors similar to those exhibited by subatomic particles.

This would be to mathematically describe the reaction-diffusion system in terms of fluid dynamics which includes vortices. One idea is to include Navier-Stokes equations or something similar for the convection process, i.e. reaction-diffusion-convection (non-linear) or vortical convection of the hypothetical sub-particle reactants. How would one go about this?

The equations for the Brusselator type system (Model G) in their present form, without vortical motion, are shown below. The figure shows the Brusselator type Model G, as a schematic. The 2nd image shows the five kinetic equations, followed by their PDEs.

$$$$\frac{\partial G(x,y,z,t)}{\partial t}=D_G \nabla^2 G - (k_{-1} + k_2)G + k_{-2} X + k_1 A \quad \textrm{(4a)}$$$$

$$$$\frac{\partial X(x,y,z,t)}{\partial t}=D_X \nabla^2 X + k_2 G - (k_{-2} + k_3 B + k_5) X + k_{-3}ZY - k_{-4}X^3 + k_4 X^2Y +k_{-5} \Omega \quad \textrm{(4b)}$$$$

$$$$\frac{\partial Y(x,y,z,t)}{\partial t}=D_Y \nabla^2 Y + k_3 B X - k_{-3} Z Y + k_{-4} X^3 - k_4 X^2 Y \quad \textrm{(4c)}$$$$

Letter symbols A and B denote the concentrations of the initial particle reactants; G, X, and Y denote the concentrations of the intermediate reactants; and Z and Ω denote the concentrations of the final reaction products. $$D_i$$ is a constant of proportionality called the diffusion coefficient.

The forward reaction rate constants, $$k_{i}$$ , above each arrow, and the reverse reaction rate constants, $$k_{-i}$$ , below each arrow specify the rate at which reactants (back of arrow) transform into products (front of arrow).