Numerical modeling of three dimensional Brusselator reaction diffusion system

In many mathematical models, positivity is one of the attributes that must be possessed by the continuous systems. For instance, the unknown quantities in the Brusselator reaction-diffusion model represent the concentration of two reactant species. The negative values of concentration produced by any numerical methods is meaningless. This work is concerned with the investigation of a novel unconditionally positivity preserving finite difference (FD) scheme to be used for the solution of three dimensional Brusselator reaction-diffusion system. Von Neumann stability method and Taylor series expansion is applied to verify unconditional stability and consistency of the proposed FD scheme. Results are compared against well-known forward Euler FD scheme and some results reported in the literature.In many mathematical models, positivity is one of the attributes that must be possessed by the continuous systems. For instance, the unknown quantities in the Brusselator reaction-diffusion model represent the concentration of two reactant species. The negative values of concentration produced by any numerical methods is meaningless. This work is concerned with the investigation of a novel unconditionally positivity preserving finite difference (FD) scheme to be used for the solution of three dimensional Brusselator reaction-diffusion system. Von Neumann stability method and Taylor series expansion is applied to verify unconditional stability and consistency of the proposed FD scheme. Results are compared against well-known forward Euler FD scheme and some results reported in the literature.

[1]  Faiz Faizullah.  On the pth moment estimates of solutions to stochastic functional differential equations in the G-framework , 2016, SpringerPlus.

[2]  M. O. Ahmad,et al.  Numerical Modeling of SEIR Measles Dynamics with Diffusion , 2018 .

[3]  Hristo V. Kojouharov,et al.  An unconditionally positivity preserving scheme for advection-diffusion reaction equations , 2013, Math. Comput. Model..

[4]  Existence of Solutions for G-SFDEs with Cauchy-Maruyama Approximation Scheme , 2014 .

[5]  Faiz Faizullah.  A Note on the Carathéodory Approximation Scheme for Stochastic Differential Equations under G-Brownian Motion , 2012 .

[6]  R. C. Mittal,et al.  A modified cubic B-spline differential quadrature method for three-dimensional non-linear diffusion equations , 2017 .

[7]  Faiz Faizullah,et al.  On existence and approximate solutions for stochastic differential equations in the framework of G-Brownian motion , 2017 .

[8]  Edward H. Twizell,et al.  Chaos-free numerical solutions of reaction-diffusion equations , 1990, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[9]  Faiz Faizullah,et al.  Existence and stability of solutions to non-linear neutral stochastic functional differential equations in the framework of G-Brownian motion , 2017, Advances in Difference Equations.

[10]  Abdul-Qayyum M. Khaliq,et al.  The locally extrapolated exponential time differencing LOD scheme for multidimensional reaction-diffusion systems , 2015, J. Comput. Appl. Math..

[11]  Faiz Faizullah.  A Note on pth Moment Estimates for Stochastic Functional Differential Equations in the Framework of G-Brownian Motion , 2017 .