Space-Time Fractional Diffusion-Advection Equation with Caputo Derivative

and Applied Analysis 3 where C is the concentration, D is the diffusion coefficient, and λ is the drift velocity; this equation considered only the distribution of one cartesian component x. To be consistent with dimensionality and following [24] we introduce an auxiliary parameter σ x and σ t as follows:


Introduction
The Diffusion-Advection Equation (DAE) describes the evolution of a concentration profile due to diffusion and advection simultaneously; this equation describes physical phenomena where concentrations as mass, energy, or other physical quantities are transferred inside a physical system due to two contributions: diffusion and convection, in this equation the concentration-dependent diffusion coefficient.The fractional calculus (FC) is the generalization of derivatives and integrals to noninteger orders; the mathematical formulation was developed for Fourier, Liouville, Abel, Riemann, Lacroix, Grünwald, Riesz, among many others.Many problems in physical science electromagnetism, electrochemistry, diffusion, and general transport theory can be solved by the fractional calculus approach [1][2][3][4][5][6][7][8][9].The FC may be approached via the theory of linear differential equations.FC involves nonlocal operators which can be applied to physical systems yielding new information about their behavior.Modeling as fractional order proves to be useful particulary for systems where the memory plays a significant role, in comparison with the ordinary calculus models; this is the main advantage [10][11][12][13].The process of sedimentation of particles dispersed in a fluid is one of great practical importance, but it has always proved extremely difficult to examine theoretically.The sedimentation phenomena describe the response of the system to the action of an external force (usually centrifugal).
The hydrodynamical problem of one particle falling through a fluid has been solved by Einstein, for Smoluchowski and many others [14,15].A random walk is a mathematical formalization of a path that consists of a succession of random steps.Random walks are related to the diffusion models and within the fractional approach it is possible to include external fields; considerations of transport in the phase space are possible within the same approach [16].A Lévy flight, also referred to as Lévy motion, is a random walk in which the step-lengths have a probability distribution that is heavytailed.When defined as a walk in a space of dimension greater than one, the steps are defined in terms of the steplengths, which have a certain probability distribution; the steps made are in isotropic random directions [17].Scher and Montroll [18] in their description of anomalous transit-time dispersion in amorphous solids considered a system where the traditional methods proved to fail.The authors presented a stochastic transport model for the current; the dynamics considered continuous time random walk approach for a variety of physical quantities in numerous experimental realizations.Jespersen et al. in [19] considered Lévy flights subject to external force fields; the authors presented a Riesz/Weyl form of the DAE; the corresponding Fokker-Planck equation contains a fractional derivative in space.Mainardi et al. in [20] present the interpretation of the corresponding Green function as a probability density for the particular cases of space fractional, time fractional, and neutral fractional diffusion; the fundamental space-time fractional diffusion equation is obtained from the standard diffusion equation by replacing the second-order space derivative with a Riesz-Feller derivative of order 0 <  ≤ 2 and the first-order time derivative with a Caputo derivative of order 0 <  ≤ 2. Recently, the fractional kinetic equations of the diffusion, DAE, and Fokker-Planck type are presented in [21]; the equations are derived from basic random walk models.Metzler and Compte in [22] consider advection processes following anomalous statistics on an effective transport level within the framework of fractional kinetic equations.Liu et al. [23] present the space fractional Fokker-Planck equation, the authors use the Riemann-Liouville and Grünwald-Letnikov definitions of fractional derivatives, the Fokker-Planck equation is transformed into a system of ordinary differential equations, and numerical results are presented.In the works of the authors mentioned above the pass from an ordinary derivative to a fractional one is direct, to be consistent with the dimensionality of the fractional differential equations (FDE) in the work [24]; the authors have proposed a systematic way to construct FDE for the physical systems analyzing the dimensionality of the ordinary derivative operator and trying to bring it to a fractional derivative operator consistently.Following [24], in this work we propose an alternative procedure for constructing the fractional DAE; the order of the derivative being considered is 0 < ,  ≤ 1 for space-time domain, respectively.The main purpose in this paper was to show an alternative solution of the DAE using the fractional derivative of Caputo type; this representation preserves the physical units of the system for any value taken by the exponent of the fractional derivative.This alternative solution in the range 0 < ,  ≤ 1 describes Levy flights (non-Markovian version) and the phenomena of subdiffusion for space-time domain, respectively.The paper is organized as follows: second section introduction to fractional calculus, third section the fractional DAE, and the fourth section the conclusions.
The RL definition of the fractional derivative for ( > 0) is The GL definition is formulated as () where  is the time increment.
The CFD for a function () is defined as follows: where  = 1, 2, . . .∈  and −1 <  ≤ .In this case, 0 <  ≤ 1 is the order of the CFD.In this definition the derivative of a constant is zero and the initial conditions for the fractional order differential equations take the same form as for the integer-order ones taking a known physical interpretation, which is very suitable in engineering [4,6,9,24].In this paper we will use the CFD.Laplace transform to CFD is given by [3] The Mittag-Leffler function has caused extensive interest among physicists due to its vast potential of applications describing realistic physical systems with memory and delay.The inverse Laplace transform requires the introduction of the Mittag-Leffler function defined by the series expansion as [3] where Γ is the gamma function, when ,  = 1; from (5) we obtain   , the exponential function as a special case of the Mittag-Leffler function.This function will often appear with the argument  = −  ; its Laplace transform then is given as Some common Laplace transforms are

Fractional Diffusion-Advection Equation
The DAE [21] is represented by where  is the concentration,  is the diffusion coefficient, and  is the drift velocity; this equation considered only the distribution of one cartesian component .
To be consistent with dimensionality and following [24] we introduce an auxiliary parameter   and   as follows: this is true if the parameter   has dimensions of length (inverse meters) and   has the dimension of time (inverse seconds).
The fractional representation of ( 8) is The order of the derivative being considered is 0 < ,  ≤ 1 for the fractional DAE in space and time domain, respectively..Considering (11) and assuming that the space derivative is fractional (9) and the time derivative is ordinary, the spatial fractional equation is
When  = 1, from (15) we have (16) represents the classical case in the underdamped case.
Due to the condition  < 2 √  we can choose an example So, the solution (15) takes its final form Figure 1 shows the concentration in the space for different values of  arbitrarily chosen.

Take into account the fact that the relation between 𝛽 and
The solution (21) takes the form where x = (( 2 /4 2 ) − (/)) 1/2  is a dimensionless parameter.
Due to the condition  > 2 √ , we have the following range of values: Then, the solution ( 21) can be written in its final form Figure 2 shows the concentration in the space for different values of  arbitrarily chosen..Considering (11) and assuming that the time derivative is fractional (10) and the space derivative is ordinary, the temporal fractional equation is

Fractional Time Diffusion-Advection Equation
A particular solution of this equation may be found in the form where ⃗  is the wave vector in the  direction and  0 is a constant.
Substituting (28) into ( 27) we obtain is the natural frequency with real and imaginary parts  and , we have where is the natural frequency in the medium in presence of fractional time components, and ⃗  2 is the natural frequency without its presence.
Substituting (33) in (29) gives The solution is written as The particular solution of the equation ( 35) as follows where   (− ⃗ ω2   ) is the Mittag-Leffler function.When  = 1, we have, ⃗  =  − , substituting this expression in (30) we have where solving for  we obtain and for (40) substitutitng ( 40) into (39) we have Equations ( 40) and (41) describe the real and imaginary part of the natural frequency in terms of the wave vector , the diffusion coefficient , and the drift velocity .From (36) we have The solution to (27) which follows from (42) is given by where ⃗  =  − , is the natural frequency,  and  is given by ( 40) and (41), respectively.The equation (42) represent the classical case.
In these case exists a physical relation between the natural frequency ⃗ , the parameter   and the period  0 given by the order  of the fractional differential equation We can use this relation to write (43) as where t = ⃗  is a dimensionless parameter.
When  = 1/2, ⃗ ω2 = ⃗  2  1/2  , from (36) we have where and the erfc() denotes the error function defined as For large values of , the error function can be approximated as substituting ( 48) into (46) leads to the solution At asymptotically large times and using (49) we have The solution (52) represents plane waves with time decaying amplitude.The asymptotic behavior of Mittag-Leffler function [25], for  = − ⃗ ω2   is Then, substituting (53) in (36) gives where is the natural frequency in the medium in presence of fractional time components, the parameter   characterizes these structures (components that show an intermediate behavior between a system conservative and dissipative) of the fractional time operator.In this case, (52) and (54) represent the time evolution of the concentration and the amplitude which exhibits an algebraic decay for  → ∞.The fractional differentiation with respect to time can be interpreted as an existence of memory effects which correspond to intrinsic dissipation in the system [6].Figure 3 shows the simulation of (45) for  = 0.9 arbitrarily chosen and Figure 4 shows the concentration in the time for different values of  also chosen arbitrarily.

Fractional Space-Time Diffusion-Advection Equation.
Now considering (11) and assuming that the space and time derivatives are fractional, the order of the time-space fractional differential equation is 0 < ,  ≤ 1,   has dimension of length, and   has dimension of time.Figures 5 and 6 show the simulation of ( 20)-( 45) and ( 26)-( 45) for different values of  and  arbitrarily chosen.

Conclusions
We present the analysis of the DAE from the point of view of fractional calculus.FDE have been examined separately; the fractional spatial derivative and the temporal fractional derivative have finally shown the numerical simulations where both derivatives are taken in simultaneous form.
The parameters   and   are introduced characterizing x, t  = 1,  = 1  = 0.9,  = 0.9  = 0.8,  = 0.8 the existence of the fractional space and time components, respectively; these parameters represent components that show an intermediate behavior between a conservative and dissipative system.The general solutions of the fractional DAE depending only on the parameters  and  are given in the form of the multivariate Mittag-Leffler functions preserving the physical units of the system studied.
It was shown in the spatial case that the solution (15) corresponds to the spatial generalized solution of the DAE, in the case  = 1; (16) represents the classical case in the underdamped case.In the overdamped case, (24) for  = 1 represents the overdamped case.
For the temporal case (36) shows the solution for the fractional time DAE.When  = 1, the natural frequency ⃗  is represented by ⃗  =  − , where  and  are given by ( 40) and (41), respectively, and describes the real and imaginary part of the natural frequency in terms of the wave vector , the diffusion coefficient , and the drift velocity ; (43) represents the classical case.In this case exists a physical relation between the parameter   and the period  0 given by the order  of the fractional differential equation, relation (44).An important measure for sedimentation processes is the speed with which the major part of the tracer is setting down; in the subdiffusive case, the drift velocity  is only an effective velocity and it is explicitly time dependent; (52) and (54) represent the time evolution of the concentration; these equations involve the drift velocity and the amplitude which exhibit an algebraic decay for  → ∞ and exhibit anomalous slow diffusion (subdiffusion).
with spatial-temporal-decaying amplitude with respect to time  and the space .The numerical simulations represent a nonlocal concentration interpreted as an existence of memory effects which correspond to intrinsic dissipation characterized by the exponent of the fractional derivatives  and  in the system and related to concentration in a fractal space-time geometry.
This alternative solution of DAE using the fractional derivative of Caputo type in the range 0 < ,  ≤ 1 may contribute to the study and interpretation of the electrokinetic phenomena, the models of porous electrodes, the groundwater dynamics, and the description of anomalous complex processes (mass transport by diffusion and convection considering the formation of different intermediate species).
In the case of Caputo's derivative it is defined in the range of 1 < ,  ≤ 2 which allows describing superdiffusion phenomena (including ballistic diffusion); the Levy flights model and subwave phenomena will be made in a future paper.