Article ID Journal Published Year Pages File Type
4627106 Applied Mathematics and Computation 2015 11 Pages PDF
Abstract

Reaction diffusion operators have been used to model many engineering and biological systems. In this study we consider a reaction diffusion system modeling various engineering and life science problems. There are many algorithms to approximate such mathematical models. Most of the algorithms are conditionally stable and convergent. For a big time step size a Krylov subspace type solver for such models converges slowly or oscillates because of the presence of the diffusion term. Here we study a multigrid preconditioned generalized minimal residual method (GMRES) for such a model. We start with a five point scheme for the spatial integration and a method of lines for the temporal integration of the system of PDEs. Then we implement a multigrid iterative algorithm for the full discrete model, and show some numerical results to demonstrate the dominance of the solver. We analyze the convergence rate of such a multigrid iterative preconditioning algorithm. Reaction diffusion systems arise in many mathematical models and thus this study has many applicabilities.

Related Topics
Physical Sciences and Engineering Mathematics Applied Mathematics
Authors
,