Modeling and Numerical Simulation of Material Science
Vol.05 No.04(2015), Article ID:63077,8 pages

Two-Temperature Generalized Thermoelasticity without Energy Dissipation of Infinite Medium with Spherical Cavity Thermally Excited by Time Exponentially Decaying Laser Pulse

Eman A. N. Al-Lehaibi

Mathematics Department, College of Science and Arts―Sharoura, Najran University, Najran, KSA

Copyright © 2015 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 14 December 2015; accepted 24 January 2016; published 27 January 2016


This work is dealing with two-temperature generalized thermoelasticity without energy dissipation infinite medium with spherical cavity when the surface of this cavity is subjected to laser heating pulse. The closed form solutions for the two types of temperature, strain, and the stress distribution due to time exponentially decaying laser pulse are constructed. The Laplace transformation method is employed when deriving the governing equations. The inversion of Laplace transform will be obtained numerically by using the Riemann-sum approximation method. The results have been presented in figures to show the effect of the time exponentially decaying laser pulse and the two temperature parameter on all the studied fields.


Generalized Thermoelasticity, Two-Temperature, Energy Dissipation, Laser Pulse

1. Introduction

The two temperatures theory of thermoelasticity was introduced by Gurtin and Williams [1] , Chen and Gurtin [2] , and Chen et al. [3] , [4] , in which the classical Clausius-Duhem inequality was replaced by another one depending on two temperatures; the conductive temperature and the thermodynamic temperature, the first is due to the thermal processes, and the second is due to the mechanical processes inherent between the particles and the layers of elastic material, this theory was also investigated by Ieşan [5] .

Among the authors who contribute to developing this theory, Quintanilla studied existence, structural stability, convergence and spatial behavior for this theory [6] , Youssef constructed the generalized Fourier’s law to the two-temperature theory of thermoelasticity and proved its uniqueness of solution for homogeneous isotropic material [7] . Puri and Jordan studied the propagation of plane harmonicwaves, recently [8] , Magaña and Quintanilla [9] have studied the uniqueness and growth solutions for the model proposed by Youssef [7] . A new theory of generalized thermoelasticity has been constructed based on two-temperature generalized thermo- elasticity theory for anisotropic and homogeneous body without energy dissipation by Youssef [10] . This new theorem has been constructed in the context of Green and Naghdi model of type II of linear thermoelasticity. Also, a theorem of general uniqueness is proved for two-temperature generalized thermoelasticity without energy dissipation [10] .

The present paper is devoted to a study of the induced temperature and stress fields in aninfinite elastic medium with aspherical cavity under the purview of two-temperature thermoelasticity without energy dissipation. The medium is considered to be an isotropic homogeneous thermoelastic material. The bounding plane surface of the cavity is thermally loaded bytime exponentially decaying laser pulse. An exact solutions of the problem is obtained in Laplace transformdomain, and the inversions of the Laplace transforms have been culculated numerically. The derived formulations are computed numerically for copper, and the results are presented in graphical form.

2. The Governing Equations

We will consider perfectly conducting, elastic, isotropic, and homogeneous medium and the governing equations will be taken in the context of two-temperature generalized thermoelasticity without energy dissipation.

According to Youssef model, the heat conduction equation takes the form [10] :

, (1)

The conduction-dynamical heat equation takes the form [10] :


The equations of motion take the form


The constitutive equations take the form




where Lame’s constants, density, specific heat at constant strain, , coefficient of linear thermal expansion, is the time, is the temperature, is the reference temperature, is the thermodynamical temperature increment such that, is the conductive heat

temperature, are the components of stress tensor, is the cubic dilatation, is the displacement components, is the characteristic of Green-Naghdi theorem, a is non-negative parameter (two- temperature parameter), and Q is the heat source per unit mass.

3. Formulation of the Problem

We will consider perfectly conducting elastic infinite bodies with spherical cavity occupy the region of an isotropic homogeneous medium whose state can be written in terms of the space variable r and the time variable t such that all of the field functions vanish at infinity. We use a spherical system of coordinates. Due to its symmetric geometry, the problem is one-dimensional with all the functions considered depending on the radial distance r and the time t and the displacement. It is assumed that there are nobody forces in the medium and it is initially quiescent.

Thus, the field equations in spherical one-dimensional case can be put as:




The non-Fourier heat transfer equation due to a laser heating pulse decaying exponentially in time can be written as [7] :


where is the power intensity of surface reflection, I0 is laser peak power intensity, is reflection coefficient, laser pulse parameter and is absorption coefficient.




The constitutive equations will take the following forms







We shall use the following non-dimensional variablesfor convenience [10] :

, , , , (14)

where and.

Equation (1) and Equations (4)-(9) assume the form (where the primes are suppressed for simplicity)


By using Equation (13) into Equation (15), we get


also, we have






where, , , , , and.

4. The Solution in the Laplace Transform Domain

We use the Laplace transform of both sides of the last equations defined as:


Hence, we obtain








Eliminating between Equations (23) and (24), we have


where, and.

From Equations (28) and (24), we obtain


From Equations (22) and (28), we have


where, and.

By eliminating between Equations (28) and (30), we obtain


where, and.

By eliminating between Equations (28) and (30), we get



The bounded solutions of the Equations (31) and (32) take the forms




where are the roots of the equation


By using Equations (33) and (34) into Equation (30), we obtain


Hence, we get


To get the constants, we have to apply the boundary conditions on the surface of the cavity when (R is the radius of the cavity). We will consider that the heat conduction and the strain have zero value on the surface of the cavity, which gives that


Thus, the system of the equations on (34) and (37) gives the following linear equations




By solving the above system, we get


Those complete the solutions as following



Substituting from Equations (41) and (42) in (28), (27) and (25) we get, and respectively.

5. Numerical Inversion of the Laplace Transform

To determine the solutions in the time domain, the Riemann-sum approximation method is used to obtain the numerical results. In this method, any function in Laplace domain can be inverted to the time domain as:


where Re is the real part and is imaginary number unit. For faster convergence, numerous numerical experiments have shown that the value of satisfies the relation Tzou [9] .

6. Numerical Results and Discussion

We now consider a numerical example for which computational results are given. For this purpose, copper is taken as the thermoelastic material for which we take the following values of the different physical constants [11] :

, , ,

, , , ,

From the above values, we get the non-dimensional values of the problem as:

, , , ,.

Figures 1-5 represent the conductive temperature distribution, the thermodynamic temperature distribution, the strain distribution, the displacement distribution, and the stress distribution respectively, in the context of one-temperature type (solid lines) and two-temperature type (dashed lines). We can notice that the two-temper- ature parameter has significant effects on all distribution. The material reaches the steady state through the two- temperature type before the one-temperature type. The peak points decrease when we use the two-temperature type.

Figure 1. The conductive temperature distribution.

Figure 2. The thermodynamic temperature distribution.

Figure 3. The strain distribution.

Figure 4. The displacement distribution.

Figure 5. The stress distribution.


I want thank Prof. Hamdy M. Youssef (Mechanics Department, Faculty of Engineering, Umm Al-Qura University, Makkah KSA) for his help and advises to me to complete this work and to choose this respected journal.

Cite this paper

Eman A. N.Al-Lehaibi, (2015) Two-Temperature Generalized Thermoelasticity without Energy Dissipation of Infinite Medium with Spherical Cavity Thermally Excited by Time Exponentially Decaying Laser Pulse. Modeling and Numerical Simulation of Material Science,05,55-62. doi: 10.4236/mnsms.2015.54006


  1. 1. Gurtin, M.E. and Williams, W.O. (1967) An Axiomatic Foundation for Continuum Thermodynamics. Archive for Rational Mechanics and Analysis, 26, 83-117.

  2. 2. Chen, P.J. and Gurtin, M.E. (1968) On a Theory of Heat Conduction Involving Two Temperatures. Zeitschrift für Angewandte Mathematik und Physik (ZAMP), 19, 614-627.

  3. 3. Chen, P.J., Gurtin, M.E. and Williams, W.O. (1969) On the Thermodynamics of Non-Simple Elastic Materials with Two Temperatures. Zeitschrift für Angewandte Mathematik und Physik (ZAMP), 20, 107-112.

  4. 4. Chen, P.J. and Williams, W.O. (1968) A Note on Non-Simple Heat Conduction. Zeitschrift für Angewandte Mathematik und Physik (ZAMP), 19, 969-970.

  5. 5. Iesan, D. (1970) On the Linear Coupled Thermoelasticity with Two Temperatures. Zeitschrift für Angewandte Mathematik und Physik (ZAMP), 21, 583-591.

  6. 6. Quintanilla, R. (2004) On Existence, Structural Stability, Convergence and Spatial Behavior in Thermoelasticity with Two Temperatures. Acta Mechanica, 168, 61-73.

  7. 7. Youssef, H.M. (2006) Theory of Two-Temperature-Generalized Thermoelasticity. IMA Journal of Applied Mathematics, 71, 383-90.

  8. 8. Puri, P. and Jordan, P.M. (2006) On the Propagation of Harmonic Plane Waves under the Two-Temperature Theory. International Journal of Engineering Science, 44, 1113-1126.

  9. 9. Magaña, A. and Quintanilla, R. (2009) Uniqueness and Growth of Solutions in Two-Temperature Generalized Thermoelastic Theories. Mathematics and Mechanics of Solids, 14, 622-634.

  10. 10. Youssef, H.M. (2011) Theory of Two-Temperature Thermoelasticity without Energy Dissipation. Journal of Thermal Stresses, 34, 138-146.

  11. 11. Youssef, H.M. and Al-Harby, A.H. (2007) State-Space Approach of Two-Temperature Generalized Thermoelasticity of Infinite Body with a Spherical Cavity Subjected to Different Types of Thermal Loading. Archive of Applied Mechanics, 77, 675-687.