The combination of the “electronic deformation” in the semiconducting medium, based on the photogeneration theory in the crystal lattice, and the “thermoelastic mechanism” owing to the integral photothermal process characterizes an important thermoelastic fact. When atoms are perturbed, they migrate from a high-density to a low-density area. The mechanical defects and internal strains make this form of transportation conceivable. When temperatures climb to the point where materials begin to melt, atomic spacing increases. In a similar vein, moisture transfer takes place when an existing concentration gradient forms as a result of the presence of unequal amounts of moisture. Temperature and humidity levels of the substance change depending on location and time. In this respect, the theories of heat transmission and moisture transfer are essentially identical. Mechanically applied stressors have the potential to significantly alter the distribution of temperature and moisture. That’s why it is important to determine exactly how mechanical deformation relates to diffusion due to changes in humidity and temperature. Many practical engineering problems involve the interplay of moisture, heat, and deformation. When heat and moisture act on solids, a phenomenon known as hygrothermoelasticity takes place. Szekeres (2000); Szekeres (2012) examined how generalized heat transmission is affected by moisture. It has been suggested by Gasch et al. (2016) that changes in humidity and temperature could do more damage than mechanical loadings. Using a fundamental comparison between heat and moisture, (Szekeres and Engelbrecht, 2000) formulated equations governing coupled hygrothermoelasticity.

Photothermal (PT) analysis of semiconductor materials’ thermal and electrical properties has grown in popularity. Semiconductors, used in sensors, solar cells, and advanced medical devices, have been extensively studied. Most renewable energy generation requires semiconductor knowledge. Semiconductors are dielectric and non-conductive. When optical energy activates semiconductor surface holes and electrons, electronic deformation occurs. The optical energy of light accelerates excited electrons to the surface, creating an electron cloud-like convective density or plasma waves. Photo-excitation and heat effect cause thermoelastic deformation. Thus, semiconductors are studied using thermoelasticity and photothermal theories. Numerous authors have developed novel methods to study photoacoustic spectroscopy’s sensitivity when a laser beam hits a semiconductor sample (Gordon et al., 1964; Kreuzer, 1971). Photothermal approaches in several physical tests verified the nano-composite semiconductor materials’ real temperatures, inner displacements, thermal diffusion, and other electrical properties (Tam, 1983; Tam, 1986; Tam, 1989; Todorovic et al., 1999; Song et al., 2010). Photo-excitation directly affects electronic deformation by causing elastic oscillations in the material’s internal structures. Hobiny and Abbas (2016) studied photothermal waves in an unbounded medium using a semiconductor-filled cylindrical cavity. Two-temperature theory and strain-stresses in a semiconducting material under photothermal waves with hydrostatic starting stress are problematic (Abo-dahab and Lotfy, 2017; Lotfy, 2017). Lotfy (2016) applied the photo-thermoelasticity hypothesis to semiconductor media thermal-plasma-elastic waves under a magnetic field and rotation.

Alterations to a material’s mechanical and thermal characteristics are often seen after a temperature change. The influence of the temperature gradient has been taken into consideration by several researchers after being ignored in previous studies of the expanded theory of thermoelasticity (Hasselman and Heller, 1980a; Youssef and Abbas, 2007a; Abouelregal and Marin, 2020; Marin et al., 2020). It is crucial to consider how temperature affects material characteristics (Marin et al., 2015) because the qualities of a material cannot be considered to have constant values under the effect of temperature change. The thermoelasticity theory emerged from discussions of coupled and uncoupled theories, both of which conflict with physical experiments. Biot (1956) introduced the unique hypothesis of linked thermoelasticity, which explained the anomaly. Infinite-speed thermal wave propagation was investigated using the CD theory. Lord and Shulman (1967) (LS) introduced a new concept when the heat equation was still in its infancy, by adding one relaxation period. Two relaxation periods were included in the heat conduction equation by Green and Lindsay (1972) (GL). The generalized thermoelasticity theory (GL) has been used by several writers (Chandrasekharaiah, 1986; Hosseini et al., 2013). Applications of the extended thermoelasticity theory including the interaction of thermal waves, electromagnetic fields, and mechanical waves in a thermoelastic solid medium are many.

Eringen (1972); Eringen and Edelen (1972) found that strain is defined as the gradient of the inner products of non-local deformations on non-local and local continuum components. This is a logical conclusion of the continuum theory of physics. Non-local thermoelastic theory (Eringen and Edelen, 1972) bridges the theoretical gap between the lattice-building theory and the classical continuum hypothesis. The classical continuum hypothesis allows for the investigation of constitutive relations between the atomic structure of lattices and the propagation of phonons. According to photomechanical waves and a moisture diffusivity process, (Alhashash et al., 2022) investigated the two-temperature thermoelasticity theory of a semiconductor model.

Physical characteristics of materials have been shown to depend on the temperature in recent studies. Deformation and thermo-mechanical behavior are both affected by the temperature dependency of these materials. The temperature-dependent thermal conductivity has a greater impact on thermal and mechanical behavior. Therefore, the thermal stress analysis is affected by the varying thermal conductivity. A thermoelastic hollow cylinder was an issue that Suhara investigated around the turn of the last century. Since then, several elastic and inelastic media issues have included the medium’s temperature-dependent physical features (Hasselman and Heller, 1980b). Youssef (2005) used the state-space technique to the solution of the generalized thermoelasticity issue involving a spherical cavity where the thermal conductivity is temperature-dependent under ramp-type heating. Youssef and El-Bary (2006) have included variable thermal conductivity in the generalized thermoelasticity with thermal shock layers. Youssef and Abbas (2007b) have solved the generalized thermoelasticity issue for an infinitely long annular cylinder with variable thermal conductivity. Nano-energy uses rely heavily on the effects of size. For the first time, a new gradient theory is presented to characterize the nonlocal equation of motion (elastic nonlocality) in nano-scale materials during the coupled between electrons and thermoelastic fields (Allen, 2014; Khamis et al., 2021). On the other hand, the nanoscale systems are very complex, therefore, we contented ourselves with studying the non-local equation of motion only.

In this research, the moisture diffusivity property is used to study the wave propagation in a non-local semiconductor medium with photo-thermoelasticity. The work is investigated in one dimension (1D) while a photothermal transport process is going on, affected by specific mechanical forces and the diffusivity upon being hit by a non-Gaussian laser pulse. The problem has been addressed at the free surface of a non-local semiconducting material when the thermal conductivity of the medium depends on the heat. By applying Laplace transforms to the time variable, the governing equations for the most important fundamental physical quantities may be solved analytically. Numerical inversion is carried out on a computer using powerful and efficient software. As a final step, normal force stress, normal displacement, carrier density, temperature distribution, and moisture concentration were all calculated numerically with the different values of thermal memory, thermal conductivity, moisture reference and nonlocality parameters.

The carrier density N r k , t (optoelectronic), moisture concentration m r k , t , temperature change T r k , t , and displacement distribution u r k , t are all used to describe the novel model developed in this study ( r k is the position and t represents the time). The non-local medium’s plasma-thermal-elastic wave equations interact with the moisture diffusion equations when the thermal conductivity is variable in the following tensor form (Sladek et al., 2020): ∂ N r i , t ∂ t = D E N , i i r i , t − N r i , t τ + κ τ T r i , t (1) D T 1 + τ θ ∂ ∂ t K m T , i , i + D T m K m m , i , i = 1 + τ q ∂ ∂ t ρ   C e ∂ T r i , t ∂ t + γ t T 0 ∂ u i , j r i , t ∂ t − E g τ N r i , t (2) D m   K m m , i , i + D m T K m T , i , i = ∂ ∂ t K m m r i , t − E g τ N r i , t + γ m m 0 D m ∂ u j , j r i , t ∂ t (3)

The non-local motion equation for diffusivity semiconductor medium may be written as: ρ 1 − ξ 2 ∇ 2 ∂ 2 u i r i , t ∂ t 2 = σ i j , j ′ (4)

The nonlocal elastic parameter in length is denoted by ξ . The strain in terms of displacement takes the form: ε i j = 1 2 u i , j + u j , i (5)

The stress equation in terms of plasma, temperature, strain and moisture concentration can be written as follows: 1 − ξ 2 ∇ 2 σ i j = σ i j ′ = C i j k l ε k l − β i j α t 1 + τ θ ∂ ∂ t T + d n N − β i j m m ,   i , j , k , l = 1,2,3 (6)

In all equations in this work a “comma” before an index implies space-differentiation and a “dot” above a symbol refers to the time-differentiation.

Take into account the moisture thermal conductivity, which is variable and may be chosen as a linear function of temperature. Non-local semiconductor material’s changing moisture thermal conductivity under the effect of a light heat source is represented as a function of temperature as shown in (Hasselman and Heller, 1980b; Alhashash et al., 2022): K m = K 0 1 + π T (7)

When the medium is independent of temperature, the constant thermal conductivity (reference) is K 0 and π represents a negative small parameter. The integral form of thermal conductivity may be derived using the Kirchhoff transform to convert the nonlinear temperature components into linear ones (Hasselman and Heller, 1980b; Alhashash et al., 2022): ϑ = 1 K 0 ∫ 0 T K m z d z (8)

Considering that the values of all physical quantities are unrelated to the y z -coordinates, all studies are performed along the x -axis (the direction of wave propagation is along the x -axis). The main Eqs. 1, 4 involving the physical variables may be recast in 1D as follows (Hasselman and Heller, 1980a): ∂ N ∂ t = D E ∂ 2 N ∂ x 2 − N τ + κ τ T (9)

The motion Eq. 4 has the following structure (Hasselman and Heller, 1980b): ρ   1 − ξ 2 ∂ 2 ∂ x 2 ∂ 2 u ∂ t 2 = 2 μ + λ ∂ 2 u ∂ x 2 − γ t 1 + τ θ ∂ ∂ t ∂ T ∂ x − δ n ∂ N ∂ x − γ m ∂ m ∂ x (10)

Where γ t , m = β   α m , T and δ n = β d n , β = 3 μ + 2 λ .

Several methods of differentiation provide Eqs. 7, 8 for the map transform, which may be used to the fundamental equations to derive the thermal conductivity, which is a variable in the original computations, as: K 0 ϑ , i = K m T , i , K 0 ϑ , i i = K m T , i , i , K 0 / K m ϑ , i = T , i . (11)

Using the same technique, the time-differentiating from the first order for diffusivity is: K 0 ∂ Μ ∂ t = K m ∂ m ∂ t (12)

As a result applying the map transform to Eq. 1, which ∂ ∂ x i may act on both sides, yields: ∂ ∂ t ∂ N ∂ x j = D E ∂ N , i i ∂ x j − 1 τ ∂ N ∂ x j + κ τ   ∂ T ∂ x j , ∂ ∂ t ∂ N ∂ x j = D E ∂ N , i i ∂ x j − 1 τ ∂ N ∂ x j + κ K 0 τ K m   ∂ ϑ ∂ x j , ∂ ∂ t ∂ N ∂ x j = D E ∂ N , i i ∂ x j − 1 τ ∂ N ∂ x j + κ τ   ∂ ϑ ∂ x j . (13)

The last term in the first half of Eq. 13 in the preceding equation may be extended as follows if the non-linear components are disregarded: κ K 0 τ K m ∂ ϑ ∂ x = κ K 0 K 0 1 + π T ∂ ϑ ∂ x = κ 1 + π T − 1 ∂ ϑ ∂ x = κ τ 1 − π T + π T 2 − . . . . . . . ) ∂ ϑ ∂ x = κ τ ∂ ϑ ∂ x − κ τ π T ∂ ϑ ∂ x + κ τ π T 2 ∂ ϑ ∂ x − . . . . . . . = κ τ ∂ ϑ ∂ x (14)

Leibniz integral rule is applied to Eq. 13, which yields: ∂ N ∂ t = D E N , i i − 1 τ N + κ τ ϑ . (15)

Rewriting the heat Eq. 2 and moisture diffusion Eq. 3 after using the map transform yields the following: D T 1 + τ θ ∂ ∂ t ϑ , i i + D T m M , i i = 1 + τ q ∂ ∂ t 1 k ∂ ϑ ∂ t + γ t T 0 K 0 ∂ u i , j ∂ t − E g K 0 τ N (16) D m M , i i + D m T ϑ , i i = ∂ M ∂ t − E g K 0 τ N + γ m m 0 D m K 0 ∂ u i , j ∂ t (17)

According to the map transform, the equation of motion (4) takes the form: ρ   1 − ξ 2 ∂ 2 ∂ x 2 ∂ 2 u ∂ t 2 = 2 μ + λ ∂ 2 u ∂ x 2 − γ t 1 + τ θ ∂ ∂ t K 0 K m ∂ ϑ ∂ x − δ n ∂ N ∂ x − γ m K 0 K m ∂ M ∂ x (18)

With neglecting the non-linear term, yields: K 0 K m ∂ ϑ ∂ x = K 0 K 0 1 + π T ∂ ϑ ∂ x = 1 + π T − 1 ∂ ϑ ∂ x = 1 − π T + π T 2 − . . . . . . . ) ∂ ϑ ∂ x = ∂ ϑ ∂ x (19)

According to relation (19), the mapped motion Eq. 18 can be rewritten as: ρ   1 − ξ 2 ∂ 2 ∂ x 2 ∂ 2 u ∂ t 2 = 2 μ + λ ∂ 2 u ∂ x 2 − γ t 1 + τ θ ∂ ∂ t ∂ ϑ ∂ x − δ n ∂ N ∂ x − γ m ∂ M ∂ x (20)

On the other hand, the constitutive equation for the non-local medium can be rewritten in 1D as: 1 − ξ 2 ∂ 2 ∂ x 2 σ x x = 2 μ + λ ∂ u ∂ x − β α t ϑ + d n N − γ m M = σ (21)

The following are additional non-dimensional values we may provide for use in simplifying expressions: x ′ , ξ ′ , u ′ = x , ξ , u C T t * , t ′ , τ q ′ , τ θ ′ = t , τ q , τ θ t * , ϑ ′ , N ′ = γ t ϑ , δ n N 2 μ + λ , σ ′ = σ μ , M ′ = M (22)

Using the dimensionless Eq. 22, the mapped Eqs. 15–17 and Eqs. 20, 21 can be reduced in the following form: ∂ 2 ∂ x 2 − q 1 − q 2 ∂ ∂ t N + ε 3   ϑ = 0 (23) 1 + τ θ ∂ ∂ t ∂ 2 ∂ x 2 − a 1 1 + τ q ∂ ∂ t ∂ ∂ t ϑ + a 2 ∂ 2 M ∂ x 2 + a 3   N − ε 1 1 + τ q ∂ ∂ t ∂ 2 u ∂ t ∂ x = 0 (24) ∂ 2 ∂ x 2 − a 4 ∂ ∂ t M + a 5 ∂ 2 ϑ ∂ x 2 + a 6   N − a 7 ∂ 2 u ∂ t ∂ x = 0 (25) ∂ 2 ∂ x 2 − 1 − ξ 2 ∂ 2 ∂ x 2 ∂ 2 ∂ t 2 u − 1 + τ θ ∂ ∂ t ∂ ϑ ∂ x − ∂ N ∂ x − a 8 ∂ M ∂ x = 0 (26) σ = a 9 ∂ u ∂ x − ϑ + N − a 10 M (27)

Where q 1 = k t * D E ρ τ C e ,   q 2 = k D E ρ C e ,   a 1 = C T 2 t * k D T ,   a 2 = D T m γ t D T 2 μ + λ , ε 2 = α T E g t * τ d n ρ C e ,   a 3 = ε 2 a 1 ε 1 = γ t 2 T 0 t * K 0 ρ ,   a 4 = C T 2 t * D m ,   a 5 = D m T 2 μ + λ D m γ t ,   a 6 = E g 2 μ + λ t * a 4 K 0 δ n τ ,   a 7 = γ m m 0 C T 2 t * K 0 ,   a 8 = γ m 2 μ + λ ε 3 = d n K 0 κ t * α T ρ τ C e D E ,   a 9 = 2 μ + λ μ ,   a 10 = γ m μ C T 2 = 2 μ + λ ρ ,   δ n = 2 μ + 3 λ d n ,   D T = k ρ C e ,   t * = k ρ C e C T 2

To aid in finding a mathematical solution when the system is initially at rest, below are the initial conditions: u x , t t = 0 = ∂ u x , t ∂ t t = 0 = 0 ,   T x , t t = 0 = ∂ T x , t ∂ t t = 0 = 0 , m x , t t = 0 = ∂ m x , t ∂ t t = 0 = 0   σ x , t t = 0 = ∂ σ x , t ∂ t t = 0 = 0 , N x , t t = 0 = ∂ N x , t ∂ t t = 0 = 0 (28)

To simplify the definition of partial differential equations, the Laplace transform is used, which is defined as: L Χ x , t = Χ ¯ x , s = ∫ 0 ∞ exp − s t Χ x , t   d   t (29)

Eq. 29 is used to convert Eqs. 23–27 as follows: D 2 − α 1 N ¯ + ε 3   ϑ ¯ = 0 (30) D 2 − α 2   ϑ ¯ + a 2 ′ D 2   M ¯ + a 3 ′   N ¯ − α 3   D u ¯ = 0 (31) D 2 − α 4   M ¯ + a 5 D 2 ϑ ¯ + a 6   N ¯ − α 5   D u ¯ = 0 (32) D 2 − Q u ¯ − Z D (   ℘ ϑ ¯ + N − R D M ¯ = 0 (33) σ ¯ = a 9 D u ¯ − ϑ ¯ + N ¯ − a 10 M ¯ (34) where, D = d d x , R = a 8 1 + ξ 2 s 2 , Z = 1 1 + ξ 2 s 2 , α 1 = q 1 + q 2 s , Q = s 2 1 + ξ 2 s 2 , α 2 = a 1 1 + τ q s s 1 + τ θ s , a 2 ′ = a 2 1 + τ θ s , α 4 = a 4 s , ℘ = 1 + τ θ s , α 5 = a 7 s , α 3 = ε 1 1 + τ q s s 1 + τ θ s a 3 ′ = a 3 1 + τ θ s .

Eliminating technique is used for the quantities ϑ ¯ , u ¯ , N ¯ and M ¯ using Mathematica softwear; yields D 8 − Θ 1 D 6 + Θ 2 D 4 − Θ 3 D 2 − Θ 4 M ¯ , N ¯ , ϑ ¯ , u ¯ x , s = 0 (35)

By conducting calculations using computer programs such as Mathematica, the basic coefficients can be calculated in Eq. 35, whose values are deduced as follows: Θ 1 = − − Q a 2 a 5 − a 2 ′ a 5 α 1 − a 5 R α 3 + Q − a 2 α 5 + R α 5 + Z ℘ α 1 + α 2 + α 3 + α 4 a 3 ′ a 2 ′ − 1 , Θ 2 = Q a 2 ′ a 5 α 1 + a 5 R α 1 α 3 − Q α 1 + α 2 + α 4 − a 2 ′ a 7 ε 3 + a 2 ′ α 1 α 5 + a 2 α 5 ε 3 − R α 1 α 5 + α 2 α 5 + Z a 3 ε 3 − α 1 α 2 − α 1 α 3 − Z ℘ α 1 α 4 − α 2 α 4 − α 3 α 4 − α 3 ε 3 a 3 ′ a 2 ′ − 1 , Θ 3 = − Q a 2 ′ a 7 ε 3 − a 3 ′ ε 3 + Z α 1 α 2 + α 1 α 4 + α 2 α 4 − a 3 ′ R α 5 ε 3 + a 7 R α 3 ε 3 + R α 1 α 2 α 5 − a 3 ′ α 4 ε 3 + α 1 α 2 α 4 + Z ℘ α 1 α 3 α 4 + α 3 α 4 ε 3 a 3 ′ a 2 ′ − 1 , Θ 4 = Q a 3 ′ α 4 ε 3 − Z ℘ α 1 α 2 α 4 a 3 ′ a 2 ′ − 1 . (36)

Eq. 35 can be solved using the factorization method, which yields: D 2 − m 1 2 D 2 − m 2 2   D 2 − m 3 2 D 2 − m 4 2 ϑ ¯ , u ¯ , N ¯ , M ¯ x , s = 0 (37)

To acquire roots in the positive real part at x → ∞ , substitute m i 2 i = 1,2,3,4 . Since this is a linear problem, the solution to Eq. 37 takes the form below. ϑ ¯ x , s = ∑ i = 1 4 D i s   e − m i x (38)

It is possible to express the solutions in terms of the other variables as follows: N ¯ x , s = ∑ i = 1 4 D i ′ s   e − m i x = ∑ i = 1 4 H 1 i D i s   e − m i x (39) u ¯ x , s = ∑ i = 1 4 D i ″ s   exp − m i x = ∑ i = 1 4 H 2 i   D i s   exp − m i x (40) M ¯ x , s = ∑ i = 1 4 D i ‴ s   exp − m i x = ∑ i = 1 4 H 3 i   D i s   exp − m i x (41) σ ¯   x , s = ∑ i = 1 4 D i 4 s   exp − m i x = ∑ i = 1 4 H 4 i   D i s   exp − m i x (42)

The quantities D i , D i ′ , D i ″ , and D i ‴   i = 1 . . 4 are the unknown parameters. The following relationship exists between the unknown parameters and may be determined using the basic Eqs 30–34: H 1 i = − ε 3 m i 2 − α 1 , H 2 i = a 2 − R m i 4 + − a 2 ′ α 1 − Z a 2 ′ ε 3 + R α 1 + α 2 m i 2 + R a 3 ′ ε 3 − α 1 α 2 m i ( m i 4 a 2 ′ + − Q a 2 ′ − Z a 2 ′ α 1 − R α 3 m i 2 + Q a 2 ′ α 1 + R α 1 α 3 , H 3 i = − m i 6 + − Q − Z α 1 − α 2 − α 3 m i 4 + ( Q α 1 + α 2 − a 3 ′ ε 3 + Z α 1 α 2 + α 3 α 1 + ε 3 ) m i 2 + Q a 3 ′ ε 3 − α 1 α 2 m i 2 − α 1 m i 2 a 2 ′ − Q a 2 ′ − R α 3 m i 2 , H 4 i = a 9 m i H 2 i − 1 + H 1 i − a 10 H 3 i

Assume that the elastic non-local semiconductor medium is exposed to thermal loads of the thermal ramp type under conditions of heat, plasma, and mechanical loading. These conditions are applied to the free outside surface of the non-local medium. In every scenario, Laplace transformations are used (Vlase et al., 2017; Sladek et al., 2020).

(I) Type ramp heating is used for the thermal boundary condition at x = 0 , which may be expressed according to the thermal excitation function F t :

The time of pulse heat flux is t 0 = 0.001 , Equation 43 is solved by the use of Laplace transforms and the attribute of being dimensionless, whereby: ∑ n = 1 4 D i s = F ¯ s 1 − e − s t 0 t 0 s 2 (44)

(II) A mechanical load is applied to the surface at x = 0 . This means that for a generalized non-local photo-thermoelastic medium with an orthotropic boundary, the boundary conditions are provided according to a given hydrostatic initial stress function q t by:

Laplace transformation is used, which yields: ∑ i = 1 4 a 9 m i H 2 i − 1 + H 1 i − a 10 H 3 i D i = q ¯ s (46)

(III) During the diffusion phase, the carriers’ plasma density, which can be quantified with the help of the plasma excitation function N * , may reach the sample’s surface at x = 0 due to the low chance of recombination. In this case, it is assumed that a laser excitation is supplied to a non-local semiconductor surface in the form of a non-Gaussian plasma shock loading as: N x , s = N * = N 0 * t   exp − t t p (47)

Maximum plasma shock loading caused by a laser pulse occurs at a time t p = t , denoted by N 0 * . According to Laplace transform, yields: N ¯ 0 , s = ∑ i = 1 4 H 1 i D i 0 , s = N 0 * s + 1 / t p 2 (48)

(Ⅳ) The boundary condition for displacement on the free surface at x = 0 is: u ¯ 0 , s = 0 (49)

Or: ∑ i = 1 4 H 2 i D i x , s = 0 (50)

To obtain the relation between the temperature T and the mapped temperature ϑ , using the map and integral Eqs. 7, 8, yields: ϑ = 1 K 0 ∫ 0 T K 0 1 + π T d T = T + π 2 T 2 = π 2 T + 1 π 2 − 1 2 π (51) T = 1 π 1 + 2 π ϑ − 1 (52)

In the same way, moisture diffusivity can be obtained.

The inverse Laplace-transform is calculated using the Riemann-sum approximation method; [see ref. (Honig and Hirdes, 1984)] for additional information. Therefore, the following expressions may be used to derive the field’s time-domain variables based on Fourier series expansion. For a Laplace-domain function ζ ¯ x , s , the transfer to the time domain looks like this: ζ x , t ′ = L − 1 ζ ¯ x , s = 1 2 π i ∫ n − i ∞ n + i ∞ ⁡ exp s t ′ ζ ¯ x , s d s (53) where s = n + i Μ ( n , Μ ∈ R ). Consequently, the equation with its variables changed around may be expressed as: ζ x , t ′ = exp n t ′ 2 π ∫ − ∞ ∞ ⁡ exp i β t ζ ¯ x , n + i β d β (54)

The following relationship may be found by expanding for the function ζ x , t ′ in the closed interval 0 , 2 t ′ using the Fourier series: ζ x , t ′ = e n t ′ t ′ 1 2 ζ ¯ x , n + R e ∑ k = 1 R ζ ¯ x , n + i k π t ′ − 1 n (55)

Where i and R e indicate the imaginary number unit and the real component, respectively, adequate R stands for the number of terms in the equation. The value of the parameter n may be determined based on the instructions in (Honig and Hirdes, 1984). Because of this, a method known as numerical inversion of the Laplace transform (NILT) is used (Brancik, 1999). This approach is based on the quick Fourier transform.

To further demonstrate the theoretical findings from the preceding part, we will now provide numerical data. Numerical values of physical characteristics such as thermal, displacement, plasma, moisture concentration, and the normal distribution of stress obtained from this problem under the impact of variable thermal conductivity over a short time may now be graphically represented in MATLAB for the researched physical fields. In the theoretical investigation, n-type silicon (Si) acts as the non-local semiconductor solid material in the device. Table 1 uses SI units for the following Si physical parameters: (Lotfy et al., 2017; Ezzat, 2020; Mondal and Sur, 2021; Zhao et al., 2022a; Zhao et al., 2022b).

In this study, we investigate how several external conditions, such as ramp-type heating, laser shock, moisture, and mechanical forces, might affect the propagation of photothermal-elastic waves in a solid non-local semiconducting medium. Electron-elastic deformations ignite elastic waves inside the non-local semiconductor, resulting in a fluctuating deformation potential. One may get the 1D case now. There is hope that the photo-thermoelasticity theories may focus attention on the peculiarities of wave motion in non-local semiconductor media. Considerable graphical analysis has been performed to investigate the significance of thermal memory, the variable thermal conductivity effect, nonlocality, and the moisture reference. As a result of its clarity and precision, the photothermal theory may be used to explain de-excitation in materials and the absorption of light. The plots show that the thermal memory, thermal conductivity, moisture content, and nonlocality of the medium have a major impact on all the considered domains. The research might help scientists learn more about how waves behave in a wide range of environments and temperatures. There is a clear correlation between the thermal conductivity, moisture content, and nonlocality of the field quantities, as the amplitude of these values varies (increasing or decreasing) with the thermal memory. Recent research has shown that semiconductors may be utilized to convert solar energy into electricity while also withstanding exposure to laser pulses, demonstrating the vital role that semiconductors play in today’s cutting-edge technologies. Modern technology relies heavily on semiconductors in a variety of devices, such as solar cells, displays, and transistors. Several fields of mechanical and electrical engineering employ them as nanomaterials. Thermomechanical, sensor, resonator, medical, and accelerometer researchers should all be able to make use of the study’s findings in their own ongoing and future investigations. At the opposite end of the spectrum, microwave and radio frequency emitters make it possible for people to communicate wirelessly. Visible and infrared diode lasers are at the center of the information technology industry. The method presented here has potential application to a wide range of photo-thermoelasticity and thermodynamic issues.