Next Article in Journal
Iterative Numerical Methods for a Fredholm–Hammerstein Integral Equation with Modified Argument
Next Article in Special Issue
A Soliton Solution for the Kadomtsev–Petviashvili Model Using Two Novel Schemes
Previous Article in Journal
Fast Calculation of the Derivatives of Bessel Functions with Respect to the Parameter and Applications
Previous Article in Special Issue
Crossover Dynamics of Rotavirus Disease under Fractional Piecewise Derivative with Vaccination Effects: Simulations with Real Data from Thailand, West Africa, and the US
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Fractional Differential Equations Using Two Different Methods

1
Department of Mathematics, Clarkson University, Potsdam, NY 13699, USA
2
Department of Computer Science and Mathematics, Lebanese American University, Beirut 1102 2801, Lebanon
3
Department of Mathematics, Art and Science Faculty, Siirt University, Siirt 56100, Türkiye
4
Mathematics Research Center, Department of Mathematics, Near East University, Near East Boulevard, Nicosia 99138, Türkiye
5
École Nationale Supérieure des Sciences de la Mer et de l’Aménagement du Littoral (ENSSMAL), Campus Universitaire de Dely Ibrahim, Bois des Cars, B.P. 19, Alger 16320, Algeria
6
Department of Electricity and Electronics, Institute of Research and Development of Processes, Faculty of Science and Technology, University of the Basque Country, 48940 Leioa, Spain
7
Department of Computer Engineering, Biruni University, Istanbul 34010, Türkiye
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(1), 65; https://doi.org/10.3390/sym15010065
Submission received: 17 November 2022 / Revised: 3 December 2022 / Accepted: 20 December 2022 / Published: 26 December 2022

Abstract

:
Numerical methods play an important role in modern mathematical research, especially studying the symmetry analysis and obtaining the numerical solutions of fractional differential equation. In the current work, we use two numerical schemes to deal with fractional differential equations. In the first case, a combination of the group preserving scheme and fictitious time integration method (FTIM) is considered to solve the problem. Firstly, we applied the FTIM role, and then the GPS came to integrate the obtained new system using initial conditions. Figure and tables containing the solutions are provided. The tabulated numerical simulations are compared with the reproducing kernel Hilbert space method (RKHSM) as well as the exact solution. The methodology of RKHSM mainly relies on the right choice of the reproducing kernel functions. The results confirm that the FTIM finds the true solution. Additionally, these numerical results indicate the effectiveness of the proposed methods.

1. Introduction

The fractional calculus’ sense (FC) is presented after classical calculus, but after identifying the limitations of the classical one, many researchers weighed the notions of fractional calculus to comprehend the character systematically. Plenty of mathematicians promoted the vital establishment with the aid of new attributes and corresponding outcomes for FC [1,2,3,4,5,6]. Specifically, the particular functions are proposed to create novel non-integer integral and differential operators. These latter are presented by many people to investigate and symbolize different equations linked with phenomena [7,8,9,10,11,12,13,14,15]. The symmetric and anti-symmetric solitons of the fractional Schrödinger equation have been studied in [16]. In [17], the authors extended the Lie symmetry analysis to the time fractional generalized KdV equations. The Adomian decomposition technique for investigating the fractional KdV–Burgers equation was applied in [18]. Since the Chebyshev collocation technique is implemented for investigating the time-fractional nonlinear Klein–Gordon equation [19], a geometric method is applied for the Korteweg–de Vries equation [20]. A combination of FTIM and geometric method is applied for the fractional Burger–Huxley equation [21]. Additionally, in [22], a lie group approach is implemented for solving the fractional equation. Some other applications can be seen in [23,24,25,26,27,28,29,30,31].
Consider the following equation:
c q C D 0 + , t α u = λ u x x + g ( x , t ) , ( x , t ) Ω ,
u ( x , t ) = H ( x , t ) , o n Γ ,
where c denotes the specific heat, c is the density, λ describes the thermal conductivity coefficient, and C D 0 + , t α shows the Caputo fractional derivative. Equations of this type are used to describe the transport processes with a long memory.The fractional heat equation is one of the most well-known fractional partial differential equations that describe the physical phenomenon. In recent years, solving the fractional heat equation magnetized the attention of mathematicians because of its importance. Many methods have been worked to solve this problem. The higher-order numerical method is used for investigating the fractional heat equation [32]. Additionally, the Laplace homotopy technique is worked to solve this equation [33]. In [34], authors used one-step backward-forward algorithms for multi-dimensional backward heat conduction problems.
Motivated by the above works, in this paper, two numerical methods are worked to solve this equation. One is a combination of a specific type of fictitious time integration technique and the Runge–Kutta method. The other is the RKHSM. The main paper’s contributions are as follows:
  • We present new results on the numerical simulation for the considered equation.
  • We apply two effective numerical methods to obtain these new accurate results.
  • The convergence analysis that confirms the theoretical parts of both methods is discussed.
Variable transformation of a time integration method, namely FTIM, was suggested by Atluri and Liu. Researchers used it for solving linear or nonlinear algebraic equations by defining the fictitious time and using it to derive a system of nonautonomous first-order ordinary differential equations that is equivalent to the original algebraic equations in an n-dimensional space. Some applications of this technique can be seen in [35,36,37,38,39].
In another aspect of this paper, as we mentioned before, we applied the RKHSM for solving the proposed equation.
Recently, the RKHSM has achieved great popularity and success. It became a powerful tool in treating different types of FPDEs, such as the fractional Bloch–Torrey equations [40] and fractional differential equations, including the ABC derivative [41], to name a few. See also [42,43,44,45,46,47] for more research about this method. The RKHSM has many advantages, such as its simplicity and flexibility in treating many fractional differential systems and the fact that it is a mesh-free method. The rest of the paper consists of the following: Section 2 recalls some essential concepts about fractional calculus and reproducing kernel theories. Section 3 and Section 4 are where we see the main theory of the FTIM and RKHSM, respectively, to build a numerical solution for the considered problem. Before finishing with the conclusion part, we validate the proposed methods through two examples.

2. Basic Definitions

Definition 1. 
The left-sided Riemannian–Liouville fractional integral of order μ R + of f C α , α 1 , is
I μ f ( x ) = 1 Γ ( μ ) 0 x f ( η ) ( x η ) 1 μ d η , 0 < μ , x > 0 , I 0 f ( x ) = f ( x ) .
Definition 2. 
We write
f ( x ) C α m , x > 0 , m N 0 ,
provided
f ( m ) C α .
Definition 3. 
Suppose a real function f (with 0 < x ) is in the space C α , α R . By p ( > α ) such that f ( x ) = x p f 1 ( x ) C [ 0 , ] . Obviously, C α C β if β α .
Definition 4. 
Suppose f C 1 m , m N , the Caputo derivative of f is
D μ f ( x ) = [ I m μ ] f ( x ) , μ ( m 1 , m ] , d m d t m f ( x ) , μ = m .
I μ I ν f = I μ + ν , 0 μ , 0 ν , f C α , 0 α .
I μ x ξ f = Γ ( γ + 1 ) Γ ( γ + μ + 1 ) x μ + ξ , 0 < μ , 1 < ξ , 0 < x .
Lemma 1. 
Assume α ( m 1 , m ] and f L 1 [ a , b ] . Then,
J a μ f ( x ) = 1 Γ ( α ) a x ( x t ) 1 μ f ( t ) d t , D a α J a α f ( x ) = f ( x ) ,
and
D a α J a α f ( x ) = f ( x ) k = 1 m 1 f ( k ) ( 0 ) ( x a ) k k ! , x > 0 .
Definition 5. 
The fractional derivative of f in the Caputo sense is
D α f ( x ) = J m α D m f ( x ) = 1 Γ ( m α ) 0 x ( x t ) m α 1 f m ( t ) d t .
Definition 6. 
The Caputo time-fractional derivative operator is
C D t α u ( x , t ) = 1 Γ ( m α ) 0 t ( t z ) m α 1 m u ( x , z ) z m d z , α ( m 1 , m ) , m u ( x , t ) t m , α = m ,
for m to be the smallest integer that exceeds α . The space-fractional derivative with β > 0 is described by
C D x β u ( x , t ) = 1 Γ ( m β ) 0 x ( x s ) m β 1 m u ( s , t ) s m d s , m 1 < β < m , m u ( x , t ) t m , β = m .

Notations

(i)
We will write
A C [ a , b ] = f : [ a , b ] R | f is absolutely continuous on [ a , b ]
to denote the collection of all absolutely continuous functions on [ a , b ] .
(ii)
We write CCF to mean a completely continuous function.
Definition 7. 
Define the function space W 2 3 [ a , b ] by
W 2 3 [ a , b ] = { f ( x ) | f ( j ) A C [ a , b ] , j = 0 , 1 , 2 , f ( 3 ) L 2 [ a , b ] , and f ( a ) = f ( b ) = 0 } .
Definition 8. 
If f , g W 2 3 [ a , b ] , the inner product and norm of this space are described to be
f , g W 2 3 = i = 0 2 f ( i ) ( a ) g ( i ) ( a ) + a b f ( 3 ) ( x ) g ( 3 ) ( x ) d x ,
and
f W 2 3 = f , f W 2 3 .
Theorem 1. 
The RK function of W 2 3 [ 0 , 1 ] is the function R ζ ( x ) described as
R ζ ( x ) = r ( x , ζ ) , x ζ , r ( ζ , x ) , x > ζ ,
with
r ( x , ζ ) = 3 x ζ 13 x ζ 5 156 + 5 x ζ 4 156 5 x ζ 3 78 5 x ζ 2 26 + 21 x 2 ζ 2 104 x 2 ζ 5 624 + 5 x 2 ζ 4 624 5 x 2 ζ 3 312 5 x 2 ζ 26 + 7 x 3 ζ 2 104 x 3 ζ 5 1872 + 5 x 3 ζ 4 1872 5 x 3 ζ 3 936 5 x 3 ζ 78 x 4 ζ 104 + x 4 ζ 5 3744 5 x 4 ζ 4 3744 + 5 x 4 ζ 3 1872 + 5 x 4 ζ 2 624 x 5 ζ 5 18720 + x 5 ζ 4 3744 x 5 ζ 3 1872 x 5 ζ 2 624 x 5 ζ 156 + x 5 120 .
For the proof, see [43].
Definition 9. 
Define the function space W 2 2 [ c , d ] by
W 2 2 [ c , d ] = { f ( t ) | f ( j ) A C [ c , d ] , j = 0 , 1 , f L 2 [ c , d ] , and f ( c ) = 0 } .
Definition 10. 
If f , g W 2 2 [ c , d ] , the inner product and norm are
f , g W 2 2 = i = 0 1 f ( i ) ( c ) g ( i ) ( c ) + c d f ( t ) g ( t ) d t ,
f W 2 2 = f , f W 2 2 .
Theorem 2. 
The RK function of W 2 2 [ c , d ] is the function K η ( t ) defined by
K η ( t ) = 1 3 c 3 + 1 2 c 2 η + c 2 η c + 1 2 c 2 η c c + η t + 1 2 η t 2 1 6 t 3 , t η , 1 3 c 3 + 1 2 c 2 η 1 6 η 3 + c 2 η c + 1 2 c 2 η c + 1 2 η 2 c + η t , t > η .
For the proof, see [43].
Definition 11. 
If f , g W 2 1 [ a , b ] , the inner product and norm are
f , g W 2 1 = f ( a ) g ( a ) + a b f ( x ) g ( x ) d x ,
and
f W 2 1 = f , f W 2 1 .
Theorem 3. 
The RK function of W 2 1 [ a , b ] is the function F ζ ( x ) defined by
F ζ ( x ) = a + 1 + t , t ζ , ζ a + 1 , t > ζ .
For the proof, see [43].
Throughout Ω = [ a , b ] × [ c , d ] .
Definition 12. 
Define the binary function space W 2 ( 3 , 2 ) ( Ω ) by
W 2 ( 3 , 2 ) ( Ω ) = u | 3 x 2 t u i s C C F o n Ω , 5 x 3 t 2 u L 2 ( Ω ) , u ( x , c ) = u ( a , t ) = u ( b , t ) = 0 .
Definition 13. 
If u , v W 2 ( 3 , 2 ) ( Ω ) , the inner product and norm of this space are
u , v W 2 ( 3 , 2 ) = j = 0 1 c d 2 t 2 j x j u ( a , t ) 2 t 2 j x j v ( a , t ) d t + c d 2 t 2 u ( b , t ) 2 t 2 v ( b , t ) d t + j = 0 1 j t j u ( x , c ) , j t j v ( x , c ) W 2 3 + c d a b 3 x 3 2 t 2 u ( x , t ) 3 x 3 2 t 2 v ( x , t ) d x d t ,
and
u W 2 ( 3 , 2 ) = u , u W 2 ( 3 , 2 ) .
Theorem 4. 
The RK function of W 2 ( 3 , 2 ) ( Ω ) is the function
Y ( ζ , η ) ( x , t ) = R ζ ( x ) K η ( t ) .
Definition 14. 
Define the binary function space W 2 ( 1 , 1 ) ( Ω ) by
W 2 ( 1 , 1 ) ( Ω ) = u ( x , t ) | u ( x , t ) i s C C F i n Ω , 2 x t u ( x , t ) L 2 ( Ω ) .
Definition 15. 
If u , v W 2 ( 1 , 1 ) ( Ω ) , the inner product and norm of this space are
u , v W 2 ( 1 , 1 ) = c d t u ( a , t ) t v ( a , t ) d t + u ( x , c ) , v ( x , c ) W 2 1 + c d a b x t u ( x , t ) x t v ( x , t ) d x d t ,
and
u W 2 ( 1 , 1 ) = u , u W 2 ( 1 , 1 ) .
Theorem 5. 
The RK function of W 2 ( 1 , 1 ) ( Ω ) is the function
T ( ζ , η ) ( x , t ) = F ζ ( x ) F η ( t ) .

3. The Fictitious Time Integration Method (FTIM)

Here, the FTIM is implemented. Consider the following equation:
c q C D 0 + , t α u + b ( x ) u x + u x x + g ( x , t ) = 0 , ( x , t ) Ω .
Using (10), 0 < α < 1 and 0 < β 1 for Equation (28), we have
c q Γ ( 1 α ) 0 t u σ ( x , σ ) ( t σ ) α d σ b ( x ) u x + u x x + g ( x , t ) = 0 .
To increase the stability of the technique, we propose a fictitious damping coefficient μ in Equation (29) by
μ c q Γ ( 1 α ) 0 t u σ ( x , σ ) ( t σ ) α d σ μ b ( x ) u x + μ u x x + μ g ( x , t ) = 0 .
Imposing, in Equation (30), the transformation
ω ( x , t , ξ ) = ( 1 + θ ) λ u ( x , t ) , 0 < λ 1 ,
yields
μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t ω σ ( x , σ , θ ) ( t σ ) α d σ + b ( x ) ω x ( x , t , θ ) + ω x x ( x , t , θ ) ) + μ g ( x , t ) = 0 .
We consider
ω θ = λ ( 1 + θ ) λ 1 u ( x , t ) .
Equation (32) will be
ω θ = μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t ω σ ( x , σ , θ ) ( t σ ) α d σ + b ( x ) ω x ( x , t , θ ) + ω x x ( x , t , θ ) ) + μ g ( x , t ) + λ ( 1 + θ ) λ 1 u .
Equation (34) can be converted to a new class of PDE for ω by choosing u = ω / ( 1 + θ ) λ :
ω θ = μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t ω σ ( x , σ , θ ) ( t σ ) α d σ + b ( x ) ω x ( x , t , θ ) + ω x x ( x , t , θ ) ) + μ g ( x , t ) + λ ω ( x , t , θ ) 1 + θ .
Using
θ ω ( 1 + θ ) λ = ω θ ( 1 + θ ) λ λ ω ( 1 + θ ) λ + 1 ,
Implementing 1 / ( 1 + θ ) λ for Equation (35), one obtains
θ ω ( 1 + θ ) λ = μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t ω σ ( x , σ , θ ) ( t σ ) α d σ + b ( x ) ω x ( x , t , θ ) + ω x x ( x , t , θ ) ) + μ g ( x , t ) .
By u = ω ( 1 + θ ) λ , we obtain
u θ = μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t u σ ( x , σ , θ ) ( t σ ) α d σ + b ( x ) u x ( x , t , θ ) + u x x ( x , t , θ ) ) + μ g ( x , t ) .
Suppose that u i j ( θ ) : = u ( x i , t j , θ ) as the discrete values of u at a grid point ( x i , t j ) , and Equation (38) converts to
d d ξ u i j ( θ ) = μ ( 1 + θ ) λ ( c q Γ ( 1 α ) 0 t j u σ ( x i , σ , θ ) ( t j σ ) α d σ + b ( x i ) u x ( x i , t j , θ ) + u x x ( x i , t j , θ ) ) + μ g ( x i , t j ) ,
where
0 t j u σ ( x i , σ , θ ) ( t j σ ) α d σ l = 1 j 1 u ( x i , t l + 1 , θ ) u ( x i , t l , θ ) Δ t ( t j t l ) α ,
where t j = j Δ t , x i = a + i Δ x , and Δ t = T n .
u = ( u 1 1 , u 1 2 , , u m n ) T , Equation (39) can be abstracted by
u = Q ( u , ξ ) , u R m × n , ξ R , M = m × n ,
where Q R M is a vector-valued function of u and θ and u is an M-dimensional vector. Now, we implement the group-preserving scheme (GPS) [48] to solve Equation (39) as
u s + 1 = u s + cosh Δ θ Q s u s 1 Q s . u s + sinh Δ θ Q s u s u s Q s Q s 2 Q s = u s + Ξ s Q s .
Now, we employ the GPS by taking the initial value of u i j ( 0 ) to solve Equation (39) from the initial fictitious time θ = 0 to a selected final fictitious time θ f . Additionally, the terminating criterion for this method is
i , j = 1 m , n [ u i j ( s + 1 ) u i j ( s ) ] 2 ε ,
where ε is a picked convergence criterion. The solution of u will be obtained by
u i j = u i j ( θ 0 ) ( 1 + θ 0 ) λ ,
where θ 0 ( θ f ) satisfies the above criterion.

4. The Application of RKHSM

4.1. Methodology for RKHSM

In this part, the RKHSM is used to solve Problems (1) and (2), using the following way:
Step 1: Considering the following transformation
v = u P ,
where P ( x , t ) = f 2 ( 0 ) f 2 ( t ) x b f 1 ( 0 ) f 1 ( t ) 1 x b + u 0 ( x ) for which u ( 0 , t ) = f 1 ( t ) , u ( b , t ) = f 2 ( t ) , and u ( x , 0 ) = u 0 ( x ) .
Consequently, the new form of (1) and (2) is as follows:
c q C D 0 + , t α v λ v x x = ( x , t ) , 0 x b , 0 t d ,
With
v ( x , t ) = 0 o n Γ ,
where x , t = λ P x x ( x , t ) c q C D 0 + , t α P ( x , t ) + g ( x , t ) .
Step 2: Defining a linear operator T : W 2 ( 3 , 2 ) ( Ω ) W 2 ( 1 , 1 ) ( Ω ) as follows
T : W 2 ( 3 , 2 ) ( Ω ) W 2 ( 1 , 1 ) ( Ω ) v c q C D 0 + , t α v λ v x x .
Lemma 2. 
The operator T is a bounded linear.
Proof. 
We begin by checking directly that T is bounded. So, we must show that
T v W 2 ( 1 , 1 ) C v W 2 ( 3 , 2 ) , with C > 0 .
We have
T v ( x , t ) W 2 ( 1 , 1 ) 2 = T v ( x , t ) , T v ( x , t ) W 2 ( 1 , 1 ) = c d t T v ( a , t ) 2 d t + T v ( x , c ) , T v ( x , c ) W 2 1 + c d a b x t T v ( x , t ) 2 d x d t = T v ( a , c ) 2 + a b x T v ( x , c ) 2 d x + c d t T v ( a , t ) 2 d t + c d a b x t T v ( x , t ) 2 d x d t .
In view of reproducing the property,
v ( x , t ) = v ( , ) , Y ( x , t ) ( , ) W 2 ( 3 , 2 ) .
In a similar way, we deduce
i x i j t j T v ( x , t ) = v ( , ) , i x i j t j T Y ( x , t ) ( , ) W 2 ( 3 , 2 ) , i , j { 0 , 1 } .
Applying the Schwarz inequality, we discover
i x i j t j T v ( x , t ) = v ( , ) , i x i j t j T Y ( x , t ) ( , ) W 2 ( 3 , 2 ) v W 2 ( 3 , 2 ) i x i j t j T Y ( x , t ) ( , ) W 2 ( 3 , 2 ) .
Since Y ( x , t ) ( , ) is continuous, we consequently have
i x i j t j T v ( x , t ) C i , j v W 2 ( 3 , 2 ) , i , j { 0 , 1 } .
Hence
T v ( x , t ) W 2 ( 1 , 1 ) 2 C 0 , 0 2 v W 2 ( 3 , 2 ) 2 + a b C 1 , 0 2 v W 2 ( 3 , 2 ) 2 d ξ + c d C 0 , 1 2 v W 2 ( 3 , 2 ) 2 d τ + c d a b C 1 , 1 2 v W 2 ( 3 , 2 ) 2 d x d t C 0 , 0 2 + C 1 , 0 2 b a + C 0 , 1 2 d c + C 1 , 1 2 b a d c v W 2 ( 3 , 2 ) 2 .
Therefore,
T v ( x , t ) W 2 ( 1 , 1 ) 2 C v W 2 ( 3 , 2 ) 2 ,
where C = C 0 , 0 2 + C 1 , 0 2 b a + C 0 , 1 2 d c + C 1 , 1 2 b a d c .
We apply, next, the operator T in order to reformulate the problem (46) and (47) to be
T v ( x , t ) = x , t , ( x , t ) Ω , v ( x , t ) = 0 on Γ ,
where x , t = λ P x x ( x , t ) c q C D 0 + , t α P ( x , t ) + g ( x , t ) .
Step 3: Construct the { Θ ¯ i } i = 1 on W 2 ( 3 , 2 ) ( Ω ) , providing this by using the Gram–Schmidt process:
Θ ¯ i ( x , t ) = k = 1 i i k Θ k ( x , t ) , 0 < i i , i = 1 , 2 , ,
where
  • Θ i ( x , t ) = T * ϱ i ( x , t ) , in which T * denotes the adjoint of T and ϱ i ( x , t ) = T ( x i , t i ) ( x , t ) where T ( x i , t i ) ( x , t ) is given by (27).
  • The countable set { x i , t i } i = 1 is dense in Ω .
  • { Θ i } i = 1 is a function system in W 2 ( 3 , 2 ) ( Ω ) and the following shows the way that we can construct it:
    Θ i ( x , t ) = T * ϱ i ( x , t ) = T * ϱ i ( ζ , η ) , Y ( x , t ) ( ζ , η ) W 2 ( 3 , 2 ) = ϱ i ( ζ , η ) , T ( ζ , η ) Y ( x , t ) ( ζ , η ) W 2 ( 1 , 1 ) = T ( ζ i , η i ) ( ζ , η ) , T ( ζ , η ) Y ( x , t ) ( ζ , η ) W 2 ( 1 , 1 ) = T ( ζ , η ) Y ( x , t ) ( ζ , η ) | ( ζ , η ) = ( x i , t i ) = T ( ζ , η ) Y ( ζ , η ) ( x , t ) | ( ζ , η ) = ( x i , t i ) = c q C D 0 + , η α Y ( ζ , η ) ( x , t ) λ ζ 2 2 Y ( ζ , η ) ( x , t ) ( ζ , η ) = ( x i , t i ) .
  • i k is the orthogonalization coefficients which are defined by
    i j = 1 Θ 1 , for i = j = 1 , 1 ς i , for i = j 1 , 1 ς i k = j i 1 C i k k j , for i > j ,
    where ς i = Θ i 2 k = 1 i 1 C i k 2 , C i k = Θ i , Θ ¯ k W 2 ( 3 , 2 ) .
Theorem 6. 
Assume { x i , t i } i = 1 is dense; therefore, { Θ i } i = 1 is the complete system of W 2 ( 3 , 2 ) ( Ω ) .
Proof. 
Clearly, Θ i ( x , t ) W 2 ( 3 , 2 ) ( Ω ) . Thus, for v ( x , t ) W 2 ( 3 , 2 ) ( Ω ) ,
v ( x , t ) , Θ i ( x , t ) W 2 ( 3 , 2 ) = 0 , i = 1 , 2 , .
as
v , Θ i W 2 ( 3 , 2 ) = v ( x , t ) , T * ϱ i ( x , t ) W 2 ( 3 , 2 ) = T v ( x , t ) , ϱ i ( x , t ) W 2 ( 1 , 1 ) = T v ( x i , t i ) = 0 ,
and due to the density of { x i , t i } i = 1 in Ω :
T v ( x , t ) = 0 .
by applying T 1 ,
v ( x , t ) = 0 .
Step 4: The solution’s representation is given by
Theorem 7. 
Assume { x i , t i } i = 1 is a dense set on Ω and (57) has a unique solution on W 2 ( 3 , 2 ) ( Ω ) , then
v ( x , t ) = i = 1 k = 1 i i k ( x k , t k ) Θ ¯ i ( x , t )
is the solution of (57), and the solution of (1) and (2) is
u ( x , t ) = i = 1 k = 1 i i k ( x k , t k ) Θ ¯ i ( x , t ) P ( x , t ) .
Proof. 
We know that the basis Θ ¯ i ( x , t ) i = 1 is a complete orthonormal system in the space W 2 ( 3 , 2 ) ( Ω ) , then
v ( x , t ) = i = 1 v ( x , t ) , Θ ¯ i ( x , t ) W 2 ( 3 , 2 ) Θ ¯ i ( x , t ) = i = 1 k = 1 i i k v ( x , t ) , Θ k ( x , t ) W 2 ( 3 , 2 ) Θ ¯ i ( x , t ) = i = 1 k = 1 i i k v ( x , t ) , T * ρ k ( x , t ) W 2 ( 3 , 2 ) Θ ¯ i ( x , t ) = i = 1 k = 1 i i k T v ( x , t ) , ϱ k ( x , t ) W 2 ( 1 , 1 ) Θ ¯ i ( x , t ) = i = 1 k = 1 i i k T v ( x , t ) , T ( x k , t k ) ( x , t ) W 2 ( 1 , 1 ) Θ ¯ i ( x , t ) = i = 1 k = 1 i i k ( x k , t k ) Θ ¯ i ( ξ , τ ) ,
with ( x k , t k ) = T v ( x k , t k ) .
On the other hand, (66) follows directly from u = v P ( x , t ) .

Remarks

1.
We have
v n ( x , t ) = i = 1 n k = 1 i i k ( x k , t k ) Θ ¯ i ( x , t ) .
2.
W 2 ( 3 , 2 ) ( Ω ) is a Hilbert space. Then, we deduce
i = 1 k = 1 i i k ( x k , t k ) Θ ¯ i ( x , t ) < .

4.2. Convergence Analysis

The approximate solution of (65) takes the form
v n ( x , t ) = i = 1 n Λ i Θ ¯ i ( x , t ) ,
where
Λ i = k = 1 i i k ( x k , t k ) .
Here, by letting ( x 1 , t 1 ) = ( a , c ) , it is possible to know the values of v ( x 1 , t 1 ) from the initial and boundary conditions. In addition, v 0 ( x 1 , t 1 ) = v ( x 1 , t 1 ) .
Theorem 8. 
Suppose that v n W 2 ( 3 , 2 ) is bounded in (57), { ( x i , t i ) } i = 1 is dense on Ω , and the solution of (70) is unique. Then, v n converges to v and
v n ( x , t ) = i = 1 n Λ i Θ ¯ i ( x , t ) .
Proof. 
(i)
From (70), we know
v n + 1 = v n + Λ n + 1 Θ ¯ n + 1 .
then,
v n + 1 W 2 ( 3 , 2 ) 2 = v n + Λ n + 1 Θ ¯ n + 1 W 2 ( 3 , 2 ) 2 = v n + Λ n + 1 Θ ¯ n + 1 , v n + Λ n + 1 Θ ¯ n + 1 W 2 ( 3 , 2 ) = v n , v n W 2 ( 3 , 2 ) + v n , Λ n + 1 Θ ¯ n + 1 W 2 ( 3 , 2 ) + Λ n + 1 Θ ¯ n + 1 , v n W 2 ( 3 , 2 ) + Λ n + 1 Θ ¯ n + 1 , Λ n + 1 Θ ¯ n + 1 W 2 ( 3 , 2 ) = v n W 2 ( 3 , 2 ) 2 + v n , Λ n + 1 Θ ¯ n + 1 W 2 ( 3 , 2 ) + Λ n + 1 Θ ¯ n + 1 , v n W 2 ( 3 , 2 ) + , Λ n + 1 2 Θ ¯ n + 1 , Θ ¯ n + 1 W 2 ( 3 , 2 ) .
Now, the orthogonality of Θ ¯ i ( x , t ) i = 1 implies
v n + 1 W 2 ( 3 , 2 ) 2 = v n W 2 ( 3 , 2 ) 2 + Λ n + 1 2 = v n 1 W 2 ( 3 , 2 ) 2 + Λ n 2 + Λ n + 1 2 = v n 2 W 2 ( 3 , 2 ) 2 + Λ n 1 2 + Λ n 2 + Λ n + 1 2 = v 1 W 2 ( 3 , 2 ) 2 + Λ 2 2 + Λ 3 2 + + Λ n 2 + Λ n + 1 2 = v 0 W 2 ( 3 , 2 ) 2 + i = 1 n + 1 Λ i 2 ,
Hence,
v n W 2 ( 3 , 2 ) v n + 1 W 2 ( 3 , 2 ) .
The convergence of v n W 2 ( 3 , 2 ) follows directly from the boundedness of v n W 2 ( 3 , 2 ) . So, there exists ϝ such that
i = 1 Λ i 2 = ϝ ,
where the constant ϝ is positive.
As a result,
Λ i 2 i = 1 2 .
As v m v m 1 v n + 1 v n and for n < m , we write
v m v n W 2 ( 3 , 2 ) 2 = v m v m 1 + v m 1 + v n + 1 v n W 2 ( 3 , 2 ) 2 = v m v m 1 W 2 ( 3 , 2 ) 2 + v m 1 v m 2 W 2 ( 3 , 2 ) 2 + + v n + 1 v n W 2 ( 3 , 2 ) 2 .
Furthermore,
v m v m 1 W 2 ( 3 , 2 ) 2 = Λ m 2 .
Thus,
v m v n W 2 ( 3 , 2 ) 2 = p = n + 1 m Λ p 2 0 , as n , m .
The completeness of W 2 ( 3 , 2 ) ( Ω ) allows us to deduce that v n v ˜ as n .
(ii)
To prove this, let us take the limits in (70)
v ˜ ( x , t ) = i = 1 Λ i Θ ¯ i ( x , t ) .
We apply the linear operator T to (81)
T v ˜ ( x , t ) = i = 1 Λ i T Θ ¯ i ( x , t ) ,
Hence,
T v ˜ ( x p , t p ) = i = 1 Λ i T Θ ¯ i ( x , t ) , ϱ p ( x , t ) W 2 ( 1 , 1 ) = i = 1 Λ i Θ ¯ i ( x , t ) , T * ϱ p ( x , t ) W 2 ( 3 , 2 ) = i = 1 Λ i Θ ¯ i ( x , t ) , Θ p ( x , t ) W 2 ( 3 , 2 ) .
Thus,
j p T v ˜ ( x p , t p ) = j p i = 1 Λ i Θ ¯ i ( x , t ) , Θ p ( x , t ) W 2 ( 3 , 2 ) ,
and we take the summation p = 1 j to deduce
p = 1 j j p T v ˜ ( x p , t p ) = i = 1 Λ i Θ ¯ i ( x , t ) , p = 1 j j p Θ p ( x , t ) W 2 ( 3 , 2 ) = i = 1 Λ i Θ ¯ i ( x , t ) , Θ ¯ j ( x , t ) W 2 ( 3 , 2 ) = Λ j .
Observe then from (71) that
T v ˜ ( x p , t p ) = ( x p , t p ) .
For all ( ζ , η ) Ω , it exists { ( x q j , t q j ) } j = 1 such that ( x q j , t q j ) ( ζ , η ) , as j .
It is well-known to us that
T v ˜ ( x q j , t q j ) = ( x q j , t q j ) .
Using the continuity of and letting j allows us to
T v ˜ ( ζ , η ) = ( ζ , η ) .

5. Numerical Experiments

We apply the proposed methods to solve some problems. In Example 1, we use RKHSM to solve the considered equation, and the GPS is considered for Example 2 to deal with the fractional convection–diffusion equation. Now, how to apply the RKHSM can be summarized in the following procedure:
  • Step 1: Setting n = p × q ;
  • Step 2: Setting Θ i ( x i , t i ) = T ( ζ , η ) Y ( x , t ) ( ζ , η ) | ( ζ , η ) = ( x i , t i ) ;
  • Step 3: Calculating the orthogonalization coefficients i j using (60);
  • Step 4: Setting Θ ¯ i ( x i , t i ) = k = 1 i i k Θ k ( x i , t i ) , i = 1 , 2 , , n ;
  • Step 5: Choosing an initial guess u 0 ( x 1 , t 1 ) ;
  • Step 6: Setting i = 1 ;
  • Step 7: Setting Λ i = k = 1 i i k ( x k , t k ) ;
  • Step 8: u i ( x i , t i ) = = 1 i Λ Θ ¯ ( x , t ) ;
  • Step 9: If i < n , set i = i + 1 . Go to step 7. Else stop, Where x i = i p , i = 1 , 2 , , p and t j = j q , j = 1 , 2 , , q . n is the grid points’ number.
Example 1. 
Considering the following problem with the fractional order α = 0.5 :
C D 0 + , t α u = u x x + g ( x , t ) , ( x , t ) [ 0 , 1 ] 2 ,
where
u ( x , t ) = e x p ( x ) x 2 ( 1 x ) 2 t α ,
and
g ( x , t ) = e x x x 3 + 6 x 2 + x 8 + 2 t α π ( x 1 ) 2 x 2 csc ( π α ) Γ ( α ) .
In this example, the RKHSM is tested with the standard grid points x i = i p , i = 1 , , p and t j = j q , j = 1 , , q with p × q = n = 100 . The comparison of (89) with (1) and (2) shows Ω = [ 0 , 1 ] × [ 0 , 1 ] , c q = 1 , λ = 1 , and u ( 0 , t ) = u ( 1 , t ) = u ( x , 0 ) = 0 . Therefore, as we see in Section 4, the approximate solution of (89) takes the form
v n ( x , t ) = i = 1 n k = 1 i i k g ( x k , t k ) Θ ¯ i ( x , t ) .
In Table 1, a numerical comparison between the obtained results via RKHSM with the exact solution, for α = 0.9, 0.8, 0.75, is given. These results clearly show that the approximate solution (using the RKHSM) converges to the exact solution. The results are in good agreement with each other, and this confirms the effectiveness of the RKHSM to solve this type of equation.
Example 2. 
Consider the problem (28) with b ( x ) = x and fractional order α = 0.01 where
u ( x , t ) = x 2 t 2 α 2 Γ ( α + 1 ) Γ ( 2 α + 1 ) ,
and
g ( x , t ) = 2 t 2 + 2 x 2 + 2 .
We solve this problem by using the GPS with taking Δ θ = 1 e 10 , m = n = 20 , and λ = 1 . Additionally, an initial guess of u i j ( 0 ) = 1 e 1 is taken. Figure 1 shows the numerical solution, exact solution, absolute error, and absolute error’s contour. Indeed, we present in Table 2 the values of absolute errors between the numerical solution (using the GPS) and the exact solution for Example 2. From this table’s results, it is clear that the error estimate confirms the accuracy of this new method, and Figure 1 shows that both graphs are very similar in their behavior.

6. Conclusions

In the current work, we successfully implemented two numerical schemes to gain approximate solutions to the considered problems. One is the FITM, which converted the original problem into a new one with one extra dimension. After that, we used GPS to solve the problem. The other is the RKHSM, which was used for the mentioned problem. The main steps for applying this method are defining an appropriate bounded linear operator and constructing an orthonormal function system of the appropriate RKHS. Indeed, the both methods are shown to have good convergence. Two examples were employed to show the capacity and reliability of the FITM and RKHSM. Our obtained results are compared with exact results and they are found to be in good agreement with each other. From the numerical results, it can be observed the suitability, ease, and effectiveness of the proposed approaches for solving such types of fractional partial differential equation. This research opens the way for the use of the two proposed methods to study the mentioned problem for various new fractional derivatives. As part of our purpose, we plan to apply the FTIM and RKHSM to multidimensional fractional partial differential equations, which will be new in the literature.

Author Contributions

Conceptualization, E.K.A.; Methodology, M.P. and N.A.; Software, A.A.; Formal analysis, N.A.; Investigation, M.P.; Resources, M.D.l.S.; Data curation, A.A.; Writing—original draft, N.A.; Writing—review & editing, N.A.; Visualization, A.A. and E.K.A.; Supervision, M.B.; Project administration, A.A., M.D.l.S. and M.B.; Funding acquisition, M.D.l.S. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

Basque Government, Grants IT1555-22 and KK-2022/00090 MCIN/AEI 269.10.13039/ 501100011033, Grant PID2021-1235430B-C21/C22.

Data Availability Statement

Data are included within this research.

Acknowledgments

We are thankful to the anonymous reviewers for carefully reading the paper and their helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liouville, J. Mémoire sur quelques questions de géométrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces questions. JEC Polytech. 1832, 13, 1–69. [Google Scholar]
  2. Riemann, G.F.B. Versuch Einer Allgemeinen Auffassung der Integration und Differentiation; Gesammelte Mathematische Werke: Leipzig, Germany, 1896. [Google Scholar]
  3. Caputo, M. Elasticita e Dissipazione; Zanichelli: Bologna, Italy, 1969. [Google Scholar]
  4. Miller, K.S.; Ross, B. An Introduction to Fractional Calculus and Fractional Diffrential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  5. Podlubny, I. Fractional Diffrential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  6. Singh, H.; Kumar, D.; Baleanu, D. Methods of Mathematical Modelling: Fractional Diffrential Equations; CRC Press Taylor & Francis Group: New York, NY, USA, 2019. [Google Scholar]
  7. Baleanu, D.; Guvenc, Z.B.; Tenreiro Machado, J.A. New Trends in Nanotechnology and Fractional Calculus Applications; Springer: Dordrecht, The Netherlands, 2010. [Google Scholar]
  8. Atangana, A.; Alkahtani, B.T. Analysis of non-homogenous heat model with new trend of derivative with fractional order. Chaos Solit. Fract. 2016, 89, 566–571. [Google Scholar]
  9. Kumar, D.; Singh, J.; Tanwar, K.; Baleanu, D. A new fractional exothermic reactions model having constant heat source in porous media with power, exponential and Mittag-Leffler laws. Int. Heat Mass Transf. 2019, 138, 1222–1227. [Google Scholar] [CrossRef]
  10. Baleanu, D.; Wu, G.C.; Zeng, S.D. Chaos analysis and asymptotic stability of generalized Caputo fractional diffrential equations. Chaos Solit. Fract. 2017, 102, 99–105. [Google Scholar] [CrossRef]
  11. Veeresha, P.; Prakasha, D.G.; Baskonus, H.M. New numerical surfaces to the mathematical model of cancer chemotherapy effct in Caputo fractional derivatives. Chaos 2019, 29, 013119. [Google Scholar] [CrossRef]
  12. Gao, W.; Veeresha, P.; Baskonus, H.M.; Prakasha, D.G.; Kumar, P. A new study of unreported cases of 2019-nCOV epidemic outbreaks. Chaos Solit. Fract. 2020, 138, 109929. [Google Scholar] [CrossRef]
  13. Kumar, D.; Singh, J.; Baleanu, D. On the analysis of vibration equation involving a fractional derivative with Mittag-Leffler law. Math. Meth. Appl. Sci. 2020, 43, 443–457. [Google Scholar] [CrossRef]
  14. Veeresha, P.; Prakasha, D.G.; Baleanu, D. Analysis of fractional Swift-Hohenberg equation using a novel computational technique. Math. Meth. Appl. Sci. 2020, 43, 1970–1987. [Google Scholar] [CrossRef]
  15. Caputo, M.; Fabrizio, M. A new defiition of fractional derivative without singular kernel. Prog. Fract. Diffr. Appl. 2015, 1, 73–85. [Google Scholar]
  16. Cao, Q.-H.; Dai, C.-Q. Symmetric and anti-symmetric solitons of the fractional second- and third-order nonlinear Schrödinger equation. Chin. Phys. Lett. 2021, 38, 090501. [Google Scholar] [CrossRef]
  17. Chen, C.; Jiang, Y.-L.; Wang, X.-T. Lie symmetry analysis of the time fractional generalized KdV equations with variable coefficients. Symmetry 2019, 11, 1281. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, Q. Numerical solution for fractional KdV-Burgers equation by Adomian decomposition method. Appl. Math. Comput. 2006, 182, 1048–1055. [Google Scholar]
  19. Nagy, A.M. Numerical solution of time fractional nonlinear Klein-Gordon equation using Sinc Chebyshev collocation method. Appl. Math. Comput. 2017, 310, 139–148. [Google Scholar] [CrossRef]
  20. Parto-Haghighi, M.; Akinlar, M.A.; Chu, Y. New numerical solutions of fractional-order Korteweg-de Vries equation. Results Phys. 2020, 19, 103326. [Google Scholar]
  21. Inc, M.; Partohaghighi, M.; Akinlar, M.A.; Agarwal, P.; Chu, Y.-M. New solutions of fractional-order Burger-Huxley equation. Results Phys. 2020, 18, 103290. [Google Scholar] [CrossRef]
  22. Hashemi, M.S.; Baleanu, D.; Parto-Haghighi, M. A lie group approach to solve the fractional Poisson equation. Rom. J. Phys. 2015, 60, 1289–1297. [Google Scholar]
  23. Akgül, A.; Ahmad, H. Reproducing kernel method for Fangzhu’s oscillator for water collection from air. Math. Meth. Appl. Sci. 2020, 1–10. [Google Scholar] [CrossRef]
  24. Ahmad, H.; Akgül, A.; Khan, T.A.; Stanimirović, P.S.; Chu, Y.-M. New perspective on the conventional solutions of the nonlinear time-fractional partial differential equations. Complexity 2020, 2020, 8829017. [Google Scholar] [CrossRef]
  25. Inc, M.; Khan, M.N.; Ahmad, I.; Yao, S.W.; Ahmad, H.; Thounthong, P. Analysing time-fractional exotic options via efficient local meshless method. Results Phys. 2020, 19, 103385. [Google Scholar] [CrossRef]
  26. Ahmad, I.; Khan, M.N.; Inc, M.; Ahmad, H.; Nisar, K.S. Numerical simulation of simulate an anomalous solute transport model via local meshless method. Alex. Eng. J. 2020, 59, 2827–2838. [Google Scholar] [CrossRef]
  27. Hussain, M.; Haq, S. Weighted meshless spectral method for the solutions of multi-term time fractional advection-diffusion problems arising in heat and mass transfer. Int. J. Heat Mass Transf. 2019, 129, 1305–1316. [Google Scholar] [CrossRef]
  28. Ahmad, I.; Ahmad, H.; Inc, M.; Yao, S.-W.; Almohsen, B. Application of local meshless method for the solution of two term time fractional-order multi-dimensional PDE arising in heat and mass transfer. Therm. Sci. 2020, 24, 95–105. [Google Scholar] [CrossRef]
  29. Pitolli, F.; Sorgentone, C.; Pellegrino, E. Approximation of the Riesz-Caputo derivative by cubic splines. Algorithms 2022, 15, 69. Available online: https://www.mdpi.com/1999-4893/15/2/69 (accessed on 15 February 2022). [CrossRef]
  30. Izadi, M.; Srivastava, H.M. A discretization approach for the nonlinear fractional logistic equation. Entropy 2020, 22, 1328. Available online: https://www.mdpi.com/1099-4300/22/11/1328 (accessed on 16 November 2020). [CrossRef]
  31. Shymanskyi, V.; Sokolovskyy, Y. Finite element calculation of the linear elasticity problem for biomaterials with fractal structure. Open Bioinform. J. 2021, 14, 114–122. [Google Scholar] [CrossRef]
  32. Priya, G.S.; Prakash, P.; Nieto, J.J.; Kayar, Z. Higher-order numerical scheme for the fractional heat equation with Dirichlet and Neumann boundary conditions. Numer. Heat Transf. 2013, 63, 540–559. [Google Scholar] [CrossRef]
  33. Mehmet, Y.; Necati, O. Numerical inverse Laplace homotopy technique for fractional heat equations. Therm. Sci. 2018, 22, 185–194. [Google Scholar]
  34. Chen, Y.-W. A highly accurate backward-forward algorithm for multi-dimensional backward heat conduction problems in fictitious time domains. Int. J. Heat Mass Transf. 2018, 120, 499–514. [Google Scholar] [CrossRef]
  35. Partohaghighi, M.; Akinlar, M.A.; Weber, G.W. New solutions of hyperbolic telegraph equation. J. Dyn. Games 2021, 8, 129–138. [Google Scholar]
  36. Partohaghighi, M.; Bayram, M.; Baleanu, D. On numerical solution of the time fractional advection-diffusion equation involving Atangana-Baleanu-Caputo derivative. Open Phys. 2019, 17, 816–822. [Google Scholar] [CrossRef]
  37. Chen, Y.-W.; Liu, C.-S.; Chang, Y.-S.; Chang, J.-R. Constraint-type fictitious time integration method for solving non-linear multi-dimensional elliptic partial differential equations. J. Mar. Sci. Technol. 2020, 28, 168–178. [Google Scholar]
  38. Chen, Y.-W. High order implicit and explicit Lie-group schemes for solving backward heat conduction problems. Int. J. Heat Mass Transf. 2016, 101, 1016–1029. [Google Scholar] [CrossRef]
  39. Liu, C.S. Cone of non-linear dynamical system and group preserving schemes. Int. J. Non Linear Mech. 2021, 36, 1047–1068. [Google Scholar] [CrossRef]
  40. Akgül, E.K. A novel method for the space and time fractional Bloch-torrey equations. Therm. Sci. 2018, 22, S253–S258. [Google Scholar] [CrossRef]
  41. Akgül, A. A novel method for a fractional derivative with non-local and non-singular kernel. Chaos Solit. Fract. 2018, 114, 478–482. [Google Scholar] [CrossRef]
  42. Hemati, F.; Ghasemi, M.; Khoshsiar Ghaziani, R. Numerical solution of the multiterm time-fractional diffusion equation based on reproducing kernel theory. Numer. Methods Partial Differ. Equ. 2021, 37, 44–68. [Google Scholar] [CrossRef]
  43. Attia, N.; Akgül, A.; Seba, D.; Nour, A. Reproducing kernel Hilbert space method for the numerical solutions of fractional cancer tumor models. Math. Meth. Appl. Sci. 2020, 1–22. [Google Scholar] [CrossRef]
  44. Abu Arqub, O.; Shawagfeh, N. Solving optimal control problems of Fredholm constraint optimality via the reproducing kernel Hilbert space method with error estimates and convergence analysis. Math. Methods Appl. Sci. 2021, 44, 7915–7932. [Google Scholar] [CrossRef]
  45. Abu Arqub, O.; Shawagfeh, N. Application of reproducing kernel algorithm for solving Dirichlet time-fractional diffusion-Gordon types equations in porous media. J. Porous Media 2019, 22, 411–434. [Google Scholar] [CrossRef]
  46. Djennadi, S.; Shawagfeh, N.; Abu Arqub, O. A numerical algorithm in reproducing kernel-based approach for solving the inverse source problem of the time-space fractional diffusion equation. Partial Differ. Equ. Appl. Math. 2021, 4, 100164. [Google Scholar] [CrossRef]
  47. Abu Arqub, O. Reproducing kernel algorithm for the analytical-numerical solutions of nonlinear systems of singular periodic boundary value problems. Math. Prob. Eng. 2015, 2015, 518406. [Google Scholar] [CrossRef]
  48. Cui, M.; Lin, Y. Nonlinear Numerical Analysis in the Reproducing Kernel Space; Nova Science Publishers Inc.: New York, NY, USA, 2009. [Google Scholar]
Figure 1. Solution under applying GPS for the fractional convection-diffusion equation in Example 2.
Figure 1. Solution under applying GPS for the fractional convection-diffusion equation in Example 2.
Symmetry 15 00065 g001
Table 1. Absolute errors of the RKHSM solution for Example 1.
Table 1. Absolute errors of the RKHSM solution for Example 1.
RKHSM-Absolute Error
(x,t) α = 0 . 9 α = 0 . 8 α = 0 . 75
(0.1, 0.1) 1.3016 × 10 3 1.3446 × 10 3 1.9199 × 10 3
(0.3, 0.3) 1.0876 × 10 4 3.8615 × 10 4 8.1880 × 10 4
(0.5, 0.5) 5.9252 × 10 5 8.7214 × 10 4 1.3792 × 10 3
(0.7, 0.7) 5.9929 × 10 4 1.0365 × 10 3 1.6629 × 10 3
(0.9,0.9) 9.5379 × 10 4 1.6982 × 10 4 4.2782 × 10 4
Table 2. Comparison between the exact solution and FITM solution for Example 2.
Table 2. Comparison between the exact solution and FITM solution for Example 2.
(x, t)ApproximateExactAbsolute Error
(0.1, 0.1)0.02150.0213 1.5308 × 10 4
(0.2, 0.2)0.08930.0864 2.9000 × 10 3
(0.3, 0.3)0.20160.1960 5.6000 × 10 3
(0.4, 0.4)0.31330.3076 5.7000 × 10 3
(0.5, 0.5)0.49580.4958 3.9629 × 10 8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Partohaghighi, M.; Akgül, A.; Akgül, E.K.; Attia, N.; De la Sen, M.; Bayram, M. Analysis of the Fractional Differential Equations Using Two Different Methods. Symmetry 2023, 15, 65. https://doi.org/10.3390/sym15010065

AMA Style

Partohaghighi M, Akgül A, Akgül EK, Attia N, De la Sen M, Bayram M. Analysis of the Fractional Differential Equations Using Two Different Methods. Symmetry. 2023; 15(1):65. https://doi.org/10.3390/sym15010065

Chicago/Turabian Style

Partohaghighi, Mohammad, Ali Akgül, Esra Karatas Akgül, Nourhane Attia, Manuel De la Sen, and Mustafa Bayram. 2023. "Analysis of the Fractional Differential Equations Using Two Different Methods" Symmetry 15, no. 1: 65. https://doi.org/10.3390/sym15010065

APA Style

Partohaghighi, M., Akgül, A., Akgül, E. K., Attia, N., De la Sen, M., & Bayram, M. (2023). Analysis of the Fractional Differential Equations Using Two Different Methods. Symmetry, 15(1), 65. https://doi.org/10.3390/sym15010065

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