Next Article in Journal
Algorithms for Multi-Customer Scheduling with Outsourcing
Next Article in Special Issue
Generalized Thermoelastic Interaction in a Half-Space under a Nonlocal Thermoelastic Model
Previous Article in Journal
A Safe and Efficient Lane Change Decision-Making Strategy of Autonomous Driving Based on Deep Reinforcement Learning
Previous Article in Special Issue
Two-Scale Homogenization of Piezoelectric Perforated Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermoelastic Plane Waves in Materials with a Microstructure Based on Micropolar Thermoelasticity with Two Temperature and Higher Order Time Derivatives

1
Department of Mathematics, College of Science and Arts, Jouf University, Al-Qurayyat 77455, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
3
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500036 Brașov, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(9), 1552; https://doi.org/10.3390/math10091552
Submission received: 30 March 2022 / Revised: 25 April 2022 / Accepted: 26 April 2022 / Published: 5 May 2022
(This article belongs to the Special Issue Applied Mathematics and Continuum Mechanics)

Abstract

:
The study of the effect of the microstructure is important and is most evident in elastic vibrations of high frequency and short-wave duration. In addition to deformation caused by temperature and acting forces, the theory of micropolar thermoelasticity is applied to investigate the microstructure of materials when the vibration of their atoms or molecules is increased. This paper addresses a two-dimensional problem involving a thermoelastic micro-polar half-space with a traction-free surface and a known conductive temperature at the medium surface. The problem is treated in the framework of the concept of two-temperature thermoelasticity with a higher-order time derivative and phase delays, which takes into consideration the impact of microscopic structures in non-simple materials. The normal mode technique was applied to find the analytical formulas for thermal stresses, displacements, micro-rotation, temperature changes, and coupled stress. The numerical results are graphed, and the effect of the discrepancy indicator and higher-order temporal derivatives is examined. There are also some exceptional cases that are covered.
MSC:
74A15; 74J15; 80M05; 80M50

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:
σ i j = λ ε k k δ i j + ( μ + α ) ε i j + μ ε j i 2 α ε i j k ω k γ θ δ i j ,
m i j = ε ω k , k δ i j + ( υ + β ) ω j , i + ( υ β ) ω i , j .
ε i j = 1 2 ( u j , i + u i , j ) .
Here, σ i j is the stress tensor components, u i represents the displacement components, ω k alludes to the components of the microrotation vector, m i j gives the components of the couple stress tensor, λ , μ , α   υ , β and ε are material constants, δ i j is known as the Kronecker delta function, ε i j describes the components of the electrode small stress tensor, ω i shows the components of rotation, θ = T T 0 denotes the change in temperature above the reference temperature T 0 and γ is the material constant given by γ = ( 3 λ + 2 μ + α ) α t .
The equations of motion are given by:
i j
σ i j , j = ρ u .. i ,
ε i j k σ j k + m j i , j = ρ j ω .. , i i , j , k = 1 , 2 , 3 ,
where ε i j k is the permutation symbol, ρ is the mass density and j is micro-inertia.
Now, using the constitutive Equations (1) and (2), we can remove the stresses σ i j and m i j from the equations of motion (4) and (5) to obtain:
( μ + α ) 2 u + ( λ + μ α ) ( . u ) + 2 α × ω γ θ = ρ 2 u t 2 ,
( υ + β ) 2 ω + ( ε + υ β ) ( . ω ) + 2 α × u 4 α ω = ρ j 2 ω t 2 .
The energy equation can be written as:
ρ C E θ t + γ T 0 t ( div u ) = . q + ρ Q ,
where C E is specific heat and q is the heat flux vector. The conventional Fourier law can be written as:
q ( x , t ) = K θ ( x , t ) ,
where x is the position vector and K 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:
q ( x , t ) = K φ ( x , t ) ,
where K denotes the thermal conductivity and φ represents the conductive temperature measured fulfils the relationship [16,17,18]:
φ θ = a 2 φ ,
where a > 0 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,
q ( x , t + τ q ) = K φ ( x , t + τ φ ) ,
where τ q 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 τ q , and τ φ until a sufficiently higher order of the m and n terms respectively [10,11]:
( 1 + r = 1 m τ q r r ! r t r ) q = K ( 1 + r = 1 n τ φ r r ! r t r ) φ .
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:
K ( 1 + r = 1 n τ φ r r ! r t r ) 2 φ = ( 1 + r = 1 m τ q r r ! r t r ) ( ρ C E θ t + γ T 0 t ( div u ) ρ Q ) .
The constraints on the thermomechanical parameters of an isotropic object satisfy the following inequalities:
3 λ + 2 μ > 0 ,   μ > 0 , μ + α > 0 , γ + ε > 0 , α , γ , ε > 0 .
Chirita et al. [51,52,53] show that there are certain limitations to the use of higher orders m and n , such as when m 5 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.

3. Particular Cases

From the modified higher order heat equation with two temperatures (14), it is possible to derive many models of thermoelasticity with two temperatures, both in the existence of the micropolar effect and in the absence of it. One-temperature thermoelastic models can also be derived.
Case 1:
In the absence of the micropolar influence, one/two-temperature thermoelastic theories can be obtained as follows:
  • The conventional thermoelastic theory (CTE) [1] when τ q = τ φ = 0 , a = 0 ,   θ = φ ;
  • The Lord-Shulman generalized thermoelastic model (LS) [2] by setting τ q = τ 0 > 0 ,   a = 0 , θ = φ , τ θ , τ φ 0 and taken m = 1 ;
  • The dual-phase-lag thermoelastic theory (DPL) [8,9] when a = 0 , τ q , τ φ > 0 , n = 1 , and m = 2 ;
  • The dual-phase-lag two-temperature thermoelastic model (2DPL) [22] if a > 0 , τ θ = τ φ , m = n = 1 ;
  • The generalized two-temperature model with one relaxation time (2LS) [20] if a > 0 , τ φ = 0 , τ q = τ 0 > 0 ,   m = n = 1 ;
  • The generalized dual-phase-lag two-temperature thermoelastic model with high-order (2HDPL) [10] is obtained when a > 0 , τ q , τ φ > 0 , n , m 1 ;
  • The generalized dual-phase-lag one-temperature thermoelastic model with high-order (1HDPL) [11] is obtained when a = 0 , τ q , τ φ > 0 , n , m 1 .
Case 2:
In the case of a micropolar effect, one or two-temperature thermoelastic theories can also be obtained. It will be the same as in the previous cases, but the new models will be referred to with the same abbreviations with the addition of the letter “M” to them to become, respectively, 2MLS, 2MDPL, 2MHDPL, MTTLS, and MHTTE.

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 T 0 . It was also hypothesized that the boundary of the medium y = 0 is traction-free and exposed to a heat source, which would decrease over time and affect a small 2 L bandwidth surrounding the x -axis. To study the problem, we will take the Cartesian coordinate system ( x , y , z ), provided that the origin of the coordinates is on the upper surface of the plane y = 0 , and the y -axis usually indicates the depth of the half-space. When considering plane waves in a plane, all particles on a line parallel to the z axis are shifted equally. As a result, all considered fields are only functions of the x , y , and t variables and do not depend on the z -coordinate. On the basis of these assumptions, the components of the displacement and microrotation vector will be in the following forms:
u x = ( x , y , t ) , u y = v ( x , y , t ) , u z = 0 , ω = ( 0 , ω , 0 ) .
Then in x - y plane, cubical dilatation e can be written as:
e = ε k k = div u = u x + v y .
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:
( λ + μ α ) e x + ( μ + α ) 2 u γ θ x + 2 α ω y = ρ 2 u t 2 ,
( λ + μ α ) e y + ( μ + α ) 2 v γ θ y 2 α ω x = ρ 2 v t 2 ,
( υ + β ) 2 ω + 2 α ( v x u y ) 4 α ω = ρ j 2 ω t 2 .
The higher order heat conduction Equation (14) with two temperatures can be expressed as:
K ( 1 + r = 1 n τ φ r r ! r t r ) 2 φ = ( 1 + r = 1 m τ q r r ! r t r ) ( ρ C E θ t + γ T 0 ε k k t ) .
In the x - y plane, the constitutive relations (1) and (2) are:
σ x x = ( λ + 2 μ ) u x + λ v y γ θ ,
σ y y = ( λ + 2 μ ) v y + λ u x γ θ ,
σ z z = λ ( u x + v y ) γ θ ,
σ x y = ( μ + α ) v x + ( μ α ) u y 2 α ω ,
σ y x = ( μ α ) v x + ( μ + α ) u y + 2 α ω ,
m z x = ( υ β ) ω x = ( υ β ) m x z ,
m z y = ( υ β ) ω y = ( υ β ) m y z ,
φ θ = a 2 φ , 2 = 2 x 2 + 2 y 2 .
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 2 L -wide band inclosing the x -axis. Accordingly, the following boundary conditions will be taken into consideration:
  • The mechanical boundary conditions:
    σ y y = σ x y = 0 a t y = 0 ,
    m y z = 0 a t y = 0 .
  • Thermal boundary condition:
    φ ( x , y , t ) = φ 0 H ( L | x | ) exp ( b t ) a t y = 0 ,
    where φ 0 is a constant and H ( . ) is the known Heaviside unit step function. This relationship also indicates that the width of the applied thermal mechanical shock is 2 L on the surface of the x -axis half space and its value is zero elsewhere.

5. Solution Methodology

For convenience, the following dimensionless quantities can be defined:
{ x , y , u , v } = c 1 η { x , y , u , v } , { θ , φ } = γ ρ c 1 2 { θ , φ } , { t , τ q , τ θ } = c 1 2 η { t , τ q , τ θ } σ i j = σ i j α + μ ,   m i j = α   η   m i j c 1 ( α + μ ) ( υ + β ) , ω = α ω α + μ , c 1 2 = ( λ + 2 μ ) ρ ,   η = k ρ C E .
After eliminating primes for convenience in the nomenclature and using dimensionless forms (33), the motion and balance of momentum Equations (18)–(21) and the modified Equation (21) can be expressed as:
β 1 2 2 u x 2 + β 2 2 2 v x y + 2 u y 2 β 1 2 θ x + 2 ω y = β 1 2 2 u t 2 ,
β 1 2 2 v y 2 + β 2 2 2 u x y + 2 v x 2 β 1 2 θ y 2 ω x = β 1 2 2 v t 2 ,
2 ω + g 1 ( v x u y ) g 2 ω = g 3 2 ω t 2 ,
( 1 + r = 1 n τ φ r r ! r t r ) 2 φ = ( 1 + r = 1 m τ q r r ! r t r ) ( θ t + g ε k k t ) ,
where
β 1 2 = λ + 2 μ α + μ , β 2 2 = ( β 1 2 1 ) , g 1 = 2 α 2 η 2 c 1 2 ( α + μ ) ( υ + β ) g 2 = 4 α η 2 c 1 2 ( υ + β ) , g 3 = ρ j c 1 2 υ + β , g = η γ 2 T 0 ρ c 1 2 K .
In the same way, using the non-dimensional constitutive Equations (22)–(29), they become:
σ x x = β 1 2 u x + δ 2 v y β 1 2 θ ,
σ y y = β 1 2 v y + δ 2 u x β 1 2 θ ,
σ z z = δ 2 ( u x + v y ) β 1 2 θ ,
σ x y = v x + δ 3 u y 2 ω ,
σ y x = u x + δ 3 v y + 2 ω ,
m z x = δ 4 ω x = δ 4 m x z ,
m z y = δ 4 ω y = δ 4 m y z ,
φ θ = a 1 2 φ ,
where
δ 2 = λ α + μ , δ 3 = μ α α + μ , δ 4 = υ β υ + β , a 1 = c 1 2 η 2 a .
Due to the large number of differential equations as well as field variables, the potential functions Φ ( x , y , t ) and Ψ ( x , y , t ) can be introduced, which are defined by the relations:
u = Φ + × Ψ , Ψ = ( 0 , 0 , Ψ ) .
Using the above relationship, the components of the displacement u ( x , y , t ) and v ( x , y , t ) can be expressed as:
u = Φ x + Ψ y ,
v = Φ y Ψ x .
Entering (49) and (50) into (34)–(37) yields the following results:
2 ( β 1 2 Φ x + Ψ y ) β 1 2 θ x + 2 ω y = β 1 2 2 t 2 ( Φ x + Ψ y ) ,
2 ( β 1 2 Φ y Ψ x ) β 1 2 θ y 2 ω x = β 1 2 2 t 2 ( Φ y Ψ x ) ,
( 2 g 3 2 t 2 g 2 ) ω = g 1 2 Ψ ,
( 1 + r = 1 n τ φ r r ! r t r ) 2 φ = ( 1 + r = 1 m τ q r r ! r t r ) ( θ t + g 2 t 2 2 Φ ) .
Equations (51) and (52) can be simplified to the following equations.
( 2 2 t 2 ) Φ = θ ,
( 2 β 1 2 2 t 2 ) Ψ = 2 ω .

6. Normal Mode Solution

To solve equations from (51) to (56), the solutions will be imposed as follows:
( θ , φ , u , v , σ i j , Φ , m i j , Ψ , ω ) = ( θ ¯ , φ ¯ , u ¯ , v ¯ , σ ¯ i j , Φ ¯ , m ¯ i j , Ψ ¯ , ω ¯ ) ( y ) exp ( i ζ x + Ω t ) ,
where Ω denotes the complex constant angular frequency and ζ denotes the wave number along the x -axis. The variables θ ¯ , φ ¯ , u ¯ , v ¯ , σ ¯ i j , Φ ¯ , Ψ ¯ and ω ¯ are indefinite amplitude functions in the variable y only and are completely separate from the time t and coordinate x .
Substituting Equation (57) in Equations (51)–(56), the following equations can be obtained:
( D 2 ε 1 ) Φ ¯ = θ ¯
( D 2 ζ 2 ) φ ¯ = ε 4 ( θ ¯ + g ( D 2 ζ 2 ) Φ ¯ )
φ ¯ θ ¯ = a 1 ( D 2 ζ 2 ) φ ¯
( D 2 ε 2 ) Ψ ¯ = 2 ω ¯
( D 2 ε 3 ) ω ¯ = g 1 ( D 2 ζ 2 ) Ψ ¯ ,
where
D = d d y , ε 1 = ζ 2 + Ω 2 , ε 2 = ζ 2 + β 1 2 Ω 2 , ε 3 = ζ 2 + g 2 + g 3 Ω 2 , ε 4 = Ω L q L φ , L q = 1 + r = 1 m Ω r τ q r r ! ,   L φ = 1 + r = 1 n Ω r τ φ r r ! .
By removing θ ¯ , φ ¯ , Φ ¯ from Equations (58)–(60) and ω ¯ , Ψ ¯ from Equations (61) and (62), we obtain:
( D 4 η 1 D 2 + γ 1 ) { θ ¯ , φ ¯ , Φ ¯ } ( y ) = 0 ,
( D 4 η 2 D 2 + γ 2 ) { ω ¯ , Ψ ¯ } ( y ) = 0 ,
where
η 1 = ( ζ 2 + ε 1 + ε 4 [ a 1 ( ε 1 + g ζ 2 ) + ( a 1 ζ 2 + 1 ) ( g + 1 ) ) ] ) ( a 1 ε 4 ( g + 1 ) + 1 ) , η 2 = ε 2 + ε 3 g 1 , γ 1 = ( ζ 2 ε 1 + ε 4 ( a 1 ζ 2 + 1 ) ( ε 1 + g ζ 2 ) ) ( a 1 ε 4 ( g + 1 ) + 1 ) , γ 2 = ε 2 ε 3 2 g 1 ζ 2 .
The solutions of Equations (64) and (65) fulfill the consistency criteria such that the functions θ ¯ , φ ¯ , u ¯ , v ¯ , σ ¯ i j , Φ ¯ , Ψ ¯ and ω ¯ trend to zero as y goes to infinity may be expressed as:
Φ ¯ ( y ) = j = 1 2 A j   e λ j y ,
Ψ ¯ ( y ) = j = 1 2 B j   e μ j y ,
θ ¯ ( y ) = j = 1 2 A j   e λ j y ,
φ ¯ ( y ) = j = 1 2 A j   e λ j y ,
ω ¯ ( y ) = j = 1 2 B j   e μ j y ,
where the parameters A j ,   A j , A j , B j and B j , ( i = 1 , 2 , 3 ) are some parameters depending on ζ and Ω . The parameters λ i , μ j ,   ( i = 1 , 2 ) are the positive roots of the equations:
λ 4 η 1 λ 2 + γ 1 = 0 ,
μ 4 η 2 μ 2 + γ 2 = 0 .
Substituting Equations (67)–(71) into Equations (58)–(62), the following relations can be obtained:
A j = ( λ j 2 ε 1 ) A j , A j = ( λ j 2 ε 1 ) 1 a 1 ( λ j 2 ζ 2 ) A j , B j = 1 2 ( μ j 2 ε 2 ) B j ,   j = 1 , 2 .
Introducing relations (74) into Equations (67)–(71), then we have:
θ ¯ = j = 1 2 ( λ j 2 ε 1 ) A j   e λ j y ,
φ ¯ = j = 1 2 λ j 2 ε 1 1 a 1 ( λ j 2 ζ 2 ) A j   e λ j y ,
ω ¯ = 1 2 j = 1 2 ( μ j 2 ε 2 ) B j   e μ j y .
To derive the displacement components u ¯ and v ¯ , we first substitute Equation (57) into Equations (49) and (50) and then use Equations (67) and (68):
u ¯ = i ζ j = 1 2 λ j A j   e λ j y j = 1 2 λ j B j   e μ j y ,
v ¯ = j = 1 2 λ j A j   e λ j y i ζ j = 1 2 λ j B j   e μ j y .
After applying the normal mode solution method, the thermal and couple stresses components σ ¯ i j and m ¯ i j can be written in the forms:
σ ¯ x x = j = 1 2 α 1 j A j   e λ j y + j = 1 2 β 1 j B j   e μ j y ,
σ ¯ y y = j = 1 2 α 2 j A j   e λ j y + j = 1 2 β 2 j B j   e μ j y ,
σ ¯ z z = j = 1 2 α 3 j A j   e λ j y + j = 1 2 β 3 j B j   e μ j y ,
σ ¯ x y = j = 1 2 α 4 j A j   e λ j y + j = 1 2 β 4 j B j   e μ j y ,
σ ¯ y x = j = 1 2 α 5 j A j   e λ j y + j = 1 2 β 5 j B j   e μ j y ,
m ¯ x z = i ζ 2 j = 1 2 ( μ j 2 ε 2 ) B j   e μ j y ,
m ¯ y z = 1 2 j = 1 2 μ j ( μ j 2 ε 2 ) B j   e μ j y ,
where
α 1 j = ( δ 2 β 1 2 ) λ j 2 + β 1 2 ε 1 ζ 2 β 1 2 , β 1 j = i ζ μ j ( δ 2 δ 1 ) , α 2 j = β 1 2 ε 1 ζ 2 δ 2 , β 2 j = i ζ μ j ( β 1 2 δ 2 ) , α 3 j = β 1 2 ( λ j 2 ε 1 ) ζ 2 δ 2 , β 3 j = i ζ μ j δ 2 ,   α 4 j = 2 i ζ λ j , β 4 j = ζ 2 ε 2 + μ j 2 ( 1 + δ 3 ) , α 5 j = i ζ λ j ( 1 + δ 3 ) , β 5 j = δ 3 ζ 2 + ε 2 .
The boundary conditions (30)–(32) can be written after using Equation (57) as:
σ ¯ y y = 0 , σ ¯ x y = 0 , m ¯ y z = 0 at y = 0 ,
φ ¯ = φ 0 exp ( [ i ζ x + ( Ω + b ) t ] ) = f ¯ ( t , ζ , Ω ) at y = 0 .
Using the boundary conditions (88) and (89), we get the following results:
j = 1 2 α 2 j A j + j = 1 2 β 2 j B j = 0 ,
j = 1 2 α 4 j A j + j = 1 2 β 4 j B j = 0 ,
j = 1 2 μ j ( μ j 2 ε 2 ) B j = 0
j = 1 2 λ j 2 ε 1 1 a 1 ( λ j 2 ζ 2 ) A j = f ¯ ( t , ζ , Ω )
j = 1 2 α y j A j + j = 1 2 β y j A j 3 = 0 ,
j = 1 2 α y x j A j + j = 1 2 β y x j A j 3 = 0 ,
1 2 j = 1 2 E j A j 3 = 0 ,
j = 1 2 E j A j = f ¯ ,
E j = λ j ( λ j 2 ε 2 ) , j = 3 , 4 E j = ( λ j 2 ε 1 ) 1 a 1 ( λ j 2 ζ 2 ) , j = 1 , 2 .
After solving the system of equations above, the constants A j , B j , ( j = 1 , 2 ) 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 T 0 = 298   K , the values of the physical parameters are [48]:
λ = 9.4 × 10 10 ( N / m 2 ) , μ = 4.5 × 10 10 ( N / m 2 ) , α = 0.5 × 10 10 ( N / m 2 ) , υ + β = 0.779 × 10 5 ( N / m 2 ) , K = 1.7 × 10 2 ( J / msK ) , C E = 1.04 × 10 3 ( J / kgK ) , ρ = 1.74 × 10 3 kg / m 3 , j = 0.2 × 10 19 m 2 , γ = 0.779 × 10 9 N .
We may take Ω as real ( Ω = 2 ) for small values of time. The other physical values are assumed to be [48]: ζ = 2 , L = 2 , a 1 = 0.2 , φ 0 = 1 , and b = 1 .
The calculations are performed at non-dimensional time t = 0.2 , in the surface x = 0.1 during the interval 0 y 3 , 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 y -axis, the distributions of conductive and dynamic temperatures ( φ and θ ), thermal stresses ( σ x x , σ y y , and σ y x ), the components of the displacement vector ( u and v ), micro-rotation ω and the components of the couple stress m y z will be graphically illustrated as a result of changing some effective parameters such as the discrepancy coefficient a 1 as well as the higher orders of time derivatives m and n . Some comparisons were also made between the different models in the presence or absence of the micropolar.
The graphs show that the discrepancy coefficient a 1 and the higher orders of time derivatives m and n 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 a 1 , 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 m and n . They demonstrated that m > 5 or n > 5 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 y 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 α = β = α = υ = ε = j = 0 . One-temperature thermoelastic theories can also be obtained in the case of a 1 = 0 . From the tables and the results obtained in references [10,12,14], it was found that it is sufficient to put m = 3 and n = 2 to get close to each other’s results.
Figure 2 displays the effect of micropolarity, the higher orders of time derivatives m and n , and the discrepancy coefficient a 1 on the temperature θ profile versus y . 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 m and n 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 a 1 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 a 1 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 m and n 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 u and v 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 a 1 . As shown in Figure 4 and Figure 5, the temperature parameter a 1 as well as the higher orders of the partial derivatives of time m and n , have a major influence on all displacements u and v . 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 σ x x , σ y y and σ y x are illustrated in Figure 6, Figure 7 and Figure 8 under the influence of the discrepancy coefficient a 1 and higher orders m and n , in addition to the micropolar. From the figures, it is clear that the discrepancy coefficient a 1 has a considerable impact on all fields. We can also see that the effect of two temperature coefficient a 1 and higher orders m and n 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 σ x x 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 σ x x . 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 σ y y 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 σ y x . From Figure 7 and Figure 8, it can be observed that the tangential thermal stress σ y y and the shear stress σ y x 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 σ y y 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 σ y x . From Figure 7 and Figure 8 it can be observed that the tangential thermal stress σ y y and shear stress σ y x , 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 m y z 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 a 1 as well as the higher orders of the partial derivatives m and n .

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.

Author Contributions

Conceptualization: A.E.A., M.M. and F.A.; methodology: A.E.A. and F.A.; validation: A.E.A., F.A. and M.M.; formal analysis: A.E.A., F.A. and M.M.; investigation: A.E.A., M.M. and F.A.; resources: F.A.; data curation: A.E.A., F.A. and M.M.; writing—original draft preparation: A.E.A., F.A. and M.M.; writing—review and editing: A.E.A.; visualization: F.A. and M.M.; supervision: A.E.A., F.A. and M.M.; project administration: A.E.A. All authors have read and agreed to the published version of the manuscript.

Funding

The Deanship of Scientific Research at Jouf University, Saudi Arabia funded this project under grant No. (DSR-2021-03-0376).

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Acknowledgments

The authors extend their appreciation to the core member of Scientific Research at Jouf University for funding this work through research. We would also like to extend our sincere thanks to the College of Science and Arts in Al-Qurayyat for its technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Biot, M. Thermoelasticity and irreversible thermodynamics. J. Appl. Phys. 1956, 27, 240–253. [Google Scholar] [CrossRef]
  2. Lord, H.W.; Shulman, Y.H. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  3. Green, A.E.; Naghdi, P.M. A re-examination of the basic results of thermomechanics. Proc. Math. Phys. Sci. 1991, 432, 171–194. [Google Scholar]
  4. Green, A.E.; Lindsay, K.A. Thermoelasticity. J. Elast. 1972, 2, 1–7. [Google Scholar] [CrossRef]
  5. Green, A.E.; Naghdi, P.M. On undamped heat waves in an elastic solid. J. Therm. Stress. 1992, 15, 252–264. [Google Scholar] [CrossRef]
  6. Green, A.E.; Naghdi, P.M. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  7. Tzou, D.Y. A unified filed approach for heat conduction from macro to macroscales. ASME J. Heat Transf. 1995, 117, 8–16. [Google Scholar] [CrossRef]
  8. Tzou, D.Y. The generalized lagging response in small-scale and high-rate heating. Int. J. Heat Mass Transf. 1995, 38, 3231–3234. [Google Scholar] [CrossRef]
  9. Tzou, D.Y. Experimental support for the lagging behavior in heat propagation. J. Thermophys. Heat Transf. 1995, 9, 686–693. [Google Scholar] [CrossRef]
  10. Abouelregal, A.E. Two-temperature thermoelastic model without energy dissipation including higher order time-derivatives and two phase-lags. Mater. Res. Express 2019, 6, 116535. [Google Scholar] [CrossRef]
  11. Abouelregal, A.E. On Green and Naghdi thermoelasticity model without energy dissipation with higher order time differential and phase-lags. J. Appl. Comput. Mech. 2020, 6, 445–456. [Google Scholar]
  12. Abouelregal, A.E. A novel generalized thermoelasticity with higher-order time-derivatives and three-phase lags. Multidiscip. Model. Mater. Struct. 2019, 16, 689–711. [Google Scholar] [CrossRef]
  13. Abouelregal, A.E. A novel model of nonlocal thermoelasticity with time derivatives of higher order. Math. Methods Appl. Sci. 2020, 43, 6746–6760. [Google Scholar] [CrossRef]
  14. Abouelregal, A.E. Three-phase-lag thermoelastic heat conduction model with higher-order time-fractional derivatives. Indian J. Phys. 2020, 94, 1949–1963. [Google Scholar] [CrossRef]
  15. Choudhuri, S.R. On a thermoelastic three-phase-lag model. J. Therm. Stress. 2007, 30, 231–238. [Google Scholar] [CrossRef]
  16. Chen, P.J.; Gurtin, M.E. On a theory of heat conduction involving two temperatures. Z. Angew. Math. Phys. 1968, 19, 614–627. [Google Scholar] [CrossRef]
  17. Chen, P.J.; Williams, W.O. A note on non-simple heat conduction. Z. Angew. Math. Phys. 1968, 19, 969–970. [Google Scholar] [CrossRef]
  18. Chen, P.J.; Gurtin, M.E.; Williams, W.O. On the thermodynamics of non-simple elastic materials with two temperatures. Z. Angew. Math. Phys. 1969, 20, 107–112. [Google Scholar] [CrossRef]
  19. Quintanilla, R. On existence, structural stability, convergence and spatial behavior in thermoelasticity with two temperatures. Acta Mech. 2004, 168, 61–73. [Google Scholar] [CrossRef]
  20. Youssef, H. Theory of two-temperature-generalized thermoelasticity. IMA J. Appl. Math. 2006, 71, 383–390. [Google Scholar] [CrossRef]
  21. Ezzat, M.A.; El-Karamany, A.S. Two temperature theory in generalized magneto thermoelasticity with two relaxation times. Meccanica 2011, 46, 785–794. [Google Scholar] [CrossRef]
  22. Mukhopadhyay, S.; Prasad, R.; Kumar, R. On the theory of two-temperature thermoelasticity with two phase-lags. J. Therm. Stress. 2011, 34, 352–365. [Google Scholar] [CrossRef]
  23. Mukhopadhyay, S.; Kumar, R. Thermoelastic Interactions on Two-Temperature Generalized Thermoelasticity in an Infinite Medium with a Cylindrical Cavity. J. Therm. Stress. 2009, 32, 341–360. [Google Scholar] [CrossRef]
  24. Fernández, J.R.; Quintanilla, R. Uniqueness and exponential instability in a new two-temperature thermoelastic theory. AIMS Math. 2021, 6, 5440–5451. [Google Scholar] [CrossRef]
  25. Sarkar, N.; Mondal, S. Two-dimensional problem of two-temperature generalized thermoelasticity using memory-dependent heat transfer: An integral transform approach. Indian J. Phys. 2020, 94, 1965–1974. [Google Scholar] [CrossRef]
  26. Hobiny, A.; Alzahrani, F.; Abbas, I.; Marin, M. The effect of fractional time derivative of bioheat model in skin tissue induced to laser irradiation. Symmetry 2020, 12, 602. [Google Scholar] [CrossRef]
  27. Hassanpour, S.; Heppler, G.R. Micropolar elasticity theory: A survey of linear isotropic equations, representative notations, and experimental investigations. Math. Mech. Solids 2017, 22, 224–242. [Google Scholar] [CrossRef]
  28. Eringen, A.C. Linear theory of micropolar elasticity. J. Appl. Math. Mech. 1966, 15, 909–923. [Google Scholar]
  29. Abouelregal, A.E.; Marin, M. The size-dependent thermoelastic vibrations of nanobeams subjected to harmonic excitation and rectified sine wave heating. Mathematics 2020, 8, 1128. [Google Scholar] [CrossRef]
  30. Nowacki, W. Theory of Asymmetric Elasticity; Pergamon Press: Oxford, NY, USA, 1986. [Google Scholar]
  31. Eringen, A.C. Foundations of Micropolar Thermoelasticity; International Centre for Mechanical Science, Udine Course and Lectures 23; Springer: Berlin, Germany, 1970. [Google Scholar]
  32. Tauchert, T.R.; Claus, W.D.; Ariman, T. The linear theory of micropolar thermoelasticity. Int. J. Eng. Sci. 1968, 6, 37–47. [Google Scholar] [CrossRef]
  33. Dost, S.; Tabarrok, B. Generalized micropolar thermoelasticity. Int. J. Eng. Sci. 1978, 16, 173. [Google Scholar] [CrossRef]
  34. Chandrasekhariah, D.S. Heat flux dependent micropolar elasticity. Int. J. Eng. Sci. 1986, 24, 1389–1395. [Google Scholar] [CrossRef]
  35. El-Karamany, A.S.; Ezzat, M.A. On the three-phase-lag linear micropolar thermoelasticity theory. Eur. J. Mech.-A Solids 2013, 40, 198–208. [Google Scholar] [CrossRef]
  36. Alharbi, A.M.; Said, S.M.; Abd-Elaziz, E.M.; Othman, M.I.A. Mathematical model for a magneto-thermoelastic micropolar medium with temperature-dependent material moduli under the effect of mechanical strip load. Acta Mech. 2021, 232, 2331–2346. [Google Scholar] [CrossRef]
  37. Marin, M.; Othman, M.I.A.; Abbas, I.A. An extension of the domain of influence theorem for generalized thermoelasticity of anisotropic material with voids. J. Comput. Theor. Nanosci. 2015, 12, 1594–1598. [Google Scholar] [CrossRef]
  38. Sharma, H.; Kumari, S.; Kumar, A. Study of micropolar thermo-elasticity. Adv. Math. Sci. Appl. 2020, 19, 929–941. [Google Scholar]
  39. Hilal, M.I.M.; Abd-Elaziz, E.M.; Hanoura, S.A. Reflection of plane waves in magneto-micropolar thermoelastic medium with voids and one relaxation time due to gravity and two-temperature theory. Indian J. Phys. 2021, 95, 915–924. [Google Scholar] [CrossRef]
  40. Kumar, R.; Prasad, R.; Kumar, R. Thermoelastic interactions on hyperbolic two-temperature generalized thermoelasticity in an infinite medium with a cylindrical cavity. Eur. J. Mech.-A Solids 2020, 8, 104007. [Google Scholar] [CrossRef]
  41. Lianngenga, R.; Singh, S.S. Reflection of coupled dilatational and shear waves in the generalized micropolar thermoelastic materials. J. Vib. Control 2020, 26, 1948–1955. [Google Scholar] [CrossRef]
  42. Othman, M.I.A.; Said, S.; Marin, M. A novel model of plane waves of two-temperature fiber-reinforced thermoelastic medium under the effect of gravity with three-phase-lag model. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 4788–4806. [Google Scholar] [CrossRef]
  43. Abouelregal, A.E.; Zenkour, A.M. Two-temperature thermoelastic surface waves in micropolar thermoelastic media via dual-phase-lag model. Adv. Aircr. Spacecr. Sci. 2017, 4, 711–727. [Google Scholar]
  44. Guesmia, A.; Muñoz Rivera, J.E.; Sepúlveda Cortés, M.A.; Vera Villagrán, O. Well-posedness and stability of a generalized micropolar thermoelastic body with infinite memory. Q. J. Math. 2021, 72, 1495–1515. [Google Scholar] [CrossRef]
  45. Marin, M. Harmonic vibrations in thermoelasticity of microstretch materials. J. Vib. Acoust. Trans. ASME 2010, 132, 044501. [Google Scholar] [CrossRef]
  46. Kumar, R.; Abbas, I.A. Deformation due to thermal source in micropolar thermoelastic media with thermal and conductive temperatures. J. Comput. Theor. Nanosci. 2013, 10, 2241–2247. [Google Scholar] [CrossRef]
  47. Shaw, S.; Mukhopadhyay, B. Moving heat source response in micropolar half-space with two-temperature theory. Contin. Mech. Thermodyn. 2013, 25, 523–535. [Google Scholar] [CrossRef]
  48. Ezzat, M.A.; Awad, E.S. Constitutive relations, uniqueness of solution and thermal shock application in the linear theory of micropolar generalized thermoelasticity involving two temperatures. J. Therm. Stress. 2010, 33, 226–250. [Google Scholar] [CrossRef]
  49. Quintanilla, R. Exponential stability and uniqueness in thermoelasticity with two temperatures, Dynamics Continous. Discret. Impulsive Sys. Ser. A Math. Anal. 2004, 11, 57–68. [Google Scholar]
  50. Quintanilla, R. A well posed problem for the Dual-Phase-Lag heat conduction. J. Therm. Stress. 2008, 31, 260–269. [Google Scholar] [CrossRef]
  51. Chiriţă, S. On the time differential dual-phase-lag thermoelastic model. Meccanica 2017, 52, 349–361. [Google Scholar] [CrossRef]
  52. Chiriţă, S.; Ciarletta, M.; Tibullo, V. On the thermomechanic consistency of the time differential dual-phase-lag models of heat conduction. Int. J. Heat Mass Transf. 2017, 114, 277–285. [Google Scholar] [CrossRef]
  53. Chiriţă, S.; Ciarletta, M.; Tibullo, V. The wave propagation in the time differential dual-phase-lag thermoelastic model. Proc. R. Soc. A 2015, 471, 20150400. [Google Scholar] [CrossRef]
  54. Chiriţă, S. High-order approximations of three-phase-lag heat conduction model: Some qualitative results. J. Therm. Stress. 2018, 41, 608–626. [Google Scholar] [CrossRef]
  55. Chiriţă, S.; D’Apice, C.; Zampoli, V. The time differential three-phase-lag heat conduction model: Thermodynamic compatibility and continuous dependence. Int. J. Heat Mass Transf. 2016, 102, 226–232. [Google Scholar] [CrossRef]
  56. Praveen, A.; Kumar, S.S.; Devinder, P. A two dimensional fibre reinforced micropolar thermoelastic problem for a half-space subjected to mechanical force. Theoret. Appl. Mech. 2015, 42, 11–25. [Google Scholar]
Figure 1. Schematic of the half-space under the influence of external heat source.
Figure 1. Schematic of the half-space under the influence of external heat source.
Mathematics 10 01552 g001
Figure 2. Dynamic temperature fluctuations θ for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 2. Dynamic temperature fluctuations θ for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g002
Figure 3. Conductive temperature fluctuations φ for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 3. Conductive temperature fluctuations φ for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g003
Figure 4. Normal displacement fluctuations u for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 4. Normal displacement fluctuations u for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g004
Figure 5. Tangential displacement fluctuations v for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 5. Tangential displacement fluctuations v for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g005
Figure 6. The normal stress variation σ x x for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 6. The normal stress variation σ x x for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g006
Figure 7. The normal stress variation σ y y for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 7. The normal stress variation σ y y for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g007
Figure 8. The shear stress variation σ y x for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 8. The shear stress variation σ y x for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g008
Figure 9. The variations of the micro-rotation ω for various models of thermoelasticity in the presence and absence of micropolarity.
Figure 9. The variations of the micro-rotation ω for various models of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g009
Figure 10. The variations of the tangential couple stresses m y z for different theories of thermoelasticity in the presence and absence of micropolarity.
Figure 10. The variations of the tangential couple stresses m y z for different theories of thermoelasticity in the presence and absence of micropolarity.
Mathematics 10 01552 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abouelregal, A.E.; Marin, M.; Alsharari, F. Thermoelastic Plane Waves in Materials with a Microstructure Based on Micropolar Thermoelasticity with Two Temperature and Higher Order Time Derivatives. Mathematics 2022, 10, 1552. https://doi.org/10.3390/math10091552

AMA Style

Abouelregal AE, Marin M, Alsharari F. Thermoelastic Plane Waves in Materials with a Microstructure Based on Micropolar Thermoelasticity with Two Temperature and Higher Order Time Derivatives. Mathematics. 2022; 10(9):1552. https://doi.org/10.3390/math10091552

Chicago/Turabian Style

Abouelregal, Ahmed E., Marin Marin, and Fahad Alsharari. 2022. "Thermoelastic Plane Waves in Materials with a Microstructure Based on Micropolar Thermoelasticity with Two Temperature and Higher Order Time Derivatives" Mathematics 10, no. 9: 1552. https://doi.org/10.3390/math10091552

APA Style

Abouelregal, A. E., Marin, M., & Alsharari, F. (2022). Thermoelastic Plane Waves in Materials with a Microstructure Based on Micropolar Thermoelasticity with Two Temperature and Higher Order Time Derivatives. Mathematics, 10(9), 1552. https://doi.org/10.3390/math10091552

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop