Next Article in Journal
Optimizing the Neural Structure and Hyperparameters of Liquid State Machines Based on Evolutionary Membrane Algorithm
Next Article in Special Issue
Non-Relativistic Treatment of the 2D Electron System Interacting via Varshni–Shukla Potential Using the Asymptotic Iteration Method
Previous Article in Journal
A Study of the Complexity of Problems Posed by Talented Students in Mathematics
Previous Article in Special Issue
The Non-Linear Fokker–Planck Equation in Low-Regularity Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials

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, Cairo 11566, Egypt
2
Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(11), 1843; https://doi.org/10.3390/math10111843
Submission received: 18 April 2022 / Revised: 17 May 2022 / Accepted: 24 May 2022 / Published: 27 May 2022

Abstract

:
Herein, a spectral Galerkin method for solving the fractional Rayleigh–Stokes problem involving a nonlinear source term is analyzed. Two kinds of basis functions that are related to the shifted sixth-kind Chebyshev polynomials are selected and utilized in the numerical treatment of the problem. Some specific integer and fractional derivative formulas are used to introduce our proposed numerical algorithm. Moreover, the stability and convergence accuracy are derived in detail. As a final validation of our theoretical results, we present a few numerical examples.

1. Introduction

The importance of non-Newtonian fluids in science and industrial applications has piqued the interest of numerous researchers. There are many examples of non-Newtonian fluids such as in natural substances (lava, magma, gums, honey), in biology (semen, mucus, synovia, blood), in industry (molten polymer, lubricant, paint, ink, glue), in food products (ketchup, butter, mustard, chocolate, mayonnaise, cheese), and in cosmetics (cream, silicone, toothpaste, nail polish, soap solution). In this regard, the fractional Rayleigh–Stokes equation (FRSE) plays an important role in describing the dynamic behavior of some non-Newtonian fluids [1,2,3,4].
The nonlinear FRSE [5] is as follows:
u t ( x , t ) D t α [ a u x x ( x , t ) ] b u x x ( x , t ) = f ( u ( x , t ) , x , t ) , 0 < α < 1 ,
with, respectively, the following homogeneous initial and boundary conditions:
u ( x , 0 ) = 0 , 0 < x < ,
and
u ( 0 , t ) = u ( , t ) = 0 , 0 < t τ ,
where a and b are two positive constants and the nonlinear source term f ( u ( x , t ) , x , t ) satisfies the global Lipschitz condition with respect to u ( x , t ) . The symbol D t α is the Caputo fractional derivative operator of order α that describes the viscoelastic behavior of the flow. Some researchers have investigated and proposed a few methods for the solution of FRSE. In Ref. [6], the authors proposed the finite element method for the numerical solution of FRSE. In Ref. [7], the authors applied the radial basis function-generated finite difference method for the solution of the FRSE, while in [8], the authors solved the FRSE by using the spectral Jacobi–Galerkin method. Furthermore, an improved tau method for the multi-dimensional FRSE for a heated generalized second grade fluid was developed in [9]. Some other studies regarding the Rayleigh–Stokes problem can be found in [10,11].
Chebyshev polynomials (CPs) play significant roles in numerical analysis and approximation theory. There are well-known four kinds of CPs, which are specific types of Jacobi polynomials. These kinds of polynomials have been extensively used in a variety of papers related to numerical analysis; see, for instance, ref. [12,13,14]. The other two kinds of Chebyshev polynomials, namely, the fifth and sixth kinds of Chebyshev polynomials were investigated in [15]. These two classes are symmetric like the first and second kinds of Chebyshev polynomials. In fact, they are particular polynomials of the so-called “generalized ultraspherical polynomials” (see, for example, ref. [16,17]). Regrading these polynomials, their theoretical, as well as practical aspects, have attracted a significant attention from several authors. In this respect, Xu et al. [18] treated the fractional optimal control problems using sixth-kind Chebyshev wavelets. Moreover, Babaei et al. [19] employed Chebyshev polynomials of the sixth kind for solving the variable-order fractional nonlinear quadratic integro-differential equations. In addition, Jafari et al. [20] developed a spectral collocation method for treating the inverse reaction-diffusion–convection based on Chebyshev polynomials of the sixth kind. Some other contributions regarding sixth-kind Chebyshev polynomials can be found in [21,22], while for some contributions regarding Chebyshev polynomials of the fifth kind, one can referred to [23,24].
Since obtaining an exact solution is very computationally expensive for fractional differential equations, it is therefore impossible or extremely difficult to analytically solve such models. As a consequence, it has become an active research pursuit to analyze and implement high-efficient numerical techniques such as spectral methods for the simulation of solutions to these models. Spectral methods are based on the idea that approximate solutions to differential equations can be expressed as a series of truncated special functions. Three main spectral methods are employed, namely, the collocation, tau, and Galerkin methods. Readers interested in this subject can consult [25,26,27] for detailed explanations and applications of these techniques.
The following is a brief summary of the principal aims of this article:
  • Construct and develop a new method for solving the nonlinear FRSE through shifted CPs of the sixth-kind by the application of the Galerkin method;
  • Discuss the convergence and error analysis of the presented method;
  • Present some numerical results to examine the applicability and accuracy of the algorithm.
The structure of the paper is as follows. Section 2 displays a few fundamental concepts related to Caputo fractional calculus. A few definitions and formulas concerning sixth-kind shifted CPs are also displayed. Section 3 discusses the Galerkin approach for the numerical treatment of the FRSE. The proposed double Chebyshev expansion is examined for convergence and error analysis in Section 4. Section 5 contains some numerical examples and comparisons between our numerical results and those produced by other approaches. A few conclusions are summarized in Section 6.

2. Preliminaries and Essential Relations

Essential definitions and formulas are included in this section.

2.1. Some Definitions and Properties of the Fractional Calculus

Definition 1
([28]). On the typical Lebesgue space L 1 [ 0 , 1 ] , the Riemann–Liouville fractional integral operator of order ρ is defined as
I ρ h ( y ) = 1 Γ ( ρ ) 0 y ( y t ) ρ 1 h ( t ) d t , ρ > 0 , h ( y ) , ρ = 0 .
Definition 2
([28]). The Caputo definition of the fractional-order derivative is:
D y ρ h ( y ) = 1 Γ ( m ρ ) 0 y ( y t ) m ρ 1 h ( m ) ( t ) d t , ρ > 0 , y > 0 ,
where m 1 ρ < m , m N .
The operator D y ρ fulfills the accompanying properties for m 1 ρ < m , m N ,
( i ) ( D y ρ I ρ h ) ( y ) = h ( y ) , ( i i ) ( I ρ D y ρ h ) ( y ) = h ( y ) k = 0 m 1 h ( k ) ( 0 + ) Γ ( k + 1 ) ( y a ) k , y > 0 , ( i i i ) D y ρ y k = Γ ( k + 1 ) Γ ( k + 1 ρ ) y k ρ , k N , k ρ ,
where ρ is the smallest integer greater than or equal to ρ .

2.2. Some Basic Formulas and Properties of Sixth-Kind CPs and Their Shifted Ones

Sixth-kind Chebyshev polynomials Y i ( t ) [15] are orthogonal polynomials with respect to the weight function w ˜ ( t ) = t 2 1 t 2 . The orthogonality relation of these polynomials is given by [21]
1 1 t 2 1 t 2 Y i ( t ) Y j ( t ) d t = h i , if   i = j , 0 , if   i j ,
where
h i = π 2 2 i + 3 1 , if   i even , i + 3 i + 1 , if   i odd .
These polynomials may be constructed using the following recursive formula:
Y i ( t ) = t Y i 1 ( t ) i ( i + 1 ) + ( 1 ) i ( 2 i + 1 ) + 1 4 i ( i + 1 ) Y i 2 ( t ) , Y 0 ( t ) = 1 , Y 1 ( t ) = t , i 2 .
In [21], the authors also provide trigonometric representations of sixth-kind Chebyshev polynomials as follows:
Y 2 i ( cos θ ) = sin ( 2 i + 2 ) θ 4 i sin ( 2 θ ) , Y 2 i + 1 ( cos θ ) = 2 ( i + 2 ) cos θ sin ( ( 2 i + 3 ) θ ) sin ( ( 2 i + 4 ) θ ) 4 i + 1 ( i + 1 ) sin ( 2 θ ) cos θ .
Now, we define the shifted orthogonal polynomials Y i * ( t ) on [ 0 , τ ] as:
Y i * ( t ) = Y i 2 t τ 1 .
The following recurrence relation:
Y i * ( t ) = 2 t τ 1 Y i 1 * ( t ) α i Y i 2 * ( t ) , Y 0 * ( t ) = 1 , Y 1 * ( t ) = 2 t τ 1 , i 2 ,
generates the sequence of the shifted sixth-kind CPs Y i * ( t ) on [ 0 , τ ] , with
α i = i ( i + 1 ) + ( 1 ) i ( 2 i + 1 ) + 1 4 i ( i + 1 ) .
These polynomials are orthogonal on [ 0 , τ ] with respect to ω ( t ) = ( 2 t τ ) 2 t τ t 2 .
More precisely, we have the following orthogonality relation (see, [21]):
0 τ ω ( t ) Y i * ( t ) Y j * ( t ) d t = h τ , i δ i , j ,
where δ i , j is the well-known Kronecker delta function and
h τ , i = π τ 4 2 2 i + 5 1 , if   i even , i + 3 i + 1 , if   i odd .
The power form representation of Y j * ( t ) is given by [21]
Y j * ( t ) = r = 0 j B r , j t r ,
where
B r , j = 2 2 r j τ r ( 2 r + 1 ) ! k = r + 1 2 j 2 ( 1 ) j 2 + k + r ( 2 k + r + 1 ) ! ( 2 k r ) ! , if j even , 2 j + 1 k = r 2 j 1 2 ( 1 ) j + 1 2 + k + r ( k + 1 ) ( 2 k + r + 2 ) ! ( 2 k r + 1 ) ! , if j odd .
Another important formula of the Y j * ( t ) is its inversion formula [21]
t j = i = 0 j Q i , j Y i * ( t ) ,
where
Q i , j = τ j ( 2 j + 1 ) ! 2 i 2 j + 2 ( j i ) ! ( i + j + 4 ) ! ( i + 2 ) i ( i + 4 ) + j 2 + j + 3 , if   i even , ( i + 1 ) ( i ( i + 4 ) + j ( j + 3 ) + 6 ) , if i odd .
For more properties about Y j * ( t ) , see [21,22].
Theorem 1.
The first derivative of Y i * ( t ) is given by [29]
d Y i * ( t ) d t = r = 0 i 1 M r , i Y r * ( t ) , i 1 ,
and the coefficients M r , i are given by
M r , i = 2 2 i + r τ r + 1 , i even , r odd , i ( r + 2 ) i + 1 , i odd , i r 1 2 even , ( i + 4 ) ( r + 2 ) i + 1 , i odd , i r 3 2 even , 0 , otherwise .

3. Galerkin Approach for Treating the FRSE

We begin by selecting our basis functions in this section. Then, using the spectral Galerkin approach, we present a numerical solution for solving the FRSE with homogeneous initial and boundary conditions.

3.1. Basis Functions Selection

The following are the basis functions that we choose:
ϕ i ( t ) = t Y i * ( t ) , ψ i ( x ) = x ( x ) Y i * ( x ) .
The orthogonality relations of ϕ i ( t ) and ψ i ( x ) are respectively given by:
0 τ ( 2 t τ ) 2 t τ t 2 t 2 ϕ i ( t ) ϕ j ( t ) d t = h τ , i δ i , j ,
and
0 ( 2 x ) 2 x 3 2 ( x ) 3 2 ψ i ( x ) ψ j ( x ) d x = h , i δ i , j ,
where h τ , i is as given in (13).
Theorem 2.
The second-order derivative of ψ i ( x ) can be written as [29]:
d 2 ψ i ( x ) d x 2 = j = 0 i λ j , i Y j * ( x ) ,
where
λ j , i = 1 2 i j 1 ( i + 1 ) ( i + 2 ) 2 , i f   i = j , j + 2 , i f   i , j even a n d i j + 2 2 odd , 3 ( j + 2 ) , i f   i , j even a n d i j + 2 2 even , ( j + 1 ) ( i + 2 j + 6 ) i + 1 , i f   i , j odd a n d i j + 2 2 even , ( j + 1 ) ( i + 2 j + 2 ) i + 1 , i f   i , j odd a n d i j + 2 2 odd , 0 , otherwise .
Theorem 3.
The first-order derivative of ϕ i ( t ) is given by
d ϕ i ( t ) d t = r = 0 i A r , i Y r * ( t ) , i 1 ,
where the coefficients A r , i are given by
A r , i = τ 2 M r 1 , i + M r , i + α r + 2 M r + 1 , i , 1 r i 2 , 2 τ + M r 1 , i , r = i , M r , i + M r 1 , i , r = i 1 , M 0 , i + α 2 M 1 , i , r = 0 .
Theorem 4.
The following approximation formula holds for 0 < α < 1
D t α ϕ j ( t ) k = 0 N σ k , r , j , α Y k * ( t ) ,
where
σ k , r , j , α = r = 0 j ( r + 1 ) ! B r , j ρ k , r + 1 α ( r + 1 α ) ! .
Remark 1.
The proofs of Theorems 3 and 4 are given in the Appendix A at the end of this paper.

3.2. Galerkin Solution for the FRSE

Now, consider the following two spaces:
Λ = span { ψ i ( x ) ϕ i ( t ) : i , j = 0 , 1 , N } , Δ = { u ( x , t ) Λ : u ( 0 , t ) = u ( , t ) = u ( x , 0 ) = 0 , 0 < x < , 0 < t τ } ,
then, any function u ˜ ( x , t ) Δ may be written as
u ˜ ( x , t ) = i = 0 N j = 0 N c i j ψ i ( x ) ϕ j ( t ) .
Thanks to Theorems 2–4 along with the recurrence relation (3), we have the following expressions:
u ˜ t ( x , t ) = 2 4 i = 0 N j = 0 N r = 0 j A r , j c i j [ Y i + 2 * ( x ) + 1 α i + 1 α i + 2 Y i * ( x ) α i α i + 1 Y i 2 * ( x ) ] Y r * ( t ) , u ˜ x x ( x , t ) = τ 2 i = 0 N j = 0 N s = 0 i λ s , i c i j Y s * ( x ) Y j + 1 * ( t ) + Y j * ( t ) + α j + 1 Y j 1 * ( t ) , D t α [ u ˜ x x ( x , t ) ] τ 2 i = 0 N j = 0 N s = 0 i k = 0 N λ s , i c i j σ k , r , j , α Y s * ( x ) Y k * ( t ) .
Now, the residual of Equation (1) may be written in the following form:
R ( x , t ) = u ˜ t ( x , t ) D t α [ a u ˜ x x ( x , t ) ] b u ˜ x x ( x , t ) f ( u ˜ ( x , t ) , x , t ) = 2 4 i = 0 N j = 0 N r = 0 j A r , j c i j [ Y i + 2 * ( x ) + 1 α i + 1 α i + 2 Y i * ( x ) α i α i + 1 Y i 2 * ( x ) ] Y r * ( t ) b τ 2 i = 0 N j = 0 N s = 0 i λ s , i c i j Y s * ( x ) Y j + 1 * ( t ) + Y j * ( t ) + α j + 1 Y j 1 * ( t ) a τ 2 i = 0 N j = 0 N s = 0 i k = 0 N λ s , i c i j σ k , r , j , α Y s * ( x ) Y k * ( t ) f ( u ˜ ( x , t ) , x , t ) .
The following system of equations can be obtained using the Galerkin method as follows:
0 τ 0 R ( x , t ) ψ m ( x ) ϕ n ( t ) w ( x , t ) d x d t = 0 , 0 m , n N ,
where w ( x , t ) = ( 2 x ) 2 ( 2 t τ ) 2 t τ t 2 t x ( x ) .
By virtue of Equation (8), we can rewrite Equation (9) as
2 4 i = 0 N j = 0 N r = 0 j A r , j c i j [ h , i + 2 δ i + 2 , m + 1 α i + 1 α i + 2 h , i δ i , m α i α i + 1 h , i 2 δ i 2 , m ] δ r , n h τ , r b τ 2 i = 0 N j = 0 N s = 0 i λ s , i c i j δ s , m h , s h τ , j + 1 δ j + 1 , n + h τ , j δ j , n + α j + 1 h τ , j 1 δ j 1 , n a τ 2 i = 0 N j = 0 N s = 0 i λ s , i c i j σ n , r , j , α δ s , m h , s h τ , n f ¯ m , n = 0 , 0 m , n N ,
where
f ¯ m , n = 0 τ 0 f ( u ˜ ( x , t ) , x , t ) ψ m ( x ) ϕ n ( t ) w ( x , t ) d x d t .
Equation (10) constructs a system of non-linear algebraic equations with unknown expansion coefficients c i j of dimension ( N + 1 ) 2 , which can be solved using the well-known Newton’s iterative approach with zero initial approximations, and thus an approximation of the solution can be obtained.

3.3. Transformation to the Homogeneous Initial and Boundary Conditions

By virtue of the following transformation:
u ( x , t ) : = v ( x , t ) + v ¯ ( x , t ) ,
where
v ¯ ( x , t ) = 1 x u ( 0 , t ) u ( 0 , 0 ) + x u ( , t ) u ( , 0 ) + u ( x , 0 ) .
the FRSE (1) with non-homogeneous initial and boundary conditions can be transformed into the following form:
v t ( x , t ) D t α [ a v x x ( x , t ) ] b v x x ( x , t ) = f ^ ( v ( x , t ) , x , t ) , 0 < α < 1 ,
with homogeneous initial and boundary conditions
v ( x , 0 ) = 0 , 0 < x < , v ( 0 , t ) = v ( , t ) = 0 , 0 < t τ ,
where
f ^ ( v ( x , t ) , x , t ) = f ( u ( x , t ) , x , t ) v ¯ t ( x , t ) + D t α [ a v ¯ x x ( x , t ) ] + b v ¯ x x ( x , t ) .

4. Convergence Analysis

We present an upper estimate for the truncation error as well as the stability of the proposed approximate solution in this section.
Theorem 5.
Consider the function: u ( x , t ) = γ 1 ( x ) γ 2 ( t ) L w ( x , t ) 2 , with γ 1 ( x ) and γ 2 ( t ) having bounded third derivatives that satisfy the expansion:
u ( x , t ) = i = 0 j = 0 c i j ψ i ( x ) ϕ j ( t ) .
Then, the above series (11) is uniformly convergent to u ( x , t ) and the expansion coefficients c i j satisfy the inequality:
| c i j | 1 i 3 j 3 , i , j > 3 ,
where means that a generic constant d exists such that | c i j | d i 3 j 3 .
Proof. 
The orthogonality relations of ψ i ( x ) and ϕ j ( t ) enable one to write
c i j = 1 h , i h τ , j 0 τ 0 u ( x , t ) ψ i ( x ) ϕ j ( t ) w ( x , t ) d x d t .
By the hypotheses of the theorem, we obtain
c i j = 1 h , i h τ , j 0 ( 2 x ) 2 x x 2 γ 1 ( x ) Y i * ( x ) d x 0 τ ( 2 t τ ) 2 t τ t 2 γ 2 ( t ) Y j * ( t ) d t .
By virtue of the two substitutions:
x = 2 ( 1 + cos θ 1 ) , t = τ 2 ( 1 + cos θ 2 ) ,
the last equation turns into the form
c i j = 1 4 h i 0 π γ 1 2 ( 1 + cos θ 1 ) Y i ( cos θ 1 ) sin 2 ( 2 θ 1 ) d θ 1 × 1 4 h j 0 π γ 2 τ 2 ( 1 + cos θ 2 ) Y j ( cos θ 2 ) sin 2 ( 2 θ 2 ) d θ 2 .
Now, we have four cases:
(i)
If i , j are even
Integrating Equation (13) by parts three times as followed in Theorem 5 in [21] and making use of the trigonometric representations (2) leads to the following estimation:
| c i j | 1 i 3 j 3 , i , j > 3 .
We can also deduce the following estimations in the following cases
(ii)
If i , j are odd
| c i j | 1 i 3 j 3 , i , j > 3 .
(iii)
If i is even ; j is odd
| c i j | 1 i 3 j 3 , i , j > 3 .
(iv)
If i is odd ; j is even
| c i j | 1 i 3 j 3 , i , j > 3 .
This completes the proof of Theorem 5. □
Theorem 6.
If u ( x , t ) satisfies the hypothesis of Theorem 5 and if u N ( x , t ) = i = 0 N j = 0 N c i j ψ i ( x ) ϕ j ( t ) , then the following estimate of truncation error is fulfilled:
u ( x , t ) u N ( x , t ) L w ( x , t ) 2 N 3 2 .
Proof. 
From the definition of u ( x , t ) and u N ( x , t ) , we obtain
u ( x , t ) u N ( x , t ) L w ( x , t ) 2 i = 0 N j = N + 1 c i j ψ i ( x ) ϕ j ( t ) L w ( x , t ) 2 + i = N + 1 j = 0 c i j ψ i ( x ) ϕ j ( t ) L w ( x , t ) 2 = i = 0 N j = N + 1 | c i j | ( 1 α i + 1 α i + 2 ) h , i h τ , j + i = N + 1 j = 0 | c i j | ( 1 α i + 1 α i + 2 ) h , i h τ , j = j = N + 1 | c 0 j | b 0 + | c 1 j | b 1 + | c 2 j | b 2 + | c 3 j | b 3 h τ , j + i = 4 N j = N + 1 | c i j | ( 1 α i + 1 α i + 2 ) h , i h τ , j + i = N + 1 | c i 0 | b 4 + | c i 1 | b 5 + | c i 2 | b 6 + | c i 3 | b 7 ( 1 α i + 1 α i + 2 ) h , i + i = N + 1 j = 4 | c i j | ( 1 α i + 1 α i + 2 ) h , i h τ , j ,
where
b r = ( 1 α r + 1 α r + 2 ) h , r , if   r = 0 , 1 , 2 , 3 , h τ , r , if   r = 4 , 5 , 6 , 7 .
Now, following steps similar to those given in Theorem 5, we obtain
| c 0 j | 1 j 3 , | c 1 j | 1 j 3 , | c 2 j | 1 j 3 , | c 3 j | 1 j 3 , j > 3 ,
and
| c i 0 | 1 i 3 , | c i 1 | 1 i 3 , | c i 2 | 1 i 3 , | c i 3 | 1 i 3 , i > 3 .
Inserting Equations (12), (15) and (16) into Equation (14) and using
( 1 α i + 1 α i + 2 ) < 1 , i 0 , a n d | h , j | 1 2 2 j , j 0 ,
we obtain
u ( x , t ) u N ( x , t ) L w ( x , t ) 2 j = N + 1 b ^ 0 j 3 2 j + i = 4 N j = N + 1 b ^ 1 i 3 j 3 2 i + j + i = N + 1 b ^ 2 i 3 2 i + i = N + 1 j = 4 b ^ 3 i 3 j 3 2 i + j ,
where b ^ r , r = 0 , 1 , 2 , 3 are constants.
However, for all j > 0 , we have 1 j 3 2 j < 1 j 5 2 ; thus
u ( x , t ) u N ( x , t ) L w ( x , t ) 2 j = N + 1 b ^ 0 j 5 2 + i = 4 N j = N + 1 b ^ 1 i 5 2 j 5 2 + i = N + 1 b ^ 2 i 5 2 + i = N + 1 j = 4 b ^ 3 i 5 2 j 5 2 .
and hence, the application of the integral test [30] enables us to write Equation (18) as
u ( x , t ) u N ( x , t ) L w ( x , t ) 2 N 3 2 .
Theorem 7.
Under the assumptions of Theorem 5, we have
u N + 1 ( x , t ) u N ( x , t ) L w ( x , t ) 2 2 N .
Proof. 
We have
u N + 1 ( x , t ) u N ( x , t ) L w ( x , t ) 2 = i = 0 N c i , N + 1 ψ i ( x ) ϕ N + 1 ( t ) + j = 0 N + 1 c N + 1 , j ψ N + 1 ( x ) ϕ j ( t ) L w ( x , t ) 2 i = 0 N c i , N + 1 ψ i ( x ) ϕ N + 1 ( t ) ) L w ( x , t ) 2 + j = 0 N + 1 c N + 1 , j ψ N + 1 ( x ) ϕ j ( t ) L w ( x , t ) 2 < i = 0 N | c i , N + 1 | h , i h τ , N + 1 + j = 0 N + 1 | c N + 1 , j | h , N + 1 h τ , j = h τ , N + 1 | c 0 , N + 1 | h , 0 + | c 1 , N + 1 | h , 1 + | c 2 , N + 1 | h , 2 + | c 3 , N + 1 | h , 3 + h , N + 1 | c N + 1 , 0 | h τ , 0 + | c N + 1 , 1 | h τ , 1 + | c N + 1 , 2 | h τ , 2 + | c N + 1 , 3 | h τ , 3 + h τ , N + 1 i = 4 N | c i , N + 1 | h , i + h , N + 1 j = 4 N + 1 | c N + 1 , j | h τ , j .
Based on Equations (15)–(17) and Theorem 5 we obtain
| c 0 , N + 1 | 1 N 3 , | c 1 , N + 1 | 1 N 3 , | c 2 , N + 1 | 1 N 3 , | c 3 , N + 1 | 1 N 3 , | c N + 1 , 0 | 1 N 3 , | c N + 1 , 1 | 1 N 3 , | c N + 1 , 2 | 1 N 3 , | c N + 1 , 3 | 1 N 3 , | h , j | 1 2 2 j , | h , N + 1 | 1 2 2 N , | h τ , j | 1 2 2 j , | h τ , N + 1 | 1 2 2 N , | c i , N + 1 | 1 i 3 N 3 , | c N + 1 , i | 1 i 3 N 3 .
Now, inserting Equation (20) into Equation (19) and using the inequality: 1 j 3 2 j 1 2 j , j > 0 , along with the following approximation
i = a + 1 b f ( i ) x = a b f ( x ) d x ,
where f is decreasing function, the desired result can be obtained. □

5. Illustrative Examples

In this section, the technique presented in Section 3 is applied to solve the nonlinear FRSE. Three illustrative examples are used to demonstrate the effectiveness and applicability of the proposed technique.
Example 1.
Consider the FRSE of the form
u t ( x , t ) D t α [ u x x ( x , t ) ] u x x ( x , t ) = f ( u ( x , t ) , x , t ) , 0 < α < 1 ,
where
f ( u ( x , t ) , x , t ) = sin ( u ) + 2 t x ( x ) 5 x 2 5 x + 2 6 Γ ( 3 α ) t 1 α + 3 t x 2 ( x ) 2 sin x 3 ( x ) 3 t 2 ,
along with the following initial and boundary conditions:
u ( x , 0 ) = 0 , 0 < x < , u ( 0 , t ) = 0 , u ( , t ) = 0 , 0 < t τ .
The exact solution of this problem is u ( x , t ) = x 3 ( x ) 3 t 2 .
In Table 1, we reported the computational time (CPU time) and compared the L 2 errors of the present method with method in [5] at = τ = 1 . We see in this table that the results are accurate for small choices of N. Table 2 lists the L errors for different values of α at N = 6 when = 2 and τ = 3 and the CPU time. Figure 1 illustrates the L error for α = 0.3 (left) and α = 0.7 (right) at N = 4 when = 5 and τ = 10 . We can see from Table 1 and Table 2 and Figure 1 that the proposed method is appropriate and effective. This demonstrates the advantage of our method compared to some other numerical methods.
Example 2.
Consider the FRSE of the form
u t ( x , t ) D t α [ u x x ( x , t ) ] u x x ( x , t ) = f ( u ( x , t ) , x , t ) , 0 < α < 1 ,
where
f ( u ( x , t ) , x , t ) = u 2 + sin ( π x ) 2 π 2 Γ ( 3 α ) t 2 α + π 2 t 2 + 2 t t 4 sin 2 ( π x ) ,
along with the following initial and boundary conditions:
u ( x , 0 ) = 0 , 0 < x < , u ( 0 , t ) = 0 , u ( , t ) = 0 , 0 < t τ .
The exact solution of this problem is u ( x , t ) = t 2 sin ( π x ) .
Table 3 presents the CPU time and a comparison of the L 2 errors between our proposed method and the method in [5] at = τ = 1 . It can be found that the obtained results of the presented method are more accurate than the method in [5]. Moreover, Figure 2 sketches the L error for different values of α at N = 16 when = 2 and τ = 3 . This figure show that the numerical and exact solutions are almost identical. In Table 4. we list the absolute error for α = 0.5 at N = 18 when = 5 and τ = 10 . As can be seen, the proposed method presents better accuracy.
Example 3.
Consider the FRSE of the form
u t ( x , t ) D t α [ u x x ( x , t ) ] u x x ( x , t ) = f ( u ( x , t ) , x , t ) , 0 < α < 1 ,
where
f ( u ( x , t ) , x , t ) = u 2 + e x ( 3 α ) t 2 α Γ ( 4 α ) Γ ( 4 2 α ) t 3 2 α t 3 α e 2 x t 6 2 α ,
along with the following initial and boundary conditions:
u ( x , 0 ) = 0 , 0 < x < 1 , u ( 0 , t ) = t 3 α , u ( 1 , t ) = e t 3 α , 0 < t 1 .
The exact solution of this problem is u ( x , t ) = e x t 3 α .
In Table 5, the absolute errors for the case corresponding to N = 18 , α = 0.1 , and α = 0.9 are displayed. This table confirms that the presented method has high performance and produces accurate results. In addition, Figure 3 illustrates the L error for α = 0.5 at N = 18 . The results show good agreement between the approximate solution and the exact one.

6. Concluding Remarks

The nonlinear FRSE was treated numerically by applying the spectral Galerkin method using some polynomials related to shifted sixth-kind Chebyshev polynomials as basis functions. The proposed problem is reduced to a nonlinear system of algebraic equations that can be solved using Newton’s iterative method. The resulting approximate solutions using the suggested method are extremely close to the exact ones, indicating that our proposed algorithm can efficiently solve the problem. To demonstrate the validity and enormous potential of the algorithm, comparisons are performed between our proposed approximate solutions and those developed by other methods in the literature. In this paper, Wolfram Mathematica 11.2 was used for all calculations. In future work, we think that the theoretical results in this paper will be useful for other types of differential equations. In addition, we think that we can derive other derivative formulas for some polynomials related to Chebyshev polynomials of the sixth kind, in order to handle types of fractional differential equations that involve terms of other high-order derivatives.

Author Contributions

Formal analysis, Y.H.Y.; Investigation, A.G.A.; Methodology, W.M.A.-E. and Y.H.Y.; Software, A.G.A. and Y.H.Y.; Supervision, G.M.M. and Y.H.Y.; Validation, W.M.A.-E.; Writing—original draft, A.G.A.; Writing—review & editing, W.M.A.-E. All authors contributed to the preparation of the paper. All authors have read and agreed to the submitted version of the manuscript.

Funding

The authors received no funding for this study.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proofs of Theorem 3 and 4

Proof of Theorem 3:
Proof. 
From (7), we have
d ϕ i ( t ) d t = t d Y i * ( t ) d t + Y i * ( t ) .
By virtue of Theorem 1, one obtains
d ϕ i ( t ) d t = r = 0 i 1 M r , i [ t Y r * ( t ) ] + Y i * ( t ) .
Based on the recurrence relation (3), we can write
t Y i * ( t ) = τ 2 Y i + 1 * ( t ) + Y i * ( t ) + α i + 1 Y i 1 * ( t ) .
Inserting Equation (A2) into the relation (A1), we obtain
d ϕ i ( t ) d t = τ 2 r = 0 i 1 M r , i Y r + 1 * ( t ) + Y r * ( t ) + α r + 1 Y r 1 * ( t ) + Y i * ( t ) .
The last formula after expanding and rearranging terms leads to the following formula:
d ϕ i ( t ) d t = r = 0 i A r , i Y r * ( t ) ,
and the coefficients A r , i are given by
A r , i = τ 2 M r 1 , i + M r , i + α r + 2 M r + 1 , i , 1 r i 2 , 2 τ + M r 1 , i , r = i , M r , i + M r 1 , i , r = i 1 , M 0 , i + α 2 M 1 , i , r = 0 .
This finalizes the theorem’s proof. □
Proof of Theorem 4:
Proof. 
If the operator D t α is applied to ϕ j ( t ) , then based on Equation (6), we obtain
D t α ϕ j ( t ) = r = 0 j B r , j ( r + 1 ) ! ( r + 1 α ) ! t r + 1 α .
Now, t r + 1 α is approximated in terms of Y k * ( t ) as
t r + 1 α k = 0 N ρ k , r + 1 α Y k * ( t ) .
to find ρ k , r + 1 α . Based on the orthogonality relation of Y k * ( t ) in (12), we obtain
ρ k , r + 1 α = 1 h τ , k 0 τ t r + 1 α Y k * ( t ) ω ( t ) d t = 1 h τ , k m = 0 k B m , k 0 τ t r + m α + 1 ( 2 t τ ) 2 t τ t 2 = 1 h τ , k m = 0 k B m , k 0 τ 4 t r + m α + 7 2 + τ 2 t r + m α + 3 2 4 τ t r + m α + 5 2 τ t d t = 1 h τ , k m = 0 k B m , k τ r + m α + 5 4 β r + m α + 9 2 , 3 2 + β r + m α + 5 2 , 3 2 4 β r + m α + 7 2 , 3 2 = 1 h τ , k m = 0 k B m , k τ r + m α + 5 Γ ( 3 2 ) 4 Γ ( r + m α + 9 2 ) Γ ( r + m α + 6 ) + Γ ( r + m α + 5 2 ) Γ ( r + m α + 4 ) 4 Γ ( r + m α + 7 2 ) Γ ( r + m α + 5 ) = m = 0 k π B m , k τ r + m α + 5 α 2 α ( 2 m + 2 r + 3 ) + ( m + r ) 2 + 3 m + 3 r + 5 Γ m + r α + 5 2 2 h τ , k Γ ( m + r α + 6 ) ,
where β ( . ) and Γ ( . ) are the well known beta and gamma functions, respectively.
Now, inserting Equation (A4) into Equation (A3), we obtain the result of Theorem 4. □

References

  1. Fetecu, C.; Fetecu, C. The Rayleigh–Stokes problem for heated second grade fluids. Int. J. Non-Linear Mech. 2002, 37, 1011–1015. [Google Scholar] [CrossRef]
  2. Shen, F.; Tan, W.; Zhao, Y.; Masuoka, T. The Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivative model. Nonlinear Anal. Real World Appl. 2006, 7, 1072–1080. [Google Scholar] [CrossRef]
  3. Abelman, S.; Hayat, T.; Momoniat, E. On the Rayleigh problem for a Sisko fluid in a rotating frame. Appl. Math. Comput. 2009, 215, 2515–2520. [Google Scholar] [CrossRef]
  4. Gabriele, A.; Spyropoulos, F.; Norton, I.T. Kinetic study of fluid gel formation and viscoelastic response with kappa-carrageenan. Food Hydrocoll. 2009, 23, 2054–2061. [Google Scholar] [CrossRef]
  5. Guan, Z.; Wang, X.; Ouyang, J. An improved finite difference/finite element method for the fractional Rayleigh–Stokes problem with a nonlinear source term. J. Appl. Math. Comput. 2021, 65, 451–479. [Google Scholar] [CrossRef]
  6. Dehghan, M.; Abbaszadeh, M. A finite element method for the numerical solution of Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives. Eng. Comput. 2017, 33, 587–605. [Google Scholar] [CrossRef]
  7. Nikan, O.; Golbabai, A.; Machado, J.A.T.; Nikazad, T. Numerical solution of the fractional Rayleigh–Stokes model arising in a heated generalized second-grade fluid. Eng. Comput. 2021, 37, 1751–1764. [Google Scholar] [CrossRef]
  8. Hafez, R.M.; Zaky, M.A.; Abdelkawy, M.A. Jacobi spectral Galerkin method for distributed-order fractional Rayleigh–Stokes problem for a generalized second grade fluid. Front. Phys. 2020, 7, 240. [Google Scholar] [CrossRef]
  9. Zaky, M.A. An improved tau method for the multi-dimensional fractional Rayleigh–Stokes problem for a heated generalized second grade fluid. Comput. Math. Appl. 2018, 75, 2243–2258. [Google Scholar] [CrossRef]
  10. Xue, C.; Nie, J. Exact solutions of the Rayleigh–Stokes problem for a heated generalized second grade fluid in a porous half-space. Appl. Math. Model. 2009, 33, 524–531. [Google Scholar] [CrossRef]
  11. Chen, C.; Liu, F.; Anh, V. Numerical analysis of the Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives. Appl. Math. Comput. 2008, 204, 340–351. [Google Scholar] [CrossRef] [Green Version]
  12. Abd-Elhameed, W.M.; Doha, E.H.; Youssri, Y.H.; Bassuony, M.A. New Tchebyshev-Galerkin operational matrix method for solving linear and nonlinear hyperbolic telegraph type equations. Numer. Methods Partial Differ. Equ. 2016, 32, 1553–1571. [Google Scholar] [CrossRef]
  13. Zhou, F.; Xu, X. The third kind Chebyshev wavelets collocation method for solving the time-fractional convection diffusion equations with variable coefficients. Appl. Math. Comput. 2016, 280, 11–29. [Google Scholar] [CrossRef]
  14. 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]
  15. Masjed-Jamei, M. Some New Classes of Orthogonal Polynomials and Special Functions: A Symmetric Generalization of Sturm-Liouville Problems and Its Consequences. Ph.D. Thesis, University of Kassel, Department of Mathematics, Kassel, Germany, 2006. [Google Scholar]
  16. Xu, Y. An integral formula for generalized Gegenbauer polynomials and Jacobi polynomials. Adv. Appl. Math. 2002, 29, 328–343. [Google Scholar] [CrossRef] [Green Version]
  17. Draux, A.; Sadik, M.; Moalla, B. Markov—Bernstein inequalities for generalized Gegenbauer weight. Appl. Numer. Math. 2011, 61, 1301–1321. [Google Scholar] [CrossRef]
  18. Xu, X.; Xiong, L.; Zhou, F. Solving fractional optimal control problems with inequality constraints by a new kind of Chebyshev wavelets method. J. Comput. Sci. 2021, 54, 101412. [Google Scholar] [CrossRef]
  19. Babaei, A.; Jafari, H.; Banihashemi, S. Numerical solution of variable order fractional nonlinear quadratic integro-differential equations based on the sixth-kind Chebyshev collocation method. J. Comput. Appl. Math. 2020, 377, 112908. [Google Scholar] [CrossRef]
  20. Jafari, H.; Babaei, A.; Banihashemi, S. A novel approach for solving an inverse reaction–diffusion–convection problem. J. Optim. Theory Appl. 2019, 183, 688–704. [Google Scholar] [CrossRef]
  21. Abd-Elhameed, W.M.; Youssri, Y.H. Sixth-kind Chebyshev spectral approach for solving fractional differential equations. Int. J. Nonlinear Sci. Numer. Simul. 2019, 20, 191–203. [Google Scholar] [CrossRef]
  22. 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]
  23. 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. 2022, 99, 966–992. [Google Scholar] [CrossRef]
  24. 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]
  25. 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]
  26. Faghih, A.; Mokhtary, P. An efficient formulation of Chebyshev tau method for constant coefficients systems of multi-order FDEs. J. Sci. Comput. 2020, 82, 6. [Google Scholar] [CrossRef]
  27. 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] [CrossRef]
  28. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  29. 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]
  30. Stewart, J. Single Variable Essential Calculus: Early Transcendentals; Cengage Learning: Boston, MA, USA, 2012. [Google Scholar]
Figure 1. The L error for Example 1.
Figure 1. The L error for Example 1.
Mathematics 10 01843 g001
Figure 2. The L error for Example 2.
Figure 2. The L error for Example 2.
Mathematics 10 01843 g002
Figure 3. Th L error for Example 3.
Figure 3. Th L error for Example 3.
Mathematics 10 01843 g003
Table 1. The L 2 errors for Example 1.
Table 1. The L 2 errors for Example 1.
Method in [5]Presented Method
α h = 1 5000 , T = 1 128 T = 1 5000 , h = 1 128 N = 6 CPU Time
0.1 1.1552 × 10 6 1.4408 × 10 6 8.5446 × 10 16 35.703
0.5 1.0805 × 10 6 1.4007 × 10 6 3.0349 × 10 15 36.109
0.9 8.1511 × 10 7 1.3682 × 10 6 7.5113 × 10 16 38.594
Table 2. The L errors for Example 1.
Table 2. The L errors for Example 1.
α 0.10.50.9
L error5.12335 ×   10 13 7.70939 ×   10 13 9.18376 ×   10 13
CPU time46.02943.81344.406
Table 3. The L 2 errors for Example 2.
Table 3. The L 2 errors for Example 2.
Method in [5]Presented Method
α h = 1 5000 , T = 1 128 T = 1 5000 , h = 1 128 N = 12 CPU Time
0.1 9.1909 × 10 5 5.1027 × 10 5 2.53148 × 10 12 94.218
0.5 8.4317 × 10 5 4.4651 × 10 5 1.25044 × 10 12 95.718
0.9 6.2864 × 10 5 4.0543 × 10 5 8.47489 × 10 13 94.78
Table 4. The absolute errors for Example 2.
Table 4. The absolute errors for Example 2.
xAbsolute Error ( t = 3 )Absolute Error ( t = 6 )Absolute Error ( t = 9 )
0.50.0004593950.002025540.00477932
10.0003488950.001698370.00419428
1.50.0002759780.001473360.00378444
20.0002347180.001342310.00352795
2.50.0002213720.001297960.00323485
30.0002347090.001330060.00170189
3.50.0002759190.001397840.00634153
40.0003485640.001261730.05467222
4.50.0004576830.000072230.21660543
Table 5. The absolute errors for Example 3.
Table 5. The absolute errors for Example 3.
α = 0.1 α = 0.9
x t = 1 10 t = 5 10 t = 9 10 t = 1 10 t = 5 10 t = 9 10
0.1 7.51088 × 10 9 3.37811 × 10 10 1.94761 × 10 10 6.06613 × 10 7 3.43728 × 10 7 2.23961 × 10 7
0.2 1.42794 × 10 8 6.40815 × 10 10 3.11168 × 10 10 1.11602 × 10 6 6.33681 × 10 7 4.13436 × 10 7
0.3 1.96572 × 10 8 8.81388 × 10 10 3.47166 × 10 10 1.51426 × 10 6 8.60691 × 10 7 5.63444 × 10 7
0.4 2.31269 × 10 8 1.03444 × 10 9 1.34724 × 10 8 1.78783 × 10 6 1.01624 × 10 6 6.84415 × 10 7
0.5 2.43513 × 10 8 1.04202 × 10 9 1.32679 × 10 7 1.92288 × 10 6 1.09187 × 10 6 1.04649 × 10 6
0.6 2.32053 × 10 8 5.93031 × 10 10 9.68495 × 10 7 1.90517 × 10 6 1.07976 × 10 6 3.98943 × 10 6
0.7 1.97889 × 10 8 2.20531 × 10 9 5.43109 × 10 6 1.72122 × 10 6 9.76103 × 10 7 2.04777 × 10 5
0.8 1.44202 × 10 8 1.35044 × 10 8 1.94111 × 10 5 1.35221 × 10 6 7.49924 × 10 7 5.23308 × 10 5
0.9 7.60732 × 10 9 4.18981 × 10 8 5.25662 × 10 5 7.86145 × 10 7 8.31597 × 10 7 6.27135 × 10 5
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. A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials. Mathematics 2022, 10, 1843. https://doi.org/10.3390/math10111843

AMA Style

Atta AG, Abd-Elhameed WM, Moatimid GM, Youssri YH. A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials. Mathematics. 2022; 10(11):1843. https://doi.org/10.3390/math10111843

Chicago/Turabian Style

Atta, Ahmed Gamal, Waleed Mohamed Abd-Elhameed, Galal Mahrous Moatimid, and Youssri Hassan Youssri. 2022. "A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials" Mathematics 10, no. 11: 1843. https://doi.org/10.3390/math10111843

APA Style

Atta, A. G., Abd-Elhameed, W. M., Moatimid, G. M., & Youssri, Y. H. (2022). A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials. Mathematics, 10(11), 1843. https://doi.org/10.3390/math10111843

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