3-Point block backward differentiation formula with an off-step point for the solutions of stiff chemical reaction problems

A major challenge in simulating chemical reaction processes is integrating the stiff systems of Ordinary Differential Equations (ODEs) describing the chemical reactions due to stiffness. Thus, it would be of interest to search systematically for stiff solvers that are close to optimal for such problems. This paper presents an implicit 3-Point Block Backward Differentiation Formula with one off-step point (3POBBDF) for the solutions of first-order stiff chemical reaction problems. In deriving the method, the Lagrange polynomial was adopted as the basis function. The paper further analyses the basic properties of the 3POBBDF which include order of accuracy, consistence, zero-stability, and convergence. The stability region as well as the interval of instability of the method was also computed. To demonstrate the accuracy of the proposed approach, some famous stiff chemical reaction problems such as Robertson problem and Chemical AKZO were solved, and the results obtained were compared with those of some existing methods. The results obtained clearly show that the 3POBBDF performs better than the existing methods with which we compared our results.