Numerical modelling of dynamical systems in isothermal chemical reactions and morphogenesis

Abstract

Mathematical models of isothermal chemical systems in reactor problems and Turing's theory of morphogenesis with an application in sea-shell patterning are studied. The reaction-diffusion systems describing these models are solved numerically. First- and second-order difference schemes are developed, which are economical and reliable in comparison to classical numerical methods. The linearization process decouples the reaction-diffusion equations thereby allowing the use of different time steps for each differential equation, which may be large due to the excellent stability properties of the methods. The methods avoid having to solve a non-linear algebraic system at each time step. The schemes are suitable for implementation on a parallel machine.