1. Introduction
Recent research in the field of hydrodynamic and hydromagnetics stability have found the presence of a family of problems in differential equations (DEs) of a high order, and which have real mathematical interest. There are various approximate (numerical) methods in the literature that have been used for the solution of boundary value problems (BVPs). The existence and uniqueness to finding the solution of higher order BVPs are systematically examined in [
1]. The BVPs of higher order DEs have been examined due to their significance and the potential for applications in applied sciences. To find the analytical solutions of such BVPs analytically is very tough and are available in very few cases. Very few researchers have tried the numerical solution of 10
th order BVPs. Some of the approximate techniques have been established over the years to the numerical solution for these kinds of BVPs. In [
2,
3], the authors has solved 10
th and 12
th order BVPs using the Adomian decomposition method (ADM) involving Green’s function. The homotopy perturbation approach was utilized in [
4] to solve BVPs of 10
th order. When a uniform magnetic field is applied across the fluid in the direction of gravity, the instability sets now as ordinary convection and it is modeled by 10
th order BVPs as discussed in [
5]. In [
6], established approximate techniques for solving the 10
th order non-linear BVPs occurring in thermal instability.
Numerical methods for the solution of non-linear BVPs of order 2 m were found in [
7]. An effective numerical procedure DTM for solving some linear and non-linear BVPs of 10
th order is discussed in [
8]. In [
9,
10], the BVPs of 9
th and 10
th order are considered by adopting homotopy perturbation technique and the modified-variational iteration technique. Also the variational iterative technique was adopted in [
11] for solving the 10
th order BVPs. Wazwaz [
12,
13,
14,
15] proposed modified form of ADM for solving 6
th, 8
th, 10
th and 12
th order.
The study of non-polynomial spline [
16] of 11
th degree is a key element to solve 10
th order BVPs. In [
17], it is depicted that the DEs that describe the 10
th order model to incorporate a 3rd order model of enlistment machine, two equations for dynamic power control, two equations for receptive power control, and three equations for edge pitch control. A 10
th order nonlinear dynamic model was developed in [
18] to turn mobile robots that incorporate slip between the driven wheels and the ground. Based on binary six-point and eight-point approximating subdivision scheme, two collocation algorithms are constructed by [
19,
20] to find the solution of BVPs. The 4
th order linear BVPs using a new cubic B-spline were solved in [
21]. Authors explained the 10
th and 12
th order BVPs by using the Galerkin weighted residual technique in [
22]. The 5
th, 6
th and 8
th order linear and non-linear BVPs by using the cubic B-spline scheme (CBS) method were solved in [
23,
24,
25]. The higher (10
th and 11
th) degree splines were tested in [
26,
27] for solving 10
th order BVPs. In [
28] they practiced 2nd order finite difference schemes for the mathematical solutions of the 8
th, 10
th and 12
th order Eigen-value problems. Galerkin method with septic B-spline and quintic B-spline was adopted in [
29,
30] for solving 10
th order BVPs. Quintic B-spline and septic-B spline collocation methods was discussed in [
31,
32] to find solution of a 10
th order BVPs.
For discrete methods, e.g., Adomian decomposition, shooting, homotopy perturbation, finite differences and variational-iterative technique, only give discrete approximate values of the unknown . For fitting curve to data we require further data processing methods. To overcome these disadvantages, we introduced a new CBS scheme for the solution of 10th order BVPs. The algorithm developed here is not only for the approximation solutions of the 10th order boundary value problems(BVPs) employing CBS but also estimate derivatives of 1st order to 10th order (where boundary conditions (BCs) are defined) of the exact solution as well.
The rest of the paper is organized as follows. The construction of CBS is presented in
Section 2. In
Section 3, the CBS scheme is utilized as an interpolating function in the solution of 10
th order nonlinear BVPs. The results and discussion are presented in
Section 4. Also some problems are considered in this section to show the efficiency of the CBS scheme. Finally, the concluding remarks are given in the final section.
2. The Construction of CBS
In this section, we construct the CBS basis functions for solving numerically the non-linear equations arising in the study of hydrodynamics and hydro-magnetic stability problems. To find the approximate solution at nodal points defined in the region
. For an interval
=
, we divide it into
n sub-intervals
=
by the equidistant knots. For this range, we select equidistant points such that
such that
i.e.,
and
.
Assume
such that
converted to to cubic-polynomial on separately sub interval
. The basis function is defined as
for
. Considering one and all
is also a piece-wise cubic with knots at
, simultaneously
.
Assume
be linearly independent and let
Thus
is
dimensional and
=
. Let
be the cubic-B spline function interpolating at the nodal points and
. Then
can be written as
Consequently now for a function
, there happened to be a distinctive cubic-B spline
, satisfying the interpolating conditions:
for
.
The values of
, and its derivatives
,
at nodal points are required and these derivatives are tabulated in
Table 1.
Assume
and
then from
may be used to determine numerical-difference formulas for
such that
, for
such that
, for
,
such that
and
,
such that
like so the errors can be obtained by using Taylor-series
3. The 10th Order Nonlinear BVPs
In this section, we consider the 10
th order nonlinear BVPs arising in the study of hydrodynamics stability and visco-elastic flows.
with BCs
where
,
,
,
,
and
,
,
,
,
are given real constants,
and
f is continuous in interval
.
The Taylor, series for
at the preferred collocation points alongside central difference (see [
31]), we have
Equation (
9) can be written as
Substituting Equation (
13) into Equation (
12), we obtain
Since
so, Equation (
14) becomes
After some simplifications the above equation becomes
Let
be the accurate solution of non-linear 10
th order BVPs
Imposing Equations (
3), (
5), (
7), (
8) and (
9) into Equation (
17), we have
Equation (
18) we will produce a new system consisting of
linear equations
with
unknowns
where
, therefore ten further equations are required. From given BCs at
, we have five equations:
similarly from
there will be other five equations
Omitting the order of the error of terms, the exact solution
is accomplished by finding solution of the discussed above linear system of
equations in
unknowns considering the Equations (
18)–(
20).
4. Convergence Analysis
Let
be the exact solution of the Equations (
10)–(
12) and also
be the CBS approximation to
. Therefore, we have
where
Also, we have assume that
be the computed cubic B spline approximation to
, namely
To approximate the error we have to estimate error and seperately
The system of
matrix can be written as:
Now, by subtracting Equations (
22) and (
23), we obtain
where
B is an
-dimensional band matrix,
and
where
T denoting transpose.
Taking the infinity norm from Equation (
24), we obtain
The B-spline
=
satisfy the following property
Using Equations (
26) and (
27) in Equation (
28)
which proves that this method is second order convergent and