# Effect of Rotation for Two-Temperature Generalized Thermoelasticity of Two-Dimensional under Thermal Shock Problem

Authors: Kh. Lotfy; Wafaa Hassan
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The theory of two-temperature generalized thermoelasticity based on the theory of Youssef is used to solve boundary value problems of two-dimensional half-space. The governing equations are solved using normal mode method under the purview of the Lord-Şhulman (LS) and the classical dynamical coupled theory (CD). The general solution obtained is applied to a specific problem of a half-space subjected to one type of heating, the thermal shock type. We study the influence of rotation on the total deformation of thermoelastic half-space and the interaction with each other under the influence of two temperature theory. The material is homogeneous isotropic elastic half-space. The methodology applied here is use of the normal mode analysis techniques that are used to solve the resulting nondimensional coupled field equations for the two theories. Numerical results for the displacement components, force stresses, and temperature distribution are presented graphically and discussed. The conductive temperature, the dynamical temperature, the stress, and the strain distributions are shown graphically with some comparisons.


---

## Body

## 1. Introduction

The linear theory of elasticity is of paramount importance in the stress analysis of steel, which is the commonest engineering structural material. To a lesser extent, linear elasticity describes the mechanical behavior of the other common solid materials, for example, concrete, wood, and coal. However, the theory does not apply to the behavior of many of the new synthetic materials of the clastomer and polymer type, for example, polymethyl-methacrylate (Perspex), polyethylene, and polyvinyl chloride. The linear theory of micropolar elasticity is adequate to represent the behavior of such materials. For ultrasonic waves, that is, for the case of elastic vibrations characterized by high frequencies and small wavelengths, the influence of the body microstructure becomes significant; this influence of microstructure results in the development of new type of waves that are not in the classical theory of elasticity. Metals, polymers, composites, solids, rocks, and concrete are typical media with microstructures. More generally, most of the natural and manmade materials including engineering, geological, and biological media possess a microstructure.The classical coupled thermoelasticity theory proposed by Biot [1] with the introduction of the strain-rate term in the Fourier heat conduction equation leads to a parabolic-type heat conduction equation, called the diffusion equation. This theory predicts finite propagation speed for elastic waves but an infinite speed for thermal disturbance. This is physically unrealistic. To overcome such an absurdity, generalized thermoelasticity theories have been propounded by Lord and Shulman [2] as well as Green and Lindsay [3] advocating the existence of finite thermal wave speed in solids.These theories have been developed by introducing one or two relaxation times in the thermoelastic process, either by modifying Fourier’s heat conduction equation or by correcting the energy equation and Neuman-Duhamel relation. According to these generalized theories, heat propagation can be visualized as a wave phenomenon rather than a diffusion one; in the literature, it is usually referred to as the second sound effect. These two theories are structurally different from one another and one cannot be obtained as a particular case of the other. Various problems characterizing these theories have been investigated and has revealed some interesting phenomena. Brief reviews of this topic have been reported by Chandrasekharaiah and Srinath [4] and Chandrasekharaiah and Murthy [5]. The interplay of the Maxwell electromagnetic field with the motion of deformable solids is largely being undertaken by many investigators owing to the possibility of its application to geophysical problems and certain topics in optics and acoustics. Moreover, the earth is subject to its own magnetic field and the material of the earth may be electrically conducting. Thus, the magnetoelastic nature of the earth's material may affect the propagation of waves. Many authors have considered the propagation of electromagnetothermoelastic waves in an electrically and thermally conducting solid. A comprehensive review of the earlier contributions to the subject can be found in the study by Puri [6]. Among the authors who considered the generalized magnetothermoelastic equations are Nayfeh and Nemat-Nasser [7] who studied the propagation of plane waves in a solid under the influence of an electromagnetic field. They have obtained the governing equations in the general case and the solution for some particular cases. Roy Choudhuri and Mukhopdhyay [8] extended these results to rotating media. Ezzat [9] has studied the problem of generation of generalized magnetothermoelastic waves by thermal shock in a perfectly conducting half-space. Ezzat et al. [10] have established the model of two dimensional equations of generalized magnetothermoelasticity. In dealing with classical or generalized thermoelastic problems in most situations, the displacement potential function approach is used. However Bahar and Hetnarski [11, 12] outlined several disadvantages of the potential function approach. These may be summarized in the fact that the boundary and initial conditions of the problem are not related directly to the potential function, as it has no physical meaning explicitly.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. To get rid of these difficulties, Bahar and Hetnarski [13] introduced the state space formulation in thermoelastic problems. This state space approach has been further developed in  Sherief [14]  to include the effect of heat sources. Sherief and Anwar [15] surveyed a two-dimensional thermal shock problem for a semi-infinite piezoelectric rod using state space approach. Youssef and El-bary [16] put forward an analysis for a generalized thermoelastic infinite layer problem under three theories using state space approach. State space formulation to the vibration of gold nanobeam in femtoseconds scale was done by Elsibai and Youssef [17].The theory of heat conduction in a deformable body, formulated by Chen and Gurtin [18] and Chen et al. [19, 20] depends on two different temperatures, the conductive temperature and the thermo dynamical temperature. Chen et al. [21] have suggested that the difference between these two temperatures is proportional to heat supply. In absence of heat supply, these two temperatures are identical for time independent situation. However, for time dependent cases, particularly for problems related to wave propagation, the two temperatures, are in general, different, regardless of heat supply. The two temperature thermoelasticity theory has gained much attention of the researchers in the recent years. The existence, structural stability, convergence, and spatial behaviour of two temperature thermoelasticity have been provided by Quintanilla and Tien [22]. Youssef [23] has developed a new model of generalized thermoelasticity that depends on two temperatures T and  ϕ, where the difference between the two temperatures is proportional to heat supply ϕ¨,ii with a nonnegative constant a (length2).Later on, Youssef and Al-Lehaibi [24], Misra et al. [25], Singh [26], and Singh and Bala [27] investigated various problems on the basis of two temperature thermoelasticity with relaxation times and showed that the obtained results are qualitatively different as compared to those in case of one temperature thermoelasticity.Roy Choudhuri and Debnath [28, 29], Othman [30, 31] Othman and Singh [32], and Othman and Song [33] studied the effect of rotation in a micropolar-generalized thermoelastic and thermoviscoelasticity half-space under different theories. The propagation of plane harmonic waves in a rotating elastic medium without thermal field has been studied. It was shown there that the rotation causes the elastic medium to be dispersive and an isotropic. These problems are based on more realistic elastic model since earth, moon, and other plants have angular velocity.Owing to the mathematical difficulties encountered in two- and three-dimensional multifield coupled generalized heat conduction problems, the problems become too complicated to obtain an analytical solution. Instead of analytical methods, several authors have applied numerical methods such as finite difference method, finite element method and boundary value method. For solving such kind of problems, one can find several two-dimensional works based on the generalized thermoelasticity by using the normal mode analysis in the literature Ezzat and Abd Elall [34], Othman et al. [35, 36], Lotfy and Othman [37], Lotfy et al. [38–40], and Sarkar and Lahiri [41]. Using the normal mode analysis technique, we will get the solution in the Fourier transformed domain actually. To apply the normal mode analysis, we have to assume that all the relations are sufficiently smooth on the real axis such that the normal mode analysis of all these functions exists.  The normal mode analysis [34–40] was used to obtain the exact expression for the temperature distribution, thermal stresses, and the displacement components.The present paper is intended to demonstrate the use of the normal mode analysis in analyzing the propagation of thermoelastic waves in two temperature theory of thermoelasticity under the effect of rotation. Normal mode analysis is employed to obtain the exact solution. Numerical values of the field quantities are calculated by considering illustrative examples and their variations with respect to space coordinate are displayed graphically and discussed under thermal shock problem.


## 2. Formulation of the Problem

In the present paper, the authors consider the problem of a homogeneous, isotropic, elastic half-space(x≥0). The surface of the half-space is subjected initially (t=0)  to a thermal shock that is a function of y and  t. Thus all the quantities considered in this problem will be functions of the time variable t and coordinates  x, y. We can introduce the equations of the problem as follows.The heat conduction equation takes the form Youssef [23]:
(1)Kφ,ii=(∂∂t+τ0∂2∂t2)(ρCET+γT0ui,j).
The constitutive equation takes the form
(2)σij=λekkδij+2μeij-γTδij.
Since the medium is rotating uniformly with an angular velocity Ω_=Ωn_, where n_  is a unit vector representing the direction of the axis of rotation. The displacement equation of motion in the rotating frame of reference has two additional terms, centripetal acceleration, Ω_×(Ω_×u_) due to time-varying motion only and the Coriolis acceleration 2Ω_×u_˙,  where u_  is the dynamic displacement vector.The equations of motion in a rotating frame of reference in the context of generalized thermo elasticity are(3)
ρ
[
u
¨
i
+
{
Ω
_
×
(
Ω
_
×
u
_
)
}
i
+
(
2
Ω
_
×
u
˙
_
)
i
]
=
σ
i
j
,
j
,
(
i
,
j
=
1,2
,
3
)
.
The relation between the heat conduction and the dynamical heat takes the form
(4)φ-T=aφ,ii,
where a>0 two-temperature parameter, Youssef [23].Now, we will suppose elastic and homogenous half-spacex≥0 which obey (1)–(4) and initially quiescent where all the state functions are depend only, on the dimension,  x, y  and the time t.The displacement components for one dimension medium have the form(5)ux=u(x,y,t),uy=v(x,y,t),uz=0.
The strain component takes the form
(6)eij=12(ui,j+uj,i).
The heat conduction equation takes the form
(7)K(∂2φ∂x2+∂2φ∂y2)=(∂∂t+τ0∂2∂t2)ρCET+γT0(∂∂t+τ0∂2∂t2)(∂u∂x+∂v∂y).
The constitutive law equations can be written as
(8)σxx=(2μ+λ)∂u∂x+λ∂v∂y-γT,σyy=(2μ+λ)∂v∂y+λ∂u∂x-γT,σxy=μ(∂u∂y+∂v∂x).
Using the summation convection. From (8), we note that the third equation of motion in (3)  is identically satisfied and first two equations become
(9)ρ(∂2u∂t2-Ω2u+2Ωv˙)=μ∇2u+(μ+λ)∂e∂x-γ∂T∂x,ρ(∂2v∂t2-Ω2v-2Ωu˙)=μ∇2v+(μ+λ)∂e∂y-γ∂T∂y.
The relation between the heat conduction and dynamical heat takes the form
(10)φ-T=a(∂2φ∂x2+∂2φ∂y2).
To transform the above equations into nondimensional forms, we define the following non-dimensional variables:
(11)(x′,y′,u′,v′)=c0η(x,y,u,v),(t′,τ0′,υ0′)=c02η(t,τ0,υ0),(θ′,φ′)=(T,φ)-T0T0,(12)σij′=σij2μ+λ,Ω′=Ωc02η,
where η=(ρCE/K), C22=(μ/ρ) and C02=(2μ+λ)/ρ.Hence, we have (dropping the dashed for convenience)(13)∇2φ-(1+τ0∂∂t)∂θ∂t-ε(1+τ0∂∂t)∂e∂t=0,(14)φ-θ=β(∂2φ∂x2+∂2φ∂y2),
where ε=(γ/ρCE) and β=aη2c02.Assume the scalar potential functionsΠ(x,y,t) and ψ(x,y,t) defined by the relations in the nondimensional form
(15)u=∂Π∂x+∂ψ∂y,v=∂Π∂y-∂ψ∂x.
By using (15) and (12) in (9), we obtain
(16)[∇2+Ω2-∂2∂t2]Π+2Ω∂ψ∂t-a0θ=0,(17)(∇2-a1∂2∂t2+a1Ω2)ψ-2Ωa1∂Π∂t=0,
where
(18)a1=ρC02μ,a0=γT0ρC02.
The heat conduction equation (13) becomes
(19)∇2φ-(1+τ0∂∂t)∂θ∂t-ε(1+τ0∂∂t)∂Π∂t=0.


## 3. Normal Mode Analysis

The solution of the physical variable can be decomposed in terms of normal modes as the following way:(20)[Π,ψ,φ,θ,σij](x,y,t)=[Π*,ψ*,u*(x),φ*(x),θ*(x),σij*(x)]×exp(ωt+iby),
where ω is the (complex) time constant, i is the imaginary, b  be a wave number in the y-direction  and  Π*, ψ*, u*(x), φ*(x), θ*(x), and σij*(x) are the amplitude of the functions. By using   the normal mode defined in (20), (17)–(19) and (14) take the following forms:
(21)[D2-A1]Π*+A0ψ*-A2θ*=0,(22)(D2-A4)ψ*-A5Π=0,(23)[D2-A3]φ*=-β*θ*,(24)(D2-b2)φ-Aθ*-BΠ*=0,
where A0=2Ωω, A=ω(1+ωτ0), B=εA, A1=b2+ω2-Ω2, A2=a0, A3=(βb2+1)/β, β*=1/β, A4=b2+a1(ω2-Ω2), A5=a1A0, and D=d/dx.Eliminatingθ*(x), Π*(x), ψ*(x), and φ*(x) between (21) and (24), we obtain the partial differential equation satisfied by  θ*(x)(25)[D6-ED4+FD2-G]Π*(x)=0,
where A6=(β*b2+A3A)/(β*+A), A7=-B/(β*+A).Since(26)E=A1+A4+A6+A2A7,F=A1A4+A0A5+A6(A1+A4)+A2A7(A3+A4),G=A1A4A6+A2A3A4A7+A0A5A6.
In a similar manner, we get
(27)[D6-ED4+FD2-G](θ*,φ*,ψ*)(x)=0.
The above equation can be factorized
(28)(D2-k12)(D2-k22)(D2-k32)Π*(x)=0,
where kn2  (n=1,2,3) are the roots of the following characteristic equation:
(29)k6-Ek4+Fk2-G=0.
The solution of (28) which is bounded as x→∞ is given by
(30)Π*(x)=∑n=13Mn(b,ω)exp(-knx).
Similarly
(31)θ*(x)=∑n=13Mn′(b,ω)exp(-knx),(32)ψ*(x)=∑n=13Mn′′(b,ω)exp(-knx),(33)φ*(x)=∑n=13Mn′′′(b,ω)exp(-knx).
Since,
(34)u*(x)=DΠ*+ibψ*,(35)v*(x)=ibΠ*-Dψ*,(36)e*(x)=Du*+ibv*.
Using (34) and (35), in order to obtain the amplitude of the displacement components u and v, which are bounded as x→∞, then (34) and (35) become
(37)u*(x)=-∑n=13Mn(b,ω)kne-knxu*(x)=+ib∑n=13Mn′′(b,ω)exp(-knx),v*(x)=ib∑n=13Mn(b,ω)e-knxv*(x)=+∑n=13Mn′′(b,ω)kne-knx,
where Mn, Mn′, Mn′′, and Mn′′′ are some parameters depending on β, b, and ω.Substituting from (30)–(32) into (21)–(24), we have
(38)Mn′(b,ω)=H1nMn(b,ω),n=1,2,3,(39)Mn′′(b,ω)=H2nMn(b,ω),n=1,2,3,(40)Mn′′′(b,ω)=H3nMn(b,ω),n=1,2,3,
where
(41)H1n=A7(kn2-A3)kn2-A6,n=1,2,3,(42)H2n=A5(kn2-A4),n=1,2,3,(43)H3n=A7β*A6-kn2,n=1,2,3.
Thus, we have
(44)θ*(x)=∑n=13H1nMn(b,β*,ω)exp(-knx),ψ*(x)=∑n=13H2nMn(b,β*,ω)exp(-knx),φ*(x)=∑n=13H3nMn(b,β*,ω)exp(-knx).
The stress components can be calculated by using (31) and (37) in (8) as follows
(45)σxx*=∑n=13hnMn(b,ω)exp(-knx),σyy*=∑n=13hn′Mn(b,ω)exp(-knx),σxy*=∑n=13hn′′Mn(b,ω)exp(-knx).
The displacements components can be reduced by using (39) in (37) and  we get
(46)u*(x)=∑n=13(ibH2n-kn)Mn(b,ω)e-knx,v*(x)=∑n=13(ib+knH2n)Mn(b,ω)e-knx,
where
(47)hn=-[kn(ibH2n-kn)-ibλ(ib+knH2n)2μ+λ+γT0H1n2μ+λ],hn′=[ib(knH2n+ib)-λkn(ibH2n-kn)2μ+λ-γT0H1n2μ+λ],hn′′=μ[ib(ibH2n-kn)-kn(ib+knH2n)]2μ+λ.
The normal mode analysis is, in fact, to look for the solution in Fourier transformed domain. Assuming that all the field quantities are sufficiently smooth on the real line such that normal mode analysis of these functions exists.


## 4. Application

### 4.1. Thermal Shock Problem

In order to determine the constantsMn  (n=1,2,3). In the physical problem, we should suppress the positive exponentials that are unbounded at infinity. The constants  M1, M2, and M3 have to be  chosen such that the boundary conditions on the surface at x=0 take the form (1)
thermal boundary conditions that the surface of the half-space subjected to thermal shock(48)θ(0,y,t)=f(0,y,t)(2)
mechanical boundary condition that surface of the half-space is traction-free(49)σxx(0,y,t)=0(3)
mechanical boundary condition that surface of the half-space is traction-free(50)σxy(0,y,t)=0,
where  f(0,y,t) is some given function in y and t.Substituting from the expressions of the considered variables into the above boundary conditions (48)–(50), we obtain the following equations satisfyed by the parameters after some simple manipulations
(51)∑n=13H1nMn(b,β*,ω)=f*(y,t),∑n=13hnMn(b,β,ω)=0,∑n=13hn′′Mn(b,β,ω)=0.
Invoking the boundary conditions (51) at the surface x=0 of the plate, we obtain a system of three equations. After applying the inverse of matrix method, we have the values of the three constants Mj, j=1,2,3. Hence, we obtain the expressions of displacements, temperature distribution, and another physical quantity of the plate.


## 4.1. Thermal Shock Problem

In order to determine the constantsMn  (n=1,2,3). In the physical problem, we should suppress the positive exponentials that are unbounded at infinity. The constants  M1, M2, and M3 have to be  chosen such that the boundary conditions on the surface at x=0 take the form (1)
thermal boundary conditions that the surface of the half-space subjected to thermal shock(48)θ(0,y,t)=f(0,y,t)(2)
mechanical boundary condition that surface of the half-space is traction-free(49)σxx(0,y,t)=0(3)
mechanical boundary condition that surface of the half-space is traction-free(50)σxy(0,y,t)=0,
where  f(0,y,t) is some given function in y and t.Substituting from the expressions of the considered variables into the above boundary conditions (48)–(50), we obtain the following equations satisfyed by the parameters after some simple manipulations
(51)∑n=13H1nMn(b,β*,ω)=f*(y,t),∑n=13hnMn(b,β,ω)=0,∑n=13hn′′Mn(b,β,ω)=0.
Invoking the boundary conditions (51) at the surface x=0 of the plate, we obtain a system of three equations. After applying the inverse of matrix method, we have the values of the three constants Mj, j=1,2,3. Hence, we obtain the expressions of displacements, temperature distribution, and another physical quantity of the plate.


## 5. Numerical Results

In order to analyze the above problem numerically, we now consider a numerical example for which computational results are given. The results depict the variation of temperature, displacement, and stress fields in the context of two theories. To study the effect of rotation and two temperature on wave propagation. The copper material was chosen for the purpose of numerical example. The numerical constants (in SI unit) of the problem were taken as(52)λ=7.59×109N/m2,μ=3.86×1010kg/ms2,ρ=8954kg/m3,τ0=0.02s,α=-1.28×109N/m2,β=0.32×109N/m2,η=8886.73m/s2,ε=0.0168,αt=1.78×10-5K-1,k=386Wm-1K-1,b=1,CE=383.1J/(kgK)T0=293K,f*=1,ω=ω0+iξ,ω0=2,ξ=1.
Since we have ω=ω0+iξ, where i is the imaginary  unit, eωt=eω0t(cosξt+isinξt) and for small value of time, we can take ω=ω0 (real).The computations were carried out fora value of time t=0.1. The numerical technique, outlined above, was used for the distribution of the real part of the thermal temperature θ and ϕ, the displacement u, v, strain, and the stress (σxx,σyy,σxy) distribution for the problem. The field quantities, temperature, displacement components,  and stress components depend not only on space x and time t but also on the thermal relaxation time τ0. Here all the variables are taken in nondimensional forms.In the first group Figures1(a), 1(b), 1(c), 1(d), 1(e), 1(f), 1(g), and  1(h), the graph shows the two curves predicted by different theories of thermoelasticity. In these figures, the solid lines represent the solution in the Coupled theory, the dashed lines represent the solution in the generalized Lord and Şhulman theory. We notice that the results for the temperature, the displacement, and stresses distribution when the relaxation time is including in the heat equation are distinctly different from those when the relaxation time is not mentioned in heat equation, because the thermal waves in the Fourier theory of heat equation travel with an infinite speed of propagation as opposed to finite speed in the non-Fourier case. This demonstrates clearly the difference between the coupled and the theory of thermoelasticity (LS).The thermal temperatureθ, the displacement u, v  the  strain   e distribution, the stresses distribution σyy, σxx, and σxy and the thermal  temperature ϕ  distribution with   β=0.1, and  t=0.1 under CD and LS theories.
(a)
The thermodynamical heat distribution
(b)
The displacement distribution (u)
(c)
The displacement distribution (v)
(d)
The strain distribution (e)
(e)
The stressσxx distribution
(f)
The stressσyy distribution
(g)
The stressσxy distribution
(h)
The conductive heat distributionThe second group Figures2(a), 2(b), 2(c), 2(d), 2(e), 2(f), 2(g), and 2(h) show the comparison between  the thermal temperature θ and ϕ, displacement components u, v,  the force stresses components σxx, σyy, and σxy, the case of different two  values of  rotation, and constant of two temperature parameter (β=1) under LS theory. It should be noted (Figure 2(a)) in this problem. It is clear from the graph that θ sharp  decreases to minimum value at the beginning, where it experiences smooth increases (with maximum positive gradient). Graph lines for both values of rotation show different slopes. In other words, the temperature lines for Ω=0.0 has the highest gradient when compared with that of Ω=0.2 in all ranges. In addition, all lines begin to coincide when the horizontal distance x increases to reach the reference temperature of the solid. These results obey physical reality for the behaviour of copper as a polycrystalline solid. Figure 2(b), the horizontal displacement u, despite the peaks (for different values of rotation), the magnitude of the maximum displacement peak strongly depends on the rotation. It is also clear that the rate of change of u decreases with increasing the rotation. On the other hand, Figure 2(c) shows atonable increase of the vertical displacement v, near the beginning reachs minimum value and then reaching zero value at the infinity (state of particles equilibrium) when Ω=0.0. Figure 2(d) displays a comparison of the strain in two cases which show the different behaviours: when Ω=0.0 and Ω=0.2 we can say that significant difference in the strain is noticed for different values of the rotation. In addition, all lines begin to coincide when the horizontal distance x increases to reach zero at infinity. In Figure 2(e) the horizontal stresses σxx graph lines for both values of rotation show different slopes. In other words, the σxx component line for Ω=0.0 has the highest gradient when compared with that of Ω=0.2. In addition, all lines begin to coincide when the horizontal distance x is increased to reach zero after their relaxations at infinity. Variation of Ω has a serious effect on both magnitudes of mechanical stresses. These trends obey elastic and thermoelastic properties of the solid under investigation. Figure 2(f) shows that the stress component  σyy takes the different behavior. In other words, the σyy component line for Ω=0.0  has the highest gradient when compared with that of Ω=0.2. Figure 2(g) shows that the stress component σxy satisfies the boundary condition; it  sharp decreases in the start and start increases (minimum) in the context of the Ω=0.2 but when Ω=0.0 take the different behaviour. The lines for Ω=0.0 has the highest gradient when compared with that of Ω=0.2. These trends obey elastic and thermoelastic properties of the solid. Figure 2(h) displays the conductive temperature in which we observe the significant difference in the conductive temperature that is noticed for the value of rotation. The conductive temperature begins from the positive values and then decreases to arrive the minimum amplitudes in two cases of Ω, also move in the wave propagation, beyond it falls again to try to retain zero at infinity.The thermal temperatureθ, the displacement u, v the   strain   e distribution, the stresses distribution σyy, σxx, and σxy and the thermal  temperature ϕ  distribution with  different values of  rotation  when β=0.1 and t=0.1.
(a)
The thermodynamical heat distribution
(b)
The displacement distribution (u)
(c)
The displacement distribution (v)
(d)
The strain distribution (e)
(e)
The stressσxx distribution
(f)
The stressσyy distribution
(g)
The stressσxy distribution
(h)
The conductive heat distributionThe third group Figures3(a), 3(b), 3(c), 3(d), 3(e), 3(f), 3(g), and 3(h) show the comparison between the thermal temperature θ and ϕ, displacement components (u and v), strain distribution, and the stress (σxx,σyy,σxy) distribution, the case of different two values of two temperature parameter. For the value of y, namely,  y=-1 were substituted in performing the computation. Figure 3(a) exhibits the space variation of temperature distribution in which we observe the following: Significant difference in the thermodynamical temperature is noticed for different value of the nondimensional two-temperature parameter. It should be noted (Figure 3(a)). It is clear from the graph that θ has decreased to arrive the minimum value at the beginning in two cases β=0.1 and β=0.2 (two temperature). The value of temperature quantity converges to zero with increasing the distance x and satisfies the boundary conditions at x=0. Also from this figure we can see when β=0.1 and β=0.2 move in the wave function.The thermal temperatureθ, the displacement u, v the  strain   e distribution, the stresses distribution σyy, σxx, and σxy and the thermal  temperature ϕ  distribution with   different values of two-temperature parameter  at the constants Ω=0.2 and t=0.1.
(a)
The thermodynamical heat distribution
(b)
The displacement distribution (u)
(c)
The displacement distribution (v)
(d)
The strain distribution (e)
(e)
The stressσxx distribution
(f)
The stressσyy distribution
(g)
The stressσxy distribution
(h)
The conductive heat distributionIn Figure3(b) the horizontal displacement u, we see that the displacement component u always starts from the negative value when β=0.1 and β=0.2 and terminates at the zero value begins with increase (then smooth  increases) to reach its maximum magnitude. Beyond it u falls again to try to retain zero at infinity, beyond reaching zero at the infinity (state of particles equilibrium). The displacements u show the same behaviours at different values of β. In Figure 3(c), the vertical displacement v, we see that the displacement component v always starts from the positive value and terminates at the zero value to reach the minimum value, beyond reaching zero at infinity with increases of x. Figure 3(d) displays a comparison of the strain in the context of two cases which show the same behaviours when β=0.1 and β=0.2. We can say that significant difference in the strain is noticed for different values of the non-dimensional two-temperature parameter. In addition, all lines begin to coincide when the horizontal distance x increases to reach zero at infinity. The strain distribution is continuous, smooth, and moves in the wave function. These trends obey elastic and thermoelastic properties of the solid. The stress component, σxx, reach coincidence with zero value (Figure 3(e)) and reaches the maximum value in the beginning and smooth decreases then  converges to zero with increasing the distance x. Figure 3(f) shows that the stress component σyy, increases in the start and arrive to maximum in the context of the two values of β. These trends obey elastic and thermoelastic properties of the solid under investigation. In Figure 3(g) the stress component σxy satisfies the boundary condition and starts from zero. It sharp decreases in the start to arrive the minimum and then start smooth increases to maximum when β=0.1 but sharp increases in the start to arrive the maximum and then start smooth decreases to minimum when β=0.2.  Figure 3(h) displays the conductive temperature in which we observe the significant difference in the conductive temperature that is noticed for the value of the non-dimensional two temperature parameter β where the case of β=0.1 and β=0.2 indicates the new case (two-temperature). The conductive temperature begins from the positive values and then decreases to arrive the minimum amplitudes in two cases of β, also move in the wave propagation when β=0.1 and  β=0.2,  beyond it falls again to try to retain zero at infinity.The forth group Figures4(a), 4(b), 4(c), 4(d), 4(e), 4(f), 4(g), and 4(h) show the comparison between  the temperature  θ, the displacement components u and the force stresses components σxx and σxy, and the case of different two  values of  time, (namely, t=0.1 and  t=0.3) under Lord-Şhulman  (LS) theory. This group shows the effect of time on the results and we found that the curve when t=0.1 is greater than the curves when t=0.2 in all figures. These results obey physical reality for the behaviour of copper as a polycrystalline solid.The thermal temperatureθ, the displacement u, v the  strain   e distribution, the stresses distribution σyy, σxx, and σxy and the thermal temperature ϕ  distribution with different values of time at the constants  Ω=0.2  and  β=0.1.
(a)
The thermodynamical heat distribution
(b)
The displacement distribution (u)
(c)
The displacement distribution (v)
(d)
The strain distribution (e)
(e)
The stressσxx distribution
(f)
The stressσyy distribution
(g)
The stressσxy distribution
(h)
The conductive heat distribution


## 6. Conclusions

(1)
The curves of the physical quantities with (CD) theory in most of figures are lower in comparison with those under (LS) theory, due to the relaxation times.
(2)
Analytical solutions based upon normal mode analysis for thermoelastic problem in solids have been developed and utilized.
(3)
The theory of two-temperature generalized thermoelasticity describes the behavior of the particles of the elastic body more real than the theory of one-temperature generalized thermoelasticity.
(4)
In the context of the theory of two-temperature, the physical functions are continuous.
(5)
The value of all the physical quantities converges to zero with an increase in distancex and all functions are continuous.
(6)
Deformation of a body depends on the nature of forced applied as well as the type of boundary conditions.
(7)
It is clear from all the figures that all the distributions considered have nonzero value only in a bounded region of the half-space. Outside of this  region, the values vanish identically and this means that the region has not  felt thermal disturbance yet.
(8)
From the temperature distributions, we have found  a wave type heat   propagation with finite speeds in the medium. The heat wave front moves forward with a finite speed in the medium with the passage of time which proves that the generalized thermoelasticity theory with two temperature heat transfer is very close to  the behavior of the elastic materials. This is not the case for the CD theory where an infinite speeds of thermal propagation can be found and hence all the considered physical quantities have a nonzero (possibly very small) value for any point in the medium. This indicates that the generalized Fourier’s heat conduction mechanism is completely different from the classical Fourier’s law.



---

*Source: 297274-2013-11-19.xml*