Fractional Order Generalized Thermoelastic Infinite Medium with Cylindrical Cavity Subjected to Harmonically Varying Heat

In this work, a mathematical model of an elastic material with cylindrical cavity will be constructed. The governing equations will be taken into the context of the fractional order generalized thermoelasticity theory (Youssef 2010). Laplace transform and direct approach will be used to obtain the solution when the boundary of the cavity is exposed to harmonically heat with constant angular frequency of thermal vibration. The inverse of Laplace transforms will be computed numerically using a method based on Fourier expansion techniques. Some comparisons have been shown in figures to present the effect of the fractional order parameter and the angular frequency of thermal vibration on all the studied felids.


Introduction
Recently, a considerable research effort has been expended to study anomalous diffusion, which is characterized by the time-fractional diffusion-wave equation by Kimmich [1]: where  is the mass density, c is the concentration,  is the diffusion conductivity, i is the coordinate symbol which takes the values 1, 2 and 3, the subscript "," means the derivative with respect to x i and notion I  is the Riemann-Liouville fractional integral is introduced as a natural generalization of the well-known n-fold repeated integral   written in a convolution-type form as in [2,3]: where     is the gamma function.According to Kimmich [1] Equation (1) describes different cases of diffusion where 0 1 It should be noted that the term diffusion is often used in a more generalized sense including various transport phenomena.Equation ( 1) is a mathematical model of a wide range of important physical phenomena, for example, the sub-diffusive transport occur in widely different systems ranging from dielectrics and semiconductors through polymers to fractals, glasses, porous, and random media.Super diffusion is comparatively rare and has been observed in porous glasses, polymer chain, biological systems, transport of organic molecules and atomic clusters on surface [4].One might expect the anomalous heat conduction in media where the anomalous diffusion is observed.
Fujita [5,6] considered the heat wave equation for the case of 1 2 where C is the specific heat, k is the thermal conductivity and the subscript "," means the derivative with respect to the coordinates i x .Equation (3) can be obtained as a consequence of the non local constitutive equation for the heat flux components i q is in the form 1 , , 1 2 (4) Povstenko [4] used the Caputo heat wave equation defines in the form: to get the stresses corresponding to the fundamental solution of a Cauchy problem for the fractional heat conduction equation in one-dimensional and two-dimensional cases.
Some applications of fractional calculus to various problems of mechanics of solids are reviewed in the literature [7,8].

Theory of Fractional Order Generalized Thermoelasticity
The classical thermoelasticity is based on the principles of the theory of heat conduction which is called Fourier law, which relates the heat flux components i q to the temperature gradient as follows: In combination with the energy conservative law, this leads to the parabolic heat conduction equation which is considered by Povstenko [4]: where dotted above T means the derivative with respect to the time t.
Recently, in the non classical thermoelasticity theories, Fourier law (6) and heat conduction (7) are replaced by more general equations, have been formulated.The first well-known generalized of such a type of Lord and Shulman [9] and it takes the form: which leads to the hyperbolic differential equation of heat conduction of Lord and Shulman [9]: where o  is non-negative constant and is called relaxation time.
In the context of the generalized thermoelasticity, the governing equations for isotropic medium are defined as follows: The equation of motion The constitution relation where , is the increment of the dynamical temperature where o T is the reference temperature,   and ij e is the strain tensor satisfy the relations The heat flux equation The entropy increment equation per unit volume takes the form where S is the entropy increment of the material.The heat flux-entropy equation The heat equation without any heat sources By using Equations (14,15,16), we have the heat equation in the form [11]: where for weak conductivity for normal conductivity for strong conductivity

The Problem Formulation
Let us consider a perfectly conducting elastic infinite body with cylindrical cavity occupies the region R r    of an isotropic homogeneous medium whose state can be expressed in terms of the space variable r and the time variable t such that all of the state functions vanish at infinity.We will use the cylindrical system of coordinates (r,,z) with the z-axis lying along the axis of the cylinder.Due to symmetry, the problem is one-dimensional with all the functions considered depending on the radial distance r and the time t.The medium described above is considered to be quiescent and the surface of the cavity is subjected to harmonically varying heat and traction free described mathematically as follow:  is constant and  is the angular frequency of thermal vibration ( 0   for a thermal shock).
It is assumed that there are no body forces and no heat sources in the medium.Thus, the field equations ( 10), ( 11), ( 12) and (17) in cylindrical case can be set as [12]:   where For convenience, we shall use the following non-dimensional variables [12]: where 2  2 Equations (20-26) assume the form (where the primes are suppressed for simplicity) where

Formulation in the Laplace Transform Domain
Taking the Laplace transform for the both sides of the Equations (27-31), this is defined as follows: where the rule for the Laplace transform of the Riemann-Liouville fractional integral for zero initial function reads from Povstenko [4]: Then, we have 4 e where and an over bar symbol denotes its Laplace transform and s denotes the Laplace transform parameter.
Eliminating u from the Equations (26,34,35), we get Eliminating e from Equations (41,42), we obtain In a similar manner, we can show that e satisfies the equation The bounded solutions of Equations (41,42) at infinity can be written in the form where K o (.) is the modified Bessel function of the second kind of order zero.A 1 , A 2 , B 1 and B 2 are all parameters depending on the parameter s of the Laplace transform, and satisfy the relations Using Equation (41), we obtain 2 2 ,i=1,2.
Thus, we have Substituting from Equation (49) into the Laplace transform of Equation ( 26), obtain where K 1 (.) is the modified Bessel function of the second kind of order one.
In deriving Equation (50), we have used the following well-known relation of the Bessel function: Finally, substituting from Equations (45,49,50) into Equations (36-38), we obtain the stress components in the form: Using the boundary conditions (39,40), we get Then, we have where Then, we have and where Those complete the solution in the Laplace transform space.

Numerical Inversion of the Laplace Transform
In order to invert the Laplace transform, we adopt a numerical inversion method based on a Fourier series expansion [13].By this method the inverse ( ) where N is a sufficiently large integer representing the number of terms in the truncated Fourier series, chosen such that where  1 is a prescribed small positive number that corresponds to the degree of accuracy required.The parameter c is a positive free parameter that must be greater than the real part of all the singularities of   f s .The optimal choice of c was obtained according to the criteria described in [13].

Numerical Results and Discussion
With a view to illustrating the analytical procedure presented earlier, we now consider a numerical example for which computational results are given.The results depict the variation of temperature, stress, displacement and strain fields in the context of Youssef model [11].
For this purpose, copper is taken as the thermoelastic material for which we take the following values of the different physical constants [12]: From the above values we get the nondimensional values for our problem as: .We can see the significant effect of the parameter  and the angular frequency of thermal vibration  on all the studied fields.

Conclusion
We considered a perfectly conducting elastic isotropic homogeneous infinite body with cylindrical cavity in the context of the fractional order generalized thermoelasticity theory (Youssef model).The effect of the fractional parameter and the angular frequency of thermal vibration on all the studied fields are very significant.New classification of the materials must be constructed according to the fractional parameter which describes the ability of the material to conduct the heat.

2 1 p and 2 2 p
are the roots of the characteristic equation
the temperature, the stresses, the displacement and the strain are shown in Figures 1-5 respectively with wide range of non dimensional distance r from r = R = 1.0 up to r = 2.0, non dimensional time t = 0.08 and non dimensional relaxation time 0.001 o   with different values of the parameter   0.8, 1.0,1.2  which describe the three types of conductivity (weak conductivity, normal conductivity, strong conductivity), respectively and with different values of the parameter