Validation of numerical solution of diffusive part in a reaction–diffusion model

In this work, we present a methodological procedure to validate the numerical solution of the diffusive part in a reaction–diffusion model. Uniform explicit finite differences method is used to generate the solution in a confined circular domain with boundary condition of zero flux. For the validation of the numerical solution, we consider three different criteria applied to normal diffusion and sub-diffusive cases: (i) the moments of concentration, ​(ii) decay of the concentration at the origin and (iii) the mass conservation. The numerical solution fulfills the validation criteria of moments and concentration decay at the origin only in the long-term. The mass conservation criterion is fulfilled when the initial condition is imposed close to the border, whereas when it is set near to the origin a dependence on the diffusion rate appears. Pattern formation is presented after validating the numerical solution for normal diffusive case. Good agreement of stationary spatial pattern against reported results is observed.