Next Article in Journal
On Reliability Estimation of Lomax Distribution under Adaptive Type-I Progressive Hybrid Censoring Scheme
Previous Article in Journal
Machine Learning Techniques Applied to Predict Tropospheric Ozone in a Semi-Arid Climate Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermo-Optical Mechanical Waves in a Rotating Solid Semiconductor Sphere Using the Improved Green–Naghdi III Model

1
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
2
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500093 Brasov, Romania
3
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(22), 2902; https://doi.org/10.3390/math9222902
Submission received: 13 October 2021 / Revised: 10 November 2021 / Accepted: 12 November 2021 / Published: 15 November 2021
(This article belongs to the Section Engineering Mathematics)

Abstract

:
The current study investigates thermophotovoltaic interactions using a new mathematical model of thermoelasticity established on a modification of the Green–Naghdi model of type III (GN-III). The basic equations, in which the heat transfer is in the form of the Moore–Gibson–Thompson (MGT) equation, are derived by adding a single delay factor to the GN-III model. The impact of temperature and electrical elastic displacement of semiconductors throughout the excited thermoelectric mechanism can be studied theoretically using this model. The proposed model was used to investigate the interactions between the processes of thermoelastic plasma in a rotating semiconductor solid sphere that was subjected to a thermal shock and crossed to an externally applied magnetic field. The influence of rotation parameters on various photothermal characteristics of silicon solid was presented and explored using the Laplace technique.

1. Introduction

Several developments have been made in the field of ultrafast carrier dynamics in semiconductor materials during the last five decades. The motives for the study include not only basic scientific interests, but also the application of semiconductor optoelectronic and electronic devices and the growing need for faster interaction and information processing. To continuously improve microelectronic semiconductor devices, various dynamic processes in semiconductors must be understood. As a result, the excitation of semiconductors to the disequilibrium state, as well as the numerous relaxation processes that follow, must be thoroughly investigated.
The photoacoustic (PA) signal can be detected on the front side of a material that has been irradiated with a modulated optical beam or on the back side of a material that has been irradiated with an unmodulated optical beam. The heat-transmission detection system may be used for a variety of applications, including the assessment of heat propagation mechanisms in solid materials [1]. When this approach is used for semiconductors, it is possible to obtain additional information about carrier transport properties [2]. Because of the periodic generation of excess carriers into semiconductors, the carrier thermalization and recombination processes will produce thermal waves. The semiconductor materials experience mechanical stress when electron-hole pairs occur.
The sound wave is formed when photo-induced free carriers induce periodic elastic strain in the material [3]. The semiconductor coefficients are also influenced by the carrier density, however, the effect is small [4]. Some materials, such as semiconductors, offer a variety of physical features that might aid research. The concept of thermoelasticity only allows semiconductor materials to be classified as elastic materials. The significance of semiconductors in technological advancement was recently highlighted when they were utilized to create electrical energy of sunlight, even when being subjected to ultrasound [5]. Superconductors are nanoparticles utilized in electrical and electronic engineering, and have a wide range of uses in contemporary industry, including transistors, displays, and solar cells.
The concept of photovoltaic energy has only recently been introduced into semiconductor media in order to generate alternative energy sources. To determine the relationship between the photothermal equations and thermoelasticity, a number of theoretical models have been examined. Gordon et al. [6] first used electronic displacements in a photo-thermal spectroscopic. Photoacoustic spectroscopy [7] is a specific and sensitive method that uses a laser source to determine the acoustic velocity of a type of semiconductor.
In photothermal processing techniques, wave propagation is used in many applications in today’s engineering industries during electro-deformations of elastic semiconductor material [8]. Abouelregal [9] investigated the responses in a solid, rotating semiconductor cylinder because of changing heat flow. Abouelregal et al. [10,11,12] studied the impact of an additional carrier on a semiconducting spinning half-space exposed to a normal force using the Green and Naghdi models.
Due to the more realistic predictions from the theory of uncoupled or coupled thermoelasticity, interest is growing these days in theories of extended thermoelasticity. The phenomenon of the unlimited speed of heat transmission inherent in the classical coupled thermoelasticity theory is being addressed by modified generalized theories [13,14,15,16] developed during the last five decades. Green and Naghdi [17,18,19] presented Green–Naghdi models (type II and III) as alternative expanded models for thermoelasticity with or without energy dissipation. The type II model does not allow for energy dissipation and is thus regarded as a special case of type III, which does allow for energy dissipation.
Numerous scholarly publications devoted to the study and interpretation of the MGT equation have emerged in recent years. The theory also focuses on third-order governing equations, which are useful in many applications of fluid mechanics [20]. Quintanilla [21,22] is developing a new MGT heat conduction model. After adding the relaxation parameter to the GN-III model, Abouelregal et al. [23,24,25] developed the suggested modified heat equation and employed the energy equation. Since the MGT equation was introduced, the number of papers on this idea has exploded [26,27,28,29]. Other results on bodies with the microstructure can be found in the papers [30,31,32,33].
Several authors have used classical models to investigate the uncoupled and coupled systems of photo-thermoelastic equations and the impacts on semiconductors. A survey of the literature found that there is no study on the volatile evaluation of semiconductor materials subjected to heat and optical loading as well as the temperature-dependent properties of materials. Among the many applications of surface waves in organic semiconductors are nuclear sciences, manufacturing engineering, electric utilities, subsea constructions, compressed gases, aircraft, chemical vessels, and metals, to name a few.
For the first time in this paper, the photothermal MGT heat conduction model (MGTPT), which describes photo-excited carriers and acoustic waves in semiconductors, are introduced. In temporal evolutions of a mainly hyperbolic nature, these equations are of the third order. Due to the presence of the thermal relaxation coefficient in front of the third-order time derivative, the MGTPT model allows for the limited rate of heat spread. To the best of our knowledge, photothermal waves along the interface of semiconducting materials have not been studied before using the MGTPT model.
The Moore–Gibson–Thompson (MGTPT) photothermal heat conduction equation is used to investigate effects throughout photothermal processes in an isotropic, homogeneous, thermoelastic semiconductor solid sphere. The analytical solution is found using numerical inversion methods, and the governing equations, including associated plasmas, thermal processes, and elastic waves, are provided in the field of Laplace transforms. Physical field variables and some analytic comparisons are shown and graphed analytically. The findings were compared to the findings of the study conducted by the researchers.

2. Formulating Mathematical Modeling

Thermoelastic and electronic deformation processes provide the following motion equations with an external force:
σ i j , j + F i = ϱ u i ¨ .
Relationships between constitutive and strain-displacement:
σ i j = C i j k l e k l ( β i j θ + d n i j N ) . e i j = 1 2 ( u i , j + u j , i ) .  
If the medium rotates at an angular velocity Ω = Ω   n uniformly, n being the axis direction of the rotation unit vector, the equation of motion (1) in the rotating reference frame includes two further terms (forces). The centripetal acceleration force ϱ Ω × ( Ω × u ) is an extra initial term, while the second term is related to time-varying motion and Corioli’s acceleration 2 ϱ ( Ω × u ˙ ) . Then, in a rotating reference frame, the equation of motion may be expressed as [25]:
σ i j , j + F i = ϱ u i ¨ + ϱ Ω × ( Ω × u ) i + 2 ϱ ( Ω × u ˙ ) i   = ϱ u i ¨ + 2 ρ ( Ω × u ˙ ) i + ϱ [ ( Ω · u ) Ω Ω 2 u ] i .
The corresponding plasma-thermal-elastic wave equation for the increase in carrier density N is as follows [34,35,36]
( D E i j N , j ) , i = ϱ N t + 1 τ N + κ θ + G ,
Vasilev and Sandomirskii [36] note that the coupled factor with thermal activation is minimal at lower temperatures when working with a harmonic modification laser.
The conventional Fourier law predicts that propagation speeds will never be infinite. A slight change in the beginning data may be demonstrated to influence the full solution throughout the entire space using this structure. Cattaneo-Vernotte developed a broader Fourier law in [37,38], by adding thermal relaxation τ 0 to the heat flow vector q as
( 1 + τ 0   t ) q = K i j θ .
The thermal relaxation parameter (single delay) τ 0 is known as the material-dependent constant, and it is an important element of this research. In terms of physics, the parameter τ 0 is the time it takes for a volume element to achieve constant heat conductivity after being exposed to a temperature gradient. This time lag may be translated into a variety of phenomena and contexts, such as when models are used to explore lithotripsy, thermal therapy, ultrasonic cleaning, and high-frequency ultrasound sonochemistry (HFU).
Based on the GN-III type, the Fourier law may be written as [18]
q = K i j   θ K i j *   ϑ ,
where ϑ represents a thermal displacement that fulfills ϑ ˙ = θ , and K i j * signify the rates of the thermal conductivity. The equation of balance energy may be expressed as [39,40]
ϱ C E θ t + T 0 t ( β i j e i j ) = · q + Q .
The energy Equation (7), when combined with the enhanced Fourier law presented in (6), has the same flaw as Fourier’s normal model, anticipating that the heat waves will spread quickly. After including the relaxation parameter into the GN-III model, Quintanilla [21,22] and Abouelregal et al. [23,24,25] developed the suggested updated heat equation. Fourier’s law was proposed to be [21,22]
( 1 + τ 0   t ) q = K i j   θ K i j *   ϑ .
Consider the case when the semiconductor elastic medium is subjected to external light beams, and the exciting free electrons produce a carrier-free charge density with semiconductor gap energy E g . A shift in electronic deformation and elastic vibrations occurs as a result of the absorbed optical energy. The overall shape of the heat conductivity equation will be affected by thermal-elastic-plasma waves in this case.
When the recombination of electron-hole pairs is taken into account, a portion of the optical energy received is thermalized. For semiconductor materials having a plasma effect, the modified Fourier law can be stated as follows:
When considering electron-hole pair recombination, the proportion of absorbed optical energy must be taken into account as
( 1 + τ 0   t ) q = K i j   θ K i j *   ϑ E g τ N d x .
The influence of heat production by carrier volume and surface de-excitation in the sample is represented by the last term in R.H.S. of Equation (9).
The result of differentiating the previous equation with respect to the position vector x , is
( 1 + τ 0   t ) ( · q ) = · ( K i j   θ ) · ( K i j *   ϑ ) E g τ N .
Incorporating Equation (10) into Equation (7) yields the proposed heat transfer equation, which describes the interplay of thermal, plasma, and elastic waves
( 1 + τ 0   t ) [ ϱ C E 2 θ t 2 + T 0 2 t 2 ( β i j u i , j ) ϱ Q t ] = ( K i j θ ˙ , j ) , i + ( K i j * θ , j ) , i + E g τ N t .
We assumed that an initial magnetic field H pervades the neighboring free space. This generates an induced electro field E and an induced magnetic field   h , both of which satisfy Maxwell’s magnetic equations and are adequate for slow-moving media.
J = × h ,   × E = μ 0 h t ,   E = μ 0 ( h t × H ) ,   h = 0 ,
τ i j = μ 0 [ H i h j + H j h i H k h k δ i j ]
where μ 0 is the magnetic permeability, J is the current density, and τ i j is the Maxwell stress tensor.

3. Application Problem

In this section, we will present a theoretical application to clarify the model proposed in the previous section. We investigate an isotropic, homogeneous, and perfect electrically conductive solid sphere of the radius R with a traction-free outer surface and time-dependent changing heat. Furthermore, we assume that the body is devoid of any sources of heat. We use a spherical coordinate system ( r , ϑ , ϕ ) with the values 0 r R , 0 ϑ 2 π , and 0 ϕ 2 π . Due to symmetry, all of the functions investigated are considered to be dependent on the distance r and the time t . The following are the components of the displacement vector in displacement-strain relationships:
u ρ = u ( r , t ) ,   u ϕ ( r , t ) = 0 = u ϑ ( r , t ) ,   e r r = u r ,   e ϕ ϕ = e ϑ ϑ = u r .  
The dilatation e may be calculated as follows:
e = e r r + e ϕ ϕ + e ϑ ϑ = u r + 2 u r = 1 r 2 ( r 2 u ) r
The photothermal stresses (2) will be in the form
σ r r = ( λ + 2 μ ) u r + 2 λ u r ( 3 λ + 2 μ ) ( α t θ + δ n N ) , σ ϑ ϑ = σ ϕ ϕ = λ u r + 2 ( μ + λ ) u r ( 3 λ + 2 μ ) ( α t θ + δ n N ) ,
where α t is the linear thermal expansion coefficient, δ n is the electronic deformation coefficient, and λ , μ are the Lame’s constants.
If the rotation is about the axial axis of the cylinder, i.e., Ω = ( 0 , 0 , Ω ) , the equation of motion can be expressed in the following terms: the impact of the magnetic field F r and the force due to the rotation of the body ϱ Ω 2 u . In this case, Equation (11) can be written as:
σ r r r + 2 r ( σ r r σ ϑ ϑ ) + F r = ϱ 2 u t 2 ϱ Ω 2 u
Assume the surface of the sphere is immersed in a magnetic field that acts in the ϕ -direction of constant strength H 0 = ( 0 , 0 , H 0 ) . According to Equation (12), we obtain
E = ( 0 , μ 0 H 0 u t , 0 ) ,   J = ( 0 , r ( 1 r 2 ( r 2 u ) r ) , 0 ) , h = ( 0 , 0 , 1 r 2 ( r 2 u ) r ) .
The Lorentz force F r can be expressed as
F r = μ 0 ( J × H 0 ) r
Thus, we have F r and Maxwell’s stress τ r r from Equations (18) and (19) as
F r = μ 0 H 0 2 r ( 1 r 2 ( r 2 u ) r ) , τ r r = μ 0 H 0 2 r 2 ( r 2 u ) r
After plugging Equations (16) and (20) into Equation (17), we get
( λ + 2 μ + μ 0 H 0 2 ) r ( 1 r 2 ( r 2 u ) r ) γ θ r d n N r = ϱ 2 u t 2 ϱ Ω 2 u .
where { γ ,   d n } = ( 3 λ + 2 μ ) { α t , δ n } .
Equation (21) can be rewritten as
( λ + 2 μ + μ 0 H 0 2 ) 2 e γ 2 θ d n 2 N = ϱ 2 e t 2 ρ Ω 2 e .
The modified photothermal MGTE heat Equation (11) is as follows in the absence of any heat sources ( Q = 0 ):
( 1 + τ 0   t ) [ ϱ C E 2 θ t 2 + γ T 0 2 e t 2 ] = K 2 θ ˙ + K * 2 θ + E g τ N t .
where 2 = 2 r 2 + 2 r r = 1 r 2 r ( r 2 r ) .
The non-dimensional variables of the partial differential equation are easy to construct. As a result, the following non-dimensional quantities are described:
{ r ,   u } = v 0 η { r , u } ,   { t , τ 0 , τ } = v 0 2 η { t , τ 0 , τ } , { θ , N } = 1 ρ v 0 2 { γ θ , d n n } , η = ρ C E K , Ω = Ω c 0 2 η ,   { σ i j , τ r r } = 1 ρ v 0 2 { σ i j , τ r r } ,   v 0 2 = v 1 2 + v a 2 ,   v 1 = λ + 2 μ ϱ , v a = μ 0 H 0 2 ϱ  
If the primes are omitted, the system of equations may be expressed by the following formula:
( 1 + τ 0   t ) [ 2 θ t 2 + ε 1 2 e t 2 ] = ( t + ω * ) 2 θ + ε 2 N t ,
2 e 2 θ 2 Ν = 2 e t 2 Ω 2 e ,
2 Ν = g 1 Ν t + g 2 Ν + g 3 θ ,
σ r r = β 2 u r + ( 1 β 2 ) e θ Ν , σ ϑ ϑ = β 2 u r + ( 1 β 2 ) e θ Ν ,
where
β 2 = 2 μ λ + 2 μ ,   ε 1 = γ 2 T 0 ϱ 2 C E c 0 2 ,   ω * = K * v 0 2 η K , ε 2 = γ E g v 0 2 τ d n K ,   g 1 = ϱ D E η ,   g 2 = 1 D E η τ 1 ,   g 3 = κ d n γ η 2 D E c 0 2 .
The initial conditions can be assumed to be
u ( ρ , 0 ) = 0 = u ( r , 0 ) r ,   Ν ( r , 0 ) = 0 = Ν ( r , 0 ) r ,   θ ( r , 0 ) = 0 = θ ( r , 0 ) r .
We also assume that the following boundary conditions are hold:
θ ( a , t ) = θ 0 H ( t ) ,   t > 0
σ r r ( R , t ) = 0
D E N ρ = s v N   at   ρ = ρ 0 ,
where θ 0 is a constant, s v is the surface recombination speed, and H ( t ) is the Heaviside function.

4. Solution Using the Laplace Transform

When solving differential equations with constant coefficients, the Laplace transform can be applied. For example, in the design of control mechanisms, the Laplace transform method is essential. The properties of both the Laplace transform and the inverse Laplace transform are also used to evaluate the dynamic control scheme. The following equation defines the Laplace transform of a function g ( t ) , often known as   L [ g ( t )   ] or by g ¯ ( s )
L [ g ( t )   ] = g ¯ ( s ) = 0 g ( t ) e s t d t ,   s > 0 .
We can derive the following findings by using the Laplace transform method to the governing Equations (25)–(28):
( 2 ψ ) θ ¯ = ψ ε 1 e ¯ ε 2 s N ¯ ,
( 2 s 2 + Ω 2 ) e ¯ = 2 θ ¯ + 2 Ν ¯ ,
( 2 g 4 ) Ν ¯ = g 3 θ ¯ ,
σ ¯ r r = β 2 d u ¯ d r + ( 1 β 2 ) e ¯ θ ¯ Ν ¯ ,
σ ¯ ϑ ϑ = β 2 u ¯ r + ( 1 β 2 ) e ¯ θ ¯ Ν ¯ ,
where ψ = s 2 ( 1 + τ 0   s ) / ( s + ω * ) .
When we decouple Equations (35) and (37), we obtain
( 6 α 2 4 + α 1 2 α 0 ) { θ ¯ , N ¯ , e ¯ } = 0 ,
where α 2 , α 1 and α 0 are defined as follows:
α 2 = s 2 Ω 2 + g 7 + g 6 g 3 , α 1 = g 7 ( s 2 Ω 2 ) + g 8 + g 6 g 5 g 3 , α 0 = g 8 ( s 2 Ω 2 ) , g 4 = s g 1 + g 2 ,   g 5 = g 4 g 3 ,   g 6 = g 3 ψ ε 1 , g 7 = g 4 + ψ ,   g 8 = g 4 ψ + s g 3 ε 2 .
If the parameters   λ 1 2 , λ 2 2 , and λ 3 2 are the roots of the equation
λ 6 α 2 λ 3 + α 1 λ 2 α 0 = 0 ,
then Equation (40) can be expressed as
( 2 λ 1 2 ) ( 2 λ 2 2 ) ( 2 λ 3 2 ) { e ¯ , θ ¯ , N ¯ } = 0 ,
The roots of Equation (42) can be calculated as
λ 1 2 = 1 3 [ 2 β 0 sin ( γ 0 ) + α 2 ] , λ 2 2 = 1 3 β 0 [ sin ( γ 0 ) + 3 cos ( γ 0 ) ] + 1 3 α 2 , λ 2 2 = 1 3 β 0 [ 3 cos ( γ 0 ) sin ( γ 0 ) ] + 1 3 α 2 ,
with
β 0 = α 2 2 3 α 1 ,   γ 0 = 1 3 sin 1 ( 2 α 2 3 9 α 2 α 1 + 27 α 0 2 β 0 3 ) .
Equation (43) has a general bounded solution that is can be expressed as
{ e ¯ ,   θ ¯ , Ν ¯ } = 1 r i = 1 3 { 1 , L i , H i } A i I 1 / 2 ( λ i r ) .
In Equation (46), the second type of modified Bessel function of order n is indicated by I n (   .   ) . Three integral parameters are dependent on s : A i ,   ( i = 1 ,   2 ,   3 ) are assumed. Moreover, A i is related to two different variables, L i and M i , by plugging in Equation (46) into Equations (35)–(40) to get
H i = g 3 ( λ i 2 s 2 ) λ i 4 g 5 λ i 2 ,   L n = ( λ i 2 s 2 ) ( λ i 2 g 4 ) λ i 4 g 5 λ i 2 ,   i = 1 ,   2 ,   3 .
The following well-known Bessel function can be used to get the displacement u ¯ :
x 3 2 I 1 2 ( x ) d x = x 3 2   I 3 2 ( x )
Then, we have
u ¯ = 1 r i = 1 3 1 λ i A i I 3 / 2 ( λ i r ) .
According to the following relationships, for each positive x , modified Bessel I n meets the following relatives:
I 1 2 ( x ) = 2 π x sin h ( x ) ,   I 3 2 ( x ) = 2 π x ( cos h ( x ) sin h ( x ) x ) .
If we insert relation (50) to Equations (46) and (49), then we get
{ e ¯ ,   θ ¯ , Ν ¯ } = 2 π r i = 1 3 { 1 , L i , H i } A i λ i r sin h ( λ i r )
u ¯ = 2 π i = 1 2 A i r λ i 3 2 ( cos h ( λ i r ) sin h ( λ i r ) λ i r )
This leads to a closed form of the final solutions for thermal stresses:
σ ¯ r r = 2 π i = 1 3 A i r ( λ i r ) 5 / 2 ( 2 β 2 λ i r cos h ( λ i r ) + ( 2 β 2 + Ω i λ i 2 r 2 ) sin h ( λ i r ) ) ,
σ ¯ ϑ ϑ = 2 π i = 1 3 A i r ( λ i r ) 5 / 2 ( ( 2 + 3 β 2 ) λ i r cos h ( λ i r ) + ( 2 + Ω i λ i 2 r 2 + β 2 ( 3 + λ i 2 r 2 ) ) sin h ( λ i r ) )
where Ω i = 1 L i H i .
When we consider the non-dimensional Maxwell’s stress τ r r , we get the following:
τ ¯ r r = v a 2 v 0 2 2 π r i = 1 3 A i λ i r sin h ( λ i r )
After applying the Laplace transform, the boundary conditions (31)–(33) take on the following forms.
θ ¯ = θ 0 s ,   σ ¯ r r ( R , s ) = 0 ,   D E N ¯ r | r = R = s f N ¯ ( R , s ) .
Equation (56) is replaced with Equations (46) and (53) to produce
2 π R i = 1 3 L i A i λ i R sin h ( λ i R ) = θ 0 s ,
i = 1 3 A i R ( λ i R ) 5 / 2 ( ( 2 β 2 + Ω i λ i 2 R 2 ) sin h ( λ i R ) 2 β 2 λ i R cos h ( λ i R ) ) = 0 ,
2 π i = 1 3 A i R ( λ i R ) 5 / 2 ( ( 2 + λ i 2 R 2 ) sin h ( λ i R ) 2 λ i R cos h ( λ i R ) ) = s f D E 2 π R i = 1 3 H i A i λ i R sin h ( λ i R ) .
To obtain the parameter values A i ,   ( i = 1 ,   2 ,   3 ) , we solve the system (57)–(59). This research used a numerical technique based on Fourier series expansion to invert Laplace transforms [41,42]. The following approach may be used to invert any function in the Laplace domain to the time domain:
Γ ( ρ , t ) = e c t t ( 1 2 Γ ¯ ( ρ , c ) + R e j = 1 N f Γ ¯ ( ρ , c + i j π t ) ( 1 ) j )
where N f stands for the number of terms, R e stands for the real part, and i is for the imaginary number unit. Many previous studies and all numerical experiments have shown that the parameter c in the previous relation can achieve the formula ≅ 4.7/t, allowing for faster convergence of numerical results [43].

5. Numerical Results and Discussion

We present a numerical example in which we need to derive mathematical solutions for several physical quantities whose behavior is largely influenced by the assumptions stated in the topic. Using Mathematica’s software, the influence of the Modified Moore–Gibson–Thompson (MGTPT) heat equation on the physical fields under consideration is now shown in the form of graphical presentations and tables. For the theoretical study, isotropic silicon (Si) is utilized as the solid semiconductor material. The physical parameters used are as follows [12]:
λ = 3.64 × 10 10   kg   m 1 s 2 ,   μ = 5.46 × 10 10   kg   m 1 s 2 ,   ρ = 2330   kg   m 3 , K = 1.51   W   m 1 K 1 ,   C E = 6.95 × 10 2   J   kg   K 1 ,   d n = 9 × 10 31   m 3 ,   T 0 = 300   K   E g = 1.11   eV ,   D E = 2.5 × 10 3   m 2   s 1 ,   s f = 2   m   s 1 ,   τ = 5 × 10 5   s .
The numerical method given in (60) was used to compute the numerical distribution of non-dimensional field variables along the radial direction of the cylinder. The numerical results for t = 0.12 and R = 1 are obtained. It is possible to study the behavior of all field variables in three different scenarios.

5.1. Comparison of Several Thermoelasticity Models

The Moore–Gibson–Thompson photothermal (MGTPT) model suggested in the presented work is a generality of numerous previous photothermal elasticity models (CPTE, PLS, PGN-II, PGN-III, and MGTPT). The intent of presenting them as explained earlier is not only to generalize but also to resolve some of the physical inconsistencies found in some of the earlier models. In Section 1 and Section 2, these discrepancies are discussed as well as how to address them. In this section, a new photothermal model is compared to earlier models in order to verify them. Before we make the comparisons, we will first explain how to derive the previous models (CPTE, PLS, PGN-II, and PGN-III) from the modified MGTPT model given in Equation (11) in some cases.
It is possible to get the coupling photo-thermoelasticity theory (CPTE) by ignoring the thermal parameters ( K i j * = τ 0 = 0 ). When the Green and Naghdi parameter K i j * = 0 , the generalized and photothermal-elasticity model with a relaxation time (PLS) is produced. The generalized photothermal model (PGN-III) without a relaxation parameter based on Green and Naghdi theory II may be developed when τ 0 = 0 . It is possible to use the photothermal Green and Naghdi model (PGN-II) if we remove the term that contains the parameter K i j in the heat Equation (11). When the parameters K i j * , τ 0 > 0 , the modified Moore–Gibson–Thompson photothermal model (MGTPT) is applied.
To clarify the study and compare the different models in more depth, as well as for the practical purposes of researchers in this field, the numerical results will be presented in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 and Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. In this case, we take τ = 0.01 , t = 0.15 , and Ω = 5 .
First, we will review several key findings that apply to all distributions.
  • As seen in the graphs and tables, the thermal parameters τ 0 and K* have a substantial influence on all studied field variables.
  • The phenomenon of restricted heat transfer in the case of the enhanced generalized thermoelastic theory (MGTPT) is clearly shown in all figures and tables. In other words, thermal and mechanical waves propagate at finite speeds within the medium, contrary to the classical theory.
  • For all physical domains, all curves for the different models converge within the mean, where r tends to be zero to satisfy the regularity requirement.
  • The introduction of the relaxation time into the thermal conductivity Equation (11) in the Green and Naghdi model of the third type (PGN-III) reduced the propagation speed of optical and mechanical thermal waves within the medium. This important conclusion demonstrates the importance of the proposed model, in which the problem of the physical inconsistency found in the CPTE and PGN-III models has been solved. This is fully consistent with the claims of Quintanilla [21] and Abouelregal et al. [23,24,25,26].
  • There is convergence of the results in the case of the LS and MGTE models, due to the inclusion of the relaxation coefficient, which led to a slower decrease in temperature.
  • On the surface of the semiconducting solid sphere, where the boundary conditions manifest, the coupled photothermal model (CPTE) and modified and generalized photothermal models (PLS, PGNII, PGNIII, and MGTPT) provide values that differ in magnitude but behave similarly.
Now, we will study some observations for each of the studied field variables.
Table 1 and Figure 1 illustrate the temperature variance θ vs. r for all photo-thermo-elasticity models across a wide range of 0 < r 1 . As can be seen in this graph, the magnitude of the temperature is larger for the PGN-III and CPTE models than for the other models (PLS, GN-II, and MGTPT). It is worth noting that the rate of temperature propagation is limited and corresponds to the physical behavior of photo-thermoelastic materials. The temperature takes its largest value at the surface of the sphere, where the thermal shock effect is present, and it gradually decreases as we head inside the body, in contrast to the increasing radius.
The different models for the distribution of displacement u across radius r are compared in Figure 2 and Table 2. The numerical data and graphs show that the displacement u starts with negative values, then rapidly increases until it reaches its maximum positive value, then gradually fades to zero. It can also be seen that the displacement curves converge in the case of all the different models, despite the different values. Thus, it can be concluded that the effect of thermal parameters is weak on the displacement distributions.
The five curves predicted by various photo-thermoelasticity models with and without thermal parameters for the carrier density N against the radial distance r are shown in Figure 3 and Table 3. Due to the recombination of charge carriers, the carrier density peaks at the surface of the sphere r = 1 and progressively falls as the radial distance increases until it approaches the steady state. The results in the case of the GN-III model are characterized by a larger carrier density N than the MGTPT, with a similar behavior to the PLS and PGN-II models.
The fluctuation of the radial stress σ r r with respect to the distance for different models is illustrated in Figure 4 and Table 4. It should be noted that the pressure always starts from zero to respect the boundary requirements, and then gradually decreases until it reaches the lower value at a certain position near the surface of the sphere, and then gradually increases to zero. Figure 5 and Table 5 illustrate the fluctuation of hoop stress σ ϑ ϑ versus r for several models of photo-thermoelasticity. In all cases, the hoop stress is initially negative and progressively reduces until it reaches zero. The stresses emphasize the location of the medium next to the spherical surface with time, which is consistent with the data in [44]. In addition, when the radial width increases, the value and amplitude of fields measured at the surface increase. Because of this, ref. [45] provides reasons for this phenomenon.
Figure 6 and Table 6 demonstrate the Maxwell’s stress τ r r fluctuation vs. r over a wide range of 0 < r 1 for different models of photo-thermoelasticity. The speed of propagation of electromagnetic waves is limited and corresponds to the physical behavior of photo-thermal materials. We also note that the behavior of the curves is the same and that there are small differences between the different models.

5.2. The Influence of Angular Velocity of Rotation

The analysis of the propagation of plane photo-thermoelastic waves in a spinning solid appears to have received little attention. Because most large bodies have an angular motion, such as the earth, moon, and other planets, studying the propagation of planar thermoelastic or magneto-photo-thermoelastic waves in a rotating medium with thermal relaxation appears to be more realistic.
Table 7 explains the variations of different photothermal fields with different rotation parameter Ω under different photo-thermoelastic models (CPTE, PLS, PGN-II, PGN-III, and MGTPT). The numerical results were calculated at a single value of the distance r , which is r = 0.9 . In the case of all the different models, it was found that the differences in the magnitudes in all the studied fields decreased with the increase in the angular velocity, except for the displacement. We note that the displacement distribution decreases with the angular velocity in the case of the PLS, PGN-II, and MGTPT models, while it increases in the case of the CPTE and PGN-III models.
In this section, we explore the response of distinct photothermal fields as the rotation parameter Ω changes only in the case of the generalized elastic photothermal model (MGTPT). For comparative purposes, we will choose three numbers for the amount of rotational speed Ω . The rotating example is chosen with Ω = 5 ,   10 , while the nonrotating case is chosen with Ω = 0 . Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 illustrate the distribution of the data for the various fields studied.
As seen in Figure 7, Figure 8 and Figure 12, the parameter Ω has a relatively little influence on temperature θ , displacement u , and Maxwell’s stress τ r r . This phenomenon is consistent with the results reported in the references [44,45]. We can also see in Figure 7, Figure 8 and Figure 12 that as the rotation parameter Ω increases, the field variables θ , u , and τ r r decrease. Furthermore, as shown in Figure 9, Figure 10 and Figure 11, the rotational parameter Ω has a great influence on the carrier density N and the thermal stresses σ r r and σ ϑ ϑ .
Figure 10 illustrates that increasing the amount of N increases the carrier density numerical values. Figure 10 shows the change in thermal stress σ r r as a function of radius r and parameter Ω . We can see that the rotational parameter Ω has a significant impact on the radial stress profile σ r r . As the parameter Ω is increased, the stress σ r r increases. Furthermore, the amplitude of the stress σ ϑ ϑ distribution increases as Ω increases.

6. Conclusions

A novel photo-thermoelastic model, including the Moore–Gibson–Thompson equation, is presented in this paper. The type III Green–Naghdi model, as well as the Lord and Shulman equations, are all included in the improved photothermal model. The interplay of heat, plasma, and elastic waves in semiconductor materials is explored in this expanded model, which is significant. Furthermore, unlike traditional models, this model permits heat-elastic optical waves to propagate at restricted speeds. So far, only a few articles on our topic have been written in the scientific literature. As special instances of the suggested model, several thermoelastic and photothermal models may be constructed.
From the quantitative data, we may conclude that the thermal relaxation time reduces the mechanical and thermo-optical waves along the axes. It is also noted that the behavior of waves in different physical fields depends on the speed of rotation of the medium. The new photo-thermoelasticity model allows heat to go through a medium at a finite velocity rather than an infinity-speed wave as in the case of the CPTE model. Compared to the PGN-III model, the extended MGTPT model more accurately describes the photo-thermoelastic process. In contrast to previous generalized thermoelasticity models, the PGN-III model’s results show convergence with the conventional elasticity model (CPTE), which does not disappear quickly within the body.
Scientists working in fields, such as physics, material design, thermal efficiency, and geophysics will benefit greatly from the information presented in this study. The approach employed in this investigation can be used to address a variety of thermodynamic challenges.

Author Contributions

Formal analysis, M.M.; Investigation, S.A.; Project administration, A.E.A. The results were discussed by all authors, and the final paper was reviewed and approved by all of them. All authors have read and agreed to the published version of the manuscript.

Funding

Research Supporting Project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

With regard to the research, writing, and publishing of this paper, the author(s) reported no possible conflicts of interest.

Abbreviations

In the governing equations, σ i j is the stress tensor, ϱ is the density of the material, u i denotes the components of the displacement vector, F i indicates the body forces components and i , j , k = 1 ,   2 ,   3 , e i j being the strain tensor, e k k = e   d n i j = d n i δ i j are the difference in deformation potential of the conduction and valence bands, and C i j k l are the elastic constants, β i j = β i δ i j are the stress-temperature coefficients. Furthermore, θ = T T 0 represents the temperature change, T 0 is the location temperature, N is the carrier density, D E i j denotes the coefficients of the diffusion, κ is the thermal activation coupling parameter, τ is the lifetime of photogenerated electron-hole pairs, and G symbolizes the carrier source, K i j denotes the tensor of thermal conductivity, C E is the specific heat and Q is the heat supply.

References

  1. Adams, M.J.; Kirkbright, G.F. Thermal diffusivity and thickness measurements for solid samples utilising the optoacoustic effect. Analyst 1977, 102, 678–682. [Google Scholar] [CrossRef]
  2. Vargas, H.; Miranda, L.C.M. Photoacoustic and Re1ated Phototherma1 Technique. Phys. Rep. 1988, 161, 43–101. [Google Scholar] [CrossRef]
  3. Ferreira, S.O.; Ying An, C.; Bandeira, I.N.; Miranda, L.C.M.; Vargas, H. Photoacoustic measurement of the thermal diffusivity ofPb1−xSnxTe alloys. Phys. Rev. B 1989, 39, 7967–7970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Othman, M.I.A.; Eraki, E.E.M. Effect of gravity on generalized thermoelastic diffusion due to laser pulse using dual-phase-lag model. Multidiscip. Modeling Mater. Struct. 2018, 14, 457–481. [Google Scholar] [CrossRef]
  5. Stearns, R.G.; Kino, G.S. Effect of electronic strain on photoacoustic generation in silicon. Appl. Phys. Lett. 1985, 47, 1048–1050. [Google Scholar] [CrossRef]
  6. Gordon, J.P.; Leite, R.C.C.; Moore, R.S.; Porto, S.P.S.; Whinnery, J.R. Long-Transient Effects in Lasers with Inserted Liquid Samples. J. Appl. Phys. 1965, 36, 3–8. [Google Scholar] [CrossRef]
  7. Todorovic, D.M.; Nikolic, P.M.; Bojicic, A.I. Photoacoustic frequency transmission technique: Electronic deformation mechanism in semiconductors. J. Appl. Phys. 1999, 85, 7716. [Google Scholar] [CrossRef]
  8. Song, Y.; Todorovic, D.M.; Cretin, B.; Vairac, P. Study on the generalized thermoelastic vibration of the optically excited semiconducting microcantilevers. Int. J. Solids Struct. 2010, 47, 1871–1875. [Google Scholar] [CrossRef] [Green Version]
  9. Abouelregal, A.E. Magnetophotothermal interaction in a rotating solid cylinder of semiconductor silicone material with time dependent heat flow. Appl. Math. Mech. Engl. Ed. 2021, 42, 39–52. [Google Scholar] [CrossRef]
  10. Abouelregal, A.E.; Sedighi, H.M.; Shirazi, A.H. The Effect of Excess Carrier on a Semiconducting Semi-Infinite Medium Subject to a Normal Force by Means of Green and Naghdi Approach. Silicon 2021. Available online: https://link.springer.com/article/10.1007/s12633-021-01289-9#citeas (accessed on 13 October 2021).
  11. Abouelregal, A.E.; Ahmad, H.; Elagan, S.K.; Alshehri, N.A. Modified Moore–Gibson–Thompson photo-thermoelastic model for a rotating semiconductor half-space subjected to a magnetic field. Int. J. Mod. Phys. C 2021. [Google Scholar] [CrossRef]
  12. Hobiny, A.D.; Alzahrani, F.S.; Abbas, I.A. A study on photo-thermo-elastic wave in a semi-conductor material caused by ramp-type heating. Alex. Eng. J. 2021, 60, 2033–2040. [Google Scholar] [CrossRef]
  13. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  14. Green, A.E.; Lindsay, K.A. Thermoelasticity. J. Elast. 1972, 2, 1–7. [Google Scholar] [CrossRef]
  15. Tzou, D.Y. Experimental support for the lagging behaviour in heat propagation. J. Thermophys. Heat Transf. 1995, 9, 686–693. [Google Scholar] [CrossRef]
  16. Tzou, D.Y. A unified approach for heat conduction from macro to microscale. J. Heat Transf. 1995, 117, 8–16. [Google Scholar] [CrossRef]
  17. Green, A.E.; Naghdi, P.M. A Re-Examination of the Basic Postulates of Thermomechanics. Proc. R. Soc. A Math. Phys. Eng. Sci. 1991, 432, 171–194. [Google Scholar]
  18. Green, A.E.; Naghdi, P.M. On undamped heat waves in an elastic solid. J. Ther. Stresses 1992, 15, 253–264. [Google Scholar] [CrossRef]
  19. Green, A.E.; Naghdi, P.M. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  20. Lasiecka, I.; Wang, X. Moore–Gibson–Thompson equation with memory, part II: General decay of energy. J. Diff. Equ. 2015, 259, 7610–7635. [Google Scholar] [CrossRef]
  21. Quintanilla, R. Moore-Gibson-Thompson thermoelasticity. Math. Mech. Solids 2019, 24, 4020–4031. [Google Scholar] [CrossRef]
  22. Quintanilla, R. Moore-Gibson-Thompson thermoelasticity with two temperatures. Appl. Eng. Sci. 2020, 1, 100006. [Google Scholar] [CrossRef]
  23. Jangid, K.; Gupta, M.; Mukhopadhyay, S. On propagation of harmonic plane waves under the Moore–Gibson–Thompson thermoelasticity theory. Waves Random Complex Media 2021. [Google Scholar] [CrossRef]
  24. Abouelregal, A.E.; Ahmad, H.; Nofal, T.A.; Abu-Zinadah, H. Moore–Gibson–Thompson thermoelasticity model with temperature-dependent properties for thermo-viscoelastic orthotropic solid cylinder of infinite length under a temperature pulse. Phys. Scr. 2021, 96, 105201. [Google Scholar] [CrossRef]
  25. Aboueregal, A.E.; Sedighi, H.M. The effect of variable properties and rotation in a visco-thermoelastic orthotropic annular cylinder under the Moore Gibson Thompson heat conduction model. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2021, 235, 1004–1020. [Google Scholar] [CrossRef]
  26. Florea, O.A.; Bobe, A. Moore–Gibson–Thompson thermoelasticity in the context of double porous materials. Contin. Mech. Thermodyn. 2021, 33, 2243–2252. [Google Scholar] [CrossRef]
  27. Abouelregal, A.E.; Ersoy, H.; Civalek, Ö. Solution of Moore–Gibson–Thompson equation of an unbounded medium with a cylindrical hole. Mathematics 2021, 9, 1536. [Google Scholar] [CrossRef]
  28. Kaltenbacher, B.; Lasiecka, I.; Marchand, R. Wellposedness and exponential decay rates for the Moore–Gibson–Thompson equation arising in high intensity ultrasound. Control Cybern. 2011, 40, 971–988. [Google Scholar]
  29. Bazarra, N.; Fernández, J.R.; Quintanilla, R. Analysis of a Moore-Gibson-Thompson thermoelastic problem. J. Comput. Appl. Math. 2021, 382, 113058. [Google Scholar] [CrossRef]
  30. Marin, M. A domain of influence theorem for microstretch elastic materials. Nonlinear Anal. Real World Appl. 2010, 11, 3446–3452. [Google Scholar] [CrossRef]
  31. Marin, M.; Agarwal, R.P.; Mahmoud, S.R. Modeling a microstretch thermo-elastic body with two temperatures. Abstr. Appl. Anal. 2013, 2013, 583464. [Google Scholar] [CrossRef] [Green Version]
  32. Abbas, I.; Marin, M. Analytical Solutions of a Two-Dimensional Generalized Thermoelastic Diffusions Problem Due to Laser Pulse. Iran. J. Sci. Technol. Trans. Mech. Eng. 2018, 42, 57–71. [Google Scholar] [CrossRef]
  33. Marin, M.; Othman, M.I.A.; Seadawy, A.R.; Carstea, C. A domain of influence in the Moore–Gibson–Thompson theory of dipolar bodies. J. Taibah Univ. Sci. 2020, 14, 653–660. [Google Scholar] [CrossRef]
  34. Song, Y.Q.; Bai, J.T.; Ren, Z.Y. Study on the reflection of photothermal waves in a semiconducting medium under generalized thermoelastic theory. Acta Mech. 2012, 223, 1545–1557. [Google Scholar] [CrossRef]
  35. Todorovic, D.M. Plasma, thermal, and elastic waves in semiconductors. Rev. Sci. Instrum. 2003, 74, 582. [Google Scholar] [CrossRef]
  36. Vasilev, A.N.; Sandomirskii, V.B. Photoacoustic effects in finite semiconductors. Sov. Phys. Semicond. 1984, 18, 1095. [Google Scholar]
  37. Cattaneo, C. A form of heat-conduction equations which eliminates the paradox of instantaneous propagation. Comptes Rendus 1958, 247, 431–433. [Google Scholar]
  38. Vernotte, P. Some possible complications in the phenomena of thermal conduction. Compte Rendus 1961, 252, 2190–2191. [Google Scholar]
  39. 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]
  40. Kumar, H.; Mukhopadhyay, S. Thermoelastic damping analysis in microbeam resonators based on Moore–Gibson–Thompson generalized thermoelasticity theory. Acta Mech. 2020, 231, 3003–3015. [Google Scholar] [CrossRef]
  41. Roychoudhuri, S.K. On a thermoelastic three-phase-lag model. J. Therm. Stress. 2007, 30, 231–238. [Google Scholar] [CrossRef]
  42. Honig, G.; Hirdes, U. A method for the numerical inversion of Laplace transform. J. Comp. Appl. Math. 1984, 10, 113–132. [Google Scholar] [CrossRef] [Green Version]
  43. Tzou, D.Y. Macro-To Micro-Scale Heat Transfer: The Lagging Behavior; Taylor & Francis: Abingdon, UK, 1997. [Google Scholar]
  44. Singh, B.; Mukhopadhyay, S. Galerkin-type solution for the Moore–Gibson–Thompson thermoelasticity theory. Acta Mech. 2021, 232, 1273–1283. [Google Scholar] [CrossRef]
  45. Trajkovski, D.; Čukić, R. A coupled problem of thermoelastic vibrations of a circular plate with exact boundary conditions. Mech. Res. Commun. 1999, 26, 217–224. [Google Scholar] [CrossRef]
Figure 1. The temperature θ under different photothermal models.
Figure 1. The temperature θ under different photothermal models.
Mathematics 09 02902 g001
Figure 2. The displacement variation u under different photothermal models.
Figure 2. The displacement variation u under different photothermal models.
Mathematics 09 02902 g002
Figure 3. The carrier density N under different photothermal models.
Figure 3. The carrier density N under different photothermal models.
Mathematics 09 02902 g003
Figure 4. The radial stress variation σrr under different photothermal models.
Figure 4. The radial stress variation σrr under different photothermal models.
Mathematics 09 02902 g004
Figure 5. The hoop stress variation σϑϑ under different photothermal models.
Figure 5. The hoop stress variation σϑϑ under different photothermal models.
Mathematics 09 02902 g005
Figure 6. The Maxwell’s stress τrr for different models of photo-thermoelasticity.
Figure 6. The Maxwell’s stress τrr for different models of photo-thermoelasticity.
Mathematics 09 02902 g006
Figure 7. The temperature variation θ for different rotation parameter Ω.
Figure 7. The temperature variation θ for different rotation parameter Ω.
Mathematics 09 02902 g007
Figure 8. The displacement variation u for different rotation parameter Ω.
Figure 8. The displacement variation u for different rotation parameter Ω.
Mathematics 09 02902 g008
Figure 9. The carrier density N variation for different rotation parameter Ω.
Figure 9. The carrier density N variation for different rotation parameter Ω.
Mathematics 09 02902 g009
Figure 10. The radial stress variation σrr for different rotation parameter Ω.
Figure 10. The radial stress variation σrr for different rotation parameter Ω.
Mathematics 09 02902 g010
Figure 11. The hoop stress variation σϑϑ for different rotation parameter Ω.
Figure 11. The hoop stress variation σϑϑ for different rotation parameter Ω.
Mathematics 09 02902 g011
Figure 12. The Maxwell’s stress τrr for different rotation parameter Ω.
Figure 12. The Maxwell’s stress τrr for different rotation parameter Ω.
Mathematics 09 02902 g012
Table 1. The temperature variation θ under photothermal models.
Table 1. The temperature variation θ under photothermal models.
r CPTEPLSPGN-IIPGN-IIIMGTPT
00.0008855590.009658490.0123413−0.002379010.00799097
0.10.0008855590.009658490.0123413−0.002379010.00799097
0.20.009801540.009725340.01284940.01006030.00802299
0.30.02615540.009869040.0139960.03228610.00808409
0.40.0523780.01019840.01648530.06661030.00820664
0.50.09227290.01110820.02204360.1164050.00853831
0.60.1514650.0140610.03492040.1862150.00977179
0.70.2380630.02458420.06569620.2818850.0152332
0.80.3637370.06398180.1411140.4107920.0412536
0.90.5469080.2167080.3309430.5837260.1702640
10.8368970.8368970.8368970.8368970.8368970
Table 2. The displacement variation u under photothermal models.
Table 2. The displacement variation u under photothermal models.
r CPTEPLSPGN-IIPGN-IIIMGTPT
00.000234702.58 × 10−71.01 × 10−40.0005106−8.58 × 10−8
0.10.00033472.88 × 10−71.02 × 10−40.00053064−9.58 × 10−8
0.20.0007351.96 × 10−62.44 × 10−40.001130111.70 × 10−7
0.30.001275278.89 × 10−64.82 × 10−40.001870182.09 × 10−6
0.40.00204863.47 × 10−60.0009060.002828111.14 × 10−5
0.50.003178010.00012910.00167590.004087695.43605 × 10−5
0.60.004826520.0004730.00307580.00573410.000249641
0.70.00715880.001692660.00556860.007793620.00110423
0.80.009705470.005472010.00932120.009572520.00434249
0.90.003770530.009961880.006679050.002109280.00978329
1−0.1162520−0.0897956−0.1092720−0.1194910−0.08472120
Table 3. The carrier density variation N under photothermal models.
Table 3. The carrier density variation N under photothermal models.
r
CPTEPLSPGN-IIPGN-IIIMGTPT
0−0.00256381−0.000144986−0.000431833−0.00403677−5.1282 × 10−5
0.1−0.00256381−0.000144986−0.000431833−0.00403677−5.12826 × 10−5
0.2−0.00297915−0.000243067−0.000625476−0.00446883−0.000100656
0.3−0.0037235−0.00048108−0.00104769−0.00520596−0.000233686
0.4−0.00488043−0.00103429−0.00190478−0.00627114−0.000584929
0.5−0.00657716−0.00233155−0.00363361−0.00768978−0.00153045
0.6−0.00898314−0.00541422−0.00714027−0.00946911−0.00412116
0.7−0.0121715−0.0127064−0.0141794−0.0114332−0.0112066
0.8−0.0143739−0.0284206−0.0266344−0.0114592−0.0291648
0.90.00826869−0.0417587−0.02726810.0145790−0.0549796
10.336636000.22328600.267491000.34493800.1716060
Table 4. The radial stress variation σrr under photothermal models.
Table 4. The radial stress variation σrr under photothermal models.
r
CPTEPLSPGN-IIPGN-IIIMGTPT
0−0.000656461−0.00393592−0.001913720.00018923−0.00395531
0.1−0.000756463−0.00493594−0.003913780.000389237−0.00495535
0.2−0.00442793−0.0053169−0.00514714−0.00469574−0.00517990
0.3−0.0111599−0.00615662−0.00762273−0.013776−0.00569276
0.4−0.0219495−0.00790313−0.0121283−0.0277876−0.00682064
0.5−0.0383551−0.0115677−0.0201443−0.0480931−0.00936626
0.6−0.0626674−0.0194178−0.0343838−0.0765139−0.0152988
0.7−0.0980382−0.0364582−0.0596874−0.1152370−0.0293835
0.8−0.1472930−0.0725051−0.1033440−0.1652930−0.0619620
0.9−0.1954270−0.1312310−0.1600040−0.2086930−0.1202020
100000
Table 5. The hoop stress variation σϑϑ under photothermal models.
Table 5. The hoop stress variation σϑϑ under photothermal models.
r CPTEPLSPGN-IIPGN-IIIMGTPT
00.001082703−0.00406557−0.0020835010.00302438−0.003656021
0.10.00182730−0.00465574−0.002835810.00402431−0.00465609
0.2−0.00136573−0.00486959−0.0037108−0.000561801−0.00470202
0.3−0.00723706−0.00532409−0.00546289−0.00877838−0.00478966
0.4−0.0166821−0.00623259−0.00864235−0.0215143−0.00494770
0.5−0.0311049−0.00807528−0.0142834−0.0400714−0.00525571
0.6−0.0525984−0.0119419−0.0242999−0.0662255−0.00594862
0.7−0.0843308−0.0204470−0.0423518−0.1024190−0.00784639
0.8−0.1325300−0.0412831−0.0768912−0.1534260−0.0150706
0.9−0.2258680−0.1120510−0.1633210−0.2450880−0.0599187
1−0.6360190−0.5511620−0.5971990−0.6420980−0.4730900
Table 6. The Maxwell’s stress τrr under different photothermal models.
Table 6. The Maxwell’s stress τrr under different photothermal models.
r CPTEPLSPGN-IIPGN-IIIMGTPT
0−0.00267372−0.000424817−0.00092268−0.00328662−0.000160729
0.1−0.00367374−0.000736172−0.00140181−0.00398662−0.00032458
0.2−0.00428631−0.00127631−0.00213414−0.00447240−0.000653511
0.3−0.00497257−0.00220043−0.00323770−0.00497870−0.00130324
0.4−0.00564912−0.00372845−0.00483778−0.00540708−0.00254364
0.5−0.00599225−0.00605572−0.00692594−0.00541727−0.00474740
0.6−0.00485913−0.00885152−0.00878874−0.00385297−0.00803334
0.70.00167694−0.00919275−0.006918230.00320529−0.01037210
0.80.027010100.006135150.012338700.029062100.000327701
0.90.116790000.08812070.098234500.118982000.076368400
10.426786000.4217000.423225000.427323000.420114000
Table 7. The variation of different photothermal fields with different rotation parameter Ω.
Table 7. The variation of different photothermal fields with different rotation parameter Ω.
r   =   0.9 . Ω θ u N σ r r σ z z τ r r
CPTE00.5799610.003014770.002452210−0.207303−0.2330910.1199820
50.5046340.003578270.001947390−0.181736−0.2028360.1086730
100.4324170.003935900.000754857−0.157142−0.1730950.0969131
PLS00.3805520.00839796−0.03440070−0.137101−0.1177670.1044720
50.3295050.00822557−0.02989960−0.120069−0.1020620.0947231
100.2806480.00789131−0.02610120−0.103682−0.0864960.0845178
PGN-II00.4390590.00712821−0.02102290−0.157300−0.1541780.1099490
50.3807770.00714752−0.01877670−0.137766−0.1333740.0996657
100.3249840.00699024−0.01693120−0.118985−0.1131270.0889246
PGN-III00.6215300.001568960.01446110−0.222413−0.2614530.1222750
50.5414390.002303690.01131410−0.195131−0.2266840.1107220
100.4645670.002826720.00831416−0.168858−0.1933050.0987226
MGTPT00.3447230.00902147−0.0437438−0.124892−0.0943120.100625
50.2980950.00874568−0.0375427−0.109368−0.0820450.0912424
100.2534940.00831745−0.0323135−0.0944313−0.0695830.0814059
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.; Askar, S. Thermo-Optical Mechanical Waves in a Rotating Solid Semiconductor Sphere Using the Improved Green–Naghdi III Model. Mathematics 2021, 9, 2902. https://doi.org/10.3390/math9222902

AMA Style

Abouelregal AE, Marin M, Askar S. Thermo-Optical Mechanical Waves in a Rotating Solid Semiconductor Sphere Using the Improved Green–Naghdi III Model. Mathematics. 2021; 9(22):2902. https://doi.org/10.3390/math9222902

Chicago/Turabian Style

Abouelregal, Ahmed E., Marin Marin, and Sameh Askar. 2021. "Thermo-Optical Mechanical Waves in a Rotating Solid Semiconductor Sphere Using the Improved Green–Naghdi III Model" Mathematics 9, no. 22: 2902. https://doi.org/10.3390/math9222902

APA Style

Abouelregal, A. E., Marin, M., & Askar, S. (2021). Thermo-Optical Mechanical Waves in a Rotating Solid Semiconductor Sphere Using the Improved Green–Naghdi III Model. Mathematics, 9(22), 2902. https://doi.org/10.3390/math9222902

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