Next Article in Journal
Adsorption on Fractal Surfaces: A Non Linear Modeling Approach of a Fractional Behavior
Next Article in Special Issue
Multi-Model Selection and Analysis for COVID-19
Previous Article in Journal
An Inverse Source Problem for the Generalized Subdiffusion Equation with Nonclassical Boundary Conditions
Previous Article in Special Issue
Analytic Solution of the Langevin Differential Equations Dominated by a Multibrot Fractal Set
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator

1
Department of Mathematical Problems of Control and Cybernetics, Yuriy Fedkovych Chernivtsi National University, Chernivtsi Region, 58012 Chernivtsi, Ukraine
2
Department of Mathematical Modeling, Yuriy Fedkovych Chernivtsi National University, Chernivtsi Region, 58012 Chernivtsi, Ukraine
3
Department of Mathematics, Rajdhani College, University of Delhi, Delhi 110015, India
4
Department of Mathematics, AMITY School of Applied Science, AMITY University, Jaipur 302002, India
*
Author to whom correspondence should be addressed.
Submission received: 24 May 2021 / Revised: 14 June 2021 / Accepted: 21 June 2021 / Published: 2 July 2021

Abstract

:
In this manuscript, the time-fractional diffusion equation in the framework of the Yang–Abdel–Cattani derivative operator is taken into account. A detailed proof for the existence, as well as the uniqueness of the solution of the time-fractional diffusion equation, in the sense of YAC derivative operator, is explained, and, using the method of α -HATM, we find the analytical solution of the time-fractional diffusion equation. Three cases are considered to exhibit the convergence and fidelity of the aforementioned α -HATM. The analytical solutions obtained for the diffusion equation using the Yang–Abdel–Cattani derivative operator are compared with the analytical solutions obtained using the Riemann–Liouville (RL) derivative operator for the fractional order γ = 0.99 (nearby 1) and with the exact solution at different values of t to verify the efficiency of the YAC derivative operator.
MSC:
Primary 92B05, 92C60; Secondary 26A33

1. Introduction

A parabolic partial differential equation that delineates the movement of energy and matter in a medium is called a diffusion equation. The diffusion process of heat or mass can be described using the concept of ordinary and partial derivatives. In today’s era, fractional calculus [1,2,3,4,5,6] is emerging as an efficient and powerful tool in the field of science and technology. It is the branch of mathematics pertaining to the derivatives and integrals of arbitrary order and it is fruitful in explaining the concepts of damping, wave propagation and diffusion, biology, genetic algorithms, control systems, economy and finance, signal processing, robotics, system identification, electromagnetism, heat transfer, and many more. The literature is brimming with developments made in the field of fractional calculus. The most widely accepted definition including the singular kernel was proposed by Riemann and Liouville–Caputo. The next classification of fractional derivatives are made on the basis of the non-singular kernel. Caputo–Fabrizio proposed the first definition in the field of fractional calculus with a non-singular kernel [7]. Further development and success in the definitions of fractional derivatives were introduced by Prabhakar, Sonine, Wiman, Miller–Ross, Gorenglo, Mainardi, Mittag–Leffler, Atangana–Baleanu, Yang–Abdel–Cattani, and a lot more, who propose that the non-singular kernels are special functions, such as the Mittag–Leffler function, Miller–Ross function, Wiman function, Kohlrausch–William–Watts function, Rabotnov function, Prabhakar function, etc., see [8,9,10,11,12,13]. The fractional derivatives in which the non-singular kernels are the special functions are called general fractional-order derivatives. Further applications of fractional derivatives can be seen in [14,15,16,17,18,19,20,21,22,23,24,25,26,27]. In this work, we will find the analytical solution of the following nonlinear fractional diffusion equation in the framework of the Yang–Abdel–Cattani (YAC) derivative operator.
0 Y A C D t γ ( w ( η , t ) ) = η w a w ( η , t ) η w ( η , 0 ) = w 0 ,
where 0 Y A C D t γ ( w ( η , t ) ) represents the Yang–Abdel–Cattani (YAC) fractional derivative of w ( η , t ) , w is the density of the diffusing medium at point η and at time t. This generalized fractional derivative was presented by Yang et al. with the Rabotnov exponential function as the non-singular kernel. In this work, we will present a detailed proof for the existence as well as the uniqueness of the solution of the time-fractional diffusion equation in the framework of the YAC derivative operator and, using the α -homotopy analysis transform method, we will find the analytical solution of the fractional diffusion equation in the sense of the YAC derivative operator. The analytical solutions obtained for the fractional diffusion equation using the YAC derivative operator are compared to the analytical solution obtained using the Riemann–Liouville (RL) derivative operator for the fractional order γ = 0.99 (nearby 1) and with the exact solution at different values of t to verify the efficiency of the YAC fractional derivative operator. Graphical representations of the analytical solutions are also given for a better understanding of the α -HATM.

2. Definitions

Definition 1.
[13] We define γ , χ R + ; the following series defines the Rabotnov exponential function of order γ
Ψ ( χ u γ ) = s = 0 χ s u ( s + 1 ) ( γ + 1 ) 1 Γ ( ( s + 1 ) ( γ + 1 ) ) , u C
Definition 2.
[13] For k on L 1 [ a , b ] , t > 0 , χ R + 0 < γ 1 , the following defines the Yang–Abdel–Cattani fractional derivative of order γ
0 Y A C D t γ ( k ( t ) ) = 0 t Ψ γ ( χ ( t ρ ) γ ) k ( ρ ) d ρ .
where, Ψ γ represents the Rabotnov exponential function of order γ.
Definition 3.
[13] The following defines the Laplace transform for the Yang–Abdel–Cattani fractional derivative
L ( 0 Y A C D t γ ( k ( t ) ) ) = 1 p γ + 1 p L [ k ( t ) ] k ( 0 ) 1 + χ p ( γ + 1 )
Definition 4.
[13] For k∈ L 1 [ a , b ] , t > 0 , 0 < γ 1 , χ R + , the following defines the fractional integral with Rabotnov fractional exponential function, of order γ
I Y A C γ k ( t ) = 0 t Ψ γ ( χ ( t ρ ) γ ) k ( ρ ) d ρ .
Definition 5.
[13] The following defines the Laplace transform for the Yang–Abdel–Cattani fractional integral
L ( I Y A C γ k ( t ) ) = 1 p γ + 1 L ( k ( t ) ) 1 + χ p ( γ + 1 )

3. Existence of Solution of Fractional Diffusion Equation Using Yang–Abdel–Cattani Derivative Operator

Theorem 1.
Let us assume that the function f ( η , t , w , w , w ) satisfies the Lipschitz condition as
| f ( η , t , w , w , w ) f ( η , t , w 1 , w 1 , w 1 ) | M 1 | w w 1 | + M 2 | w w 1 | + M 3 | w w 1 | .
We also assume that
| w w 1 | k 1 | w w 1 |
| w w 1 | k 2 | w w 1 |
where k 1 , k 2 R + then there exists a unique solution for the following time-fractional differential equation.
0 Y A C D t γ ( w ( η , t ) ) = η w a w ( η , t ) η
Proof. 
We define
Φ ( w , η ) = f ( η , t , w , w , w ) = η w a w ( η , t ) η
We first show that Φ ( w , η ) satisfies Lipschitz condition. Consider
| | Φ ( w , η ) Φ ( w 1 , η ) | | = | | f ( η , t , w , w , w ) f ( η , t , w 1 , w 1 , w 1 ) | | M 1 | | w w 1 | | + M 2 | | w w 1 | | + M 3 | | w w 1 | | M 1 | | w w 1 | | + M 2 k 1 | | w w 1 | | + M 3 k 2 | | w w 1 | | = ( M 1 + M 2 k 1 + M 3 k 2 ) | | w w 1 | | .
We also define M 1 + M 2 k 1 + M 3 k 2 = M . So finally, we have
| | Φ ( w , η ) Φ ( w 1 , η ) | | = | | f ( η , t , w , w , w ) f ( η , t , w 1 , w 1 , w 1 ) | |
                                                                                           M | w w 1 |
Using Picard’s theorem, we obtain
w ( η , t ) = w ( η , 0 ) + 0 t Ψ γ ( χ ( t ρ ) γ ) Φ ( w , η ( ρ ) ) d ρ .
For convenience, we write
0 t Ψ γ ( χ ( t ρ ) γ ) Φ ( w , η ( ρ ) ) d ρ = I Y A C γ Φ ( w , η )
Finally, we have
(14) w ( η , t ) = w ( η , 0 ) + I Y A C γ Φ ( w , η ) (15) w ( η , t ) w ( η , 0 ) = I Y A C γ Φ ( w , η ) (16) | | w ( η , t ) w ( η , 0 ) | | = | | I Y A C γ Φ ( w , η ) | | (17) = 0 t Ψ γ ( χ ( t ρ ) γ ) Φ ( w , η ( ρ ) ) d ρ . (18) 0 t | | Ψ γ ( χ ( t ρ ) γ ) | | | | Φ ( w , η ( ρ ) ) | | d ρ . (19) | | Φ ( w , η ( ρ ) ) | | 0 t Ψ γ ( χ ( t ρ ) γ ) (20) = I Y A C γ ( 1 ) | | Φ ( w , η ( ρ ) ) | |
As we proved that Φ ( w , η ) satisfies Lipschitz condition, so the following holds
| | Φ ( w , η ) | | K
Hence
| | w ( η , t ) w ( η , 0 ) | | K I Y A C γ ( 1 )
Finally, we consider
(23) | | w ( η , t ) w 1 ( η , t ) | | = | | I Y A C γ Φ ( w , η ) I Y A C γ Φ ( w 1 , η ) | | (24) I Y A C γ ( 1 ) | | Φ ( w , η ) Φ ( w 1 , η ) | | (25) I Y A C γ M | | w w 1 | |
For the above map to be a contraction, we must have
(26) M I Y A C γ 1 (27) I Y A C γ 1 M
Hence the existence and the uniqueness of the solution follows as a consequence of the Banach fixed point theorem. □

4. α -HATM Solution of Nonlinear Time-Fractional Diffusion Equation

Consider the nonlinear fractional diffusion equation, given as:
0 Y A C D t γ ( w ( η , t ) ) = η w a w ( η , t ) η
We rewrite the above equation as
0 Y A C D t γ ( w ( η , t ) ) a w a 1 w η 2 + w a 2 w η 2 = 0
Taking the Laplace transform on the two sides of Equation (29),
1 p γ + 1 p L [ w ( η ) ] w ( η , 0 ) 1 + χ p ( γ + 1 ) L a w a 1 w η 2 + w a 2 w η 2 = 0
By simplifying, we obtain
L [ w ( η , t ) ] w ( η , 0 ) p p γ + 1 + χ p L a w a 1 w η 2 + w a 2 w η 2 = 0
Let N be a nonlinear operator, defined as
N [ β ( η , t , α ) ] = L [ β ( η , t , α ) ] w ( η , 0 ) p p γ + 1 + χ p L a β a 1 β η 2 + β a 2 β η 2
where β ( η , t , α ) is a function in η , t , α , and α [ 0 , 1 / r ] is an embedding parameter. Now construct the homotopy as
( 1 r α ) L [ β ( η , t , α ) w 0 ( η , t ) ] = α k N [ β ( η , t , α ) ]
where β ( η , t , α ) is a function of η , t and α , w 0 ( η , t ) is an initial guess of w ( η , t ) , L is the Laplace transform, and k 0 is an auxiliary parameter. From above equation, we see that
when α = 0 , β ( η , t , 0 ) = w 0 ( η , t )
when α = 1 r , β ( η , t , 1 / r ) = w ( η , t )
This shows that as α varies from 0 to 1/r, the solution β ( η , t , α ) changes from w 0 ( η , t ) the initial guess to w ( η , t ) the exact solution.
Expand β ( η , t , α ) with respect to α using the Taylor series, we obtain
β ( η , t , α ) = w 0 η , t + i = 1 w i η , t α i
where
w i η , t = 1 i ! i β ( η , t ; α ) α i α = 0
Let w 0 η , t , k , r be selected appropriately, the series defined in Equation (35) converges at α = 1 / r , hence
w ( η , t ) = w 0 η , t + i = 1 w i η , t 1 r i
Defining the vectors w ¯ n = w 0 , w 1 , , w n , and differentiating Equation (33) i-times with respect to α , and substituting α = 0 , and lastly dividing them by i ! , we obtain:
L w i η , t ζ i w i 1 η , t = k i ( w i 1 η , t ) ,
where
i ( w i 1 ( η , t ) ) = 1 ( i 1 ) ! i 1 N ( η , t ; α ) α i 1 α = 0
and
ζ i = 0 , i 1 r , o t h e r w i s e .
Using Equations (32) and (38), we obtain
i ( w i 1 ( η , t ) ) = L [ w i 1 ] w ( η , 0 ) p 1 ζ i r + L a w i 1 a 1 w i 1 η 2 + w i 1 a 2 w i 1 η 2
Lastly, take the inverse Laplace transform on both sides of Equation (37),
w m η , t = ζ i w i 1 η , t + k L 1 i ( w i 1 ) .
Finally, opting for the suitable values of k and r, the α -HATM series solution is obtained, which is given as
w ( η , t ) = lim N i = 0 w i ( η , t ) 1 r i
We now consider different cases of the above nonlinear diffusion equation.

4.1. Case 1

We will now find the analytical solution of the following non-linear fractional diffusion equation in sense of YAC derivative operator using above mentioned α HATM.
0 Y A C D t γ ( v ( η , t ) ) = η v ( η , t ) η , 0 < t < 1 ; 0 < γ < 1 , v ( η , 0 ) = s i n ( π η )
Using the α -HATM, the series solution is given as
v 0 Y A C ( η , t ) = s i n ( π η )
v 1 Y A C ( η , t ) = k π 2 t χ + t γ Γ ( 1 γ ) s i n ( π η )
v 2 Y A C ( η , t ) = r v 1 Y A C ( η , t ) + 1 2 k 2 π 2 2 π 2 t 2 γ Γ ( 1 2 γ ) + 2 t γ Γ ( 1 γ ) t χ 2 + π 2 t χ + 4 π 2 t χ Γ ( 2 γ ) s i n ( π η )
The following diffusion equation’s
v ( η , t ) t = 2 v ( η , t ) η 2
exact solution is given as
s i n ( π η ) e π 2 t
In Table 1, we will compare the α -HATM solutions obtained for Case 1 in terms of the YAC fractional derivative operator with the analytical solutions obtained using the Riemann–Liouville fractional derivative operator [28] and with the considered diffusion equation’s exact solution at different values of t for η = 0.25 and γ = 0.99 ( γ close to 1) k = 0.01 , r = 3 and χ = 3 .
In Figure 1, the 3-D plots of the α -HATM solution for the YAC operator are compared to the exact solution for γ = 0.99 ( γ nearby 1).

4.2. Case 2

We will now find the analytical solution of the following non-linear fractional diffusion equation in sense of YAC derivative operator using above mentioned α HATM.
0 Y A C D t γ ( v ( η , t ) ) = η v 2 v ( η , t ) η , 0 < γ < 1 ; 0 < t < 1 , v ( η , 0 ) = η + b c
Using the α -HATM, the series solution is given as
v 0 Y A C ( η , t ) = η + b c
v 1 Y A C ( η , t ) = k ( b + η ) 4 c 3 t χ + t γ Γ ( 1 γ )
v 2 Y A C ( η , t ) = r v 1 Y A C ( η , t ) + k 2 t 2 γ ( b + η ) 4 c 5 Γ ( 1 2 γ ) Γ ( 1 γ ) Γ ( 2 γ ) 2 Γ ( 1 γ ) Γ ( 2 γ ) + t γ Γ ( 1 2 γ ) ( c 2 Γ ( 1 2 γ ) + t χ Γ ( 1 γ ) [ 4 t γ ( c 2 t χ ) Γ ( 2 γ ) ] )
The following diffusion equation’s
v ( η , t ) t = η v 2 v ( η , t ) η
exact solution is given as
η + b 2 c 2 t , t c 2
In Table 2, we will compare the α -HATM solutions obtained for Case 2 in terms of the YAC fractional derivative operator with the analytical solutions obtained in terms of the Riemann–Liouville fractional derivative operator [28] and with the considered diffusion equation’s exact solution at different values of t for k = 0.19 , r = 1 , χ = 1 , b = 2 , c = 3 , η = 0.25 and γ = 0.99 .
In Figure 2, the 3-D plots of the α -HATM solution for the YAC operator are compared to the exact solution for γ = 0.99 ( γ nearby 1).

4.3. Case 3

We will now find the analytical solution of the following non-linear fractional diffusion equation in sense of YAC derivative operator using above mentioned α HATM.
0 Y A C D t γ ( v ( η , t ) ) = η v 2 v ( η , t ) η , 0 < t < 1 ; 0 < γ < 1 , v ( η , 0 ) = 1 1 + η 2
Using the α -HATM, the series solution is given as
v 0 Y A C ( η , t ) = 1 1 + η 2
v 1 Y A C ( η , t ) = k η 2 ( 1 + η 2 ) 3 / 2 1 1 + η 2
The following diffusion equation’s
v ( η , t ) t = η v 2 v ( η , t ) η
exact solution is
1 η 2 + e 2 t
In Table 3, we will compare the α -HATM solutions obtained for Case 3 in terms of the YAC fractional derivative operator with the analytical solutions obtained in terms of the Riemann–Liouville fractional derivative operator [28] and with the considered diffusion equation’s exact solution at different values of t for k = 0.6 , r = 1 , χ = 1 , η = 0.25 , and γ = 0.99 .
In Figure 3, the 3-D plots of the α -HATM solution for the YAC operator are compared to the exact solution for γ = 0.99 ( γ nearby 1).

5. Conclusions

In this paper, the time-fractional nonlinear diffusion equation is taken into consideration in regards to the Yang–Abdel–Cattani fractional derivative operator. The α -HATM is used to find the analytical solution of the nonlinear fractional diffusion equation. The analytical solutions obtained from the α -HATM, in the sense of the YAC derivative operator, are compared to the analytical solutions obtained via the Riemann–Liouville derivative operator and with the exact solution for all the three cases at distinct values of time t, and we observe that the analytical solutions obtained using the YAC derivative operator coincide with the exact solution more closely as compared to the Riemann–Liouville derivative operator when the value of fractional order is close to 1, i.e., γ = 0.99 . In Figure 1, Figure 2, Figure 3 we plotted the 3-D representations of the solution v ( η , t ) with the exact solution for case 1, case 2, and case 3, respectively. In Figure 1a, Figure 2a, Figure 3a, the 3-D representation of the analytical solution attained using α -HATM along with the exact solution is given. In Figure 1b, Figure 2b, Figure 3b, the α -HATM solutions are plotted for different values of k along with the exact solution and we see that the α -HATM solution is in best fit with the exact solution for a suitable value of k and hence k = 0.01 works as an optimal value for case 1, k = 0.19 works as an optimal value for case 2, and k = 0.6 works as an optimal value for case 3. In Figure 1c, Figure 2c, Figure 3c, the α -HATM solutions are plotted for different values of r along with the exact solution and we see that the α -HATM solution is in best fit with the exact solution for a suitable value of r and hence r = 3 works as an optimal value for case 1 and r = 1 works as an optimal value for cases 2 and 3. In Figure 1d, Figure 2d, Figure 3d, the α -HATM solutions are plotted for different values of χ along with the exact solution and we see that the α -HATM solution is in best fit with the exact solution for a suitable value of χ and hence χ = 3 works as an optimal value for case 1 and χ = 1 works as optimal value for case 2 and 3. Hence, we see the efficient role of the parameters k, r provided in the α -HATM, and χ , given in the YAC definition in regulating the convergence of the solution. Hence we see that YAC fractional derivative operator provides better results than the Riemann–Liouville fractional derivative operator.

Author Contributions

Conceptualization, I.V.M., M.G., S.S. and R.S.D.; methodology, S.S. and R.S.D.; software, S.S. and R.S.D.; validation, I.V.M., M.G., A.C. and R.S.D.; formal analysis, I.V.M., M.G., A.C. and R.S.D.; writing-original draft preparation, S.S. and R.S.D.; writing–review and editing, S.S., R.S.D. and A.C.; supervision, I.V.M., M.G., A.C. and R.S.D.; funding acquisition, I.V.M. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Igor V. Malyk.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our sincere thanks to the reviewers for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; A Wiley-Interscience Publication; John Wiley and Sons, Inc.: New York, NY, USA, 1993; p. xvi+366. ISBN 0-471-58884-9. [Google Scholar]
  2. Oldham, K.B.; Spanier, J. The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order; With an Annotated Chronological Bibliography by Bertram Ross; Academic Press: New York, NY, USA; London, UK,, 1974; Volume 111, p. xiii+234. ISBN 9780125255509. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006; Volume 204, p. xvi+523. ISBN 978-0-444-51832-3. [Google Scholar]
  4. 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, Inc.: San Diego, CA, USA, 1999; Volume 198, p. xxiv+340. ISBN 0-12-558840-2. [Google Scholar]
  5. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2014; p. x+293. ISBN 978-981-4579-89-6. [Google Scholar] [CrossRef]
  6. Yang, X.J. Advanced Local Fractional Calculus and Its Applications; World Science: New York, NY, USA, 2012; ISBN 13 9781938576010. [Google Scholar]
  7. Caputo, M.; Fabrizio, M. A New Definition of Fractional Derivative without Singular Kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar] [CrossRef]
  8. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and nonsingular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  9. Khalil, R.; Horani, M.A.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. [Google Scholar] [CrossRef]
  10. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst.-S 2020, 13, 709–722. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, X.J. General Fractional Derivatives: Theory, Methods and Applications; Chapman and Hall/ CRC: Boca Raton, FL, USA, 2019; ISBN 9781138336162. [Google Scholar]
  12. Yang, X.J.; Gao, F.; Ju, Y. General Fractional Derivatives with Applications in Viscoelasticity; Academic Press: Cambridge, MA, USA, 2020; ISBN 9780128172087. [Google Scholar]
  13. Yang, X.J.; Abdel-Aty, M.; Cattani, C. A new general fractional order derivative with Rabotnov fractional-exponential kernel applied to model the anomalous heat. Therm. Sci. 2019, 23, 1677–1681. [Google Scholar] [CrossRef] [Green Version]
  14. Atangana, A. On the new fractional derivative and application to non-linear Fisher’s reaction-diffusion equation. Appl. Math. Comput. 2016, 273, 948–956. [Google Scholar] [CrossRef]
  15. Aguilar, J.P.; Korbel, J.; Luchko, Y. Application of the Fractional Diffusion Equation to Option Pricing and Risk Calculations. Mathematics 2019, 7, 796. [Google Scholar] [CrossRef] [Green Version]
  16. Sontakke, B.R.; Shelke, A.S. Appropriate Scheme for Time fractional Diffusion Equation and Its Applications. Glob. J. Pure Appl. Math. 2019, 13, 4333–4345. [Google Scholar]
  17. Constantinescu, C.D.; Ramirez, J.M.; Zhu, W.B. An application of fractional differential equations to risk theory. Math. Stochastics 2019, 23, 1001–1024. [Google Scholar] [CrossRef] [Green Version]
  18. Ali, K.K.; Osman, M.S.; Baskonus, H.C.; Elazabb, N.S.; İlhan, E. Analytical and numerical study of the HIV-1 infection of CD4+ T-cells conformable fractional mathematical model that causes acquired immunodeficiency syndrome with the effect of antiviral drug therapy. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  19. Gao, W.; Ghanbari, B.; Baskonus, H.C. New numerical simulations for some real world problems with Atangana-Baleanu fractional derivative. Chaos Solitons Fractals 2019, 128, 34–43. [Google Scholar] [CrossRef]
  20. Ghanbari, B.; Gómez-Aguilar, J.F. Modelling the dynamics of nutrient-phytoplankton-zooplankton system with variable-order fractional derivatives. Chaos Solitons Fractals 2018, 116, 114–120. [Google Scholar] [CrossRef]
  21. Kumar, S.; Ghosh, S.; Samet, B.; Goufo, E.F.D. An analysis for heat equations arises in diffusion process using new Yang-Abdel-Aty-Cattani fractional operator. Math. Methods Appl. Sci. 2020, 1–19. [Google Scholar] [CrossRef]
  22. Jleli, M.; Kumar, S.; Kumar, R.; Samet, B. Analytical approach for time fractional wave equations in the sense of Yang-Abdel-Aty-Cattani via the homotopy perturbation transform method. Alex. Eng. J. 2019. [Google Scholar] [CrossRef]
  23. Zhang, X.; Chen, Y.Q. Admissibility and robust stabilization of continuous linear singular fractional order systems with the fractional order α: The 0 < α < 1 case. ISA Trans. 2018, 82, 42–50. [Google Scholar] [CrossRef] [PubMed]
  24. Belmor, S.; Ravichandran, C.; Jarad, F. Nonlinear generalized fractional differential equations with generalized fractional integral conditions. J. Taibah Univ. Sci. 2020. [Google Scholar] [CrossRef] [Green Version]
  25. Tarasov, V.E. Fractional nonlinear dynamics of learning with memory. Non-Linear Dyn. 2020, 100, 1231–1242. [Google Scholar] [CrossRef]
  26. Shone, T.T.; Patra, A. Solution for non-linear fractional partial differential equations using fractional complex transform. Int. J. Appl. Comput. Math. 2019, 5, 90. [Google Scholar] [CrossRef]
  27. Tarasov, V.E. Exact solutions of Bernoulli and Logistic fractional differential equations with power law coefficients. Mathematics 2020, 8, 2231. [Google Scholar] [CrossRef]
  28. Dubey, R.S.; Goswami, P. Analytical solution of the nonlinear diffusion equation. Eur. Phys. J. Plus 2018, 133. [Google Scholar] [CrossRef]
Figure 1. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 1 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with exact solution for distinct values of k, Yellow (k = −0.01), Blue (k = −0.04), Green (k = −0.08), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 3), Blue ( χ = 5), Green ( χ = 7), Red(Exact Solution).
Figure 1. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 1 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with exact solution for distinct values of k, Yellow (k = −0.01), Blue (k = −0.04), Green (k = −0.08), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 3), Blue ( χ = 5), Green ( χ = 7), Red(Exact Solution).
Fractalfract 05 00064 g001
Figure 2. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 2 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with the exact solution for distinct values of k, Yellow (k = −0.19), Blue (k = −1), Green (k = −2), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 1), Blue ( χ = 5), Green ( χ = 10), Red (Exact Solution).
Figure 2. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 2 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with the exact solution for distinct values of k, Yellow (k = −0.19), Blue (k = −1), Green (k = −2), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 1), Blue ( χ = 5), Green ( χ = 10), Red (Exact Solution).
Fractalfract 05 00064 g002
Figure 3. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 3 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with exact solution for distinct values of k, Yellow (k = −0.6), Blue (k = −0.8), Green (k = −1), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 1), Blue ( χ = 2), Green ( χ = 3), Red (Exact Solution).
Figure 3. 3-D Plot of α -HATM solution for YAC operator with exact solution for case 3 when γ = 0.99 . (a) α -HATM solution with exact solution. (b) α -HATM solution with exact solution for distinct values of k, Yellow (k = −0.6), Blue (k = −0.8), Green (k = −1), Red (Exact Solution). (c) α -HATM solution with exact solution for distinct values of r, Yellow (r = 1), Blue (r = 2), Green (r = 3), Red (Exact Solution). (d) α -HATM solution with exact solution for distinct values of χ , Yellow ( χ = 1), Blue ( χ = 2), Green ( χ = 3), Red (Exact Solution).
Fractalfract 05 00064 g003aFractalfract 05 00064 g003b
Table 1. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 1.
Table 1. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 1.
t v ( RL ) v ( YAC ) v ( Exact ) Error ( RL ) Error ( YAC )
0.10.1965710.2634640.263544 6.6973 × 10 2 8 × 10 5
0.20.9020910.09816880.098225 8.0159 × 10 3 5.62 × 10 5
0.30.03030860.03613510.0366092 6.3006 × 10 3 4.741 × 10 4
0.40.01653240.01317740.0136445 2.8879 × 10 3 4.671 × 10 4
0.50.004639990.005079240.00508543 4.4544 × 10 4 6.19 × 10 6
0.60.001341190.001882480.00189538 5.5419 × 10 4 1.29 × 10 5
0.70.0005917540.0006453950.000706423 1.14669 × 10 4 6.1028 × 10 5
0.80.0002086320.0002403790.000263289 5.4657 × 10 5 2.291 × 10 5
0.90.00009355880.00009716670.00009813 4.5712 × 10 6 9.633 × 10 7
1.00.00003009470.00003515930.0000365738 6.4791 × 10 6 1.4145 × 10 6
Table 2. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 2.
Table 2. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 2.
t v ( RL ) v ( YAC ) v ( Exact ) Error ( RL ) Error ( YAC )
0.10.3629160.3770390.377101 1.4185 × 10 2 6.2 × 10 5
0.20.3754160.3777590.379237 3.821 × 10 3 1.478 × 10 3
0.30.3794740.3810230.381411 1.937 × 10 3 3.88 × 10 4
0.40.3801590.3830750.383622 3.463 × 10 3 5.47 × 10 4
0.50.3824430.3842920.3858872 3.4442 × 10 3 1.58 × 10 3
0.60.3875660.3880410.388162 5.96 × 10 4 1.21 × 10 4
0.70.3887690.3904070.390493 1.724 × 10 3 8.6 × 10 5
0.80.390340.3926720.392867 2.527 × 10 3 1.72 × 10 4
0.90.3907570.3950340.395285 4.528 × 10 3 2.51 × 10 4
1.00.4054650.3977220.397748 7.717 × 10 3 2.6 × 10 5
Table 3. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 3.
Table 3. Comaprison of α HATM Solution in sense of YAC derivative operator, RL derivative operator and exact solution for Case 3.
t v ( RL ) v ( YAC ) v ( Exact ) Error ( RL ) Error ( YAC )
0.10.8425890.8615150.882539 3.995 × 10 2 2.1024 × 10 2
0.20.7804120.8053780.802101 2.1689 × 10 2 3.277 × 10 3
0.30.6984540.7233380.728431 2.9972 × 10 2 5.093 × 10 3
0.40.6400230.641230.661101 2.1078 × 10 2 1.9871 × 10 2
0.50.5898740.5590940.599676 9.802 × 10 3 4.058 × 10 3
0.60.5211760.5536650.543718 2.254 × 10 2 9.947 × 10 3
0.70.4623970.4842870.492802 3.0405 × 10 2 8.515 × 10 3
0.80.4600850.458750.446521 1.3564 × 10 2 1.23 × 10 2
0.90.4287530.3948340.404486 2.4267 × 10 2 9.652 × 10 3
1.00.3875620.3309260.366333 2.1229 × 10 2 3.24 × 10 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Malyk, I.V.; Gorbatenko, M.; Chaudhary, A.; Sharma, S.; Dubey, R.S. Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator. Fractal Fract. 2021, 5, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030064

AMA Style

Malyk IV, Gorbatenko M, Chaudhary A, Sharma S, Dubey RS. Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator. Fractal and Fractional. 2021; 5(3):64. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030064

Chicago/Turabian Style

Malyk, Igor V., Mykola Gorbatenko, Arun Chaudhary, Shivani Sharma, and Ravi Shanker Dubey. 2021. "Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator" Fractal and Fractional 5, no. 3: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030064

Article Metrics

Back to TopTop