Next Article in Journal
Spectral Entropy, Empirical Entropy and Empirical Exergy for Deterministic Boundary-Layer Structures
Next Article in Special Issue
Learning Entropy: Multiscale Measure for Incremental Learning
Previous Article in Journal / Special Issue
Time Eigenstates for Potential Functions without Extremal Points
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion

by
Yuriy Povstenko
1,2
1
Institute of Mathematics and Computer Science, Jan Długosz University in Częstochowa, Armii Krajowej 13/15, Częstochowa 42-200, Poland
2
Department of Computer Science, European University of Informatics and Economics (EWSIE) Białostocka 22, Warsaw 03-741, Poland 
Entropy 2013, 15(10), 4122-4133; https://doi.org/10.3390/e15104122
Submission received: 27 August 2013 / Revised: 22 September 2013 / Accepted: 22 September 2013 / Published: 27 September 2013
(This article belongs to the Special Issue Dynamical Systems)

1. Introduction

The standard heat conduction (diffusion) equation for temperature T
Entropy 15 04122 i006
is obtained from the balance equation for energy
Entropy 15 04122 i007
where ρ is the mass density, C is the specific heat capacity, q is the heat flux vector, and the classical Fourier law which states the linear dependence between the heat flux vector q and the temperature gradient
Entropy 15 04122 i008
with k being the thermal conductivity. In the heat conduction Equation (1) Entropy 15 04122 i009 is the heat diffusivity coefficient.
To describe heat conduction in media with complex internal structure, the standard parabolic Equation (1) is no longer accurate enough. In nonclassical theories, the Fourier law Equation (3) and the parabolic heat conduction Equation (1) are replaced by more general equations (see [1,2,3,4,5,6]). The time-nonlocal dependence between the heat flux vector q and the temperature gradient [7,8]
Entropy 15 04122 i010
results in the heat conduction with memory [7,8]
Entropy 15 04122 i011
Several particular cases of choice of the memory kernel Entropy 15 04122 i012 were analyzed in [9,10,11,12]. The time-nonlocal dependence between the heat flux vector q and the temperature gradient with the long-tail power kernel [9,10,11,12]
Entropy 15 04122 i013
Entropy 15 04122 i014
where Entropy 15 04122 i015 is the gamma function, can be interpreted in terms of fractional calculus:
Entropy 15 04122 i016
Entropy 15 04122 i017
where Entropy 15 04122 i018 and Entropy 15 04122 i019 are the Riemann–Liouville fractional integral and derivative of the order respectively [13,14,15,16]:
Entropy 15 04122 i020
Entropy 15 04122 i021
The balance Equation (2) and the constitutive Equations (8) and (9) yield the time-fractional equation
Entropy 15 04122 i022
with the Caputo fractional derivative
Entropy 15 04122 i023
The details of obtaining the time-fractional heat conduction Equation (12) from the balance Equation (2) and the constitutive Equations (8) and (9) can be found in [17].
Equations with fractional derivatives, in particular the time-fractional heat conduction equation (diffusion-wave equation), describe many important physical phenomena in different media (see [9,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], among many others). Fractional calculus plays a significant part in studies of entropy [33,34,35,36,37,38]. It should be noted that entropy is also used in analysis of anomalous diffusion processes and fractional diffusion equation [39,40,41,42,43,44,45].
Different kinds of boundary conditions for Equation (12) in a bounded domain were analyzed in [46,47]. It should be emphasized that due to the generalized constitutive equations for the heat flux (8) and (9) the boundary conditions for the time-fractional heat conduction equation have their traits in comparison with those for the standard heat conduction equation. The Dirichlet boundary condition specifies the temperature over the surface of a body
Entropy 15 04122 i024
For time-fractional heat conduction Equation (12) two types of Neumann boundary condition can be considered: the mathematical condition with the prescribed boundary value of the normal derivative of temperature
Entropy 15 04122 i025
and the physical condition with the prescribed boundary value of the heat flux
Entropy 15 04122 i026
Entropy 15 04122 i027
Here n is the outer unit normal the boundary surface. Similarly, the mathematical Robin boundary condition is a specification of a linear combination of the values of temperature and the values of its normal derivative at the boundary of the domain
Entropy 15 04122 i028
with some nonzero constants c1 and c2, while the physical Robin boundary condition specifies a linear combination of the values of temperature and the values of the heat flux at the boundary of the domain. For example, the Newton condition of convective heat exchange between a body and the environment with the temperature TE
Entropy 15 04122 i029
where h is the convective heat transfer coefficient, leads to
Entropy 15 04122 i030
Entropy 15 04122 i031
If the surfaces of two solids are in perfect thermal contact, the temperatures on the contact surface and the heat fluxes through the contact surface are the same for both solids, and the boundary conditions of the fourth kind are obtained:
Entropy 15 04122 i032
Entropy 15 04122 i033
where subscripts 1 and 2 refer to the first and second solid, respectively, and n is the common unit normal at the contact surface. In fractional calculus, where integrals and derivatives of arbitrary (not only integer) order are considered, there is no sharp boundary between integration and differentiation. For this reason, some authors [15,25] do not use a separate notation for the fractional integral Entropy 15 04122 i034 The fractional integral of the order Entropy 15 04122 i035 is denoted as Entropy 15 04122 i036. In the equation of perfect thermal contact (23) Entropy 15 04122 i037 and Entropy 15 04122 i038, are understood in this sense.
Starting from the pioneering papers [48,49,50,51,52], considerable interest has been shown in solutions to time-fractional heat conduction equation. In the literature, there are only a few papers in which the fractional heat conduction equation (fractional diffusion-wave equation) is studied in composite medium [47,53,54]. In the present paper, the problem of fractional heat conduction in a composite medium consisting of a spherical inclusion Entropy 15 04122 i039 and a matrix Entropy 15 04122 i040 being in perfect thermal contact at Entropy 15 04122 i041 is considered. The heat conduction in each region is described by the time-fractional heat conduction equation with the Caputo derivative of fractional order Entropy 15 04122 i042 and Entropy 15 04122 i043, respectively.

2. Statement of the Problem

Consider the time-fractional heat conduction equations in a spherical inclusion
Entropy 15 04122 i044
and in a matrix
Entropy 15 04122 i045
under the initial conditions
Entropy 15 04122 i046
Entropy 15 04122 i047
Entropy 15 04122 i048
Entropy 15 04122 i049
and the boundary condition of perfect thermal contact
Entropy 15 04122 i050
Entropy 15 04122 i051
The boundedness condition at the origin and the zero condition at infinity are also assumed:
Entropy 15 04122 i052
The limitations on α and β in Equations (26–29) express the fact that if Entropy 15 04122 i053 or Entropy 15 04122 i054, then the additional condition on the first time derivative should be also imposed.
In what follows we restrict ourselves to the particular case when a sphere Entropy 15 04122 i055 is at initial uniform temperature T0 and the matrix Entropy 15 04122 i056 is at initial zero temperature
Entropy 15 04122 i057
Entropy 15 04122 i058
Entropy 15 04122 i059
Entropy 15 04122 i060
The Laplace transform with respect to time t applied to Equations (24) and (25) leads to two ordinary differential equations
Entropy 15 04122 i061
Entropy 15 04122 i062
having the solutions
Entropy 15 04122 i063
Entropy 15 04122 i064
It follows from conditions at the origin and at infinity Equation (32) that
Entropy 15 04122 i065
The integration constants B1 and B2 are obtained from the perfect thermal contact boundary conditions Equations (30) and (31)
Entropy 15 04122 i066
Entropy 15 04122 i067
Hence, the solution is written as
Entropy 15 04122 i068
Entropy 15 04122 i069
Now we will investigate the approximate solution of the considered problem for small values of time. In the case of classical heat conduction this method was described in [55,56]. Based on Tauberian theorems for the Laplace transform (see, for example [57]), for small values of time t (the large values of the transform variable S) we can neglect the exponential term in comparison with 1,
Entropy 15 04122 i070
thus obtaining
Entropy 15 04122 i071
Entropy 15 04122 i072
In the following particular cases Entropy 15 04122 i073, Entropy 15 04122 i074; Entropy 15 04122 i075 the denominator in Equations (47) and (48) can be treated as a cubic equation and the decomposition into the sum of partial fractions can be obtained similar to that used in [58].
Now we will consider another particular case when Entropy 15 04122 i076
To invert the Laplace transform the following formula will be used [14,15,16]
Entropy 15 04122 i077
where Entropy 15 04122 i078 is the generalized Mittag-Leffler function in two parameters
Entropy 15 04122 i079
Additionally [51,52,59,60,61]
Entropy 15 04122 i080
Entropy 15 04122 i081
Entropy 15 04122 i082
Here Entropy 15 04122 i083 is the Wright function [1,51,52,62]
Entropy 15 04122 i084
whereas Entropy 15 04122 i085 is the Mainardi function [15,51,52]
Entropy 15 04122 i086
From Equations (47) and (48) we get:
Entropy 15 04122 i087
Entropy 15 04122 i088
where
Entropy 15 04122 i089
It should be emphasized that the solution is expressed in terms of the Mainardi function Entropy 15 04122 i090 and the Wright function Entropy 15 04122 i091. The limitation Entropy 15 04122 i092 in Equations (51–53) means that Entropy 15 04122 i093 in Equations (56) and (57).

4. Conclusions

We have obtained the approximate solution to the time-fractional heat conduction equations in a composite body consisting of a matrix and spherical inclusion with different thermophysical properties. The conditions of perfect thermal contact have been assumed: the temperatures at the boundary surface are equal and the heat fluxes through the contact surface are the same. The Laplace integral transform allows us to obtain the ordinary differential equations for temperatures. Inversion of the Laplace transform has been carried out analytically for small values of time.

Acknowledgments

The author thanks the anonymous reviewers for their helpful suggestions.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Petrov, N.; Vulchanov, N. A note on the non-classical heat condition. Bulg. Acad. Sci. Theor. Appl. Mech. 1982, 13, 35–39. [Google Scholar]
  2. Chandrasekharaiah, D.S. Thermoelasticity with second sound: a review. Appl. Mech. Rev. 1986, 39, 355–376. [Google Scholar] [CrossRef]
  3. Joseph, D.D.; Preziosi, L. Heat waves. Rev. Mod. Phys. 1989, 61, 41–73. [Google Scholar] [CrossRef]
  4. Tamma, K.; Zhou, X. Macroscale and microscale thermal transport and thermo-mechanical inter-actions: some noteworthy perspectives. J. Thermal Stresses 1998, 21, 405–449. [Google Scholar] [CrossRef]
  5. Chandrasekharaiah, D.S. Hyperbolic thermoelasticity: a review of recent literature. Appl. Mech. Rev. 1998, 51, 705–729. [Google Scholar] [CrossRef]
  6. Ignaczak, J.; Ostoja-Starzewski, M. Thermoelasticity with Finite Wave Speeds; Oxford University Press: London, UK, 2009. [Google Scholar]
  7. Nigmatullin, R.R. To the theoretical explanation of the “universal response”. Phys. Stat. Sol. (b) 1984, 123, 739–745. [Google Scholar] [CrossRef]
  8. Nigmatullin, R.R. On the theory of relaxation for systems with “remnant” memory. Phys. Stat. Sol. (b) 1984, 124, 389–393. [Google Scholar] [CrossRef]
  9. Povstenko, Y.Z. Fractional heat conduction equation and associated thermal stresses. J. Thermal Stresses 2005, 28, 83–102. [Google Scholar] [CrossRef]
  10. Povstenko, Y.Z. Thermoelasticity which uses fractional heat conduction equation. J. Math. Sci. 2009, 162, 296–305. [Google Scholar] [CrossRef]
  11. Povstenko, Y.Z. Theory of thermoelasticity based on the space-time-fractional heat conduction equation. Phys. Scr. 2009. [Google Scholar] [CrossRef]
  12. Povstenko, Y.Z. Fractional Cattaneo-type equations and generalized thermoelasticity. J. Thermal Stresses 2011, 34, 97–114. [Google Scholar] [CrossRef]
  13. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publisher: New York, NY, USA, 1993. [Google Scholar]
  14. 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-Verlag: New York, NY, USA, 1997; pp. 223–276. [Google Scholar]
  15. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  16. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  17. 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]
  18. 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]
  19. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A: Math. Gen. 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  20. Zaslavsky, G.M. Chaos, fractional kinetics, and anomalous transport. Phys. Rep. 2002, 371, 461–580. [Google Scholar] [CrossRef]
  21. Rabotnov, Yu.N. Creep Problems in Structural Members; North-Holland Publishing Company: Amsterdam, The Netherlands, 1969. [Google Scholar]
  22. Mainardi, F. Applications of fractional calculus in mechanics. In Transform Methods and Special Functions; Rusev, P., Dimovski, I., Kiryakova, V., Eds.; Bulgarian Academy of Sciences: Sofia, Bulgaria, 1998; pp. 309–334. [Google Scholar]
  23. Rossikhin, Yu.A.; Shitikova, M.V. Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids. Appl. Mech. Rev. 1997, 50, 15–67. [Google Scholar] [CrossRef]
  24. West, B.J.; Bologna, M.; Grigolini, P. Physics of Fractals Operators; Springer-Verlag: New York, NY, USA, 2003. [Google Scholar]
  25. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers, Inc.: Redding, CA, USA, 2006. [Google Scholar]
  26. Sabatier, J.; Agrawal, O.P.; Tenreiro Machado, J.A. (Eds.) Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering; Springer-Verlag: Dordrecht, The Netherlands, 2007.
  27. Gafiychuk, V.; Datsko, B. Mathematical modeling of different types of instabilities in time fractional reaction-diffusion systems. Comput. Math. Appl. 2010, 59, 1101–1107. [Google Scholar] [CrossRef]
  28. Baleanu, D.; Güvenç, Z.B.; Tenreiro Machado, J.A. (Eds.) New Trends in Nanotechnology and Fractional Calculus Applications; Springer-Verlag: New York, NY, USA, 2010.
  29. Rossikhin, Y.A.; Shitikova, M.V. Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results. Appl. Mech. Rev. 2010, 63, 010801. [Google Scholar] [CrossRef]
  30. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  31. Datsko, B.; Gafiychuk, V. Complex nonlinear dynamics in subdiffusive activator–inhibitor systems. Commun. Nonlinear Sci. Numer. Simulat. 2012, 17, 1673–1680. [Google Scholar] [CrossRef]
  32. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers; Springer-Verlag: Berlin, Germany, 2013. [Google Scholar]
  33. Zunino, L.; Pérez, D.G.; Martín, M.T.; Garavaglia, M.; Plastino, A.; Rosso, O.A. Permutation entropy of fractional Brownian motion and fractional Gaussian noise. Phys. Lett. A 2008, 372, 4768–4774. [Google Scholar] [CrossRef]
  34. Ubriaco, M.R. Entropies based on fractional calculus. Phys. Lett. A 2009, 373, 2516–2519. [Google Scholar] [CrossRef]
  35. Tenreiro Machado, J.A. Entropy analysis of integer and fractional dynamical systems. Nonlinear Dyn. 2010, 62, 371–378. [Google Scholar] [CrossRef]
  36. Tenreiro Machado, J.A. Fractional dynamics of a system with particles subjected to impacts. Commun. Nonlinear Sci. Numer. Simulat. 2011, 16, 4596–4601. [Google Scholar] [CrossRef]
  37. Jumarie, G. Path probability of random fractional systems defined by white noises in coarse-grained time applications of fractional entropy. Fract. Diff. Calc. 2011, 1, 45–87. [Google Scholar] [CrossRef]
  38. Tenreiro Machado, J.A. Shannon information and power law analysis of the chromosome code. Abstr. Appl. Anal. 2012. [Google Scholar] [CrossRef]
  39. Essex, C.; Schulzky, C.; Franz, A.; Hoffmann, K.H. Tsallis and Rényi entropies in fractional diffusion and entropy production. Physica A 2000, 284, 299–308. [Google Scholar] [CrossRef]
  40. Cifani, S.; Jakobsen, E.R. Entropy solution theory for fractional degenerate convection–diffusion equations. Ann. Inst. Henri Poincare (C) Nonlinear Anal. 2011, 28, 413–441. [Google Scholar] [CrossRef] [Green Version]
  41. Magin, R.; Ingo, C. Entropy and information in a fractional order model of anomalous diffusion. In Proceedings of the 16th IFAC Symposium on System Identification, Brussels, Belgium, 11–13 July 2011; Kinnaert, M., Ed.; International Federation of Automatic Control: Brussels, Belgium, 2012; pp. 428–433. [Google Scholar]
  42. Magin, R.; Ingo, C. Spectral entropy in a fractional order model of anomalous diffusion. In Proceedings of the 13th International Carpathian Control Conference, High Tatras, Slovakia, 28–31 May 2012; Petraš, I., Podlubny, I., Kostúr, J., Kačur, J., Mojžišová, A, Eds.; Institute of Electrical and Electronics Engineers: Košice, Slovakia, 2012; pp. 458–463. [Google Scholar]
  43. Prehl, J.; Essex, C.; Hoffmann, K.H. Tsallis relative entropy and anomalous diffusion. Entropy 2012, 14, 701–706. [Google Scholar] [CrossRef]
  44. Prehl, J.; Boldt, F.; Essex, C.; Hoffmann, K.H. Time evolution of relative entropies for anomalous diffusion. Entropy 2013, 15, 2989–3006. [Google Scholar] [CrossRef]
  45. Magin, R.L.; Ingo, C.; Colon-Perez, L.; Triplett, W.; Mareci, T.H. Characterization of anomalous diffusion in porous biological tissues using fractional order derivatives and entropy. Microporous Mesoporous Mater. 2013, 178, 39–43. [Google Scholar] [CrossRef] [PubMed]
  46. Povstenko, Y. Different kinds of boundary conditions for time-fractional heat conduction equation. In Proceedings of the 13th International Carpathian Control Conference, High Tatras, Slovakia, 28–31 May 2012; Petraš, I., Podlubny, I., Kostúr, J., Kačur, J., Mojžišová, A, Eds.; Institute of Electrical and Electronics Engineers: Košice, Slovakia, 2012; pp. 588–591. [Google Scholar]
  47. Povstenko, Y.Z. Fractional heat conduction in infinite one-dimensional composite medium. J. Thermal Stresses 2013, 36, 351–363. [Google Scholar] [CrossRef]
  48. Wyss, W. The fractional diffusion equation. J. Math. Phys. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  49. Schneider, W.R.; Wyss, W. Fractional diffusion and wave equations. J. Math. Phys. 1989, 30, 134–144. [Google Scholar] [CrossRef]
  50. Fujita, Y. Integrodifferential equation which interpolates the heat equation and the wave equation. Osaka J. Math. 1990, 27, 309–321. [Google Scholar]
  51. Mainardi, F. The fundamental solutions for the fractional diffusion-wave equation. Appl. Math. Lett. 1996, 9, 23–28. [Google Scholar] [CrossRef]
  52. Mainardi, F. Fractional relaxation-oscillation and fractional diffusion-wave phenomena. Chaos, Solitons Fractals 1996, 7, 1461–1477. [Google Scholar] [CrossRef]
  53. Chen, S.; Jiang, X. Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus. Physica A 2012, 392, 3865–3874. [Google Scholar] [CrossRef]
  54. Povstenko, Y. Fundamental solutions to time-fractional heat conduction equations in two joint half-lines. Cent. Eur. J. Phys. 2013. [Google Scholar] [CrossRef]
  55. Luikov, A.V. Analytical Heat Diffusion Theory; Academic Press: New York, NY, USA, 1968. [Google Scholar]
  56. Özişik, M.N. Heat Conduction; John Wiley: New York, NY, USA, 1980. [Google Scholar]
  57. Debnath, L.; Bhatta, D. Integral Transforms and Their Applications; Chapman & Hall/CRC: Boca Raton, FL, USA, 2007. [Google Scholar]
  58. Povstenko, Y. Time-fractional heat conduction in an infinite medium with a spherical hole under Robin boundary condition. Fract. Calc. Appl. Anal. 2013, 16, 354–369. [Google Scholar] [CrossRef]
  59. Mikusiński, J. On the function whose Laplace transform is Entropy 15 04122 i094. Stud. Math. 1959, 18, 191–198. [Google Scholar]
  60. Stanković, B. On the function of E.M. Wright. Publ. Inst. Math. 1970, 10, 113–124. [Google Scholar]
  61. Gajić, Lj.; Stanković, B. Some properties of Wright’s function. Publ. Inst. Math. 1976, 20, 91–98. [Google Scholar]
  62. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; Volume 3, McGraw-Hill: New York, NY, USA, 1955. [Google Scholar]

Share and Cite

MDPI and ACS Style

Povstenko, Y. Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion. Entropy 2013, 15, 4122-4133. https://doi.org/10.3390/e15104122

AMA Style

Povstenko Y. Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion. Entropy. 2013; 15(10):4122-4133. https://doi.org/10.3390/e15104122

Chicago/Turabian Style

Povstenko, Yuriy. 2013. "Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion" Entropy 15, no. 10: 4122-4133. https://doi.org/10.3390/e15104122

APA Style

Povstenko, Y. (2013). Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion. Entropy, 15(10), 4122-4133. https://doi.org/10.3390/e15104122

Article Metrics

Back to TopTop