1. Introduction
Waves propagation in a thermoelastic and thermoplastic homogeneous or non-homogeneous media has considerable importance in diverse topics, such as engineering, earthquakes, acoustics, and seismology, due to the non-homogeneities presence in the different layers of the earth’s crust. A lot of materials, e.g., polymers, composites of solids, metals, and rocks, nearly have the same media properties as microstructures.
The strain and thermal fields coupling gave rise to the coupled thermoelasticity theory (CT theory). Duhamel [
1,
2] was the first postulated the originator of the thermal stresses theory who introduced the term of dilatation in equation of heat conduction. Biot [
3] proposed the coupled thermoelasticity theory by inserting the modified Fourier heat conduction equation strain-rate term (diffusion equation) approaching the equation of heat conduction in a parabolic-type, which predicts a propagation with finite speed for the waves in elastic media, although the thermal disturbance has infinite speed. This is physically unrealistic and conflicts with practical results. The theory of Lord-Shulman is the first generalization of the coupled theory. The model basis proposed by Lord and Shulman [
4] modified the Fourier’s law of the equation of heat conduction, presenting a new physical concept (acceleration needs a relaxation time for the heat flow) and denoted by (LS) theory. The heat equation that concerns the LS theory of these types of waves ensures finite speeds of propagation of elastic and heat waves. Lord and Șhulman [
4] and Green and Lindsay [
5] presented generalized thermoelasticity theories advocating the finite wave speed due to the thermal signals in the solids. Chen and Gurtin [
6] and Chen et al. [
7] presented the theory of heat conduction in an elastic body depending on two different temperatures (conductivity temperature and thermodynamic temperature). Puri [
8] pointed out the magneto-thermoelasticity effect on the plane waves’ propagation.
Relaxation times during the process of generalized thermoelastic waves have been inserted to develop these theories made by modifying the conduction heat equation of Fourier′s or reconsidering the Neuman-Duhamel relation (energy equation). Many researchers considered the electromagneto-thermoelastic wave propagation in a thermally conducting and electrically solid. The plane wave propagation in a solid considering the effect of the electro-magnetic field explained by Nayfeh and Nemat-Nasser [
9]. Green and Naghdi [
10,
11,
12] considered the concern (GN, types I, II, and III) theories; the linearized model-I version approaches to the model of thermoelastic in a classical form. Concerns (GN, type II) model, the production rate internal of entropy is identical to zero and implying no thermal energy dissipation. Also, this theory admits that the thermoelastic waves are undamped in a thermoelastic media and denoted by the thermoelasticity model with no energy dissipation. Finally, the GN-III model includes the two previous models as special cases admitting, in general, the dissipation of energy.
Chandrasekharaiah and Srinath [
13] investigated the problem of thermoelastic waves without dissipation of energy due to an external point heat source. The waves’ propagation phenomenon in solids, as a globe of elastic, was first discussed by Bromwich [
14] who considered the gravity effect. Love [
15] considered the effect of gravity parameter on Rayleigh waves and concluded that the Rayleigh waves velocity increased considerably from the gravity field for the large values of wavelengths. The gravity effect on propagation of waves in an elastic layer has been discussed [
16]. Ezzat et al. [
17] presented a magneto-thermoelasticity generalized model in a two-dimensional form. The problem of generalized thermoelasticity considering three theories and the state-space approach has been presented by Youssef and El-bary [
18]. Youssef [
19] formatted a new model of generalized thermoelasticity under two temperature dependence.
The heat conduction model of three-phase-lag has attracted considerable interest. Roy Choudhuri [
20] considered the addition phase lag associated with the variable displacement of temperature to the two relaxation times presented by Tzou (please, see Ref. [
21] and the therein references) that was related to the heat flux vector and the temperature gradient. The above-mentioned displacement variable of thermal is attributable to Green and Naghdi [
10,
11,
12]. Youssef and Al-Lehaibi [
22] pointed out some new problems with two temperature thermoelasticity considering relaxation time. The formulation of the state space vibration in the femtoseconds scale of nano-beam gold was studied by Elsibai and Youssef [
23].
Abd-Alla and Abo-Dahab [
24] pointed out the initial stress, diffusion, and rotation effect on a magneto-thermoelastic generalized problem with a spherical cavity. Abouelregal and Abo-Dahab [
25] discussed the dual-phase-lag (DPL) model on non-homogeneous an infinite spherical cavity solid with a magnetic field. The initial stress and gravitational effect on generalized microstretch-magneto-thermoelastic solid medium for different models were presented by Othman et al. [
26]. Lotfy [
27] investigated interactions in an elastic generalized magneto-thermoelastic medium under two temperatures and three models. Lotfy and Hassan [
28] explored thermal shock and two-temperature problems on the generalized thermoelasticity using normal mode techniques. The theory of two-temperature generalized thermoelasticity considering Youssef’s theory for solving the boundary value problems of two dimensional solid with different types of heating on its boundary was highlighted by Lotfy and Abo-Dahab [
29]. The interaction between the magnetic field, thermal field, and elasticity considering heat transfer fractional derivative for a fibre-reinforced rotation thermoelastic has been investigated by Lotfy and Abo-Dahab [
30]. Abo-Dahab et al. [
31] discussed a generalized two dimensional thermoelasticity considering the magnetic field and rotation.
Recently, Lotfy et al. [
32] have investigated the response of the thermomechanical model for semiconductor medium on a reflection diffusion photothermal waves. Abo-Dahab et al. [
33] discussed the mechanical changes in a microstretch thermoelastic medium considering two-temperature due to pulse heating. Said and Othman [
34] investigated the generalized electro–magneto-thermoelasticity considering three theories under the internal heat source and two temperatures in a finite conducting medium. Othman and Abd-Elaziz [
35] applied three thermoelastic theories considering the micro-temperatures and gravity influence on the porous thermoelastic medium. Carini and Zampoli [
36] studied in detail the three-phase-lag theory in porous matrices assuming three delay times in linearity of thermoelasticity. Marin et al. [
37] discussed some results on the dipolar structure of bodies in the context of Green–Lindsay thermoelasticity. Kumar et al. [
38] discussed the interactions of thermomechanical in magnetothermoelastic and isotropic transversely considering GN type II theory under rotation.
The Earth’s electromagnetic impact through seismic propagation, the machine design of different elements, electromagnetic radiations emissions from plasma physics, nuclear devices, etc. has been discussed by more researchers. Marin and Craciun [
39] investigated the uniqueness concerns the model composite materials boundary value problem in dipolar thermoelasticity.
In this paper, two-dimensional generalized thermoelasticity theory under two temperatures is considered for solve the boundary value problems half-space under the initial stress, gravity, and magnetic field. The governing equations have been solved using new mathematical methods considering Lord-Şhulman (LS), Green-Naghdi theory of type III (GN type III), and three-phase-lag (3PHL) theories to investigate displacement and stresses components and temperature distribution. Comparisons have been made with the predicted results by the three theories; (LS), (GN III), and (3PHL) with the absence and the presence of the magnetic field, initial stress, and gravity. The obtained results were numerically calculated and presented graphically to figure out the physical meaning of the phenomenon. A comparison has been made between the present results in the absence and the presence of the external considered parameters and with the previously obtained results by other researchers.
2. Formulation of the Problem
Considering an isotropic semi-infinite elastic solid, Oxyz is a Cartesian orthogonal coordinate system, any point O of the boundary of the plane, and Oy vertically downward to the medium, as displayed in
Figure 1.
We formulate the fundamental equation of the problem, as follows:
(i) The constitutive equation (stress–strain relation) considering initial stress takes the following form:
(ii) The equation of heat conduction assuming three thermoelastic theories forms is as follows (Kar and Kanoria [
40]):
(iii) The motion equation with body force and heat source absent takes the following form:
(iv) The conductive and thermodynamic temperatures relation takes the form [
19]
Consider that the displacement current is absent; the equations of linearized Maxwell’s that govern the magnetic field for a moving solid slowly having perfect electrical conductivity in the following form
where we used
Using Equation (5), we obtain
The stress of Maxwell produced from the magnetic field can be formed as
which is reduced to
where all notations used have been defined in the “nomenclature” at the end of the paper.
Equation (2) is the generalized thermoelastic solid equation field, applicable to the following:
(LS) theory:
(GN type II) theory:
(3PHL) theory:
The non-dimensional variables take the following form
where
and
.
When we substitute from Equation (8) into Equations (2)–(4), we get
where
The equations of motion approach
where
.
Assuming the scalar potential and vector potential functions
and
Substituting from Equation (13) into Equations (11) and (12), we get
where
The temperature Equation (9) approaches
5. Numerical Results
To display the figure out for the physical meaning of the obtained thee earlier analytical procedure, we assume a numerical example for obtaining computational results in
Table 1. The graphs show the influence of thermoelastic theories, initial stress, electro-magnetic field, and gravity with respect to distance z on the displacement components u and w, components of stress
as well as temperature T.
Additionally, we used the constant values , , , , , s, and .
Considering these physical constants, the MATLAB package is used for calculating the numerical results.
From all the analytical results and numerical calculations via graphs obtained, we obvious that all physical quantities satisfy the boundary conditions at the origin point considered at the starting of the movement of the waves and approach to zero as the axial x approaches to infinity. This indicates that if these quantities vanish, the waves far from the origin point at the wall or the boundary agree with the practical results, and the previous results obtained by others.
First, from
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14,
Figure 15,
Figure 16,
Figure 17,
Figure 18,
Figure 19,
Figure 20,
Figure 21,
Figure 22,
Figure 23,
Figure 24 and
Figure 25, we can discuss the influence of initial stress P on the procedure calculated in the previous obtained results.
Figure 2 and
Figure 3 show that the components of displacement components u and w have a strong effect on neglecting initial stress at the interval (0, 1.15). However, they have a strong effect on the presence of initial stress at the interval (1.15, 1.7). The decrease or increase approaches to zero as the distance z approaches to infinity. It is shown that the normal stresses displayed in
Figure 4 and
Figure 5 behave in the same manner during the intervals (0, 0.3). While they have a strong impact with the initial stress, with the absence of P at (0.3, 0.8) they have a strong effect. After that, the increase or decrease periodically approaches zero as the distance approaches to infinity. While the absence of P has a strong effect on the shear stress
at the interval (0, 0.65) that schematics in
Figure 6, it has a strong effect on (0.65, 1.4) and periodically decreases or increases approach zero as z approaches infinity.
Figure 7 illustrates that Maxwell’s normal stress
has a strong effect on the interval
in the absence of initial stress, taking an inverse behavior for (0.35, 1). After that, it periodically decreases or increases as large values of the distance z approaches zero. The total normal mechanical and Maxwell stresses have a strong behavior in the presence of P for
. In the absence of P, it decreases in (0.3, 0.65) and it is interrupted periodically tending to zero as z approaches to infinity (see,
Figure 8). Finally,
Figure 9 displays the temperature distribution concerns to z considering the absence and presence of P.
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8 and
Figure 9 illustrate the influence of theories with respect to distance z on the displacement components u and w; the components of stress
; as well as the temperature T. From these figures, we noticed the strong effect of Green–Naghdi theory (G–N) compared with the slight effect of Lord–Shulman (L-S), and the three-phase-lag model (3PHL) on u; w;
; and T. We display this effect in detail. In
Figure 2, for
, the displacement component u decreases for z in (1.25, 1.7). Then, u periodically decreases or increases because the large values of z approach zero. The displacement component w behaves the same as u but in a different interval of the distance z, as shown in
Figure 3.
Figure 4 and
Figure 5 show that the components of stress
and
have the same manner concerning z, since they increase for
and decrease for
. After that, they periodically increase or decrease because the large values of z approach zero. However,
, as shown in
Figure 6, decreases when z increases in
, increases in (0.65, 1.4), and periodically decreases or increases approaching zero as z approaches infinity. Likewise,
decreases for
and increases for
. Then,
decreases or increases periodically approaching zero as z approaches infinity (see,
Figure 7).
Figure 8 illustrates that
have an inverse behavior when compared with
. The temperature distribution T decreases for
and increases for
. After that, T decreases or increases periodically approaching zero as z approaches infinity (see
Figure 9).
The influence of the electric field on the displacement components u and w, the components of stress
; and the temperature T concerning distance z appears in
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14,
Figure 15,
Figure 16 and
Figure 17. In
Figure 10, the displacement component u increases with respect to z as the increased of electric field and approaches to zero as z approaches to infinity. On the contrary, w decreases with respect to z as the electric field increases and approaches to zero as z approaches to infinity, as shown in
Figure 11. Moreover,
Figure 12 illustrates that the stress component
decreases with respect to z as the electric field increases and approaches zero as z approaches infinity. However,
decreases for
and increases for
as the electric field increases and approaches zero as z approaches infinity, as seen in
Figure 13. While the electric field increases in
Figure 14, it decreases when
z takes the interval
and approaches zero as z approaches infinity. However,
decreases for
and increases for
as the electric field increases and approaches zero as z approaches infinity, as shown in
Figure 15. Likewise,
decreases for
and increases for
as the electric field increases to zero as z approaches infinity, as shown in
Figure 16. The temperature distribution T, on the contrary to
, increases for
and decreases for
as the electric field increases and approaches zero as z approaches infinity as shown in
Figure 17.
The
Figure 18,
Figure 19,
Figure 20,
Figure 21,
Figure 22,
Figure 23,
Figure 24 and
Figure 25 show the impact of the magnetic field concerning distance z on the components u, w,
,
, and T. In
Figure 18, the component of displacement u increases with respect to z as the magnetic field increases and approaches zero as z approaches infinity. On the contrary, w decreases concerning z as the magnetic field increases and approaches zero as z approaches infinity, as shown in
Figure 19.
Figure 20 illustrates that the stress component
decreases with respect to z as the magnetic field increases and approaches zero as z approaches infinity. However,
decreases for
and increases for
as the magnetic field increases and approaches to zero as z approaches to infinity, see
Figure 21. As the magnetic field increases in
Figure 22,
decreases with respect to z and approaches to zero as z approaches infinity. However,
increases for
and decreases for
as the magnetic field increases and approaches zero as z approaches infinity (see,
Figure 23).
Figure 24 shows that
decreases for
but increases for
as the magnetic field increases and approaches zero as z approaches infinity. As shown in
Figure 25, the temperature distribution T, on the contrary to
increases for
and decreases for
as the magnetic field increases and approaches to zero as z approaches infinity.
Figure 26,
Figure 27,
Figure 28,
Figure 29,
Figure 30,
Figure 31,
Figure 32 and
Figure 33 highlight the influence of the gravity with respect to distance z on the components of displacement u and w; the components of stress
; and the temperature T. In
Figure 26, the displacement component u increases for z as the gravity increases and approaches zero as z approaches infinity.
Figure 27 shows that as the gravity increases, w decreases for
, increases for
, and approaches zero as z approaches infinity. Additionally,
Figure 28 shows that the component of stress
increases concerning z as the gravity increases and approaches zero as z approaches to infinity. However,
increases for
and decreases for
because the gravity increases and approaches zero as z approaches infinity, see
Figure 29. As gravity increases in
Figure 30,
increases for
, decreases for
, and approaches zero as z approaches infinity. Likewise,
increases for
and decreases for
because the gravity increases and approaches zero as z approaches infinity, as displayed in
Figure 31. As the gravity increases,
decrease for
and increase for
and approach zero as z approaches infinity, as shown in
Figure 32. The gravity effect on the temperature distribution T concerning z is shown in
Figure 33. Physically, it is obvious that all dependent components approach zero as the distance approaches infinity. This agrees with the physical meaning of the waves phenomenon if it is far from the origin point of propagation.
Finally, when the initial stress and gravity parameters are vanishes, the obtained results are deduced to the study of Said and Othman [
34] which considered GN (type II) theory and internal heat source as special case from the present study.
6. Conclusions
In this work, we presented an analytical solution based upon the Lame potentials and the normal mode technique for the problem of thermoelastic in a solid medium has been developed, utilized, and compared via graphs. The strong effect of Green-Naghdi theory (GN) compared with the slight effect between Lord–Shulman (LS) and three-phase-lag model (3PHL). The physical quantities converge to zero with an increasing of the distance z and continuously satisfying the assumed boundary conditions. The body deformation depends on the nature of external forces applied (electromagnetic field, two-temperature, initial stress, and gravity), also, the type of thermoelastic theories and boundary conditions. The time parameter, as well as relaxation time, gravity, and electromagnetic field have a strong effect and play a significant strong role in all the physical quantities obtained for the components of stresses, components of displacement, and temperature decreasingly or increasingly. Therefore, the presence of the field of electro-magnetic, gravity, two-temperature, initial stress, and relaxation times in the present model is of significance. The considered method is interesting and applicable to a wide range of phenomena in thermodynamics, thermoelasticity, and magneto-thermoelasticity. The transient behaviors of field variables are studied in detail, and the influences of variations in field variables on each other are discussed. Thus, they provide useful information for practical scientists/ technologists/ researchers/ seismologists/ engineers work in this experimental field on propagation of waves. This paper introduced the effect of gravity, initial stress, electromagnetic field, and two temperatures dependence on the components of displacement, temperature, and components of stress that indicate to their significant effects. Finally, the results provide a significant motivation to study the magneto-thermoelectric conducting materials as a new applicable class of electro-magneto-thermoelectric solids and should prove the useful for the researchers in the material science, designers of new materials, physicists, engineers, and those working on the electro-magneto-thermoelasticity development and in practical situations, especially in optics, geomagnetic, geophysics, acoustics, and oil prospecting.