Next Article in Journal
A Newly Secure Solution to MIMOME OFDM-Based SWIPT Frameworks: A Two-Stage Stackelberg Game for a Multi-User Strategy
Next Article in Special Issue
An Efficient Computational Technique for Fractal Vehicular Traffic Flow
Previous Article in Journal
Using Entropy in Web Usage Data Preprocessing
Previous Article in Special Issue
Chaos in a Cancer Model via Fractional Derivatives with Exponential Decay and Mittag-Leffler Law
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anomalous Advection-Dispersion Equations within General Fractional-Order Derivatives: Models and Series Solutions

1
State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, China
2
School of Mechanics and Civil Engineering, China University of Mining and Technology, Xuzhou 221116, China
3
College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
4
Institute of Geotechnical Engineering, Xi’an University of Technology, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Submission received: 3 November 2017 / Revised: 18 January 2018 / Accepted: 19 January 2018 / Published: 22 January 2018
(This article belongs to the Special Issue Wavelets, Fractals and Information Theory III)

Abstract

:
In this paper, an anomalous advection-dispersion model involving a new general Liouville–Caputo fractional-order derivative is addressed for the first time. The series solutions of the general fractional advection-dispersion equations are obtained with the aid of the Laplace transform. The results are given to demonstrate the efficiency of the proposed formulations to describe the anomalous advection dispersion processes.

1. Introduction

Fractional differential equations have been widely applied to describe the anomalous phenomena in multiple scientific fields, such as physical chemistry, environmental engineering, biology, etc. [1,2,3,4]. In general, the fractional advection-dispersion models with the Caputo fractional derivative [5] may match the real observation better than the classical advection-dispersion models [6,7], which were applied to describe the transport of chemical pollutants in shale gas exploitation [8]. For example, the time- [9] and space- [10] fractional advection-dispersion models have been verified to be able to capture some non-Fickian transport. The tempered advection-dispersion models, such as the promising models, can capture the scale-dependent dispersion (see [5]) and predict the truncated power-law breakthrough curves very nicely (see [11]). In fact, the spatial evolution of the conservative solute molecules with the complex distribution is closely related to the physical and chemical interactions between them and the porous media [12,13,14]. To describe the characteristics of the solute molecules, the entropy solution was used to discuss a class of the fractional degenerate advection-dispersion models in [15]. The fractional transient advection-dispersion model for the entropy density of a reactive plume was discussed in [16].
Recently, a great many new fractional calculus operators have been proposed in [17,18,19,20,21]. For example, the Prabhakar fractional operator via the three-parameter Mittag-Leffler function was suggested in [17]. More recently, a family of the general fractional calculus operators containing the Mittag-Leffler functions were presented in [18]. In particular, Giusti and Colombaro [17] presented how these general fractional calculus operators can reduce to the classical fractional calculus operators. The new general fractional calculus operators of Mittag-Leffler type had been employed to characterize the anomalous relaxation behaviors [19], linear viscoelastic system [20] and anomalous diffusion [21]. The new general Liouville–Caputo fractional derivative operator of Wiman type was considered to describe the anomalous relaxation behavior in [18]. However, the new general Liouville–Caputo fractional derivative operator of Wiman type has not been applied to model the advection-dispersion processes.
In view of the above, the principal objective of this paper is to explore the anomalous advection-dispersion equation describing the pollutant transport in shale gas exploitation within the general Liouville–Caputo fractional derivative operator of Wiman type.
The remainder of the present work is arranged as follows. In Section 2, brief reviews of the definitions of the Mittag-Leffler functions and generations and the new general Liouville–Caputo fractional operator of Wiman type are presented. In Section 3, the series solutions of an anomalous advection-dispersion models via the general Liouville–Caputo fractional derivative operator of Wiman type are obtained. Finally, the conclusions are summarized in Section 4.

2. Mittag-Leffler Function and a New General Liouville–Caputo Fractional Derivative of Wiman Type

In this section, we introduce the family of the Mittag-Leffler function and the new general Liouville–Caputo fractional derivative operator of Wiman type.

2.1. Mittag-Leffler Functions

In this subsection, the definitions of the Mittag-Leffler functions are called (see [1,2,17,18,20,21,22]).
Let C , R and N be the sets of complex numbers, real numbers and positive integers, respectively.
Definition 1.
The Mittag-Leffler function with one parameter, proposed by Gösta Magnus Mittag-Leffler in 1903, is defined as (see [23]):
E ϑ φ = λ = 0 φ λ Γ ϑ λ + 1
where φ , ϑ C , R e ϑ > 0 , λ N and Γ · is the Gamma function.
Definition 2.
In 1905, Wiman extended the Mittag-Leffler function with two parameters, given as [24]:
E ϑ , ϖ φ = λ = 0 φ λ Γ ϑ λ + ϖ
where φ , ϑ , ϖ C , R e ϑ > 0 and λ N .
Definition 3.
A further extension of (1) and (2) with three complex parameters, proposed by Prabhakar in 1971, is defined as [25]:
E ϑ , ϖ τ φ = λ = 0 τ λ Γ ϑ λ + ϖ φ λ λ !
where φ , ϑ , ϖ , τ C , R e ϑ , τ > 0 , λ N , and the Pochhammer symbol is
τ λ = τ τ + 1 ( τ + λ - 1 ) = 1 , λ = 0 Γ τ + λ Γ τ , λ 1
Definition 4.
The Laplace transform of the function f t is defined by [26]:
L f t = f s = 0 f t e - s t d t
where L is the Laplace transform operator.
The Laplace transforms of the family of the Mittag-Leffler functions are listed in Table 1 (see [18,22,27]).

2.2. A New General Liouville–Caputo Fractional-Order Derivative of Wiman Type

Definition 5.
Let ϑ , β , ϖ R , 0 < ϑ < 1 and 0 < β + ϖ < 1 . A new general Liouville–Caputo fractional-order derivative of Wiman type is defined as [18]:
0 L C D t ϑ Ω t = 0 γ γ - t β + ϖ - 1 E ϑ , β + ϖ γ - t ϑ Ω 1 t d t
where
Ω 1 t = d Ω t d t
The Laplace transform of Equation (5) is given as [18]:
L 0 L C D t ϑ Ω t = s - β + ϖ 1 - s - ϑ - 1 s Ω ( s ) - Ω 0
Remark.
If there exists [17]
lim η 0 t β + ϖ - 1 E ϑ , β + ϖ η t ϑ = t β + ϖ - 1 Γ β + ϖ
then the Liouville–Caputo fractional derivative is the special case of the general Liouville–Caputo fractional-order derivative of Wiman type.

3. The Anomalous Advection-Dispersion Model with General Liouville–Caputo Fractional-Order Derivative of Wiman Type

3.1. The Model Background

The unconventional oil or gas (e.g., shale gas) development benefits many countries economically while bringing serious environmental pollution [8,28]. In this process, the fracturing fluids containing many chemical additives will diffuse into the aquifers through the porous media, resulting in the formation contamination [8]. The mathematical model of the advection dispersion of the chemical pollutants in shale gas extraction is shown in Figure 1.
The mathematical model of anomalous advection-dispersion process with general Liouville–Caputo fractional-order derivative of Wiman type is given by
0 L C D t ϑ u x , t = k 2 u x , t x 2 - ρ u x , t x x , t > 0
with the initial value condition
u x , 0 = 0
and the boundary value conditions
u 0 , t = c u , t = B
where u x , t is the concentration of chemical pollutants in the aquifers, k is the dispersion coefficient of the aquifers, ρ is the seepage velocity of chemical pollutants, c is the concentration of chemical pollutants at the lower boundary of aquifers, and B is bounded.

3.2. The Series Solutions for General Fractional Advection-Dispersion Equation

Now, we find the series solutions for the general fractional advection-dispersion model within general Liouville–Caputo fractional-order derivative of Wiman type.
On performing the Laplace transform of Equation (8), we obtain
s - β + ϖ 1 - s - ϑ - 1 × s u ( x , s ) - u x , 0 = k d 2 u ( x , s ) d x 2 - ρ d u ( x , s ) d x
where
L u x , t = u ( x , s )
In a similar manner, the corresponding boundary-value conditions can be written as:
u 0 , s = c s u , s = B
Substituting Equation (9) into Equation (11), we have
M u ( x , s ) = k d 2 u ( x , s ) d x 2 - ρ d u ( x , s ) d x
where M = s 1 - β + ϖ 1 - s - ϑ - 1 .
Making use of the eigenvalue method [8], we obtain the general solution of Equation (13), given as:
u ( x , s ) = m 1 e ρ + ρ 2 + 4 k M 2 k x + m 2 e ρ - ρ 2 + 4 k M 2 k x
With the aid of Equation (12), we have
u ( x , s ) = c s e ρ - ρ 2 + 4 k M 2 k x = c s e ρ 2 k x e - ρ 2 + 4 k M 2 k x
In order to obtain the solution of Equation (8) in the series form, we present
u ( x , s ) = c s e ρ 2 k x e - ρ 2 + 4 k M 2 k x = c s e ρ 2 k x 1 - ρ 2 + 4 k M 4 k 2 1 1 2 2 x + 1 2 ρ 2 + 4 k M 4 k 2 x 2 - 1 6 ρ 2 + 4 k M 4 k 2 3 3 2 2 x 3 + 1 24 ρ 2 + 4 k M 4 k 2 2 x 4 + = c s e ρ 2 k x + c s e ρ 2 k x n = 1 - x n k n n 2 2 n ! Θ n M
where
Θ n M = ρ 2 4 k + M n n 2 2
Furthermore, Θ n M can be expanded as:
Θ n M = ρ 2 4 k + M n n 2 2 = ρ 2 4 k n n 2 2 + n 2 ρ 2 4 k n n 2 2 - 1 M + 1 2 ! n 2 n 2 - 1 ρ 2 4 k n n 2 2 - 2 M 2 + + 1 N ! n 2 n 2 - 1 n 2 - N + 1 ρ 2 4 k n n 2 2 - N M N + = ρ 2 4 k n n 2 2 + N = 1 n 2 n 2 - 1 n 2 - N + 1 N ! ρ 2 4 k n n 2 2 - N M N
The substitution of Equation (17) into Equation (16) results in
u ( x , s ) = c s e ρ - ρ 2 + 4 k M 2 k x = c s e ρ 2 k x + c s e ρ 2 k x n = 1 - x n k n n 2 2 n ! ρ 2 4 k n n 2 2 + N = 1 A M N = c s e ρ 2 k x + c e ρ 2 k x n = 1 - x n k n n 2 2 n ! s - 1 ρ 2 4 k n n 2 2 + N = 1 A s - 1 M N
where
A = n 2 n 2 - 1 n 2 - N + 1 N ! ρ 2 4 k n n 2 2 - N
With the use of Table 1, we have
L - 1 s - 1 M N = t - N 1 - β - ϖ E ϑ , 1 - N 1 - β - ϖ N t ϑ
where L - 1 · represents the inverse Laplace transform operator.
Finally, substituting Equation (19) into Equation (18), we have the series solution of Equation (8) as:
u x , t = c e ρ 2 k x + c e ρ 2 k x n = 1 - x n k n n 2 2 n ! ρ 2 4 k n n 2 2 + N = 1 A t - N 1 - β - ϖ E ϑ , 1 - N 1 - β - ϖ N t ϑ
and the corresponding plots of the changes of the concentration for the different parameters are displayed in Figure 2.
The plot of the changes of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the spaces x = 0.5 , x = 0.8 , x = 0.7 and x = 0.8 is demonstrated in Figure 3.
The plot of the changes of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the times t = 2.5 , t = 2.6 , t = 2.7 , and t = 2.8 is shown in Figure 4.

4. Conclusions and Remarking Comments

In the current work, a novel anomalous advection-dispersion model with general Liouville–Caputo fractional-order derivative of Wiman type was proposed. The series solution of the anomalous advection-dispersion equation was obtained with the aid of the Laplace transform. The fractional advection-dispersion equation within fractional Liouville–Caputo fractional derivative is a special case due to the fact that the Liouville–Caputo fractional derivative is obtained by the general Liouville–Caputo fractional-order derivative of Wiman type. The proposed model is more efficient for the description of the anomalous advection-dispersion process than the classical model with Liouville–Caputo fractional-order derivative [9], which reduces to the classical advection-dispersion [6,7]. The results show that the general Liouville–Caputo fractional-order derivative of Wiman type is important for us to model power-law behaviors in nature.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China (51574219) and the Initiation Fund of Doctor’s Research (No. 107-451117008). We would like to express thanks for the anonymous reviewers to improve the manuscript.

Author Contributions

All authors have worked equally in this manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, X.J.; Machado, J.A.T. A new fractional operator of variable order: Application in the description of anomalous diffusion. Phys. A Stat. Mech. Appl. 2017, 481, 276–283. [Google Scholar] [CrossRef]
  2. Gao, F. General fractional calculus in nonsingular power-law kernel applied to model anomalous diffusion phenomena in heat-transfer problems. Therm. Sci. 2017, 21, S11–S18. [Google Scholar] [CrossRef]
  3. Ionescu, C.; Lopes, A.; Copot, D.; Machado, J.A.T.; Bates, J.H.T. The role of fractional calculus in modelling biological phenomena: A review. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 141–159. [Google Scholar] [CrossRef]
  4. Singh, M.K.; Chatterjee, A.; Singh, V.P. Solution of one-dimensional time fractional Advection dispersion equation by homotopy analysis method. J. Eng. Mech. 2017, 143, 1–16. [Google Scholar] [CrossRef]
  5. Garrard, R.M.; Zhang, Y.; Wei, S.; Sun, H.G.; Qian, J.Z. Can a Time Fractional-Derivative Model Capture Scale-Dependent Dispersion in Saturated Soils? Groundwater 2017, 55, 857–870. [Google Scholar] [CrossRef] [PubMed]
  6. Huang, K.; Toride, N.; Genuchten, M.T.V. Experimental investigation of solute transport in large, homogeneous and heterogeneous, saturated soil columns. Transp. Porous Media 1995, 18, 283–302. [Google Scholar] [CrossRef]
  7. Danckwerts, P.V. Continuous flow systems: Distribution of residence times. Chem. Eng. Sci. 1953, 2, 1–13. [Google Scholar] [CrossRef]
  8. Liang, X.; Liu, G.N.; Su, S.J. Applications of a novel integral transform to the convection-dispersion equations. Therm. Sci. 2017, 21, S233–S240. [Google Scholar] [CrossRef]
  9. Zhang, H.; Liu, F.; Phanikumar, M.S.; Meerschaert, M.M. A novel numerical method for the time variable fractional order mobile–immobile advection—Dispersion model. Comput. Math. Appl. 2013, 66, 693–701. [Google Scholar] [CrossRef]
  10. Benson, D.A.; Wheatcraft, S.W.; Meerschaert, M.M. Application of a fractional advection-dispersion equation. Water Resour. 2000, 36, 1403–1412. [Google Scholar] [CrossRef]
  11. Meerschaert, M.M.; Zhang, Y.; Baeumer, B. Tempered anomalous diffusion in heterogeneous systems. Geophys. Res. Lett. 2008, 35, L17403. [Google Scholar] [CrossRef]
  12. Gorenflo, R.; Mainardi, F.; Moretti, D.; Pagnini, G.; Paradisi, P. Fractional diffusion: Probability distributions and random walk models. Phys. A Stat. Mech. Appl. 2002, 305, 106–112. [Google Scholar] [CrossRef]
  13. Vot, F.L.; Abad, E.; Yuste, S.B. Continuous-time random-walk model for anomalous diffusion in expanding media. Phys. Rev. E 2017, 96, 032117. [Google Scholar] [CrossRef] [PubMed]
  14. Magin, R.L.; Ingo, C. Entropy and information in a fractional order model of anomalous diffusion. IFAC Proc. Vol. 2012, 45, 428–433. [Google Scholar] [CrossRef]
  15. Povstenko, Y. Generalized boundary conditions for the time-fractional advection diffusion equation. Entropy 2015, 17, 4028–4039. [Google Scholar] [CrossRef]
  16. Chiogna, G.; Rolle, M. Entropy-based critical reaction time for mixing-controlled reactive transport. Water Resour. Res. 2017, 53, 7488–7498. [Google Scholar] [CrossRef]
  17. Giusti, A.; Colombaro, I. Prabhakar-like fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 138–143. [Google Scholar] [CrossRef]
  18. Yang, X.J. General fractional calculus operators containing the generalized Mittag-Leffler functions applied to anomalous relaxation. Therm. Sci. 2017, 21, S317–S326. [Google Scholar] [CrossRef]
  19. Yang, X.J. Fractional derivatives of constant and variable orders applied to anomalous relaxation models in heat-transfer problems. Therm. Sci. 2017, 21, 1161–1171. [Google Scholar] [CrossRef]
  20. Yang, X.J. New general fractional-order rheological models with kernels of Mittag-Leffler functions. Rom. Rep. Phys. 2017, 69, 118. [Google Scholar]
  21. Yang, X.J.; Machado, J.A.; Baleanu, D. Anomalous diffusion models with general fractional derivatives within the kernels of the extended Mittag-Leffler type functions. Rom. Rep. Phys. 2017, 69, 115. [Google Scholar]
  22. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: New York, NY, USA, 2014; pp. 7–16, 83,98. [Google Scholar]
  23. Mittag-Leffler, G.M. Sur la nouvelle function Eαx. C. R. Acad. Sci. 1903, 137, 554–558. [Google Scholar]
  24. Wiman, A. Über den Fundamentalsatz in der Teorie der Funktionen Eαx. Acta Math. 1905, 29, 191–201. [Google Scholar] [CrossRef]
  25. Prabhakar, T.R. A Singular Integral Equation with a Generalized Mittag-Leffler Function in the Kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  26. Liang, X.; Gao, F.; Su, S.J.; Wang, Z.; Yang, X.J. Classifications and duality relations for several integral transforms. J. Nonlinear Sci. Appl. 2017, 10, 6324–6332. [Google Scholar] [CrossRef]
  27. Mathai, A.M.; Haubold, H.J. Special Functions for Applied Scientists; Springer: New York, NY, USA, 2008; p. 95. [Google Scholar]
  28. Torres, L.; Yadav, O.P.; Khan, E. A review on risk assessment techniques for hydraulic fracturing water and produced water management implemented in onshore unconventional oil and gas production. Sci. Total Environ. 2016, 539, 478–493. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The mathematical model of the advection dispersion of the chemical pollutants in shale gas extraction involving general Liouville–Caputo fractional-order derivative of Wiman type.
Figure 1. The mathematical model of the advection dispersion of the chemical pollutants in shale gas extraction involving general Liouville–Caputo fractional-order derivative of Wiman type.
Entropy 20 00078 g001
Figure 2. The concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 .
Figure 2. The concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 .
Entropy 20 00078 g002
Figure 3. The changes of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the spaces x = 0.5 , x = 0.6 , x = 0.7 and x = 0.8 .
Figure 3. The changes of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the spaces x = 0.5 , x = 0.6 , x = 0.7 and x = 0.8 .
Entropy 20 00078 g003
Figure 4. The change plot of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the times t = 2.5 , t = 2.6 , t = 2.7 , and t = 2.8 .
Figure 4. The change plot of the concentration of chemical pollutants in the aquifers for the values n = 1 , n = 2 , n = 3 , n = 4 and 5, and the parameters N = 5 , k = 0.8 , c = 0.5 , ρ = 2 , β = 0.4 , ϖ = 0.5 , ϑ = 0.2 in the times t = 2.5 , t = 2.6 , t = 2.7 , and t = 2.8 .
Entropy 20 00078 g004
Table 1. The Laplace transforms of Mittag-Leffler functions with power-law functions.
Table 1. The Laplace transforms of Mittag-Leffler functions with power-law functions.
Mittag-Leffler Functions with Power-Law FunctionsLaplace Transforms
t β - 1 E ϑ , β t ϑ R e ϑ , β > 0 s - β 1 - s - ϑ - 1
t β - 1 E ϑ , β τ t ϑ R e ϑ , β , τ > 0 s - β 1 - s - ϑ - τ
t β + ϖ - 1 E ϑ , β + ϖ t ϑ R e ϑ , β , ϖ > 0 s - β + ϖ 1 - s - ϑ - 1
t β + ϖ - 1 E ϑ , β + ϖ τ t ϑ R e ϑ , β , ϖ , τ > 0 s - β + ϖ 1 - s - ϑ - τ

Share and Cite

MDPI and ACS Style

Liang, X.; Yang, Y.-G.; Gao, F.; Yang, X.-J.; Xue, Y. Anomalous Advection-Dispersion Equations within General Fractional-Order Derivatives: Models and Series Solutions. Entropy 2018, 20, 78. https://0-doi-org.brum.beds.ac.uk/10.3390/e20010078

AMA Style

Liang X, Yang Y-G, Gao F, Yang X-J, Xue Y. Anomalous Advection-Dispersion Equations within General Fractional-Order Derivatives: Models and Series Solutions. Entropy. 2018; 20(1):78. https://0-doi-org.brum.beds.ac.uk/10.3390/e20010078

Chicago/Turabian Style

Liang, Xin, Yu-Gui Yang, Feng Gao, Xiao-Jun Yang, and Yi Xue. 2018. "Anomalous Advection-Dispersion Equations within General Fractional-Order Derivatives: Models and Series Solutions" Entropy 20, no. 1: 78. https://0-doi-org.brum.beds.ac.uk/10.3390/e20010078

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