Next Article in Journal
Coarse Graining Shannon and von Neumann Entropies
Next Article in Special Issue
Meromorphic Non-Integrability of Several 3D Dynamical Systems
Previous Article in Journal
Kinetics of Interactions of Matter, Antimatter and Radiation Consistent with Antisymmetric (CPT-Invariant) Thermodynamics
Previous Article in Special Issue
Entropy in Natural Time and the Associated Complexity Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Diffusion in a Solid with Mass Absorption

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
3
Institute of Preschool and School Education, Faculty of Pedagogy, Jan Długosz University in Czȩstochowa, Waszyngtona 4/8, 42-200 Czȩstochowa, Poland
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(5), 203; https://doi.org/10.3390/e19050203
Submission received: 28 March 2017 / Revised: 24 April 2017 / Accepted: 29 April 2017 / Published: 2 May 2017
(This article belongs to the Special Issue Complex Systems, Non-Equilibrium Dynamics and Self-Organisation)

Abstract

:
The space-time-fractional diffusion equation with the Caputo time-fractional derivative and Riesz fractional Laplacian is considered in the case of axial symmetry. Mass absorption (mass release) is described by a source term proportional to concentration. The integral transform technique is used. Different particular cases of the solution are studied. The numerical results are illustrated graphically.

1. Introduction

The conventional theory of diffusion is based on the classical Fick law, which relates the matter flux j to the concentration gradient
j = k grad c ,
where k is the diffusion conductivity.
In combination with the balance equation for mass, the classical Fick law leads to the standard diffusion equation
c t = a Δ c
with a being the diffusivity coefficient.
Mass transport in a medium with a first order chemical reaction is described by an additional linear source term in the diffusion equation [1]:
c t = a Δ c b c ,
where the values of the coefficient b > 0 and b < 0 correspond to mass absorption and mass release, respectively. Similarly, Equation (3) can also represent heat conduction with heat release or absorption [2]. This equation also governs heat transfer in a thin plate whose lateral surfaces exchange heat with the ambient medium having constant temperature; the values b > 0 and b < 0 correspond to the temperature of the medium greater and less than that of the plate, respectively [3]. Similar equations appear in the theory of bio-heat transfer [4,5,6] and in the survival probability (see [7] and references therein).
Nonclassical theories, in which the Fick law and the standard diffusion equation are replaced by more general equations, constantly attract the attention of researchers. From a mathematical point of view, the Fick law in the theory of diffusion, the Fourier law in the theory of heat conduction, and the Darcy law in the theory of fluid flow through a porous medium are identical. Some of the generalized theories are formulated in terms of diffusion, others in terms of heat conduction or fluid flow through a porous medium.
The time-nonlocal dependence between the matter flux j and the concentration gradient grad c with the long-tail power kernel (see [8,9,10,11]) can be interpreted in terms of fractional integrals and derivatives and results in the time-fractional diffusion equation
α c t α = a Δ c ,
where α c / t α is the Caputo fractional derivative [12,13,14]:
α f ( t ) t α = 1 Γ ( n α ) 0 t ( t τ ) n α 1 d n f ( τ ) d τ n d τ , n 1 < α < n ,
and Γ ( α ) is the gamma function.
The Caputo fractional derivative has the following Laplace transform rule
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, s is the Laplace transform variable.
Space nonlocal generalizations of the Fick law with the power kernel can also be interpreted in terms of fractional calculus and give the space-fractional diffusion equation
c t = a Δ β / 2 c ,
where the positive powers of the Laplace operator Δ β / 2 , β > 0 , are called the Riesz derivatives.
The space-fractional heat conduction equation in the case of one spatial coordinate was considered by Gorenflo and Mainardi [15], in the case of higher dimensions by Hanyga [16]. The definitions of space-fractional differential operators can be found, for example, in [14,17,18,19,20].
The cumbersome aspects of space-fractional differential operators disappear when one computes their Fourier transforms:
F Δ β / 2 f ( x ) = | ξ | β F f ( x ) , β > 0 ,
where ξ is the transform-variable vector.
It is obvious that Equation (8) 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 ) .
If the considered function f ( x , y ) depends only on the radial coordinate r = x 2 + y 2 1 / 2 , then the two-fold Fourier transform
f ˜ ( ξ , η ) = 1 2 π f ( x , y ) e i ( x ξ + y η ) d x d y
can be simplified (see [11,21]) and is expressed as
F f ( x , y ) = f ˜ ( ξ , η ) = H f ( r ) = f ^ ( ϱ ) = 0 r f ( r ) J 0 ( r ϱ ) d r ,
F 1 f ˜ ( ξ , η ) = f ( x , y ) = H 1 f ^ ( ϱ ) = f ( r ) = 0 ϱ f ^ ( ρ ) J 0 ( r ϱ ) d ϱ ,
where ϱ = | ξ | = ξ 2 + η 2 . Hence, in the case of axial symmetry the two-fold Fourier transform with respect to the Cartesian coordinates x and y is reduced to the Hankel transform with respect to the radial coordinate r. Here J 0 ( r ) is the Bessel function of the first kind of the zeroth order, the tilde marks the Fourier transform, the hat denotes the Hankel transform.
Equation (8) for the Fourier transform of the fractional Laplace operator in the case of axial symmetry also simplifies:
F Δ β / 2 f ( x , y ) = H Δ β / 2 f ( r ) = ϱ β f ^ ( ϱ ) .
If the time nonlocality is accompanied with the spatial nonlocalty, then the general space-time-fractional heat conduction equation is obtained (see [18,19,22,23]):
α c t α = a Δ β / 2 c .
It should be emphasized that fractional calculus (the theory of derivatives and integrals of non-integer order) has many important applications in description of processes in media with complex internal structure: in fractional dynamics [24,25,26,27], fractional kinetics [28,29,30], fractional thermoelasticity [8,9,31,32], biology [33,34], fractional control [35,36,37], description of physical processes in colloid, glassy and porous materials [38,39,40], comb structures [41,42], dielectrics and semiconductors [30,43,44], among others.
The fractional counterpart of Equation (3) has the following form:
α c t α = a Δ β / 2 c b c .
Physical aspects of fractional reaction-diffusion were extensively studied in the literature (see, for example, [45,46,47,48,49] and references therein).
In this paper, we study the fundamental solutions to the Cauchy problem and to the source problem for Equation (15) in the case of axial symmetry. The integral transform technique is used. The numerical results are illustrated graphically.

2. Fundamental Solution to the Cauchy Problem

We consider the space-time-fractional diffusion equation with mass absorption (mass release) in the axisymmetric case:
α c t α = a Δ β / 2 c b c , 0 < α 1 , 1 β 2 ,
under the initial condition
t = 0 : c = p 0 2 π r δ ( r R ) ,
where δ ( r ) is the Dirac delta function.
The constant multiplier p 0 has been introduced in Equation (17) to get the nondimensional quantities used in numerical calculations (see Equation (30)).
The Laplace transform with respect to time t and the Hankel transform with respect to the radial coordinate r give the solution in the transform domain
c ^ * = p 0 2 π s α 1 s α + a ϱ β + b J 0 R ϱ .
After inversion of the integral transforms, we get:
c ( r , t ) = p 0 2 π 0 E α a ϱ β + b t α J 0 R ϱ J 0 r ϱ ϱ d ϱ ,
where E α ( z ) is the Mittag-Leffler function in one parameter α [12,14,50]:
E α z = n = 0 z n Γ ( α n + 1 ) , α > 0 , z C ,
and the following formula for the inverse Laplace transform
L 1 s α 1 s α + b = E α b t α
has been used.
The asymptotic behavior of the solution is determined by the asymptotic behavior of the Mittag-Leffler function E α ( x ) . For 0 < α < 1 , we have [51,52,53] :
E α a ϱ β + b t α exp a ϱ β + b t α Γ ( 1 + α ) for small values of t ,
E α a ϱ β + b t α 1 Γ ( 1 α ) 1 a ϱ β + b t α for large values of t .
Consider several particular cases of the obtained solution corresponding to different particular cases of the Mittag-Leffler function.

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

Since E 1 ( x ) = e x , taking into account integral (A1) from Appendix A, we get
c ( r , t ) = p 0 4 π a t I 0 r R 2 a t exp r 2 + R 2 4 a t b t .

2.2. Localized Diffusion ( α 0 , β = 2 )

In this case E 0 ( x ) = 1 1 + x 2 and using Equation (A2), we obtain
c ( r , t ) = p 0 2 π a I 0 r ( 1 + b ) a K 0 R ( 1 + b ) a , 0 r R , I 0 R ( 1 + b ) a K 0 r ( 1 + b ) a , R r < .

2.3. Subdiffusion with α = 1 / 2 and β = 2

The Mittag-Leffler function E 1 / 2 ( x ) can be written in the following integral form [11]:
E 1 / 2 ( x ) = 2 π 0 exp u 2 2 u x d u .
From Equations (19) and (26) with taking into account (A1), we have
c ( r , t ) = p 0 4 π 3 / 2 a t 1 / 2 0 I 0 r R 4 a t 1 / 2 u × exp u 2 2 b t 1 / 2 u r 2 + R 2 8 a t 1 / 2 u 1 u d u .

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

For the so called Cauchy diffusion [22], which is characterized by the values α = β = 1 , evaluating the corresponding integral (see Equation (A3) from Appendix A), we arrive at
c ( r , t ) = p 0 a t 8 π 2 e b t k 3 E ( k ) ( 1 k 2 ) ( r R ) 3 / 2 ,
where E ( k ) is the complete elliptic integral of the first kind,
k = 2 r R ( a t ) 2 + ( r + R ) 2 .
Results of numerical calculations are shown in Figure 1, Figure 2, Figure 3 and Figure 4. We have introduced the following nondimensional quantities
c ¯ = 2 π R 2 p 0 c , r ¯ = r R , t ¯ = a t α R β , b ¯ = t α b .

3. Fundamental Solution to the Source Problem

In this case, we study the equation
α c t α = a Δ β / 2 c b c + q 0 2 π r δ ( t ) δ ( r R ) , 0 < α 1 , 1 β 2 ,
under zero initial condition
t = 0 : c = 0 .
The Laplace transform with respect to time t and the Hankel transform with respect to the radial coordinate r give the solution in the transform domain
c ^ * = q 0 2 π 1 s α + a ϱ β + b J 0 R ϱ .
After inversion of the integral transforms, we arrive at
c ( r , t ) = q 0 t α 1 2 π 0 E α , α a ϱ β + b t α J 0 R ϱ J 0 r ϱ ϱ d ϱ ,
where E α , γ ( z ) is the Mittag-Leffler function in two parameters α and γ [12,14,50]:
E α , γ z = n = 0 z n Γ ( α n + γ ) , α > 0 , γ > 0 , z C ,
and
L 1 s α γ s α + b = t γ 1 E α , γ b t α .

Subdiffusion with α = 1 / 2 and β = 2

The Mittag-Leffler function E 1 / 2 , 1 / 2 ( x ) is presented as [11]:
E 1 / 2 , 1 / 2 ( x ) = 2 π 0 exp u 2 2 u x u d u .
Hence
c ( r , t ) = p 0 4 π 3 / 2 a t 1 / 2 0 I 0 r R 4 a t 1 / 2 u × exp u 2 2 b t 1 / 2 u r 2 + R 2 8 a t 1 / 2 u d u .
The results of numerical calculations are shown in Figure 5 and Figure 6, where the nondimensional concentration is defined as
c ¯ = 2 π R 2 q 0 t α 1 c ,
other nondimensional quantities are the same as in (30).

4. Discussion

We have considered the fundamental solutions to the Cauchy problem and to the source problem for the space-time fractional diffusion equation with the linear source term. Time nonlocality deals with memory effects, whereas space nonlocality describes the long-range interaction. It is seen from Figures that for positive values of the parameter b (mass absorption) the mass concentration is decreased, for negative values of the parameter b (mass release) the mass concentration is increased. For α = 1 , the shape of curves depends on the value of nondimensional time t ¯ (Figure 1 presents typical results for small values of t ¯ , Figure 2—for t ¯ 1 ). For fractional value of α , the fundamental solution to the Cauchy problem has a cusp at r ¯ = 1 . Such a cusp also appears in the curves describing the fundamental solution to the source problem for decreasing value of the order of fractional Laplacian (compare Figure 4 and Figure 6). To calculate the Mittag-Leffler functions E α ( z ) in the fundamental solution (19) and E α , β ( z ) in the fundamental solution (34) we have used the algorithms proposed in [54].

Author Contributions

Yuriy Povstenko and Tamara Kyrylych wrote the paper; Grażyna Rygał carried out numerical calculations for the fundamental solution to the source problem and prepared the corresponding Figures; Tamara Kyrylych performed numerical calculations for the Cauchy problem and prepared the corresponding Figures. All the authors have equally contributed in 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.

Appendix A. Integrals

In Appendix, we present integrals used in the paper. Equations (A1) and (A2) are taken from [55], the Lipshitz-Hankel integral (A3) was evaluated in [56].
0 e a x 2 J 0 ( b x ) J 0 ( c x ) x d x = 1 2 a exp b 2 + c 2 4 a I 0 b c 2 a , a > 0 , b > 0 , c > 0 .
0 x x 2 + a 2 J 0 ( b x ) J 0 ( c x ) d x = I 0 ( a b ) K 0 ( a c ) , 0 < b < c , I 0 ( a c ) K 0 ( b c ) , 0 < c < b , a > 0 .
0 e a x J 0 ( b x ) J 0 ( c x ) x d x = a k 3 4 π ( 1 k 2 ) ( b c ) 3 / 2 E ( k ) , a > 0 , b > 0 , c > 0 ,
where E ( k ) is is the complete elliptic integral of the first kind,
k = 2 b c a 2 + ( b + c ) 2 .

References

  1. Crank, J. The Mathematics of Diffusion, 2nd ed.; Oxford University Press: Oxford, UK, 1975. [Google Scholar]
  2. Carslaw, H.S.; Jaeger, J.C. Conduction of Heat in Solids, 2nd ed.; Oxford University Press: Oxford, UK, 1959. [Google Scholar]
  3. Polyanin, A.D. Handbook of Linear Partial Differential Equations for Engineers and Scientists; Chapman & Hall/CRC: Boca Raton, FL, USA, 2002. [Google Scholar]
  4. Pennes, H.H. Analysis of tissue and arterial blood temperatures in the resting human forearm. J. Appl. Physiol. 1948, 1, 93–122. [Google Scholar] [PubMed]
  5. Nyborg, W.L. Solutions of the bio-heat transfer equation. Phys. Med. Biol. 1988, 33, 785–792. [Google Scholar] [CrossRef] [PubMed]
  6. Lakhssassi, A.; Kengne, E.; Semmaoui, H. Investigation of nonlinear temperature distribution in biological tissues by using bioheat transfer equation of Pennes’ type. Nat. Sci. 2010, 2, 131–138. [Google Scholar] [CrossRef]
  7. Abad, E.; Yuste, S.B.; Lindenberg, K. Survival probability of an immobile target in a sea of evanescent diffusive or subdiffusive traps: A fractional equation approach. Phys. Rev. E 2012, 86, 061120. [Google Scholar] [CrossRef] [PubMed]
  8. Povstenko, Y. Fractional heat conduction equation and associated thermal stresses. J. Therm. Stresses 2005, 28, 83–102. [Google Scholar] [CrossRef]
  9. Povstenko, Y. Theory of thermoelasticity based on the space-time-fractional heat conduction equation. Phys. Scr. 2009, 136, 014017. [Google Scholar] [CrossRef]
  10. Povstenko, Y. Non-axisymmetric solutions to time-fractional diffusion-wave equation in an infinite cylinder. Fract. Calc. Appl. Anal. 2011, 14, 418–435. [Google Scholar] [CrossRef]
  11. Povstenko, Y. Linear Fractional Diffusion-Wave Equation for Scientists and Engineers; Birkhäuser: New York, NY, USA, 2015. [Google Scholar]
  12. 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]
  13. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  14. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  15. Gorenflo, R.; Mainardi, F. Random walk models for space-fractional diffusion processes. Fract. Calc. Appl. Anal. 1998, 1, 167–191. [Google Scholar]
  16. Hanyga, A. Multidimensional solutions of space-fractional diffusion equations. Proc. R. Soc. Lond. A 2001, 457, 2993–3005. [Google Scholar] [CrossRef]
  17. 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]
  18. Saichev, A.I.; Zaslavsky, G.M. Fractional kinetic equations: Solutions and applications. Chaos 1997, 7, 753–764. [Google Scholar] [CrossRef] [PubMed]
  19. 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]
  20. Matignion, D. Diffusive representations for fractional Laplacian: System theory framework and numerical issues. Phys. Scr. T 2009, 136, 014009. [Google Scholar] [CrossRef]
  21. Sneddon, I.N. The Use of Integral Transforms; McGraw-Hill: New York, NY, USA, 1972. [Google Scholar]
  22. 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]
  23. Hanyga, A. Multidimensional solutions of space-time-fractional diffusion equations. Proc. R. Soc. Lond. A 2002, 458, 429–450. [Google Scholar] [CrossRef]
  24. 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]
  25. Datsko, B.; Gafiychuk, V. Complex nonlinear dynamics in subdiffusive activator-inhibitor systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1673–1680. [Google Scholar] [CrossRef]
  26. Baleanu, D.; Tenreiro Machado, J.A.; Luo, A.C.J. (Eds.) Fractional Dynamics and Control; Springer: New York, NY, USA, 2012. [Google Scholar]
  27. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Higher Education Press: Bejing, China; Springer: Berlin, Germany, 2010. [Google Scholar]
  28. Herrmann, R. Fractional Calculus: An Introduction for Physicists, 2nd ed.; World Scientific: Singapore, 2014. [Google Scholar]
  29. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin, Germany, 2013. [Google Scholar]
  30. Uchaikin, V.; Sibatov, R. Fractional Kinetics in Solids: Anomalous Charge Transport in Semiconductors, Dielectrics and Nanosystems; World Scientific: Hackensack, NJ, USA, 2013. [Google Scholar]
  31. Povstenko, Y. Thermoelasticity based on fractional heat conduction equation. In Proceedings of the 6th International Congress on Thermal Stresses, Vienna, Austria, 26–29 May 2005; Ziegler, F., Heuer, R., Adam, C., Eds.; Vienna University of Technology: Vienna, Austria, 2005; Volume 2, pp. 501–504. [Google Scholar]
  32. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  33. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers, Inc.: Redding, CA, USA, 2006. [Google Scholar]
  34. Nonnenmacher, T.F.; Metzler, R. Applications of fractional calculus techniques to problems in biophysics. In Applications of Fractional Calculus in Physics; Hilfer, R., Ed.; World Scientific: Singapore, 2000; pp. 377–428. [Google Scholar]
  35. Caponetto, R.; Dongola, G.; Fortuna, L.; Petráš, I. Fractional Order Systems. Modeling and Control Applications; World Scientific: Hackensack, NJ, USA, 2010. [Google Scholar]
  36. 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]
  37. Valério, D.; Sá da Costa, J. An Introduction to Fractional Control; The Institution of Engineering and Technology: London, UK, 2013. [Google Scholar]
  38. Weeks, E.R.; Weitz, D.A. Subdiffusion and the cage effect studied near the colloidal glass transition. Chem. Phys. 2002, 284, 361–377. [Google Scholar] [CrossRef]
  39. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  40. Kimmich, R. Strange kinetics, porous media, and NMR. Chem. Phys. 2002, 284, 253–285. [Google Scholar] [CrossRef]
  41. Arkhincheev, V.E. Anomalous diffusion and charge relaxation on comb model: Exact solutions. Phys. A Stat. Mech. Appl. 2000, 280, 304–314. [Google Scholar] [CrossRef]
  42. Arkhincheev, V.E. Diffusion on random comb structure: Effective medium approximation. Phys. A Stat. Mech. Appl. 2002, 307, 131–141. [Google Scholar] [CrossRef]
  43. Nigmatullin, R.R. To the theoretical explanation of the “universal response”. Phys. Status Solidi B 1984, 123, 739–745. [Google Scholar] [CrossRef]
  44. Nigmatullin, R.R. On the theory of relaxation with remnant temperature. Phys. Status Solidi B 1984, 124, 389–393. [Google Scholar] [CrossRef]
  45. Sokolov, I.M.; Schmidt, M.G.W.; Sagués, F. Reaction-subdiffusion equations. Phys. Rev. E 2006, 73, 031102. [Google Scholar] [CrossRef] [PubMed]
  46. Henry, B.I.; Langlands, T.A.M.; Wearne, S.L. Anomalous diffusion with linear reaction dynamics: From continuous time random walks to fractional reaction-diffusion equations. Phys. Rev. E 2006, 74, 031116. [Google Scholar] [CrossRef] [PubMed]
  47. Abad, E.; Yuste, S.B.; Lindenberg, K. Reaction-subdiffusion and reaction-superdiffusion equations for evanescent particles performing continuous-time random walks. Phys. Rev. E 2010, 81, 031115. [Google Scholar] [CrossRef] [PubMed]
  48. Méndez, V.; Fedotov, S.; Horsthemke, W. Reaction-Transport Systems: Mesoscopic Foundations, Fronts, and Spatial Instabilities; Springer: Berlin, Germany, 2010. [Google Scholar]
  49. Yuste, S.B.; Abad, E.; Lindenberg, K. Reactions in subdiffusive media and associated fractional equations. In Fractional Dynamics. Recent Advances; Klafter, J., Lim, S.C., Metzler, R., Eds.; World Scientific: Hackensack, NJ, USA, 2012; pp. 77–106. [Google Scholar]
  50. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  51. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1955; Volume 3. [Google Scholar]
  52. 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]
  53. 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]
  54. 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]
  55. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Volume 2: Special Functions; Gordon and Breach: Amsterdam, The Netherlands, 1986. [Google Scholar]
  56. Eason, G.; Noble, B.; Sneddon, I.N. On certain integrals of Lipschitz-Hankel type involving products of Bessel functions. Philos. Trans. R. Soc. Lond. Ser. A 1955, 247, 529–551. [Google Scholar] [CrossRef]
Figure 1. Dependence of the fundamental solution to the Cauchy problem on distance for α = 1 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Figure 1. Dependence of the fundamental solution to the Cauchy problem on distance for α = 1 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Entropy 19 00203 g001
Figure 2. Dependence of the fundamental solution to the Cauchy problem on distance for α = 1 , β = 2 , t ¯ = 1 and various values of b ¯ .
Figure 2. Dependence of the fundamental solution to the Cauchy problem on distance for α = 1 , β = 2 , t ¯ = 1 and various values of b ¯ .
Entropy 19 00203 g002
Figure 3. Dependence of the fundamental solution to the Cauchy problem on distance for α = 0 . 5 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Figure 3. Dependence of the fundamental solution to the Cauchy problem on distance for α = 0 . 5 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Entropy 19 00203 g003
Figure 4. Dependence of the fundamental solution to the Cauchy problem on distance for α = 0 . 5 , t ¯ = 0 . 25 , b ¯ = 0 . 5 and various values of β .
Figure 4. Dependence of the fundamental solution to the Cauchy problem on distance for α = 0 . 5 , t ¯ = 0 . 25 , b ¯ = 0 . 5 and various values of β .
Entropy 19 00203 g004
Figure 5. Dependence of the fundamental solution to the source problem on distance for α = 0 . 5 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Figure 5. Dependence of the fundamental solution to the source problem on distance for α = 0 . 5 , β = 2 , t ¯ = 0 . 25 and various values of b ¯ .
Entropy 19 00203 g005
Figure 6. Dependence of the fundamental solution to the source problem on distance for α = 0 . 5 , t ¯ = 0 . 25 , b ¯ = 0 . 5 and various values of β .
Figure 6. Dependence of the fundamental solution to the source problem on distance for α = 0 . 5 , t ¯ = 0 . 25 , b ¯ = 0 . 5 and various values of β .
Entropy 19 00203 g006

Share and Cite

MDPI and ACS Style

Povstenko, Y.; Kyrylych, T.; Rygał, G. Fractional Diffusion in a Solid with Mass Absorption. Entropy 2017, 19, 203. https://doi.org/10.3390/e19050203

AMA Style

Povstenko Y, Kyrylych T, Rygał G. Fractional Diffusion in a Solid with Mass Absorption. Entropy. 2017; 19(5):203. https://doi.org/10.3390/e19050203

Chicago/Turabian Style

Povstenko, Yuriy, Tamara Kyrylych, and Grażyna Rygał. 2017. "Fractional Diffusion in a Solid with Mass Absorption" Entropy 19, no. 5: 203. https://doi.org/10.3390/e19050203

APA Style

Povstenko, Y., Kyrylych, T., & Rygał, G. (2017). Fractional Diffusion in a Solid with Mass Absorption. Entropy, 19(5), 203. https://doi.org/10.3390/e19050203

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