Developing Higher-Order Unconditionally Positive Finite Difference Methods for the Advection Diffusion Reaction Equations

This study introduces the higher-order unconditionally positive finite difference (HUPFD) methods to solve the linear, nonlinear, and system of advection-diffusion-reaction (ADR) equations. The stability and consistency of the developed methods are analyzed, which are necessary and sufficient for the numerical approach to converge to the exact solution. The problem under consideration is of the Cauchy type, and hence, Von Neumann stability analysis is used to analyze the stability of the proposed schemes. The HUPFD's efficacy and efficiency are investigated by calculating the error, convergence rate, and computing time. For validation purposes, the higher-order unconditionally positive finite difference solutions are compared to analytical calculations. The numerical results demonstrate that the proposed methods produce accurate solutions to solve the advection diffusion reaction equations. The results also show that increasing the order of the unconditionally positive finite difference leads an implicit scheme that is conditionally stable and has a higher order of accuracy with respect to time and space.