Abstract
The objective of this article is to present the computable solution of space-time advection-dispersion equation of fractional order associated with Hilfer-Prabhakar fractional derivative operator as well as fractional Laplace operator. The method followed in deriving the solution is that of joint Sumudu and Fourier transforms. The solution is derived in compact and graceful forms in terms of the generalized Mittag-Leffler function, which is suitable for numerical computation. Some illustration and special cases of main theorem are also discussed.
Introduction
In the last decade, considerable interest in fractional differential equations has been stimulated due to their numerous applications in the areas of physics, biology, engineering, and other areas. Several numerical and analytical methods have been developed to study the solutions of nonlinear fractional partial differential equations, for details, refer to the work in [1–6]. Fractional equations have enabled the investigation of the nonlocal response of multiple phenomena such as diffusion processes, electrodynamics, fluid flow, elasticity, and many more. Nowadays, fractional derivatives have gained a significant development to model some real life phenomena in the form of partial differential equations or the ordinary equations. Several researchers have performed the numerical simulation for fractional problems and revealed their applications in different directions include [7–12] and references therein. The exchange of heat, mass and momentum are considered to be the fundamental transfer phenomena in the universe. The mathematical framework for heat and mass transfer are of same kind, basically encompass by advection-dispersion equation. In recent work many authors have demonstrated the depth of mathematics and related physical issues of advection-dispersion equations. Schumer et al. [13] gave physical interpretation of space-time fractional advection-dispersion equation. Space-time fractional advection-dispersion equations are generalizations of classical advection-dispersion equations. The use of Hilfer-Prabhakar fractional derivative operator is gaining importance in physics because of their specific properties. The objective of this paper is to derive the solution of Cauchy type generalized fractional advection dispersion equation (18), associated with the Hilfer-Prabhakar fractional derivative. This paper provides an elegant extension of results, given earlier by Haung and Liu [14], Haubold et al. [15], Saxena et al. [16], and Agarwal et al. [17].
Results Required in the Sequel
In early 90s, Watugala [18] introduced Sumudu transform, which is defined as,
for all real t ≥ 0 the Sumudu transform of function f(t) ∈ A is defined as,
inversion formula of (2), is given by
where γ being a fixed real number.
Among others, the Sumudu transform was shown to have units preserving properties, and hence may be used to solve problems without resorting to the frequency domain. Further details and properties about this transform can be found in Belgacem [19], Belgacem et al. [20], and Katatbeh and Belgacem [21].
For a function u(x, t), the Fourier transform of with respect to x is defined by
and for the function u*(η, t), inverse Fourier transform with respect to η is given by the formula
For more details of Fourier transform, see [Debnath and Bhatta [22]].
Mittag-Leffler function of two parameters is studied by Wiman [23] as
Mittag-Leffler function of three parameter introduced by Prabhakar [24] as
Riemann-Liouville fractional integral (right-sided) of order α is defined in [25]
The right sided Riemann-Liouville fractional derivative of order α defined as
here [x] is the integral part of x.
Caputo [26], introduced fractional derivative of order R(α) > 0 as
The Sumudu transform of (10) is given in [27], as
where ū(x, s) is the Sumudu transform of u(x, t).
Hilfer [28], gave a fractional derivative operator of two parameters μ and ν, which is generalization of (9) and (10), in the form
For ν = 0, equation (12) reduces into (9) and for ν = 1, equation (12) reduces into (10).
The Sumudu transform of (12) is given in [29], as
Where the initial value term involves the Riemann-Liouville fractional derivative operator of order (1 − ν)(1 − μ) as t → 0 +.
A generalization of Hilfer derivate is given in [30], known as Hilfer-Prabhakar derivative, is defined as:
Let μ ∈ (0, 1), ν ∈ [0, 1], and let f. belongs to the set of locally integrable real valued functions i.e., The Hilfer-Prabhakar derivative is defined by
where γ, ω ∈ R, ρ > 0, and where We observe that (14) reduces to the Hilfer derivative for γ = 0. The Sumudu transform of this derivative operator (14) is given in [31], in the form:
For details of this derivative, refer to the work in [30, 31].
Brockmann and Sokolov [32], defined a fractional Laplace operator as:
where the operators are defined by
and
The Fourier transform of is given in [32], as
Inverse Sumudu transform of the following function is directly applicable in this sequel:
In the complex plane C, for any R(α) > 0, R(β) > 0, and ω ∈ C
Space-time Fractional Advection-Dispersion Equation
Here we will find, the solution of the generalized space-time Advection-Dispersion equation (18) under the conditions given in (19) and (20). Our main findings in the form of the following Theorem 3.1 and Corollary 3.2.
Theorem 3.1. Consider the generalized fractional order space-time advection-dispersion equation of Cauchy type
where λ ∈ (0, 2] x ∈ R, t ∈ R+, μ ∈ (0, 1), ν ∈ [0, 1],
with initial condition,
and boundary condition
where is the Laplace operator of fractional order λ, λ ∈ (0, 2]. The positive constant η represent the average fluid velocity and ς (positive constant) represent the dispersion coefficient. Subject to the above constraints, solution of equation (18), is
Proof: First, take the Fourier transform of equation (18) with respect to the space variable x, then
u*(k, t) represent Fourier transform of u(x, t). Again, apply Sumudu transform on (22) with respect to time variable t, we get
where S[u(k, t);s] = ū(k, s).
Solve equation (23), by using conditions (19)-(20), we get
On taking inverse Sumudu transform of equation (24), and after little simplification, apply result (17), it gives
Taking inverse Fourier transform of (25), get our required result (21).
This completes the proof of the theorem 3.1.
On taking in Theorem 3.1, we arrive at:
Corollary 3.2. Consider the following one dimensional space-time Schrödinger equation of fractional order, for a free nature particle of mass m is
where λ ∈ (0, 2], x ∈ R, t ∈ R+, μ ∈ (0, 1), ν ∈ [0, 1],
with initial condition
and boundary condition
where is same as we defined earlier and h = 6.625 × 10−27ergs = 4.21 × 10−21Mev s is the Planck constant. Subject to the above constraints, solution of equation (26), is
Proof: For obtaining the solution of Corollary 3.2, we follow same procedure, as we used in the proof of Theorem 3.1, and after little simplification, finally we obtain the desired result (29).
Illustration
Example 4.1. To describe solute transport in aquifers, consider the following generalized fractional advection dispersion equation
with initial condition
and boundary condition
where and we consider a dimensionless parameter, called Peclet number, where L is the packing length. The Peclet number determines the nature of the problem, that is, the Peclet number is low for dispersion-dominated problems and is large for advective dominated problems, d is the dispersion coefficient [L2T−1] and ν′ is the Darcy velocity [LT−1].
Our interest is in the solution of (30), for this we follow same procedure, as we applied in the proof of Theorem 3.1, and after little simplification, finally we obtain
Here u(x, t) represent the analytical expression of solute concentration and
Example 4.2. Consider the generalized fractional order space-time advection-dispersion equation
with the initial condition
Here δ(x) is Dirac-delta function and boundary condition
The solution of (34) can be obtained by same technique as we applied in proof of Theorem 3.1
Special Cases
Some interesting special cases of Theorem 3.1 are enumerated below:
If we set γ = 0, in (14), then Hilfer-Prabhakar derivative reduces to Hilfer derivative (12), and the Theorem 3.1 reduces to:
(I). Consider the generalized fractional order space-time advection-dispersion equation of Cauchy type
where (0 < λ ≤ 2), x ∈ R, t ∈ R+, μ ∈ (0, 1), ν ∈ [0, 1],
with initial condition
and boundary condition
For obtaining the solution of (38), follow same procedure as we used in the proof of theorem 3.1, and use (13), after little simplification, obtain the following
Again, use convolution theorem of the Fourier transform to (41), then we get solution of (38), in term of Green's function as
Here Green's function is given as
If we set ν = 1 in (12), then Hilfer fractional derivative reduces to Caputo fractional derivative operator (10) and the equation (38), yields the following:
(II). Consider the generalized fractional order space-time advection-dispersion equation of Cauchy type
where (0 < λ ≤ 2), x ∈ R, t ∈ R+, μ ∈ (0, 1),
with initial condition
and boundary condition
For obtaining the solution of (42), follow same procedure as we used in the proof of theorem 3.1, and use (11), after little simplification, obtain the following
Again, use convolution theorem of the Fourier transform to (45) then we get solution of (42), in term of Green's function as
Here Green's function is given as
(III). On giving suitable value to the parameters involved in Theorem 3.1, we can obtained same results, earlier given by Haung and Liu [14], Haubold et al. [15], Saxena et al. [16], and Agarwal et al. [17].
Conclusion
In this paper, we have presented a solution of generalized space-time fractional advection-dispersion equation. The solution has been developed in terms of Mittag-Leffler function with the help of Sumudu transform and Fourier transform. We can develop the efficient numerical techniques to find solution of various fractional partial differential equations arising in various fields by considering these analytic solutions as base. For future research, the methodology presented in this paper can serve as a good working template to solve any fractional advection-dispersion equations in higher dimensions.
Statements
Author contributions
VG, JS, and YS designed the study, developed the methodology, collected the data, performed the analysis, and wrote the manuscript.
Acknowledgments
The authors are grateful to referees for their suggestions and useful comments on this paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
AtanganaAGómez-AguilarJF. Numerical approximation of Riemann-Liouville definition of fractional derivative: from Riemann-Liouville to Atangana-Baleanu. Numerical Methods for Partial Differential Equations (2018) 34:1502–23. 10.1002/num.22195
2.
GillVModiKSinghY. Analytic solutions of fractional differential equation associated with RLC electrical circuit. J Stat Manage Syst. (2018) 21:575–82. 10.1080/09720510.2018.1466966
3.
Gómez-AguilarJFYépez-MartínezHTorres-JimenezJCordova-FragaTEscobar-JimenezRFOlivares-PeregrinoVH. Homotopy perturbation transform method for nonlinear differential equations involving to fractional operator with exponential kernel. Adv Diff Equations (2017) 2017:68. 10.1186/s13662-017-1120-7,1–18
4.
Morales-DelgadoVFGómez-AguilarJFYépez-MartínezHBaleanuDEscobar-JimenezRFOlivares-PeregrinoVH. Laplace homotopy analysis method for solving linear partial differential equations using a fractional derivative with and without kernel singular. Adv Diff Eq. (2016) 2016:164. 10.1186/s13662-016-0891-6,1–17
5.
Yépez-MartínezHGómez-AguilarJFSosaaIOReyesaJMTorres-JimenezJ. The Feng's first integral method applied to the nonlinear mKdV space-time fractional partial differential equation. Rev Mex Fís. (2016) 62:310–6.
6.
Yépez-MartínezHGómez-AguilarJF. Numerical and analytical solutions of nonlinear differential equations involving fractional operators with power and Mittag-Leffler kernel. Math Modell Nat Phenomena (2018) 13:13. 10.1051/mmnp/2018002
7.
KumarDSinghJBaleanuD. A new analysis of Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler type kernel. Eur J Phys Plus (2018) 133:70. 10.1140/epjp/i2018-11934-y
8.
OwolabiKM. Mathematical analysis and numerical simulation of patterns in fractional and classical reaction-diffusion systems. Chaos Solitons Fractals (2016) 93:89–98. 10.1016/j.chaos.2016.10.005
9.
OwolabiKM. Robust and adaptive techniques for numerical simulation of nonlinear partial differential equations of fractional order. Commun Nonlinear Sci Num Simulat. (2017) 44:304–17. 10.1016/j.cnsns.2016.08.021
10.
OwolabiKMAtanganaA. Analysis of mathematics and numerical pattern formation in superdiffusive fractional multicomponent system. Adv Appl Math Mech. (2017) 9:1438–60. 10.4208/aamm.OA-2016-0115
11.
KumarDSinghJBaleanuDRathoreS. Analysis of a fractional model of Ambartsumian equation, Eur J Phys Plus (2018) 133:259. 10.1140/epjp/i2018-12081-3
12.
SinghJKumarDHammouchZAtanganaA. A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl Math Comput. (2018) 316:504–15. 10.1016/j.amc.2017.08.048
13.
SchumerRBensonDAMeerschaertMWheatcraftSW. Eulerian derivation of the fractional advection-dispersion equation. J Contam Hydrol. (2001) 48:69–88. 10.1016/S0169-7722(00)00170-4
14.
HaungFLiuF. The fundamental solution of the space-time fractional advection- disperssion equation. J Appl Math Comput. (2015) 18:339–50. 10.1007/BF02936577
15.
HauboldHJMathaiAMSaxenaRK. Solution of reaction-diffusion equations in terms of the H- function. Bull Astro Soc India (2007) 35:681–9.
16.
SaxenaRKSaxenaRKallaSL. Solution of the space-time fractional Schrdinger equation occurring in quantum mechanics. Fract Calc Appl Anal. (2010) 13:177–99.
17.
AgarwalRJainSAgarwalRP. Analytic solution of generalized space time advection-dispersion equation with fractional Laplace operator. J Nonlinear Sci Appl. (2016) 9:3445–54. 10.22436/jnsa.009.06.09
18.
WatugalaGK. Sumudu Transform: a new integral transform to solve differential equations and control engineering problems. Int J Math Educ Sci Technol. (1993) 24:35–43. 10.1080/0020739930240105
19.
BelgacemFBM. chapter 6, Applications of the sumudu transform to indefinite periodic parabolic equations. In: Proceedings of the 6th International Conference on Mathematical Problems and Aerospace Sciences (ICNPAA '06), Cambridge: Cambridge Scientific (2007). pp. 51–60.
20.
BelgacemFBMKaraballiAAKallaSL. Analytical investigations of the sumudu transform and applications to integral production equations. Math Problem Eng. (2003) 2003:103–18. 10.1155/S1024123X03207018
21.
KatatbehQDBelgacemFBM. Applications of the Sumudu transform to fractional differential equations. Nonlinear Studies (2011) 18:99–112.
22.
DebnathLBhattaD. Integral Transforms and Their Applications. Boca Raton, FL: CRC Press (1995).
23.
WimanA. ber den Fundamentals in der Teorie der Funktionen Ea(x), (German). Acta Math. (1905) 29:191–201.
24.
PrabhakarTR. A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math J. (1971) 19:7–15.
25.
SamkoSGKilbasAAMarichevOI. Fractional Integrals and Derivatives-Theory and Applications. Linghorne: CRC Press (1993).
26.
CaputoM. Elasticita e dissipazione.Bologana: Zani-chelli (1969).
27.
ChaurasiaVBLSinghJ. Application of sumudu transform in Schrdinger equation occurring in quantum mechanics. Appl Math Sci. (2010) 4:2843–50.
28.
HilferR. Applications of Fractional Calculus in Physics.Singapore: World Scientific (2000). 10.1142/3779
29.
AkahtaniBGulatiVKilicmanA. Application of sumudu transform in generalized fractional reaction-diffusion equation. Int J Appl Comput Math. (2016) 2:387–94. 10.1007/s40819-015-0066-2
30.
GarraRGorenfloRPolitoFTomovskiZ. Hilfer-Prabhakar derivatives and some applications. Appl Math Comput. (2014) 242:576–89. 10.1016/j.amc.2014.05.129
31.
PanchalSKKhandagaleADDolePV. Sumudu transform of Hilfer–Prabhakar fractional derivatives with applications. In: Proceeding of Natural Conference on Recent Trends in Mathematics, Vol.1. pp. 60–66.
32.
BrockmannDSokolovIM. Levy flights in external force fields: from model to equations. Chem. Phys. (2002) 284:409–21. 10.1016/S0301-0104(02)00671-7
Summary
Keywords
space-time fractional advection-dispersion equation, Fourier transforms, Sumudu transforms, Hilfer-Prabhakar fractional derivative, fractional laplacian operator, Mittag-Leffler function
Citation
Gill V, Singh J and Singh Y (2019) Analytical Solution of Generalized Space-Time Fractional Advection-Dispersion Equation via Coupling of Sumudu and Fourier Transforms. Front. Phys. 6:151. doi: 10.3389/fphy.2018.00151
Received
01 October 2018
Accepted
12 December 2018
Published
08 January 2019
Volume
6 - 2018
Edited by
Dumitru Baleanu, University of Craiova, Romania
Reviewed by
Francisco Gomez, Centro Nacional de Investigación y Desarrollo Tecnológico, Mexico; Kolade Matthew Owolabi, Federal University of Technology, Nigeria
Updates
Copyright
© 2019 Gill, Singh and Singh.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jagdev Singh jagdevsinghrathore@gmail.com
This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.