Computational analysis of glycolytic reaction in open spatial reactor

We consider the computational analysis of processes within the spatially-distributed model simulating the glycolytic reaction taking place in the one-side fed open chemical reactor. The main point of the simulation is the decomposition of the reaction–diffusion system into unidirectional reaction in a bulk supplied by feedback terms stated as boundary conditions on the lower boundary of the reactor, i.e. the unique plane where an exchange with an outer medium is possible within the real experimental situation. Analysis of the curvature of the reagents distribution curves proves kinematic character of the observed lateral waves corresponding to the picture of experimentally observed glycolytic traveling waves. At the same time, their origin relates to diffusion of the reagents in a vertical cross-section of the reactor. Study of the solutions for the concerned reaction–diffusion model in the case of stochastically different diffusion coefficients reveals the Turing structures.