Journal of Mathematics

Volume 2013 (2013), Article ID 489863, 11 pages

http://dx.doi.org/10.1155/2013/489863

Research Article

The Effect of Magnetic Field and Initial Stress on Fractional Order Generalized Thermoelastic Half-Space

Sunita Deswal, Sandeep Singh Sheoran, and Kapil Kumar Kalkal

Department of Mathematics, G. J. University of Science and Technology, Haryana, Hisar 125001, India

Received 26 December 2012; Revised 5 March 2013; Accepted 12 March 2013

Academic Editor: Petr Ekel

Copyright © 2013 Sunita Deswal 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.

The aim of this paper is to study magneto-thermoelastic interactions in an initially stressed isotropic homogeneous half-space in the context of fractional order theory of generalized thermoelasticity. State space formulation with the Laplace transform technique is used to obtain the general solution, and the resulting formulation is applied to the ramp type increase in thermal load and zero stress. Solutions of the problem in the physical domain are obtained by using a numerical method of the Laplace inverse transform based on the Fourier expansion technique, and the expressions for the displacement, temperature, and stress inside the half-space are obtained. Numerical computations are carried out for a particular material for illustrating the results. Results obtained for the field variables are displayed graphically. Some comparisons have been shown in figures to present the effect of fractional parameter, ramp parameter, magnetic field, and initial stress on the field variables. Some particular cases of special interest have been deduced from the present investigation.

1. Introduction

Biot [1] developed the coupled theory of thermoelasticity to overcome the paradox inherent in the uncoupled theory that elastic changes have no effect on temperature. In this theory, the equations of elasticity and heat conduction are coupled. But, it shares the defect of the uncoupled theory in which it predicts infinite speed of propagation for heat waves. Generalized thermoelastic theories have been developed with the objective of removing the defect of coupled theory. At present, mainly two different models of generalized thermoelasticity are being extensively used, one proposed by Lord and Shulman [2] and the other proposed by Green and Lindsay [3]. The Lord-Shulman theory only modified the Fourier heat conduction equation and suggested one relaxation time, whereas the Green-Lindsay theory modified both the energy equation and the equation of motion and suggested two relaxation times. Dhaliwal and Sherief [4] extended the Lord-Shulman theory to include anisotropic case. Hetnarski and Ignaczak [5] presented a survey article of various representative theories in the range of generalized thermoelasticity. Youssef [6] studied two-dimensional generalized thermoelasticity problem for half-space subjected to ramp type heating.

The theory of magneto-thermoelasticity has received attention of many researchers due to its application in widely diverse fields such as geophysics for understanding the effect of earth's magnetic field on seismic waves, damping of acoustic waves, emission of electromagnetic radiations from nuclear devices, optics, and so forth. The theory of magneto-thermoelasticity was introduced by Knopoff [7] and Chadwick [8] and developed by Kaliski and Petykiewicz [9]. The theoretical outline of the development of magneto-thermoelasticity was discussed by Paria [10]. Paria studied the propagation of plane magneto-thermoelastic waves in an isotropic unbounded medium under the influence of a magnetic field acting transversely to the direction of propagation. Nayfeh and Nemat-Nasser [11] studied the propagation of plane waves in a solid under the influence of electro-magnetic field. Sherief and Ezzat [12] discussed a thermal shock problem in magneto-thermoelasticity with thermal relaxation. Sherief and Helmy [13] illustrated a two-dimensional half-space problem in the context of electromagneto-thermoelasticity theory subjected to a nonuniform thermal shock. Ezzat and Youssef [14] constructed a problem on generalized magneto-thermoelasticity in a perfectly conducting medium. Baksi et al. [15] examined magneto-thermoelastic problem with thermal relaxation and heat source in three-dimensional, infinite rotating elastic medium.

The development of initial stress in the medium is due to many reasons such as the process of quenching, resulting from difference of temperatures, slow process of creep, differential external forces, and gravity variations. The earth is supposed to be under high initial stress. The researchers have shown much interest to study the effect of these stresses on the propagation of waves. Biot [16] solved the dynamic problem of elastic medium under initial stress. Chattopadhyay et al. [17] studied the reflection of elastic waves under initial stress at a free surface. Montanaro [18] studied the isotropic linear thermoelasticity with hydrostatic initial stress by using Biot's linearization of the constitutive law for stress. Othman and Song [19] investigated the reflection of plane waves from an elastic solid half-space under hydrostatic initial stress without energy dissipation. Singh [20] explored the effect of hydrostatic initial stress on waves in a thermoelastic half-space.

The theory of fractional derivative and integral was established in the second half of nineteenth century. The first application of fractional derivative was given by Abel who applied fractional calculus in the solution of an integral equation that arises in the formulation of tautochrone problem. In the recent years, fractional calculus has been applied successfully in various areas to modify many existing models of physical processes such as heat conduction, diffusion, viscoelasticity, wave propagation, and electronics. Caputo and Mainardi [21, 22] and Caputo [23] have established the relation between fractional derivative and theory of linear viscoelasticity. 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 appeared in [24–27]. One can refer to Podlubny [28] for a survey of applications of fractional calculus. Povstenko [29] has proposed a quasistatic uncoupled theory of thermoelasticity based on fractional heat conduction equation. Sherief et al. [30] introduced a new model of thermoelasticity using fractional calculus, proved a uniqueness theorem, and derived a reciprocity relation and a variational principle. In this model, heat conduction equation takes the form as

where are the components of the heat flux vector, is the temperature, is the thermal relaxation time parameter, is thermal conductivity tensor, and is a fractional parameter such that . The above heat conduction equation reduces to the Maxwell-Cattaneo law in the limiting case when . It should be mentioned here that the Maxwell-Cattaneo law has been employed by Lord and Shulman [2] to develop first generalized theory of thermoelasticity. Youssef [31] constructed another model of thermoelasticity in the context of a new consideration of heat conduction with a fractional order and proved the uniqueness theorem. In this model, Youssef described different cases of conductivity: corresponds to weak conductivity, corresponds to normal conductivity, and corresponds to superconductivity. Ezzat [32, 33] established a model of fractional heat conduction equation by using the new Taylor series expansion of time-fractional order developed by Jumarie [34]. El-Karamany and Ezzat [35] 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. The two-temperature dynamic coupled, Lord-Shulman and fractional coupled thermoelasticity theories, result as limit cases. For fractional thermoelasticity not involving two-temperatures, El-Karamany and Ezzat [36] established the uniqueness, reciprocal theorems and convolution variational principle. The dynamic coupled and the Green-Naghdi thermoelasticity theories result as limit cases. The reciprocity relation in case of quiescent initial state is found to be independent of the order of differintegration [35, 36]. Fractional order theory of a perfect conducting thermoelastic medium not involving two temperatures was investigated by El-Karamany and Ezzat [37]. Kothari and Mukhopadhyay [38] studied a half-space problem under fractional order theory of thermoelasticity and analyzed the effect of the fractional order parameter on the field variables.

In the present paper, we study the effect of magnetic field and initial stress under fractional order theory of thermoelasticity proposed by Sherief et al. [30]. We employ a state space approach developed by Bahar and Hetnarski [39] on the formulation. The Laplace transform technique is used to obtain the general solution. The inverse Laplace transform is carried out using a numerical inversion method developed by Honig and Hirdes [40]. Finally, the effect of fractional parameter, ramp parameter, magnetic field, and initial stress on field variables is displayed graphically.

2. Governing Equations

The governing equations in the context of fractional order theory of generalized thermoelasticity with initial stress and magnetic field for isotropic and homogeneous elastic medium are considered as

(i) the equation of motion

 where ,(ii) heat conduction equation

(iii) constitutive relations

We take the linearized Maxwell equations governing the electromagnetic field for a perfectly conducting medium as

where are the components of displacement vector , , is the absolute temperature, is the reference temperature assumed to obey the inequality , is the thermal relaxation time, are the components of the stress tensor, are the components of strain tensor, is the Kronecker delta function, is the cubical dilation, is the density of the medium, is the specific heat, is the thermal conductivity, , is the coefficient of linear thermal expansion, and are the lame constants, is the initial stress, are the components of Lorentz's body force vector , is the magnetic permeability, is the electric permittivity, is the applied magnetic field, is the induced magnetic field, is the induced electric field, and is the current density vector.

3. Problem Formulation

We consider a perfectly conducting isotropic homogeneous and fractional order generalized thermoelastic half-space with hydrostatic initial stress subjected to a constant magnetic field which produces an induced magnetic field and induced electric field . We assume one-dimensional motion for which all the field quantities are functions of and .

The displacement components take the form

The strain component becomes

The components of magnetic field vectors are

The electric intensity vector is parallel to current density vector . Hence, components of and are given as

Now, the Maxwell equation (5) provides the following results:

Using (8) and (10) into the relation , we obtain

The governing equations for one-dimensional case become

Now, we will use the following nondimensional variables:

where

Expressing (12)–(14) in terms of the nondimensional variables given by (15) and dropping the prime sign for convenience, we have the following forms:

where

These equations will be supplemented with appropriate boundary conditions relevant to the particular application under consideration.

Taking the Laplace transform of (17) by using homogeneous initial conditions defined and denoted as

we obtain

Eliminating the value offrom (20) and (21) by using (22), we obtain

where

4. State-Space Formulation

Now, choosing the temperature of heat conduction and the stress component in -direction as state variables, one can write (23) in matrix form as

where

The formal solution of (25) can be written in the form

where

In the above solution, we have cancelled the part of exponential having positive power to get bounded solution for large .

Now, we will use the Cayley-Hamilton theorem to find the form of the matrix , and, for this, we proceed as follows.

The characteristic equation of the matrix is obtained as

where the roots of (29), namely, , and satisfy the following relations:

Now, we write the spectral decomposition of matrix as

where and are called the projectors of and satisfy the following conditions

Then, we have

where

Finally, we get

The characteristic equation of matrix can be written as

where the roots of (36), namely, and , can be written as

Now, the Taylor series expansion of yields

Using the Cayley-Hamilton theorem, we can express and higher orders of matrix in terms of and , where is the identity matrix of second order, Bahar and Hetnarski [39].

Therefore, the infinite series in (38) can be expressed as

where and are coefficients depending on and . By the Cayley-Hamilton theorem, the characteristic roots and of the matrix must satisfy (39), so we have

By solving the above system of equations and using (37), we get

Plugging the values of and in (39), we have

where the entries are given as

The solution of (25) can be written in the following form:

Substituting the values of , , and into (44) and performing the necessary matrix operations, we obtain the values of and as

where

Now, considering (12) along with (15), and the Laplace transform, the displacement component is evaluated as

Substituting (46) into (48), we get

5. Application

Problem: the ramp type boundary temperature of an elastic half-space.

We consider a homogeneous isotropic thermoelastic solid occupying the half-space . The boundary of half-space is affected by ramp type heating. In mathematical notations, the boundary conditions can be denoted as

where is constant temperature and is defined as

where is ramping parameter.

We employ nondimensional variables given in (15) on (50) and taking the Laplace transform, we get

Substituting the values of and from (52) into (45)-(46) and (49), we find

where

6. Limiting Cases

(1)The mathematical expressions for the field variables studied in the context of generalized magneto-thermoelasticity theory with initial stress under applied boundary condition can be obtained by applying in (13).(2)Neglecting initial stress effect, by substituting in (14) and in (58), we obtain the expressions of field variables under fractional order generalized magneto-thermoelasticity.(3)The expressions for studied fields in the context of the fractional order generalized thermoelasticity theory with initial stress can be deduced by setting in (12).

7. Numerical Inversion of the Laplace Transforms

We will now outline the numerical inversion method used to find the solution in the physical domain. The inversion formula of the Laplace transform is defined as

where is the Laplace transform of function .

In order to invert the Laplace transforms in the equations given in Section 5, we apply a numerical inversion method based on the Fourier series expansion explained by Honig and Hirdes [40]. In this method, the inverse transform of the Laplace transform is approximated by the relation as

where is a sufficiently large integer representing the number of terms in the truncated Fourier series chosen such that

where is a prescribed small positive value that corresponds to the degree of accuracy to be achieved and is a positive constant and must be greater than the real parts of all the singularities of . The optimal choice of was obtained according to the criteria described by Honig and Hirdes [40].

8. Numerical Results and Discussions

To illustrate and compare the theoretical results obtained in the Section 5, we now present some numerical results which depict the variations of displacement, temperature, and stress component. The material chosen for the purpose of numerical evaluations is copper, for which we take the following values of the different physical constants

 kgm⁻¹ s⁻²,  Fm⁻¹,  Wm⁻¹ K⁻¹,  K,  kgm⁻³,  s, ,  K⁻¹,  Jkg⁻¹ K⁻¹,  Am⁻¹,  Hm⁻¹.

The general Lame constants and are given as

where is the initial stress parameter, is Young's modulus; and is Poisson ratio. For isotropic elastic medium with no initial stress, we take .

The computations are carried out for , , , and . The numerical technique, outlined in previous section, was used to invert the Laplace transform in (53), providing the displacement, temperature, and stress distributions in the physical domain. The results are represented graphically for different positions of .

Figures 1, 2, and 3 exhibit the space variations of the field quantities in the context of fractional order theory of thermoelasticity with magnetic field and initial stress for different values of fractional parameter .

Figure 1: Displacement distribution for different values of at .

Figure 2: Temperature distribution for different values of at .

Figure 3: Stress distribution for different values of at .

Figure 1 displays the variations of displacement component for different values of , and it is noticed that the magnitude of displacement component decreases with the increase in the value of fractional parameter . In both the cases (i.e., and ), the displacement component attains maximum value at the boundary of half-space, and then continuously decreases to zero. Hence, displacement component has similar trend for both the values of .

Figure 2 depicts the variations of temperature with distance for different values of , and it is noticed that in both the cases (i.e., and ), maximum value of is 2 which is on the boundary of half-space. We observe from the figure that the difference is negligible in the beginning, and with the increase in , the difference is much pronounced up to ; both the series approach to zero. The trends of both the series are alike only up to .

Figure 3 shows the variations of stress component with distance for different values of . It is evident from the figure that both the series have similar trend, that is, first increases to a maximum value and then decreases to a minimum value. The value of increases with the increase of fractional parameter when , while the trend of change is adverse when . The difference is significant in the range .

Figures 4, 5, and 6 depict the variations of displacement component, temperature, and stress component for three different cases defined as (a) ISMT (initially stressed magneto-thermoelasticity), , ; (b) MT (magneto-thermoelasticity), , ; and (c) IST (initially stressed thermoelasticity), , , . In all the above cases (a), (b), and (c), the value of fractional order parameter is taken to be .

Figure 4: Displacement distribution at , .

Figure 5: Temperature distribution at , .

Figure 6: Stress distribution at , .

Figure 4 shows the variations of displacement component under the three cases discussed in (a), (b), and (c). It is observed that the displacement component in case of ISMT and IST has maximum positive values but in case of MT, the displacement component has maximum negative value on the boundary of half-space. The trend for displacement component in case of ISMT and IST is same, that is, both the series are continuously decreasing to zero, but in case of MT, the series rapidly approaches to zero. Figure 4 also exhibits that in the absence of magnetic field (IST) and initial stress (MT) decreases the value of displacement. It is also observed that as compared to ISMT, the difference is significant in case of IST but much pronounced in case of MT.

Figure 5 displays the variations of temperature under the ISMT, MT, and IST theories. It is observed that the trend of the series in case of ISMT, MT, and IST is similar, and the difference is significant. It is found that the absence of initial stress (MT) increases the value of temperature component, but the absence of magnetic field (IST) decreases the value of temperature as compared to the general case (ISMT).

Figure 6 exhibits the variations of stress under ISMT, MT, and IST theories. We found that the behavior of stress component in all the three cases is alike. It also observed that as compared to ISMT theory the stress component has large values in IST theory, but has small values in MT theory. The difference is significant in both IST and MT theories compared with ISMT theory, but the difference is much pronounced in MT theory as compared to IST theory.

Figures 7, 8, and 9 exhibit the variations of displacement component, temperature, and stress component in the context of the fractional order theory of thermoelasticity with magnetic field and initial stress for different values of ramp parameter. The values of ramp parameter are taken as 0.1, 0.3, and 0.5, and the value of fractional order parameter is taken to be for all the three cases.

Figure 7: Displacement distribution for different values of ramp parameter at , .

Figure 8: Temperature distribution for different values of ramp parameter at , .

Figure 9: Stress distribution for different values of ramp parameter at , .

Figure 7 demonstrates that as we increase the value of ramp parameter, then the value of displacement component also increases. It is also observed that trend for all the cases is similar and effect of change in ramp parameter is significant in first two cases and much pronounced in last two cases.

Figure 8 depicts that as we increase the value of ramp parameter, then the value of the temperature rapidly decreases. It is also observed that trend for the series in all the three cases is similar and the effect of change in ramp parameter is much pronounced in the initial range of distance .

Figure 9 displays that as we increase the value of ramp parameter, then the value of the stress component decreases. It is also observed that trend for the series in all the three cases is similar (i.e., all the series have same initial value, and first increases to a maximum value then decreases to a minimum value) and the effect of change in ramp parameter is much pronounced. It is also apparent from the figure that as we increase the value of ramp parameter from 0.1 to 0.3, the values of stress component decrease rapidly as compared to increase in the value of ramp parameter from 0.3 to 0.5.

9. Summary

We consider a perfectly conducting elastic homogeneous half-space in the context of fractional order generalized thermoelasticity theory with magnetic field and initial stress. The method of the matrix exponential, which constitutes the basis of the state space approach of modern theory, is applied to the nondimensional equations. The importance of state space approach is recognized in the fields, where the time behaviour of physical process is of interest. The state space approach is more general than the classical Laplace and Fourier transform technique. Consequently, state space is applicable to all systems that can be analyzed by integral transform in time and also is applicable to many systems for which transform theory breaks down [41]. The potential function approach is often used to solve problems of thermoelasticity theory. This, however, has several disadvantages as outlined in [39]. These may be summarized in the fact that the boundary conditions for physical problems are related directly to the physical quantities under consideration not to the potential functions. Secondly, more stringent assumptions must be made on the behaviour of potential functions than on the actual physical quantities. Last of all, it was found that many integral representations of physical quantities are convergent in the classical sense, while their potential function representations only converge in the mean. All these reasons have led many authors to avoid the use of potential functions. Among the alternatives is the state space formulation. This approach enables one to use the methodology of modern control theory in solving problems of thermoelasticity.

The main conclusions due to the influence of magnetic field, initial stress, fractional parameter, and ramp parameter can be summarized as follows:

(1)the phenomenon of finite speed of propagation is preserved for all the field variables except for stress component due to the presence of initial stress, and all results are in agreement with the generalized theory of thermoelasticity;(2)the fractional parameter has a significant effect on all the studied fields. The displacement component decreases with the increase in the value of the fractional parameter; (3)the thermodynamic temperature first decreases then increases, but the stress component first increases then decreases to a minimum value with the increase in the value of fractional parameter;(4)the effect of magnetic field is silent as compared to initial stress on studied fields. Also, all the field variables for the cases ISMT, MT, and IST behave alike except the displacement component for MT;(5)in the context of fractional order generalized thermoelasticity theory with magnetic field and initial stress, the increase in the value of ramp parameter increases the magnitude of displacement component but decreases the magnitude of temperature distribution and stress component;(6)displacement component , temperature distribution , and stress component show almost similar pattern for different values of fractional parameter and ramp parameter .

Acknowledgment

One of the authors, Sandeep Singh Sheoran, is thankful to U. G. C., New Delhi, for the financial support Vide Letter no. F. 17-11/2008 (SA-1).

References

M. A. Biot, “Thermoelasticity and irreversible thermodynamics,” Journal of Applied Physics, vol. 27, pp. 240–253, 1956.

H. W. Lord and Y. 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.

R. S. Dhaliwal and H. H. Sherief, “Generalized thermoelasticity for anisotropic media,” Quarterly of Applied Mathematics, vol. 38, no. 1, pp. 1–8, 1980.

R. B. Hetnarski and J. Ignaczak, “Generalized thermoelasticity,” Journal of Thermal Stresses, vol. 22, no. 4-5, pp. 451–476, 1999.

H. M. Youssef, “Two-dimensional generalized thermoelasticity problem for a half-space subjected to ramp-type heating,” European Journal of Mechanics A/Solids, vol. 25, no. 5, pp. 745–763, 2006.

L. Knopoff, “The interaction between elastic wave motion and a magnetic field in electrical conductors,” Journal of Geophysical Research, vol. 60, pp. 441–456, 1955.

P. Chadwick, “Elastic wave propagation in a magnetic field,” in Proceedings of the International Congress of Applied Mechanics, pp. 143–153, Brussels, Belgium, 1957.

S. Kaliski and J. Petykiewicz, “Equation of motion coupled with the field of temperature in a magnetic field involving mechanical and electrical relaxation for anisotropic bodies,” Proceedings of Vibration Problems, vol. 4, pp. 1–12, 1959.

G. Paria, “On magneto-thermo-elastic plane waves,” Proceedings of the Cambridge Philosophical Society, vol. 58, pp. 527–531, 1962.

A. H. Nayfeh and S. Nemat-Nasser, “Electromagneto-thermoelastic plane waves in solids with thermal relaxation,” Journal of Applied Mechanics, Transactions ASME, vol. 39, no. 1, pp. 108–113, 1972.

H. H. Sherief and M. A. Ezzat, “A thermal-shock problem in magneto-thermoelasticity with thermal relaxation,” International Journal of Solids and Structures, vol. 33, no. 30, pp. 4449–4459, 1996.

H. H. Sherief and K. A. Helmy, “A two-dimensional problem for a half-space in magneto-thermoelasticity with thermal relaxation,” International Journal of Engineering Science, vol. 40, no. 5, pp. 587–604, 2002.

M. A. Ezzat and H. M. Youssef, “Generalized magneto-thermoelasticity in a perfectly conducting medium,” International Journal of Solids and Structures, vol. 42, no. 24-25, pp. 6319–6334, 2005.

A. Baksi, R. K. Bera, and L. Debnath, “A study of magneto-thermoelastic problems with thermal relaxation and heat sources in a three-dimensional infinite rotating elastic medium,” International Journal of Engineering Science, vol. 43, no. 19-20, pp. 1419–1434, 2005.

M. A. Biot, Mechanics of Incremental Deformation, John Wiley & Sons, New York, NY, USA, 1965.

A. Chattopadhyay, S. Bose, and M. Chakraborty, “Reflection of elastic waves under initial stress at a free surface: P and SV motion,” Journal of the Acoustical Society of America, vol. 72, no. 1, pp. 255–263, 1982.

A. Montanaro, “On singular surfaces in isotropic linear thermoelasticity with initial stress,” Journal of the Acoustical Society of America, vol. 106, no. 3, pp. 1586–1588, 1999.

M. I. A. Othman and Y. Song, “Reflection of plane waves from an elastic solid half-space under hydrostatic initial stress without energy dissipation,” International Journal of Solids and Structures, vol. 44, no. 17, pp. 5651–5664, 2007.

B. Singh, “Effect of hydrostatic initial stresses on waves in a thermoelastic solid half-space,” Applied Mathematics and Computation, vol. 198, no. 2, pp. 494–505, 2008.

M. Caputo and F. Mainardi, “A new dissipation model based on memory mechanism,” Pure and Applied Geophysics, vol. 91, no. 1, pp. 134–147, 1971.

M. Caputo and F. Mainardi, “Linear models of dissipation in anelastic solids,” La Rivista del Nuovo Cimento, vol. 1, no. 2, pp. 161–198, 1971.

M. Caputo, “Vibrations of an infinite viscoelastic layer with a dissipative memory,” Journal of the Acoustical Society of America, vol. 56, no. 3, pp. 897–904, 1974.

R. Gorenflo and F. Mainardi, “Fractional calculus: integral and differential equations of fractional orders,” in Fractals and Fractional Calculus in Continuum Mechanics, vol. 378, pp. 223–276, Springer, Vienna, Austria, 1997.

R. Hilfer, Application of Fraction Calculus in Physics, World Scientific Publishing, Singapore, 2000.

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

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

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

Y. Z. Povstenko, “Fractional heat conduction equation and associated thermal stress,” Journal of Thermal Stresses, vol. 28, no. 1, pp. 83–102, 2005.

H. H. Sherief, A. M. A. El-Sayed, and A. M. Abd El-Latief, “Fractional order theory of thermoelasticity,” International Journal of Solids and Structures, vol. 47, no. 2, pp. 269–275, 2010.

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

M. A. Ezzat, “Thermoelectric MHD non-Newtonian fluid with fractional derivative heat transfer,” Physica B, vol. 405, no. 19, pp. 4188–4194, 2010.

M. A. Ezzat, “Magneto-thermoelasticity with thermoelectric properties and fractional derivative heat transfer,” Physica B, vol. 406, no. 1, pp. 30–35, 2011.

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. S. El-Karamany and M. A. Ezzat, “Fractional order theory of a perfect conducting thermoelastic medium,” Canadian Journal of Physics, vol. 89, no. 3, pp. 311–318, 2011.

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.

L. Y. Bahar and R. B. Hetnarski, “State space approach to thermoelasticity,” Journal of Thermal Stresses, vol. 1, pp. 135–147, 1978.

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.

D. Wiberg, Theory and Problems of State Space and Linear Systems. Schaum's Outline Series in Engineering, Mc Graw-Hill, New York, NY, USA, 1971.