Numerical treatment of Gray-Scott model with operator splitting method

Abstract

This article focuses on the numerical solution of a classical, irreversible Gray Scott reaction-diffusion system describing the kinetics of a simple autocatalytic reaction in an unstirred ow reactor. A novel finite element numerical scheme based on B-spline collocation method is developed to solve this model. Before applying finite element method, ``strang splitting'' idea especially popularized for reaction-diffusion PDEs has been applied to the model. Then, using the underlying idea behind finite element approximation, the domain of integration is partitioned into subintervals which is sought as the basis for the B-spline approximate solution. Thus, the partial derivatives are transformed into a system of algebraic equations. Applicability and accuracy of this method is justified via comparison with the exact solution and calculating both the error norms \begin{document}$ L_2 $\end{document} and \begin{document}$ L_\infty $\end{document} . Numerical results arising from the simulation experiments are also presented.