Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 8605056, 7 pages

http://dx.doi.org/10.1155/2016/8605056

## Time-Fractional Heat Conduction in a Half-Line Domain due to Boundary Value of Temperature Varying Harmonically in Time

Institute of Mathematics and Computer Science, Faculty of Mathematics and Natural Sciences, Jan Dlugosz University in Czestochowa, Armii Krajowej 13/15, 42-200 Czestochowa, Poland

Received 23 August 2016; Accepted 23 October 2016

Academic Editor: Filippo de Monte

Copyright © 2016 Yuriy Povstenko. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The Dirichlet problem for the time-fractional heat conduction equation in a half-line domain is studied with the boundary value of temperature varying harmonically in time. The Caputo fractional derivative is employed. The Laplace transform with respect to time and the sin-Fourier transform with respect to the spatial coordinate are used. Different formulations of the considered problem for the classical heat conduction equation and for the wave equation describing ballistic heat conduction are discussed.

#### 1. Introduction

In the paper [1] and later on in the book [2] Nowacki studied the classical parabolic heat conduction equation with a heat source term varying harmonically as a function of timein the domain . Here is the thermal diffusivity coefficient, is the Dirac delta function, and denotes the frequency.

Nowacki’s solution of (1) is based on the assumption that temperature can be expressed as a product of the auxiliary function and the time harmonic termIn this case, there are no initial and boundary conditions (excepting the zero condition at ), and the problem is reduced to solving the corresponding equation for the auxiliary function The final result reads The square root of the imaginary unit is defined as

If the heat conduction equationis considered in a half-line domain , then the boundary condition at should be imposed. For example, we can assume the Dirichlet boundary condition varying harmonically in time:Similar analysis can be also carried out in the case of the boundary value of heat flux varying harmonically in time (the physical Neumann boundary condition). Boundary conditions varying harmonically in time describe various situations, in particular, thermal processing of materials using pulsed lasers or collection of solar energy [3].

Under Nowacki’s assumption (2), there is no initial condition, and for the auxiliary function we obtainTo compare with the subsequent results it is worthwhile to solve (6) under boundary conditions (7) and (8) using the sin-Fourier transform with respect to the spatial coordinate . The solution has a form Hence,If the surface temperature is described by the dependencethen the solution becomes [3]

Many experimental and theoretical investigations testify that in media with complex internal structure the standard heat conduction equation is no longer sufficiently accurate. This results in formulation of nonclassical theories, in which the parabolic heat conduction equation is replaced by more general one (see [4–12] and the references therein).

For example, Green and Naghdi [7] proposed the theory of thermoelasticity without energy dissipation based on the wave equation for temperature. In the framework of this theory, the following boundary value problem can be studied:Under the assumption (2), using the sin-Fourier transform, we get the solution (see (A.3) from appendix)

#### 2. Time-Fractional Heat Conduction

The time-nonlocal generalization of the Fourier law with the “long-tail” power kernel [11, 13–15] can be interpreted in terms of fractional calculus (theory of integrals and derivatives of noninteger order) and results in the time-fractional heat conduction equationwith the Caputo fractional derivative of order defined as [16–18]and having the following Laplace transform rule:where the asterisk denotes the transform, is the transform variable, and is the gamma function. The Caputo fractional derivative is a regularization in the time origin for the Riemann-Liouville fractional derivative by incorporating the relevant initial conditions [19]. The major utility of the Caputo fractional derivative is caused by the treatment of differential equations of fractional order for physical applications, where the initial conditions are usually expressed in terms of a given function and its derivatives of integer (not fractional) order, even if the governing equation is of fractional order [17, 20]. Additional discussion on the use of the Caputo and Riemann-Liouville fractional derivatives can be found in [21] (see Section 3.4 “Which type of fractional derivative? Caputo or Riemann-Liouville?” in this book).

Equations with fractional derivatives describe many important physical phenomena in amorphous, colloid, glassy, and porous materials, in fractals, comb structures, polymers, and random and disordered materials, in viscoelasticity and hereditary mechanics of solids, in biological systems, and in geophysical and geological processes (see, e.g., [22–30] and the references therein). Important applications of fractional calculus can be found in such fields as fractional dynamics [31–35], fractional kinetics [36–38], and fractional thermoelasticity [11, 12, 39–41].

Equation (15) describes the whole spectrum from localized heat conduction (the Helmholtz equation for ) through the standard heat conduction () to the ballistic heat conduction (the wave equation when ).

The interested reader is referred to the book [15], which systematically presents solutions to different initial and boundary value problems for the time-fractional diffusion-wave equation (15) in Cartesian, cylindrical, and spherical coordinates. In [42, 43], this equation was considered in unbounded domains with the source term varying harmonically in time.

In the present paper, we study the Dirichlet problem for the time-fractional heat conduction equation in a half-line domain with the surface value of temperature varying harmonically in time. The integral transform technique is used. Different formulations of the considered problem for the classical heat conduction equation () and for the wave equation describing ballistic heat conduction () are discussed.

#### 3. Formulation of the Problem

We consider the time-fractional heat conduction equation in a half-line:under the harmonic boundary conditionand zero condition at infinity

For the Caputo derivative of the exponential function we haveSubstituting (with ) gives the final resultwhere is the incomplete gamma function [44]:

Hence, for fractional (noninteger) values of the order of derivative,Therefore, the Nowacki assumption (2) cannot be used for the time-fractional heat conduction equation and the corresponding initial conditions should be imposed (see also [42, 43]). In the present paper we assume zero initial conditions:

#### 4. Solution to the Problem

Application of the Laplace transform with respect to time and the sin-Fourier transform with respect to the spatial coordinate to (18) under the initial conditions (25) and the boundary conditions (19), (20) giveswhere the tilde denotes the sin-Fourier transform and is the transform variable.

At first, we analyze the standard heat conduction equation corresponding to :The inverse Laplace transform gives (see (A.6) from appendix):Integrals (A.2) and (A.4) allow us to invert the sin-Fourier transform and to obtain the solutionThe first term in (29) coincides with solution (10) and describes the quasi-steady-state oscillations; the second and third ones describe the transient process.

Starting in (27) from the inversion of the sin-Fourier transform, we get

Taking into account (A.10) and using the convolution theorem for the Laplace transform allow us to obtain an alternative form of the solution to the standard heat conduction equation:

Another particular case of solution (26) in the transform domain corresponds to the ballistic heat conduction ():Inversion of the sin-Fourier transform (see (A.2)) results inTaking into account (A.9), we obtainIt should be noted that solution (34) describes the wavefront at

In applications there often appears the value (see investigations of diffusion on fractals [45] and comb structures [46, 47]). Using (A.7), we get

Now we return to the analysis of the time-fractional heat conduction equation and its solution in the transform domain (26). The inverse sin-Fourier transform givesThe inverse Laplace transform of is expressed in terms of the Mainardi function (see (A.11) from appendix). The solution has the form

The particular case of the Mainardi function is reduced to the exponential function (see [48, 49]): and solution (37) for coincides with solution (31).

Starting from the inverse Laplace transform of (26), we havewhere is the Mittag-Leffler function in two parameters and (see appendix). In parallel with (37), the inverse Fourier transform of (39) leads to another form of the solutionComparison of (37) and (40) allows us to establish the relation between the Mainardi function and the Mittag-Leffler function in the form of sin-Fourier transform (see also [15, 50], where the similar relations were obtained in terms of the cos-Fourier transform).

Figures 1 and 2 present the dependence of solution on distance in the case of the boundary condition for different values of the order of fractional derivative and different values of time. In numerical calculations we have used the following nondimensional quantities:To evaluate the Mittag-Leffler function , the algorithm suggested in the paper [51] has been used.