Impact Factor 2.638 | CiteScore 2.3
More on impact ›

# Frontiers in Physics

## ORIGINAL RESEARCH article

Front. Phys., 08 January 2019 | https://doi.org/10.3389/fphy.2018.00151

# Analytical Solution of Generalized Space-Time Fractional Advection-Dispersion Equation via Coupling of Sumudu and Fourier Transforms

Vinod Gill1, Jagdev Singh2* and Yudhveer Singh3
• 1Department of Mathematics, Government Post Graduate College Hisar, Hisar, India
• 2Department of Mathematics, JECRC University, Jaipur, India
• 3Amity Institute of Information Technology, Amity University, Jaipur, India

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 [16]. 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 [712] 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 ${E}_{\rho ,0,\omega ,0+}^{0}f=f.$ 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

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] xR, tR+, μ ∈ (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], xR, tR+, μ ∈ (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 ${\mu }^{\prime }=\frac{d}{{\nu }^{\prime }L}$ and we consider a dimensionless parameter, called Peclet number, $Pe=\frac{1}{{\mu }^{\prime }}$ 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 $g\left(k\right)=\frac{1}{\sqrt{2\pi }}\left[\frac{{e}^{-\left(1+ik\right)}-1}{1+ik}\right].$

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), xR, tR+, μ ∈ (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), xR, tR+, μ ∈ (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.

## Author Contributions

VG, JS, and YS designed the study, developed the methodology, collected the data, performed the analysis, and wrote the manuscript.

## Conflict of Interest Statement

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.

## Acknowledgments

The authors are grateful to referees for their suggestions and useful comments on this paper.

## References

1. Atangana A, Gómez-Aguilar JF. 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. doi: 10.1002/num.22195

2. Gill V, Modi K, Singh Y. Analytic solutions of fractional differential equation associated with RLC electrical circuit. J Stat Manage Syst. (2018) 21:575–82. doi: 10.1080/09720510.2018.1466966

3. Gómez-Aguilar JF, Yépez-Martínez H, Torres-Jimenez J, Cordova-Fraga T, Escobar-Jimenez RF, Olivares-Peregrino VH. Homotopy perturbation transform method for nonlinear differential equations involving to fractional operator with exponential kernel. Adv Diff Equations (2017) 2017:68. doi: 10.1186/s13662-017-1120-7,1–18

4. Morales-Delgado VF, Gómez-Aguilar JF, Yépez-Martínez H, Baleanu D, Escobar-Jimenez RF, Olivares-Peregrino VH. 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. doi: 10.1186/s13662-016-0891-6,1–17

5. Yépez-Martínez H, Gómez-Aguilar JF, Sosaa IO, Reyesa JM, Torres-Jimenez J. 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ínez H, Gómez-Aguilar JF. Numerical and analytical solutions of nonlinear differential equations involving fractional operators with power and Mittag-Leffler kernel. Math Modell Nat Phenomena (2018) 13:13. doi: 10.1051/mmnp/2018002

CrossRef Full Text

7. Kumar D, Singh J, Baleanu D. A new analysis of Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler type kernel. Eur J Phys Plus (2018) 133:70. doi: 10.1140/epjp/i2018-11934-y

8. Owolabi KM. Mathematical analysis and numerical simulation of patterns in fractional and classical reaction-diffusion systems. Chaos Solitons Fractals (2016) 93:89–98. doi: 10.1016/j.chaos.2016.10.005

9. Owolabi KM. Robust and adaptive techniques for numerical simulation of nonlinear partial differential equations of fractional order. Commun Nonlinear Sci Num Simulat. (2017) 44:304–17. doi: 10.1016/j.cnsns.2016.08.021

10. Owolabi KM, Atangana A. Analysis of mathematics and numerical pattern formation in superdiffusive fractional multicomponent system. Adv Appl Math Mech. (2017) 9:1438–60. doi: 10.4208/aamm.OA-2016-0115

11. Kumar D, Singh J, Baleanu D, Rathore S. Analysis of a fractional model of Ambartsumian equation, Eur J Phys Plus (2018) 133:259. doi: 10.1140/epjp/i2018-12081-3

12. Singh J, Kumar D, Hammouch Z, Atangana A. A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl Math Comput. (2018) 316:504–15. doi: 10.1016/j.amc.2017.08.048

13. Schumer R, Benson DA, Meerschaert M, Wheatcraft SW. Eulerian derivation of the fractional advection-dispersion equation. J Contam Hydrol. (2001) 48:69–88. doi: 10.1016/S0169-7722(00)00170-4

14. Haung F, Liu F. The fundamental solution of the space-time fractional advection- disperssion equation. J Appl Math Comput. (2015) 18:339–50. doi: 10.1007/BF02936577

15. Haubold HJ, Mathai AM, Saxena RK. Solution of reaction-diffusion equations in terms of the H- function. Bull Astro Soc India (2007) 35:681–9.

16. Saxena RK, Saxena R, Kalla SL. Solution of the space-time fractional Schrdinger equation occurring in quantum mechanics. Fract Calc Appl Anal. (2010) 13:177–99.

17. Agarwal R, Jain S, Agarwal RP. Analytic solution of generalized space time advection-dispersion equation with fractional Laplace operator. J Nonlinear Sci Appl. (2016) 9:3445–54. doi: 10.22436/jnsa.009.06.09

18. Watugala GK. Sumudu Transform: a new integral transform to solve differential equations and control engineering problems. Int J Math Educ Sci Technol. (1993) 24:35–43. doi: 10.1080/0020739930240105

19. Belgacem FBM. 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. Belgacem FBM, Karaballi AA, Kalla SL. Analytical investigations of the sumudu transform and applications to integral production equations. Math Problem Eng. (2003) 2003:103–18. doi: 10.1155/S1024123X03207018

21. Katatbeh QD, Belgacem FBM. Applications of the Sumudu transform to fractional differential equations. Nonlinear Studies (2011) 18:99–112.

22. Debnath L, Bhatta D. Integral Transforms and Their Applications. Boca Raton, FL: CRC Press (1995).

23. Wiman A. ber den Fundamentals in der Teorie der Funktionen Ea(x), (German). Acta Math. (1905) 29:191–201.

24. Prabhakar TR. A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math J. (1971) 19:7–15.

25. Samko SG, Kilbas AA, Marichev OI. Fractional Integrals and Derivatives-Theory and Applications. Linghorne: CRC Press (1993).

26. Caputo M. Elasticita e dissipazione. Bologana: Zani-chelli (1969).

27. Chaurasia VBL, Singh J. Application of sumudu transform in Schrdinger equation occurring in quantum mechanics. Appl Math Sci. (2010) 4:2843–50.

28. Hilfer R. Applications of Fractional Calculus in Physics. Singapore: World Scientific (2000). doi: 10.1142/3779

CrossRef Full Text

29. Akahtani B, Gulati V, Kilicman A. Application of sumudu transform in generalized fractional reaction-diffusion equation. Int J Appl Comput Math. (2016) 2:387–94. doi: 10.1007/s40819-015-0066-2

30. Garra R, Gorenflo R, Polito F, Tomovski Z. Hilfer-Prabhakar derivatives and some applications. Appl Math Comput. (2014) 242:576–89. doi: 10.1016/j.amc.2014.05.129

31. Panchal SK, Khandagale AD, Dole PV. 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. Brockmann D, Sokolov IM. Levy flights in external force fields: from model to equations. Chem. Phys. (2002) 284:409–21. doi: 10.1016/S0301-0104(02)00671-7

Keywords: space-time fractional advection-dispersion equation, Fourier transforms, Sumudu transforms, Hilfer-Prabhakar fractional derivative, fractional laplacian operator, Mittag-Leffler function, 2010 Mathematics Subject Classification:, 26A33, 33E12, 34A08, 42A38, 49K20

2010 Mathematics Subject Classification: 26A33, 33E12, 34A08, 42A38, 49K20.

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.

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

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