This article proposes an implicit operator splitting nonstandard finite difference (OS-NSFD) scheme for numerical treatment of two species in three space dimensions reaction-diffusion glycolysis model. Since, the unknown state variables exhibiting the concentrations of species in glycolysis models and they cannot be negative and obtaining their positive solutions is a challenging task. The established theoretical result ensures that our proposed OS-NSFD scheme is unconditionally convergent at equilibrium point and fulfills the condition of positivity of solutions on contrary to other methods. Further, we analyze the existence and uniqueness of the solution obtained for the underlying system. To highlight the effectiveness of OS-NSFD scheme we compare the simulation results of OS-NSFD scheme with three well-known existing operators splitting finite difference (FD) schemes, namely, forward Euler explicit, backward Euler implicit and Crank Nicolson splitting schemes. Many existing techniques provide with the restricted positive solutions which do not work always. These techniques are only applicable if certain conditions on the discretized parameters are considered otherwise; they produce negative solutions, which is not the physical feature of the real system. The current work bridges this gap by catering the unconditional positive solutions to the reaction diffusion models.
Read full abstract