Next Article in Journal
Analytical Results for Classical Axially Symmetric States of Relativistic Hydrogenic Ions in Collinear Electric and Magnetic Fields of Arbitrary Strengths
Previous Article in Journal
Fractional-View Analysis of Fokker-Planck Equations by ZZ Transform with Mittag-Leffler Kernel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Effective Projection Method for Solving a Coupled System of Fractional-Order Bagley–Torvik Equations via Fractional Shifted Legendre Polynomials

by
Saeed Althubiti
1,† and
Abdelaziz Mennouni
2,*,†
1
Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
2
Department of Mathematics, LTM, University of Batna 2, Mostefa Ben Boulaïd, Fesdis, Batna 05078, Algeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(8), 1514; https://doi.org/10.3390/sym14081514
Submission received: 26 June 2022 / Revised: 13 July 2022 / Accepted: 21 July 2022 / Published: 25 July 2022
(This article belongs to the Section Mathematics)

Abstract

:
This work proposes a valuable and successful strategy for approximating the solutions to the Bagley–Torvik system, which plays an essential role in fractional calculus. The Caputo sense is used to derive the basic conformable fractional. The Bagley–Torvik problem is numerically solved in this study using an effective symmetric projection method. From this symmetry, there are some interesting original results. The proposed approach has two key benefits. We began by converting the connected fractional Bagley–Torvik equations into two fractional-order Bagley–Torvik equations, which we then solved using the current method. Second, two linear equation systems are solved to obtain approximate solutions.

1. Introduction

Because of its applications, fractional calculus has received a lot of attention from mathematicians and scientists in recent decades. Numerous investigators have shown the ramifications of fractional order derivative as a different physical phenomenon. Nonlinear dynamics, Maclaurin theory, theoretical physics, biomedicine, polymeric materials, dynamical mechanics, electron microscopy, transportation system, discharge and viscoelasticity, confusion and geometrisc, identification and modeling, filtration, financial system, and finance are all implementations of fractional differential problems.
The Bagley–Torvik equation plays a critical part in fractional differential equations, because it appears in the mathematical model of the movement of a big thin plate immersed in a Newtonian fluid. Torvik and Bagley introduced this fractional differential equation for the first time in 1984. The researchers stated that the fractional derivative is usually noticed in the behavior of real materials (see [1]). The authors of [2] investigated the stability conditions for linear matrix inequality in fractional order systems. In [3], a novel approach for fractional-order inputted systems is presented that is deduced from the guardian map and described as a family of monoparametric pseudo-polynomials. Reference [4] discussed a numerical approximation to a class of variable-order fractional functional boundary value problems.
The Bagley–Torvik equation has been investigated by a number of researchers using various approximate numerical methods. The goal of [5] is to propose and study a novel operational Tau approach with a high order to achieve the numerical solution to a Bagley–Torvik problem that uses Chebyshev polynomials as basis functions. Some derivatives of these equations’ solutions have a singularity at their origin, according to the author. The authors of [6] proposed a new numerical method for approximating the solution to the Bagley–Torvik problem. The fundamental goal of their strategy is to obtain a more general approximation solution in the form of Müntz–Legendre polynomials. The authors of [7] proposed an effective numerical technique founded on shifted Chebyshev polynomials for deriving numerical solutions to the Bagley–Torvik equation. The authors of [8] considered restructuring the Bagley–Torvik equation as a system of fractional differential equations of the order 1 / 2 . An approximate method for solving the Bagley–Torvik problem over the range [ 0 , 1 ] was shown by the authors of [9], via the second-order Chebyshev wavelet.
Projection approximation methods are at the heart of approximation theory and can be used for a wide range of things, including solving integral and integro-differential equations (see [10,11]).
This paper serves two main purposes. On one hand, we extend projection methods to numerically solve a system of Bagley–Torvik equations. We offer a new symmetric projection approach based on the fractional shifted Legendre polynomials. This symmetry leads to some interesting original results. On the other hand, we use a new procedure to solve the system of Bagley–Torvik equations.
There are two main advantages to the suggested method. We started by converting the problem into a system of two separable fractional-order Bagley–Torvik equations, which we then solved with the current method. Second, two systems of linear equations are solved to achieve approximate solutions.

2. Preliminaries

We will start by reviewing certain basic fractional calculus vocabulary and mathematical notations in this part.
Let Γ denote the Euler Gamma function, which is fundamental in fractional calculus theory.
Definition 1.
The left sided Riemann–Liouville fractional integral of order μ > 0 of an integrable function φ : ( 0 , ) R is described by
J 0 + μ φ ( s ) : = 1 Γ ( μ ) 0 s φ ( τ ) ( s τ ) 1 μ d τ , s > 0 .
Remark 1.
The above integral can be represented in convolution form as follows:
J 0 + μ φ ( s ) = ( ψ μ φ ) ( s ) ,
where
ψ μ ( s ) : = s μ 1 / Γ ( μ ) , s > 0 , 0 s 0 .
Definition 2.
The left sided Riemann–Liouville fractional derivative of order 0 < μ < 1 of a continuous function φ : ( 0 , ) R is defined by
D 0 + μ φ ( s ) : = 1 Γ ( 1 μ ) d d s 0 s φ ( τ ) ( s τ ) μ d τ .
Remark 2.
For 0 < μ < 1 , we have
D 0 + μ φ ( s ) = d d s J 0 + 1 μ φ ( s ) .
J 0 + μ D 0 + μ φ ( s ) φ ( 0 ) = φ ( s ) φ ( 0 ) .
Definition 3.
For an absolutely continuous function ϕ ( . ) , the Caputo fractional derivative of order 0 < α < 1 is defined by
c D 0 + α ϕ ( s ) = J 0 + 1 α d d s ϕ ( s ) = 1 Γ ( 1 α ) 0 s ( s t ) α ϕ ( t ) d t .
Remark 3.
1. 
We have
c D 0 + μ φ ( s ) = J 0 + 1 μ d d s φ ( s ) ,
2. 
For a continuous function φ, the relationship between the Caputo and Riemann–Liouville fractional derivatives is provided by
c D 0 + μ φ ( s ) = D 0 + μ ( φ ( s ) φ ( 0 ) ) , 0 < μ < 1 ,
and
c D 0 + μ φ ( s ) = D 0 + μ φ ( s ) k = 0 r 1 s k k ! φ ( k ) ( 0 ) , s > 0 , r 1 < μ < r .
Reference [12] demonstrates the fundamentals of fractional differential equations.

3. System of Coupled Fractional Bagley–Torvik Equations

Let us consider the following system of coupled fractional Bagley–Torvik equations:
A 1 D 2 u ( s ) + A 2 D ν u ( s ) + A 3 v ( s ) = f ( s ) , 0 s R , 1 < ν < 2 , A 1 D 2 v ( s ) + A 2 D ν v ( s ) + A 3 u ( s ) = g ( s ) ,
depending on the initial conditions (ICs):
u ( 0 ) = α 0 , u ( 0 ) = α 1 , v ( 0 ) = δ 0 , v ( 0 ) = δ 1 ,
or the boundary conditions (BCs):
u ( 0 ) = β 0 , u ( R ) = β 1 , v ( 0 ) = γ 0 , v ( R ) = γ 1 .
Our first goal is to convert the above system of connected fractional Bagley–Torvik equations into two fractional-order Bagley–Torvik equations. As in [13], let
ζ : = u + v , F : = f g ;
ξ : = u v , G : = f + g .
Theorem 1.
The problem (1) can be represented as follows
A 1 D 2 ζ ( s ) + A 2 D ν ζ ( s ) + A 3 ζ ( s ) = G ( s ) , A 1 D 2 ξ ( s ) + A 2 D ν ξ ( s ) A 3 ξ ( s ) = F ( s ) ,
depending on the initial conditions (ICs):
ξ ( 0 ) = ρ 0 , ξ ( 0 ) = ρ 1 , ζ ( 0 ) = η 0 , ζ ( 0 ) = η 1 ,
or the boundary conditions (BCs):
ξ ( 0 ) = λ 0 , ξ ( R ) = λ 1 , ζ ( 0 ) = σ 0 , ζ ( R ) = σ 1 .
Proof. 
It is obvious that
u = ζ + ξ 2 , f = G + F 2 ;
v = ζ ξ 2 , g = G F 2 .
Substituting this into (1) produces
A 1 D 2 ζ ( s ) + ξ ( s ) + A 2 D ν ζ ( s ) + ξ ( s ) + A 3 ζ ( s ) ξ ( s ) = G ( s ) + F ( s ) ,
A 1 D 2 ζ ( s ) ξ ( s ) + A 2 D ν ζ ( s ) ξ ( s ) + A 3 ζ ( s ) + ξ ( s ) = G ( s ) F ( s ) .
By adding the two Equations (3) and (4) together and by subtracting (4) from (3) we obtain (2). Moreover, for the initial conditions, we have
ρ 0 : = α 0 δ 0 , ρ 1 : = α 1 δ 1 , η 0 : = α 0 + δ 0 , η 1 : = α 1 + δ 1 .
For the boundary conditions, we have
λ 0 : = β 0 γ 0 , λ 1 : = β 1 γ 1 , σ 0 : = β 0 + γ 0 , σ 1 : = β 1 + γ 1 .
 □

4. The Fractional-Order Legendre Polynomials

Legendre polynomials are one of physics the most important functions in physics. They are crucial in the theory of approximation, specifically in the solution of integro-differential and integral equations.
The classic Legendre polynomials defined in [ 1 , 1 ] are famously defined by the following recurrence relation:
L 0 ( s ) = 1 , L 1 ( s ) = s , L j + 1 ( s ) = 2 j + 1 j + 1 s L j ( s ) j j + 1 L j 1 ( s ) , j = 1 , 2 , 3 ,
The following are the first Legendre polynomials:
L 0 ( s ) = 1 ; L 1 ( s ) = s ; L 2 ( s ) = 1 2 3 s 2 1 ; L 3 ( s ) = 1 2 5 s 3 3 s ; L 4 ( s ) = 1 8 35 s 4 30 s 2 + 3 ; L 5 ( s ) = 1 8 63 s 5 70 s 3 + 15 s .
Let the shifted Legendre polynomials L i ( 2 τ T 1 ) be denoted by L i S ( τ ) , τ 0 , T , from which the following recurrence relation can be achieved.
L 0 S ( τ ) = 1 , L 1 S ( τ ) = 2 τ T 1 , L j + 1 S ( τ ) = 2 j + 1 j + 1 2 τ T 1 L j S ( τ ) j j + 1 L j 1 S ( τ ) , j = 1 , 2 , 3 ,
We recall that
0 T L n S ( τ ) L m S ( τ ) d τ = T 2 n + 1 n = m , 0 n m .
Regarding the formulation of the shifted Legendre polynomials, the fractional shifted Legendre polynomials (FSLPs) can be described by shifting the variable τ = t α for α > 0 . Let L n , α S ( . ) signifies the n t h FSLP; that is,
L n , α S ( t ) = L n S ( t α ) , for   all t 0 , T .
It can be provided using the fractional shifted Legendre polynomials recurrence formula as follows:
L 0 , α S ( t ) = 1 , L 1 , α S ( t ) = 2 t α T 1 , L j + 1 , α S ( t ) = 2 j + 1 j + 1 2 t α T 1 L j , α S ( t ) j j + 1 L j 1 , α S ( t ) , j = 1 , 2 , 3 ,
Furthermore, the analytical version of the FSLP L n , α S ( . ) is as follows:
L n , α S ( t ) = j = 0 n c n , j t α j , n = 0 , 1 , 2 , 3 , ,
where
c n , p : = ( 1 ) n + p ( n + p ) ! ( n p ) ! ( p ! ) 2 T p , p = 0 , , n .
Letting ω α ( t ) = t α 1 . Define the following weight inner product
ϑ , ψ ω α : = 0 T 1 α ϑ ( t ) ψ ( t ) ω α ( t ) d t .
Via the orthogonality property (5), the sequence of FSLPs is orthogonal with regard to the weight function ω ( . ) in the interval [ 0 , T 1 α ] .
0 T 1 α L m , α S ( τ ) L n , α S ( τ ) ω α ( τ ) d τ = δ m n T α ( 2 n + 1 ) .
Denoting the corresponding normalized sequence by N n , α S . That is,
N n , α S ( t ) : = L m , α S ( t ) T α ( 2 n + 1 ) .
Hence
0 T 1 α N m , α S ( t ) N n , α S ( t ) ω α ( t ) d t = δ m n .

5. Galerkin Projection Method

Let ( π n S ) n 0 denote the sequence of bounded finite rank orthogonal projections described by
π n S f : = j = 0 n 1 f , N j , α S ω α N j , α S .
Any function φ ( . ) defined over the range [ 0 , 1 ] can be expanded in terms of normalized fractional shifted Legendre polynomials (NFSLPs) as follows:
φ ( s ) φ n ( s ) = π n S φ ( s ) ,
or equivalently
φ ( s ) φ n ( s ) = C t Ψ α ( s ) ,
where
C : = [ c 0 , c 2 , , c n 1 ] t , Ψ α ( s ) = N 0 , α S ( s ) , N 1 , α S ( s ) , , N n 1 , α S ( s ) .
Several results for solving compact operator equations using the Galerkin method have been established over the last two decades. The Galerkin method is used in this section to approximate the solution of the fractional Bagley–Torvik equations. For D 2 n : = π n S D 2 π n S , D ν n : = π n S D ν π n S , we obtain the following system of Legendre–Galerkin approximate equations
A 1 π n S D 2 ζ n ( s ) + A 2 π n S D ν ζ n ( s ) + A 3 ζ n ( s ) = π n S G ( s ) , A 1 π n S D 2 ξ n ( s ) + A 2 π n S D ν ξ n ( s ) A 3 ξ n ( s ) = π n S F ( s )
which has a unique solution ( ζ n , ξ n ) , given by
ζ n = j = 0 n a n , j N j , α S , ξ n = j = 0 n b n , j N j , α S ,
for some scalars a n , j and b n , j . Equation (6) reads as
j = 0 n a n , j A 1 π n S D 2 N j , α S + A 2 π n S D ν N j , α S + A 3 N j , α S = π n S G , j = 0 n b n , j A 1 π n S D 2 N j , α S + A 2 π n S D ν N j , α S A 3 N j , α S = π n S F ,
so that
j = 0 n a n , j [ A 1 i = 0 n D 2 N j , α S , N i , α S ω α N i , α S + A 2 i = 0 n a n , j D ν N j , α S , N i , α S ω α N i , α S + A 3 N j , α S ] = i = 0 n G , N i , α S ω α N i , α S ,
j = 0 n b n , j [ A 1 i = 0 n D 2 N j , α S , N i , α S ω α N i , α S + A 2 i = 0 n a n , j D ν N j , α S , N i , α S ω α N i , α S A 3 N j , α S ] = i = 0 n F , N i , α S ω α N i , α S .
Hence
j = 0 n a n , j [ A 1 i = 0 n D 2 N j , α S , N i , α S ω α N i , α S , N k , α S ω α + A 2 i = 0 n a n , j D ν N j , α S , N i , α S ω α N i , α S , N k , α S ω α + A 3 N j , α S , N k , α S ω α ] = i = 0 n G , N i , α S ω α N i , α S , N k , α S ω α , j = 0 n b n , j [ A 1 i = 0 n D 2 N j , α S , N i , α S ω α N i , α S , N k , α S ω α + A 2 i = 0 n a n , j D ν N j , α S , N i , α S ω α N i , α S , N k , α S ω α A 3 N j , α S , N k , α S ω α ] = i = 0 n F , N i , α S ω α N i , α S , N k , α S ω α .
Thus,
j = 0 n a n , j [ A 1 D 2 N j , α S , N k , α S ω α + A 2 D ν N j , α S , N k , α S ω α ] + A 3 a n , k = G , N k , α S ω α , j = 0 n b n , j [ A 1 D 2 N j , α S , N k , α S ω α + A 2 D ν N j , α S , N k , α S ω α ] A 3 b n , k = F , N k , α S ω α .
To put it differently, the coefficients a n , j and b n , j are found by solving the linear systems below:
( A 1 Θ n + A 2 Λ n + A 3 I ) a n = P n , ( A 1 Θ n + A 2 Λ n A 3 I ) b n = M n ,
where
Θ n ( k , j ) : = D 2 N j , α S , N k , α S ω α , Λ n ( k , j ) : = D ν N j , α S , N k , α S ω α , P n ( k ) : = G , N k , α S ω α , M n ( k ) : = F , N k , α S ω α .
or equivalently,
Θ n ( k , j ) : = 0 T 1 α D 2 N j , α S ( τ ) N k , α S ( τ ) ω α ( τ ) d τ , Λ n ( k , j ) : = 0 T 1 α D ν N j , α S ( τ ) N k , α S ( τ ) ω α ( τ ) d τ , P n ( k ) : = 0 T 1 α G ( τ ) N k , α S ( τ ) ω α ( τ ) d τ , M n ( k ) : = 0 T 1 α F ( τ ) N k , α S ( τ ) ω α ( τ ) d τ .

6. Convergence Analysis

Assume that
( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 γ 1 for   some   positive   constant γ 1 , ( A 1 D 2 n + A 2 D ν n A 3 I ) 1 γ 2 for   some   positive   constant γ 2 .
Theorem 2.
The following estimates hold:
ζ n ζ ω α γ ( n + 1 ) ! 2 2 n + 1 f o r   s o m e   p o s i t i v e   c o n s t a n t γ , ξ n ξ ω α γ ( n + 1 ) ! 2 2 n + 1 f o r   s o m e   p o s i t i v e   c o n s t a n t γ .
Proof. 
We have
ζ n ζ = ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 π n S G ( A 1 D 2 + A 2 D ν + A 3 I ) 1 G = ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 π n S G ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 G + ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 G ( A 1 D 2 + A 2 D ν + A 3 I ) 1 G = ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 ( π n S I ) G + ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 [ ( A 1 D 2 + A 2 D ν + A 3 I ) ( A 1 D 2 n + A 2 D ν n + A 3 I ) ] ( A 1 D 2 + A 2 D ν + A 3 I ) 1 G = ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 ( π n S I ) G + A 1 ( D 2 D 2 n ) ζ + A 2 ( D ν D ν n ) ζ .
But
( D 2 D 2 n ) ζ = ( D 2 π n S D 2 π n S ) ζ = ( I π n S ) D 2 ζ + π n S D 2 ( I π n S ) ζ ,
and
( D ν D ν n ) ζ = ( D ν π n S D ν π n S ) ζ = ( I π n S ) D ν ζ + π n S D ν ( I π n S ) ζ .
Thus,
ζ n ζ = ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 [ ( π n S I ) G + A 1 ( I π n S ) D 2 ζ + π n S D 2 ( I π n S ) ζ + A 2 ( I π n S ) D ν ζ + π n S D ν ( I π n S ) ζ ] .
Following [14], it is straightforward to show that for all continuous function φ , we have
φ π n S φ α c φ ( n + 1 ) ! 2 2 n + 1 , for   some   positive   constant c φ .
Since π n S = 1 , and since
ζ n ζ ω α ( A 1 D 2 n + A 2 D ν n + A 3 I ) 1 [ ( π n S I ) G ω α + A 1 ( I π n S ) D 2 ζ ω α + D 2 ( I π n S ) ζ ω α + A 2 ( I π n S ) D ν ζ ω α + D ν ( I π n S ) ζ ω α ] .
Hence, there exist some positive constants c G , c ζ , c ζ and c ζ such that
ζ n ζ ω α γ 1 ( n + 1 ) ! 2 2 n + 1 c G + A 1 c ζ + D 2 c ζ + A 2 c ζ + D ν c ζ .
Letting
γ : = γ 1 c G + A 1 c ζ + D 2 c ζ + A 2 c ζ + D ν c ζ ,
we achieve the desired result.
The second estimate can be proven in the same way as the first one provided. □
Theorem 3.
The following estimates hold:
u n u ω α γ + γ ( n + 1 ) ! 2 2 n + 2 , v n v ω α γ + γ ( n + 1 ) ! 2 2 n + 2 .
Proof. 
Since
u n = ζ n + ξ n 2 , v n = ζ n ξ n 2 ,
we obtain
u n u = ζ n ζ + ξ n ξ 2 , v n v = ζ n ζ + ξ ξ n 2 .
We achieve the intended outcome by applying the preceding theorem. □

7. Numerical Example

The system of Bagley–Torvik equations can be used to simulate the motion of a plate in a Newtonian fluid. In this section, we consider the Bagley–Torvik equation system (1), with the following exact solution
u ( s ) = 1 2 s 2 + s 5 , v ( s ) = 1 2 s 2 s 5 ,
and the following initial conditions:
ρ 0 = 0 , ρ 1 = 0 , η 0 = 0 , η 1 = 0 ,
where
A 1 = 1 , A 2 = 2 , A 3 = 1 , ν = 1.5 .
To demonstrate our method, consider the following results for n = 5 obtained by solving linear systems (7), for the unknowns a 5 , 0 , , a 5 , 5 and b 5 , 0 , , b 5 , 5 , respectively.
a 5 , 0 = 0.34994688584237489950 , a 5 , 1 = 0.24242181863767043857 , a 5 , 2 = 0.24082750802218830743 × 10 1 , a 5 , 3 = 0.35772806140104335114 × 10 2 , a 5 , 4 = 0.10641147349488199344 × 10 2 , a 5 , 5 = 0.47311441812443250151 × 10 3 .
b 5 , 0 = 0.18842256637045510705 , b 5 , 1 = 0.20397238982845960939 , b 5 , 2 = 0.97015810181886046239 × 10 1 , b 5 , 3 = 0.20870155082626365279 × 10 1 , b 5 , 4 = 0.94639740185993567875 × 10 3 , b 5 , 5 = 0.77426268074446688191 × 10 4 .
In this case, the approximate solution ( ζ 5 , ξ 5 ) is
ζ 5 ( s ) = 0.228203152181668994 × 10 2 + 0.35783798758255339381 s 3 2 + 1.4989320129175203445 s 3 1.8554152043491273285 s 9 2 + 1.4844204098105792136 s 6 0.48429356253136554893 s 15 2 ,
ξ 5 ( s ) = 0.6770852731039787 × 10 4 + 0.68959666818352111152 s 9 2 0.079255760895919951972 s 15 2 + 0.0037925790697463120 s 3 2 0.0556763054987918832 s 3 + 0.44154942872731887767 s 6 .
The ways of connecting errors for this example are shown in Table 1.
Figure 1 and Figure 2 compare the exact and approximate solutions at n = 25.

8. Conclusions

The purpose of this work was to develop an efficient projection method and to analyze its convergence for a fractional-order Bagley–Torvik system. We employed fractional shifted Legendre polynomials as the basis for our development. Two significant advantages exist with the proposed method. We began by converting the problem into a system of two separable fractional-order Bagley–Torvik equations, which we then solved using the present method. Second, approximate solutions were obtained by solving two linear equation systems. To demonstrate the method’s accuracy and performance, numerical results are provided. The current method can be updated to use polynomials other than Legendre’s polynomials. In the future, we will investigate this problem using a family of fractional Chebyshev polynomials.

Author Contributions

Conceptualization, S.A. and A.M.; Methodology, S.A. and A.M.; Investigation, S.A. and A.M.; Resources, S.A. and A.M.; Data curation, S.A. and A.M.; Writing—original draft preparation, S.A. and A.M.; writing—review and editing, S.A. and A.M.; supervision, S.A. and A.M.; funding acquisition, S.A.; Writing—original draft, S.A. and A.M.; Writing—review, S.A. and A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The research of the first author was funded by Taif University Researchers Supporting Project number (TURSP-2020/320), Taif University, Taif, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the reviewers and the editor for their insightful comments and suggestions, which significantly improved this original manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bagley, R.L.; Torvik, P.J. On the appearance of the fractional derivative in the behavior of real materials. ASME Trans. J. Appl. Mech. 1984, 51, 294–298. [Google Scholar]
  2. López-Rentería, J.A.; Aguirre-Hernández, B.; Fernández-Anaya, G. LMI stability test for initialized fractional order control systems. Appl. Comput. Math. 2019, 18, 50–61. [Google Scholar]
  3. López-Rentería, J.A.; Aguirre-Hernández, B.; Fernández-Anaya, G. A new guardian map and boundary theorems applied to the stabilization of initialized fractional control systems. Math. Methods Appl. Sci. 2022, 45, 7832–7844. [Google Scholar] [CrossRef]
  4. Dehghan, R. A numerical solution of variable order fractional functional differential equation based on the shifted Legendre polynomials. SeMA J. 2019, 76, 217–226. [Google Scholar] [CrossRef]
  5. Mokhtary, P. Numerical treatment of a well-posed Chebyshev Tau method for Bagley-Torvik equation with high-order of accuracy. Numer. Algor. 2016, 72, 875–891. [Google Scholar] [CrossRef]
  6. Rahimkhani, P.; Ordokhani, Y. Application of Müntz–Legendre polynomials for solving the Bagley–Torvik equation in a large interval. SeMA J. 2018, 75, 517–533. [Google Scholar] [CrossRef]
  7. Ji, T.; Hou, J.; Yang, C. Numerical solution of the Bagley-Torvik equation using shifted Chebyshev operational matrix. Adv. Differ. Equ. 2020, 2020, 648. [Google Scholar] [CrossRef]
  8. Diethelm, K.; Ford, J. Numerical Solution of the Bagley-Torvik Equation. BIT Numer. Math. 2002, 42, 490–507. [Google Scholar] [CrossRef]
  9. Setia, A.; Liu, Y.; Vatsala, A.S. The Solution of the Bagley-Torvik Equation by Using Second Kind Chebyshev Wavelet. In Proceedings of the 11th International Conference on Information Technology: New Generations, Las Vegas, NV, USA, 7–9 April 2014; pp. 443–446. [Google Scholar] [CrossRef]
  10. Mennouni, A. A projection method for solving Cauchy singular integro-differential equations. Appl. Math. Lett. 2012, 25, 986–989. [Google Scholar] [CrossRef] [Green Version]
  11. Mennouni, A. Airfoil polynomials for solving integro-differential equations with logarithmic kernel. Appl. Math. Comput. 2012, 218, 11947–11951. [Google Scholar] [CrossRef]
  12. Samko, S.G.; Kilbas, A.A.; Marichev, D.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993. [Google Scholar]
  13. Mennouni, A. A new efficient strategy for solving the system of Cauchy integral equations via two projection methods. Transylv. J. Math. Mech. 2022, 14, 63–71. [Google Scholar]
  14. Mohammadi, F.; Mohyud-Din, S.T. A fractional-order Legendre collocation method for solving the Bagley-Torvik equations. Adv. Differ. Equ. 2016, 2016, 269. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The approximate solution ζ n for n = 25 .
Figure 1. The approximate solution ζ n for n = 25 .
Symmetry 14 01514 g001
Figure 2. The approximate solution ξ n for n = 25 .
Figure 2. The approximate solution ξ n for n = 25 .
Symmetry 14 01514 g002
Table 1. Numerical Example.
Table 1. Numerical Example.
n ζ ζ n 2 ξ ξ n 2
31.2108 × 10−39.7890 × 10−4
52.5271 × 10−41.3877 × 10−4
78.5470 × 10−5115295 × 10−6
93.7174 × 10−51.8479 × 10−7
111.8893 × 10−54.2876 × 10−8
131.0679 × 10−51.2666 × 10−8
156.5213 × 10−64.4412 × 10−9
174.2231 × 10−61.7705 × 10−9
192.8642 × 10−67.8008 × 10−10
212.0163 × 10−63.7235 × 10−10
231.4638 × 10−61.8982 × 10−10
251.0904 × 10−61.0221 × 10−10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Althubiti, S.; Mennouni, A. An Effective Projection Method for Solving a Coupled System of Fractional-Order Bagley–Torvik Equations via Fractional Shifted Legendre Polynomials. Symmetry 2022, 14, 1514. https://doi.org/10.3390/sym14081514

AMA Style

Althubiti S, Mennouni A. An Effective Projection Method for Solving a Coupled System of Fractional-Order Bagley–Torvik Equations via Fractional Shifted Legendre Polynomials. Symmetry. 2022; 14(8):1514. https://doi.org/10.3390/sym14081514

Chicago/Turabian Style

Althubiti, Saeed, and Abdelaziz Mennouni. 2022. "An Effective Projection Method for Solving a Coupled System of Fractional-Order Bagley–Torvik Equations via Fractional Shifted Legendre Polynomials" Symmetry 14, no. 8: 1514. https://doi.org/10.3390/sym14081514

APA Style

Althubiti, S., & Mennouni, A. (2022). An Effective Projection Method for Solving a Coupled System of Fractional-Order Bagley–Torvik Equations via Fractional Shifted Legendre Polynomials. Symmetry, 14(8), 1514. https://doi.org/10.3390/sym14081514

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