Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion

The problem of fractional heat conduction in a composite medium consisting of a spherical inclusion (0< r < R) and a matrix (R < r < ∞) being in perfect thermal contact at r = R is considered. The heat conduction in each region is described by the time-fractional heat conduction equation with the Caputo derivative of fractional order 0 < a ≤ 2 and 0 < β ≤ 2, respectively. The Laplace transform with respect to time is used. The approximate solution valid for small values of time is obtained in terms of the Mittag-Leffler, Wright, and Mainardi functions.


Abstract:
The problem of fractional heat conduction in a composite medium consisting of a spherical inclusion ) 0 ( R r   and a matrix ) (    r R being in perfect thermal contact at R r  is considered.The heat conduction in each region is described by the time-fractional heat conduction equation with the Caputo derivative of fractional order 2 0    and , 2 0 

Introduction
The standard heat conduction (diffusion) equation for temperature is obtained from the balance equation for energy 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 with k being the thermal conductivity.In the heat conduction Equation ( 1) 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] results in the heat conduction with memory [7,8]   .
Several particular cases of choice of the memory kernel   were analyzed in [9][10][11][12].The timenonlocal dependence between the heat flux vector q and the temperature gradient with the long-tail power kernel [9][10][11][12]  where     is the gamma function, can be interpreted in terms of fractional calculus: where are the Riemann-Liouville fractional integral and derivative of the order ,  respectively [13][14][15][16]: The balance Equation (2) and the constitutive Equations ( 8) and ( 9) yield the time-fractional equation with the Caputo fractional derivative 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].
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  .
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 and the physical condition with the prescribed boundary value of the heat flux 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 with some nonzero constants 1 c and 2 c , 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 E T  , where h is the convective heat transfer coefficient, leads to 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: (23) 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 The fractional integral of the order 0   is denoted as 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

Statement of the Problem
Consider the time-fractional heat conduction equations in a spherical inclusion (24) and in a matrix and the boundary condition of perfect thermal contact     ., , : The boundedness condition at the origin and the zero condition at infinity are also assumed: The limitations on  and  in Equations (26)(27)(28)(29) express the fact that if In what follows we restrict ourselves to the particular case when a sphere R r   0 is at initial uniform temperature 0 T and the matrix The Laplace transform with respect to time t applied to Equations ( 24) and ( 25) leads to two ordinary differential equations having the solutions It follows from conditions at the origin and at infinity Equation ( 32) that .0 , 0 The integration constants 1 B and 2 B are obtained from the perfect thermal contact boundary conditions Equations ( 30) and ( 31) Hence, the solution is written as 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, , cosh sinh 1 In the following particular cases 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 .

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.
conduction in each region is described by the time- fractional heat conduction equation with the Caputo derivative of fractional order on the first time derivative should be also imposed.