Next Article in Journal
Synergy Prediction Model of Information Entropy Based on Zone Safety Degree and Stope Roof Weighting Step Analysis
Next Article in Special Issue
Best Constant in Ulam Stability for the Third Order Linear Differential Operator with Constant Coefficients
Previous Article in Journal
Ratio-Type Estimator for Estimating the Neutrosophic Population Mean in Simple Random Sampling under Intuitionistic Fuzzy Cost Function
Previous Article in Special Issue
New Results for the Investigation of the Asymptotic Behavior of Solutions of Nonlinear Perturbed Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies

by
Ezekiel Olaoluwa Omole
1,
Emmanuel Oluseye Adeyefa
1,
Victoria Iyadunni Ayodele
2,
Ali Shokri
3 and
Yuanheng Wang
4,*
1
Department of Mathematics, Federal University Oye-Ekiti, P.M.B. 373, Oye-Ekiti 370112, Ekiti State, Nigeria
2
Department of Computer Science and Mathematics, Nigeria Police Academy, Wudil-Kano 713101, Kano State, Nigeria
3
Department of Mathematics, Faculty of Sciences, University of Maragheh, Maragheh 83111-55181, Iran
4
College of Mathematics and Computer Science, Zhejiang Normal University, Jinhua 321004, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(9), 891; https://doi.org/10.3390/axioms12090891
Submission received: 19 July 2023 / Revised: 12 September 2023 / Accepted: 12 September 2023 / Published: 18 September 2023
(This article belongs to the Special Issue Differential Equations and Related Topics)

Abstract

:
A computational approach with the aid of the Linear Multistep Method (LMM) for the numerical solution of differential equations with initial value problems or boundary conditions has appeared several times in the literature due to its good accuracy and stability properties. The major objective of this article is to extend a multistep approach for the numerical solution of the Partial Differential Equation (PDE) originating from fluid mechanics in a two-dimensional space with initial and boundary conditions, as a result of the importance and utility of the models of partial differential equations in applications, particularly in physical phenomena, such as in convection-diffusion models, and fluid flow problems. Thus, a multistep collocation formula, which is based on orthogonal polynomials is proposed. Ninth-order Multistep Collocation Formulas (NMCFs) were formulated through the principle of interpolation and collocation processes. The theoretical analysis of the NMCFs reveals that they have algebraic order nine, are zero-stable, consistent, and, thus, convergent. The implementation strategies of the NMCFs are comprehensively discussed. Some numerical test problems were presented to evaluate the efficacy and applicability of the proposed formulas. Comparisons with other methods were also presented to demonstrate the new formulas’ productivity. Finally, figures were presented to illustrate the behavior of the numerical examples.

1. Introduction

This paper focuses on the numerical solution of a second-order partial differential equation arising from the fluid dynamics. It describes how convection effects interact with diffusion transports. Thus, we are interested in the convection diffusion reaction equation, which is given below:
a 2 u x 2 + b 2 u y 2 + p u x + q u t + r u y + k u = ν
subject to the given initial conditions:
s 1 ( x ) = u ( x , 0 ) , x [ a , b ] , s 2 ( x ) = u y ( x , 0 ) , x [ c , d ]
and the following boundary conditions:
t 1 ( x ) = u ( a , y ) , y 0 t 2 ( x ) = u ( b , t ) , t 0 t 3 ( x ) = u ( b , y ) , y 0 t 4 ( x ) = u ( b , t ) , t 0
or any of the conditions below:
v 1 ( x ) = u ( a , t ) , t 0 v 2 ( x ) = u ( b , y ) , y 0 v 3 ( x ) = u ( a , t ) , t 0 v 4 ( x ) = u ( b , y ) , y 0 .
Symbolically, in ( 1 ) , u ( x , y ) symbolizes the dependent variable. On the other hand, x, y, and t stand for the diffusion function; likewise, ν is the non-homogeneous function, also known as the prescribed source, and a 1 , b 1 , p 1 , q 1 , and r 1 are the constants. Model Equation ( 1 ) has numerous applications in physics, mathematics, and engineering (Liu et al. [1]. According to Kamran et al. [2], it can be used to describe many real-life physical systems, among which are the Langmuir wave packet estimation in plasma, electrostatics, the non-relativistic limit of the Klein–Gordon equation, fluid flow, and bimolecular dynamics (Kilic and Celik [3], Akbarov et al. [4], Chen et al. [5], Singh et al. [6], Mirzaee et al. [7], Khan et al. [8], and Aliev et al. [9]). These classes of partial differential equations lack analytic solutions and are very difficult to approximate due to the nonlinear parameters. Accordingly, the numerical method is an alternative approach for approximating the PDEs’ solution, which is being proposed (Adeyefa et al. [10]). Because of the aforementioned applications, various numerical techniques for solving various application problems, such as time–frequency analysis, signal delay, convection–diffusion equations, nonlinear approximation, and Monte Carlo simulation, arising in fluid dynamics problems have been developed. For instance, predictor–corrector techniques (Su et al. [11], Awoyemi and Idowu [12], Iskandarov and Komartsova [13], Ashry et al. [14], Asif [15]); Galerkin methods (Guo et al. [16]); Haar wavelets methods (Aziz and Khan [17], Shiralashetti et al. [18], Saparova et al. [19]); Runge–Kutta methods (Takei and Iwata [20], Yakubu et al. [21], Zhao and Huang [22]); Newtral network model (Mall and Chakraverty [23]); multigrid technique (Ghaffar et al. [24], Ge [25], Gupta et al. [26]); finite difference methods (Mulla et al. [27]); and finite element methods (Harari and Hughes [28]).
In spite of the success recorded by the above-mentioned techniques, there is still a need for improvement in convergence and accuracy due to the fact that the majority of the methods discussed above are not self-starting and do not use second-derivative values in their derivations. A number of authors have proposed numerical self-starting block algorithms for solving differential equations with either boundary or initial conditions. Sunday et al. [29] constructed a class of hybrid algorithm block methods with variable steps and applied the developed method to solve application problems, namely the Kepler problem. Ramos and Vigo-Aguiar in reference [30] presented a BDF-type approach that is L-stable in nature for the numerical approximation of stiff problems through the method of line. Meanwhile, Ngwane and Jator [31] considered a numerical solution of an oscillatory second-order initial and Hamiltonian system of equations using a trigonometrically fitted block method. In reference [32], Modebei et al. proposed a class of numerical solvers for approximating a fourth-order partial differential equation via a block approach with a uniform order. In the same vein, Jator in [33] developed a block-like unification scheme for the numerical solution of a class of elliptic, telegraph, and sine-Gordon partial differential equations with an emphasis on accuracy. Likewise, Olaiya et al. in [34] presented a numerical approach for simulating the Black–Scholes partial differential equation via a two-step off-grip block of algorithms of algebraic order seven. Lastly, the work of Familua et al. in [35] considered a higher-order block technique for the numerical simulation of third-order boundary value problems with applications. The theoretical analysis of the methods was investigated and discussed comprehensively. It has been established in the literature that the block approach was first initiated by Milne [36]. Runge–Kutta’s self-starting nature, high-convergence rates, and unique property of yielding approximations of the solution at various points are all still present in the block approach. The convection–diffusion equation and the Helmholtz equation, to name a couple, are examples of nonlinear partial differential equations that arise in the sciences and engineering. In light of this, ninth-order multi-step collocation formulas (NMCFs) are proposed with the potential to be a useful tool for solving these equations. This study aims to develop a method that solves nonlinear second-order PDEs more precisely and with higher convergence than other methods already described in the literature.

2. Design of NMCFS

To develop a numerical method for the solution of ( 1 ) , the orthogonal polynomial of the form:
U x = r = 0 p + q 1 ζ r Ω r x ,
is adopted as the basis function. ζ r are the coefficients to be determined and Ω r is defined on the interval ( , ) , with the help of the recurrence formula:
Ω n + 1 x = x Ω n ( x ) Ω n ( x ) .
The polynomials are orthonormal with respect to the weight function e x 2 .
The first five sets of orthogonal polynomials (Hermite polynomials), as contained in Adeyefa et al. [10], are as follows:
Ω 0 ( x ) = 1 Ω 1 ( x ) = x Ω 2 ( x ) = x 2 1 Ω 3 ( x ) = x 3 3 x Ω 4 ( x ) = x 4 6 x 2 + 3 .
On the other hand, we differentiate ( 5 ) twice to obtain:
U x = r = 0 p + q 1 ζ r Ω r x .
Set p = 2 , q = 8 , with k as the step number. It follows that p is the number of distinct interpolation points, which must coincide with the order of the d.e in ( 1 ) , meanwhile q denotes the number of collocation points chosen. Thus, ( 5 ) and ( 6 ) reduces to
U x = r = 0 10 ζ r Ω r x ,
U x = r = 0 10 ζ r Ω r x .
Now, interpolate the approximate solution to ( 7 ) at x n + w , w = 0 ( 1 ) and also collocate ( 8 ) at x n + w , w = 0 ( 1 ) 8 , which produces seven non-singular equations, which can be written as a system in matrix form as
D J = L
where
D = Ω 0 ( x n ) Ω 1 ( x n ) Ω 2 ( x n ) Ω 3 ( x n ) Ω k + 2 ( x n ) Ω 0 ( x n + 1 ) Ω 1 ( x n + 1 ) Ω 2 ( x n + 1 ) Ω 0 ( x n + 1 ) Ω k + 2 ( x n + 1 ) Ω 0 ( x n ) Ω 1 ( x n ) Ω 2 ( x n ) Ω 3 ( x n ) Ω k + 2 ( x n ) Ω 0 ( x n + 1 ) Ω 1 ( x n + 1 ) Ω 2 ( x n + 1 ) Ω 3 ( x n + 1 ) Ω k + 2 ( x n + 1 ) Ω 0 ( x n + k ) Ω 1 ( x n + k ) Ω 2 ( x n + k ) Ω 3 ( x n + k ) Ω k + 2 ( x n + k )
J = ζ 0 , ζ 1 , ζ 2 , ζ 3 , , ζ 7 T , L = u m , n , u m + 1 , n , g m , n , g m + 1 , n , , g m + 8 , n T .
Using computer software, like Maple, where J = D 1 L , one can solve the matrix Equation ( 9 ) for the unknown values of ζ i , i = 0 ( 1 ) 10 . The values obtained are subsequently plugged into the Formula ( 7 ) with the setting x = ϕ h + x n + 7 to obtain the continuous function of the form
u m + j , n ϕ = Ψ 0 u m , n + Ψ 1 u m + 1 , n + h 2 j = 0 8 Δ j ϕ g m + j , n , j = 0 1 8 ,
where the coefficients of the continuous function ( 10 ) are given in (Appendix A).
The main formulas are produced by evaluating ( 10 ) at ϕ = 5 , 4 , 3 , 2 , 1 , 0, and 1, which gives the following discrete formulas.
u m + 2 , n = u m , n + 2 u m + 1 , n + 33953 518400 h 2 g m , n + 424759 453600 h 2 g m + 1 , n 81629 453600 h 2 g m + 2 , n + 11143 28350 h 2 g m + 3 , n 27533 72576 h 2 g m + 4 , n + 110563 453600 h 2 g m + 5 , n 23017 226800 h 2 g m + 6 , n + 5627 226800 h 2 g m + 7 , n 9829 3628800 h 2 g m + 8 , n ,
u m + 3 , n = 2 u m , n + 3 u m + 1 , n + 155171 1209600 h 2 g m , n + 9421 4800 h 2 g m + 1 , n + 144847 302400 h 2 g m + 2 , n + 252101 302400 h 2 g m + 3 , n 5701 8064 h 2 g m + 4 , n + 135901 302400 h 2 g m + 5 , n 56473 302400 h 2 g m + 6 , n + 4601 100800 h 2 g m + 7 , n 6029 1209600 h 2 g m + 8 , n ,
u m + 4 , n = 3 u m , n + 4 u m + 1 , n + 12869 67200 h 2 g m , n + 225469 75600 h 2 g m + 1 , n + 94001 75600 h 2 g m + 2 , n + 3271 1575 h 2 g m + 3 , n 11279 12096 h 2 g m + 4 , n + 49313 75600 h 2 g m + 5 , n 3449 12600 h 2 g m + 6 , n + 2537 37800 h 2 f m + 7 , n 4439 604800 h 2 g m + 8 , n ,
u m + 5 , n = 4 u m , n + 5 u m + 1 , n + 18481 72576 h 2 g m , n + 363173 90720 h 2 g m + 1 , n + 181289 90720 h 2 g m + 2 , n + 311363 90720 h 2 g m + 3 , n 13261 36288 h 2 g m + 4 , n + 17551 18144 h 2 g m + 5 , n 33583 90720 h 2 g m + 6 , n + 1163 12960 h 2 g m + 7 , n 3547 362880 h 2 g m + 8 , n ,
u m + 6 , n = 5 u m , n + 6 u m + 1 , n + 76859 241920 h 2 g m , n + 16883 3360 h 2 g m + 1 , n + 83207 30240 h 2 g m + 2 , n + 9043 1890 h 2 g m + 3 , n + 2449 8064 h 2 g m + 4 , n + 63047 30240 h 2 g m + 5 , n 5461 15120 h 2 g m + 6 , n + 533 5040 h 2 g m + 7 , n 407 34560 h 2 g m + 8 , n ,
u m + 7 , n = 6 u m , n + 7 u m + 1 , n + 7319 19200 h 2 g m , n + 261023 43200 h 2 g m + 1 , n + 152107 43200 h 2 g m + 2 , n + 87827 14400 h 2 g m + 3 , n + 3541 3456 h 2 g m + 4 , n + 140401 43200 h 2 g m + 5 , n + 7009 14400 h 2 g m + 6 , n + 9143 43200 h 2 f m + 7 , n 2849 172800 h 2 g m + 8 , n ,
u m + 8 , n = 7 u m , n + 8 u m + 1 , n + 57281 129600 h 2 g m , n + 114769 16200 h 2 g m + 1 , n + 67861 16200 h 2 g m + 2 , n + 15506 2025 h 2 g m + 3 , n + 3541 2592 h 2 g m + 4 , n + 77893 16200 h 2 g m + 5 , n + 9353 8100 h 2 g m + 6 , n + 10157 8100 h 2 f m + 7 , n + 5741 129600 h 2 g m + 8 , n ,
The first derivative of ( 10 ) is given below:
u m + j , n t = Ψ 0 u m , n + Ψ 1 u m + 1 , n + h 2 j = 0 8 Δ j ϕ f m + j , n , j = 0 1 8 .
The coefficients of the continuous function ( 18 ) and the first derivative of ( 10 ) are given in (Appendix B).
The additional discrete formulas are generated by evaluating ( 18 ) at the points ϕ = 7 , 6 , 5 , 4 , 3 , 2 , 1 , 0, and 1, which gives the following first derivative discrete formulas:
u m , n = 1 7257600 h ( 1624505 h 2 g m , n + 4124232 h 2 g m + 1 , n 5225624 h 2 g m + 2 , n + 6488192 h 2 g m + 3 , n 5888310 g m + 4 , n h 2 + 3698920 g m + 5 , n h 2 1522672 h 2 g m + 6 , n + 369744 h 2 g m + 7 , n 40187 h 2 g m + 8 , n + 7257600 u m , n 7257600 u m + 1 , n ) ,
u m + 1 , n = 1 7257600 h ( 515529 h 2 g m , n + 4809956 h 2 g m + 1 , n 3983564 h 2 g m + 2 , n + 4702524 h 2 g m + 3 , n 4177930 g m + 4 , n h 2 + 2593756 g m + 5 , n h 2 1059756 h 2 g m + 6 , n + 256004 h 2 g m + 7 , n 27719 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 2 , n = 1 7257600 h ( 447623 h 2 g m , n + 7561144 h 2 g m + 1 , n + 2506008 h 2 g m + 2 , n + 1197440 h 2 g m + 3 , n 1543370 g m + 4 , n h 2 + 1083672 g m + 5 , n h 2 471184 h 2 g m + 6 , n + 118192 h 2 g m + 7 , n 13125 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 3 , n = 1 7257600 h ( 462217 h 2 g m , n + 7361892 h 2 g m + 1 , n + 5782580 h 2 g m + 2 , n + 6461116 h 2 g m + 3 , n 3209610 g m + 4 , n h 2 + 1879388 g m + 5 , n h 2 755372 h 2 g m + 6 , n + 181380 h 2 g m + 7 , n 19591 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 4 , n = 1 7257600 h ( 455751 h 2 g m , n + 7434680 h 2 g m + 1 , n + 5350552 h 2 g m + 2 , n + 10280832 h 2 g m + 3 , n + 1239350 g m + 4 , n h 2 + 1027864 g m + 5 , n h 2 502800 h 2 g m + 6 , n + 129968 h 2 g m + 7 , n 14597 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 5 , n = 1 7257600 h ( 460745 h 2 g m , n + 7383268 h 2 g m + 1 , n + 5603124 h 2 g m + 2 , n + 9429308 h 2 g m + 3 , n + 5688310 g m + 4 , n h 2 + 4847580 g m + 5 , n h 2 934828 h 2 g m + 6 , n + 202756 h 2 g m + 7 , n 21063 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 6 , n = 1 7257600 h ( 454279 h 2 g m , n + 7446456 h 2 g m + 1 , n + 5318936 h 2 g m + 2 , n + 10225024 h 2 g m + 3 , n + 4022070 g m + 4 , n h 2 + 10111256 g m + 5 , n h 2 + 2341744 h 2 g m + 6 , n + 3504 h 2 g m + 7 , n 6469 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 7 , n = 1 7257600 h ( 468873 h 2 g m , n + 7308644 h 2 g m + 1 , n + 5907508 h 2 g m + 2 , n + 8714940 h 2 g m + 3 , n + 6656630 g m + 4 , n h 2 + 6606172 g m + 5 , n h 2 + 8831316 h 2 g m + 6 , n + 2754692 h 2 g m + 7 , n 74375 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) ,
u m + 8 , n = 1 7257600 h ( 400967 h 2 g m , n + 7934392 h 2 g m + 1 , n + 3325080 h 2 g m + 2 , n + 15007616 h 2 g m + 3 , n 3409610 g m + 4 , n h 2 + 17796888 g m + 5 , n h 2 377872 h 2 g m + 6 , n + 11688880 h 2 g m + 7 , n + 2065659 h 2 g m + 8 , n 7257600 u m , n + 7257600 u m + 1 , n ) .

Implementation Strategies of the NMCFs

Here, we combined the discrete Formulas ( 11 ) ( 17 ) and its derivatives ( 19 ) in matrix form below:
P U m = Q η 0 + R η 1 + h 2 [ S η 2 + T η 3 ] ,
where
P = 7257600 3628800 0 0 0 0 0 0 3628800 0 1209600 0 0 0 0 0 2419200 0 0 604800 0 0 0 0 1814400 0 0 0 362880 0 0 0 1451520 0 0 0 0 241920 0 0 1209600 0 0 0 0 0 172800 0 1036800 0 0 0 0 0 0 129600 7257600 0 0 0 0 0 0 0 , U m = u m + 1 , n u m + 2 , n u m + 3 , n u m + 4 , n u m + 5 , n u m + 6 n u m + 7 , n u m + 8 , n
Q = 0 0 0 0 0 0 0 3628800 0 0 0 0 0 0 0 2419200 0 0 0 0 0 0 0 1814400 0 0 0 0 0 0 0 1451520 0 0 0 0 0 0 0 1209600 0 0 0 0 0 0 0 1036800 0 0 0 0 0 0 0 907200 0 0 0 0 0 0 0 7257600 , η 0 = u m 1 , n u m 2 , n u m 3 , n u m 4 , n u m 5 , n u m 6 n u m 7 , n u m , n
R = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7257600 , η 1 = u m 1 , n u m 2 , n u m 3 , n u m 4 , n u m 5 , n u m 6 n u m 7 , n u m , n
S = 0 0 0 0 0 0 0 237671 0 0 0 0 0 0 0 155171 0 0 0 0 0 0 0 115821 0 0 0 0 0 0 0 92405 0 0 0 0 0 0 0 76859 0 0 0 0 0 0 0 65871 0 0 0 0 0 0 0 57281 0 0 0 0 0 0 0 1624505 , η 2 = g m 1 , n g m 2 , n g m 3 , n g m 4 , n g m 5 , n g m 6 n g m 7 , n g m , n
T = 3398072 653032 1426304 1376650 884504 368272 90032 9829 2374092 579388 1008404 855150 543604 225892 55212 6029 1803752 752008 1256064 563950 394504 165552 40592 4439 1452692 725156 1245452 132610 351020 134332 32564 3547 1215576 665656 1157504 73470 504376 87376 25584 2849 1044092 608428 1053924 177050 561604 84108 36572 2849 918152 542888 992384 177050 623144 149648 162512 5741 4124232 5225624 6488192 5888310 3698920 1522672 369744 40187
η 3 = g m + 1 , n g m + 2 , n g m + 3 , n g m + 4 , n g m + 5 , n g m + 6 n g m + 7 , n g m + 8 , n .
The NMCFs are obtained by multiplying the matrix Equation ( 28 ) by the inverse of P:
U m = Q ¯ η 0 + R ¯ η 1 + h 2 [ S ¯ η 2 + T ¯ η 3 ]
U m = 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 , Q ¯ = 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1
R ¯ = h 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 7 0 0 0 0 0 0 0 8 , S ¯ = 0 0 0 0 0 0 0 324901 1451520 0 0 0 0 0 0 0 58193 113400 0 0 0 0 0 0 0 71661 89600 0 0 0 0 0 0 0 15406 14175 0 0 0 0 0 0 0 56975 41472 0 0 0 0 0 0 0 93 56 0 0 0 0 0 0 0 2019731 1036800 0 0 0 0 0 0 0 31648 14175
T ¯ = 8183 14400 653203 907200 50689 56700 196277 241920 92473 181440 95167 453600 7703 151200 5741 1036800 29384 14175 81 50 30916 14175 22703 11340 5968 4725 14773 28350 1796 14175 521 37800 1467 400 4707 2800 225 64 28143 8960 11079 5600 9141 11200 2223 11200 387 17920 24832 4725 928 567 80128 14175 188 45 38144 14175 15776 14175 256 945 418 14175 248375 36288 19375 12096 143375 18144 641875 145152 225 64 12875 9072 3125 9072 3625 96768 1476 175 549 350 1776 175 639 140 36 7 81 50 72 175 9 200 216433 21600 98441 64800 1601467 129600 160867 34560 55223 8100 127253 129600 8183 14400 57281 1036800 23552 2025 7424 4725 41984 2835 14528 2835 41984 4725 7424 14175 23552 14175 0 .
The matrix equation in ( 29 ) could be written explicitly as follows:
u m + 1 , n = u m , n + u m , n h + 324901 1451520 h 2 g m , n + 8183 14400 h 2 g m + 1 , n 653203 907200 h 2 g m + 2 , n + 50689 56700 h 2 g m + 3 , n 196277 241920 h 2 g m + 4 , n + 92473 181440 h 2 g m + 5 , n 95167 453600 h 2 g m + 6 , n + 7703 151200 h 2 g m + 7 , n 5741 1036800 h 2 g m + 8 , n ,
u m + 2 , n = u m , n + 2 h u m , n + 58193 113400 h 2 g m , n + 29384 14175 h 2 g m + 1 , n 81 50 h 2 g m + 2 , n + 30916 14175 h 2 g m + 3 , n 22703 11340 h 2 g m + 4 , n + 5968 4725 h 2 g m + 5 , n 14773 28350 h 2 g m + 6 , n + 1796 14175 h 2 g m + 7 , n 521 37800 h 2 g m + 8 , n ,
u m + 3 , n = u m , n + 3 h u m , n + 71661 89600 h 2 g m , n + 1467 400 h 2 g m + 1 , n 4707 2800 h 2 g m + 2 , n + 225 64 h 2 g m + 3 , n 28143 8960 h 2 f m + 4 , n + 11079 5600 h 2 g m + 5 , n 9141 11200 h 2 g m + 6 , n + 2223 11200 h 2 g m + 7 , n 387 17920 h 2 g m + 8 , n ,
u m + 4 , n = u m , n + 4 h u m , n + 15406 14175 h 2 g m , n + 24832 4725 h 2 g m + 1 , n 928 567 h 2 g m + 2 , n + 80128 14175 h 2 g m + 3 , n 188 45 h 2 f m + 4 , n + 38144 14175 h 2 g m + 5 , n 15776 14175 h 2 g m + 6 , n + 256 945 h 2 g m + 7 , n 418 14175 h 2 g m + 8 , n ,
u m + 5 , n = u m , n + 5 h u m , n + 56975 41472 h 2 g m , n + 248375 36288 h 2 g m + 1 , n 19375 12096 h 2 g m + 2 , n + 143375 18144 h 2 g m + 3 , n 641875 145152 h 2 g m + 4 , n + 225 64 h 2 g m + 5 , n 12875 9072 h 2 g m + 6 , n + 3125 9072 h 2 f m + 7 , n 3625 96768 h 2 g m + 8 , n ,
u m + 6 , n = u m , n + 6 h u m , n + 93 56 h 2 g m , n + 1476 175 h 2 g m + 1 , n 549 350 h 2 g m + 2 , n + 1776 175 h 2 g m + 3 , n 639 140 h 2 g m + 4 , n + 36 7 h 2 g m + 5 , n 81 50 h 2 f m + 6 , n + 72 175 h 2 g m + 7 , n 9 200 h 2 g m + 8 , n ,
u m + 7 , n = u m , n + 7 h u m , n + 2019731 1036800 h 2 g m , n + 216433 21600 h 2 g m + 1 , n 98441 64800 h 2 g m + 2 , n + 1601467 129600 h 2 g m + 3 , n 160867 34560 h 2 g m + 4 , n + 55223 8100 h 2 g m + 5 , n 127253 129600 h 2 g m + 6 , n + 8183 14400 h 2 g m + 7 , n 57281 1036800 h 2 g m + 8 , n ,
u m + 8 , n = u m , n + 8 h u m , n + 31648 14175 h 2 g m , n + 23552 2025 h 2 g m + 1 , n 7424 4725 h 2 g m + 2 , n + 41984 2835 h 2 g m + 3 , n 14528 2835 h 2 g m + 4 , n + 41984 4725 h 2 g m + 5 , n 7424 14175 h 2 g m + 6 , n + 23552 14175 h 2 g m + 7 , n .
In order to obtain the first derivatives of the NMCFs, substitute ( 30 ) ( 37 ) into ( 20 ) ( 27 ) , which gives
u m + 1 , n = u m , n + 1070017 3628800 h g m , n + 2233547 1814400 h g m + 1 , n 2302297 1814400 h g m + 2 , n + 2797679 1814400 h g m + 3 , n 31457 22680 h g m + 4 , n + 1573169 1814400 h g m + 5 , n 645607 1814400 h g m + 6 , n + 156437 1814400 h g m + 7 , n 33953 3628800 h g m + 8 , n ,
u m + 2 , n = u m , n + 32377 113400 h g m , n + 22823 14175 h g m + 1 , n 21247 56700 h g m + 2 , n + 15011 14175 h g m + 3 , n 2903 2835 h g m + 4 , n + 9341 14175 h g m + 5 , n 15577 56700 h f m + 6 , n + 953 14175 h g m + 7 , n 119 16200 h g m + 8 , n ,
u m + 3 , n = u m , n + 12881 44800 h g m , n + 35451 22400 h g m + 1 , n + 1719 22400 h g m + 2 , n + 39967 22400 h g m + 3 , n 351 280 h g m + 4 , n + 17217 22400 h g m + 5 , n 7031 22400 h g m + 6 , n + 243 3200 h g m + 7 , n 369 44800 h g m + 8 , n ,
u m + 4 , n = u m , n + 4063 14175 h g m , n + 22576 14175 h g m + 1 , n + 244 14175 h g m + 2 , n + 32752 14175 h g m + 3 , n 1816 2835 h g m + 4 , n + 9232 14175 h g m + 5 , n 3956 14175 h g m + 6 , n + 976 14175 h g m + 7 , n 107 14175 h g m + 8 , n ,
u m + 5 , n = u m , n + 41705 145152 h g m , n + 115075 72576 h g m + 1 , n + 3775 72576 h g m + 2 , n + 159175 72576 h g m + 3 , n 125 4536 h g m + 4 , n + 85465 72576 h g m + 5 , n 24575 72576 h f m + 6 , n + 5725 72576 h g m + 7 , n 175 20736 h f m + 8 , n ,
u m + 6 , n = u m , n + 401 1400 h g m , n + 279 175 h g m + 1 , n + 9 700 h f m + 2 , n + 403 175 h g m + 3 , n 9 35 h g m + 4 , n + 333 175 h f m + 5 , n + 79 700 h g m + 6 , n + 9 175 h g m + 7 , n 9 1400 h g m + 8 , n ,
u m + 7 , n = u m , n + 149527 518400 h g m , n + 408317 259200 h g m + 1 , n + 24353 259200 h g m + 2 , n + 542969 259200 h g m + 3 , n + 343 3240 h g m + 4 , n + 368039 259200 h g m + 5 , n + 261023 259200 h f m + 6 , n + 111587 259200 h g m + 7 , n 8183 518400 h g m + 8 , n ,
u m + 8 , n = u m , n + 3956 14175 h g m , n + 23552 14175 h g m + 1 , n 3712 14175 h g m + 2 , n + 41984 14175 h g m + 3 , n 3632 2835 h g m + 4 , n + 41984 14175 h g m + 5 , n 3712 14175 h f m + 6 , n + 23552 14175 h g m + 7 , n + 3956 14175 h f m + 8 , n .
Remark 1.
The matrix D in Equation ( 9 ) must be a square matrix. Otherwise, the computation will be indeterminable.
Remark 2.
The matrix D in Equation ( 9 ) must be a nonsingular matrix. Otherwise, the determinant will be equal to zero, which implies that the solution will not exist.
Remark 3.
The Equations ( 30 ) ( 37 ) and ( 38 ) ( 44 ) formed the NMCFs required to simultaneously solve the second-order PDEs.
Remark 4.
The Equations ( 30 ) ( 37 ) and ( 38 ) ( 44 ) possessed a uniform order of accuracy as shown below.

3. Theoretical Analysis of the NMCFs

3.1. Introduction of the Analysis of the NMCFs

The order, error constants, zero-stability, and consistency of the NMCFs, will be examined in this section in accordance with (Familua et al. [35], Jain et al. [36]).
The Formula ( 18 ) together with their associated derivatives ( 29 ) can be expressed by the linear operator below:
L [ u ( x ) ; h ] = U m Q ¯ η 0 R ¯ η 1 h 2 [ S ¯ η 2 + T ¯ η 3 ] ,
where u m is continuously differentiable, and U m , Q ¯ η 0 , R ¯ η 1 , S ¯ η 2 , and T ¯ η 3 have their usual meaning as stated above. Expanding U m , η 2 , and η 3 in ( 46 ) , respectively, in Taylor series about x n , and collecting their like terms in powers of h and u yields
L [ u ( x ) ; h ] = C 0 u x + C 1 h u x + C 2 h 2 u x + + C q h q u ( q ) x ,
where C q , q = 1 , 2 ,
Definition 1
(Yakubu et al. [21]). The NMCFs ( 18 ) and their linear operators are said to have order p if C 0 = C 1 = = C p = 0 , C p + 2 0 .
Definition 2
(Sunday et al. [29]). The term C p + 2 in Definition 1 is referred to as the error constants, which indicate the local truncation error ( 18 ) given as
L T E = C p + 2 h p + 2 u ( p + 2 ) x n + O h ( p + 3 ) .
Definition 3
(Modebei et al. [32]). Any LMM class with an order greater than or equal to one is said to be consistent.
Definition 4
(Olaiya et al. [34]). If the roots of any class of LMM do not exceed the order of the differential equations considered, the class is said to be zero-stable.
Definition 5
(Lambert [37]). If an LMM class is zero-stable and consistent, it is said to be convergent.
Definition 6
(Lambert [37]). A matrix whose determinant is zero is called a singular matrix.
Definition 7
(Jain et al. [38]). If the Linear Multistep Method has order 2 k , where k is even, and order 2 k 1 , where k is odd, it is said to be of maximal order.
Definition 8
(Jain et al. [38]). A matrix whose determinant is not equal to zero is called a nonsingular matrix.
Definition 9
(Awoyemi and Idowu [12]). An LMM is said to be P-stable if its periodicity interval is ( 0 , ) .
Definition 10
(Jain et al. [38]). A matrix that has one in the leading diagonal and zero elsewhere.
Definition 11
(Henricin [39]). An LMM is said to be A-stable if its periodicity interval is ( , 0 ) .

3.1.1. Order and Error Constant of the NMCFs

The order and error constant of the NMCFs are analysed following the approach and procedure discussed in Definition 1. Each of ( 30 ) ( 37 ) , which make up the NMCFs, is analysed.
Hence, the NMCFs is of order p = [ 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 , 9 ] T , with error constants C p + 2 =
C 11 = 4.718068 × 10 3 , 1.179146 × 10 2 , 1.851157 × 10 2 , 2.531131 × 10 2 , T
3.211106 × 10 2 , 3.883117 × 10 2 , 4.590456 × 10 2 , 5.062263 × 10 2 T ,
together with its derivatives ( 38 ) ( 44 ) :
7.892554 × 10 3 , 6.428571 × 10 3 , 6.975446 × 10 3 , 6.631393 × 10 3 , T
6.975446 × 10 3 , 6.428571 × 10 3 , 7.892554 × 10 3 , 0 T .

3.1.2. Consistency of the NMCFs (Omole et al. [40,41])

Applying the definition 3 of consistency to the NMCFs ( 30 ) ( 37 ) , it is said to be consistent if it has an order of more than or equal to one. Therefore, the NMCFs are consistent according to Lambert [37].

3.1.3. Zero-Stability of the NMCFs (Fatunla [42])

Similarly, the zero-stability of the NMCFs can be obtained using the first characteristics polynomial of the NMCFs given by
Π ( z ) = det z U m Q ¯ = 0 .
Thus,
Π ( z ) = z 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 = 0
Π ( z ) = z 7 ( z 1 ) = 0 ,
Solving ( 50 ) for the values of z, z = 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 .
Hence, the NMCFs are zero-stable.

3.1.4. Convergence of the NMCFs (Henrici [39], Rufai et al. [43])

According to Definition 5, consistency and zero-stability are all that are required for the Linear Multistep Method to be convergent, so the NMCFs’ consistency and zero-stability imply that they are convergent at all points, concluding the prove.

3.1.5. Region of Absolute Stability of the NMCFs (Yakubu et al. [21], Lambert [37])

Finally, the stability of the NMCFs is examined and discussed using the procedure described in Yakubu et al. [21] and Lambert [37].
M ( z ) = V + z B ( I z A ) 1 U ,
In addition to the stability function,
p ( n , z ) = d e t ( n I M ( z ) ) ,
For the stability properties, the Formulas ( 30 ) ( 37 ) were formulated as
Y Y i + 1 = A U B V h 2 f ( u ) Y i 1 ,
where
Y i 1 = u m + 1 , n u m , n , Y i + 1 = u m + 1 , n u m + 8 , n , V = 0 1 0 1 , I = 1 0 0 1 ,
The stability polynomial ( 54 ) and its first derivative ( 55 ) are then obtained by changing the values of A, B, U, V, M, and I (see Appendix C) in Equations ( 51 ) and ( 52 ) . This is then coded in the MATLAB (R2012a) environment. Figure 1 in the text below depicts the NMCFs’ stability nature.
f z = η + 963008 42525 z 2 + 390304 14175 z 1 η 7
f z = 32 42525 60188 z + 36591 η 7
The NMCFs’ absolute stability region is P-stable because it consists of the complex plane outside the enclosed figure and its periodicity interval lies between ( 0.069 , 0 ) , which falls within the periodicity interval for P-stability ( , 0 ) .

4. Implementation Strategy

In this section, the implementation strategy is discussed in detail. It follows that variable y is discretized as follows:
h = b a M , x i = a + i h , i = 0 , 1 , , M . , h = d c N , t i = c + i h , i = 0 , 1 , , N .
For a fixed x in the interval [ a , b ] , i = 0 , , M , and for a fixed t in the interval [ c , d ] , i = 0 , , N .
N is the number of subintervals or iterations. The spatial derivative is approximated by the difference operator and replaced accordingly:
u t u ( x m + 1 , t ) u ( x m 1 , t ) 2 h ,
u y u ( x , y m + 1 ) u ( x , y m 1 ) 2 h ,
2 u y 2 u ( x , y m + 1 ) 2 u ( x , y m ) + u ( x , y m 1 , y ) ( h ) 2 ,
u ( x , y m + 1 ) is termed the numerical approximation to u ( x , y m + 1 ) and, consequently, ( 1 ) has been converted to a system of ordinary differential equations and ( 1 ) has the semi-discretized form given below:
d 2 u m , n d x 2 = 1 a m , n [ b m , n u ( x , y m + 1 ) 2 u ( x , y m ) + u ( x , y m 1 , y ) ( h ) 2 p m , n d u m , n d x q m , n u ( x m + 1 , t ) u ( x m 1 , t ) 2 h r m , n u ( x , y m + 1 ) u ( x , y m 1 ) 2 h k m , n u m , n + g m , n ] ,
where it follows that g are the non-homogeneous terms, h is the step size, and W is the tridiagonal matrix generated from ( 60 ) (shown in Appendix D).
The proposed formulas, namely the NMCFs, are then applied to solve the resulting equations of the ODEs with initial or boundary conditions ( 60 ) , with the aid of Mathematica 11.0, with the features Nsolve for linear and Findroot for nonlinear, with an HP Laptop G250, 8GB RAM, ITERABYTE.

5. Numerical Examples

The accuracy and convergence of the NMCFs are presented in this section. Five numerical examples including nonlinear partial differentials originating from fluid dynamics were resolved from the literature. A comparison of the numerical solution produced by the NMCFs, and the exact results were made and also compared with the errors produced by other existing methods in the literature. To highlight the accuracy of the NMCFs and their benefits over existing techniques, the results are presented in tabular form.
For instance, the absolute error of the approximate solutions is computed alongside the exact solution and compared with the results from the other existing methods, particularly those proposed by Yagider and Karabacak [44], Lima et al. [45], Biala and Jator [46], Xu and Wang [47], and Volkov et al. [48]. The results from the methods are also discussed here.
The absolute errors (AEs) are given by = Max | u ( x m , y m ) u m ( y n ) |. u ( x m , y m ) symbolises the exact solution, and u m ( y n ) denotes the approximate solution at the mesh point ( x m , y m ) .

5.1. Test Problem One

Consider the following Helmholtz equation:
2 u x 2 ( x , y ) + 2 u y 2 ( x , y ) + 8 u ( x , y ) = 0 , 0 < x , y < 1 ,
subjected to the following conditions:
u ( 0 , y ) = sin ( 2 y ) , for 0 x 1 , u x ( 0 , y ) = 0 , for 0 y 1 ,
the theoretical solution of ( 61 ) is given as
u ( x , y ) = cos ( 2 x ) sin ( 2 y ) .
Source: Yagider and Karabacak [44].
The prominent Helmholtz equation, which is an example of the elliptic PDEs considered in test problem one, was solved using the newly formulated formulas named NMCFs. The numerical solution is demonstrated in Table 1. The problem had been solved earlier by Yagider and Karabacak [44]. The numerical results were generated at x-values ranging between 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 , and 1.0 . The NMCFs results obtained show some level of improvement over the method constructed by Yagider and Karabacak [44], who proposed a multivariate Padé approximation as displayed in Table 2. Seemingly, the results generated by the NMCFs compared favorably with Yagider and Karabacak [44], despite using the same mesh size. A comparison of the errors is also displayed in Figure 2 below to show the behaviour of the performance of the NMCFs versus Yagider and Karabacak [44].

5.2. Test Problem Two

To demonstrate the accuracy of various numerical methods and gain a conceptual understanding of physical flows, the modified Burgers–Fisher equation, which combines the reaction, convection, and diffusion mechanism, has been taken into consideration by a number of authors. The Burgers–Fisher equation refers to this equation because it combines characteristics of the Fisher equation for diffusion, transport, and reactions with those of the Burgers–Fisher equation for the convective phenomenon. We are also interested in analyzing the numerical solution of the modified Burgers–Fisher equation because it exhibits relatively quick convergence and accuracy, and this will help to show the accuracy of the suggested method. If we take into account the modified Burgers–Fisher equation formed below:
u t + u u x = 2 u x 2 + u ( 1 u ) , ( x , t ) u [ 0 , T ] × ( 0 , T ] , T > 0 ,
subjected to the following initial conditions:
u ( x , 0 ) = 1 2 1 2 t a n h ( x 4 ) , x Ω
and the boundary conditions given by
u ( 0 , t ) = 1 2 + 1 2 t a n h ( 5 t 8 ) , t ( 0 , T ] , x δ Ω , u ( 1 , t ) = 1 2 + 1 2 t a n h ( 5 t 8 1 4 ) , t ( 0 , T ] , x δ Ω
the theoretical solution of ( 64 ) is given as
u ( x , y ) = 1 2 1 2 t a n h ( 5 t 8 x 4 ) .
Source: Lima et al. [45].
The second test problem considered belongs to a class of convection–diffusion reaction equations known as the modified Burgers–Fisher equation. The numerical solution is taken into consideration as a means of applying the proposed formulas (NMCFs). This problem was earlier solved by Lima et al. [45], who proposed a finite element method. In Table 3, the numerical results were presented, while, in Table 4, the comparison of the absolute errors in the NMCFs against that of Lima et al. [45] is illustrated within the interval of integration. The results revealed that the performance of the NMCFs is far better than the method of Lima et al. [45] in terms of accuracy and give better convergence. The comparison of the errors is displayed in Figure 3 for better interpretation and reader understanding.

5.3. Test Problem Three

The Burgers equation is a convection–diffusion equation that can explain the evolutionary process by which a convective phenomenon can maintain balance with a diffusive behavior in a variety of areas of applied mathematics to comprehend the Navier–Stokes equations’ fundamental characteristics. The pressure can be ignored in this simple equation, but the nonlinear and viscous terms still have an impact. The Reynolds number in the Navier–Stokes equations is the proportion of a flow’s advective to viscous contributions. It is crucial to take into account an accurate and reliable numerical method for the simulation in order to obtain the evolution of this flow. The behavior of the modified Burgers equation in the following form was examined by the authors in Lima et al. [45].
u t + u u x = 1 4 2 u x 2 , ( x , t ) u [ 0 , T ] × ( 0 , T ] , T > 0 ,
subjected to the following initial conditions:
u ( x , 0 ) = 1 ( 1 + e 2 x ) x Ω ,
and the boundary conditions given by
u ( 0 , t ) = 1 ( 1 + e t ) , t ( 0 , T ] , x δ Ω , u ( 1 , t ) = 1 ( 1 + e 2 t ) , t ( 0 , T ] , x δ Ω
the theoretical solution of (68) is given as
u ( x , y ) = 1 ( 1 + e ( 2 x t ) 4 μ 0 ) .
Source: Lima et al. [45].
Next, in test problem three, the second-order nonlinear Navier–Stokes equations, which are a classical example of convection–diffusion equations, were put into consideration and solved using the NMCFs. The numerical solution, which comprises the exact results, NMCFs results, and the absolute errors, are shown in Table 5. The problem had been solved earlier by Lima et al. [45]. The numerical solution was executed at x-values ranging between 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 , and 1.0 . The NMCFs results generated show some level of outstanding performance over the method proposed by Lima et al. (2021), who constructed a finite element method with error correlation. Apparently, the results generated by the NMCFs give a minimal error against the method proposed by Lima et al. [45], as seen in Table 6. From the results, one can see that the proposed method has demonstrated excellently, with good performance with varying mesh points. The behavior of the performance of test problem three with the NMCFs and the method proposed by the authors in Lima et al. [45] is presented in plot form for better clarification in Figure 4.

5.4. Test Problem Four

Consider the Laplace equation given below:
2 u x 2 + 2 u y 2 = 0 , x , y [ 0 , 1 ] ,
subjected to the conditions on the boundary domain given below:
u ( x , 0 ) = e x , u y ( 0 , x ) = cos ( y ) ,
the exact solution of ( 72 ) is given as
u ( x , y ) = e x cos ( y ) .
Source: Xu and Wang [47].
Similarly, in test problem four, the Laplace equation was solved using the NMCFs. The numerical solution, which comprises the exact results, NMCFs results. and the absolute errors are shown in Table 7 and Table 8. Before now, Biala et al. [46] and Xu and Wang [47] had solved the problem with their own approach. For instance, Biala et al. [46] proposed a block unification algorithm of order two via a shifted Chebychev’s polynomial. Meanwhile, Xu and Wang [47] constructed a method named a parallel iterative algorithm. The numerical solution was implemented at N-values ranging between 32 , 40 , and 48. The NMCFs obtained show better performances in terms of accuracy and convergence over the method proposed by the existing methods. The behavior of the performance of test problem four with the NMCFs and the method proposed by the authors in Biala et al. [46] and Xu and Wang [47] are presented in a curve for better interpretation in Figure 5.

5.5. Test Problem Five

Finally, the PDE with the Dirichlet boundary condition as follows is considered:
2 u x 2 + 2 u y 2 = 2 π 2 π y 2 2 π y e π ( 1 y ) sin ( π x ) , x , y [ 0 , 1 ] ,
subjected to the conditions on the boundary domain given below:
u ( x , 0 ) = u ( x , 1 ) ,
the exact solution of ( 74 ) is given as
u ( x , y ) = e π x sin ( π x ) + e π ( 1 y ) sin ( π x ) .
Source: Volkov et al. [48].
Lastly, test problem five is solved using the NMCFs. The numerical results containing the exact results, NMCFs results, and the absolute errors are shown in Table 9 and Table 10. Biala and Jator [46] and Volkov et al. [48] had solved the problem with their own techniques before now. In the work of Biala and Jator [46], a block unification algorithm of order two via a shifted Chebychev’s polynomial was constructed, analyzed, and implemented. Likewise, Volkov et al. [48] presented a parallel iterative algorithm. The numerical solution was implemented at N-values ranging between 16 , 32 , 64 , and 128. The NMCFs obtained show superiority in performance in terms of accuracy over the method proposed by the existing methods. The behavior of the performance of test problem five with the NMCFs and the method proposed by the authors in Biala and Jator [46] and Volkov et al. [48] are presented in a figure for better usage in Figure 6.

6. Conclusions

This article presents a ninth-order multistep collocation formula that was designed using orthogonal polynomial collocation. The formulas derived from the continuous function were combined in a step-by-step block approach algorithm. The proposed formulas (NMCFs) were used to solve a class of partial differential equations ranging from the Helmholtz equation to the convection diffusion reaction equations resulting from the semi-discretization of the problems studied. Because the NMCFs were implemented block by block, they do not require the starting values and predictors that are associated with the predictor-corrector method. The numerical results demonstrate a significant improvement over other methods in the literature.
When compared to existing methods in the literature, the numerical experiments presented in this paper clearly show that the NMCFs have a reasonably wide stability region displayed in Figure 1 and enjoy good accuracy and fast convergence advantages. Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 show the comparison of errors for test problems 1–5. Because they are more accurate, provide better convergence, and are computationally faster, the NMCFs can be extended to solve problems that model real-life phenomena with numerous applications of higher-order PDEs in the future.

Author Contributions

The final manuscript was written by all the authors, who also read and approved it. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No.: 12171435).

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their gratitude to the anonymous referees for their thorough review of the manuscript and for their insightful comments and suggestions, which helped to improve the manuscript’s presentation.

Conflicts of Interest

The authors state that there are no conflicting interests.

Appendix A

Ψ 0 ( ϕ ) = ( ϕ + 6 ) Ψ 1 ( ϕ ) = ( ϕ + 7 ) Δ 0 ( ϕ ) = h 2 ϑ 7257600 2 ϕ 8 + 24 ϕ 7 + 99 ϕ 6 + 105 ϕ 5 189 ϕ 4 693 ϕ 3 + 1287 ϕ 2 9225 ϕ + 65871 Δ 1 ( ϕ ) = h 2 ϑ 1814400 4 ϕ 8 + 53 ϕ 7 + 223 ϕ 6 + 275 ϕ 5 593 ϕ 4 439 ϕ 3 5327 ϕ 2 + 37289 ϕ 261023 Δ 2 ( ϕ ) = h 2 ϑ 1814400 14 ϕ 8 + 203 ϕ 7 + 913 ϕ 6 + 1145 ϕ 5 2033 ϕ 4 4777 ϕ 3 + 67 ϕ 2 11917 ϕ + 152107 Δ 3 ( ϕ ) = h 2 ϑ 1814400 28 ϕ 8 + 441 ϕ 7 + 2181 ϕ 6 + 2925 ϕ 5 4551 ϕ 4 12405 ϕ 3 2913 ϕ 2 + 29679 ϕ 263481 Δ 4 ( ϕ ) = h 2 ϑ 725760 14 ϕ 8 + 238 ϕ 7 + 1313 ϕ 6 + 2095 ϕ 5 2875 ϕ 4 7523 ϕ 3 3631 ϕ 2 + 10369 ϕ + 17705 Δ 5 ( ϕ ) = h 2 ϑ 1814400 28 ϕ 8 + 511 ϕ 7 + 3161 ϕ 6 + 6445 ϕ 5 5791 ϕ 4 25937 ϕ 3 + 2063 ϕ 2 + 4135 ϕ 140401 Δ 6 ( ϕ ) = h 2 ϑ 1814400 14 ϕ 8 + 273 ϕ 7 + 1893 ϕ 6 + 5025 ϕ 5 393 ϕ 4 21729 ϕ 3 14757 ϕ 2 + 46059 ϕ + 21027 Δ 7 ( ϕ ) = h 2 ϑ 1814400 4 ϕ 8 + 83 ϕ 7 + 643 ϕ 6 + 2195 ϕ 5 + 2167 ϕ 4 5827 ϕ 3 17183 ϕ 2 13567 ϕ 9143 Δ 8 ( ϕ ) = h 2 ϑ 7257600 2 ϕ 8 + 44 ϕ 7 + 379 ϕ 6 + 1625 ϕ 5 + 3571 ϕ 4 + 3515 ϕ 3 + 343 ϕ 2 889 ϕ 2849 ϑ = ( ϕ + 6 ) ( ϕ + 7 )

Appendix B

Ψ 0 ( ϕ ) = 1 , Ψ 1 ( ϕ ) = 1 Δ 0 ( ϕ ) = h 2 7257600 20 ϕ 9 + 450 ϕ 8 + 3960 ϕ 7 + 16800 ϕ 6 + 32004 ϕ 5 + 6300 ϕ 4 62640 ϕ 3 64800 ϕ 2 + 468873 Δ 1 ( ϕ ) = h 2 1814400 40 ϕ 9 + 945 ϕ 8 + 8640 ϕ 7 + 37800 ϕ 6 + 74088 ϕ 5 + 17010 ϕ 4 143760 ϕ 3 151200 ϕ 2 1827161 Δ 2 ( ϕ ) = h 2 1814400 7140 ϕ 9 + 3465 ϕ 8 + 33120 ϕ 7 + 150780 ϕ 6 + 307188 ϕ 5 + 84420 ϕ 4 589680 ϕ 3 635040 ϕ 2 + 1476877 Δ 3 ( ϕ ) = h 2 1814400 280 ϕ 9 + 7245 ϕ 8 + 72720 ϕ 7 + 348600 ϕ 6 + 750456 ϕ 5 + 256410 ϕ 4 1421280 ϕ 3 1587600 ϕ 2 2178735 Δ 4 ( ϕ ) = h 2 725760 140 ϕ 9 + 3780 ϕ 8 + 39960 ϕ 7 + 204120 ϕ 6 + 477036 ϕ 5 + 215460 ϕ 4 888720 ϕ 3 1058400 ϕ 2 + 665663 Δ 5 ( ϕ ) = h 2 1814400 280 ϕ 9 + 7875 ϕ 8 + 87840 ϕ 7 + 483000 ϕ 6 + 1264536 ϕ 5 + 847350 ϕ 4 2313360 ϕ 3 3175200 ϕ 2 1651543 Δ 6 ( ϕ ) = h 2 1814400 140 ϕ 9 + 4095 ϕ 8 + 48240 ϕ 7 + 287700 ϕ 6 + 866628 ϕ 5 + 921060 ϕ 4 1254960 ϕ 3 3175200 ϕ 2 + 2207829 Δ 7 ( ϕ ) = h 2 1814400 40 ϕ 9 + 1215 ϕ 8 + 15120 ϕ 7 + 98280 ϕ 6 + 346248 ϕ 5 + 572670 ϕ 4 7680 ϕ 3 1445040 ϕ 2 1814400 ϕ 688673 Δ 8 ( ϕ ) = h 2 7257600 20 ϕ 9 + 630 ϕ 8 + 8280 ϕ 7 + 58800 ϕ 6 + 243684 ϕ 5 + 590940 ϕ 4 + 784080 ϕ 3 + 453600 ϕ 2 74375

Appendix C

A = 0 0 0 0 0 0 0 0 0 324901 1451520 1169 2 653203 907200 50689 56700 196277 241920 92473 181440 95167 453600 7703 151200 5741 1036800 58193 113400 14692 7 81 50 30916 14175 22703 11340 5968 4725 14773 28350 1796 14175 521 37800 71661 89600 1467 400 4707 2800 225 64 28143 8960 11079 5600 9141 11200 2223 11200 387 17920 15406 14175 24832 4725 928 567 80128 14175 188 45 38144 14175 15776 14175 256 945 418 14175 56975 41472 248375 36288 19375 12096 143375 18144 641875 145152 225 64 12875 9072 3125 9072 3625 96768 93 56 1476 175 549 350 1776 175 639 140 36 7 81 50 72 175 9 200 2019731 1036800 216433 21600 98441 64800 1601467 129600 160867 34560 55223 8100 127253 129600 8183 14400 57281 1036800 31648 14175 23552 2025 7424 4725 41984 2835 14528 2835 41984 4725 7424 14175 23552 14175 0 , U = 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 ,
B = 324901 1451520 1169 2 653203 907200 50689 56700 196277 241920 92473 181440 95167 453600 7703 151200 5741 1036800 58193 113400 14692 7 81 50 30916 14175 22703 11340 5968 4725 14773 28350 1796 14175 521 37800 ,

Appendix D

W = 2 ( h ) 2 1 ( h ) 2 0 0 0 0 1 ( h ) 2 2 ( h ) 2 1 ( h ) 2 0 0 0 0 1 ( h ) 2 2 ( h ) 2 0 0 0 0 0 0 2 ( h ) 2 1 ( h ) 2 0 0 0 0 1 ( h ) 2 2 ( h ) 2 1 ( h ) 2 0 0 0 0 1 ( h ) 2 2 ( h ) 2

References

  1. Liu, X.; Ahsan, M.; Ahmad, M.; Nisar, M.; Liu, X.; Ahmad, I.; Ahmad, H. Applications of Haar Wavelet-Finite Difference Hybrid Method and Its Convergence for Hyperbolic Nonlinear Schrodinger Equation with Energy and Mass Conversion. Energies 2021, 14, 7831. [Google Scholar] [CrossRef]
  2. Kamran, M.; Abbas, M.; Majeed, A.; Emadifar, H.; Nazir, T. Numerical Sim-ulation of Time Fractional BBM-Burger Equation Using Cubic B-Spline Functions. Hindawi J. Funct. Spaces 2022, 2022, 2119416. [Google Scholar]
  3. Kilic, S.S.S.; Celik, E. Complex Solutions to the Higher-Order Nonlinear Boussinesq Type Wave Equation Transform. Ric. Mat. 2022, 1–8. [Google Scholar]
  4. Akbarov, S.D.; Mehdiyev, M.A.; Zeynalov, A.M. Dynamic of the moving ring-load acting in the interior of the bi-layered hollow cylinder with imperfect contact between the layers. TWMS J. Pure Appl. Math. 2021, 12, 223–242. [Google Scholar]
  5. Chen, W.; Wu, W.X.; Teng, Z.D. Complete dynamics in a nonlocal dispersal two-strain SIV epidemic model with vaccination and latent delays. Appl. Comput. Math. 2020, 19, 360–391. [Google Scholar]
  6. Singh, J.; Kumar, D.; Kumar, S. An efficient computational method for local fractional transport equation occurring in fractal porous media. Comput. Appl. Math. 2020, 39, 137. [Google Scholar] [CrossRef]
  7. Mirzaee, F.; Rezaei, S.; Samadyar, N. Numerical solution of two-dimensional stochastic time-fractional sine-Gordon equation on non-rectangular domains using finite difference and meshfree methods. Eng. Anal. Bound. Elem. 2021, 127, 53–63. [Google Scholar] [CrossRef]
  8. Khan, A.; Khan, H.; Gómez-Aguilar, J.F.; Abdeljawad, T. Existence and Hyers-Ulam stability for a nonlinear singular fractional differential equations with Mittag-Leffler kernel. Chaos Solitons Fractals 2019, 127, 422–427. [Google Scholar] [CrossRef]
  9. Aliev, F.A.; Aliyev, N.A.; Hajiyeva, N.S.; Mahmudov, N.I. Some mathematical problems and their solutions for the oscillating systems with liquid dampers: A review. Appl. Comput. Math. 2021, 20, 339–365. [Google Scholar]
  10. Adeyefa, E.O.; Omole, E.O.; Shokri, A.; Yao, S.-W. Hermite Fitted Block Integrator for Solving Second-Order Anisotropic Elliptic Type PDEs. Fractal Fract. 2022, 6, 497. [Google Scholar] [CrossRef]
  11. Su, X.; Zhou, Y. A Fast High-Order Predictor–Corrector Method on Graded Meshes for Solving Fractional Differential Equations. Fractal Fract. 2022, 6, 516. [Google Scholar] [CrossRef]
  12. Awoyemi, D.O.; Idowu, M.O. A class of hybrid collocation methods for third order ordinary differential equations. Int. J. Comput. Math. 2005, 82, 1–7. [Google Scholar] [CrossRef]
  13. Iskandarov, S.; Komartsova, E. On the influence of integral perturbations on the boundedness of solutions of a fourth-order linear differential equation. TWMS J. Pure Appl. Math. 2022, 13, 3–9. [Google Scholar]
  14. Ashry, H.; Abd-Elhameed, W.M.; Moatimid, G.M.; Youssri, Y.H. Spectral Treatment of One and Two Dimensional Second-Order BVPs via Certain Modified Shifted Chebyshev Polynomials. Int. J. Appl. Comput. Math. 2021, 7, 248. [Google Scholar] [CrossRef]
  15. Asif, M. Haar Wavelet Approach for Numerical Solution of Three-dimensional Partial Differential Equations. Ph.D. Dissertation, University of Peshawar, Peshawar, Pakistan, 2017. [Google Scholar]
  16. Guo, L.; Xu, Y. Energy conserving local discontinuous Galerkin methods for the nonlinear Schrödinger equation with wave operator. J. Sci. Comput. 2015, 65, 622–647. [Google Scholar] [CrossRef]
  17. Aziz, I.; Khan, F. A comparative study of numerical integra-tion based on haar wavelets and hybrid functions. Comput. Math. Appl. 2010, 59, 2026–2036. [Google Scholar]
  18. Shiralashetti, S.C.; Kantli, M.H.; Deshi, A.B. Biorthogonal wavelet-based full-approximation schemes for the numerical solution of elasto-hydrodynamic lubrication problems. J. Math. Model. 2018, 6, 105–122. [Google Scholar]
  19. Saparova, B.; Mamytova, R.; Kurbanbaeva, N.; Ahmedov, A.A. A Haar Wavelet Series Solution of Heat Equation with Involution. J. Adv. Res. Fluid -Chanics Therm. Sci. 2021, 86, 50–55. [Google Scholar] [CrossRef]
  20. Takei, Y.; Iwata, Y. Numerical Scheme Based on the Implicit Runge-Kutta Method and Spectral Method for Calculating Nonlinear Hyperbolic Evolution Equations. Axioms 2022, 11, 28. [Google Scholar] [CrossRef]
  21. Yakubu, D.G.; Aminu, M.; Tumba, P.; Abdulhameed, M. An efficient family of second derivative Runge-Kutta collocation methods for oscillatory systems. J. Niger. Math. Soc. 2018, 37, 111–138. [Google Scholar]
  22. Zhao, E.; Huang, J. Boundary treatment of implicit-explicit Runge-Kutta method for hyperbolic systems with source terms. J. Comput. Phys. 2020, 423, 109828. [Google Scholar] [CrossRef]
  23. Mall, S.; Chakraverty, S. Single layer Chebyshev neural network model for solving elliptic partial differential equations. Neural Process. Lett. 2017, 45, 825–840. [Google Scholar] [CrossRef]
  24. Ghaffar, F.; Badshah, N.; Khan, M.A.; Islam, S. Multigrid method for 2D Helmholtz equation using higher order finite difference scheme accelerated by Krylov subspace. J. Appl. Environ. Biol. Sci. 2014, 4, 169–179. [Google Scholar]
  25. Ge, Y. Multigrid method and fourth-order compact difference discretization scheme with unequal meshsizes for 3D Poisson equation. J. Comput. Phys. 2010, 229, 6381–6391. [Google Scholar] [CrossRef]
  26. Gupta, M.M.; Kouatchou, J.; Zhang, J. Comparison of second and fourth order discretizations for multigrid Poisson solver. J. Comput. Phys. 1997, 132, 226–232. [Google Scholar] [CrossRef]
  27. Mulla, M.A.M.; Gaweash, A.M.A.; Bakur, H.Y.I. Numerical Solution of Parabolic in Partial Differential Equations (PDEs) in One and Two Space Variable. J. Appl. Math. Phys. 2022, 10, 311–321. [Google Scholar] [CrossRef]
  28. Harari, I.; Hughes, T.J. Finite element methods for the Helmholtz equation in an exterior domain model problem. Comput. Methods Appl. Mech. Eng. 1991, 87, 59–96. [Google Scholar] [CrossRef]
  29. Sunday, S.J.; Shokri, A.; Marian, D. Variable Step Hybrid Block Method for the Approximation of Kepler Problem. Fractal Fract. 2022, 6, 343. [Google Scholar] [CrossRef]
  30. Ramos, H.; Vigo-Aguiar, J. An almost L-stable BDF-type method for the numerical solution of stiff ODEs arising from the method of lines. Inter. J. Numer. Methods Partial. Differ. Equations 2006, 23, 1110–1121. [Google Scholar] [CrossRef]
  31. Ngwane, F.F.; Jator, S.N. A Trigonometrically Fitted Block Method for Solving Oscillatory Second-Order Initial Value Problems and Hamiltonian Systems. Inter. J. Diff. Eqn. 2017, 2017, 9293530. [Google Scholar] [CrossRef]
  32. Modebei, M.I.; Adeniyi, R.B.; Jator, S.N. Numerical approximation of fourth-order PDEs using block unification method. J. Niger. Math. Soc. 2020, 39, 47–68. [Google Scholar]
  33. Jator, S.N. Block Unification Scheme for Elliptic, Telegraph, and Sine-Gordon Partial Differential Equations. Am. J. Comput. Math. 2015, 5, 175–185. [Google Scholar] [CrossRef]
  34. Olaiya, O.O.; Oduwole, H.K.; Odeyemi, J.K. Numerical solution of Black-Scholes Partial Differential Equation Using Direct Solution of Second-Order Ordinary Differential Equation With Two-Step Hybrid Block Method Of Order Seven. Sci. World J. 2019, 14, 23–29. [Google Scholar]
  35. Familua, A.B.; Omole, E.O.; Ukpebor, L.A. A Higher-order Block Method for Numerical Approximation of Third-order Boundary Value Problems in ODEs. J. Niger. Soc. Phys. Sci. 2022, 4, 706. [Google Scholar] [CrossRef]
  36. Milne, W.E. Numerical Solution of Differential Equations; John Wiley and Sons: Hoboken, NJ, USA, 1953. [Google Scholar]
  37. Lambert, J.D. Computational Methods in Ordinary Differential Equations; John Wiley & Sons Inc.: Hoboken, NJ, USA, 1973. [Google Scholar]
  38. Jain, M.K.; Iyengar, S.K.; Jain, R.K. Numerical Methods for Scientific and Engineering Computation, 5th ed.; New Age International: Dhaka, Bangladesh, 2007; pp. 282–283. Available online: https://books.google.co.jp/books?hl=en&lr=&id=5XappvcENCMC&oi=fnd&pg=PA1&dq=Jain,+M.K.%3B+Iyengar,+S.K.%3B+Jain,+R.K.+Numerical+Methods+for+Scientific+and+Engineering+Computation,+5th+ed.,+2007%3B+pp.+282%E2%80%93283.&ots=RWHgTT0nnz&sig=BY3hSRTPK_x_M-a6G7H73_LbKUk#v=onepage&q&f=false (accessed on 19 July 2023).
  39. Henrici, P. Discrete Variable Method in Ordinary Differential Equations; John Wiley & Sons: New York, NY, USA, 1972. [Google Scholar]
  40. Omole, E.O.; Jeremiah, O.A.; Adoghe, L.O. A Class of Continuous Im-plicit Seventh-eight method for solving y = f(x,y) using power series as basic function. Int. J. Chem. Math. Phys. (IJCMP) 2020, 4, 39–50. [Google Scholar] [CrossRef]
  41. Omole, E.O.; Obarhua, F.O.; Familua, A.B.; Shokri, A. Algorithms of algebraic order nine for numerically solving second-order boundary and initial value problems in ordinary differential equations. Int. J. Math. Oper. Res. 2022, 25, 343–368. [Google Scholar] [CrossRef]
  42. Fatunla, S.O. Numerical Methods for Initial Value Problems in Ordinary Differential Equations; Academic Press Inc.: Cambridge, MA, USA; Harcourt Brace Jovanovich Publishers: New York, NY, USA, 1988. [Google Scholar]
  43. Rufai, M.A.; Shokri, A.; Omole, E.O. A One-Point Third-Derivative Hy-brid Multistep Technique for Solving Second-Order Oscillatory and Periodic Problems. Hindawi J. Math. 2023, 2023, 2343215. [Google Scholar] [CrossRef]
  44. Yigider, M.; Karabacak, M. The numerical solution of helmholtz equation via multivariate pade approximation. Int. J. Res. Eng. Technol. 2014, 4, 670–678. [Google Scholar]
  45. Lima, S.A.; Kamrujjaman, M.; Islam1, M.S. Numerical solution of con-vection diffusion reaction equations by a finite element method with error correlation. AIP Adv. 2021, 11, 085225. [Google Scholar] [CrossRef]
  46. Biala, T.A.; Jator, S.N. Block unification algorithm for 2d and 3d elliptic PDEs. J. Niger. Math. Soc. 2017, 36, 319–335. [Google Scholar]
  47. Xu, Q.; Wang, W. A New Parallel Iterative Algorithm for Solving 2D Poisson Equation. Numer. Methods Partial. DiffErential Equations 2011, 27, 749–1001. [Google Scholar] [CrossRef]
  48. Volkov, E.A.; Dosiyev, A.A.; Buranay, S.C. On the solution of a nonlocal problem. Comput. Math. Appl. 2013, 66, 330–338. [Google Scholar] [CrossRef]
Figure 1. Region of absolute stability of the NMCFs.
Figure 1. Region of absolute stability of the NMCFs.
Axioms 12 00891 g001
Figure 2. Comparison of errors in NMCFs versus Yagider and Karabacak [44] for test problem 5.1.
Figure 2. Comparison of errors in NMCFs versus Yagider and Karabacak [44] for test problem 5.1.
Axioms 12 00891 g002
Figure 3. Comparison of errors in NMCFs versus Lima et al. [45] for test problem 5.2.
Figure 3. Comparison of errors in NMCFs versus Lima et al. [45] for test problem 5.2.
Axioms 12 00891 g003
Figure 4. Comparison of errors in NMCFs versus Lima et al. [45] for test problem 5.3.
Figure 4. Comparison of errors in NMCFs versus Lima et al. [45] for test problem 5.3.
Axioms 12 00891 g004
Figure 5. Comparison of errors in NMCFs versus Biala et al. [46] and Xu and Wang [47] for test problem 5.4.
Figure 5. Comparison of errors in NMCFs versus Biala et al. [46] and Xu and Wang [47] for test problem 5.4.
Axioms 12 00891 g005
Figure 6. Comparison of errors in NMCFs versus Biala and Jator [46] and Volkov et al. [48] for test problem 5.5.
Figure 6. Comparison of errors in NMCFs versus Biala and Jator [46] and Volkov et al. [48] for test problem 5.5.
Axioms 12 00891 g006
Table 1. Showing the numerical results for test problem 5.1.
Table 1. Showing the numerical results for test problem 5.1.
xExact ResultsNMCFs ResultsAEs in NMCFs
0.1 0.1947083556680159200 0.1947091711543252300 8.1548630 × 10 07
0.2 0.1829813709520361500 0.1829865712999870800 5.2003480 × 10 06
0.3 0.1639638132942570800 0.1639688742954361300 5.0610010 × 10 06
0.4 0.1384113212579074500 0.1384142557064305700 2.9344490 × 10 06
0.5 0.1073395314304028500 0.1073414975338517800 1.9661030 × 10 06
0.6 0.0719877982915558400 0.0719893725902818500 1.5742990 × 10 06
0.7 0.0337668596353071160 0.0337672585371394250 3.9890180 × 10 07
0.8 0.005800535872526165 0.005801049555132515 5.1368260 × 10 07
0.9 0.045138469521280010 0.045138088107911740 3.8141340 × 10 07
1.0 0.082675613529302500 0.082675613529302500 0.0000000
Table 2. Comparison of errors in NMCFs and Yagider and Karabacak [44] for test problem 5.1.
Table 2. Comparison of errors in NMCFs and Yagider and Karabacak [44] for test problem 5.1.
xAEs in NMCFsAEs in Yagider and Karabacak [44]
0.1 8.1548630 × 10 07 2.7000000 × 10 09
0.2 5.2003480 × 10 06 1.6800000 × 10 07
0.3 5.0610010 × 10 06 1.8551000 × 10 06
0.4 2.9344490 × 10 06 1.0017300 × 10 05
0.5 1.9661030 × 10 06 3.6361800 × 10 05
0.6 1.5742990 × 10 06 1.0226850 × 10 04
0.7 3.9890180 × 10 07 2.4043250 × 10 04
0.8 5.1368260 × 10 07 4.9437380 × 10 04
0.9 3.8141340 × 10 07 9.1526730 × 10 04
1.0 0.0000000 × 10 00 1.5558685 × 10 03
Table 3. Showing the computational results for test problem 5.2.
Table 3. Showing the computational results for test problem 5.2.
xExact ResultsNMCFs ResultsAEs in NMCFs
0.1 0.49106345176865770 0.491062975092073030 4.8 × 10 07
0.2 0.47732805754796130 0.477326955296469500 1.1 × 10 06
0.3 0.46225942047677904 0.462257661946340660 1.8 × 10 06
0.4 0.45270432190392235 0.452702242491803130 2.1 × 10 06
0.5 0.44047118471204494 0.440468940064822230 2.2 × 10 06
0.6 0.42696371653650966 0.426961857234471400 1.9 × 10 06
0.7 0.42696185723447140 0.416233741164731440 8.9 × 10 07
0.8 0.40425905953324964 0.404260159925367000 1.1 × 10 06
0.9 0.39108589473886196 0.391090346399732170 4.5 × 10 06
1.0 0.38195699940974250 0.381964156062750360 7.2 × 10 06
Table 4. Comparison of AEs in NMCFs against Lima et al. [45] for test problem 5.2.
Table 4. Comparison of AEs in NMCFs against Lima et al. [45] for test problem 5.2.
xAEs in NMCFsAEs in Lima et al. [45]
0.1 4.8 × 10 07 1.0 × 10 04
0.2 1.1 × 10 06 8.0 × 10 04
0.3 1.8 × 10 06 6.0 × 10 04
0.4 2.1 × 10 06 7.0 × 10 04
0.5 2.2 × 10 06 8.0 × 10 04
0.6 1.9 × 10 06 8.0 × 10 04
0.7 8.9 × 10 07 9.0 × 10 04
0.8 1.1 × 10 06 8.0 × 10 04
0.9 4.5 × 10 06 1.4 × 10 03
1.0 7.2 × 10 06 7.0 × 10 04
Table 5. Showing the results for test problem 5.3.
Table 5. Showing the results for test problem 5.3.
xExact ResultsNMCFs ResultsAEs in NMCFs
0.1 0.45338575868105696 0.4533846995902536 1.1 × 10 06
0.2 0.40492152613386595 0.40491952790249813 2.0 × 10 06
0.3 0.3582425862347542 0.35823985347997894 2.7 × 10 06
0.4 0.31410440319529837 0.3141011779927165 3.2 × 10 06
0.5 0.273090250573869 0.27308678350641447 3.5 × 10 06
0.6 0.2316531587086934 0.23164969524452944 3.4 × 10 06
0.7 0.1982925977359844 0.19828934296755002 3.3 × 10 06
0.8 0.1686815943978195 0.16867869569239505 2.9 × 10 06
0.9 0.1427053084191849 0.14270286639846103 2.4 × 10 06
1.0 0.1201511074867948 0.12014918056539994 1.9 × 10 06
Table 6. Comparison of errors in NMCFs and Lima et al. [45] for test problem 5.3.
Table 6. Comparison of errors in NMCFs and Lima et al. [45] for test problem 5.3.
xAEs in NMCFsAEs in Lima et al. [45]
0.1 1.1 × 10 06 1.9 × 10 03
0.2 2.0 × 10 06 2.4 × 10 03
0.3 2.7 × 10 06 2.1 × 10 03
0.4 3.2 × 10 06 2.0 × 10 03
0.5 3.5 × 10 06 1.9 × 10 03
0.6 3.4 × 10 06 1.7 × 10 03
0.7 3.3 × 10 06 1.5 × 10 03
0.8 2.9 × 10 06 1.3 × 10 03
0.9 2.4 × 10 06 1.3 × 10 03
1.0 1.9 × 10 06 3.0 × 10 04
Table 7. Showing the results for test problem 5.4.
Table 7. Showing the results for test problem 5.4.
NExact ResultsNMCFs ResultsAEs in NMCFs
16 1.0643645158455162 1.0643645183902837 2.5448 × 10 09
24 1.0419993959674516 1.0419993959681841 7.3253 × 10 09
32 1.0316174618605911 1.0316174648212710 2.9607 × 10 09
40 1.0253134960263788 1.0253134960268170 4.3832 × 10 13
48 1.0216581918640595 1.0216581918644914 4.3188 × 10 13
Table 8. Comparison of errors in NMCFs with Xu and Wang [47] and Biala and Jator [46] for test problem 5.4.
Table 8. Comparison of errors in NMCFs with Xu and Wang [47] and Biala and Jator [46] for test problem 5.4.
NAEs in NMCFsAEs in Biala and Jator [46]AEs in Xu and Wang [47]
16 2.5448 × 10 09 1.2800 × 10 06 3.9000 × 10 05
24 7.3253 × 10 09 2.6800 × 10 07 1.7400 × 10 05
32 2.9607 × 10 09 8.7600 × 10 08 9.7700 × 10 06
40 4.3832 × 10 13 3.6700 × 10 08 6.2600 × 10 06
48 4.3188 × 10 13 1.8000 × 10 08 4.3500 × 10 06
Table 9. Showing the results for test problem 5.5.
Table 9. Showing the results for test problem 5.5.
NExact ResultsNMCFs ResultsAEs in NMCFs
16 0.21242359899585922 0.2124235991762159 1.8036 × 10 10
32 0.11276006163890931 0.1127600617361654 9.7256 × 10 11
64 0.06256009950548104 0.06256009955588866 5.0408 × 10 11
128 0.24947430690102754 0.24947430707441792 1.7339 × 10 12
Table 10. Comparison of errors in NMCFs with Volkov et al. [48] and Biala and Jator [46] for test problem 5.5.
Table 10. Comparison of errors in NMCFs with Volkov et al. [48] and Biala and Jator [46] for test problem 5.5.
NAEs in NMCFsAEs in Biala and Jator [46]AEs in Volkov et al. [48]
16 1.8036 × 10 10 1.2860 × 10 03 3.2660 × 10 02
32 9.7256 × 10 11 3.1530 × 10 04 8.2100 × 10 03
64 5.0408 × 10 11 7.9130 × 10 05 2.053 × 10 03
128 1.7339 × 10 12 1.9750 × 10 05 5.1280 × 10 04
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Omole, E.O.; Adeyefa, E.O.; Ayodele, V.I.; Shokri, A.; Wang, Y. Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies. Axioms 2023, 12, 891. https://doi.org/10.3390/axioms12090891

AMA Style

Omole EO, Adeyefa EO, Ayodele VI, Shokri A, Wang Y. Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies. Axioms. 2023; 12(9):891. https://doi.org/10.3390/axioms12090891

Chicago/Turabian Style

Omole, Ezekiel Olaoluwa, Emmanuel Oluseye Adeyefa, Victoria Iyadunni Ayodele, Ali Shokri, and Yuanheng Wang. 2023. "Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies" Axioms 12, no. 9: 891. https://doi.org/10.3390/axioms12090891

APA Style

Omole, E. O., Adeyefa, E. O., Ayodele, V. I., Shokri, A., & Wang, Y. (2023). Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies. Axioms, 12(9), 891. https://doi.org/10.3390/axioms12090891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop