Next Article in Journal
The Legendre Transform in Non-Additive Thermodynamics and Complexity
Next Article in Special Issue
A Novel Numerical Approach for a Nonlinear Fractional Dynamical Model of Interpersonal and Romantic Relationships
Previous Article in Journal
Bodily Processing: The Role of Morphological Computation
Previous Article in Special Issue
Analytical Approximate Solutions of (n + 1)-Dimensional Fractal Heat-Like and Wave-Like Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Approaches to Obtaining the Space-Time Fractional Advection-Diffusion Equation

1
Institute of Mathematics and Computer Sciences, Faculty of Mathematical and Natural Sciences, Jan Długosz University in Czȩstochowa, al. Armii Krajowej 13/15, 42-200 Czȩstochowa, Poland
2
Institute of Law, Administration and Management, Faculty of Philology and History, Jan Długosz University in Czȩstochowa, Zbierskiego 2/4, 42-200 Czȩstochowa, Poland
*
Author to whom correspondence should be addressed.
Submission received: 3 May 2017 / Revised: 13 June 2017 / Accepted: 21 June 2017 / Published: 23 June 2017
(This article belongs to the Special Issue Complex Systems and Fractional Dynamics)

Abstract

:
Two approaches resulting in two different generalizations of the space-time-fractional advection-diffusion equation are discussed. The Caputo time-fractional derivative and Riesz fractional Laplacian are used. The fundamental solutions to the corresponding Cauchy and source problems in the case of one spatial variable are studied using the Laplace transform with respect to time and the Fourier transform with respect to the spatial coordinate. The numerical results are illustrated graphically.

1. Introduction

The Fokker–Planck equation with constant coefficients has the form [1,2,3,4,5]:
w t = a Δ w - v · grad w .
Here, a > 0 is the diffusion coefficient, and v is the drift vector.
The Fokker–Planck equation has many different physical applications: statistics of laser light [1], superionic conductors under the influence of an additional external field [1], diffusion in an external field [3], transport processes in porous media [6,7], groundwater hydrology [8], quantum optics [9], etc. The terms “diffusion coefficient” and “drift vector” do not predetermine their various physical interpretations [3].
Equation (1) can be written in the form of the continuity equation for probability density:
w t = - div J
and the constitutive equation for the probability current [1,3,10,11]:
J = - a grad w + v w .
Another way of obtaining Equation (1) consists of considering the conservation law (for example, for mass concentration) [7,12,13,14]:
D c D t = - div j ,
where:
D D t = t + v ·
is the material derivative, v is a velocity vector and stands for the gradient operator.
Using the Fick law:
j = - a grad c ,
we get:
c t = a Δ c - v · grad c .
Investigation of different physical phenomena in media with complex internal structure has led to considering mathematical models involving differential equations with derivatives of fractional order (see [15,16,17,18,19,20,21,22,23,24,25,26,27], among others). The space fractional [10,28,29,30,31,32,33,34,35,36,37,38,39], time fractional [10,28,29,34,35,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58] and space-time-fractional [59,60,61,62,63,64,65,66] generalizations of the advection-diffusion equation were studied by many authors. Both physical aspects of the fractional advection-diffusion equation and numerical schemes used for their solution were extensively studied in many papers (an extensive literature review can be found in [67]).
Effective implicit numerical methods for the solution of the space-time fractional Fokker–Planck equation and fractional advection diffusion equation were proposed in [68,69]; their stability and convergence were studied in [62] (see also [70,71]).
We can also mention the implicit difference method based on the shifted Grünwald–Letnikov approximation [36], transformation of fractional differential equation into a system of ordinary differential equation [37], the random walk algorithms [38,39], the spectral regularization method [52], the Crank–Nicholson difference scheme [53], Adomian’s decomposition [50], a spatial and temporal discretization [64], the fractional variational iteration method [54], the homotopy perturbation method [51,63,72] and the Jacobi collocation method [73].
In this paper, we discuss two possibilities of obtaining the space-time fractional generalization of the advection-diffusion equation. In the case of the time-fractional advection-diffusion equation, for these possibilities, the terms “Galilei variant” and “Galilei invariant” equations are used [34,42,44]. The Caputo time-fractional derivative and Riesz fractional Laplacian are employed. The fundamental solutions to the corresponding Cauchy and source problems in the case of one spatial variable are studied. It should be emphasized that in the case of the classical advection-diffusion equation ( α = 1 , β = 2 ), the fundamental solutions to the Cauchy problem and to the source problem coincide for t > 0 ; in the case of the fractional advection-diffusion equation, they are substantially different. The properties of the fundamental solution to the Cauchy problem for the space-time fractional advection-diffusion equation in the case of the first approach were investigated in [74]. In that paper, the explicit representation of the fundamental solution for the space fractional advection-diffusion equation ( α = 1 ) was obtained. In the present paper, we supplement the findings of [74] for the Cauchy problem in the case of the first approach by analysis of several particular cases and by the results of numerical calculations. The analytical form of the fundamental solutions for the Cauchy problem in the case of the second approach, as well as of the fundamental solutions to the source problems are obtained in the present paper for the first time. Most attention is concentrated on the solutions of equations with the value of the Caputo derivative α = 1 / 2 , which allows us to obtain solutions in the form of integrals amenable for numerical treatment. The numerical results are illustrated graphically.

2. Mathematical Preliminaries

The Riemann–Liouville integral of fractional order α is defined as [75,76,77]:
I α f ( t ) = 1 Γ ( α ) 0 t ( t - τ ) α - 1 f ( τ ) d τ , α > 0 ,
where Γ ( α ) is the gamma function.
The Riemann–Liouville derivative of fractional order α is introduced as a left-inverse to the fractional integral I α :
D R L α f ( t ) = d n d t n 1 Γ ( n - α ) 0 t ( t - τ ) n - α - 1 f ( τ ) d τ , n - 1 < α < n ,
whereas the Caputo fractional derivative has the form [75,76,77]:
α f ( t ) t α = 1 Γ ( n - α ) 0 t ( t - τ ) n - α - 1 d n f ( τ ) d τ n d τ , n - 1 < α < n .
The introduced fractional operators have the following Laplace transform rules:
L I α f ( t ) = 1 s α f * ( s ) ,
L D R L α f ( t ) = s α f * ( s ) - k = 0 n - 1 D k I n - α f ( 0 + ) s n - 1 - k , n - 1 < α < n ,
L α f t α = s α f * ( s ) - k = 0 n - 1 f ( k ) ( 0 + ) s α - 1 - k , n - 1 < α < n .
Here, the asterisk denotes the transform, and s is the Laplace transform variable.
The one-dimensional Riesz derivative can be defined by its Fourier transform rule [77]:
F d β f ( x ) d | x | β = - | ξ | β f ˜ ( ξ ) , 0 < β 2 ,
where the tilde marks the Fourier transform and ξ is the transform variable. For β = 2 , the standard formula is obtained:
F d 2 f ( x ) d x 2 = - ξ 2 f ˜ ( ξ ) .
In the case of several spatial variables, the positive powers of the Laplace operator, - - Δ β / 2 , β > 0 , are also called the Riesz derivatives and are defined by their Fourier transforms [77,78,79,80]:
F - Δ β / 2 f ( x ) = | ξ | β F f ( x ) , β > 0 ,
where ξ is the transform-variable vector.
Equation (16) is a fractional generalization of the standard formula for the Fourier transform of the Laplace operator corresponding to β = 2 :
F - Δ f ( x ) = | ξ | 2 F f ( x ) .

3. Fractional Advection-Diffusion Equations

3.1. The First Approach

Gurtin and Pupkin [81] proposed the general time-nonlocal dependence between the heat flux and the temperature gradient. Nigmatullin [82,83] considered the following general form of such an equation:
q ( t ) = - k 0 t K ( t - τ ) grad T ( τ ) d τ
resulting in the heat conduction equation with memory:
T t = a T 0 t K ( t - τ ) Δ T ( τ ) d τ .
The constitutive Equation (3) can be generalized in a similar way. We will restrict our consideration to the case v = const . For example, the constitutive equation:
J = 0 t K ( t - τ ) - a grad w ( τ ) + v w ( τ ) d τ
leads to the advection-diffusion equation with the general memory kernel [27,56]:
w t = 0 t K ( t - τ ) a Δ w ( τ ) - v · grad w ( τ ) d τ
(see also [1], where similar equations were obtained for the Fokker–Planck equation in the case of one spatial variable, and the point - was chosen as a starting point (a lower limit) in the integral describing non-Markovian process).
The time-nonlocal constitutive equations for the probability current with the long-tail power kernel [27,56]:
J = D R L 1 - α - a grad w + v w , 0 < α 1 ,
give the time-fractional advection-diffusion equation:
α w t α = a Δ w - v · grad w , 0 < α 1 .
The space-time-fractional advection-diffusion equation takes the form:
α w t α = - a - Δ β / 2 w - v · grad w , 0 < α 1 , 1 β 2 .
The first term in the right-hand side of Equation (24) represents the long-range interaction and provides an attempt to extend the continuum approach to smaller length scales and to link some aspects of lattice mechanics to continuum theory. The left-right side term with the Caputo fractional derivative describes the history effect on the concentration, but there is no memory effect on the drift.

3.2. The Second Approach

In this case, we have the conservation law (4), and instead of the Fick law (6), we assume its time-nonlocal counterpart:
j = - a 0 t K ( t - τ ) grad c ( τ ) d τ ,
which leads to the general equation:
c t + v · grad c = a 0 t K ( t - τ ) Δ c ( τ ) d τ .
For the time-fractional generalization of the Fick law:
j = - a D R L 1 - α grad c , 0 < α 1 ,
we arrive at:
c t + v · grad c = a D R L 1 - α Δ c , 0 < α 1 ,
or:
α c t α = a Δ c - I 1 - α v · grad c , 0 < α 1 .
The corresponding space-time-fractional equation reads:
α c t α = - a - Δ β / 2 c - I 1 - α v · grad c , 0 < α 1 , 1 β 2 .
In Equation (30), the space-fractional Laplace operator (the Riesz operator) describes the long-range interaction, whereas the time-fractional operators refer to the memory effects (both on the concentration and the drift in contrast to Equation (24).

4. Fundamental Solutions to the Cauchy Problems

4.1. The First Approach

Consider the space-time-fractional advection-diffusion Equation (24) in the domain - < x < :
α w t α = a β w | x | β - v w x , 0 < α 1 , 1 β 2 ,
under initial condition:
t = 0 : w = p 0 δ ( x ) .
The Laplace transform with respect to time t and exponential Fourier transform with respect to the spatial coordinate x give:
w ˜ * ( ξ , s ) = p 0 2 π s α - 1 s α + a | ξ | β - i v ξ .
The inverse integral transforms lead to:
w ( x , t ) = p 0 2 π - E α - a | ξ | β - i v ξ t α e - i x ξ d ξ ,
where the following formula has been used [75,76,77]:
L - 1 s α - γ s α + b = t γ - 1 E α , γ - b t α
with E α , γ ( z ) being the Mittag–Leffler function in two parameters α and γ :
E α , γ z = n = 0 z n Γ ( α n + γ ) , α > 0 , γ > 0 , z C ,
and E α ( z ) E α , 1 ( z ) .
When α = 1 , E 1 ( z ) = e z and Equation (34) simplifies:
w ( x , t ) = p 0 π 0 exp - a t ξ β cos ( x - v t ) ξ d ξ .

4.1.1. Standard Diffusion ( α = 1 , β = 2 )

In this case, the solution is well known (see, for example, [1,9]):
w ( x , t ) = p 0 2 π a t exp - ( x - v t ) 2 4 a t .

4.1.2. Cauchy Diffusion with α = 1 , β = 1

For the so-called Cauchy diffusion [84], which is characterized by the values α = β = 1 , we get:
w ( x , t ) = p 0 π a t a 2 t 2 + ( x - v t ) 2 .

4.1.3. Subdiffusion with α = 1 / 2

The Mittag–Leffler function E 1 / 2 ( - z ) has the following integral representation [85]:
E 1 / 2 ( - z ) = 2 π 0 exp - u 2 - 2 u z d u ,
and the solution reads:
w ( x , t ) = p 0 π 3 / 2 0 - exp - u 2 - 2 u a t | ξ | β cos x - 2 u v t ξ d ξ d u .
Taking in (41) β = 2 , we obtain [27]:
w ( x , t ) = p 0 π 2 a t 1 / 4 0 1 u exp - u 2 - x - 2 u v t 2 8 a u t d u .
For β = 1 , we have:
w ( x , t ) = 4 a t p 0 π 3 / 2 0 e - u 2 u 4 a 2 t u 2 + x - 2 u v t 2 d u .
The results of numerical calculations are shown in Figure 1, Figure 2, Figure 3 and Figure 4 for different values of the orders of derivatives α and β and the drift parameter v. In the calculations, we have introduce the following nondimensional quantities:
w ¯ = a 1 / β t α / β p 0 w , x ¯ = 1 a 1 / β t α / β x , v ¯ = t α - α / β a 1 / β v .

4.2. The Second Approach

In this case, we consider the space-time-fractional advection-diffusion equation:
α c t α = a β c | x | β - I 1 - α v c x , 0 < α 1 , 1 β 2 ,
under initial condition:
t = 0 : c = p 0 δ ( x ) .
The integral transform technique allows us to obtain the solution in the transform domain:
c ˜ * ( ξ , s ) = p 0 2 π 1 s + a | ξ | β s 1 - α - i v ξ .
In what follows, we confirm ourselves to the case α = 1 / 2 :
c ˜ * ( ξ , s ) = p 0 2 π 1 s + a | ξ | β s - i v ξ .
The partial-fraction decomposition of (48) gives:
c ˜ * ( ξ , s ) = p 0 2 π 1 a 2 | ξ | 2 β + 4 i v ξ × 1 s + 0 . 5 a | ξ | β - a 2 | ξ | 2 β + 4 i v ξ - 1 s + 0 . 5 a | ξ | β + a 2 | ξ | 2 β + 4 i v ξ
and:
c ( x , t ) = p 0 2 π t - e - i x ξ 1 a 2 | ξ | 2 β + 4 i v ξ × { E 1 / 2 , 1 / 2 - 0 . 5 a | ξ | β - a 2 | ξ | 2 β + 4 i v ξ t - E 1 / 2 , 1 / 2 - 0 . 5 a | ξ | β + a 2 | ξ | 2 β + 4 i v ξ t } d ξ .
Taking into account the integral representation of the Mittag–Leffler function E 1 / 2 , 1 / 2 ( - z ) [85]:
E 1 / 2 , 1 / 2 ( - z ) = 2 π 0 exp - u 2 - 2 u z u d u ,
we get:
c ( x , t ) = p 0 π 3 / 2 t - 0 u r exp - u 2 - a t u | ξ | β × { exp u r t cos ( φ / 2 ) cos x ξ + φ / 2 - u r t sin ( φ / 2 ) - exp - u r t cos ( φ / 2 ) cos x ξ + φ / 2 + u r t sin ( φ / 2 ) } d u d ξ ,
where:
r = a 4 | ξ | 4 β + 16 v 2 ξ 2 , φ = arctan 4 v ξ a 2 | ξ | 2 β .
The solution (52) is shown in Figure 5 and Figure 6 for β = 2 and β = 1 . 5 , respectively. The nondimensional quantities are introduced as:
c ¯ = a 1 / β t α / β p 0 c , w ¯ x ¯ = 1 a 1 / β t α / β x , v ¯ = t 1 - α / β a 1 / β v .

5. Fundamental Solution to the Source Problem

5.1. The First Approach

Consider the space-time-fractional advection-diffusion Equation (24) with the source term in the domain - < x < :
α w t α = a β w | x | β - v w x + q 0 δ ( x ) δ ( t ) , 0 < α 1 , 1 β 2 ,
under zero initial condition:
t = 0 : w = 0 .
The Laplace transform with respect to time t and exponential Fourier transform with respect to the spatial coordinate x give:
w ˜ * ( ξ , s ) = q 0 2 π 1 s α + a | ξ | β - i v ξ .
The inverse integral transforms result in the solution:
w ( x , t ) = q 0 t α - 1 2 π - E α , α - a | ξ | β - i v ξ t α e - i x ξ d ξ .

Subdiffusion with α = 1 / 2

It follows from Equations (51) and (58) that:
w ( x , t ) = q 0 π 3 / 2 t 0 - u exp - u 2 - 2 u a t | ξ | β cos x - 2 u v t ξ d ξ d u
with two particular cases corresponding to β = 2 [27,56]:
w ( x , t ) = q 0 2 a π t 3 / 4 0 exp - u 2 - x - 2 u v t 2 8 a t u u d u
and β = 1 :
w ( x , t ) = 4 a q 0 π 3 / 2 0 e - u 2 u 2 4 a 2 t u 2 + x - 2 u v t 2 d u ,
respectively.
The results of numerical calculations are shown in Figure 7, Figure 8 and Figure 9 for different values of the orders of derivatives and the drift parameter v. In the calculations, we have introduce the following nondimensional quantities:
w ¯ = a 1 / β t α / β - α + 1 q 0 w , x ¯ = 1 a 1 / β t α / β x , v ¯ = t α - α / β a 1 / β v .

5.2. The Second Approach

Consider the space-time fractional advection-diffusion Equation (45) with the source term:
α c t α = a β c | x | β - I 1 - α v c x + q 0 δ ( x ) δ ( t ) , 0 < α 1 , 1 β 2 ,
under zero initial condition:
t = 0 : c = 0 .
The integral transform technique allows us to obtain the solution in the transform domain:
c ˜ * ( ξ , s ) = q 0 2 π s 1 - α s + a | ξ | β s 1 - α - i v ξ .
In what follows, we confirm ourselves to the case α = 1 / 2 :
c ˜ * ( ξ , s ) = q 0 2 π s s + a | ξ | β s - i v ξ .
The partial-fraction decomposition of (66) yields:
c ˜ * ( ξ , s ) = q 0 2 π 1 2 a 2 | ξ | 2 β + 4 i v ξ × a 2 | ξ | 2 β + 4 i v ξ - a | ξ | β s + 0 . 5 a | ξ | β - a 2 | ξ | 2 β + 4 i v ξ + a 2 | ξ | 2 β + 4 i v ξ + a | ξ | β s + 0 . 5 a | ξ | β + a 2 | ξ | 2 β + 4 i v ξ
and:
c ( x , t ) = q 0 4 π t - e - i x ξ 1 a 2 | ξ | 2 β + 4 i v ξ × { a 2 | ξ | 2 β + 4 i v ξ - a | ξ | β E 1 / 2 , 1 / 2 - 0 . 5 a | ξ | β - a 2 | ξ | 2 β + 4 i v ξ t + a 2 | ξ | 2 β + 4 i v ξ + a | ξ | β E 1 / 2 , 1 / 2 - 0 . 5 a | ξ | β + a 2 | ξ | 2 β + 4 i v ξ t } d ξ .
After accounting for the integral representation of the Mittag–Leffler function E 1 / 2 , 1 / 2 ( - x ) (51), we arrive at:
c ( x , t ) = q 0 2 π 3 / 2 t - 0 u r exp - u 2 - a t u | ξ | β × { r exp u r t cos ( φ / 2 ) cos x ξ - u r t sin ( φ / 2 ) + r exp - u r t cos ( φ / 2 ) cos x ξ + u r t sin ( φ / 2 ) - a | ξ | β exp u r t cos ( φ / 2 ) cos x ξ + φ / 2 - u r t sin ( φ / 2 ) + a | ξ | β exp - u r t cos ( φ / 2 ) cos x ξ + φ / 2 + u r t sin ( φ / 2 ) } d u d ξ .
The solution (69) is presented in Figure 10 and Figure 11 with the nondimensional quantities:
c ¯ = a 1 / β t α / β - α + 1 q 0 c , x ¯ = 1 a 1 / β t α / β x , v ¯ = t 1 - α / β a 1 / β v .

6. Discussion

We have considered two approaches to deriving the space-time fractional advection-diffusion equation. In the case of one spatial dimension, we have studied the fundamental solutions to the Cauchy and source problems for the obtained equations. As is seen from the figures, the solutions corresponding to the first and second approach have different behaviors in the direction of drift ( x > 0 in the figures). For α = 1 and 1 β 2 , the solution has no cusp at x = 0 , the quantity v only causes a drift of the maximum value of the solution in the x-direction ( x - v t in the solution (37); the typical curves are shown in Figure 1). For fractional values of α and 1 < β 2 , the fundamental solutions to the Cauchy problems have a cusp at x = 0 . For fractional values of α and β = 1 , the fundamental solution to the Cauchy problem has singularity at x = 0 , and drift caused by the quantity v is less noticeable as is seen from Figure 4.
For the source problem, the quantity v causes drift of the maximum value of the solution in the x-direction with the significant difference between the solutions in the first and second approaches: in the first approach, the maximum value of the solution decreases with the increasing v, whereas in the second approach, the maximum value of the solution increases with the increasing v. The obtained solutions can also be used for testing numerical algorithms for solving the fractional advection-diffusion equation. The reader interested in evaluation of the Mittag–Leffler functions is referred to the paper [86] and the MATLAB program elaborated by Igor Podlubny [87] that implements the algorithms suggested in [86].

Acknowledgments

The support of Jan Długosz University in Czȩstochowa is gratefully acknowledged.

Author Contributions

Yuriy Povstenko wrote the paper. Tamara Kyrylych performed numerical calculations and prepared the corresponding figures. Both authors have equally contributed to the discussion and overall preparation of the manuscript, as well as read and improved the final version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Risken, H. The Fokker-Planck Equation. Methods of Solution and Applications, 2nd ed.; Springer: Berlin, Germany, 1989. [Google Scholar]
  2. Frank, T.D. Nonlinear Fokker-Planck Equations. Fundamentals and Applications; Springer: Berlin, Germany, 2005. [Google Scholar]
  3. Van Kampen, N.G. Stochastic Processes in Physics and Chemistry, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  4. Feller, W. An Introduction to Probability Theory and Its Applications, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1968. [Google Scholar]
  5. Kaviany, M. Heat Transfer Physics; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  6. Kaviany, M. Principles of Heat Transfer in Porous Media, 2nd ed.; Springer: New York, NY, USA, 1995. [Google Scholar]
  7. Nield, D.A.; Bejan, A. Convection in Porous Media, 3rd ed.; Springer: New York, NY, USA, 2006. [Google Scholar]
  8. Rushton, K.R. Groundwater Hydrology. Conceptual and Computational Models; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  9. Carmichael, H.J. Statistical Methods in Quantum Optics. Vol. 1 Master Equations and Fokker-Planck Equations; Springer: Berlin, Germany, 1999. [Google Scholar]
  10. Sokolov, I.M. Thermodynamics and fractional Fokker-Planck equations. Phys. Rev. E 2001, 63, 056111. [Google Scholar] [CrossRef] [PubMed]
  11. Mehrer, H. Diffusion in Solids. Fundamentals, Methods, Materials, Diffusion-Controlled Processes; Springer: Berlin, Germany, 2007. [Google Scholar]
  12. Bejan, A. Convection Heat Transfer, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  13. Bejan, A.; Kraus, A.D. Heat Transfer Handbook; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  14. Brenn, G. Analytical Solutions for Transport Processes. Fluid Mechanics, Heat and Mass Transfer; Springer: Berlin, Germany, 2017. [Google Scholar]
  15. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers, Inc.: Redding, CA, USA, 2006. [Google Scholar]
  16. Povstenko, Y. Fractional heat conduction equation and associated thermal stresses. J. Therm. Stress. 2005, 28, 83–102. [Google Scholar] [CrossRef]
  17. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  18. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  19. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu-Batlle, V. Fractional-Order Systems and Controls. Fundamentals and Applications; Springer: London, UK, 2010. [Google Scholar]
  20. Datsko, B.; Gafiychuk, V. Complex nonlinear dynamics in subdiffusive activator-inhibitor systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1673–1680. [Google Scholar] [CrossRef]
  21. Baleanu, D.; Tenreiro Machado, J.A.; Luo, A.C.J. (Eds.) Fractional Dynamics and Control; Springer: New York, NY, USA, 2012. [Google Scholar]
  22. Valério, D.; Sá da Costa, J. An Introduction to Fractional Control; The Institution of Engineering and Technology: London, UK, 2013. [Google Scholar]
  23. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin, Germany, 2013. [Google Scholar]
  24. Uchaikin, V.; Sibatov, R. Fractional Kinetics in Solids: Anomalous Charge Transport in Semiconductors, Dielectrics and Nanosystems; World Scientific: New Jersey, NJ, USA, 2013. [Google Scholar]
  25. Atanacković, T.M.; Pilipović, S.; Stanković, B.; Zorica, D. Fractional Calculus with Applications in Mechanics: Vibrations and Diffusion Processes; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  26. Herrmann, R. Fractional Calculus: An Introduction for Physicists, 2nd ed.; World Scientific: Singapore, 2014. [Google Scholar]
  27. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  28. Compte, A. Continuous time random walks on moving fluids. Phys. Rev. E 1997, 55, 6821–6831. [Google Scholar] [CrossRef]
  29. Compte, A.; Cáceres, M.O. Fractional dynamics in random velocity fields. Phys. Rev. Lett. 1998, 81, 3140–3143. [Google Scholar] [CrossRef]
  30. Chaves, A.S. A fractional diffusion equation to describe Lévy flights. Phys. Lett. A 1998, 239, 13–16. [Google Scholar] [CrossRef]
  31. Jespersen, S.; Metzler, R.; Fogedby, H.S. Lévy flights in external force fields: Langevin and fractional Fokker-Planck equations and their solutions. Phys. Rev. E 1999, 59, 2736–2745. [Google Scholar] [CrossRef]
  32. Yanovsky, V.C.; Chechkin, A.V.; Schertzer, D.; Tur, A.V. Lévy anomalous diffusion and fractional Fokker-Planck equation. Physica A 2000, 282, 13–34. [Google Scholar] [CrossRef]
  33. Benson, D.A.; Wheatcraft, S.W.; Meerschaert, M.M. Application of a fractional advection-dispersion equation. Water Resour. Res. 2000, 36, 1403–1412. [Google Scholar] [CrossRef]
  34. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  35. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A Math. Gen. 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  36. Meerschaert, M.M.; Tadjeran, C. Finite difference approximations for fractional advection-dispersion flow equations. J. Comput. Appl. Math. 2004, 172, 65–77. [Google Scholar] [CrossRef]
  37. Liu, F.; Anh, V.; Turner, I. Numerical solution of the space fractional Fokker-Planck equation. J. Comput. Appl. Math. 2004, 166, 209–219. [Google Scholar] [CrossRef]
  38. Zhang, Y.; Benson, D.A.; Meerschaert, M.M.; Scheffler, H.-P. On using random walks to solve the space-fractional advection-dispersion equations. J. Stat. Phys. 2006, 123, 89–110. [Google Scholar]
  39. Liu, Q.; Liu, F.; Turner, I.; Anh, V. Approximation of the Lévy-Feller advection-dispersion process by random walk and finite difference method. J. Comput. Phys. 2007, 222, 57–70. [Google Scholar] [CrossRef]
  40. Jumarie, G. A Fokker-Planck equation of fractional order with respect to time. J. Math. Phys. 1992, 33, 3536–3542. [Google Scholar] [CrossRef]
  41. Metzler, R.; Klafter, J.; Sokolov, I.M. Anomalous transport in external fields: Continuous time random walks and fractional diffusion equations extended. Phys. Rev. E 1998, 58, 1621–1633. [Google Scholar] [CrossRef]
  42. Metzler, R.; Barkai, E.; Klafter, J. Anomalous transport in disordered systems under the influence of external fields. Physica A 1999, 266, 343–350. [Google Scholar] [CrossRef]
  43. Metzler, R.; Barkai, E.; Klafter, J. Anomalous diffusion and relaxation close to thermal equilibrium: A fractional Fokker-Planck equation approach. Phys. Rev. Lett. 1999, 82, 3563–3567. [Google Scholar] [CrossRef]
  44. Metzler, R.; Compte, A. Generalized diffusion-advection schemes and dispersive sedimentation: A fractional approach. J. Phys. Chem. B 2000, 104, 3858–3865. [Google Scholar] [CrossRef]
  45. Barkai, E.; Metzler, R.; Klafter, J. From continuous time random walks to the fractional Fokker-Planck equation. Phys. Rev. E 2000, 61, 132–138. [Google Scholar] [CrossRef]
  46. Metzler, R.; Klafter, J. The fractional Fokker-Planck equation: Dispersive transport in an external force field. J. Mol. Liquids 2000, 86, 219–228. [Google Scholar] [CrossRef]
  47. Barkai, E. Fractional Fokker-Planck equation, solution, and application. Phys. Rev. E 2001, 63, 046118. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, F.; Anh, V.; Turner, I.; Zhuang, P. Time-fractional advection-dispersion equation. J. Appl. Math. Comput. 2003, 13, 233–245. [Google Scholar] [CrossRef]
  49. Huang, F.; Liu, F. The time fractional diffusion equation and the advection-dispersion equation. ANZIAM J. 2005, 46, 317–330. [Google Scholar] [CrossRef]
  50. Momani, S. An algorithm for solving the fractional convection–diffusion equation with nonlinear source term. Commun. Nonlinear Sci. Numer. Simul. 2007, 12, 1283–1290. [Google Scholar] [CrossRef]
  51. Momani, S.; Yıldırım, A. Analytical approximate solutions of the fractional convection-diffusion equation with nonlinear source term by He’s homotopy perturbation method. Int. J. Comput. Math. 2010, 87, 1057–1065. [Google Scholar] [CrossRef]
  52. Zheng, G.H.; Wei, T. Spectral regularization method for a Cauchy problem of the time fractional advection-dispersion equation. J. Comput. Appl. Math. 2010, 233, 2631–2640. [Google Scholar] [CrossRef]
  53. Karatay, I.; Bayramoglu, S.R. An efficient scheme for time fractional advection dispersion equations. Appl. Math. Sci. 2012, 6, 4869–4878. [Google Scholar]
  54. Merdan, M. Analytical approximate solutions of fractional convection-diffusion equation with modified Riemann-Liouville derivative by means of fractional variational iteration method. Iran. J. Sci. Technol. 2013, 1, 83–92. [Google Scholar]
  55. Povstenko, Y. Fundamental solutions to time-fractional advection diffusion equation in a case of two space variables. Math. Probl. Eng. 2014, 705364. [Google Scholar] [CrossRef]
  56. Povstenko, Y. Theory of diffusive stresses based on the fractional advection-diffusion equation. In Fractional Calculus: Applications; Abi Zeid Daou, R., Xavier, M., Eds.; NOVA Science Publishers: New York, NY, USA, 2015; pp. 227–242. [Google Scholar]
  57. Povstenko, Y. Generalized boundary conditions for the time-fractional advection diffusion equation. Entropy 2015, 17, 4028–4039. [Google Scholar] [CrossRef]
  58. Povstenko, Y.; Klekot, J. The Dirichlet problem for the time-fractional advection-diffusion equation in a line segment. Bound. Value Probl. 2016, 2016, 89. [Google Scholar] [CrossRef]
  59. Zaslavsky, G.M.; Edelman, M.; Niyazov, B.A. Self-similarity, renormalization, and phase space nonuniformity of Hamiltonian chaotic dynamics. Chaos 1997, 7, 159–181. [Google Scholar] [CrossRef] [PubMed]
  60. Saichev, A.I.; Zaslavsky, G.M. Fractional kinetic equations: solutions and applications. Chaos 1997, 7, 753–764. [Google Scholar] [CrossRef] [PubMed]
  61. Zaslavsky, G.M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: New York, NY, USA, 2005. [Google Scholar]
  62. Liu, F.; Zhuang, P.; Anh, V.; Turner, I.; Burrage, K. Stability and convergence of the difference methods for the space-time fractional advection-diffusion equation. Appl. Math. Comput. 2007, 191, 12–20. [Google Scholar] [CrossRef]
  63. Yıldırım, A.; Koçak, H. Homotopy perturbation method for solving the space-time fractional advection-dispersion equation. Adv. Water Res. 2009, 32, 1711–1716. [Google Scholar] [CrossRef]
  64. Abdel-Rehim, E.A. Explicit approximation solutions and proof of convergence of space-time fractional advection dispersion equations. Appl. Math. 2013, 4, 1427–1440. [Google Scholar] [CrossRef]
  65. Metzler, R.; Jeon, J.-H. Anomalous diffusion and fractional transport equations. In Fractional Dynamics. Recent Advances; Klafter, J., Lim, S.-C., Metzler, R., Eds.; World Scientific: Hackensack, NJ, USA, 2012; pp. 3–32. [Google Scholar]
  66. Povstenko, Y. Space-time-fractional advection diffusion equation in a plane. In Advances in Modelling and Control of Non-Integer Order Systems. Lecture Notes in Electrical Engineering; Latawiec, K.J., Lukaniszyn, M., Stanisławski, R., Eds.; Springer: New York, NY, USA, 2015; Volume 320, pp. 275–284. [Google Scholar]
  67. Zhang, Y.; Benson, D.A.; Reeves, D.M. Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications. Adv. Water Res. 2009, 32, 561–581. [Google Scholar] [CrossRef]
  68. Zhuang, P.; Liu, F.; Anh, V.; Turner, I. Numerical treatment for the fractional Fokker-Planck equation. ANZIAM J. 2007, 48, C759–C774. [Google Scholar] [CrossRef]
  69. Chen, C.; Liu, F.; Turner, I.; Anh, V. Implicit difference approximation of the Galilei invariant fractional advection diffusion equation. ANZIAM J. 2007, 48, C775–C789. [Google Scholar] [CrossRef]
  70. Liu, F.; Zhuang, P.; Burrage, K. Numerical methods and analysis for a class of fractional advection-dispersion models. Comput. Math. Appl. 2012, 64, 2990–3007. [Google Scholar] [CrossRef]
  71. Shen, S.; Liu, F.; Anh, V. Numerical approximations and solution techniques for the space-time Riesz-Caputo fractional advection-diffusion equation. Numer. Algorithms 2011, 56, 383–403. [Google Scholar] [CrossRef]
  72. Panday, R.K.; Singh, O.P.; Baranwal, V.K. An analytic algorithm for the space-time fractional advection-dispersion equation. Comput. Phys. Commun. 2011, 182, 1134–1144. [Google Scholar] [CrossRef]
  73. Parvizi, M.; Eslahchi, M.R.; Dehghan, M. Numerical solution of fractional advection-diffusion equation with a nonlinear source term. Numer. Algorithms 2015, 68, 601–629. [Google Scholar]
  74. Huang, F.; Liu, F. The fundamental solution of the space-time fractional advection-dispersion equation. J. Appl. Math. Comput. 2005, 18, 339–350. [Google Scholar] [CrossRef]
  75. Gorenflo, R.; Mainardi, F. Fractional calculus: Integral and differential equations of fractional order. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: Wien, Austria, 1997; pp. 223–276. [Google Scholar]
  76. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  77. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  78. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives, Theory and Applications; Gordon and Breach: Amsterdam, The Netherlands, 1993. [Google Scholar]
  79. Gorenflo, R.; Mainardi, F.; Moretti, D.; Pagnini, G.; Paradisi, P. Discrete random walk models for space-time fractional diffusion. Chem. Phys. 2002, 284, 521–541. [Google Scholar] [CrossRef]
  80. Matignion, D. Diffusive representations for fractional Laplacian: System theory framework and numerical issues. Phys. Scr. 2009, 136, 014009. [Google Scholar] [CrossRef]
  81. Gurtin, M.E.; Pipkin, A.C. A general theory of heat conduction with finite wave speeds. Arch. Ration. Mech. Anal. 1968, 31, 113–126. [Google Scholar] [CrossRef]
  82. Nigmatullin, R.R. To the theoretical explanation of the “universal response”. Phys. Status Solidi (b) 1984, 123, 739–745. [Google Scholar] [CrossRef]
  83. Nigmatullin, R.R. On the theory of relaxation with remnant temperature. Phys. Status Solidi (b) 1984, 124, 389–393. [Google Scholar] [CrossRef]
  84. Mainardi, F.; Luchko, Y.; Pagnini, G. The fundamental solution of the space-time fractional diffusion equation. Fract. Calc. Appl. Anal. 2001, 4, 153–192. [Google Scholar]
  85. Povstenko, Y. Linear Fractional Diffusion-Wave Equation for Scientists and Engineers; Birkhäuser: New York, NY, USA, 2015. [Google Scholar]
  86. Gorenflo, R.; Loutchko, J.; Luchko, Y. Computation of the Mittag-Leffler function and its derivatives. Fract. Calc. Appl. Anal. 2002, 5, 491–518. [Google Scholar]
  87. Matlab File Exchange 2005, Matlab-Code that Calculates the Mittag-Leffler Function with Desired Accuracy. Available online: www.mathworks.com/matlabcentral/fileexchange/8738-Mittag-Leffler-function (accessed on 17 October 2005).
Figure 1. The fundamental solution to the Cauchy problem; Equation (37); α = 1 , β = 1 . 5 .
Figure 1. The fundamental solution to the Cauchy problem; Equation (37); α = 1 , β = 1 . 5 .
Entropy 19 00297 g001
Figure 2. The fundamental solution to the Cauchy problem; the first approach, Equation (42); α = 0 . 5 , β = 2 .
Figure 2. The fundamental solution to the Cauchy problem; the first approach, Equation (42); α = 0 . 5 , β = 2 .
Entropy 19 00297 g002
Figure 3. The fundamental solution to the Cauchy problem; the first approach, Equation (41); α = 0 . 5 , β = 1 . 5 .
Figure 3. The fundamental solution to the Cauchy problem; the first approach, Equation (41); α = 0 . 5 , β = 1 . 5 .
Entropy 19 00297 g003
Figure 4. The fundamental solution to the Cauchy problem; the first approach, Equation (43); α = 0 . 5 , β = 1 .
Figure 4. The fundamental solution to the Cauchy problem; the first approach, Equation (43); α = 0 . 5 , β = 1 .
Entropy 19 00297 g004
Figure 5. The fundamental solution to the Cauchy problem; the second approach, Equation (52); α = 0 . 5 , β = 2 .
Figure 5. The fundamental solution to the Cauchy problem; the second approach, Equation (52); α = 0 . 5 , β = 2 .
Entropy 19 00297 g005
Figure 6. The fundamental solution to the Cauchy problem; the second approach, Equation (52); α = 0 . 5 , β = 1 . 5 .
Figure 6. The fundamental solution to the Cauchy problem; the second approach, Equation (52); α = 0 . 5 , β = 1 . 5 .
Entropy 19 00297 g006
Figure 7. The fundamental solution to the source problem; the first approach, Equation (60); α = 0 . 5 , β = 2 .
Figure 7. The fundamental solution to the source problem; the first approach, Equation (60); α = 0 . 5 , β = 2 .
Entropy 19 00297 g007
Figure 8. The fundamental solution to the source problem; the first approach, Equation (59); α = 0 . 5 , β = 1 . 5 .
Figure 8. The fundamental solution to the source problem; the first approach, Equation (59); α = 0 . 5 , β = 1 . 5 .
Entropy 19 00297 g008
Figure 9. The fundamental solution to the source problem; the first approach, Equation (61); α = 0 . 5 , β = 1 .
Figure 9. The fundamental solution to the source problem; the first approach, Equation (61); α = 0 . 5 , β = 1 .
Entropy 19 00297 g009
Figure 10. The fundamental solution to the source problem; the second approach, Equation (69); α = 0 . 5 , β = 2 .
Figure 10. The fundamental solution to the source problem; the second approach, Equation (69); α = 0 . 5 , β = 2 .
Entropy 19 00297 g010
Figure 11. The fundamental solution to the source problem; the second approach, Equation (69); α = 0 . 5 , β = 1 . 5 .
Figure 11. The fundamental solution to the source problem; the second approach, Equation (69); α = 0 . 5 , β = 1 . 5 .
Entropy 19 00297 g011

Share and Cite

MDPI and ACS Style

Povstenko, Y.; Kyrylych, T. Two Approaches to Obtaining the Space-Time Fractional Advection-Diffusion Equation. Entropy 2017, 19, 297. https://0-doi-org.brum.beds.ac.uk/10.3390/e19070297

AMA Style

Povstenko Y, Kyrylych T. Two Approaches to Obtaining the Space-Time Fractional Advection-Diffusion Equation. Entropy. 2017; 19(7):297. https://0-doi-org.brum.beds.ac.uk/10.3390/e19070297

Chicago/Turabian Style

Povstenko, Yuriy, and Tamara Kyrylych. 2017. "Two Approaches to Obtaining the Space-Time Fractional Advection-Diffusion Equation" Entropy 19, no. 7: 297. https://0-doi-org.brum.beds.ac.uk/10.3390/e19070297

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