Next Article in Journal
Integral Involving Bessel Functions Arising in Propagation Phenomena
Next Article in Special Issue
Fractional Calculus as a Simple Tool for Modeling and Analysis of Long Memory Process in Industry
Previous Article in Journal
Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect
Previous Article in Special Issue
Adaptive Pinning Synchronization of Fractional Complex Networks with Impulses and Reaction–Diffusion Terms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact

1
Faculty of Mathematics and Applied Physics, Rzeszow University of Technology, Powstancow Warszawy 8, 35-959 Rzeszow, Poland
2
Institute for Applied Problems of Mechanics and Mathematics NAS of Ukraine, 79060 Lviv, Ukraine
3
BERG Faculty, Technical University of Kosice, B. Nemcovej 3, 04200 Kosice, Slovakia
4
Faculty of Mathematical and Natural Sciences, Jan Dlugosz University in Czestochowa, Armii Krajowej 13/15, 42-200 Czestochowa, Poland
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(5), 433; https://doi.org/10.3390/math7050433
Submission received: 24 April 2019 / Revised: 10 May 2019 / Accepted: 12 May 2019 / Published: 16 May 2019
(This article belongs to the Special Issue Fractional Order Systems)

Abstract

:
The time-fractional diffusion equation with mass absorption in a sphere is considered under harmonic impact on the surface of a sphere. The Caputo time-fractional derivative is used. The Laplace transform with respect to time and the finite sin-Fourier transform with respect to the spatial coordinate are employed. A graphical representation of the obtained analytical solution for different sets of the parameters including the order of fractional derivative is given.

1. Introduction

The classical parabolic diffusion equation with heat or mass absorption [1]
u t = a Δ u b u
also describes bioheat transfer, lateral surface mass or heat exchange in a thin plate, heating of tissue during laser treatment irradiation, etc. (see, for example [2,3,4,5]). The Klein-Gordon equation
2 u t 2 = a Δ u b u
is used in solid state physics, classical mechanics, nonlinear optics, and quantum field theory [6,7].
The time-fractional equation
α u t α = a Δ u b u , 0 < α 2 ,
can be considered as the extension of the parabolic Equation (1) and hyperbolic Equation (2) and was studied in several publications [8,9,10,11,12,13,14].
It should be noted that such a generalization of many classical differential equations with integer derivatives has numerous applications in rheology, geology, physics, plasma physics, chemistry, geophysics, engineering, biology, bio-engineering, finance, and medicine (see [15,16,17,18,19,20,21,22,23,24,25,26,27,28,29], among many others). There is a great variety of inhomogeneous media where transport phenomena exhibit anomalous properties, the investigation of which is essential to refine understanding the basic characteristics of complex systems widely met in nature. Therefore, studying fractional equations has generated increasing attention of scientists in many disciplines. At present, the fractional diffusion-wave equation is generally used to describe a large class of systems at different scales (from the molecular [30] to the space one [31]) which cover media of the diverse nature (from plasma physics [29] to living tissue [3]). The study of this equation is also of interest from the point of view of understanding the complex spatio-temporal dynamics in nonlinear systems of fractional order [32,33].
In Equation (3) and further in this paper, for more concise notation, d α f ( t ) d t α denotes the Caputo fractional derivative [16,34]
d α f ( t ) d t α = 1 Γ ( n α ) 0 t ( t τ ) n α 1 d n f ( τ ) d τ n d τ , n 1 < α < n ,
and Γ ( α ) denotes the gamma function.
Ångström was the first to investigate the standard parabolic heat conduction equation under harmonic impact and laid the foundations for the new area of study known as “oscillatory diffusion” or “diffusion-waves” (see [35,36,37] and references therein). Periodic solutions of the bioheat equation were investigated in [38]. The harmonic point source in the bioheat equation was used in therapeutic hypotermia [39,40]; applications of the time-harmonic impact in ultrasound surgery were studied in [41].
As a rule, in the previous studies of diffusion or heat conduction equation the quasi-steady-state oscillations were investigated when the solution u ( x , t ) was represented as a product of a function of the spatial coordinates U ( x ) and the time-harmonic term e i ω t with the angular frequency ω
u ( x , t ) = U ( x ) e i ω t
without consideration of the initial conditions.
The use of assumption (5) is based on the well known formula for the derivative of the integer order n of the exponential function
d n e λ t d t n = λ n e λ t .
In the event of the non-integer order of time derivative, the assumption (5) cannot be used since [42]
d α e λ t d t α = λ α e λ t γ ( n α ) , λ t ) Γ ( n α ) λ α e λ t , n 1 < α < n ,
with γ ( a , x ) being the incomplete gamma function [43]
γ ( a , x ) = 0 x e u u a 1 d u .
It is worthy of notice that for the Riemann-Liouville fractional derivative [16,34] with the lower limit of integration at 0
D R L α f ( t ) = d n d t n 1 Γ ( n α ) 0 t ( t τ ) n α 1 f ( τ ) d τ , n 1 < α < n ,
we also have [16]
D R L α e λ t = t α E 1 , 1 α λ t λ α e λ t .
Here E α , β ( z ) is the Mittag-Leffler function in two parameters α and β [16,34]
E α , β z = n = 0 z n Γ ( α n + β ) , α > 0 , β C , z C .
In this paper, the initial-boundary-value problem for Equation (3) is studied in a spherical domain for the case of central symmetry under the Dirichlet boundary condition varying harmonically in time. The present paper develops and extends the results of the previous investigations [44,45], where the corresponding problems for line and half-line domains were investigated.

2. Statement of the Problem

The time-fractional diffusion equation with mass absorption (mass release) is examined in a sphere
α u t α = a 2 u r 2 + 2 r u r b u , 0 < r < R , 0 < t < , 0 < α 2 ,
under zero initial conditions
t = 0 : u = 0 , 0 < α 2 ,
t = 0 : u t = 0 , 1 < α 2 ,
and harmonic impact on the surface of a sphere
r = R : u = u 0 e i ω t .
As in the case of classical diffusion equation (when α = 1 ) and the wave equation (when α = 2 ) the boundedness condition at the origin is also adopted:
r = 0 : u .
In what follows, the integral transform technique will be used. Recall the Laplace transform rule for the Caputo derivative
L d α f ( t ) d t α = s α f * ( s ) k = 0 n 1 f ( k ) ( 0 + ) s α 1 k , n 1 < α < n ,
where the transform is marked by the asterisk, and s is the Laplace transform variable.
The following finite sin-Fourier transform is amenable to the central symmetric problem in a spherical domain 0 r R [46]. For the Dirichlet boundary condition:
F { f ( r ) } = f ˜ ( ξ k ) = 0 R r f ( r ) sin ( r ξ k ) ξ k d r ,
F 1 { f ˜ ( ξ k ) } = f ( r ) = 2 R k = 1 ξ k f ˜ ( ξ k ) sin ( r ξ k ) r ,
where the transform is marked by the tilde, and
ξ k = k π R
is the Fourier transform variable.
For the central symmetric Laplace operator
F d 2 f ( r ) d r 2 + 2 r d f ( r ) d r = ξ k 2 f ˜ ( ξ k ) + ( 1 ) k + 1 R f ( R ) .
Applying to the problems (12)–(16) the Laplace transform with respect to time t and the finite sin-Fourier transform (18) with respect to the radial coordinate r, we get in the transform domain
u ˜ * ( ξ k , s ) = ( 1 ) k + 1 a R u 0 1 s α + a ξ k 2 + b 1 s i ω .
The solution is obtained after inversion of the integral transforms:
u ( r , t ) = 2 a u 0 r k = 1 ( 1 ) k + 1 ξ k sin r ξ k 0 t τ α 1 E α , α a ξ k 2 + b τ α e i ω ( t τ ) d τ ,
where E α , β ( z ) is the Mittag-Leffler function (11), and the convolution theorem as well as the following equation for the inverse Laplace transform [16]
L 1 s α β s α + γ = t β 1 E α , β γ t α , α > 0 , β > 0 ,
have been used.
In numerical calculations, the nondimensional quantities are used:
u ¯ = u u 0 , r ¯ = r R , t ¯ = a 1 / α R 2 / α t , ω ¯ = R 2 / α a 1 / α ω , b ¯ = R 2 a b , η k = R ξ k = k π .
Hence, using in integral in (23) the substitution τ = t w , for the real part of the solution we get
u ¯ ( r ¯ , t ¯ , b ¯ , ω ¯ ) = 2 t ¯ α r ¯ k = 1 ( 1 ) k + 1 η k sin r ¯ η k 0 1 w α 1 E α , α η k 2 + b ¯ t ¯ α w α cos ω ¯ t ¯ ( 1 w ) d w ,
where we can see that the solution u ¯ depends not only on time and spatial coordinate, but also on the parameters b ¯ and ω ¯ .
To simplify calculations, it would be worthwhile to introduce a substitution z = w α . Hence,
u ¯ ( r ¯ , t ¯ , b ¯ , ω ¯ ) = 2 t ¯ α r ¯ k = 1 ( 1 ) k + 1 η k sin r ¯ η k 0 1 E α , α η k 2 + b ¯ t ¯ α z cos ω ¯ t ¯ 1 z 1 / α d z .
To evaluate the Mittag-Leffler function the algorithms suggested in [47] were used; see also the MATLAB function [48] that implements these algorithms.
The numerical results are shown in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. Calculations were carried out for the grid size step Δ t ¯ = 0.1 , Δ r ¯ = 0.05 . A graphical representation of the solution (27) makes it possible to analyze not only the limiting cases of the problem (see Section 3), but also to understand the influence of the main parameters of the problem (including the value of the order of fractional derivative) on the spatial-temporal evolution of the solution.
It is seen from Figures that time oscillations of the solution are governed by the harmonic term e i ω t . The increasing of absorption parameter b decreases the oscillation amplitude (Figure 2b and Figure 3d), whereas increasing α increases it (Figure 2a–d and Figure 3d). As the Mittag-Leffler function
E 2 , 2 x 2 = sin x x ,
the space oscillations of the solution depending on the order of fractional derivative appear for α 1.5 (Figure 1d and Figure 2d) and become well-marked for α approaching 2. The influence of both factors is evident from Figure 4 and Figure 5.

3. Analysis of the Quasi-Steady-State Oscillations

Now, we shall investigate two particular cases of the problem studied in the previous section corresponding to the integer values of the order of time derivative. For α = 1 , we have
u ˜ * ( ξ k , s ) = ( 1 ) k + 1 a R u 0 1 s + a ξ k 2 + b 1 s i ω .
Taking into account that [49]
L 1 1 s + p s + q = e q t e p t p q ,
we arrive at the solution to the bioheat equation
u ( r , t ) = 2 a u 0 r k = 1 1 k + 1 ξ k sin r ξ k a ξ k 2 + b + i ω e i ω t 2 a u 0 r k = 1 1 k + 1 ξ k sin r ξ k a ξ k 2 + b + i ω e a ξ k 2 + b t .
Similarly, for α = 2 ,
u ˜ * ( ξ k , s ) = ( 1 ) k + 1 a R u 0 1 s 2 + a ξ k 2 + b 1 s i ω .
Taking into consideration that [49]
L 1 1 s 2 + p 2 s + q = 1 p 2 + q 2 e q t cos p t + q p sin p t ,
we obtain the solution to the Klein-Gordon equation
u ( r , t ) = 2 a u 0 r k = 1 1 k + 1 ξ k sin r ξ k a ξ k 2 + b ω 2 e i ω t 2 a u 0 r k = 1 1 k + 1 ξ k sin r ξ k a ξ k 2 + b ω 2 cos a ξ k 2 + b t + i ω a ξ k 2 + b sin a ξ k 2 + b t .
For integer α , we can assume that
u ( r , t ) = U ( r ) e i ω t .
For α = 1 , the function U ( r ) fulfills the equation
d 2 U d r 2 + 2 r d U d r b + i ω a U = 0 ,
under the boundary condition
r = R : U ( r ) = u 0 ,
and for b > 0 has the solution bounded at the origin
U ( r ) = R u 0 r sinh r ( b + i ω ) / a sinh R ( b + i ω ) / a .
Therefore,
u ( r , t ) = R u 0 r sinh r ( b + i ω ) / a sinh R ( b + i ω ) / a e i ω t
(for negative value of b, sinh will be substituted by sin).
The first term in the solution (31) can be evaluated analytically using the following formula [50]
k = 1 1 k + 1 k k 2 + p 2 sin k r = π 2 sinh ( r p ) sinh ( π p ) , π < r < π , .
k = 1 1 k + 1 k k 2 p 2 sin k r = π 2 sin ( r p ) sin ( π p ) , π < r < π .
It was emphasized in [51] that Equation (41) is also valid for complex values of p and hence turns into Equation (40) for imaginary p.
Taking into account Equation (40), we obtain that the first term in the solution (31) coincides with the quasi-steady-state solution (39), whereas the second term in Equation (31) describes the transient process.
The similar analysis can be carried out for α = 2 based on the assumption (35). In this case, the function U ( r ) fulfills the equation
d 2 U d r 2 + 2 r d U d r b ω 2 a U = 0 ,
under the boundary condition
r = R : U ( r ) = u 0 ,
and for b > ω 2 has the solution bounded at the origin
U ( r ) = R u 0 r sinh r ( b ω 2 ) / a sinh R ( b ω 2 ) / a ,
whereas for b < ω 2
U ( r ) = R u 0 r sin r ( ω 2 b ) / a sin R ( ω 2 b ) / a .
Hence, for b > ω 2
u ( r , t ) = R u 0 r sinh r ( b ω 2 ) / a sinh R ( b ω 2 ) / a e i ω t .
and for b < ω 2
u ( r , t ) = R u 0 r sin r ( ω 2 b ) / a sin R ( ω 2 b ) / a e i ω t .
The first term in the solution (34) after accounting for Equations (40) and (41) coincides with the quasi-steady-state solutions (46) and (47), respectively, whereas the second term in (34) describes the transient process.

4. Conclusions

The time-fractional diffusion-wave equation with the Caputo fractional derivative of the order 0 < α 2 with mass absorption was studied in a spherical domain under the Dirichlet boundary condition varying harmonically in time. The Caputo derivative of the exponential function has a much more complicated form than the corresponding derivative of the integer order. Hence, the assumption that the solution of the problem can be represented as a product of a function of the spatial coordinate and the time-harmonic term without consideration of the initial conditions cannot be used. The solution is obtained using the Laplace transform with respect to time and the finite sin-Fourier transform specifically adapted for a spherical domain and is expressed in terms of the Mittag-Leffler function. A graphical representation of the obtained analytical solution demonstrates the influence of the main parameters of the problem including the value of the order of fractional derivative on the spatial-temporal evolution of the solution.

Author Contributions

All authors have equally contributed to this work. All authors read and approved the final manuscript.

Funding

See acknowledgements of support below. The APC was funded by APVV-14-0892.

Acknowledgments

The first author is thankful for support from Rzeszow University of Technology (grant DS.FD.18.001). The second author acknowledges support provided by grants APVV-14-0892, APVV-18-0526, VEGA 1/0365/19, SK-SRB-18-0011, SK-AT-2017-0015, ARO W911NF-15-1-0228, COST CA15225. The third author would like to acknowledge the support of Jan Dlugosz University in Czestochowa (grant DS/WMP/6011/2018).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Crank, J. The Mathematics of Diffusion, 2nd ed.; Clarendon Press: Oxford, UK, 1975. [Google Scholar]
  2. Pennes, H.H. Analysis of tissue and arterial blood temperatures in the resting human forearm. J. Appl. Physiol. 1948, 1, 93–122. [Google Scholar] [CrossRef] [PubMed]
  3. Gafiychuk, V.V.; Lubashevsky, I.A.; Datsko, B.Y. Fast heat propagation in living tissue caused by branching artery network. Phys. Rev. E 2005, 72, 051920. [Google Scholar] [CrossRef]
  4. Datsko, B.Y.; Gafiychuk, V.V.; Lubashevsky, I.A.; Priezzhev, A.V. Self-localization of laser-induced tumor coagulation limited by heat diffusion through active tissue. J. Med. Eng. Technol. 2006, 30, 390–396. [Google Scholar] [CrossRef]
  5. Polyanin, A.D. Handbook of Linear Partial Differential Equations for Engineers and Scientists; Chapman & Hall/CRC: Boca Raton, FL, USA, 2002. [Google Scholar]
  6. Gravel, P.; Gauthier, C. Classical applications of the Klein-Gordon equation. Am. J. Phys. 2011, 79, 447–453. [Google Scholar] [CrossRef]
  7. Wazwaz, A.-M. Partial Differential Equations and Solitary Waves Theory; Higher Education Press: Beijing, China; Springer: Berlin, Germany, 2009. [Google Scholar]
  8. Abuteen, A.; Freihat, A.; Al-Smadi, M.; Khalil, H.; Khan, R.A. Approximate series solution of nonlinear, fractional Klein-Gordon equations using fractional reduced differential transform method. J. Math. Stat. 2016, 12, 23–33. [Google Scholar] [CrossRef]
  9. Damor, R.S.; Kumar, S.; Shukla, A.K. Solution of fractional bioheat equation in terms of Fox’s H-Function. SpringerPlus 2016, 5, 1–10. [Google Scholar] [CrossRef]
  10. Ferrás, L.L.; Ford, N.J.; Morgado, M.L.; Nóbrega, J.M.; Rebelo, M.S. Fractional Pennes’ bioheat equation: theoretical and numerical studies. Fract. Calc. Appl. Anal. 2015, 18, 1080–1106. [Google Scholar] [CrossRef]
  11. Golmankhaneh, A.K.; Golmankhaneh, A.K.; Baleanu, D. On nolinear fractional Klein-Gordon equation. Signal Process. 2011, 91, 446–451. [Google Scholar] [CrossRef]
  12. Kheiri, H.; Shahi, S.; Mojaver, A. Analytical solutions for the fractional Klein-Gordon equation. Comput. Meth. Diff. Equ. 2014, 2, 99–114. [Google Scholar]
  13. Qin, Y.; Wu, K. Numerical solution of fractional bioheat equation by quadratic spline collocation method. J. Nonlinear Sci. Appl. 2016, 9, 5061–5072. [Google Scholar] [CrossRef]
  14. Vitali, S.; Castellani, G.; Mainardi, F. Time fractional cable equation and applications in neurophysiology. Chaos Solitons Fractals 2017, 102, 467–472. [Google Scholar] [CrossRef] [Green Version]
  15. Gorenflo, R.; Mainardi, F. Fractional calculus: Integral and differential equations of fractional order. In Fractals and Fractional Calculus in Continuum Mechanics; Springer: Wien, Austria, 1997; pp. 223–276. [Google Scholar]
  16. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  17. Povstenko, Y. Fractional heat conduction equation and associated thermal stresses. J. Therm. Stress. 2005, 28, 83–102. [Google Scholar] [CrossRef]
  18. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers, Inc.: Redding, CA, USA, 2006. [Google Scholar]
  19. Gafiychuk, V.V.; Datsko, B.Y. Spatiotemporal pattern formation in fractional reaction-diffusion systems with indices of different order. Phys. Rev. E 2008, 77, 066210. [Google Scholar] [CrossRef]
  20. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  21. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  22. Datsko, B.; Luchko, Y.; Gafiychuk, V. Pattern formation in fractional reaction-diffusion systems with multiple homogeneous states. Int. J. Bifurcat. Chaos 2012, 22, 1250087. [Google Scholar] [CrossRef]
  23. Datsko, B.; Gafiychuk, V. Complex nonlinear dynamics in subdiffusive activator-inhibitor systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1673–1680. [Google Scholar] [CrossRef]
  24. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin, Germany, 2013. [Google Scholar]
  25. 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]
  26. Herrmann, R. Fractional Calculus: An Introduction for Physicists, 2nd ed.; World Scientific: Singapore, 2014. [Google Scholar]
  27. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  28. Datsko, B.; Gafiychuk, V.; Podlubny, I. Solitary travelling auto-waves in fractional reaction–diffusion systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 23, 378–387. [Google Scholar] [CrossRef] [Green Version]
  29. Anderson, J.; Moradi, S.; Rafiq, T. Non-linear Langevin and fractional Fokker-Planck equations for anomalous diffusion by Lévy stable processes. Entropy 2018, 20, 760. [Google Scholar] [CrossRef]
  30. Weiss, M.; Nilsson, T. In a mirror dimly: Tracing the movements of molecules in living cells. Trends Cell Biol. 2004, 14, 267–273. [Google Scholar] [CrossRef]
  31. Zelenyi, L.M.; Milovanov, A.V. Fractal topology and strange kinetics: From percolation theory to problems in cosmic electrodynamics. Phys. Uspekhi 2004, 47, 809–852. [Google Scholar] [CrossRef]
  32. Gafiychuk, V.; Datsko, B. Different types of instabilities and complex dynamics in reaction-diffusion systems with fractional derivatives. J. Comp. Nonlinear Dyn. 2012, 7, 031001. [Google Scholar] [CrossRef]
  33. Datsko, B.; Gafiychuk, V. Complex spatio-temporal solutions in fractional reaction-diffusion systems near a bifurcation point. Fract. Calc. Appl. Anal. 2018, 21, 237–253. [Google Scholar] [CrossRef]
  34. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  35. Mandelis, A. Diffusion waves and their uses. Phys. Today 2000, 53, 29–33. [Google Scholar] [CrossRef]
  36. Mandelis, A. Diffusion-Wave Fields: Mathematical Methods and Green Functions; Springer: New York, NY, USA, 2001. [Google Scholar]
  37. Vrentas, J.S.; Vrentas, C.M. Diffusion and Mass Transfer; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  38. Lakhssassi, A.; Kengne, E.; Semmaoui, H. Modifed Pennes’ equation modelling bio-heat transfer in living tissues: analytical and numerical analysis. Natl. Sci. 2010, 2, 1375–1385. [Google Scholar] [CrossRef]
  39. Kengne, E.; Lakhssassi, A.; Vaillancourt, R. Temperature distributions for regional hypothermia based on nonlinear bioheat equation of Pennes type: Dermis and subcutaneous tissues. Appl. Math. 2012, 3, 217–224. [Google Scholar] [CrossRef]
  40. Fasano, A.; Sequeira, A. Hemomath. The Mathematics of Blood; Springer: Cham, Switzerland, 2017. [Google Scholar]
  41. Malinen, M.; Huttunen, T.; Kaipio, J.P. Thermal dose optimization method for ultrasound surgery. Phys. Med. Biol. 2003, 48, 745–762. [Google Scholar] [CrossRef]
  42. Povstenko, Y. Fractional heat conduction in a space with a source varying harmonically in time and associated thermal stresses. J. Therm. Stress. 2016, 39, 1442–1450. [Google Scholar] [CrossRef]
  43. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables; Dover: New York, NY, USA, 1972. [Google Scholar]
  44. Povstenko, Y.; Kyrylych, T. Time-fractional diffusion with mass absorption under harmonic impact. Fract. Calc. Appl. Anal. 2018, 21, 118–133. [Google Scholar] [CrossRef]
  45. Povstenko, Y.; Kyrylych, T. Time-fractional diffusion with mass absorption in a half-line domain due to boundary value of concentration varying harmonically in time. Entropy 2018, 19, 346. [Google Scholar] [CrossRef]
  46. Povstenko, Y. Linear Fractional Diffusion-Wave Equation for Scientists and Engineers; Birkhäuser: New York, NY, USA, 2015. [Google Scholar]
  47. 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]
  48. Podlubny, I. Mittag-Leffler Function; Calculates the Mittag-Leffler Function with Desired Accuracy, MATLAB Central File Exchange, File ID 8738. Available online: www.mathworks.com/matlabcentral/fileexchange/8738 (accessed on 17 April 2019).
  49. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Tables of Integral Transforms; McGraw-Hill: New York, NY, USA, 1954; Volume 1. [Google Scholar]
  50. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Volume 1: Elementary Functions; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1986. [Google Scholar]
  51. Magnus, W.; Oberhettinger, F. Formeln und Sätze für die Speziellen Funkttionen der Mathematischen Physik, 2nd ed.; Springer: Berlin, Germany, 1948. [Google Scholar]
Figure 1. Evolution of the solution for the problems (12)–(15) under constant impact. The results of computer simulation of the formula (27) for the parameters b ¯ = 4 , ω ¯ = 0 and different values of α : α = 0.5 —(a); α = 0.75 —(b); α = 1.25 —(c); α = 1.5 —(d).
Figure 1. Evolution of the solution for the problems (12)–(15) under constant impact. The results of computer simulation of the formula (27) for the parameters b ¯ = 4 , ω ¯ = 0 and different values of α : α = 0.5 —(a); α = 0.75 —(b); α = 1.25 —(c); α = 1.5 —(d).
Mathematics 07 00433 g001
Figure 2. Evolution of the solution for the problems (12)–(15) under harmonic impact. The results of computer simulation of the formula (27) for the parameters b ¯ = 4 , ω ¯ = 4 and different values of α : α = 0.5 —(a); α = 0.75 —(b); α = 1.25 —(c); α = 1.5 —(d).
Figure 2. Evolution of the solution for the problems (12)–(15) under harmonic impact. The results of computer simulation of the formula (27) for the parameters b ¯ = 4 , ω ¯ = 4 and different values of α : α = 0.5 —(a); α = 0.75 —(b); α = 1.25 —(c); α = 1.5 —(d).
Mathematics 07 00433 g002
Figure 3. Evolution of the solution for the problems (12)–(15) with increasing the frequency ω ¯ in the sub-diffusive case. The results of computer simulation of the formula (27) for the parameters b ¯ = 1 , α = 0.75 and different values of ω ¯ : ω ¯ = 0.0 —(a); ω ¯ = 1.0 —(b); ω ¯ = 2.0 —(c); ω ¯ = 4.0 —(d).
Figure 3. Evolution of the solution for the problems (12)–(15) with increasing the frequency ω ¯ in the sub-diffusive case. The results of computer simulation of the formula (27) for the parameters b ¯ = 1 , α = 0.75 and different values of ω ¯ : ω ¯ = 0.0 —(a); ω ¯ = 1.0 —(b); ω ¯ = 2.0 —(c); ω ¯ = 4.0 —(d).
Mathematics 07 00433 g003
Figure 4. Evolution of the solution for the problems (12)–(15) with increasing the frequency ω ¯ in the sub-wave case. The results of computer simulation of the formula (27) for the parameters b ¯ = 1 , α = 1.75 and different values of ω ¯ : ω ¯ = 0.0 —(a); ω ¯ = 1.0 —(b); ω ¯ = 2.0 —(c); ω ¯ = 4.0 —(d).
Figure 4. Evolution of the solution for the problems (12)–(15) with increasing the frequency ω ¯ in the sub-wave case. The results of computer simulation of the formula (27) for the parameters b ¯ = 1 , α = 1.75 and different values of ω ¯ : ω ¯ = 0.0 —(a); ω ¯ = 1.0 —(b); ω ¯ = 2.0 —(c); ω ¯ = 4.0 —(d).
Mathematics 07 00433 g004
Figure 5. Evolution of the solution for the problems (12)–(15) for different orders of fractional derivative. The results of computer simulation of the formula (27) for the parameters b ¯ = 2 , ω ¯ = 2.0 and different values of α : α = 0.75 —(a); α = 0.95 —(b); α = 1.75 —(c); α = 1.95 —(d).
Figure 5. Evolution of the solution for the problems (12)–(15) for different orders of fractional derivative. The results of computer simulation of the formula (27) for the parameters b ¯ = 2 , ω ¯ = 2.0 and different values of α : α = 0.75 —(a); α = 0.95 —(b); α = 1.75 —(c); α = 1.95 —(d).
Mathematics 07 00433 g005

Share and Cite

MDPI and ACS Style

Datsko, B.; Podlubny, I.; Povstenko, Y. Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact. Mathematics 2019, 7, 433. https://doi.org/10.3390/math7050433

AMA Style

Datsko B, Podlubny I, Povstenko Y. Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact. Mathematics. 2019; 7(5):433. https://doi.org/10.3390/math7050433

Chicago/Turabian Style

Datsko, Bohdan, Igor Podlubny, and Yuriy Povstenko. 2019. "Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact" Mathematics 7, no. 5: 433. https://doi.org/10.3390/math7050433

APA Style

Datsko, B., Podlubny, I., & Povstenko, Y. (2019). Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact. Mathematics, 7(5), 433. https://doi.org/10.3390/math7050433

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