International Scholarly Research Notices
Volume 2014 (2014), Article ID 646049, 13 pages
http://dx.doi.org/10.1155/2014/646049

Research Article

Fractional Order Two-Temperature Dual-Phase-Lag Thermoelasticity with Variable Thermal Conductivity

Sudip Mondal,¹ Sadek Hossain Mallik,² and M. Kanoria³

¹Bhatkunda High School, Burdwan 713153, India
²Department of Mathematics, Aliah University, Kolkata 700091, India
³Department of Applied Mathematics, University of Calcutta, Kolkata 700009, India

Received 10 May 2014; Revised 4 August 2014; Accepted 4 August 2014; Published 29 October 2014

Academic Editor: Abdelouahed Tounsi

Copyright © 2014 Sudip Mondal et al. 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.

A new theory of two-temperature generalized thermoelasticity is constructed in the context of a new consideration of dual-phase-lag heat conduction with fractional orders. The theory is then adopted to study thermoelastic interaction in an isotropic homogenous semi-infinite generalized thermoelastic solids with variable thermal conductivity whose boundary is subjected to thermal and mechanical loading. The basic equations of the problem have been written in the form of a vector-matrix differential equation in the Laplace transform domain, which is then solved by using a state space approach. The inversion of Laplace transforms is computed numerically using the method of Fourier series expansion technique. The numerical estimates of the quantities of physical interest are obtained and depicted graphically. Some comparisons of the thermophysical quantities are shown in figures to study the effects of the variable thermal conductivity, temperature discrepancy, and the fractional order parameter.

1. Introduction

During the last five decades, nonclassical thermoelasticity theories involving hyperbolic type heat transport equations admitting finite speed of thermal signals have been formulated. According to these theories, heat propagation is to be viewed as a wave phenomenon rather than a diffusion phenomenon.

In order to overcome the paradox of an infinite speed of thermal wave inherent in CTE and CCTE (classical coupled theory of thermoelasticity), efforts were made to modify coupled thermoelasticity, on different grounds, to obtain a wave-type heat conduction equation by different researchers. Lord and Shulman [1] formulated the generalized thermoelasticity theory introducing one relaxation time in Fourier’s law of heat conduction equation and thus transforming the heat conduction equation into a hyperbolic type.

Green and Lindsay [2] used the theory of two different relaxation times in the constitutive relations for the stress tensor and the entropy equation. Later Green and Naghdi [3–5] have proposed three models, labeled as types I, II, and III. When they are linearized, type I is the same as the classical heat equation (based on Fourier’s law) whereas the linearized versions of type-II and type-III theories permit propagation of thermal waves at finite speed. The entropy flux vector in type-II (i.e., thermoelasticity without energy dissipation) and type-III (i.e., thermoelasticity with energy dissipation) models is determined in terms of the potential that also determines stresses. When Fourier conductivity is dominant, then the temperature equation reduces to classical Fourier’s law of heat conduction and when the effect of conductivity is negligible, then the equation has undamped thermal wave solutions without energy dissipation.

Tzou [6] introduced two-phase-lag models to both the heat flux vector and the temperature gradient. According to this model, classical Fourier’s law has been replaced by , where the temperature gradient at a point of the material at time corresponds to the heat flux vector at the same point at time . Here is the thermal conductivity of the material. The delay time is interpreted as that caused by the microstructural interactions and is called the phase-lag of the temperature gradient. The other delay time is interpreted as the relaxation time due to the fast transient effects of thermal inertia and is called the phase-lag of the heat flux. The case corresponds to classical Fourier’s law. If and , Tzou refers to the model as single-phase-lag model. Recently Several researchers have attempted to solve their problems on the basis of the theory of dual-phase-lag model. Roychoudhuri [7] has studied one-dimensional thermoelastic wave propagation in an elastic half-space in the context of dual-phase-lag model. The exponential stability [8] and condition of the delay parameters in the dual-phase-lag theory [9] under this model have been studied by Quintanilla. Wang et al. [10, 11] have studied the well-posedness and solution structure of the dual-phase-lag heat conduction equation. Wang and Mingtian [12] have studied the thermal oscillation and resonance in dual-phase-lag heat conduction equation. Ailawalia and Budhiraja have studied a problem of dual-phase-lag model with internal heat source [13].

Gurtin and Williams [14, 15] have suggested that there is no a priori ground for assuming that the second law of thermodynamics for continuous bodies involves only a single temperature; that is, it is more logical to assume a second law in which the entropy contribution due to heat conduction is governed by one temperature, that of the heat supply by another.

Chen and Gurtin [16] and Chen et al. [17, 18] have formulated a theory of heat conduction in deformable bodies, which depends on two distinct temperatures: the conductive temperature and the thermodynamic temperature . For time-independent situations, the differences between these two temperatures are proportional to the heat supply, and in the absence of any heat supply, the two temperatures are identical [18]. For time-dependent problems, however, and for wave propagation problems in particular, the two temperatures are, in general, different, independent of the presence of a heat supply. The key element that sets the two-temperature thermoelasticity (2TT) apart from the classical theory of thermoelasticity (CTE) is the material parameter , called the temperature discrepancy [18]. Specifically, if , then and the field equations of the 2TT reduce to those of CTE.

The linearized version of the two-temperature theory (2TT) has been studied by many authors. Warren and Chen [19] have investigated the wave propagation in the two-temperature theory of thermoelasticity. Lesan [20] has established uniqueness and reciprocity theorems for the 2TT. Puri and Jordan [21] have studied propagation of plane waves under the 2TT. The existence, structural stability, and spatial behavior of the solution in 2TT have been discussed by Quintanilla [22].

It should be pointed out that both CTE and 2TT suffer from the so-called paradox of heat conduction, that is, the prediction that a thermal disturbance at some point in a body is felt instantly, but unequally, throughout the body. Mondal et al. [23, 24] have studied problems on two-temperature Green Naghdi III and dual-phase-lag model with variable thermal conductivity. Pal et al. [25], Islam et al. [26], Das and Kanoria [27], and Banik and Kanoria [28] have studied on two-temperature generalized thermoelasticity. Kumar et al. [29, 30] have established variational and reciprocal principles and some theorems in two-temperature generalized thermoelasticity. Recently Ailawalia et al. have solved a dynamic problem on Green-Naghdi (Type III) half-space with two-temperature [31] theory.

In this work, we have studied the thermoelastic stress, strain, displacement, conductive temperature, and the thermodynamic temperature in an infinite, isotropic, homogeneous elastic half-space under thermal shock using two-temperature dual-phase-lag generalized thermoelasticity in the context of fractional heat conduction equation. The governing equations of two-temperature generalized thermoelasticity theory are formed in the Laplace transform domain which are then solved by state space approach. The inversion of the transformed solutions is carried out numerically, applying a method based on a Fourier series expansion technique [32]. Finally, numerical estimate for different thermophysical quantities are obtained for copper materials. A complete and comprehensive analysis of the results is presented and the effects of fractional order parameter, two-type temperature, and variable thermal conductivity are discussed.

2. Development of Fractional Order Theory

Differential equations of fractional order have been the focus of many studies due to their frequent appearance in various applications in fluid mechanics, viscoelasticity, biology, physics, and engineering. The most important advantage of using fractional differential equations in these and other applications is their nonlocal property. It is well-known that the integer order differential operator is a local operator but the fractional order differential operator is nonlocal. This means that the next state of a system depends not only upon its current state, but also upon all of its historical states. This is more realistic, and this is one reason why fractional calculus has become more and more popular [33–35].

Fractional calculus has been used successfully to modify many existing models of physical processes. One can state that the whole theory of fractional derivatives and integrals was established in the second half of the nineteenth century. The first application of fractional derivatives was given by Abel who applied fractional calculus in the solution of an integral equation that arises in the formulation of the tautochrone problem. The generalization of the concept of derivative and integral to a noninteger order has been subjected to several approaches, and some various alternative definitions of fractional derivatives have appeared [36–39]. In the last few years, fractional calculus has been applied successfully in various areas to modify many existing models of physical processes, for example, chemistry, biology, modeling and identification, electronics, wave propagation, and viscoelasticity. One can refer to Podlubny [35] for a survey of applications of fractional calculus.

Youssef [40] introduced the formula of heat conduction given by

where the notation is the Riemann-Liouville fractional integral, introduced as a natural generalization of the well-known -fold repeated integral written in a convolution-type form as in [41] which is written as follows:

where is the Gamma function. A uniqueness theorem also has been proved. Sur and Kanoria have employed the theory to study problems on functionally graded [42] and viscoelastic material [43].

Ezzat and El-Karamany [44, 45] established a new model of fractional heat conduction equation by using the new Taylor series expansion of time-fractional order, developed by Jumarie [46] as

El-Karamany and Ezzat [47] introduced two general models of fractional heat conduction law for a nonhomogeneous anisotropic elastic solid. Uniqueness and reciprocal theorems are proved, and the convolutional variational principle is established and used to prove a uniqueness theorem with no restriction on the elasticity or thermal conductivity tensors except symmetry conditions. For fractional thermoelasticity not involving two-temperature, El-Karamany and Ezzat [48] established the uniqueness, reciprocal theorems and convolution variational principle. For two-temperature theory the formula of heat conduction has been replaced by

where is the conductive temperature. Two-temperature fractional order thermoelasticity problem of LS and Green Naghdi models (types II and III) have been solved by Sur and Kanoria [49].

Several researchers have solved different problems [50–54] using fractional order generalized thermoelasticity theory. More detailed discussion on the subject is available in the books of Hetnarski and Eslami [55], Eslami et al. [56], and Ignaczak and Ostoja-Starzewski [57].

3. Basic Formulation

The constitutive equations are

where and .

In the context of two-temperature dual-phase-lag (2TDPL) generalized thermoelasticity theory, the equation of motion in the absence of body forces and the heat conduction equation in absence of heat sources for a linearly isotropic generalized thermoelastic solid based on the theory of fractional integral are, respectively, given by [58]

where is the density, and are Lamé’s constants, is thermal conductivity, ,   being the coefficient of linear thermal expansion, is the reference temperature, and ,   being the specific heat at constant strain.

The relation between conductive temperature and thermodynamic temperature is given by

where (≥0) is the two-temperature parameter, called temperature discrepancy.

We will consider the thermal conductivity as a linear function of thermodynamical temperature as follows [59]:

where is a constant which is equal to the thermal conductivity of the material when it does not depend on thermodynamical temperature and is a nonpositive small parameter.

Substituting from (8) into (9), we get

We now use the following mapping [59]:

Differentiating (11) with respect to , we get

Differentiating again the above equation, we obtain

Differentiating (12) with respect to , we get

Differentiating (12) with respect to , we get

Substituting from (14) and (16) into (7), we obtain

Neglecting the last term on right hand side of the above equation due to nonlinearity, we get

Using (15) in (6), we get

For linearity we can approximate last equation to the following form:

Now to transform (8) by using (11) and (12), we first replace the dummy variable with and then differentiating with respect to and finally multiplying by we get

Now substituting from (10) into (21) we have

For linearity we can approximate the last equation as

Retaining only the linear terms, (14) can be written as

Now substituting from (13), (15), and (24) into (23) we have

Now integrating with respect to we get

4. Formulation of the Problem

We consider a half space with -axis pointing to the medium. This half-space is subjected to thermal and mechanical loads on the bounding plane that depends on the time and is linearly quiescent. We will consider one-dimensional thermoelastic deformation of the body so that the displacement components can be taken in the following form:

The strain displacement relation is

and the constitutive relation (5) takes the form

The equation of motion, heat transport equation, and relation between conductive temperature and thermodynamic temperature can be written as

We now use the following nondimensional variables, to make the above equations nondimensional:

Then the corresponding nondimensional equations, after omitting the primes, are

where

Initial and regularity conditions for the problem are given by

5. Method of Solution

Applying the Laplace transform defined by

to both the sides of (32), we obtain

Eliminating from (37)–(39) we get

where

5.1. State Space Approach

The equations (40) can be written in the form of a vector matrix differential equations [60] as follows:

where

The formal solution of system (42) bounded at infinity can be written as

where

We will use the well-known Cayley-Hamilton theorem to find the form of the matrix

The characteristic equation for the matrix can be written as

The roots of this equation, namely, and , satisfy the following relations:

The Taylor series expansion for the matrix exponential in (44) is given by

Using Cayley-Hamilton theorem, we can express and higher powers of the matrix in terms of and , where is the unit matrix of the second order.

Thus the infinite series in (49) can be reduced to the form

where and are coefficients depending on and .

By Cayley-Hamilton theorem, the characteristic roots and of the matrix must satisfy (50); thus we have

By solving the above linear system of equations, we get

Hence from (50) we get

where

Using (53) we can write the solution in (44) in the following form:

Hence the solution for and can be obtained from (55) as follows:

Using (56), the solution for can be obtained from (39) as follows:

5.2. Application to Thermal Shock Problem

We will consider the bounding plane of the medium at subjected to thermal shock in the following nondimensional form:

where is constant. Now applying the Laplace transform to (59) we get

Using (9), (11), and (60) we get

where .

5.3. The Mechanical Boundary Condition

The mechanical boundary condition is taken in the form

This implies

Applying the boundary conditions (61) and (63) to (56)–(58) we get

Displacement component can be obtained from (28) using (65) in the following form:

Using (9), (11), and (64) the solution for can be obtained as follows:

Again using (9), (12), and (66) the solution for can be obtained in the following form:

The solution for can be obtained from (36) using (65) and (69) in the following form:

This completes the solution of the thermal shock problem in Laplace transform domain.

6. Numerical Inversion of Laplace Transform

It is difficult to find the analytical inverse of Laplace transform of the complicated solutions for the displacement, thermodynamic temperature, conductive temperature, stress, and strain in Laplace transform domain. So we have to resort to numerical computations. We now outline the numerical procedure to solve the problem. Let be the Laplace transform of a function .

Then the inversion formula for Laplace transform can be written as

where is an arbitrary real number greater than real parts of all the singularities of .

Taking , the preceding integral takes the form

Expanding the function in a Fourier series in the interval we obtain the approximate formula [32]

where

The discretization error can be made arbitrary small by choosing large enough [32]. Since the infinite series in (74) can be summed up to a finite number of terms, the approximate value of becomes

Using the preceding formula to evaluate we introduce a truncation error that must be added to the discretization error to produce total approximation error.

Two methods are used to reduce the total error. First the “Korrektur” method is used to reduce the discretization error. Next the -algorithm is used to accelerate convergence [32].

The Korrektur method uses the following formula to evaluate the function :

where the discretization error . Thus, the approximate value of becomes

where is an integer such that .

We will now describe the -algorithm that is used to accelerate the convergence of the series in (75). Let , where is a natural number, and let be the sequence of partial sum of series in (75).

We define the -sequence by

It can be shown that [32] the sequence

converges to faster than the sequence of partial sums ,   .

The actual procedure used to invert the Laplace transform consists of using (77) together with the -algorithm. The values of and are chosen according to the criterion outlined in [32].

7. Numerical Results and Discussion

To get the solution for strain , thermal displacement component , conductive temperature , thermodynamic temperature , and thermal stress in the space time domain we have applied Laplace inversion formula to (65), (67), (68), (69), and (70), respectively, which have been done numerically using a method based on Fourier series expansion technique [32]. The numerical code has been prepared using Fortran 77 programming language. For computational purpose copper material has been taken into consideration. The values of the material constants are taken as follows [61]:

Also we take time , for computational purpose.

Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, and 15 are drawn to represent the variation of said thermophysical quantities versus the space variable for different (=0,−2,−4), (=, and (=). Here , and corresponds to week conductivity, normal conductivity, and super conductivity, respectively; and corresponds to one-temperature and two-temperature theory, respectively.

Figure 1: The thermodynamic temperature for different at .

Figure 2: The conductive temperature for different at .

Figure 3: The displacement for different at .

Figure 4: The stress for different at .

Figure 5: The strain for different at .

Figure 6: The thermodynamic temperature for different at .

Figure 7: The conductive temperature for different at .

Figure 8: The displacement for different at .

Figure 9: The stress for different at .

Figure 10: The strain for different at .

Figure 11: The comparison of thermodynamic temperature between 1TT and 2TT at .

Figure 12: The comparison of conductive temperature between 1TT and 2TT at .

Figure 13: The comparison of displacement between 1TT and 2TT at .

Figure 14: The comparison of stress between 1TT and 2TT at .

Figure 15: The comparison of radial strain between 1TT and 2TT at .

Figures 1–5 show the effect of on the said five thermophysical quantities for two-temperature theory and fractional order parameter . From these figures it is clear that magnitude of all the quantities, that is, thermodynamic temperature , conductive temperature , displacement , stress component , and strain component , is greater for smaller magnitude of .

Figures 6–10 show the effect of for two-temperature theory and on those five quantities. Figures 6–9 show that the magnitude of thermodynamic temperature , conductive temperature , displacement component , and stress component has greater value for smaller magnitude of . But from Figure 10 it is observed that in the region (approximate) strain component has larger value for than for , which is again larger than for .

Figures 11–15 are drawn to compare between the results of one-temperature theory and two-temperature theory for and for five different thermophysical quantities. Figures 11, 13, 14, and 15 show that the magnitude of is greater for one-temperature case than two-temperature case. The only exception is here.

It is observed that at the boundary plane , (Figures 2, 7, and 12), and (Figures 5, 10, and 15), which satisfies our theoretical boundary condition. It ensures the correctness of the numerical code used. In Figures 6–10 the results for agree with the corresponding results of Mondal et al. [24].

Figures 16 and 17 represent the variation of the thermodynamic temperature and conductive temperature against time for different value of , namely, , when , , and . It is observed from the figures that the magnitudes of and are oscillatory in nature and the magnitude of peek of oscillation decreases with time.

Figure 16: The variation of thermodynamic temperature with time at for ; .

Figure 17: The variation of conductive temperature with time at for ; .

8. Conclusion

State space approach has been applied to solve a generalized thermoelastic problem of an isotropic, half-space with variable thermal conductivity. The boundary of the half-space is subjected to thermal and mechanical loads. Variation of thermal conductivity has been taken as linear function of temperature. The problem has been studied using the two-temperataure dual-phase-lag model of generalized thermoelasticity in consideration of fractional order heat conduction equation.

(1)The phenomenon of finite speeds of propagation is observed in all depicted figures. This is expected since the thermal wave travels with finite speed.(2)The effects of the fractional parameter on all the studied fields are very significant.(3)The value of has an essential role in changing the value of the distributions.(4)Significant differences in the physical quantities are observed between the one-temperature theory and the two-temperature theory. The two-temperature theory is more realistic than the one-temperature theory in the case of generalized thermoelasticity.(5)Initially (at ) the conductive temperature has the value of and strain has the value of . These results agree with the boundary conditions.

Nomenclature

: | Lamé’s constant
: | Density
: | Specific heat at constant strain
: | Time
: | Parameter of the Riemann-Liouville fractional integral
: | Conductive temperature
: | Thermodynamic temperature
: | Coefficient of linear thermal expansion
: | Components of stress tensor
: | Components of strain tensor
: | Components of displacement vector
: | Thermal conductivity
: | Relaxation time
: | Longitudinal wave speed
: | Thermal diffusivity
: | Two-temperature parameter
: | Dimensionless two-temperature parameter
: | Dimensionless mechanical coupling constant
: | Heat flux vector
: | Poisson’s ratio.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to Professor S. C. Bose of the Department of Applied Mathematics, University of Calcutta, for his valuable suggestions and guidance in preparation of the paper. The authors also express their sincere thanks to the reviewer for providing valuable suggestions to improve the paper.

References

H. W. Lord and Y. A. Shulman, “A generalized dynamical theory of thermoelasticity,” Journal of the Mechanics and Physics of Solids, vol. 15, no. 5, pp. 299–309, 1967.

A. E. Green and K. A. Lindsay, “Thermoelasticity,” Journal of Elasticity, vol. 2, no. 1, pp. 1–7, 1972.

A. E. Green and P. M. Naghdi, “A re-examination of the basic postulates of thermomechanics,” Proceedings of Royal Society of London, vol. 432, no. 1885, pp. 171–194, 1991.

A. E. Green and P. M. Naghdi, “On undamped heat waves in an elastic solid,” Journal of Thermal Stresses, vol. 15, no. 2, pp. 253–264, 1992.

A. E. Green and P. M. Naghdi, “Thermoelasticity without energy dissipation,” Journal of Elasticity: The Physical and Mathematical Science of Solids, vol. 31, no. 3, pp. 189–208, 1993.

D. Y. Tzou, “Unified field approach for heat conduction from macro- to micro-scales,” Journal of Heat Transfer, vol. 117, no. 1, pp. 8–16, 1995.

S. K. Roychoudhuri, “One-dimensional thermoelastic waves in elastic half-space with dual phase-lag effects,” Journal of Mechanics of Materials and Structures, vol. 2, no. 3, pp. 489–503, 2007.

R. Quintanilla, “Exponential stability in the dual-phase-lag heat conduction theory,” Journal of Non-Equilibrium Thermodynamics, vol. 27, no. 3, pp. 217–227, 2002.

R. Quintanilla, “A condition on the delay parameters in the one-dimensional dual-phase-lag thermoelastic theory,” Journal of Thermal Stresses, vol. 26, no. 7, pp. 713–721, 2003.

L. Wang, M. Xu, and X. Zhou, “Well-posedness and solution structure of dual-phase-lagging heat conduction,” International Journal of Heat and Mass Transfer, vol. 44, no. 9, pp. 1659–1669, 2001.

L. Wang and M. Xu, “Well-posedness of dual-phase-lagging heat conduction equation: higher dimensions,” International Journal of Heat and Mass Transfer, vol. 45, no. 4, pp. 1165–1171, 2002.

L. Wang and X. Mingtian, “Thermal oscillation and resonance in dual-phase-lagging heat conduction,” International Journal of Heat and Mass Transfer, vol. 45, no. 5, pp. 1055–1061, 2002.

P. Ailawalia and S. Budhiraja, “Dynamic problem in thermoelastic solid using dual-phase-lag model with internal heat source,” Journal of Mathematical Sciences and Applications, vol. 2, no. 1, pp. 10–16, 2014.

M. E. Gurtin and W. O. Williams, “On the Clausius-Duhem inequality,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 17, no. 5, pp. 626–633, 1966.

M. E. Gurtin and W. O. Williams, “An axiomatic foundation for continuum thermodynamics,” Archive for Rational Mechanics and Analysis, vol. 26, no. 2, pp. 83–117, 1967.

P. J. Chen and M. E. Gurtin, “On a theory of heat conduction involving two temperatures,” Zeitschrift für Angewandte Mathematik und Physik, vol. 19, no. 4, pp. 614–627, 1968.

P. J. Chen, M. E. Gurtin, and W. O. Williams, “A note on non simple heat conduction,” Zeitschrift für angewandte Mathematik und Physik, vol. 19, no. 6, pp. 969–970, 1968.

P. J. Chen, M. E. Gurtin, and W. O. Williams, “On the thermodynamics of non-simple elastic materials with two temperatures,” Zeitschrift für angewandte Mathematik und Physik, vol. 20, no. 1, pp. 107–112, 1969.

W. E. Warren and P. J. Chen, “Wave propagation in the two temperature theory of thermoelasticity,” Acta Mechanica, vol. 16, no. 1-2, pp. 21–33, 1973.

D. Lesan, “On the thermodynamics of non-simple elastic materials with two temperatures,” Journal of Applied Mathematics and Physics, vol. 21, pp. 583–591, 1970.

P. Puri and P. M. Jordan, “On the propagation of harmonic plane waves under the two-temperature theory,” International Journal of Engineering Science, vol. 44, no. 17, pp. 1113–1126, 2006.

R. Quintanilla, “On existence, structural stability, convergence and spatial behavior in thermoelasticity with two temperatures,” Acta Mechanica, vol. 168, no. 1-2, pp. 61–73, 2004.

S. Mondal, S. H. Mallik, and M. Kanoria, “Two temperature Green Naghdi models with variable thermal conductivity,” Review Bulletin of Calcutta Mathematical Society, vol. 20, no. 2, pp. 121–130, 2012.

S. Mondal, S. H. Mallik, and M. Kanoria, “Two temperature dual-phase-lag model with variable thermal conductivity,” Indian Journal of Theoretical Physics, vol. 61, no. 4, pp. 247–262, 2013.

P. Pal, A. Kar, and M. Kanoria, “Effect of three-phase-lag on two temperature Generalized thermoelasticity for homogeneous isotropic half space,” International Journal of Mathematical Science and Engineerig Applications, vol. 8, no. 1, pp. 81–104, 2014.

M. Islam, A. Kar, and M. Kanoria, “Two-temperature generalized thermoelasticity in a fiber-reinforced hollow cylinder under thermal shock,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 14, no. 5, pp. 367–390, 2013.

P. Das and M. Kanoria, “Two-Temperature Magneto-thermo-elastic response in a perfectly conducting medium based on GN-III model,” International Journal of Pure and Applied Mathematics, vol. 81, no. 4, pp. 199–229, 2012.

S. Banik and M. Kanoria, “Two-temperature generalized thermoelastic interactions in an infinite body with a spherical cavity,” International Journal of Thermophysics, vol. 32, no. 6, pp. 1247–1270, 2011.

R. Kumar, R. Prasad, and S. Mukhopadhyay, “Some theorems on two-temperature generalized thermoelasticity,” Archive of Applied Mechanics, vol. 81, no. 8, pp. 1031–1040, 2011.

R. Kumar, R. Prasad, and S. Mukhopadhyay, “Variational and reciprocal principles in two-temperature generalized thermoelasticity,” Journal of Thermal Stresses, vol. 33, no. 3, pp. 161–171, 2010.

P. Ailawalia, S. Budhiraja, and A. Singla, “Dynamic problem in Green-Naghdi (type III) thermoelastic half-space with two temperature,” Mechanics of Advanced Materials and Structures, vol. 21, no. 7, pp. 544–552, 2014.

G. Honig and U. Hirdes, “A method for the numerical inversion of Laplace transforms,” Journal of Computational and Applied Mathematics, vol. 10, no. 1, pp. 113–132, 1984.

M. Caputo, “Linear models of dissipation whose Q is almost frequently independent II,” Geophys Journal of Royal Astronomical Society, vol. 13, pp. 529–539, 1967.

A. Carpinteri and F. Mainardi, Eds., Fractals and Fractional Calculus in Continuum Mechanics, Springer, New York, NY, USA, 1997.

I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, New York, NY, USA, 1999.

R. Gorenflo and F. Mainardi, “Fractional calculus: integral and differential equations of fractional order,” in Fractals and Fractional Calculus in Continuum Mechanics, Springer, Vienna, Austria, 1997.

R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.

J. Ignaczak and M. Ostoja-Starzeweski, Thermoelasticity with Finite Speed, Oxford University Press, New York, NY, USA, 2010.

L. Debnath and D. Bhatta, Integral Transforms and their Applications, Chapman & Hall/CRC, Taylor & Francis Group, London, UK, 2nd edition, 2007.

H. M. Youssef, “Theory of fractional order generalized thermoelasticity,” Journal of Heat Transfer, vol. 132, no. 6, pp. 1–7, 2010.

F. Mainardi and R. Gorenflo, “On Mittag-Leffler-type functions in fractional evolution processes,” Journal of Computational and Applied Mathematics, vol. 118, no. 1-2, pp. 283–299, 2000.

A. Sur and M. Kanoria, “Fractional order generalized thermoelastic functionally graded solid with variable material properties,” Journal of Solid Mechanics, vol. 6, no. 1, pp. 54–69, 2014.

A. Sur and M. Kanoria, “Fractional heat conduction with finite wave speed in a thermo-visco-elastic spherical shell,” Latin American Journal of Solids and Structures, vol. 11, no. 7, pp. 1132–1162, 2014.

M. A. Ezzat and A. S. El-Karamany, “The uniqueness and reciprocity theorems for generalized thermoviscoelasticity for anisotropic media,” Journal of Thermal Stresses, vol. 25, no. 6, pp. 507–522, 2002.

M. A. Ezzat and A. S. El-Karamany, “On uniqueness and reciprocity theorems for generalized thermoviscoelasticity with thermal relaxation,” Canadian Journal of Physics, vol. 81, no. 6, pp. 823–833, 2003.

G. Jumarie, “Derivation and solutions of some fractional Black-Scholes equations in coarse-grained space and time: application to Merton's optimal portfolio,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1142–1164, 2010.

A. S. El-Karamany and M. A. Ezzat, “Convolutional variational principle, reciprocal and uniqueness theorems in linear fractional two-temperature thermoelasticity,” Journal of Thermal Stresses, vol. 34, no. 3, pp. 264–284, 2011.

A. S. El-Karamany and M. A. Ezzat, “On fractional thermoelasticity,” Mathematics and Mechanics of Solids, vol. 16, no. 3, pp. 334–346, 2011.

A. Sur and M. Kanoria, “Fractional order two-temperature thermoelasticity with finite wave speed,” Acta Mechanica, vol. 223, no. 12, pp. 2685–2701, 2012.

S. Kothari and S. Mukhopadhyay, “A problem on elastic half space under fractional order theory of thermoelasticity,” Journal of Thermal Stresses, vol. 34, no. 7, pp. 724–739, 2011.

M. A. Ezzat and M. A. Fayik, “Fractional order theory of thermoelastic diffusion,” Journal of Thermal Stresses, vol. 34, no. 8, pp. 851–872, 2011.

A. E. Abouelregal, “Fractional order generalized thermo-piezoelectric semi-infinite medium with temperature-dependent properties subjected to a ramp-type heating,” Journal of Thermal Stresses, vol. 34, no. 11, pp. 1139–1155, 2011.

F. Hamza, M. Abdou, and A. M. Abd El-Latief, “Generalized fractional thermoelasticity associated with two relaxation times,” Journal of Thermal Stresses, vol. 37, no. 9, pp. 1080–1098, 2014.

A. M. Zenkour and A. E. Abouelregal, “State-space approach for an infinite medium with a spherical cavity based upon two-temperature generalized thermoelasticity theory and fractional heat conduction,” Zeitschrift für angewandte Mathematik und Physik, vol. 65, no. 1, pp. 149–164, 2014.

R. B. Hetnarski and M. R. Eslami, Thermal Stressss—Advanced Theory and Application, Springer, 2008.

M. R. Eslami, R. B. Hetnarski, J. Ignaczak, N. Noda, N. Sumi, and Y. Tanigawa, Theory of Elasticity and Thermal Stress—Explanatiions, Problems and Solutions, Springer Science + Business Media, Dordrecht, The Netherlands, 2013.

J. Ignaczak and M. Ostoja-Starzewski, Thermoelasticity with Finite Wave Speeds, Oxford Science, 2010.

R. Quintanilla, “A well-posed problem for the dual-phase-lag heat conduction,” Journal of Thermal Stresses, vol. 31, no. 3, pp. 260–269, 2008.

H. M. Youssef, “State-space approach on generalized thermoelasticity for an infinite material with a spherical cavity and variable thermal conductivity subjected to ramp-type heating,” The Canadian Applied Mathematics Quarterly, vol. 13, no. 4, pp. 369–390, 2005.

L. Y. Bahar and R. B. Hetnarski, “State space approach to thermoelasticity,” in Proceeding of 6th Candian Congress of Applied Mechanics, pp. 17–18, University of British Columbia, Vancouver, BC, Canada, 1977.

H. M. Youssef and E. A. Al-Lehaibi, “State-space approach of two-temperature generalized thermoelasticity of one-dimensional problem,” International Journal of Solids and Structures, vol. 44, no. 5, pp. 1550–1562, 2007.