A Robust and Effective Method for Solving Two-Point BVP in Modelling Viscoelastic Flows ()
1. Introduction
In modelling viscoelastic flows, differential equations of elliptic-hyperbolic operator types arise. Simulations of such flows have been studied extensively lately. The main characteristics of such elliptic-hyperbolic operators can be captured in a nonlinear fifth order two-point boundary value problem in one dimension [1] [2]. The elasticity parameter in the problem is of major importance in the investigation; that is, if it is increased, then depending on the formulation of the problem.
Many researchers have discussed solutions of viscoelastic flow model, like the Galerkin method discussed in [1], collocation method discussed in [2], the Rung-Kutta discussed in [3], the homotopy method discussed in [4]. Hermitian finite elements in [5] and the shooting method in [6] [7] [8]. Some authors have discussed other methods of solution of this model in [9] [10] [11].
In this paper, we are going to introduce Chebyshev-collocation method for the numerical solution of the fifth-order non-linear two-point boundary value problem in modelling viscoelastic flow:
(1)
concerned to the posterior boundary conditions
(2)
with regard to the given positive constants
and c which act as elasticity parameter and a boundary stress, respectively. Moreover, c is equal unity in this paper.
In recent years, a lot of attention has been devoted to the study of Chebyshev methods to investigate various scientific models. Using these methods made it possible to solve Troeschs problem [12], twelfth-order boundary-value problems [13], high-order nonlinear ordinary differential equations [14], linear integro-differential equations [15], fourth-order Sturm-Liouville problems [16], Generalized Sturm-Liouville problems [17], the parabolic inverse problem [18], two-dimensional heat equation [19], fractional diffusion equation [20], elliptic partial differential equations [21], integral and integro-differential equations of the third kind [22], the constant mobility Cahn-Hilliard equation in a square domain [23], Lane-Emden problem [24]. Recently, has been made numerical comparison of sinc-collocation and Chebychev-collocation methods for determining the eigenvalues of Sturm-Liouville problems with parameter-dependent boundary conditions by El-Gamel [25].
Chebyshev methods for ordinary differential equations have many salient features due to the properties of the basis functions and the manner in which the problem is discretized. The approximating discrete system depends only on parameters of the differential equation. There are many advantages of using Chebyshev polynomials as expansion function presented in that are good representations of smooth functions. What we do here is to seek a special Chebyshev solution which also satisfies the given boundary conditions.
We organize our paper as follows. In Section 2, we present the preliminaries which we used in this paper. Method of the solution is given in Section 3. Some numerical results are presented in Section 4 to show the efficiency of the proposed method. Finally, we draw some conclusions and closing remarks.
2. Essential Relations
Chebyshev polynomial formula of the first kind of degree m is chiefly defined and bounded in interval
see [26]
or, in didactic organism,
As for the shifted Chebyshev polynomial
of the first kind on interval
.
With leading coefficient is equal to
, by Compensation in the previous equation
In addition, the definition of the collocation points is worded as follow
(3)
Moreover, the relation between Chebyshev coefficient matrix
and
in the interval
is
where
result from the characteristics of Chebyshev polynomial.
All in all, the use of half interval
is more favored in modelling viscoelastic flows. The shifted Chebyshev polynomials can also be worded as follows
This is deduced from definition of the collocation points:
(4)
Similarly, the relation between Chebyshev coefficient matrix
and
in the interval
is
3. Method of Solution
Let is assume the approximate solution
of the main problem (1) and its derivatives is
(5)
(6)
where N is chosen as any positive integer such that
. Besides, the anonymous Chebychev coefficients of
and its derivatives are
and
, respectively. The approximate solution and its derivatives in the matrix format are
(7)
whereas the definitions of ultimate matrices are:
and
(8)
for odd N.
(9)
for even N. We need the following lemma where
and K are both positive integer.
Lemma 1. [13] [14] The following relation holds
(10)
where
We need the following the theorem:
Theorem 2. If the considered approximate solution of the problem (1) is (7), so that the discrete Chebyshev system is availed by
(11)
where
Proof. Replacing each term of (1) with the approximation defined in (7) and (10), and applying the collocation points (4) to it.
The matrices for the boundary conditions subjected to Equation (2) are
(12)
Thus, in the matrix
we will replace 5th rows by the Equation (12), we have the augmented matrix
(13)
Now we will solve a linear system (13) of
equations of the
unknown coefficients
,
. So as to gain the coefficients of the approximate solution
by the Q-R method.
Algorithm
· Input (integer) N.
· Input (double) tol.
· Input (array)
(Initial approximation,
with
dimension, can be chosen so that the boundary conditions are satisfied).
·
is a linear algebraic equation system. Then solve this system to find
.
· If
then
, break (the program is finished).
· Else then
.
4. Examples and Comparison
In this section we give an illustrative example to authenticate the obtained results on Equation (1). The performance of Chebyshev method is measured by the root mean square errors
which is defined as
Example 1: [1] [2] [3] Consider the following fifth order nonlinear two-point BVP
whose exact solution is
Table 1 exhibits a comparison between the root mean square errors obtained by using the present method, Beam function using Galerkin, Beam function using Collocation, Runge-Kutta methods and Chebyfun for Example 1. This comparison shows the strength of the first scheme. Figure 1 demonstrates Chebyshev approximate solution versus the exact solution.
Example 2: [1] [2] [3] Consider the following fifth order nonlinear two-point BVP
whose exact solution is
Table 2 represents a comparison between the root mean square errors with the methods in [1] [2] [3]. Table 3 shows a comparison of between the root mean square errors of Chebyshev and Chebyfun at
and
for Example 2. The graph of Chebyshev approximate solution and the exact solution of have been plotted in Figure 2.
Table 1. Comparison of the root mean square errors at
for Example 1.
Table 2. Comparison of the root mean square errors for Example 2, for the present method, Galerkin, collocation methods and Runge-Kutta method.
Table 3. Comparison of the root mean square errors at
and
for Example 2.
Figure 1. The exact solution and approximation solution for Example 1.
Figure 2. The exact solution and approximation solution for Example 2.
Example 3: [1] [2] [3] Consider the following fifth order nonlinear two-point BVP
whose exact solution is
and
Tables 4-7 illustrate the comparison between result of Chebyshev polynomial method and result of methods in [1] and [2] at
and 103 with N from 6 to 16. Besides, Table 8 depicts
with values 10−1, 10−2 and 10−3 that has an impact on the result of Chebyshev polynomial method at
, Chebfun matlab program and the Runge-Kutta in [3] at
. Figure 3 demonstrates Chebyshev approximate solution versus the exact solution.
Table 4. Comparison of the root mean square errors at
for Example 3.
Table 5. Comparison of the root mean square errors at
for Example 3.
Table 6. Comparison of the root mean square errors at
for Example 3.
Table 7. Comparison of the root mean square errors at
for Example 3.
Table 8. Comparison of the root mean square errors for Chebyshev polynomial method and the Runge-Kutta method for Example 3.
Figure 3. The exact solution and approximation solution for Example 3.
Example 4: Now we turn to a singular problem
and
whose exact solution is
Applying L’Hospital rule to BVP in order to remove singularity. The latter form is
The computational results are summarized in Table 9.
Table 10 exhibits a comparison between the root mean square errors at
and
for Example 4 obtained by using Chebyshev method and Chebfun. Figure 4 demonstrates Chebyshev approximate solution and the exact which appears to be in good agreement with each other.
5. Conclusion
In this article, we present a method to approximate the solution of viscoelastic flows. The numerical method is based on the operational matrix of Chebychev polynomials. We present four examples, the first three examples are nonlinear fifth-order differential equation and fourth example is singular nonlinear. We compared the results of this algorithm with others and showed the accuracy and potential applicability of the given method. The proposed method is a powerful tool for obtaining novel numerical solutions of such equations. It is advisable to use it for other nonlinear differential equations.
Table 9. Comparison of the root mean square errors at
for Example 4.
Table 10. Comparison of the root mean square errors at
and
for Example 4.
Figure 4. The exact solution and approximation solution for Example 4.