1. Introduction
In 1837, Duhamel established the coupled thermoelasticity theory that takes into account the constitutive relationship between the fields of heat and strain. Biot [
1] introduced the model of coupled temperature elasticity, in which the basic equations were built using Fourier’s law, where the many theories of thermoelasticity were defined using the thermodynamics of irreversible processes. Based on this idea, the system governing the heat equation is of the parabolic type, which states that any thermal oscillation in a substance will affect all locations of the body instantly.
Several different models have been proposed by different researchers in the field of thermoelasticity to obtain hyperbolic heat conduction equations that allow finite velocities of heat waves. The first generalization in this context, known as “generalized LS theory”, was suggested by Lord and Shulman [
2]. Green and Lindsey [
3] proposed the GL model, which generalized the constitutive connections of stress and entropy by taking into account two alternative relaxation factors. Green and Naghdi [
4,
5,
6] provided an additional extension of the theory of thermoelasticity as they proposed three new thermoelastic theories for homogeneous materials, termed GN models I, II and III. The dual-phase-delay theory (DPL), which was refined by Tzou [
7,
8,
9], is the next extension of the thermoelastic theory. By incorporating dual-phase-lag into the heat flow and the temperature gradient, Tzou [
7,
8,
9] constructed a foundational equation to describe the delayed performance of heat and mass transfer in different materials and microstructural factors such as phonon-electron interplay and phonon scattering.
Abouelregal [
10,
11,
12,
13,
14] recently developed a set of new mathematical models that describe the heat transfer process in elastic bodies, involving higher-order time derivatives and phase delays (HOPL). These models are an extension of the mechanical frameworks for the heat transfer theories of Green-Naghdi [
5,
6] and Choudhuri [
15] and heat transfer models with three-phase-delay (TPL) as well as Tzou [
7,
8,
9] models. HOPL models provide a broad theoretical model of heat transfer with diverse microstructural interests, allowing researchers working in the field of heat transfer to reliably predict the thermal response of structures using a multiscale model.
The thermoelasticity model with two-temperature (2TT) was developed by Chen and Gurtin [
16] and Chen et al. [
17,
18]. The Clausius–Duhem inequality was modified in this theory by a model based on two different temperatures: conductive and thermodynamic temperatures. The first was caused by a heat process, and the second was caused by a mechanical system that involves placing between the particles and the slabs of elastic material. They argued that the distinction between these three temperatures is proportional to the amount of heat supplied and that, in the absence of heat supply, the two temperatures are equal in a time-independent scenario. There are no differences between the two temperatures in simple substances and vice versa in the second category. The main difference between this theory and the classical one is the thermal dependence. One of the advantages of this model is that it describes the thermodynamic behavior better in thermoelastic problems that involve time-dependent heat sources, as the two temperatures are proportional to the heat source. In the case where the heat source is absent, the two temperatures are equal.
Quintanilla [
19] discovered two-temperature thermoelasticity and described its existence, structural stability, convergence, and spatial behavior. Youssef [
20] continued this concept based on the thermoelastic theory of heat transfer with a relaxation factor. In this regard, Abouelregal [
10] also created a two-temperature modified thermoelastic version with higher-order time derivatives (HOPL) and three distinct phase delays. Ezzat and El-Karamany [
21] developed a state space technique to provide a model of one-dimensional equations of two temperatures, extending magnetothermoelastic theory in a perfect electric conducting medium with two relaxation periods. Mukhopadhyay et al. expanded thermoelastic theory with two temperatures and dual-phase-lag in their paper [
22]. In recent decades, researchers have devoted close attention to the concept of two-temperature thermoelasticity [
23,
24,
25,
26].
The micropolar elasticity theory, also known as the Cosserat elasticity theory or micropolar continuity mechanics, involves local point rotation in addition to the transformation assumed in the conventional model of elasticity, plus couple stress and force per unit area. In the conventional theory of elasticity, where there is no other type of stress, force stress is simply referred to as “stress”. The concept of couple stress can be traced back to Voigt’s early work on elasticity theory. Couple stress theories have recently been developed, utilizing the full range of current continuity mechanics possibilities. Because a drop in the stress concentration factor near holes and fractures is expected, generalized continuum concepts such as Cosserat elasticity are relevant to material performance. This can lead to increased hardness [
27].
In recent years, Eringen’s theory of micropolar elasticity has attracted a lot of attention because of its potential value in examining the deformation characteristics of solids for which the conventional model is insufficient [
28]. The micropolar concept is thought to be especially effective in studying materials made up of bar-like molecules with micro-rotational influences and the ability to sustain body and surface couples. A micropolar continuum is a compendium of linked particles that take the shape of tiny stiff bodies and move in both directions. The stress vector at a place on a body’s surface element entirely defines the force there [
29]. Micropolarity has a substantial impact on all of the domains covered. Micropolarity has a diminishing influence on the magnitudes of all thermo-physical fields investigated.
The linear description of micropolar thermoelasticity was created by including the temperature influence into the concept of micropolar continua. The micropolar thermoelastic theory was established by Nowacki [
30] and Eringen [
31], who integrated thermal properties into the micropolar concept. Tauchert et al. [
32] have developed the basic equations of the mathematical model of micropolar thermoelasticity. By studying the Green–Lindsay model, Dost and Tabarrok [
33] were able to derive the equations for micropolar generalized thermoelasticity. Chandrasekharaiah [
34] established an energy balance equation and a uniqueness theorem for anisotropic materials using a micropolar thermo elasticity model in which constitutive variables are dependent on heat flow. The constitutive foundations of the three-phase-lag theory of micropolar thermoelasticity were developed by El-Karamany and Ezzat [
35]. They constructed a variational principle for a linear micropolar anisotropic and heterogeneous thermoelastic solid by proving the reciprocity and uniqueness theorems. Under the influence of mechanical plate stress, Alharbi et al. [
36] proposed a mathematical model of a thermoelastic magnetic micropolar half-space with temperature-dependent material variables. Several comprehensive works have been presented that include the theory of micropolar thermoelasticity [
37,
38,
39,
40,
41,
42,
43,
44,
45].
The main object of this investigation is to present a modified model of micropolar thermoelasticity with higher-order derivatives and dual-phase delay time. In addition to the proposed model incorporating microstructural influences in the heat transfer process, the macroscopic construction was also taken into account, with the assumption that the phonon–electron responses lead to a delay in the lattice temperature growth on the macroscopic scales. The proposed model was derived by applying the Taylor series expansions of Fourier’s law and the relationship between the two temperatures while maintaining conditions in phase lags and up to suitable higher orders.
The proposed model with phase lags is considered to be an extension of the two-temperature thermoelastic theory with one relaxation time [
20] and a two-temperature model with two phase lags [
22]. Through several previous studies, it has been realized that the concepts of dual-temperature thermoelasticity may be more applicable in real-world settings. As a result, it is expected that as the practical and theoretical study progresses, these generalized ideas of thermal elasticity with higher time derivatives will be revealed to be highly relevant to many technological applications and challenges.
The topic of wave propagation on the surface of an isotropic micropolar semi-space whose boundary is traction-free was investigated using the proposed model. The expressions for thermodynamic temperature, conductive temperature, microrotation, displacements, and thermal stresses are derived. For the purpose of comparison and investigation in the presented model, the distributions of the examined variables were estimated in tables and figures.
2. Mathematical Model and Basic Equations
In this section, field equations and constitutive relationships will be presented in a micropolar solid in the case of a thermal conductivity model with two temperature biphasic delays and higher-order time derivatives. The basic equations as presented by Eringen [
31] in the absence of body forces, body couples, heat supply, and the balanced external force of the body take the following forms [
46,
47,
48]:
The constitutive equations:
Here, is the stress tensor components, represents the displacement components, alludes to the components of the microrotation vector, gives the components of the couple stress tensor, , , , and are material constants, is known as the Kronecker delta function, describes the components of the electrode small stress tensor, shows the components of rotation, denotes the change in temperature above the reference temperature and is the material constant given by .
The equations of motion are given by:
where
is the permutation symbol,
is the mass density and
is micro-inertia.
Now, using the constitutive Equations (1) and (2), we can remove the stresses
and
from the equations of motion (4) and (5) to obtain:
The energy equation can be written as:
where
is specific heat and
is the heat flux vector. The conventional Fourier law can be written as:
where
is the position vector and
is the thermal conductivity.
A non-simple substance, according to Gurtin and Williams [
17], is one in which the stress, energy, entropy, heat flow, and thermodynamic temperature at a given time all depend on the pasts of the deformation gradient, conduction temperature, and that temperature gradient up to that point.
Quintanilla [
49,
50] replaced the traditional Fourier law (9) as a consequence of the two-temperature concept to be in the form:
where
denotes the thermal conductivity and
represents the conductive temperature measured fulfils the relationship [
16,
17,
18]:
where
is the two-temperature factor.
Quintanilla [
20] and Mukhopadhyay et al. [
22] developed the heat conduction with dual-phase-lag, which included a two-temperature model that took into account microstructural influences in the heat transmission process,
where
and
the phase delays of the heat flux and the conductive temperature gradient.
The expansion of the Taylor series for both sides of Equation (12) will be applied separately in the phase delays
, and
until a sufficiently higher order of the
and
terms respectively [
10,
11]:
If Equation (13) is combined with the energy Equation (8), we get the modified equation for higher order thermal conductivity with two temperatures as follows:
The constraints on the thermomechanical parameters of an isotropic object satisfy the following inequalities:
Chirita et al. [
51,
52,
53] show that there are certain limitations to the use of higher orders
and
, such as when
produces an unstable system incapable of describing an actual physical aspect. However, when the approximation orders are less than or equal to four, the compliance with the Second Law of Thermodynamics must be investigated.
4. Problem Formulation
In this section, a half-space area of a micropolar homogeneous thermal material will be considered (see
Figure 1). Initially, the medium is not deformed, uncompressed, and at a constant temperature of
. It was also hypothesized that the boundary of the medium
is traction-free and exposed to a heat source, which would decrease over time and affect a small
bandwidth surrounding the
-axis. To study the problem, we will take the Cartesian coordinate system (
), provided that the origin of the coordinates is on the upper surface of the plane
, and the
-axis usually indicates the depth of the half-space. When considering plane waves in a plane, all particles on a line parallel to the
axis are shifted equally. As a result, all considered fields are only functions of the
,
, and
variables and do not depend on the
-coordinate. On the basis of these assumptions, the components of the displacement and microrotation vector will be in the following forms:
Then in
-
plane, cubical dilatation
can be written as:
For two-dimensional problems, two equations for motion are obtained from Equation (6) and one for momentum balance from Equation (7) and can be expressed as follows:
The higher order heat conduction Equation (14) with two temperatures can be expressed as:
In the
-
plane, the constitutive relations (1) and (2) are:
The above basic equations can be applied to any boundary condition problem. The micropolar thermoelastic semi-area with a traction-free surface exposed to a time-dependent decreasing heat source will be considered that affects a small -wide band inclosing the -axis. Accordingly, the following boundary conditions will be taken into consideration:
The mechanical boundary conditions:
Thermal boundary condition:
where
is a constant and
is the known Heaviside unit step function. This relationship also indicates that the width of the applied thermal mechanical shock is
on the surface of the
-axis half space and its value is zero elsewhere.
6. Normal Mode Solution
To solve equations from (51) to (56), the solutions will be imposed as follows:
where
denotes the complex constant angular frequency and
denotes the wave number along the
-axis. The variables
,
,
,
,
,
,
and
are indefinite amplitude functions in the variable
only and are completely separate from the time
and coordinate
.
Substituting Equation (57) in Equations (51)–(56), the following equations can be obtained:
where
By removing
,
,
from Equations (58)–(60) and
,
from Equations (61) and (62), we obtain:
where
The solutions of Equations (64) and (65) fulfill the consistency criteria such that the functions
,
,
,
,
,
,
and
trend to zero as
goes to infinity may be expressed as:
where the parameters
,
,
and
,
are some parameters depending on
and
. The parameters
are the positive roots of the equations:
Substituting Equations (67)–(71) into Equations (58)–(62), the following relations can be obtained:
Introducing relations (74) into Equations (67)–(71), then we have:
To derive the displacement components
and
, we first substitute Equation (57) into Equations (49) and (50) and then use Equations (67) and (68):
After applying the normal mode solution method, the thermal and couple stresses components
and
can be written in the forms:
where
The boundary conditions (30)–(32) can be written after using Equation (57) as:
Using the boundary conditions (88) and (89), we get the following results:
After solving the system of equations above, the constants , , can be set. Thus, the full expressions for the different studied physical fields have been obtained.
7. Discussion of the Numerical Results
We now present some discussion of the numerical results of the physical field variables to highlight the theoretical results and the derived mathematical model in the previous sections. For the purpose of comparisons, an investigation will be conducted on a magnesium crystal material. At standard temperature
, the values of the physical parameters are [
48]:
We may take
as real (
) for small values of time. The other physical values are assumed to be [
48]:
,
,
,
, and
.
The calculations are performed at non-dimensional time , in the surface during the interval , taking into account the real component of the amplitude of the physical field variables shown on the vertical axis. At different vertical positions of the -axis, the distributions of conductive and dynamic temperatures ( and ), thermal stresses (, , and ), the components of the displacement vector ( and ), micro-rotation and the components of the couple stress will be graphically illustrated as a result of changing some effective parameters such as the discrepancy coefficient as well as the higher orders of time derivatives and . Some comparisons were also made between the different models in the presence or absence of the micropolar.
The graphs show that the discrepancy coefficient and the higher orders of time derivatives and have a substantial effect on all the physical fields. It is also shown that, depending on the value of the higher order of the time derivatives as well as the discrepancy coefficient , thermal and mechanical waves achieve a stable state. As a result of this analysis, it is necessary to consider the effect of these parameters and also to distinguish between dynamic and conductivity temperatures.
We remember that Chirita et al. [
51,
52] gave interesting results with respect to some limitations when considering more information about the possibilities available for Taylor series expansion orders
and
. They demonstrated that
or
equivalent models always lead to mechanically unstable systems (see also [
53]). When expansion orders are less than or equal to four, the necessary models can be thermodynamically consistent if appropriate assumptions are made throughout the delay [
54,
55].
In this section, the differences in the distribution of the studied fields with distance
will be studied in the case of the thermoelastic micropolar model and two temperatures with one relaxation time (2MLS), in the case of the two-temperature micropolar thermoelastic model with dual-phase lag (2MDPL), and in the case of the one/two-temperature thermoelastic dual-phase lag model and higher orders of time derivatives in the presence and absence of the micropolar (2MHDPL, 2HDPL and 1HDPL), respectively. In fact,
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10 represent the amount of thermodynamic processes that occurred within the material under investigation. We note that the theories of thermoelasticity with two temperatures can be obtained if we ignore the effect of micropolarity, that is, when
. One-temperature thermoelastic theories can also be obtained in the case of
. From the tables and the results obtained in references [
10,
12,
14], it was found that it is sufficient to put
and
to get close to each other’s results.
Figure 2 displays the effect of micropolarity, the higher orders of time derivatives
and
, and the discrepancy coefficient
on the temperature
profile versus
. In the direction of growth of the vertical distance y and in the direction of wave propagation, the dynamic temperature
gradually decreases. It is detected from the figure that, in different theories, the temperature
rises rapidly to its maximum value, which is called the “peak temperature” due to the presence of decreasing thermal load. Then, in a specific area, the capacity of heat gradually decreases until it fades. The figure also confirms that in the absence of micropolarity, the temperature values are much greater than in the case of its presence. As a result, the role of micropolarity in reducing heat wave propagation emerges.
By examining the figure and comparing the curves of the theory 2MHDPL, which includes higher orders of partial derivatives, and the theory of 2MDPL, it is clear that the presence of higher derivatives
and
leads to a limitation of heat spread within the medium. Another important observation that can be extracted from the figure is the role of the two-temperature parameter
in changing the temperature distribution. By comparison, it appears that the temperature curve in the case of the 2HDPL model is much less than that of the curves of the 1HDPL theory. Micropolar thermoelasticity theory will be particularly beneficial in studying materials made up of bar-like molecules with micro-rotational effects and the ability to sustain body and surface couples. From the results obtained and as mentioned in [
33,
34], when analyzing geophysical difficulties, the micropolar elastic model is thought to be more realistic than the traditional elastic model.
As shown in
Figure 3, in all cases where the conductive temperature
begins at a constant value at the free surface of the medium, in full accordance with the thermal boundary conditions. For the generalized theory, the conductive temperature starts at its highest point at the origin (attributed to the diffusion of the thermal boundary) and gradually declines until it attains zero in the direction of the heat wave. When compared with the data for all two-temperature models, the conductive temperature
deals with smaller values in the case of the 1HDP model when the discrepancy coefficient
is omitted. This also shows the discrepancy in the results in the presence and absence of micropolarity. Furthermore, it is clear that higher-order derivatives
and
play an important role in reducing the thermal conductivity distribution. The behavior of micropolar materials depends on the material parameters, phase delay and micro-rotation, and the parameters of the higher time derivatives [
14,
29,
35].
Figure 4 and
Figure 5 exhibit the area variant in normal and transverse displacements
and
in the framework of different micropolar thermoelastic models. We notice that there is a large difference in the distribution of displacements in the presence and absence of a non-dimensional two-temperature parameter
. As shown in
Figure 4 and
Figure 5, the temperature parameter
as well as the higher orders of the partial derivatives of time
and
have a major influence on all displacements
and
. The distribution of the displacements has the same pattern in all three cases that include the micropolar effect (2MLS, 2MDPL, and 2MHDPL) and differs significantly from the two cases in which the micropolar effect disappears (2HDPL, and 1HDPL).
The fluctuations in the vertical and shear stresses
,
and
are illustrated in
Figure 6,
Figure 7 and
Figure 8 under the influence of the discrepancy coefficient
and higher orders
and
, in addition to the micropolar. From the figures, it is clear that the discrepancy coefficient
has a considerable impact on all fields. We can also see that the effect of two temperature coefficient
and higher orders
and
continues to reduce the levels of normal stresses and shear stresses. Depending on the value of the temperature difference as well as higher orders, the waves approach a steady state. In
Figure 6, in the presence of the micropolar effect, the behavior of the contrast in the normal stress
is observed to be similar over the entire scale, with differences in the degree of contrast. The figure also shows a reflection of the curves in the presence and absence of the micropolar, which highlights the prominent effect of the presence of the micropolar on the distribution of vertical stress
. Alharbi et al. [
36] investigated the influence of mechanical strip stress on a mathematical model for a magneto-thermoelastic micropolar medium. They also demonstrated how the micropolar material constants influence stress factors.
It can be seen from
Figure 7 and
Figure 8 that the absolute value of the tangential stress
in the case of micropolar theories (2MLS, 2MDPL, and 2MHDPL) is very small compared to the absolute values in the absence of micropolar (2HDPL, and 1HDPL) and vice versa in the case of shear stress
. From
Figure 7 and
Figure 8, it can be observed that the tangential thermal stress
and the shear stress
in the case of different theories fulfill the mechanical boundary conditions imposed on the problem at the free surface of the medium.
It can be seen from
Figure 7 and
Figure 8 that the absolute value of the tangential stress
in the case of micropolar theories (2MLS, 2MDPL and 2MHDPL) is very small compared to the absolute values in the absence of micropolar (2HDPL and 1HDPL) and vice versa in the case of shear stress
. From
Figure 7 and
Figure 8 it can be observed that the tangential thermal stress
and shear stress
, in the case of different theories, fulfill the mechanical boundary conditions imposed on the problem at the free surface of the medium. Because of the anisotropy and micropolarity, the values of the normal displacement, normal force stress and temperature distribution are significantly variable near the point of application of the heat source. This phenomenon is consistent with previous results, as in [
56].
Figure 9 and
Figure 10 show the variance of micro-rotation
and tangential couple stresses
in the framework of the three theories that include the micropolar effect (2MLS, 2MDPL and 2MHDPL), and the two theories in which the micropolar effect is neglected (2HDPL and 1HDPL). A large difference in the values was observed due to the presence of the non-dimensional discrepancy coefficient
as well as the higher orders of the partial derivatives
and
.
8. Conclusions
In this article, a new heat transfer model with two temperatures is presented in the field of generalized micropolar thermoelasticity. When we distinguish the concept of two temperatures, the first (dynamic temperature) comes from the mechanical process and the second comes from the thermal process (conductive temperature). In addition to the above, the heat conduction equation includes higher orders for the partial derivatives of time. From the proposed model, many previous models can be obtained as special cases, whether in the presence or absence of the micropolar effects, as long as the discrepancy coefficient is neglected. The suggested thermoelastic model was used to examine the behavior of thermal stresses, temperatures, displacements, and micro-rotation in an isotropic, homogeneous, micropolar, thermoelastic half-space by means of normal mode analysis. Calculations and discussions showed the following conclusions:
The results indicate that the discrepancy coefficient has a substantial influence on the thermoelastic distributions within the medium, but the effect on the displacements and thermal stress disturbances is very clear;
There is a large discrepancy in the results between the cases of theories that include two-temperature and those that include one-temperature. The coefficient of two-temperature works to reduce thermal and mechanical waves. Thus, the study of elastic bodies in the case of two-temperature theories is more realistic than the generalized thermoelastic theory at one temperature. Thus, the so-called conductive heat wave must be separated from the so-called thermodynamic heat wave;
Higher orders for partial derivatives play a critical role in all distributions of the investigated domain variables. Thermal and mechanical waves are reduced by using higher orders for partial derivatives. Thus, the extended theory with phase lag times and higher order derivatives may be a better option for describing thermoelasticity than the previous generalized theories as well as the traditional ones;
Micropolarity has a substantial impact on all of the domains covered. Except for tangential stress and thermodynamical temperature, micropolarity has a diminishing influence on the magnitudes of all thermo-physical fields investigated.
Finally, when examining the real behavior of some properties of materials in conjunction with the appropriate geometry of the presented model, the proposed problem takes on a whole new meaning. For example, within the earth, precious materials such as oil and liquid-like elements are found in unrefined form, while the rocks and minerals present may be granular in nature. The fields of geomechanics, seismic engineering, soil dynamics, and other fields are also considered to have practical applications for the specialization of thermal elasticity and waves.