Analysis of Convergence and Stability Properties of Diagonally Implicit 3-Point Block Backward Differentiation Formula for First Order Stiff Initial Value Problems

Authors

  • Buhari Alhassan Department of Mathematics and Statistics, College of Natural and Applied Sciences, Al-Qalam University, Katsina, Nigeria https://orcid.org/0000-0003-0784-0160
  • Hamza Yusuf Department of Mathematics, Faculty of Physical Sciences, Federal University, Dutsin-Ma, Katsina State, Nigeria
  • Abdulrahman Adamu Department of Mathematics and Statistics, Umaru Musa Yar’adua University, Katsina, Nigeria

DOI:

https://doi.org/10.56919/usci.2432.021

Keywords:

Order, Consistency, Zero-stability, A-stability, Absolute stability region, Convergence, Block Backward Differentiation Formula

Abstract

Study’s Excerpt/Novelty

  • This paper introduces a novel diagonally implicit 3-point block backward differentiation formula (BDF) for efficiently solving first-order stiff initial value problems.
  • The method, proven to have an order of accuracy of 5, meets the necessary and sufficient conditions for convergence, including consistency and zero stability.
  • Comparative numerical results highlight its superior performance in terms of maximum error and CPU time, demonstrating its advantage over existing methods and offering a robust solution for integrating stiff initial value problems.

Full Abstract

This paper comprehensively analyses the diagonally implicit 3-point block backward differentiation formula (BDF) for solving first-order stiff initial value problems.  We establish the necessary and sufficient conditions for convergence, including consistency and zero stability, and derive the method's order of accuracy, which is found to be 5.  Stability analysis reveals that the method is almost A-stable, with an absolute stability region plotted.  A C programming language code is developed using Newton's Iteration for numerical implementation and compiled in the Microsoft Dev C++ compiler environment.  Comparative numerical results demonstrate the superior performance of the proposed method over the existing fully implicit 3-point block backward differentiation formula (3BBDF) in terms of maximum error and CPU time.  Therefore, this method offers a new and efficient numerical solution for integrating stiff initial value problems.

References

Abasi N., Suleiman M. B., Abassi N. and Musa H. (2014). 2-point block BDF with off-step points for solving stiff ODEs. Journal of soft computing and Applications, 39, 1-15. https://doi.org/10.5899/2014/jsca-00039

Alhassan B., Musa H, and Abasi N. (2022). Convergence and Order of the 2-point Diagonally Implicit Block Backward Differentiation Formula with Two Off-step Points. UMYU Scientifica, 1(1), 30-38. https://doi.org/10.56919/usci.1222.004

Alhassan B., Yusuf H., Musa H. & Abasi N. (2023a). A New Fifth Order Variable Step Size Block Backward Differentiation Formula with Off-Step Points for the Numerical of Stiff Initial Ordinary Differential Equations. Applied Mathematics and Computational Intelligence, 12, 94-121. https://doi.org/10.58915/amci.v12i4.235

Alhassan B. and Musa H. (2023b). Diagonally Implicit Extended 2-Point Super Class of Block Backward Differentiation Formula for solving First Order Stiff Initial Value Problems. Applied Mathematics and Computational Intelligence, 12 (1), 101-124.

Aminikhah H. and Hemmatnezhad M. (2011): An Effective Modification of the Homotopy Perturbation Method for Stiff Systems of Ordinary Differential Equations. Applied Mathematics Letters, 24(9), 1502-1508. https://doi.org/10.1016/j.aml.2011.03.032

Butcher, C. J. (2016): Numerical Methods for Ordinary Differential Equations. Third Edition,John Wiley & Sons, Ltd. https://doi.org/10.1002/9781119121534

Bala N., Alhassan B. and Musa H. (2022): A new 3-point implicit block numerical method for solving first order stiff initial value problems. Timbou-Africa Academic Publications, International Journal of Science Research and Technology, 9 (9), ISSN: 2623-7861.

Curtis C.F. and Hirschfelder J.O. (1952): Integration of stiff Equations, Proceedings of the National Academy of Sciences of the United States of America, 38(3), 235-243. https://doi.org/10.1073/pnas.38.3.235

Dahlquist G. (1974). Problems related to the numerical treatment of stiff differential equations. International Computing symposium, 307-314.

Haziza M. J. and Ibrahim Z. B. (2019): Diagonally Implicit Block Backward Differentiation Formula with Optimal Stability Properties for Stiff Ordinary Differential Equations, Symmetry, 1-18. https://doi.org/10.3390/sym11111342:

Ibrahim, Z., B, Suleiman, M., B., Johari, R. and Ismail, F. (2003): On the stability of fully implicit block backward differentiation formulae. Jobatan Matematik, UTM, 83-89.

Ibrahim Z. B., Othman K. I. and Suleiman M. B. (2007): Implicit r-point block backward differentiation formula for first order stiff ODEs, Applied Mathematics and Computation, 186 (1), 558-565. https://doi.org/10.1016/j.amc.2006.07.116

Ibrahim Z. B. and Nasaruddin A. A. (2020): A class of hybrid multistep block methods with A-Stability for the Numerical Solution of Stiff Ordinary Differential equations. MDPI mathematics, 914(8), 1-20. https://doi.org/10.3390/math8060914:

Lambert, J. D. (1973): Computational methods in Ordinary Differential Equations. New York: John wiley& Sons.

Lambert J. D. (1991): Numerical methods for ordinary differential systems. New York: John wiley & Sons.

Musa H., Suleiman M. B., Ismail F., Senu N., Majid Z. A. and Ibrahim Z. B. (2014): A new fifth order implicit block method for solving first order stiff ODEs. Malaysian journal of Mathematical sciences 8(8), 45-59.

Musa H. and Bala N. (2019): 3-point diagonally implicit super class of block backward differentiation formula for solving stiff initial value problems. Dutse Journal of Pure and Applied Sciences (DUJOPAS), 5 (1b), 1-10. https://doi.org/10.59568/JASIC-2022-3-2-01

Musa H., Alhassan B. &Abasi N. (2022). Diagonally Implicit 2-point Block Backward Differentiation Formula with Two Off-Step Points for Solving Stiff Initial Value Problems. Nigerian Journal of Mathematics and Applications. 32 (II):108-119. https://doi.org/10.56919/usci.1222.004

Nasir N. A. A. M., Ibrahim Z. B. and Suleiman, M. B. (2011): Fifth order two-point block backward differentiation formulas for solving ordinary differential equations. Applied Mathematical Sciences. 5(71), 3505-3518.

Noor, N. M., Yatim, S. A. M. and Ibrahim, Z. B. (2024).Fractional Block Method for the Solution of Fractional Order Differential Equations. Malaysian Journal of Mathematical Sciences, 18(1), 185-205. https://dio.org/10.47836/mjms.18.1.11.

Suleiman M.B., Musa H., Ismail F. and Senu N. (2013): A new variable step size block backward differentiation formula for solving stiff initial value problems. International Journal of Computer Mathematics, https://doi.org/10.1080/00207160.2013.776677

Suleiman M. B., Musa H. and Ismail F. (2015): An implicit 2-point block extended backward differentiation formula for integration of stiff initial value problems. Malaysian journal mathematical sciences, 9(1), 33-51.

Xin-Yuan., W. (1998): A Sixth-Order A-Stable Explicit One-Step Method for Stiff Systems. Computers Mathematics Application, 35 (9), 59-64. https://doi.org/10.1016/S0898-1221(98)00057-1

Yaakub A. R. and Evans D. J. (2003): New L-Stable Modified Trapezoidal Methods for the Iinitial Value Problems. International Journal of Computer Mathematics, 80(1), 95-104. https://doi.org/10.1080/00207160304663

Yusuf H., Musa H. and Alhassan B. (2024): A New Fixed Coefficient Diagonally Implicit Block Backward Differentiation Formula for Solving Stiff Initial Value Problems. UMYU Scientifica, 3(1), 1-14. https://doi.org/10.56919/usci.2431.001

Zawawi I.S.M, Ibrahim Z.B., Ismail F. and Majid Z.A. (2012); Diagonally implicit block backward differentiation formulas for solving ODEs, Hindawi Publishing Coorperation. International journal of mathematics and mathematical sciences. 767328. https://doi.org/10.1155/2012/767328

Published

2024-06-28

How to Cite

Alhassan, B., Yusuf, H., & Adamu, A. (2024). Analysis of Convergence and Stability Properties of Diagonally Implicit 3-Point Block Backward Differentiation Formula for First Order Stiff Initial Value Problems. UMYU Scientifica, 3(2). https://doi.org/10.56919/usci.2432.021