Efficient Embedded Diagonal Implicit Runge-Kutta Method for Directly Solving Third Order ODEs

Abstract

This paper presents two important contributions to the field of numerical analysis for third-order ordinary differential equations (ODEs). First, a new class of direct implicit Runge-Kutta (RK) processes, called RKTDIO, is introduced as solutions to third-order ODEs. Secondly, it develops the ERKTDIO method, which is an embedded pairwise diagonal implicit RK method. The study begins by introducing the theory of relevant-colored trees and B-series as fundamental concepts. By utilizing the order constraints, two RKTDIO methods are derived: a fifth-order method with three stages and a sixth-order method with four stages. In addition, an embedded method called ERKTDIO6(5) is derived, which has orders six and five. The derivation of the embedded method includes strategies to ensure that the higher-order method achieves high accuracy while the lower-order method provides optimal error estimates. To evaluate the effectiveness of the proposed methods, variable step-size codes are developed and applied to a set of specific third-order problems. The numerical evaluation involves converting the problems into a system of first-order ODEs and comparing the results with existing methods in terms of accuracy and function evaluations. The numerical demonstrations emphasise the superior performance and efficiency of the new methods in solving third-order ODEs. The comparative analysis shows the accuracy achieved by the higher-order method and the improved error estimation of the lower-order method. The results validate the efficacy of the proposed approaches and their potential for practical applications in various domains.

Article Details

How to Cite
[1]
Fawzi, F.A. et al. 2024. Efficient Embedded Diagonal Implicit Runge-Kutta Method for Directly Solving Third Order ODEs. Ibn AL-Haitham Journal For Pure and Applied Sciences. 37, 2 (Apr. 2024), 388–408. DOI:https://doi.org/10.30526/37.2.3409.
Issue
Section
Mathematics

How to Cite

[1]
Fawzi, F.A. et al. 2024. Efficient Embedded Diagonal Implicit Runge-Kutta Method for Directly Solving Third Order ODEs. Ibn AL-Haitham Journal For Pure and Applied Sciences. 37, 2 (Apr. 2024), 388–408. DOI:https://doi.org/10.30526/37.2.3409.

References

Boutayeb, A.; Chetouani A. A mini-review of numerical methods for high-order problems. International Journal of Computer Mathematics. 2007, 84(4), 563-579, https://doi.org/10.1080/00207160701242250.

Wu, Xiong-Jian; Yigong Wang; W. G. Price. Multiple resonances, responses, and parametric instabilities in offshore structures. Journal of ship research. 1988, 32(04), 285-296, https://doi.org/10.5957/jsr.1988.32.4.285

Cortell, R. Application of the fourth-order Runge-Kutta method for the solution of high-order general initial value problems. Computers & structures. 1993, 49(5), 897-900, https://doi.org/10.1016/0045-7949(93)90036-D

Malek, A.; R. Shekari Beidokhti. Numerical solution for high order differential equations using a hybrid neural network-optimization method. Applied Mathematics and Computation. 2006, 183(1), 260-271, https://doi.org/10.1016/j.amc.2006.05.068.

Awoyemi, D. O. Algorithmic collocation approach for direct solution of fourth-order initial-value problems of ordinary differential equations. International Journal of Computer Mathematics. 2005, 82(3), 321-329, https://doi.org/10.1080/00207160412331296634.

Adel, F.; Senu, N.; Ismail, F.; Majid, Z. A. New efficient phase-fitted and amplification-fitted Runge-Kutta method for oscillatory problems. International J Pure Appl Math, 2016, 107(1), 69-86, https://doi.org/10.12732/ijpam.v107i1.7

Hussain, K.; Ismail, F., ; Senu, N. A new optimized Runge-Kutta-Nyström method to solve oscillation problems. World Applied Sciences Journal, 2015, 33(10), 1614-1622, https://doi.org/10.5829/idosi.wasj.2015.33.10.300.

Hussain, K. A.; Abdulnaby, Z. E.. A new two derivative FSAL Runge-Kutta method of order five in four stages. Order, 2020, 17(1). https://doi.org/10.21123/bsj.2020.17.1.0166

Lee, K. C.; Senu, N., Ahmadian, A., Ibrahim, S. N. I., & Baleanu, D. Numerical study of third-order ordinary differential equations using a new class of two derivative Runge-Kutta type methods. Alexandria Engineering Journal, 2020, 59(4), 2449-2467, https://doi.org/10.1016/j.aej.2020.03.008.

Ahmad, N. A.; Senu, N.; Ismail, F. Phase-fitted and amplification-fitted higher order two-derivative Runge-Kutta method for the numerical solution of orbital and related periodical IVPs. Mathematical Problems in Engineering, 2017. 10.1109/ICREM.2015.7357020.

Senu, N.; Ahmad, N. A., ; zarina, B. I. B. I. Numerical Study on Phase-Fitted and Amplification-Fitted Diagonally Implicit Two Derivative Runge-Kutta Method for Periodic IVPs. Sains Malaysiana, 2021, 50(6), 1799-1814, 10.17576/jsm-2021-5006-25.

Fawzi, F. A. An Embedded 5 (4) Pair of Optimized Runge-Kutta Method for the Numerical Solution of Periodic Initial Value Problems. Iraqi Journal of Science, 2022, 3889-3901, https://doi.org/10.24996/ijs.2022.63.9.21.

SENU, N.; AHMAD, N.; ISMAIL, F. Embedded 5 (4) Pair Trigonometrically-Fitted Two Derivative Runge-Kutta Method with FSAL Property for Numerical Solution of Oscillatory Problems, WSEAS TRANSACTIONS on COMPUTERS, 2017,16, 155-162.

Jikantoro, Y. D.; Ismail, F., Senu, N., ; Ibrahim, Z. B., Hybrid methods for direct integration of special third order ordinary differential equations. Applied Mathematics and Computation, 2018, 320, 452-463, https://doi.org/10.1016/j.amc.2017.10.003.

Shokri, A.; Mehdizadeh Khalsaraei, M. A new efficient high order four-step multiderivative method for the numerical solution of second-order IVPs with oscillating solutions. Mathematics Interdisciplinary Research, 2020, 5(2), 157-172, 10.22052/mir.2020.211603.1185.

Mechee, M.; Senu, N.; Ismail, F.; Nikouravan, B.; Siri, Z. A three-stage fifth-order Runge-Kutta method for directly solving special third-order differential equation with application to thin film flow problem. Mathematical Problems in Engineering. 2013, 2013, https://doi.org/10.1155/2013/795397.

Fawzi, F. A.; Jumaa M, H. The Implementations Special Third-Order Ordinary Differential Equations (ODE) for 5th-order 3rd-stage Diagonally Implicit Type Runge-Kutta Method (DITRKM). Ibn AL-Haitham Journal For Pure and Applied Sciences. 2022, 35(1), 92-101, https://doi.org/10.30526/35.1.2803.

ISMAIL F.; Al-Khassawneh, R. A. Solving Delay Differential Equations Using embedded singly diagonally Implicit runge-kutta methods. Acta Mathematica Vietnamica. 2008, 33(2), 95-105.

Jikantoro, Y. D., Ismail, F., Senu, N., & Ibrahim, Z. B,. A new integrator for special third order differential equations with application to thin film flow problem. Indian Journal of Pure and Applied Mathematics, 2018, 49, 151-167, 10.1007/s13226-018-0259-6.

Lee, K. C.; Senu, ; Ahmadian, A. ; Ibrahim, S. N. I., On two-derivative Runge–Kutta type methods for solving u‴= f (x, u (x)) with application to thin film flow problem. Symmetry, 2020, 12(6), 924, https://doi.org/10.3390/sym12060924.

Hussain, K.; Ismail, F.; Senu, N. Solving directly special fourth-order ordinary differential equations using Runge–Kutta type method. Journal of Computational and Applied Mathematics. 2016, 306, 179-199, https://doi.org/10.1016/j.cam.2016.04.002.

Ghawadri, N.; Senu, N.; Adel Fawzi, F.; Ismail, F.; Ibrahim, Z. B. Explicit Integrator of Runge-Kutta Type for Direct Solution of u^((4))= f (x,u,u',u^''). Symmetry. 2019, 11(2), 246, https://doi.org/10.3390/sym11020246.

Dormand, J. R. Numerical methods for differential equations: a computational approach. CRC press, 1996, https://doi.org/10.1201/9781351075107 .

Fawzi, F. A.; Jaleel, N. W. Exponentially Fitted Diagonally Implicit EDITRK Method for Solving ODEs. Ibn AL-Haitham Journal For Pure and Applied Sciences. 2023, 36(1), 389-399, https://doi.org/10.30526/36.1.2883 .

Saleh, O. M.; Fawzi, F. A.; Hussain, K. A.; Ghawadri, N. G. Derivation of Embedded Diagonally Implicit Methods for Directly Solving Fourth-order ODEs. Ibn AL-Haitham Journal For Pure and Applied Sciences, 2024, 37(1), 375-385, https://doi.org/10.30526/37.1.3288.

Senu, N.; Suleiman, M.; Ismail, F. An embedded explicit Runge–Kutta–Nyström method for solving oscillatory problems. Physica Scripta. 2009, 80(1), 015005. https://doi.org/10.1088/0031-8949/80/01/015005

Hairer E.; Hochbruck, M.; Iserles, A.; Lubich, C. Geometric numerical integration; Oberwolfach Reports. 2006, 3(1), 805-882. https://doi.org/10.1017/S0962492902000144.

Hairer E, Norsett SP, Wanner G. Solving Ordinary Differential Equations I: Nonstiff Problems, II: Stiff Problems. Springer, 1987. https://doi.org/10.1007/978-3-540-78862-1.

Lambert, J. D. Numerical methods for ordinary differential systems. New York: Wiley, 1991. https://doi.org/10.1002/9781119121534