Block hybrid method with off-step points for numerical solutions of classical Blasius Equation in fluid dynamics

Document Type : Original Article

Authors

1 Department of Mathematics, University of Jos, Plateau State, Nigeria.

2 Department of Statistics, Federal College of Forestry, Jos, Plateau State, Nigeria.

3 Department of Mathematics and Statistics, University of Windsor, Ontario Canada.

Abstract
Changes in conditions and parameters of equations in fluid dynamics can place a significant burden on procedures that lead to their approximate analytic solutions. While some classes of numerical methods can easily accommodate such changes, explicit linear multistep methods (LMM) have problems of instability. In this research – in order to address problem of instability associated with explicit LMM when used to solve nonlinear differential equations – a hybrid form of explicit LMM incorporating two off-step points in block form was derived using multistep collocation and matrix inversion technique. The two off-step points involved an interpolation point and a collocation point. Furthermore, the derived block method was shown to be convergent; investigation of its region of absolute stability indicated an $A(\alpha )$-stable method having large region of absolute stability with $\alpha=88.7^{0}$ and capable of handling nonlinear systems. For implementation purpose, it was tested on classical Blasius problem in fluid dynamics and results compared well with Matlab ode23s known for solving stiff and nonlinear problems.

Keywords

Subjects

[1] Abbasbandy, S. , (2007). A numerical solution of Blasius equation by Adomians decomposition method and comparison with homotopy perturbation method, Chaos, Soliton and Fractals, 31(1) , 257-260. 
[2] Ahmad, I., Bilal, M. , (2014). Numerical Solution of Blasius Equation through Neural Networks Algorithm, American J. Compu. Math. 4 , 223-232. 
[3] Ahmad, I., Bilal, M. , (2017). A novel method for the solution of Blasius equation in semi-infinite domains, An Int. J. of Opt. and Control:Theories and Apps. 7(2) , 225-233. 
[4] Akinfenwa, O.A., Jator, S.N., Yao, N.M. , (2012). On the stability of continuous block backward differentiation formula for solving stiff ordinary differential equations, J. Modern Meth. Num. Meth., 3(2) , 50-58.
[5] Al-Khaled, K.,Ajeel, M. S., Darweesh, A., & Al-Khalid, H., (2024). A reliable algorithm for solving Blasius boundary value problem, Results in Nonlinear Anal., 7(4), 1-8. 
[6] Blasius, P. H.(1908). Grenzschichten in Flussigkeiten mit kleiner Reibung. Zeitschrift fur Mathematik und Physik, 56(1).
[7] Cebeci, T., Keller, H.B., (1971). Shooting and parallel shooting methods for solving the Falkner-Skan boundary layer equation, J. Comput. Phy., 7(2) , 289-300. 
[8] Celaya, E. A., Anza, J. J. (2013). BDF-α: A multistep method with numerical damping control. Universal Journal of Computational Mathematics, 1(3), 96-108. 
[9] Chu, M.T., Hamilton, H. , (1987). Parallel solution of ODEs by multi-block methods, Siam J. Sci. Stat. Comput, 8(3) , 342-353.
[10] Cortell, R., (2005). Numerical solutions of the classical Blasius flat-plate problem, Appl. Math. Comput. 170 , 706-710. 
[11] Fatunla, S.O., (1995). A Class of Block Methods for Second Order IVPs, International J. Comput. Math. 55 , 119-133.
[12] Fatunla, S.O., (1991). Block Methods for Second Order ODEs, Int. J. Compu. Math. 41 , 55-63. 
[13] Ganji, D.D., Babazadeh, H., Noori, F., Pirouz, M.M., Janipour, M., (2009). An application of homotopy perturbation method for nonlinear blasius equation to boundary layer flow over a flat plate, International Journal of Nonlinear Science. 7 , 399-404.
[14] He, J.H., (1997). A new approach to nonlinear partial differential equations, Comm. Nonlinear Sci. Nume. Simulation. 2 , 230-235. 
[15] Henrici, P., (1962). Discrete variable methods in ODEs, John Wiley, New York. 
[16] Howarth, L., (1938). On the solution of the laminar boundary layer equations, Proc. London Math. Soc. 2 , 547-579.
[17] Ibrahim, Z.B., Othman, K.I., Suleiman, M., (2007). Implicit r-point block backward differentiation formula for solving first-order stiff ODEs, Appl. Math. Comput. 186 , 558-565. 
[18] , M.K., Molla, M.R., Sultana, S., (2011). Numerical approximations of Blasius boundary layer equation, Dhaka University J. Sci. 59 , 87-90. 
[19] Jator, S.N., Sahi, R., Akinyemi, M., Nyonna, D., (2021). Exponentially fitted block backward differentiation formulas for pricing options, Cogent Economics & Fin. 9(1), 1875565. 
[20] Liao, S.J., (1999). An explicit, totally analytic approximate solution for Blasius viscous flow problems, Int. J. Non-Linear Mech. 34 , 759-778.
[21] Liu, Y., Kurra, S.N., (2012). Solving Blasius equation using HVIM, World Acad. Sci. 65, 578-581.
[22] Miansari, M.O., Miansari, M.E., Barari, A., Domairry, G., (2010). Analysis of Blasius Equation for Flat-plate Flow with Infinite Boundary Value, Int. J. Comput. Meth. Engnr. Sci. & Mech. 11(2) , 79-84.

[23] Okuonghae, R.I., Ikhile, M.N.O., (2010). Stiffly stable second derivative linear multistep methods with two hybrid points, Num. Analysis & Appl. 8, 248-259.
[24] Okuonghae, R.I., Ikhile, M.N.O., (2013). A class of hybrid linear multistep methods with A(α)-stability properties for stiff IVPs in ODEs, J. Num. Math. 21(2), 157-172, doi: 10.1515/jnum-2013-0006.
[25] Okuonghae, R. I., & Ozobokeme, J. K., (2024). Falkner Hybrid Block Methods for Second-Order IVPs: A Novel Approach to Enhancing Accuracy and Stability Properties, J. Numer. Anal. Approx. Theory, 53(2), 324-342. 
[26] Onumanyi, P., Sirisena, U.W., Jator, S.N., (1999). Continuous finite difference approximations for solving differential equations, Int. J. Comput. Math. 72(1) , 15-27.
[27] Parand, K., Taghavi, A., (2009). Rational scaled generalized Laguerre function collocation method for solving the Blasius equation, J. Comput. & Appl. Math. 233(4) , 980-989. 
[28] Rahman, M.M., (2018). The Comparison between analytical and numerical solution of Blasius equation on boundary layer parallel flow over a flat plate, J. Sci. Engg. Res. 6(2), 1-12. 
[29] Rahman, M.M., Khan, S., Akbar, M.A., (2018). Numerical and analytical solutions of new Blasius equation for turbulent flow, Heliyon. 9(2), e14319, 1-13.
[30] Selamat, M.S., Halmi, N.A., Ayob, N.A., (2019). A semi analytic iterative method for solving two forms of Blasius equation, J. Academia. 7(2) , 76-85. 
[31] Shateyi, S., (2023). On the application of the block hybrid methods to solve linear and non-linear first order differential equations, Axioms. 12(2), 189. 
[32] Suleiman, M., Musa, H., Ismail, F., Senu, N., (2013). A new variable step size block backward differentiation formula for solving stiff initial value problems, Int. J. Comput. Math. 90, 2391-2408. 
[33] Voss, D., Abbas, S., (1997). Block predictorcorrector schemes for the parallel solution of ODEs, Comput. & Math. Appl. 33(6) , 65-72. 
[34] Wang, L., (2004). A new algorithm for solving classical Blasius equation, Appl. Math. & Comput. 157(1) , 1-9. 
[35] Wang, L., (1991). Exact solutions of the steady-state Navier-Stokes equations, Ann. Review. Fluid Mech. 23(1),159-177. 
[36] Weyl, H., (1942). On the differential equations of the simplest boundary-layer problems, Annal. Math. 43 , 381-407.
[37] White, F.M., (1991). Viscous Fluid Flow. McGraw-Hill, New York. 
[38] Xu, D., Xu, J., Xie, G., (2014). Revisiting Blasius flow by fixed point method, Abstract & Appl. Analysis, Article ID 953151, 1-9.
[39] Yasuri, A.K., (2023). An analytical self-consistent method for differential forms of the Blasius equation, Math. Meth. Appl. Sci. 46, 5836-5849. 
[40] Yao, B., Chen, J., (2009). A new analytical solution branch for the Blasius equation with a shrinking sheet, Appl. Math. & Comput. 215(3) , 1146-1153


Articles in Press, Accepted Manuscript
Available Online from 29 May 2026

  • Receive Date 07 March 2025
  • Revise Date 26 May 2026
  • Accept Date 29 May 2026