Next Article in Journal / Special Issue
Analytical and Computational Problems Related to Fractional Gaussian Noise
Previous Article in Journal
Analysis of a Hybrid Coupled System of ψ-Caputo Fractional Derivatives with Generalized Slit-Strips-Type Integral Boundary Conditions and Impulses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modal Shifted Fifth-Kind Chebyshev Tau Integral Approach for Solving Heat Conduction Equation

by
Ahmed Gamal Atta
1,
Waleed Mohamed Abd-Elhameed
2,
Galal Mahrous Moatimid
1 and
Youssri Hassan Youssri
2,*
1
Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo 11341, Egypt
2
Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(11), 619; https://doi.org/10.3390/fractalfract6110619
Submission received: 1 October 2022 / Revised: 18 October 2022 / Accepted: 20 October 2022 / Published: 22 October 2022
(This article belongs to the Special Issue Feature Papers for Numerical and Computational Methods Section)

Abstract

:
In this study, a spectral tau solution to the heat conduction equation is introduced. As basis functions, the orthogonal polynomials, namely, the shifted fifth-kind Chebyshev polynomials (5CPs), are used. The proposed method’s derivation is based on solving the integral equation that corresponds to the original problem. The tau approach and some theoretical findings serve to transform the problem with its underlying conditions into a suitable system of equations that can be successfully solved by the Gaussian elimination method. For the applicability and precision of our suggested algorithm, some numerical examples are given.

1. Introduction

The heat equation pioneered by Fourier [1] describes the distribution of heat in a given body over time [2], which is a type of second-order parabolic partial differential equation. It has many applications in diverse scientific fields. Moreover, it has been studied analytically and numerically. For example, Meyu and Koriche [3] proposed two techniques based on the separation of variables and finite-difference methods to solve the heat equation in one dimension. Liu and Chang [4] used a method of nonlocal boundary shape functions to solve a nonlinear heat equation with nonlocal boundary conditions. Tassaddiq et al. [5] introduced an approximate approach based on a cubic B-spline collocation method to solve the heat equation with classical and nonclassical boundary conditions.
It is well-known that obtaining accurate and efficient methods for solving differential equations has become an important research point. There are several analytic and numerical methods, such as the homotopy analysis method [6,7], the variational iteration method [8,9], the Adomian decomposition method [10,11], the finite-difference method [12,13,14], the finite-element method [15,16,17], and spectral methods [18,19,20,21,22]. Spectral methods have many advantages if compared with the other methods because they yield exponential rate convergence, a good accuracy, and the computational efficiency of the solutions while failing for many complicated problems with singular solutions. Thus, it is relevant to be interested in how to enlarge the adaptability of spectral methods and construct certain simple approximation schemes without a loss of accuracy for more complicated problems. Further applications of spectral methods in different disciplines may be found in [23,24,25,26,27,28,29].
Orthogonal polynomials, such as Legendre polynomials and Chebyshev polynomials, have received a lot of attention from both theoretical and practical perspectives [30,31,32]. Chebyshev polynomials have been used as an important category of basis functions to solve ordinary, partial, and fractional differential equations, see for instance [33,34,35,36,37,38]. Two major reasons for the widespread use of these polynomials are the high accuracy of the approximation and the simplicity of numerical methods established based on these polynomials. There are six types of Chebyshev polynomials, they are Chebyshev polynomials of the first, second, third, fourth, fifth, and sixth kind. All the kinds of Chebyshev polynomials have their important parts in numerical analysis and approximation theory. There are old and recent contributions regarding the first four kinds, see, for example, [39,40], while the fifth and sixth kind of Chebyshev polynomials have gained recently a fast-growing attention from many authors. For instance, Sadri and Aminikhah in [41] treated a multiterm variable-order time-fractional diffusion-wave equation using a new efficient algorithm based on the 5 C P s . Moreover, Abd-Elhameed and Youssri in [42] employed the 5 C P s for solving the convection–diffusion equation.
The following items are the main goals of this paper:
  • Deriving new theorems, corollaries, and lemmas concerned with the shifted 5 C P s that serve in the derivation of our proposed numerical scheme.
  • Presenting a new spectral tau algorithm for the numerical treatment of the heat conduction equation.
  • Investigating the convergence analysis of the proposed double-shifted Chebyshev expansion.
  • Performing some comparisons to clarify the efficiency and accuracy of our method.
To the best of our knowledge, some advantages of the proposed technique can be mentioned as follows:
  • By choosing the shifted 5 C P s as basis functions, and taking a few terms of the retained modes, it is possible to produce approximations with excellent precision. Less calculation is required. In addition, the resulting errors are small.
  • In comparison to other Chebyshev polynomials, the shifted 5 C P s are not as well-studied or used. This motivates us to find theoretical findings concerning them. Furthermore, we found that the obtained numerical results, if they are used as basis functions, are satisfactory.
We point out here that the novelty of our contribution in this paper can be listed as follows:
  • Some derivatives and integral formulas of the shifted 5 C P s are given in reduced formulas that do not involve any hypergeometric forms.
  • The employment of these basis functions to the numerical treatment of the heat conduction equation is new.
The contents of the paper are arranged as follows. Section 2 is devoted to presenting mathematical preliminaries containing some relevant properties of 5 C P s and their shifted ones. In addition, some new formulas concerning the shifted 5 C P s are derived. In Section 3, we present and implement a spectral tau method for solving the heat conduction equation based on employing the shifted 5 C P s . In Section 4, we investigate in detail the convergence and error analysis of the suggested shifted 5 C P s . In Section 5, some numerical examples are given to ensure the efficiency, simplicity, and applicability of the suggested method. Finally, conclusions are reported in Section 6.

2. An Account on the Shifted 5 C P s and Some New Useful Formulas

This section is confined to presenting an account on the 5 C P s , X j ( t ) , j 0 , and their shifted ones. In addition, building on some of their fundamental relations, we derive some new specific formulas that serve in the derivation of our proposed numerical scheme. More precisely, we establish the second-order derivative formulas of the shifted polynomials and also the corresponding integral formulas of these polynomials.

2.1. An Account on the Shifted 5 C P s

The 5 C P s ( X i ( t ) ) are a sequence of orthogonal polynomials on [ 1 , 1 ] (see, [28,43]) that satisfy the following orthogonality relation:
1 1 t 2 1 t 2 X i ( t ) X j ( t ) d t = h i δ i , j ,
where
h i = π 2 2 i + 1 1 , i f i even , i + 2 i , i f i odd ,
and
δ i , j = 1 , i f i = j , 0 , i f i j .
X i ( t ) may be generated with the aid of the following recursive formula:
X i ( t ) = t X i 1 ( t ) ϱ i X i 2 ( t ) , i 2 ,
where X 0 ( t ) = 1 , X 1 ( t ) = t and
ϱ i = ( i 1 ) 2 + i + ( 1 ) i ( 2 i 1 ) 4 i ( i 1 ) .
The shifted orthogonal 5 C P s on [ 0 , τ ] are defined as
C i ( t ) = X i 2 t τ 1 , τ > 0 ,
with the orthogonality relation
0 τ ω ( t ) C i ( t ) C j ( t ) d t = h τ , i δ i , j ,
where h τ , i = τ 2 h i and ω ( t ) = ( 2 t τ ) 2 t τ t 2 .
Lemma 1
([28]). The analytic formula of X j ( x ) may be split to the following two analytic formulas:
X 2 i ( x ) = ( 2 i + 1 ) n = 0 i ( 1 ) n 2 2 n ( 2 i n ) ! n ! ( 2 i 2 n + 1 ) ! x 2 i 2 n , i 0 ,
X 2 i + 1 ( x ) = i + 3 2 ! ( 2 i + 1 ) ! n = 0 i ( 1 ) n i i n ( 2 i n + 1 ) ! i n + 3 2 ! x 2 i 2 n + 1 , i 0 .
Theorem 1
([28]). The following two inversion formulas hold for the polynomials X j ( x ) :
x 2 j = r = 0 j 2 j + 1 r 2 2 r X 2 j 2 r ( x ) , j 0 ,
x 2 j + 1 = j + 3 2 ! r = 0 j ( j j r ) ( 2 j 2 r + 2 ) ! j r + 3 2 ! ( 2 j r + 2 ) ! X 2 j 2 r + 1 ( x ) , j 0 .

2.2. Derivation of the Second-Order Derivative Formulas of C j ( t )

The following theorem exhibits the expressions of the second-order derivatives of X j ( t ) in terms of their original ones.
Theorem 2.
The second-order derivative of the polynomials X i ( t ) can be expressed explicitly as:
D 2 X 2 i ( t ) = n = 0 i 4 n ( 1 ) n ( 2 i + 2 n + 1 ) ( 2 i + 1 ) 2 i ( n + 1 ) + n 2 + n + 1 X 2 i 2 n 2 ( t ) ,
D 2 X 2 i + 1 ( t ) = ( 2 i + 3 ) n = 0 i 2 1 2 n ( n i ) ( 1 ) n 4 ( i n ) 2 1 4 i ( i + 2 ) + 3 2 i ( n + 1 ) + n 2 2 i 2 n + 1 X 2 i 2 n 1 ( t ) .
Proof. 
First, we prove relation (6). The power-form representation of X i ( t ) in (2) enables one to express D 2 X 2 i ( t ) , in the following form:
D 2 X 2 i ( t ) = ( 2 i + 1 ) n = 0 i ( 1 ) n 2 1 2 n ( 2 i 2 n 1 ) ( i n ) ( 2 i n ) ! n ! ( 2 i 2 n + 1 ) ! t 2 i 2 n 2 ,
which can be written with the aid of the inversion formula (4) as
D 2 X 2 i ( t ) = ( 2 i + 1 ) n = 0 i ( 1 ) n 2 1 2 n ( 2 i 2 n 1 ) ( i n ) ( 2 i n ) ! n ! ( 2 i 2 n + 1 ) ! × r = 0 i n 1 2 i 2 n 1 r 2 2 r X 2 i 2 n 2 r 2 ( t ) .
The last relation after expanding and rearranging the terms can be converted into
D 2 X 2 i ( t ) = ( 2 i + 1 ) n = 0 i 4 n × r = 0 n ( 1 ) n r + 1 ( 2 i + 2 n 2 r + 1 ) ( 2 i n + r ) ! ( r ) ! ( 2 i 2 n + 2 r + 1 ) ( n r ) ! ( 2 i 2 n + r 1 ) ! X 2 i 2 n 2 ( t ) .
Now, in order to reduce the summation on the right-hand side of the last formula, set
G n , i = r = 0 n ( 1 ) n r + 1 ( 2 i + 2 n 2 r + 1 ) ( 2 i n + r ) ! ( r ) ! ( 2 i 2 n + 2 r + 1 ) ( n r ) ! ( 2 i 2 n + r 1 ) ! .
The application of Zeilberger’s algorithm mentioned in [44] enables us to get the following recurrence relation for G n , i :
4 i 2 n 6 i n 2 + 2 n 3 + 10 i 2 26 i n + 12 n 2 29 i + 25 n + 18 G n , i + 4 i 2 4 i n + 2 n 2 4 i + 6 n + 1 G n + 1 , i = 0 + 4 i 2 n + 6 i n 2 2 n 3 6 i 2 + 14 i n 6 n 2 + 9 i 7 n 3 G n + 2 , i = 0 ,
with the initial values:
G 0 , i = 2 i ( 2 i 1 ) 2 i + 1 , G 1 , i = 2 ( 4 i 2 3 ) 2 i + 1 .
The recurrence relation (8) can be exactly solved to give
G n , i = 2 i ( n + 1 ) + ( 1 ) n ( 2 i + 2 n + 1 ) 2 i + 1 n 2 n 1 ,
and therefore, relation (6) can be obtained.
Now, we prove Formula (7). Based on relation (3), we have
D 2 X 2 i + 1 ( t ) = n = 0 i ( 2 i + 3 ) ( 1 ) n + 1 4 n ( 2 i n + 1 ) ! n ! ( 2 i + 2 n 3 ) ( 2 i 2 n 1 ) ! t 2 i 2 n 1 .
Making use of Formula (5) yields
D 2 X 2 i + 1 ( t ) = n = 0 i ( 2 i + 3 ) ( 1 ) n + 1 4 n ( 2 i n + 1 ) ! i n + 1 2 ! n ! ( 2 i + 2 n 3 ) ( 2 i 2 n 1 ) ! × r = 0 i n 1 ( 2 i 2 n 2 r ) ! i n 1 i n r 1 ( 2 i 2 n r ) ! i n r + 1 2 ! X 2 i 2 n 2 r 1 ( t ) .
The last relation after expanding and rearranging the terms can be converted into
D 2 X 2 i + 1 ( t ) = n = 0 i ( 2 i + 3 ) 2 1 2 n ( n i ) 2 i 2 n + 1 × r = 0 n ( 1 ) n r ( 2 i + 2 n 2 r 1 ) ( 2 i n + r + 1 ) ! ( r ) ! ( 2 i 2 n + 2 r + 3 ) ( n r ) ! ( 2 i 2 n + r ) ! X 2 i 2 n 1 ( t ) .
Now, set
G ˜ n , i = r = 0 n ( 1 ) n r ( 2 i + 2 n 2 r 1 ) ( 2 i n + r + 1 ) ! ( r ) ! ( 2 i 2 n + 2 r + 3 ) ( n r ) ! ( 2 i 2 n + r ) ! ,
and utilize again Zeilberger’s algorithm to show that G ˜ n , i satisfies the following recurrence relation:
4 i n 2 n 2 + 6 i 2 n + 1 1 + 2 i 2 n 2 G ˜ n + 2 , i 4 4 i 2 + 8 i 1 1 + i n G ˜ n + 1 , i 4 i n 2 n 2 + 10 i 6 n 3 3 + 2 i 2 n 2 G ˜ n , i = 0 ,
with the initial values:
G ˜ 0 , i = ( 1 + 2 i ) 2 2 i + 3 , G ˜ 1 , i = 4 i ( 4 i ( i + 2 ) 1 ) 4 i ( i + 2 ) + 3 .
The recurrence relation (9) can be exactly solved to give
G ˜ n , i = ( 1 ) n 4 ( i n ) 2 1 4 i ( i + 2 ) + 3 2 i ( n + 1 ) + n 2 ,
and therefore, relation (7) can be obtained. ◻
As a result of Theorem 2, the formula expressing the derivatives of the 5 C P s can be merged to give the following result.
Corollary 1.
Let j 2 . The second-order derivative of the polynomials X j ( t ) can be expressed explicitly as:
D 2 X j ( t ) = r = 0 j 2 λ r , j X r ( t ) ,
where
λ r , j = 2 r j ( j + r + 2 ) j 2 ( j + 1 ) r + j 4 , if j even a n d j r 2 2 even , ( j r ) ( j ( j + r + 3 ) + r 2 ) , if j even a n d j r 4 2 even , ( r + 1 ) j j 2 ( j + 4 ) 8 + ( j ( j + 2 ) + 4 ) r 2 + 2 ( j ( j + 2 ) + 4 ) r j ( r + 2 ) if j odd a n d j r 2 2 even , ( j ( j + 2 ) 4 ) ( r + 1 ) ( j r ) ( j + r + 2 ) j ( r + 2 ) , if j odd a n d j r 4 2 even , 0 , o t h e r w i s e .
Now, the second-order derivatives of the shifted polynomials C j ( t ) can be easily deduced. The following corollary exhibits this result.
Corollary 2.
Let j 2 . The second-order derivative of the polynomials C j ( t ) can be expressed explicitly as:
D 2 C j ( t ) = r = 0 j 2 λ ¯ r , j C r ( t ) ,
where
λ ¯ r , j = 2 r j + 2 τ 2 ( j + r + 2 ) j 2 ( j + 1 ) r + j 4 , if j even a n d j r 2 2 even , ( j r ) ( j ( j + r + 3 ) + r 2 ) , if j even a n d j r 4 2 even , ( r + 1 ) j j 2 ( j + 4 ) 8 + ( j ( j + 2 ) + 4 ) r 2 + 2 ( j ( j + 2 ) + 4 ) r j ( r + 2 ) if j odd a n d j r 2 2 even , ( j ( j + 2 ) 4 ) ( r + 1 ) ( j r ) ( j + r + 2 ) j ( r + 2 ) , if j odd a n d j r 4 2 even , 0 , o t h e r w i s e .
Proof. 
The result is a direct consequence of Corollary 1 by replacing t by 2 t τ 1 . ◻

2.3. Derivation of Integral Formulas of C j ( x )

In this section, new integral formulas of C j ( t ) are derived in detail. For this derivation, the following two lemmas are useful.
Lemma 1.
Let i > 0 and 1 < n < i + 2 . One has
3 F 2 1 n , i n + 3 2 , 2 + 2 i n 4 + 2 i 2 n , i n + 5 2 1 = 2 ( 1 ) n ( n 1 ) ! ( 2 i + 1 ) ( 2 i + n 1 ) n 2 ,
where ( a ) n = Γ ( a + n ) Γ ( a ) .
Lemma 2.
Let i > 0 and 1 < n < i + 2 . One has
3 F 2 n , i n + 1 2 , 1 + 2 i n 3 + 2 i 2 n , i n + 3 2 1 = 4 i ( 1 ) n ( n ) ! i ( 2 i + 1 ) ( 2 i 2 n + 3 ) ( n 2 i ) n 2 .
Proof. 
The proofs of Lemmas 1 and 2 can be done through some algebraic manipulations along with Zeilberger’s algorithm [44]. ◻
Theorem 3.
For all i 0 , the following integral formulas hold:
0 t X 2 i ( z ) d z = n = 2 i 2 3 2 n ( 2 i 2 n + 1 ) n ( 2 i 2 n + 1 ) 2 ( 2 i 2 n + 3 ) 2 ( n 2 i ) n 2 X 2 i 2 n + 1 ( t ) + ( 2 i 1 ) ( 4 i ( i + 3 ) 1 ) 4 1 4 i 2 2 X 2 i 1 ( t ) + 1 2 i + 1 X 2 i + 1 ( t ) ,
0 t X 2 i + 1 ( z ) d z = ( 2 i + 3 ) 4 ( i + 1 ) 2 2 i + 3 X 2 i + 2 ( t ) 1 2 ( 2 i + 1 ) X 2 i ( t ) + 1 4 i .
Proof. 
We prove formula (11). The power-form representation (2) enables one to express 0 t X 2 i ( z ) d z as
0 t X 2 i ( z ) d z = ( 2 i + 1 ) n = 0 i ( 1 ) n ( 2 i n ) ! 2 2 n n ! ( 2 i 2 n + 1 ) ( 2 i 2 n + 1 ) ! t 2 i 2 n + 1 .
In virtue of relation (5), the last equation may be written alternatively as
0 t X 2 i ( z ) d z = ( 2 i + 1 ) n = 0 i ( 1 ) n ( 2 i n ) ! i n + 3 2 ! 2 2 n n ! ( 2 i 2 n + 1 ) ( 2 i 2 n + 1 ) ! × r = 0 i n ( 2 i 2 n 2 r + 2 ) ! i n i n r ( 2 i 2 n r + 2 ) ! i n r + 3 2 ! X 2 i 2 n 2 r + 1 ( t ) .
After rearranging and expanding the terms in the previous equation, one gets
0 t X 2 i ( z ) d z = ( 2 i + 1 ) n = 0 i 1 4 n ( 2 i n ) ! ( 2 i 2 n + 1 ) 2 ( 2 i + 2 n 3 ) 3 ( n ) ! ( 2 i 2 n ) ! × ( 2 i + 2 n 3 ) 3 3 F 2 n , i n + 1 2 , 1 + 2 i n 3 + 2 i 2 n , i n + 3 2 1 + 2 n ( 2 i + n 1 ) ( 2 i + 2 n 1 ) 3 F 2 1 n , i n + 3 2 , 2 + 2 i n 4 + 2 i 2 n , i n + 5 2 1 X 2 i 2 n + 1 ( t ) .
Thanks to Lemmas 1 and 2, we get the desired relation (11).
Relation (12) can be similarly proved through some algebraic computations. ◻
The following corollary is a direct consequence of Theorem 3.
Corollary 3.
For all i 0 , the following integrals formulas hold
0 t C 2 i ( z ) d z = n = 2 i τ 2 2 2 n ( 2 i 2 n + 1 ) n ( 2 i 2 n + 1 ) 2 ( 2 i 2 n + 3 ) 2 ( n 2 i ) n 2 C 2 i 2 n + 1 ( t ) + τ ( 2 i 1 ) ( 4 i ( i + 3 ) 1 ) 8 1 4 i 2 2 C 2 i 1 ( t ) + τ 2 ( 2 i + 1 ) C 2 i + 1 ( t ) + c ¯ ,
0 t C 2 i + 1 ( z ) d z = ( 2 i + 3 ) τ 8 ( i + 1 ) 2 2 i + 3 C 2 i + 2 ( t ) 1 2 ( 2 i + 1 ) C 2 i ( t ) + 1 4 i + c ¯ ¯ ,
where c ¯ and c ¯ ¯ are constants.

3. A Numerical Tau Approach for the Treatment of the Heat Conduction Equation

This section focuses on obtaining a new spectral solution to the heat conduction equation subject to an initial condition and homogeneous or nonhomogeneous boundary conditions with the aid of the spectral tau method.

Treatment of the Equation Subject to Homogeneous Boundary Conditions

Now, consider the following heat conduction equation [45]:
u ( x , t ) t = k 2 u ( x , t ) x 2 + f ( x , t ) , 0 x , 0 t τ ,
governed by the initial condition:
u ( x , 0 ) = ψ ( x ) ,
and by the homogeneous boundary conditions:
u ( 0 , t ) = 0 , u ( , t ) = 0 ,
where f ( x , t ) represents the source term and k is a real constant.
If we integrate Equation (13) with respect to t , then the following equation is obtained:
u ( x , t ) = k 0 t 2 u ( x , z ) x 2 d z + f ¯ ( x , t ) , 0 x , 0 t τ ,
subject to the homogeneous boundary conditions:
u ( 0 , t ) = 0 , u ( , t ) = 0 ,
where
f ¯ ( x , t ) = 0 t f ( x , z ) d z + ψ ( x ) ,
then, we can alternatively solve Equation (15) instead of Equation (13).
Now, define
P N = span C i ( x ) C j ( t ) : i , j = 0 , 1 , . . . , N ,
then, any function u ( x , t ) L ω ^ ( x , t ) 2 can be approximated by the truncated double series
u N ( x , t ) = i = 0 N j = 0 N a i j C i ( x ) C j ( t ) = C ( x ) A C T ( t ) ,
where
C ( x ) = [ C 0 ( x ) , C 1 ( x ) , . . . , C N ( x ) ] ,
A = ( a i j ) 0 i , j N is a matrix of order ( N + 1 ) × ( N + 1 ) ,   I = [ 0 , ] × [ 0 , τ ] and ω ^ ( x , t ) = ( 2 x ) 2 x x 2 ( 2 t τ ) 2 t τ t 2 .
Now, the application of the tau method implies that
u N ( x , t ) , C m ( x ) C n ( t ) ω ^ ( x , t ) k 0 t 2 u N ( x , z ) x 2 d z , C m ( x ) C n ( t ) ω ^ ( x , t ) = f ¯ ( x , t ) , C m ( x ) C n ( t ) ω ^ ( x , t ) , 1 m , n N .
Let us denote
B = ( b m i ) N × ( N + 1 ) , b m i = ( C i ( x ) , C m ( x ) ) ω ( x ) , D = d m i N × ( N + 1 ) , d m i = ( C i ( x ) , C m ( x ) ) ω ( x ) , G = ( g n j ) N × ( N + 1 ) , g n j = 0 t C j ( z ) d z , C n ( t ) ω ( t ) , F = ( f m n ) N × N , f r , s = ( f ¯ ( x , t ) , C m ( x ) C n ( t ) ) ω ^ ( x , t ) ,
where
( u ( x ) , v ( x ) ) ω ( x ) = 0 u ( x ) v ( x ) ω ( x ) d x
and
( u ( x , t ) , v ( x , t ) ) ω ^ ( x , t ) = 0 0 τ u ( x , t ) v ( x , t ) ω ^ ( x , t ) d t d x .
In matrix form, Equation (18) may be rewritten as
B A B T = k D A G T + F ,
where the nonzero elements of the matrices B , D , and G are given as in the next theorem. In addition, making use of the homogeneous boundary conditions (16) yields
C ( 0 ) A C T i + 1 N + 2 = 0 , i = 0 , 1 , N , C ( ) A C T i + 1 N + 2 = 0 , i = 0 , 1 , N 1 .
Equations (20) and (21) generate a system of algebraic equations of dimension ( N + 1 ) 2 in the unknown expansion coefficients a i j . Thanks to the Gaussian elimination technique, the required numerical solution can be obtained.
Theorem 4.
The elements of the matrices B , D , and G can be computed explicitly as follows:
b m i = h , i , if i = m , 0 , o t h e r w i s e ,
d m i = λ ¯ m , i h , m , if i m 2 , 0 , o t h e r w i s e ,
g n j = s = 2 j 2 M s , j 2 h τ , j 2 s + 1 δ j 2 s + 1 , n + M ¯ j 2 h τ , j 1 δ j 1 , n + M ¯ ¯ j 2 h τ , j + 1 δ j + 1 , n , if j even , τ ( 4 j h τ , j + 1 δ j + 1 , n ( j + 2 ) h τ , j 1 δ j 1 , n ) 8 j ( j + 1 ) , i f j odd ,
where
M s , j = ( 1 ) 2 s τ 2 2 2 s ( 2 j 2 s + 1 ) s ( 2 j 2 s + 1 ) 2 ( 2 j 2 s + 3 ) 2 ( s 2 j ) s 2 , M ¯ j = ( 2 j 1 ) ( 4 j ( j + 3 ) 1 ) τ 2 4 1 4 j 2 2 and M ¯ ¯ j = τ 2 ( 2 j + 1 ) .
Proof. 
First, we prove the two Formulae (22) and (23). From Equations (10) and (19), we have
b m i = ( C i ( x ) , C m ( x ) ) ω ( x ) = 0 ω ( x ) C i ( x ) C m ( x ) d x
and
d m i = D 2 C i ( x ) , C m ( x ) ω ( x ) = r = 0 i 2 λ ¯ r , i 0 ω ( x ) C r ( x ) C m ( x ) d x .
Making use of the orthogonality relation (1), the desired results (22) and (23) can be obtained.
Let us prove Formula (24). In virtue of Corollary 3 and Equation (19), we get
g n , 2 j = 0 t C 2 j ( z ) d z , C n ( t ) ω ( t ) = s = 2 j ( M s , j ( C 2 j 2 s + 1 ( t ) , C n ( t ) ) ω ( t ) + M ¯ j ( C 2 j 1 ( t ) , C n ( t ) ) ω ( t ) + M ¯ ¯ j ( C 2 j + 1 ( t ) , C n ( t ) ) ω ( t ) + ( c ¯ , C n ( t ) ) ω ( t ) )
and
g n , 2 j + 1 = ( 2 j + 3 ) τ 8 ( j + 1 ) 2 2 j + 3 ( C 2 j + 2 ( t ) , C n ( t ) ) ω ( t ) 1 2 ( 2 j + 1 ) ( C 2 j ( t ) , C n ( t ) ) ω ( t ) + ( 2 j + 3 ) 1 4 j τ 8 ( j + 1 ) + c ¯ ¯ , C n ( t ) ω ( t ) .
The application of the orthogonality relation (1) enables us to write Equations (25) and (26) as
g n j = s = 2 j 2 M s , j 2 h τ , j 2 s + 1 δ j 2 s + 1 , n + M ¯ j 2 h τ , j 1 δ j 1 , n + M ¯ ¯ j 2 h τ , j + 1 δ j + 1 , n , if j even , τ ( 4 j h τ , j + 1 δ j + 1 , n ( j + 2 ) h τ , j 1 δ j 1 , n ) 8 j ( j + 1 ) , if j odd ,
and this proves Formula (24). ◻
Remark 1.
Consider the heat conduction Equation (13) subject to the initial condition (14) and the nonhomogeneous boundary conditions
u ( 0 , t ) = h 1 ( t ) , u ( , t ) = h 2 ( t ) , 0 < t τ .
Then, the following transformation
v ( x , t ) : = u ( x , t ) 1 x h 1 ( t ) x h 2 ( t ) ,
enables us to convert the nonhomogeneous boundary conditions (27) to the homogeneous ones.

4. Convergence and Error Analysis

In this section, we study the convergence of the numerical solution (17) to the exact solution u ( x , t ) of Equation (13). We discuss the analysis of the convergence for the following two cases:
  • The case in which the solution u ( x , t ) is separable.
  • The case in which the solution u ( x , t ) is not separable.

4.1. The Case Where the Solution Is Separable

Theorem 5.
Assume that the function u ( x , t ) = g 1 ( x ) g 2 ( t ) L ω ^ ( x , t ) 2 ( I ) and assume that each of g 1 ( x ) and g 2 ( t ) has a bounded third derivative such that
u ( x , t ) = i = 0 j = 0 a i j C i ( x ) C j ( t ) .
Then, the above series (28) is uniformly convergent to u ( x , t ) , and the expansion coefficients a i j satisfy the inequality:
| a i j | 1 ( i j ) 3 2 i + j , i , j > 3 ,
where the expression X Y means that there exists a generic constant n independent of N and any function such that X n Y .
Proof. 
The orthogonality relation of C i ( x ) enables us to write the expansion coefficients a i j as
a i j = 1 h , i h τ , j 0 0 τ ω ^ ( x , t ) u ( x , t ) C i ( x ) C j ( t ) d t d x .
According to the hypotheses of the theorem, one can write
a i j = 1 h , i h τ , j 0 ω ( x ) g 1 ( x ) C i ( x ) d x 0 τ ω ( t ) g 2 ( t ) C j ( t ) d t .
If we make use of the two following substitutions
2 x 1 = cos ζ 1 , 2 t τ 1 = cos ζ 2 ,
then the last equation may be rewritten in the following form
a i j = 2 τ 2 h , i h τ , j 0 π g 1 2 ( 1 + cos ζ 1 ) X i ( cos ζ 1 ) cos 2 ζ 1 d ζ 1 × 0 π g 2 2 ( 1 + cos ζ 2 ) X j ( cos ζ 2 ) cos 2 ζ 2 d ζ 2 .
Bsed on the following trigonometric representations [43]
X i ( cos θ ) = 1 2 i cos ( i + 1 ) θ cos θ , if i even , ( i + 2 ) cos θ cos ( i + 1 ) θ cos ( i + 2 ) θ i cos 2 θ , if i odd ,
along with the assumptions that g 1 ( x ) and g 2 ( t ) have a bounded third derivative and following similar steps to those followed in [43], the desired result can be obtained. ◻
Theorem 6.
The following truncation error estimate is valid
| u ( x , t ) u N ( x , t ) | 1 2 2 N .
Proof. 
Based on relations (17) and (28), we get
| u ( x , t ) u N ( x , t ) | j = N + 1 ( a 0 j C 0 ( x ) + a 1 j C 1 ( x ) + a 2 j C 2 ( x ) + a 3 j C 3 ( x ) ) C ( t ) + i = N + 1 ( a i 0 C 0 ( t ) + a i 1 C 1 ( t ) + a i 2 C 2 ( t ) + a i 3 C 3 ( t ) ) C i ( x ) + i = 4 N j = N + 1 a i j C i ( x ) C j ( t ) + i = N + 1 j = 4 a i j C i ( x ) C j ( t ) .
Following similar steps to those given in [43], we get
| a 0 j | 1 j 3 2 j , | a 1 j | 1 j 3 2 j , | a 2 j | 1 j 3 2 j , | a 3 j | 1 j 3 2 j , | a i 0 | 1 i 3 2 i , | a i 1 | 1 i 3 2 i , | a i 2 | 1 i 3 2 i , | a i 3 | 1 i 3 2 i .
Inserting Equation (30) into Equation (29) and using Theorem 5 along with the identity | C n ( x ) | < n 2 n 1 , lead to the desired estimation. ◻
Lemma 3.
The following inequalities hold for the first and second derivatives of C n ( x ) :
| D C n ( x ) | 4 n 3 2 n , x [ 0 , ] , n 0 ,
and
| D 2 C n ( x ) | n 5 2 2 n , x [ 0 , ] , n 0 .
Proof. 
The proofs of (31) and (32) are similar. Now, we prove (31).
Consider the following two cases:
  • For n = 2 i , we have (see Theorem 2.4 in [28])
    C 2 i ( x ) = 2 1 2 i j = 0 i ( 1 ) i + j δ j T 2 j * ( x ) ,
    where T j * ( x ) are the shifted first-kind Chebyshev polynomials. Based on the inequality:
    | D T n * ( x ) | 2 n 2 , we get
    | D C n ( x ) | = 2 1 2 i j = 0 i ( 1 ) i + j δ 2 j D T 2 j * ( x ) 2 2 2 i j = 0 i ( 2 j ) 2 = 2 3 2 i 3 i ( i + 1 ) ( 2 i + 1 ) 4 ( 2 i ) 3 2 2 i = 4 n 3 2 n .
  • For n = 2 i + 1 , we have (see Theorem 2.4 in [28])
    C 2 i + 1 ( x ) = 2 2 i 2 i + 1 j = 0 i ( 1 ) i + j ( 2 j + 1 ) T 2 j + 1 * ( x ) .
    Now, it is easy to write
    | D C n ( x ) | 2 1 2 i ( 2 i + 1 ) j = 0 i ( 2 j + 1 ) 3 = 2 1 2 i ( 2 i + 1 ) ( i + 1 ) 2 ( 1 + 2 i ( 2 + i ) ) 4 ( 2 i + 1 ) 3 2 2 i + 1 = 4 n 3 2 n .
    The above two cases lead to the estimation
    | D C n ( x ) | 4 n 3 2 n , x [ 0 , ] , n 0 .
The Inequality in (32) can be obtained using the inequality: | D 2 T n * ( x ) | 4 n 2 ( n 2 1 ) 3 2 and imitating the previous steps. ◻
Lemma 4.
Let u ( x , t ) and u N ( x , t ) satisfy the assumptions of Theorem 5. One gets
( u ( x , t ) u N ( x , t ) ) t N 2 N 2 ,
and
2 ( u ( x , t ) u N ( x , t ) ) x 2 N 2 2 N 4 .
Proof. 
Based on Lemma 3 and following similar steps as in Theorem 6, we get the desired results. ◻
Theorem 7.
Assume that R N ( x , t ) is the residual of Equation (13), then | R N ( x , t ) | 0 as N .
Proof. 
R N ( x , t ) of Equation (13) can be written as
R N ( x , t ) = u N ( x , t ) t k 2 u N ( x , t ) x 2 f ( x , t ) ,
which may be written alternatively with the aid of Equation (13) as
| R N ( x , t ) | = ( u ( x , t ) u N ( x , t ) ) t k 2 ( u ( x , t ) u N ( x , t ) ) x 2 ( u ( x , t ) u N ( x , t ) ) t + | k | 2 ( u ( x , t ) u N ( x , t ) ) x 2 .
Now, the application of Lemma 4 leads us to deduce that:
( u ( x , t ) u N ( x , t ) ) t 0 and 2 ( u ( x , t ) u N ( x , t ) ) x 2 0 ,
and therefore, | R N ( x , t ) | 0 as N .  ◻

4.2. The Case Where the Solution Is Nonseparable

Here, we follow Sadri and Aminikhah [41] to introduce two theorems about the convergence of our spectral tau method in the two-dimensional Chebyshev-weighted Sobolev space:
CH ω ^ ( x , t ) m ( I ) = u : i + j u x i t j L ω ^ ( x , t ) 2 ( I ) , 0 i + j m , m N ,
endowed with the norm
| | u | | CH ω ^ ( x , t ) m ( I ) = | n | m | n | u x n 1 t n 2 L ω ^ ( x , t ) 2 2 1 2 ,
where n = ( n 1 , n 2 ) such that n 1 , n 2 Z + , | n | = n 1 + n 2 .
Theorem 8
([41]). Assume that 0 n m < N + 1 , 0 k m < N + 1 , and u N ( x , t ) is the shifted fifth-kind Chebyshev approximation of u ( x , t ) CH ω ^ ( x , t ) m ( I ) . Then, the following estimations are satisfied
n x n ( u ( x , t ) u N ( x , t ) ) L ω ^ ( x , t ) 2 3 γ 1 N 7 4 ( m n ) m u ( x , t ) t m L ω ^ ( x , t ) 2
and
k t k ( u ( x , t ) u N ( x , t ) ) L ω ^ ( x , t ) 2 3 γ 2 N 7 4 ( m k ) m u ( x , t ) t m L ω ^ ( x , t ) 2 ,
where γ 1 and γ 2 are positive constants independent of any function.
Corollary 4.
Assume that 0 m < N + 1 , and u N ( x , t ) is the shifted fifth-kind Chebyshev approximation of u ( x , t ) CH ω ^ ( x , t ) m ( I ) . Then, the following estimation is satisfied
u ( x , t ) u N ( x , t ) L ω ^ ( x , t ) 2 γ 3 N 7 m 4 m u ( x , t ) t m L ω ^ ( x , t ) 2 ,
where γ 3 is a positive constant independent of any function.
Proof. 
The proof of this corollary is a direct result of Theorem 8. ◻
Theorem 9.
Let u N ( x , t ) be the shifted fifth-kind Chebyshev approximation of u ( x , t ) CH ω ^ ( x , t ) m ( I ) . Then, R N ( x , t ) L ω ^ ( x , t ) 2 0 as N .
Proof. 
Based on Theorem 7, R N ( x , t ) L ω ^ ( x , t ) 2 can be written as
R N ( x , t ) L ω ^ ( x , t ) 2 ( u ( x , t ) u N ( x , t ) ) t L ω ^ ( x , t ) 2 + | k | 2 ( u ( x , t ) u N ( x , t ) ) x 2 L ω ^ ( x , t ) 2 .
Now, the application of Theorem 8 enables us to write Equation (33) as
R N ( x , t ) L ω ^ ( x , t ) 2 3 γ 2 N 7 4 ( m 1 ) m u ( x , t ) t m L ω ^ ( x , t ) 2 + 3 | k | γ 1 N 7 4 ( m 2 ) m u ( x , t ) t m L ω ^ ( x , t ) 2 ,
and hence, it is clear that R N ( x , t ) L ω ^ ( x , t ) 2 0 as N . This completes the proof of Theorem 9. ◻

5. Illustrative Examples

In order to show the convenience and validity of the presented algorithm, three numerical examples are presented accompanied by comparisons with some other methods in the literature.
Example 1.
Consider the following heat conduction equation [46]
u ( x , t ) t 2 u ( x , t ) x 2 = 0 , 0 x , 0 t τ ,
along with the following initial and boundary conditions:
u ( x , 0 ) = sin ( x ) , u ( 0 , t ) = 0 , u ( , t ) = e t sin ( ) ,
where the exact solution is: u ( x , t ) = e t sin ( x ) .
Figure 1 shows the approximate solution and the maximum absolute error ( M A E ) graphs for the case N = 12 and = τ = 1 . This figure shows that the numerical solution is close to the exact solution. In Table 1, we illustrate the absolute error (AE) for different values of t at N = 12 and = τ = 1 . This table fully shows that the expressed method has a good precision. Furthermore, the AEs for different values of t at N = 8 when = 3 and τ = 2 are shown in Table 2. We can see from Table 1 and Table 2 and Figure 1 that the proposed method is appropriate and effective.
Example 2.
Consider the following heat conduction equation [47,48]:
u ( x , t ) t 2 u ( x , t ) x 2 = 0 , 0 x , 0 t τ ,
along with the following initial and boundary conditions:
u ( x , 0 ) = e x 2 cos x 2 , u ( 0 , t ) = cos ( t ) , u ( , t ) = e 2 cos t 2 ,
where the exact solution is: u ( x , t ) = e x 2 cos ( t x 2 ) .
In Figure 2, we sketched the exact and approximate solutions for the case N = 12 and = τ = 1 . This figure shows that the numerical and exact solutions are almost identical. In Figure 3, we plotted the M A E s when = 3 and for different values of τ at N = 12 . In Table 3, we list the M A E s for different values of N and for = τ = 1 . In Table 4, we give a comparison between the M A E s obtained from the application of the numerical scheme presented in [47] and our method. The results of Table 3 and Table 4 and Figure 2 and Figure 3 show that our numerical results when taking few terms of the proposed shifted fifth-kind Chebyshev expansion are more accurate. This demonstrates the advantage of our method when compared with some other numerical methods.
Example 3.
Consider the following heat conduction equation
u ( x , t ) t 2 u ( x , t ) x 2 = ( t + 1 ) 2 e x t 2 + 1 2 , 0 x , 0 t τ ,
along with the following initial and boundary conditions:
u ( x , 0 ) = e x , u ( 0 , t ) = 1 t 2 + 1 , u ( , t ) = e t 2 + 1 ,
where the exact solution is: u ( x , t ) = e x t 2 + 1 .
In Figure 4, we illustrate the approximate solution and the M A E graphs for the case N = 12 and = τ = 1 . Table 5 shows the AEs for different values of t at N = 12 and = τ = 1 . This table fully reveals that the expressed method has a good precision. In Table 6, we report the AEs for different values of t at N = 10 when = 3 and τ = 1 . We can see from the tabulated AEs of Table 5 and Table 6 and Figure 4 that the proposed method is suitable and powerful for solving the heat conduction equation.

6. Concluding Remarks

In this paper, we treated numerically one of the well-known equations named the heat conduction equation. The shifted Chebyshev polynomials of the fifth kind were used as basis functions. Some new theoretical results concerning specific formulas of the derivatives and integrals formulas of these polynomials were established. A numerical scheme to solve this equation was analyzed and implemented in detail. The basic idea behind the proposed algorithm was built on solving the corresponding integral equation to the heat conduction equation, and after that employing the spectral tau method to convert the integral equation governed by its boundary condition into an algebraic system of equations that could be solved via a suitable numerical solver. The performance of our presented method was evaluated in terms of absolute errors and maximum absolute errors. The numerical results demonstrated the good accuracy of this scheme and the ability to simulate the exact solution well. All codes were written and debugged by Mathematica 11 on HP Z420 Workstation, Processor: Intel (R) Xeon(R) CPU E5-1620—3.6 GHz, 16 GB Ram DDR3, and 512 GB storage.

Author Contributions

Conceptualization, W.M.A.-E., Y.H.Y.; Methodology, A.G.A., W.M.A.-E. and Y.H.Y.; Software, A.G.A., Y.H.Y.; Validation, W.M.A.-E.; Formal Analysis, YHY; Investigation, W.M.A.-E.; Resources, Y.H.Y.; Data Curation, A.G.A.; Writing—Original Draft Preparation, A.G.A., W.M.A.-E., Y.H.Y.; Writing—Review & Editing, W.M.A.-E.; Visualization, A.G.A., Y.H.Y.; Supervision, W.M.A.-E., G.M.M., Y.H.Y.; Project Administration, W.M.A.-E., Y.H.Y.; Funding Acquisition, Y.H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The authors received no funding for this study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fourier, J. The Analytical Theory of Heat; The University Press: Oxford, UK, 1878. [Google Scholar]
  2. Narasimhan, T.N. Fourier’s heat conduction equation: History, influence, and connections. Rev. Geophys. 1999, 37, 151–172. [Google Scholar]
  3. Meyu, G.; Koriche, K. Analytical solution vs. numerical solution of heat equation flow through rod of length 8 units in one dimension. Int. J. Appl. Math. Theor. Phys. 2021, 7, 53. [Google Scholar] [CrossRef]
  4. Liu, C.S.; Chang, J.R. Solving a nonlinear heat equation with nonlocal boundary conditions by a method of nonlocal boundary shape functions. Numer. Heat Transf. Part B Fundam. 2021, 80, 1–13. [Google Scholar]
  5. Tassaddiq, A.; Yaseen, M.; Yousaf, A.; Srivastava, R. A cubic B-spline collocation method with new approximation for the numerical treatment of the heat equation with classical and non-classical boundary conditions. Phys. Scr. 2021, 96, 045212. [Google Scholar] [CrossRef]
  6. Dehghan, M.; Shirilord, A. The use of homotopy analysis method for solving generalized Sylvester matrix equation with applications. Eng. Comput. 2021, 38, 2699–2716. [Google Scholar] [CrossRef]
  7. Saratha, S.R.; Krishnan, G.S.S.; Bagyalakshmi, M. Analysis of a fractional epidemic model by fractional generalised homotopy analysis method using modified Riemann-Liouville derivative. Appl. Math. Model. 2021, 92, 525–545. [Google Scholar] [CrossRef]
  8. Gonzalez-Gaxiola, O.; Biswas, A.; Ekici, M.; Khan, S. Highly dispersive optical solitons with quadratic–cubic law of refractive index by the variational iteration method. J. Opt. 2021, 51, 29–36. [Google Scholar] [CrossRef]
  9. Doeva, O.; Masjedi, P.K.; Weaver, P.M. A semi-analytical approach based on the variational iteration method for static analysis of composite beams. Compos. Struct. 2021, 257, 113110. [Google Scholar]
  10. Noeiaghdam, S.; Sidorov, D.; Wazwaz, A.M.; Sidorov, N.; Sizikov, V. The Numerical Validation of the Adomian Decomposition Method for Solving Volterra Integral Equation with Discontinuous Kernels Using the CESTAC Method. Mathematics 2021, 9, 260. [Google Scholar] [CrossRef]
  11. Lu, T.T.; Zheng, W.Q. Adomian decomposition method for first order PDEs with unprescribed data. Alex. Eng. J. 2021, 60, 2563–2572. [Google Scholar]
  12. Li, P.W. Space–time generalized finite difference nonlinear model for solving unsteady Burgers’ equations. Appl. Math. Lett. 2021, 114, 106896. [Google Scholar]
  13. Bullo, T.; Duressa, G.F.; Degla, G. Accelerated fitted operator finite difference method for singularly perturbed parabolic reaction-diffusion problems. Comput. Meth. Diff. Eqs. 2021, 9, 886–898. [Google Scholar]
  14. Vargas, A.M. Finite difference method for solving fractional differential equations at irregular meshes. Math. Comput. Simul. 2022, 193, 204–216. [Google Scholar] [CrossRef]
  15. Gharibi, Z.; Dehghan, M. Convergence analysis of weak Galerkin flux-based mixed finite element method for solving singularly perturbed convection-diffusion-reaction problem. Appl. Numer. Math. 2021, 163, 303–316. [Google Scholar]
  16. Yu, B.; Hu, P.; Saputra, A.A.; Gu, Y. The scaled boundary finite element method based on the hybrid quadtree mesh for solving transient heat conduction problems. Appl. Math. Model. 2021, 89, 541–571. [Google Scholar] [CrossRef]
  17. Namala, S.; Uddin, R. Hybrid nodal integral/finite element method for time-dependent convection diffusion equation. J. Nuc. Eng. Radiat. Sci. 2022, 8, 021406. [Google Scholar] [CrossRef]
  18. Zhang, H.; Liu, F.; Jiang, X.; Turner, I. Spectral method for the two-dimensional time distributed-order diffusion-wave equation on a semi-infinite domain. J. Comput. Appl. Math. 2022, 399, 113712. [Google Scholar]
  19. Atta, A.G.; Moatimid, G.M.; Youssri, Y.H. Generalized Fibonacci Operational tau Algorithm for Fractional Bagley-Torvik Equation. Prog. Fract. Differ. Appl. 2020, 6, 215–224. [Google Scholar]
  20. Tu, H.; Wang, Y.; Lan, Q.; Liu, W.; Xiao, W.; Ma, S. A Chebyshev-Tau spectral method for normal modes of underwater sound propagation with a layered marine environment. J. Sound Vib. 2021, 492, 115784. [Google Scholar] [CrossRef]
  21. Doha, E.H.; Abd-Elhameed, W.M. Accurate spectral solutions for the parabolic and elliptic partial differential equations by the ultraspherical tau method. J. Comput. Appl. Math. 2005, 181, 24–45. [Google Scholar] [CrossRef] [Green Version]
  22. Atta, A.G.; Moatimid, G.M.; Youssri, Y.H. Generalized Fibonacci operational collocation approach for fractional initial value problems. Int. J. Appl. Comput. Math. 2019, 5, 1–11. [Google Scholar] [CrossRef]
  23. Hesthaven, J.S.; Gottlieb, S.; Gottlieb, D. Spectral Methods for Time-Dependent Problems; Cambridge University Press: Cambridge, UK, 2007; Volume 21. [Google Scholar]
  24. Shizgal, B. Spectral Methods in Chemistry and Physics; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  25. Costabile, F.; Napoli, A. A class of Birkhoff–Lagrange-collocation methods for high order boundary value problems. Appl. Numer. Math. 2017, 116, 129–140. [Google Scholar]
  26. Youssri, Y.H. Orthonormal Ultraspherical Operational Matrix Algorithm for Fractal–Fractional Riccati Equation with Generalized Caputo Derivative. Fractal Fract. 2021, 5, 100. [Google Scholar]
  27. Abd-Elhameed, W.M.; Youssri, Y.H. Explicit shifted second-kind Chebyshev spectral treatment for fractional Riccati differential equation. CMES Comput. Model. Eng. Sci. 2019, 121, 1029–1049. [Google Scholar]
  28. Atta, A.G.; Abd-Elhameed, W.M.; Moatimid, G.M.; Youssri, Y.H. Shifted fifth-kind Chebyshev Galerkin treatment for linear hyperbolic first-order partial differential equations. Appl. Numer. Math. 2021, 167, 237–256. [Google Scholar]
  29. Atta, A.G.; Abd-Elhameed, W.M.; Youssri, Y.H. Shifted fifth-kind Chebyshev polynomials Galerkin-based procedure for treating fractional diffusion-wave equation. Int. J. Mod. Phys. C 2022, 33, 2250102. [Google Scholar] [CrossRef]
  30. Khalil, H.; Khan, R.A. A new method based on Legendre polynomials for solutions of the fractional two-dimensional heat conduction equation. Comput. Math. Appl. 2014, 67, 1938–1953. [Google Scholar]
  31. Abbasbandy, S.; Kazem, S.; Alhuthali, M.S.; Alsulami, H.H. Application of the operational matrix of fractional-order Legendre functions for solving the time-fractional convection–diffusion equation. Appl. Math. Comput. 2015, 266, 31–40. [Google Scholar] [CrossRef]
  32. Abdelkawy, M.A.; Alyami, S.A. Legendre-Chebyshev spectral collocation method for two-dimensional nonlinear reaction-diffusion equation with Riesz space-fractional. Chaos Solitons Fractals 2021, 151, 111279. [Google Scholar]
  33. Atta, A.G.; Abd-Elhameed, W.M.; Moatimid, G.M.; Youssri, Y.H. Advanced shifted sixth-kind Chebyshev tau approach for solving linear one-dimensional hyperbolic telegraph type problem. Math. Sci. 2022. [Google Scholar] [CrossRef]
  34. Pourbabaee, M.; Saadatmandi, A. Collocation method based on Chebyshev polynomials for solving distributed order fractional differential equations. Comput. Meth. Differ. Eqs. 2021, 9, 858–873. [Google Scholar]
  35. Mohammadi, F.; Moradi, L.; Conte, D. Discrete Chebyshev polynomials for solving fractional variational problems. Stat. Optim. Inf. Comput. 2021, 9, 502–515. [Google Scholar] [CrossRef]
  36. Abbas, W.; Fathy, M.; Mostafa, M.; Hashem, A.M.A. Galerkin method for nonlinear Volterra-Fredholm integro-differential equations based on Chebyshev polynomials. Eng. Res. J. 2021, 170, 169–183. [Google Scholar]
  37. Panghal, S.; Kumar, M. Multilayer perceptron and Chebyshev polynomials-based functional link artificial neural network for solving differential equations. Int. J. Model. Simulat. Sci. Comput. 2021, 12, 2150011. [Google Scholar] [CrossRef]
  38. Abd-Elhameed, W.M. Novel expressions for the derivatives of sixth kind Chebyshev polynomials: Spectral solution of the non-linear one-dimensional Burgers’ Equation. Fractal Fract. 2021, 5, 53. [Google Scholar] [CrossRef]
  39. Doha, E.H.; Abd-Elhameed, W.M.; Bassuony, M.A. On the coefficients of differentiated expansions and derivatives of Chebyshev polynomials of the third and fourth kinds. Acta Math. Sci. 2015, 35, 326–338. [Google Scholar] [CrossRef]
  40. Doha, E.H.; Abd-Elhameed, W.M.; Bassuony, M.A. On using third and fourth kinds Chebyshev operational matrices for solving Lane-Emden type equations. Rom. J. Phys. 2015, 60, 281–292. [Google Scholar]
  41. Sadri, K.; Aminikhah, H. A new efficient algorithm based on fifth-kind Chebyshev polynomials for solving multi-term variable-order time-fractional diffusion-wave equation. Int. J. Comput. Math. 2021, 99, 966–992. [Google Scholar] [CrossRef]
  42. Abd-Elhameed, W.M.; Youssri, Y.H. New formulas of the high-order derivatives of fifth-kind Chebyshev polynomials: Spectral solution of the convection–diffusion equation. Numer. Methods Partial. Differ. Equ. 2021. [Google Scholar] [CrossRef]
  43. Abd-Elhameed, W.M.; Youssri, Y.H. Fifth-kind orthonormal Chebyshev polynomial solutions for fractional differential equations. Comput. Appl. Math. 2018, 37, 2897–2921. [Google Scholar] [CrossRef]
  44. Koepf, W. Hypergeometric Summation, 2nd ed.; Springer Universitext Series: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  45. Li, J.; Cheng, Y. Linear barycentric rational collocation method for solving heat conduction equation. Numer. Methods Partial. Differ. Equ. 2021, 37, 533–545. [Google Scholar] [CrossRef]
  46. Merga, F.E.; Chemeda, H.M. Modified Crank-Nicolson scheme with Richardson extrapolation for one-dimensional heat equation. Iran. J. Sci. Technol. Trans. Sci. 2021, 45, 1725–1734. [Google Scholar] [CrossRef]
  47. Hamaidi, M.; Naji, A.; Charafi, A. Space–time localized radial basis function collocation method for solving parabolic and hyperbolic equations. Eng. Anal. Bound. Elem. 2016, 67, 152–163. [Google Scholar] [CrossRef]
  48. Ku, C.Y.; Liu, C.Y.; Xiao, J.E.; Chen, M.R. Solving backward heat conduction problems using a novel space–time radial polynomial basis function collocation method. Appl. Sci. 2020, 10, 3215. [Google Scholar]
Figure 1. The approximate solution and the M A E graphs of Example 1.
Figure 1. The approximate solution and the M A E graphs of Example 1.
Fractalfract 06 00619 g001
Figure 2. The exact and approximate solutions of Example 2.
Figure 2. The exact and approximate solutions of Example 2.
Fractalfract 06 00619 g002
Figure 3. M A E graphs of Example 2.
Figure 3. M A E graphs of Example 2.
Fractalfract 06 00619 g003
Figure 4. The approximate solution and the M A E graphs of Example 3.
Figure 4. The approximate solution and the M A E graphs of Example 3.
Fractalfract 06 00619 g004
Table 1. The AEs of Example 1.
Table 1. The AEs of Example 1.
x t = 2 10 t = 4 10 t = 6 10 t = 8 10
0.1 1.07216 × 10 6 1.04375 × 10 7 1.04421 × 10 8 1.74522 × 10 9
0.2 1.01956 × 10 6 1.66463 × 10 7 1.97095 × 10 8 2.84439 × 10 9
0.3 9.15588 × 10 7 2.09069 × 10 7 2.70454 × 10 8 3.78237 × 10 9
0.4 7.82278 × 10 7 2.30326 × 10 7 3.17442 × 10 8 4.40703 × 10 9
0.5 6.37195 × 10 7 2.30262 × 10 7 3.33491 × 10 8 4.62481 × 10 9
0.6 4.92708 × 10 7 2.10443 × 10 7 3.17013 × 10 8 4.39902 × 10 9
0.7 3.55998 × 10 7 1.73636 × 10 7 2.69589 × 10 8 3.74448 × 10 9
0.8 2.29584 × 10 7 1.23449 × 10 7 1.95837 × 10 8 2.72157 × 10 9
0.9 1.12214 × 10 7 6.4063 × 10 7 1.02952 × 10 8 1.42996 × 10 9
Table 2. The AEs of Example 1.
Table 2. The AEs of Example 1.
x t = 3 10 t = 7 10 t = 11 10 t = 15 10 t = 19 10
0.29.81704 × 10 5 2.42225 × 10 4 2.24868 × 10 6 5.37087 × 10 3 3.24886 × 10 3
0.41.61729 × 10 4 5.11161 × 10 4 4.67631 × 10 4 2.53555 × 10 3 1.11043 × 10 3
0.61.84251 × 10 4 8.30112 × 10 4 8.43864 × 10 4 8.92868 × 10 4 1.84493 × 10 3
0.81.26461 × 10 4 1.21323 × 10 3 1.14871 × 10 3 1.94866 × 10 4 4.31002 × 10 3
16.98793 × 10 5 1.66147 × 10 3 1.38598 × 10 3 2.11464 × 10 4 1.64701 × 10 3
1.24.82847 × 10 4 2.15759 × 10 3 1.54812 × 10 3 7.01591 × 10 4 3.46544 × 10 2
1.41.19936 × 10 3 2.66197 × 10 3 1.62211 × 10 3 1.39606 × 10 3 5.36674 × 10 2
1.62.29473 × 10 3 3.11052 × 10 3 1.59488 × 10 3 1.99516 × 10 3 6.71459 × 10 2
1.83.79811 × 10 3 3.41681 × 10 3 1.45906 × 10 3 2.18735 × 10 3 6.83315 × 10 2
25.64132 × 10 3 3.48161 × 10 3 1.21885 × 10 3 1.69924 × 10 3 5.15014 × 10 3
2.27.58846 × 10 3 3.21358 × 10 3 8.95812 × 10 4 3.90051 × 10 4 1.43644 × 10 2
2.49.14396 × 10 3 2.56621 × 10 3 5.33517 × 10 4 1.59222 × 10 3 3.82048 × 10 2
2.69.43725 × 10 3 1.59666 × 10 3 2.00103 × 10 4 3.57282 × 10 3 8.86543 × 10 2
2.87.08014 × 10 3 5.53469 × 10 4 1.32394 × 10 5 4.00638 × 10 3 9.89625 × 10 2
Table 3. M A E s of Example 2.
Table 3. M A E s of Example 2.
N24681012
M A E 3.73882 × 10 1 1.57463 × 10 2 1.50513 × 10 2 3.26227 × 10 3 1.16386 × 10 3 1.30992 × 10 8
Table 4. Comparison between our method and the method in [47] for Example 2.
Table 4. Comparison between our method and the method in [47] for Example 2.
Our method 1.30992 × 10 8
method in [47] 2.28 × 10 5
Table 5. The AEs of Example 3.
Table 5. The AEs of Example 3.
x t = 1 10 t = 3 10 t = 5 10 t = 7 10 t = 9 10
0.1 2.53661 × 10 6 6.78466 × 10 7 4.66726 × 10 8 6.34649 × 10 9 1.09628 × 10 7
0.2 2.12895 × 10 6 8.21776 × 10 7 9.03818 × 10 8 1.21842 × 10 8 2.68961 × 10 9
0.3 1.75779 × 10 6 8.88757 × 10 7 1.24666 × 10 7 1.69281 × 10 8 1.93979 × 10 8
0.4 1.43444 × 10 6 8.86624 × 10 7 1.46208 × 10 7 2.00541 × 10 8 5.19067 × 10 9
0.5 1.15321 × 10 6 8.25147 × 10 7 1.53111 × 10 7 2.11703 × 10 8 1.08241 × 10 8
0.6 9.01903 × 10 7 7.15177 × 10 7 1.44973 × 10 7 2.01492 × 10 8 2.08671 × 10 8
0.7 6.69298 × 10 7 5.67622 × 10 7 1.22826 × 10 7 1.71658 × 10 8 3.45005 × 10 8
0.8 4.45284 × 10 7 3.92877 × 10 7 8.89726 × 10 8 1.25433 × 10 8 4.42532 × 10 9
0.9 2.25786 × 10 7 3.92877 × 10 7 4.67338 × 10 8 6.77284 × 10 9 1.53647 × 10 7
Table 6. The AEs of Example 3.
Table 6. The AEs of Example 3.
x t = 1 10 t = 4 10 t = 7 10 t = 1
0 1.75063 × 10 7 4.24278 × 10 7 2.88407 × 10 8 7.33871 × 10 5
0.3 2.16085 × 10 5 1.95217 × 10 4 2.59148 × 10 4 2.84798 × 10 4
0.6 7.64342 × 10 5 4.37391 × 10 4 5.12969 × 10 4 5.05394 × 10 4
0.9 2.26323 × 10 4 7.54558 × 10 4 7.48026 × 10 4 7.13391 × 10 4
1.2 5.74283 × 10 4 1.13808 × 10 3 9.39723 × 10 4 8.97356 × 10 4
1.5 1.24217 × 10 3 1.51991 × 10 3 1.05868 × 10 4 1.05103 × 10 3
1.8 2.28929 × 10 3 1.76226 × 10 3 1.08712 × 10 3 1.16802 × 10 3
2.1 3.55861 × 10 3 1.68673 × 10 3 1.03616 × 10 3 1.23095 × 10 3
2.4 4.47518 × 10 3 1.17726 × 10 3 9.38583 × 10 4 1.19524 × 10 3
2.7 3.87464 × 10 3 3.82554 × 10 3 7.46487 × 10 4 9.78603 × 10 4
3 5.70265 × 10 7 8.39831 × 10 8 2.23171 × 10 7 4.82824 × 10 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Atta, A.G.; Abd-Elhameed, W.M.; Moatimid, G.M.; Youssri, Y.H. Modal Shifted Fifth-Kind Chebyshev Tau Integral Approach for Solving Heat Conduction Equation. Fractal Fract. 2022, 6, 619. https://doi.org/10.3390/fractalfract6110619

AMA Style

Atta AG, Abd-Elhameed WM, Moatimid GM, Youssri YH. Modal Shifted Fifth-Kind Chebyshev Tau Integral Approach for Solving Heat Conduction Equation. Fractal and Fractional. 2022; 6(11):619. https://doi.org/10.3390/fractalfract6110619

Chicago/Turabian Style

Atta, Ahmed Gamal, Waleed Mohamed Abd-Elhameed, Galal Mahrous Moatimid, and Youssri Hassan Youssri. 2022. "Modal Shifted Fifth-Kind Chebyshev Tau Integral Approach for Solving Heat Conduction Equation" Fractal and Fractional 6, no. 11: 619. https://doi.org/10.3390/fractalfract6110619

APA Style

Atta, A. G., Abd-Elhameed, W. M., Moatimid, G. M., & Youssri, Y. H. (2022). Modal Shifted Fifth-Kind Chebyshev Tau Integral Approach for Solving Heat Conduction Equation. Fractal and Fractional, 6(11), 619. https://doi.org/10.3390/fractalfract6110619

Article Metrics

Back to TopTop