Abstract

An alternative construction for the space-time fractional diffusion-advection equation for the sedimentation phenomena is presented. The order of the derivative is considered as , for the space and time domain, respectively. The fractional derivative of Caputo type is considered. In the spatial case we obtain the fractional solution for the underdamped, undamped, and overdamped case. In the temporal case we show that the concentration has amplitude which exhibits an algebraic decay at asymptotically large times and also shows numerical simulations where both derivatives are taken in simultaneous form. In order that the equation preserves the physical units of the system two auxiliary parameters and are introduced characterizing the existence of fractional space and time components, respectively. A physical relation between these parameters is reported and the solutions in space-time are given in terms of the Mittag-Leffler function depending on the parameters and . The generalization of the fractional diffusion-advection equation in space-time exhibits anomalous behavior.

1. 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 [19]. 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 [1013]. 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 heavy-tailed. When defined as a walk in a space of dimension greater than one, the steps are defined in terms of the step-lengths, 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 and the first-order time derivative with a Caputo derivative of order . 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 , 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 , 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.

2. Fundamentals of Fractional Calculus

The most commonly used definitions in fractional calculus are Riemann-Liouville (RL), Grünwald-Letnnikov definition (GL), and Caputo fractional derivative (CFD) [3].

The RL definition of the fractional derivative for is

The GL definition is formulated as where is the time increment.

The CFD for a function is defined as follows: where and . In this case, 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 ; 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

3. 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 , for the fractional DAE in space and time domain, respectively.

3.1. Fractional Space Diffusion-Advection Equation

Considering (11) and assuming that the space derivative is fractional (9) and the time derivative is ordinary, the spatial fractional equation is

A particular solution of this equation may be found in the form where is the natural frequency and is a constant.

Substituting (13) in (12) we obtain

The solution of (14) may be obtained applying (4). Taking the solution (13) we have

For underdamped case, with , . Considering and , in (15), is the wave vector and is the damping factor.

When , from (15) we have (16) represents the classical case in the underdamped case. From (16) we see that there is a relation between and given by

Then, the solution (15) for the underdamped case or takes the form where is a dimensionless parameter.

Due to the condition 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.

In the overdamped case, or , the solution of (15) has the form

When , from (21) we have where is the initial concentration in . The solution (22) represents the change of the concentration . This represents the classical case.

Take into account the fact that the relation between and is

The solution (21) takes the form where is a dimensionless parameter.

Due to the condition , 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.

3.2. Fractional Time Diffusion-Advection Equation

Considering (11) and assuming that the time derivative is fractional (10) and the space derivative is ordinary, the temporal fractional equation is

A particular solution of this equation may be found in the form where is the wave vector in the direction and is a constant.

Substituting (28) into (27) we obtain if 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 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 is the Mittag-Leffler function.

When , we have, , substituting this expression in (30) we have where solving for we obtain and for 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 given by the order of the fractional differential equation

We can use this relation to write (43) as where is a dimensionless parameter.

When , , from (36) we have where and the 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 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 arbitrarily chosen and Figure 4 shows the concentration in the time for different values of also chosen arbitrarily.

3.3. 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 , , 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.

4. 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 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 ; (16) represents the classical case in the underdamped case. In the overdamped case, (24) for represents the overdamped case.

For the temporal case (36) shows the solution for the fractional time DAE. When , 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 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).

The solutions (15) and (36) correspond to space-time generalized solutions of the DAE. Figures 5 and 6 show the case where both derivatives are considered (time-space) simultaneous; besides, it was shown that when and are less than 1, the concentration behaves like a concentration 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 , 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 , which allows describing superdiffusion phenomena (including ballistic diffusion); the Levy flights model and subwave phenomena will be made in a future paper.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank Mayra Martínez and Professor Dumitru Baleanu for the interesting discussions. The authors acknowledge PAPIIT-DGAPA-UNAM (IN112212) and CONACYT-0167485 for the economic support. José Francisco Gómez Aguilar acknowledges the support provided by CONACYT through the assignment postdoctoral fellowship, convocatory 2013.