Fractional heat conduction with heat absorption in a sphere under Dirichlet boundary condition

The time-fractional heat conduction equation with the Caputo derivative and with heat absorption term proportional to temperature is considered in a sphere in the case of central symmetry. The fundamental solution to the Dirichlet boundary value problem is found, and the solution to the problem under constant boundary value of temperature is studied. The integral transform technique is used. The solutions are obtained in terms of series containing the Mittag-Leffler functions being the generalization of the exponential function. The numerical results are illustrated graphically.


Introduction
The classical parabolic heat conduction equation with the source term proportional to tem- was considered, e.g., in Carslaw and Jaeger (1959), Crank (1975), Nyborg (1988), Polyanin (2002). Here, T is temperature, t is time, a stands for the thermal diffusivity coefficient, denotes the Laplace operator, the coefficient b describes the heat absorption (heat release).
In the last few decades, differential equations with derivatives of non-integer order attract the attention of the researchers as such equations provide a very suitable tool for description of many important phenomena in physics, geophysics, chemistry, biology, engineering and solid mechanics (see, for example, Gafiychuk et al. 2008;Herrmann 2011;Magin 2006;Mainardi 2010;Povstenko 2015a;Sabatier et al. 2007;Tarasov 2010;Tenreiro Machado 2011;Uchaikin 2013).
In this paper, we consider the time-fractional equation where is the Caputo fractional derivative (Kilbas et al. 2006;Podlubny 1999). Equation (2) results from the time-nonlocal generalization of the Fourier law with the "long-tail" power kernel. Such a generalization can be interpreted in terms of derivatives and integrals of non-integer order. Equation (2) takes into consideration the memory effects with respect to time. The interested reader is also referred to Povstenko (2011Povstenko ( , 2015a.
In the literature, there are papers devoted to investigation of Eq.
(2) in the case of one Cartesian spatial coordinate (Damor et al. 2016;Ferrás et al. 2015;Qin and Wu 2016;Vitali et al. 2017). Here, we study Eq. (2) in a spherical domain 0 ≤ r < R in the case of central symmetry under Dirichlet boundary condition.
The paper is organized as follows. In Sect. 2, we find the fundamental solution to the Dirichlet boundary value problem using the Laplace transform with respect to time t and the finite sin-Fourier transform with respect to the spatial coordinate r . The corresponding problem under constant boundary value of temperature at the surface r = R is investigated in Sect. 3. In both cases, the numerical results are illustrated graphically. Conclusions are presented in Sect. 4.

The fundamental solution to the Dirichlet problem
We consider the time-fractional heat conduction equation with heat absorption term in spherical coordinates in the case of central symmetry in the domain 0 ≤r < R, 0< t < ∞ with a > 0 and the order of the Caputo derivative 0 <α ≤ 1. Equation (4) is studied under zero initial condition and the Dirichlet boundary condition with δ(t) being the Dirac delta function. The constant multiplier p 0 is introduced in (6) to obtain the non-dimensional quantities used in numerical calculations [see Eq. (17)].
In the general case, the Caputo fractional derivative has the following Laplace transform rule where the asterisk denotes the Laplace transform, s is the transform variable.
In what follows, we will use the finite sin-Fourier transform with respect to the spatial coordinate r in the domain 0 ≤ r ≤ R (see, for example, Povstenko 2015b): where This transform is the convenient reformulation of the sin-Fourier series and is used for Dirichlet boundary condition with the prescribed boundary value of a function, since for the Laplace operator in the case of central symmetry we have Applying the integral transforms to (4) under the initial condition (5) and boundary condition (6) gives in the transform domain: Inversion of the Laplace and finite sin-Fourier transforms results in the sought-for fundamental solution: To obtain (13) the following formula has been used, where E α,β (z) is the Mittag-Leffler function in two parameters α and β described by the series representation (Gorenflo et al. 2014;Kilbas et al. 2006;Podlubny 1999) For b = 0, the fundamental solution Eq. (13) coincides with the solution obtained in Povstenko (2008Povstenko ( , 2015b. In the particular case α = 1, the Mittag-Leffler function E 1,1 (z) = e z ; hence, the fundamental solution to the Dirichlet problem for the classical heat conduction equation with heat absorption has the form Introducing the non-dimensional quantities we get the solution in the general case and for the classical heat conduction respectively. The results of numerical calculations are shown in Figs. 1, 2, 3, and 4. Numerical results demonstrate the significant influence of the order of fractional derivative α and the absorption parameter b on the heat conduction process. When the fractional derivative of order 0 < α < 1 replaces the standard first-order time derivative in the diffusion (heat conduction) equation, this leads to slow diffusion (see, for example, Chen 2017; Kimmich 2002; Metzler and Klafter 2004). Figures show evidently the slower diffusion with decreasing α. Heat absorption also results in slower heat diffusion.

The constant boundary value of temperature
Now we solve the time-fractional heat conduction equation with the heat absorption term under zero initial condition and the Dirichlet boundary condition with the constant boundary value of temperature: As above, the Laplace transform with respect to time t and the finite sin-Fourier transform with respect to the spatial coordinate r give the solution in the transform domain: In connection with we obtain and after inverting the integral transforms we arrive at In this case E α (z) is the Mittag-Leffler function in one parameter α having the series representation Taking into account the following series (Prudnikov et al. 1986) we get The advantage of the solution (29) is that the first term satisfies the boundary condition (22), whereas the second term vanishes at r = R. In terms of non-dimensional quantities, we have where T = T /T 0 , other non-dimensional parameters are the same as in (17).

Conclusions
In this paper, we have investigated the time-fractional heat conduction equation with the Caputo derivative of order 0 < α ≤ 1 and the heat absorption term proportional to temperature. The fundamental solution to the Dirichlet boundary problem and the solution to the problem with constant boundary value of temperature have been found. The solutions have been obtained in terms of series containing the Mittag-Leffler functions being the generalization of the exponential function. To evaluate the Mittag-Leffler functions E α (z) and E α,α (z) we have used the algorithms proposed in  (see also the Matlab programs that implement these algorithms Matlab File Exchange 2005). The obtained results can be generalized in the future works for isotropic fractal media within the framework of the non-integer dimensional space approach.