Fractional Calculus and its Applications in Applied Mathematics and other Sciences
View this Special IssueResearch Article  Open Access
On the Generalized Mass Transport Equation to the Concept of Variable Fractional Derivative
Abstract
The hydrodynamic dispersion equation was generalized using the concept of variational order derivative. The modified equation was numerically solved via the CrankNicholson scheme. The stability and convergence of the scheme in this case were presented. The numerical simulations showed that, the modified equation is more reliable in predicting the movement of pollution in the deformable aquifers, than the constant fractional and integer derivatives.
1. Introduction
Anomalous dispersion phenomena are systematically pragmatic in physics, chemistry, and biology fields [1–4]. To distinguish irregular diffusion phenomena, constantorder fractional diffusion equations were initiated and have received breathtaking success. Nonetheless, it has been recognized that the constant order fractional dispersion equations are not competent of typifying some multifaceted dispersion processes, for example, dispersion process in inhomogeneous or heterogeneous medium [5]. In addition, when we think about dispersion process in porous medium, if the medium structure or external field changes with time, in these circumstances, the constantorder fractional diffusion equation model cannot be used to well illustrate such phenomenon [6, 7]. This is the case of the groundwater pollution problem; the medium through which the pollution is dispersing is heterogeneous and change with time. Motionless in some environmental science dispersion processes, the concentration of particles will determine the dispersion pattern [8, 9]. To solve the above problems, the variableorder (VO) fractional dispersion equation models have been suggested for use [10].
The leading edge work of VO operator can be traced to Samko et al. by introducing the variable order integration and RiemannLiouville derivative in [10]. It has been acknowledged as a prevailing modelling approach in the fields of viscoelasticity [10], viscoelastic deformation [11], viscous fluid [12], and anomalous transmission [13]. The problem encounter in groundwater is that the geometry of the aquifer in which the flow or the pollution takes place is not well known. In addition, the geological formation through which the contamination takes place changes in time and space. These phenomena cannot therefore accurately be described neither via the classical hydrodynamic dispersion equation nor via the fractional version of hydrodynamic dispersion equation. In this paper, we generalize the hydrodynamic dispersion equation using the concept of the variable order derivative.
2. Definition and Problem Formulation
The fractional operators (fractional derivatives and integrals) refer to the differential and integral operators of arbitrary order, and fractional differential equations refer to those containing fractional derivatives. The former are the generalization of integerorder differential and integral operators and the latter are the generalization of differential equations of integer order. The operators of variableorder, which fall into a more complex category, are the derivatives and integrals whose orders are the functions of certain variables.
2.1. Variational Order Operator
Definition 1. Left and right RiemannLiouville integrals of variable order: let for all and ; then, is called the left RiemannLiouville integral of variable fractional order , while is referred to as the right RiemannLiouville integral of variable fractional order .
Definition 2. Left and right RiemannLiouville derivatives of variable fractional order: let for all . If , then, the left RiemannLiouville derivative of variable fractional order which is given as is called the left RiemannLiouville derivative of variable fractional order , while in the same line of idea we have the following expression: which is referred to as the right RiemannLiouville derivative of variable fractional order .
Definition 3. Left and right Caputo derivatives of variable fractional order: let for all . If , then, the left RiemannLiouville derivative of variable fractional order is given as
which is called the left Caputo derivative of variable fractional order , while in the same line of idea we have the following expression:
which is referred to as the right Caputo derivative of variable fractional order .
However, this work we will throughout use the following definition.
2.2. Problem Formulation
The concern here is the generalization of advection dispersion equation by including a possible effect of heterogeneity or variability of the aquifer into the mathematical formulation. A onedimensional model consisting of an infinitely long homogeneous isotropic porous medium with a steady state uniform flow with a seepage velocity is considered here. A particular chemical is injected from one end of the model for a period of time such that the input concentration varies as an exponential function of time. The value of that chemical concentration at any time and at a distance from the injection boundary, allowing for the decay and adsorption, may be obtained from the solution of the following set of equations [14]; more details for this model can be found in [15–18]. Consider Subject to the initial and boundary conditions: where is the dispersion coefficient, is the seepage velocity, is the retardation factor, is the radioactive decay constant, is the initial concentration, is a positive constant, and is any source and sink in the system. However, in the case of the groundwater pollution, the function is always neglected because we assume that there is no source and sink in the system under investigation. Therefore, in our case we will set the function to be zero. In order to include explicitly the medium structure or external field changes with time into the mathematical formulation, we replace the standard derivative operator by the variational order derivative operator as follows: Subject to the initial and boundary conditions: where is a continuous function in that can be approximated from the field observation. Note that there is no analytical solution available for this new equation; therefore, our next concern is to establish the existence and uniqueness of such solution. This will be done in the next section.
3. Numerical Solution
Numerical methods acquiesce estimated solutions to the overriding equation all the way through the discretization of space and time. Inside the discredited problem domain, the variable internal properties, boundaries, and stresses of the system are approximated. Deterministic, distributedparameter, and numerical models can relax the rigid idealized conditions of analytical models or lumpedparameter models, and they can therefore be more realistic and flexible for simulating fields’ conditions [19–29]. The finite difference schemes for constantorder time or space fractional diffusion equations have been widely studied [19–29]. Before performing the numerical methods, we assume that (9) has a unique and sufficiently smooth solution. To establish the numerical schemes for the above equation, we let , , , , , and , and we let be the step and be the time size, and and be grid points.
3.1. CrankNicolson Scheme [30]
We bring in the CrankNicolson idea as follows. At the outset, the discretization of first and secondorder space derivative is stated as The CrankNicolson scheme for the VO time fractional diffusion model can be stated as follows: Now replacing (11), (12), (13), and (14) in (9) we obtain the following: Here , , and . It is important to point out that the sum term on the righthand side of (15) automatically vanishes when . Then, (15) can be reformulated as: The above discretization can be rewritten in the following matrix form.
3.2. Stability Analysis of the CrankNicolson Scheme
In this section, we will analyze the stability conditions of the CrankNicolson scheme for the generalized advection dispersion equation.
Let ; here is the approximate solution at the point , and in addition and the function is chosen to be Then, the function can be expressed in Fourier series as follows: It was established by [25] that Observe that, for all , ; in addition, according to the problem in point, the velocity seepage , the dispersion coefficient , the retardation factor , and the radioactive decay constant are positive constants. Then the following properties of the coefficients , , , and can be established as follows.(1), are positive for all .(2) is negative for all .(3) for all .(4), for all .
The error committed while approximating the solution of the generalized advection dispersion equation with CrankNicolson scheme can be presented as follows: If we assume that in (22) can be put in the deltaexponential form as follows: where is a real spatial wave number, replacing (23) in (22) we obtain Equation (24) can be written in the following form: Our next concern here is to show that for all the solution of (25) satisfies the following condition: To achieve this, we make use of the recurrence technique on the natural number .
For and remembering that , are positive for all , then we obtain Assuming that for the property is verified, then Making use of the triangular inequality, we obtain Using the recurrence hypothesis, we have and this completes the proof.
3.3. Convergence Analysis of the CrankNicolson Scheme
If we assume that is the exact solution of our problem at the point , by letting and and substituting this in (18), we obtain Here From (11), (12), and (14), From the above we have that where , , , and are constants. Taking into account Caputo type fractional derivative, the detailed error analysis on the above schemes can refer to the work by Diethelm et al. [31] and further work by Li and Tao [32].
Lemma 4. is true for () where and is a constant. In addition,
This can be achieved via the recurrence technique on the natural number .
When , one has the following:
Now suppose that , . Then
This completes the proof.
Theorem 5. The CrankNicolson scheme is convergent, and there exists a positive constant such that
An interested reader can find the solvability of the CrankNicholson scheme in the work done by [33]. Therefore, the details of the proof will not be presented in this paper.
4. Numerical Simulations
In this section, we present the numerical solutions of the following three cases.
The case where we assume that the variation of the geological formation in which the contamination is migrating can be modelled as function of time and space as . To include explicitly this variability in the mathematical formulation, we reformulate the hydrodynamic dispersion equation as follows: The theoretical parameters used in this simulation are given in Table 1.

Figure 1 shows the approximate solution of the main problem as function of space and time. Figure 2 shows the density plot of the approximate solution of the main problem. Figure 3 shows the variation of the concentration as function of time for a fixed space. Figures 1 and 3 show the solution of the variation order advection and advection equation, respectively, and Figure 2 shows the contour map of variational order advection equation.
5. Discussion and Conclusions
Natural geological deposits with highly contrasting permeability may form mobile and relatively immobile zones, where the potential mass exchange between mobile and immobile zones results in a wide time distribution for solute “trapping.” The transport process, groundwater, is, by its very nature, always in contact with the matrix of an aquifer. There is thus a possibility that the solutes may interact with the rock matrix and one another. A true mathematical model for groundwater pollution must therefore be able to account for interactions between the dissolved solids and matrix of the aquifer. It will thus be advantageous to look at the nature of the interactions between dissolved solids and a porous medium that may be expected in groundwater pollution. Experimental evidence indicates that when a dissolved solid comes in contact with the matrix of a porous medium it may (a) pass through the medium with no apparent effect, (b) be absorbed by the porous matrix, and (c) react with the porous matrix and other substances dissolved in the fluid. The dissolved solids encountered in porous flow are, for this reason, often classified as conservative, nonconservative, and reactive tracers. This behaviour implies that the quantity of dissolved solids in a porous medium depends not only on the flow pattern but also on the nature of the porous matrix and the solution. These situations (a), (b), and (c) can be characterized efficiently by the timenonlocal model, including the time variational order advection dispersion equation. If the highpermeable material tends to form preferential flow paths, such as the interconnected paleochannels observed in alluvial depositional systems, then the solute transport may show a heavy leading edge, which can be described by the Variational Order Advection Dispersion Equation (VOADE). Development of partial differential equations such as the advectiondispersion equation (ADE) begins with assumptions about the random behaviour of a single particle: possible velocities it may experience in a flow field and the length of time it may be immobilized. When assumptions underlying the ADE are relaxed, a fractional ADE (FADE) can arise, with a nonintegerorder derivative on time or space terms. Fractional ADEs are nonlocal; they describe transport affected by hydraulic conditions at a distance. Space fractional ADEs arise when velocity variations are heavy tailed and describe particle motion that accounts for variation in the flow field over the entire system [34]. But when the geological formation, under which this pollution moves, is changing in time and space, the VOADE can be used. Many other uses of the variableorder derivatives can be found in the following works [35–40].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgement
This work was partially financially supported by the Leon Claude Fellowships of 2014 of South Africa. Adem Kilicman gratefully acknowledges that this research was partially supported by the University Putra Malaysia under the Grant Scheme GBIBT having the project number GBIBT/2013/9420100.
References
 T. H. Solomon, E. R. Weeks, and H. L. Swinney, “Observation of anomalous diffusion and Lévy flights in a twodimensional rotating flow,” Physical Review Letters, vol. 71, no. 24, pp. 3975–3978, 1993. View at: Publisher Site  Google Scholar
 S. Bhalekar, V. DaftardarGejji, D. Baleanu, and R. Magin, “Fractional Bloch equation with delay,” Computers & Mathematics with Applications, vol. 61, no. 5, pp. 1355–1365, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. L. Magin, Fractional Calculus in Bioengineering, Begell House Publisher, Redding, Conn, USA, 2006.
 R. L. Magin, O. Abdullah, D. Baleanu, and X. J. Zhou, “Anomalous diffusion expressed through fractional order differential operators in the BlochTorrey equation,” Journal of Magnetic Resonance, vol. 190, no. 2, pp. 255–270, 2008. View at: Publisher Site  Google Scholar
 A. V. Chechkin, R. Gorenflo, and I. M. Sokolov, “Fractional diffusion in inhomogeneous media,” Journal of Physics, vol. 38, no. 42, pp. L679–L684, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Santamaria, S. Wils, E. de Schutter, and G. J. Augustine, “Anomalous diffusion in Purkinje cell dendrites caused by spines,” Neuron, vol. 52, no. 4, pp. 635–648, 2006. View at: Publisher Site  Google Scholar
 H. Sun, W. Chen, and Y. Q. Chen, “Variableorder fractional differential operators in anomalous diffusion modeling,” Physica A, vol. 388, no. 21, pp. 4586–4592, 2009. View at: Publisher Site  Google Scholar
 H. Sun, Y. Q. Chen, and W. Chen, “Randomorder fractional differential equation models,” Signal Processing, vol. 91, no. 3, pp. 525–530, 2011. View at: Publisher Site  Google Scholar
 Y. Q. Chen and K. L. Moore, “Discretization schemes for fractionalorder differentiators and integrators,” IEEE Transactions on Circuits and Systems I, vol. 49, no. 3, pp. 363–367, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 S. Umarov and S. Steinberg, “Variable order differential equations with piecewise constant orderfunction and diffusion with changing modes,” Zeitschrift für Analysis und ihre Anwendungen, vol. 28, no. 4, pp. 431–450, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Ross and S. Samko, “Fractional integration operator of variable order in the Hölder spaces ${H}^{\lambda (x)}$,” International Journal of Mathematics and Mathematical Sciences, vol. 18, no. 4, pp. 777–788, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. T. C. Pedro, M. H. Kobayashi, J. M. C. Pereira, and C. F. M. Coimbra, “Variable order modeling of diffusiveconvective effects on the oscillatory flow past a sphere,” Journal of Vibration and Control, vol. 14, no. 910, pp. 1659–1672, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Ingman and J. Suzdalnitsky, “Application of differential operator with servoorder function in model of viscoelastic deformation process,” Journal of Engineering Mechanics, vol. 131, no. 7, pp. 763–767, 2005. View at: Publisher Site  Google Scholar
 R. W. Cleary and M. J. Ungs, Groundwater Pollution and Hydrology. Mathematical Models and Computer Programs, vol. 08540, Water Resources Programs, Princeton University, Princeton, NJ, USA, 1978.
 M. T. van Genuchten and W. J. Alves, “Analytical solutions of one dimensional convectivedispersive solute transport equations,” United States Department of Agriculture, no. 1661, 1982. View at: Google Scholar
 D. K. Jaiswal, A. Kumar, N. Kumar, and R. R. Yadav, “Analytical solutions for temporally and spatially dependent solute dispersion of pulse type input concentration in onedimensional semiinfinite media,” Journal of HydroEnvironment Research, vol. 2, no. 4, pp. 254–263, 2009. View at: Publisher Site  Google Scholar
 D. K. Jaiswal, A. Kumar, N. Kumar, and M. K. Singh, “Solute transport along temporally and spatially dependent flows through horizontal semiinfinite media: dispersion proportional to square of velocity,” Journal of Hydrologic Engineering, vol. 16, no. 3, pp. 228–238, 2011. View at: Publisher Site  Google Scholar
 R. R. Yadav, D. K. Jaiswal, and G. Gulrana, “Twodimensional solute transport for periodic flow in isotropic porous media: an analytical solution,” Hydrological Processes, vol. 26, no. 22, pp. 3425–3433, 2012. View at: Publisher Site  Google Scholar
 A. Atangana and S. C. Oukouomi Noutchie, “Stability and convergence of a timefractional variable order Hantush equation for a deformable aquifer,” Abstract and Applied Analysis, vol. 2013, Article ID 691060, 8 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Zhang, “A finite difference method for fractional partial differential equation,” Applied Mathematics and Computation, vol. 215, no. 2, pp. 524–529, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Tadjeran, M. M. Meerschaert, and H. P. Scheffler, “A secondorder accurate numerical approximation for the fractional diffusion equation,” Journal of Computational Physics, vol. 213, no. 1, pp. 205–213, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for fractional advectiondispersion flow equations,” Journal of Computational and Applied Mathematics, vol. 172, no. 1, pp. 65–77, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. H. Deng, “Numerical algorithm for the time fractional FokkerPlanck equation,” Journal of Computational Physics, vol. 227, no. 2, pp. 1510–1522, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. P. Li, A. Chen, and J. Ye, “Numerical approaches to fractional calculus and fractional ordinary differential equation,” Journal of Computational Physics, vol. 230, no. 9, pp. 3352–3368, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. M. Chen, F. Liu, I. Turner, and V. Anh, “A Fourier method for the fractional diffusion equation describing subdiffusion,” Journal of Computational Physics, vol. 227, no. 2, pp. 886–897, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. B. Yuste and L. Acedo, “An explicit finite difference method and a new von Neumanntype stability analysis for fractional diffusion equations,” SIAM Journal on Numerical Analysis, vol. 42, no. 5, pp. 1862–1874, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. Podlubny, A. Chechkin, T. Skovranek, Y. Chen, and B. M. Vinagre Jara, “Matrix approach to discrete fractional calculus. II. Partial fractional differential equations,” Journal of Computational Physics, vol. 228, no. 8, pp. 3137–3153, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Atangana, “On the solution of an acoustic wave equation with variableorder derivative loss operator,” Advances in Difference Equations, vol. 2013, article 167, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 R. Lin, F. Liu, V. Anh, and I. Turner, “Stability and convergence of a new explicit finitedifference approximation for the variableorder nonlinear fractional diffusion equation,” Applied Mathematics and Computation, vol. 212, no. 2, pp. 435–445, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Crank and P. Nicolson, “A practical method for numerical evaluation of solutions of partial differential equations of the heatconduction type,” vol. 43, pp. 50–67, 1947. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 K. Diethelm, N. J. Ford, and A. D. Freed, “Detailed error analysis for a fractional Adams method,” Numerical Algorithms, vol. 36, no. 1, pp. 31–52, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. P. Li and C. X. Tao, “On the fractional Adams method,” Computers & Mathematics with Applications, vol. 58, no. 8, pp. 1573–1588, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Cui, “Compact finite difference method for the fractional diffusion equation,” Journal of Computational Physics, vol. 228, no. 20, pp. 7792–7804, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Atangana, A generic assessment of waste disposal at Douala city: principle, practices and uncertainties [Ph.D. thesis], University of the Free State South Africa, Bloemfontein, South Africa, 2012.
 D. Valério and J. Sá da Costa, “Variableorder fractional derivatives and their numerical approximations,” Advances in Fractional Signals and Systems, vol. 91, no. 3, pp. 470–483, 2011. View at: Google Scholar
 D. Valério, G. Vinagre, J. Domingues, and J. Sá da Costa, “Variableorder fractional derivatives and their numerical approximations I—real orders,” in Symposium on Fractional Signals and Systems, Caparica, 2009. View at: Google Scholar
 D. Valério, G. Vinagre, J. Domingues, and J. Sá da Costa, “Variableorder fractional derivatives and their numerical approximations II—complex orders,” in Symposium on Fractional Signals and Systems, Caparica, 2009. View at: Google Scholar
 A. Razminia, A. F. Dizaji, and V. J. Majd, “Solution existence for nonautonomous variableorder fractional differential equations,” Mathematical and Computer Modelling, vol. 55, no. 34, pp. 1106–1117, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Ma, Y. Xu, and W. Yue, “Numerical solutions of a variableorder fractional financial system,” Journal of Applied Mathematics, vol. 2012, Article ID 417942, 14 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Sheng, H. Sun, Y. Chen, and T. Qiu, “Synthesis of multifractional Gaussian noises based on variableorder fractional operators,” Signal Processing, vol. 91, no. 7, pp. 1645–1650, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Abdon Atangana and Adem Kilicman. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.