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:
Relationships between constitutive and strain-displacement:
If the medium rotates at an angular velocity
uniformly,
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
is an extra initial term, while the second term is related to time-varying motion and Corioli’s acceleration
. Then, in a rotating reference frame, the equation of motion may be expressed as [
25]:
The corresponding plasma-thermal-elastic wave equation for the increase in carrier density
is as follows [
34,
35,
36]
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
to the heat flow vector
as
The thermal relaxation parameter (single delay) is known as the material-dependent constant, and it is an important element of this research. In terms of physics, the parameter 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]
where
represents a thermal displacement that fulfills
, and
signify the rates of the thermal conductivity. The equation of balance energy may be expressed as [
39,
40]
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]
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 . 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
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
, is
Incorporating Equation (10) into Equation (7) yields the proposed heat transfer equation, which describes the interplay of thermal, plasma, and elastic waves
We assumed that an initial magnetic field
pervades the neighboring free space. This generates an induced electro field
and an induced magnetic field
, both of which satisfy Maxwell’s magnetic equations and are adequate for slow-moving media.
where
is the magnetic permeability,
is the current density, and
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
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
with the values
,
, and
. Due to symmetry, all of the functions investigated are considered to be dependent on the distance
and the time
. The following are the components of the displacement vector in displacement-strain relationships:
The dilatation
may be calculated as follows:
The photothermal stresses (2) will be in the form
where
is the linear thermal expansion coefficient,
is the electronic deformation coefficient, and
,
are the Lame’s constants.
If the rotation is about the axial axis of the cylinder, i.e.,
, the equation of motion can be expressed in the following terms: the impact of the magnetic field
and the force due to the rotation of the body
. In this case, Equation (11) can be written as:
Assume the surface of the sphere is immersed in a magnetic field that acts in the
-direction of constant strength
. According to Equation (12), we obtain
The Lorentz force
can be expressed as
Thus, we have
and Maxwell’s stress
from Equations (18) and (19) as
After plugging Equations (16) and (20) into Equation (17), we get
where
.
Equation (21) can be rewritten as
The modified photothermal MGTE heat Equation (11) is as follows in the absence of any heat sources (
):
where
.
The non-dimensional variables of the partial differential equation are easy to construct. As a result, the following non-dimensional quantities are described:
If the primes are omitted, the system of equations may be expressed by the following formula:
where
The initial conditions can be assumed to be
We also assume that the following boundary conditions are hold:
where
is a constant,
is the surface recombination speed, and
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
, often known as
or by
We can derive the following findings by using the Laplace transform method to the governing Equations (25)–(28):
where
.
When we decouple Equations (35) and (37), we obtain
where
,
and
are defined as follows:
If the parameters
,
, and
are the roots of the equation
then Equation (40) can be expressed as
The roots of Equation (42) can be calculated as
with
Equation (43) has a general bounded solution that is can be expressed as
In Equation (46), the second type of modified Bessel function of order
is indicated by
. Three integral parameters are dependent on
:
are assumed. Moreover,
is related to two different variables,
and
, by plugging in Equation (46) into Equations (35)–(40) to get
The following well-known Bessel function can be used to get the displacement
:
According to the following relationships, for each positive
, modified Bessel
meets the following relatives:
If we insert relation (50) to Equations (46) and (49), then we get
This leads to a closed form of the final solutions for thermal stresses:
where
.
When we consider the non-dimensional Maxwell’s stress
, we get the following:
After applying the Laplace transform, the boundary conditions (31)–(33) take on the following forms.
Equation (56) is replaced with Equations (46) and (53) to produce
To obtain the parameter values
, 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:
where
stands for the number of terms,
stands for the real part, and
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]:
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 and 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 (). When the Green and Naghdi parameter , 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 . It is possible to use the photothermal Green and Naghdi model (PGN-II) if we remove the term that contains the parameter in the heat Equation (11). When the parameters , 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
,
and
.
First, we will review several key findings that apply to all distributions.
As seen in the graphs and tables, the thermal parameters 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.
for all photo-thermo-elasticity models across a wide range of
. 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
across radius
are compared in
Figure 2 and
Table 2. The numerical data and graphs show that the displacement
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
against the radial distance
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
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
than the MGTPT, with a similar behavior to the PLS and PGN-II models.
The fluctuation of the radial stress
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
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
fluctuation vs.
over a wide range of
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
, which is
. 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
, while the nonrotating case is chosen with
.
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
, and Maxwell’s stress
. 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
,
, and
decrease. Furthermore, as shown in
Figure 9,
Figure 10 and
Figure 11, the rotational parameter
has a great influence on the carrier density
and the thermal stresses
and
.
Figure 10 illustrates that increasing the amount of
increases the carrier density numerical values.
Figure 10 shows the change in thermal stress
as a function of radius
and parameter
. We can see that the rotational parameter
has a significant impact on the radial stress profile
. As the parameter
is increased, the stress
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.