Axisymmetric Solutions to Time-Fractional Heat Conduction Equation in a Half-Space under Robin Boundary Conditions
Y. Z. Povstenko1,2
Academic Editor: Nikolai Leonenko
Received23 May 2012
Revised18 Jul 2012
Accepted23 Jul 2012
Published29 Aug 2012
Abstract
The time-fractional heat conduction equation with the Caputo derivative of the order
is considered in a half-space in axisymmetric case under two types of Robin
boundary condition: the mathematical one with the prescribed linear combination of
the values of temperature and the values of its normal derivative and the physical
condition with the prescribed linear combination of the values of temperature and the
values of the heat flux at the boundary.
1. Introduction
The generalized Fourier law, the time-nonlocal dependence between the heat flux vector and the temperature gradient with the βlong-tail memoryβ power kernel [1, 2] (see also [3]) as
where is the gamma function, can be interpreted in terms of the fractional calculus:
and in combination with the law of conservation of energy as
results in the time-fractional heat conduction equation with the Caputo fractional derivative
Here, is the mass density, denotes the specific heat capacity, and is the thermal diffusivity coefficient.
Recall that the Riemann-Liouville fractional integral and derivative are defined as follows (see [4β6]):
whereas the Caputo fractional derivative has the following form [5β7]:
A detailed explanation of derivation of time-fractional heat conduction equation (1.4) from the constitutive equations (1.2) and the law of conservation of energy (1.3) can be found in [8]. Here we briefly present the main idea. In the case , as a consequence of (1.2), (1.3), and (1.6), we have
or after integration with respect to time as
Applying to both sides of (1.9) the Caputo derivative , and taking into account that for [6],
we obtain the time-fractional heat conduction equation (1.4) for .
Similarly, for , we get
Applying to both sides of (1.11) gives (1.4) for as the following equality fulfills [5]
If the heat conduction equation is investigated in a bounded domain, the boundary conditions should be prescribed. 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 considered domain
with some nonzero constants and , 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 condition of convective heat exchange between a body and the environment with the temperature
where is the convective heat transfer coefficient, leads to
The literature on mathematical aspects concerning correctness of initial-boundary-value problems for time-fractional diffusion equation and form and properties of its solutions is quite extensive (see, e.g., [9β16], among others). Geometrical explanation of fractional calculus is given in [17β19].
There are only a few papers [20, 21] in which the fractional diffusion is investigated under the mathematical Robin boundary condition. In previous publications, problems for a cylinder [22] and a sphere [23] under mathematical and physical Neumann boundary conditions were considered. In the present paper, for the first time, the solutions to time-fractional heat conduction equation in a half-space are studied under both the mathematical and physical Robin boundary conditions. The Laplace integral transform with respect to time , the Hankel transform with respect to the spatial coordinate , and the sin-cos-Fourier transforms with respect to spatial coordinate are used. The solutions under the mathematical and physical Neumann boundary conditions are obtained as particular cases.
2. Mathematical Preliminaries
2.1. Laplace Transform
The Laplace transform is defined as
where is the transform variable.
The inverse Laplace transfrom is carried out according to the Fourier-Mellin formula:
where is a positive fixed number.
The Laplace transform of the Riemann-Liouville fractional integral of the order is carried out according to the formula similar to the Laplace transform of -fold primitive of a function as
The Caputo derivative for its Laplace transform requires the knowledge of the initial values of the function and its integer derivatives of the order
whereas the Riemann-Liouville derivative for its Laplace transform rule requires the knowledge of the initial values of the fractional integral and its derivatives of the order
The reader interested in applications of integral transforms in fractional calculus is referred to [4β7, 24].
2.2. Hankel Transform
The Hankel transform is used to solve problem in cylindrical coordinates in the domain and is defined as
where is the Bessel function of the order .
The following formula is fulfilled:
2.3. Sin-Cos-Fourier Transform
In the case of boundary condition of the third kind with the prescribed boundary value of linear combination of a function and its normal derivative
the following sin-cos-Fourier transform [25] is employed:
with the kernel
Application of the sin-cos-Fourier transform to the second derivative of a function gives
3. Solution to the Problem under Mathematical Robin Boundary Condition
Consider the axisymmetric time-fractional heat conduction equation in cylindrical coordinates
with zero initial conditions
and the mathematical Robin boundary condition
The zero conditions at infinity are also assumed
The solution to the initial-boundary-value problem (3.1)β(3.4) can be written as
where is the fundamental solution being the solution of the following problem:
Here, is the Dirac delta function.
The Laplace transform with respect to time , the Hankel transform with respect to the spatial coordinate , and the sin-cos-Fourier transform with respect to the spatial coordinate result in
where each of the integral transforms is denoted by the asterisk.
Invertion of the integral transforms leads to
Here, is the Mittag-Leffler function in two parameters and defined by the series representation:
The essential role of the Mittag-Leffler function in fractional calculus results from the following formula for the inverse Laplace transform
Consider several particular cases of solution (3.8).
The case corresponds to the mathematical Neumann boundary condition with the prescribed boundary value of the normal derivative of temperature, and the solution reads
In the case of classical heat conduction (), the Mittag-Leffler function:
and taking into account the following integrals [26, 27]:
where is the modified Bessel function, is the complementary error function, we get
For the wave equation (),
and after evaluation of the necessary integrals appearing in solution (3.8) (under assumption ) [26, 27]
we obtain the solution
Of particular interest is also the case for which
4. Solution to the Problem under Physical Robin Boundary Condition
Consider the following axisymmetric time-fractional heat conduction equation:
under zero initial conditions
and the physical Robin boundary condition
where .
The solution to the initial-boundary-value problem (4.1)β(4.3) has the form
where the fundamental solution fulfills the following equation
under the following conditions:
The Laplace transform with respect to time leads to the following equation:
and the following boundary condition:
In this case, the kernel of the sin-cos-Fourier transform with respect to the spatial coordinate depends on the Laplace transform variable
and in the transform domain we obtain
Inversion of the Laplace transform in (4.10) depends on the value of . For , we have
whereas for , we get
The fundamental solution under physical Neumann boundary condition is obtained for as
Formulae (4.11) and (4.12) simplify for and , respectively, taking into account (3.12) and that
where
is the Dawson integral.
5. Conclusion
We have derived the analytical solutions to time-fractional heat conduction equation in a half-space under mathematical and physical Robin boundary conditions. The integral transform technique has been used. It should be emphasized that in the case of physical Robin boundary condition, the order of integral transforms is important as the kernel of the sin-cos-Fourier transform depends on the Laplace transform variable. The limiting case corresponds to solutions of problems under mathematical and physical Neumann boundary conditions with the prescribed boundary value of the normal derivative and with the prescribed boundary value of the heat flux, respectively. The difference between mathematical and physical boundary conditions (as well as the difference between the solutions) disappears in the case of standard heat conduction equation ().
References
Y. Z. Povstenko, βFractional heat conduction equation and associated thermal stress,β Journal of Thermal Stresses, vol. 28, no. 1, pp. 83β102, 2005.
Y. Z. Povstenko, βTheory of thermoelasticity based on the space-time-fractional heat conduction equation,β Physica Scripta, vol. T136, Article ID 014017, 6 pages, 2009.
R. Gorenflo, F. Mainardi, D. Moretti, and P. Paradisi, βTime fractional diffusion: a discrete random walk approach,β Nonlinear Dynamics, vol. 29, no. 1–4, pp. 129β143, 2002.
S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993.
I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of North-Holland Mathematics Studies, Elsevier Science, Amsterdam, The Netherlands, 2006.
R. Gorenflo and F. Mainardi, βFractional calculus: integral and differential equations of fractional order,β in Fractals and Fractional Calculus in Continuum Mechanics, A. Carpinteri and F. Mainardi, Eds., vol. 378 of CISM Courses and Lectures, pp. 223β276, Springer, Vienna, 1997.
Y. Povstenko, βNon-axisymmetric solutions to time-fractional diffusion-wave equation in an infinite cylinder,β Fractional Calculus and Applied Analysis, vol. 14, no. 3, pp. 418β435, 2011.
F. Mainardi, βThe fundamental solutions for the fractional diffusion-wave equation,β Applied Mathematics Letters, vol. 9, no. 6, pp. 23β28, 1996.
F. Mainardi, βFractional relaxation-oscillation and fractional diffusion-wave phenomena,β Chaos, Solitons and Fractals, vol. 7, no. 9, pp. 1461β1477, 1996.
R. Gorenflo and F. Mainardi, βSignalling problem and Dirichlet-Neumann map for time-fractional diffusion-wave equations,β Matimyás Matematika, vol. 21, pp. 109β118, 1998.
A. Hanyga, βMultidimensional solutions of time-fractional diffusion-wave equations,β Proceedings of the Royal Society A, vol. 458, no. 2020, pp. 933β957, 2002.
S. D. Eidelman and A. N. Kochubei, βCauchy problem for fractional diffusion equations,β Journal of Differential Equations, vol. 199, no. 2, pp. 211β255, 2004.
Y. Luchko, βSome uniqueness and existence results for the initial-boundary-value problems for the generalized time-fractional diffusion equation,β Computers & Mathematics with Applications, vol. 59, no. 5, pp. 1766β1772, 2010.
J.-H. He, βApproximate analytical solution for seepage flow with fractional derivatives in porous media,β Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 57β68, 1998.
J.-H. He, S. K. Elagan, and Z. B. Li, βGeometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus,β Physics Letters A, vol. 376, no. 4, pp. 257β259, 2012.
A. K. Bazzaev and M. Kh. Shkhanukov-Lafishev, βLocally one-dimensional scheme
for fractional diffusion equation with Robin boundary conditions,β Computational Mathematics and Mathematical Physics, vol. 50, no. 7, pp. 1141β1149, 2010.
J. Kemppainen, βExistence and uniqueness of the solution for a time-fractional diffusion equation with Robin boundary condition,β Abstract and Applied Analysis, vol. 2011, Article ID 321903, 11 pages, 2011.
Y. Povstenko, βTime-fractional radial heat conduction in a cylinder and associated thermal stresses,β Archive of Applied Mechanics, vol. 82, no. 3, pp. 345β362, 2012.
Y. Z. Povstenko, βCentral symmetric solution to the Neumann problem for a time-fractional diffusion-wave equation in a sphere,β Nonlinear Analysis: Real World Applications, vol. 13, no. 3, pp. 1229β1238, 2012.