Next Article in Journal
Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being
Previous Article in Journal
On the Normalization of Interval Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation

Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Majmaah 11952, Saudi Arabia
Mathematics 2020, 8(11), 2094; https://doi.org/10.3390/math8112094
Submission received: 29 October 2020 / Revised: 18 November 2020 / Accepted: 19 November 2020 / Published: 23 November 2020
(This article belongs to the Section Mathematical Physics)

Abstract

:
In this work, the generalized theory with dual-phase-lag of thermoelasticity is employed to study the influences induced by absorbing a penetrating laser radiation inside a 2D thermoelastic semi-infinite medium. The medium’s surface is presumed to be exposed to temperature-dependent heat losses and is traction-free. The considered problem is solved using the integral transforms technique by applying the double-transformation Laplace and Hankel. A numerical fashion is applied to obtain the inverse of the Laplace transformation. The results of this problem are presented graphically for some studied fields.

1. Introduction

One of the celebrated contributions to the field of thermoelasticity is the coupled theory, which was formulated by Biot in 1956 [1], in which the author tackled one of the faults of the classical uncoupled theory of thermoelasticity by modifying the energy equation to contain a term representing the strain field. Since Biot’s model was based upon Fourier’s law, it failed to address the second defect of the classical uncoupled theory, which means that the generated heat waves still propagate with an infinite velocity. Hence, the behavior of Biot’s theory contradicts the laws of physics, especially for some problems such as utilizing laser heat sources, where their pulse duration is extremely short, and problems with large temperature gradients [2,3].
After Biot’s attempt, great efforts were made to obtain a generalization that addressed the defect found in the coupled theory; this point became the focus of attention for those interested in thermoelasticity. Authentic progress was made in 1967 when Lord and Shulman [4] formulated a new model based on Cattaneo’s approach [5,6] instead of Fourier’s Law. This theory was followed by several essential generalizations [7,8,9]; most of these theories were dependent on adjusting the heat equation. The DPL model or the generalized theory with dual-phase-lag is one of the essential generalizations developed by Tzou [10,11], where the energy equation was modified to contain two distinct phase-lags: one symbolizes the temperature gradient, while the other symbolizes the heat flux. Several authors have employed the DPL model to study thermoelastic waves under the effect of different fields [12,13,14,15,16,17,18,19].
Since the discovery of lasers, many applications in physics, mechanics, and engineering have been implemented based on their various properties, especially in materials processing, such as drilling of holes, spot welding, glazing of materials, and surface-hardening scribing [20]. Several authors have employed the generalized theories of thermoelasticity to study the thermoelastic interactions induced by a pulsed laser, considering different aspects [21,22,23,24]. Zenkour and Aboelregal [25] studied the fractional influences in a half-space solid exposed to non-Gaussian laser radiation by utilizing the two-temperature generalized model with a fractional-order energy equation. Bassiouny et al. [26] investigated the thermoelastic conduct of an elastic semi-infinite solid under the influence of a pulsed laser, by applying the theory of fractional-order strain thermoelasticity. Allam and Tayel [27] applied five different theories to study the interaction caused by surface absorption of a pulsed laser in a functionally graded thermoelastic semi-infinite medium. Abbas and Marin [28] employed the Lord and Shulman model to study the effect of the relaxation time in a 2D thermoelastic semi-infinite solid illuminated by a laser pulse. Abouelregal and Zenkour [29] employed the generalized two-temperature thermoelastic theory to discuss thermoelastic responses in a half-space, whose pounding plane was exposed to a laser pulse. Tayel and Hassan [30] and Tayel [31] studied the thermoelastic interaction induced by both surface and volumetric absorption techniques, respectively, in a one-dimensional semi-infinite medium, considering the absence and existence of cooling, by employing various thermoelasticity theories.
Integral transforms have been effectively employed for solving many problems in physics, applied mathematics, and engineering sciences. Although Laplace transform is considered one of the most essential integral transformations applied to solve coupling problems in the theory of thermoelasticity, its inverse is very difficult to obtain using the usual methods. This difficulty has prompted those interested in thermoelasticity to apply numerical methods or to use asymptotic expansion techniques that are adequate for short times. This technique was used in the context of the coupled theory by Hetnarski [32,33] and applied by several authors to some generalized thermoelastic problems [34,35,36].
The main goals of this work are to study the thermoelastic interaction induced by the volumetric absorption technique of a laser pulse in a 2D semi-infinite medium and to investigate the effects of the phase-lag parameters of the DPL theory of thermoelasticity on the existence and absence of a cooling effect. The medium is studied under the influence of a laser beam, the spatial and temporal profiles of which are considered Gaussian. The medium surface is exposed to temperature-dependent heat losses and is considered traction-free.

2. Problem Formulation

A thermoelastic, isotropic, and homogeneous 2D semi-infinite medium ( z 0 ) is considered under the influence of a laser beam that is incident-normal to a small specified region of the medium’s surface. The medium is initially at a uniform temperature T 0 and its surface ( z = 0 ) is taken as stress-free and assumed to be exposed to temperature-dependent cooling. The cylindrical coordinate system ( r ,   ϕ ,   z ) is utilized in the problem, with the z-axis oriented inward and vertically.
According to the heating process, which is considered axisymmetric, all quantities will depend on r , z and t only and, thus, the displacement vector u will possess the following components:
u r = u ( r ,   z ,   t ) ,   u ϕ = 0 ,   u z = w ( r ,   z ,   t ) .
Consequently, the nonvanishing strain components are
e r r =     u   r ,   e z z =   w   z ,   e ϕ ϕ =   u r   ,   e r z = 1 2 ( u z + w r ) .
Thus, the dilatation e takes the form
e = e r r + e ϕ ϕ + e z z = 1 r   r ( r   u ) +   w   z .
Following [10,11], the linear heat equation for a homogeneous and isotropic body in the context of the DPL model is given as
k   ( 1 + τ 1   t ) 2 T =   ( 1 + τ 2   t + τ 2 2 2 2   t 2 ) ( ρ   C E   T   t + γ   T 0   e   t Q ( r , z , t ) )   , with   2 ( ) = ( 2   r 2 + 1 r   r + 2   z 2 ) ( ) ,
where k stands for thermal conductivity; ρ for density; T for the absolute temperature, C E , for the specific heat at constant strain; T 0 for the reference temperature, which is chosen such that | ( T T 0 ) / T 0 | 1 and τ 1 , τ 2 for phase-lags of the temperature gradient and the heat flux, respectively, where ( τ 1 , τ 2 > 0 ) . Moreover, γ = ( 3 λ + 2 μ ) α T , in which λ and μ are lamé’s constants and α T is the thermal expansion.
For volumetric absorption in a 2D medium, Q   ( r , z , t ) is given as
Q   ( r , z , t ) = A 0 q 0 Q ( t )   e r 2 ω 2 β   e β   z   ,
where A 0 ,   q 0 , Q ( t ) , ω , β denote the transition coefficient of the illuminated surface, the maximum value of the laser intensity, the temporal laser pulse shape, a spatial constant, and the absorption coefficient of the material, respectively.
The axially symmetric equations of motion in terms of displacements and absence of body forces are
( λ + μ )   e   r + μ   ( 2 1 r 2 ) u γ T   r = ρ 2 u   t 2 ,
( λ + μ )     e   z + μ   2   w γ   T   z = ρ 2 w   t 2 .
Furthermore, the nonvanishing stresses are determined by
σ z z = 2 μ w z + λ e γ   θ ,
σ r r = 2 μ u r + λ e γ   θ ,
σ ϕ ϕ = 2 μ u r + λ e γ   θ ,
σ r z = μ (   u   z +   w   r ) ,
where θ = T T 0 represents the increment of the temperature.
Due to the foregoing formulation, the boundary conditions will be
σ z z | z = 0 = σ r z | z = 0 = 0 , k ( 1 + τ 1   t )     θ   z | z = 0 =   ( 1 + τ 2   t + τ 2 2 2 2   t 2 ) h   θ ( r , 0 , t ) , }
where h represents the coefficient of the heat transfer responsible for cooling, and θ ( r , 0 , t ) is the surface temperature.
Additionally, the initial conditions are
θ ( r , z , t ) | t = 0 = θ ( r , z , t ) t | t = 0 = 0 ,   u ( r , z , t ) | t = 0 = u ( r , z , t ) t | t = 0 = 0 , w ( r , z , t ) | t = 0 = w ( r , z , t ) t | t = 0 = 0 . }
Now, we introduce the following nondimensional variables
( u , w ,   r , z ,   ω ) = c   η ( u , w ,   r , z , ω ) θ = γ   θ ( λ + 2 μ ) σ i j = σ i j μ ( t , τ 1 ,   τ 2 ) = c 2 η ( t , τ 1 ,   τ 2 ) h = h c   ρ   C E β = β c   η q 0 = γ q 0 ( λ + 2 μ ) c   ρ   C E
where c 2 = λ + 2 μ ρ , η = ρ c E   k .
For seeking the simplicity primes will be dropped, and thus Equations (3), (5) and (6) in the nondimensional forms will be
( 1 + τ 1   t )   2 θ = ( 1 + τ 2 t + τ 2 2 2 2 t 2 ) ( θ t + ε   e t A 0 q 0 β   Q ( t ) e r 2 ω 2 e β z ) ,
( g     1 ) e r + ( 2   1 r 2 ) u g θ r = g 2 u t 2 ,
( g     1 ) e z + 2 w g θ z = g 2 w t 2 ,
where g = λ + 2 μ μ   , ε = T 0 γ 2 ρ   C E ( λ + 2 μ ) .
Combining Equations (15) and (16) by applying the operators 1 r (     r r . )   for (15) and z   for (16), then adding them, one gets
2 e 2 θ = 2 e   t 2 .
In the nondimensional form, the stress components become
σ r r = 2 u r + ( g 2 ) e g θ ,
σ z z = 2 w z + ( g 2 ) e g θ ,
σ ϕ ϕ = 2 u r + ( g 2 ) e g θ ,
σ r z = u z + w r .
In addition, Equation (11) will be written as
( σ z z = σ r z ) | z = 0 = 0 , ( 1 + τ 1   t )   θ   z | z = 0 = ( 1 + τ 2   t   +   τ 2 2 2   2   t 2 )   h   θ ( r , 0 , t ) . }

3. Problem Solution

The specific conditions in Equation (12) and the basic Equations (14) and (17) predict that the double-transformation Laplace and Hankel integral transforms [37] are more appropriate to be applied for t and r , respectively, in which the transformation is given as
f ¯ ˜ (   α ,   z ,   s ) = 0 0 r J n ( α   r ) f ( r , z , t )   e s t d r   d   t , f ( r , z , t ) = 1 2 π i   σ i   σ + i   0 α J n ( α   r ) f ¯ ˜ ( α , z , s ) e s t d   α   d   s , }
where s is the Laplace transform variable, the superscripts dash and tilde indicate the Laplace and Hankel transformation of f ( r , z , t ) , respectively, and J n is the Bessel function of the first kind of order n in which
n = { 0 for   f ¯ ( r ,   z ,   s ) =   { θ ¯ ,   w ¯ ,   e ¯ ,   σ ¯ z z ,   σ ¯ r r ,   σ ¯ ϕ ϕ } ( r ,   z ,   s ) 1 for   f ¯ ( r ,   z ,   s ) =   { u ¯ ,   σ ¯ r z } ( r ,   z ,   s )     .
Thus, Equations (14) and (17) become, respectively
( D 2 α 2 s   Γ ) θ ¯ ˜ ε s   Γ   e ¯ ˜ = A 0 q 0 β   Γ Q ¯ ( s ) η ˜ ( α ) e β z .
( D 2 α 2 s 2 ) e ¯ ˜ ( D 2 α 2 ) θ ¯ ˜ = 0 ,
where Γ = 1 + τ 2 s + τ 2 2 2 s 2 1 + τ 1 s , Q ¯ ( s ) is the Laplace transform of Q ( t ) , and η ˜ ( α ) = ω 2 2 e 1 4 α 2 ω 2 is the Hankel transform of order zero of e r 2 ω 2 .
The elimination of e ¯ ˜ between Equations (24) and (25) gives
( D 4 b 1 ( α , s ) D 2 + b 2 ( α , s ) ) θ ¯ ˜ = A 0 q 0 β   Γ ( β 2 α 2 s 2 ) Q ¯ ( s ) η ˜ ( α ) e β z ,
where b 1 ( α , s ) = 2 α 2 + s 2 + s Γ ( 1 + ε ) ,   b 2 ( α , s ) = α 4 + α 2 s 2 + s Γ α 2 ( 1 + ε )   + s 3 Γ .
Consequently, Equation (26) can be written as
( D 2 k 1 2 ) ( D 2 k 2 2 ) θ ¯ ˜ = A 0 q 0 β   Γ ( β 2 α 2 s 2 ) Q ¯ ( s ) η ˜ ( α ) e β z .
where k i 2 are the roots of the following auxiliary equation:
k 4 b 1 ( α , s ) k 2 + b 2 ( α , s ) = 0 .
The solution of Equation (27) after considering the behavior of θ ¯ ˜ as z is given as
θ ¯ ˜ ( α , z , s ) = i = 1 2 B i ( α , s )   e k i z + M ( α , s ) e β z ,
where M ( α , s ) represents the coefficient of the particular solution of Equation (27); using the undetermined coefficients method, one gets
M ( α , s ) = A 0 q 0 Γ β ( β 2 α 2 s 2 ) Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) ,
and thus
θ ¯ ˜ ( α , z , s ) = i = 1 2 B i ( α , s )   e k i z A 0 q 0 Γ β ( β 2 α 2 s 2 ) Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z .
In a similar fashion, the elimination of θ ¯ ˜ between Equations (24) and (25) gives
( D 2 k 1 2 ) ( D 2 k 2 2 ) e ¯ ˜ = A 0 q 0 β   Γ ( β 2 α 2 ) Q ¯ ( s ) η ˜ ( α ) e β z .
Solving (32), one can obtain
e ¯ ˜ ( α , z , s ) = i = 1 2 E i e k i z + e p ( α , s ) e β z ,
where e p is the coefficient of the particular solution of (32).
e p = A 0 q 0 Γ β ( β 2 α 2 ) Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) .
Introducing Equations (29) and (33) into Equation (25) yields
E i ( α , s ) = ( k i 2 α 2 k i 2 α 2 s 2 ) B i ( α , s ) ,
e p ( α , s ) = ( β 2 α 2 β 2 α 2 s 2 ) M ( α , s ) ,
So that
e ¯ ˜ = i = 1 2 ( k i 2 α 2 k i 2 α 2 s 2 ) B i e k i z A 0 q 0 Γ β ( β 2 α 2 ) Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z .
Applying (23) to Equation (16), one gets
( D 2 α 2 g   s 2 )   w ¯ ˜ = g D   θ ¯ ˜ ( g 1 ) D   e ¯ ˜ .
Substituting for θ ¯ ˜ and e ¯ ˜ in Equation (38), then solving the nonhomogeneous second-order ordinary differential equation, this gives
w ¯ ˜ = N ( α , s ) e ξ   z i = 1 2 ( k i k i 2 α 2 s 2 ) B i ( α , s ) e k i z   + A 0 q 0 Γ β 2 Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z ,
where ξ = α 2 + g   s 2 .
Applying (23) on the nondimensional form of Equation (2) and making use of the relation H n { r n 1 d d r ( r 1 + n f ¯ ( r , z , s ) ) } = α H n + 1 { f ¯ ( r , z , s ) } for n = 0 , it follows
e ¯ ˜ = α   u ¯ ˜ + d   w ¯ ˜ d   z ,
where H n { f ( r ) } represents the Hankel transform of f ( r ) .
Now, in order to obtain u ¯ ˜ , we shall substitute Equations (37) and (39) into Equation (40); one gets
u ¯ ˜ = 1 α ( ξ   N e ξ   z i = 1 2 ( α 2 k i 2 α 2 s 2 ) B i e k i z + A 0 q 0 Γ β α 2 Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z ) .
The stress components of the problem can be obtained by applying (23) on (18), (19), (20) and (21); this gives the following relations:
σ ¯ ˜ r r = 2 α   u ¯ ˜ + ( g 2 ) e ¯ ˜ g θ ¯ ˜ ,
σ ¯ ˜ z z = 2 d   w ¯ ˜ d z + ( g 2 ) e ¯ ˜ g θ ¯ ˜ ,
σ ¯ ˜ ϕ ϕ = 2 α   u ¯ ˜ + ( g 2 ) e ¯ ˜ g θ ¯ ˜ ,
σ ¯ ˜ r z = d   u ¯ ˜ d z α   w ¯ ˜ .
From the stresses above, only the stresses that related to the boundary conditions will be picked out to be calculated; hence, introducing Equations (31), (37), (39), and (41) into (43) and (45), it follows
σ ¯ ˜ z z = 2 ξ   N e ξ   z + i = 1 2 ( Ω k i 2 α 2 s 2 )   B i e k i z   A 0 q 0 Γ β   Ω   Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z ,
σ ¯ ˜ r z = 1 α ( Ω   N e ξ   z + i = 1 2 ( 2 α 2 k i k i 2 α 2 s 2 )   B i e k i z   2 A 0 q 0 Γ β 2 α 2   Q ¯ ( s ) η ˜ ( α ) ( β 2 k 1 2 ) ( β 2 k 2 2 ) e β z ) .
where Ω = 2 α 2 + g   s 2 .
Now making use of Equations (31), (46), and (47), Equations (22) can be written as
i = 1 2 k i   B i ( α , s ) =   M 1 + Γ   h   θ ¯ ˜ ( α , 0 , s ) ,
i = 1 2 F i   B i ( α , s ) 2   ξ   N ( α , s ) = M 2 ,
i = 1 2 F i + 2   B i ( α , s )   Ω   N ( α , s ) =   M 3 ,
where
M 1 = β   M ,   M 2 =   Ω   M ( β 2 α 2 s 2 ) ,   M 3 = 2 β   α 2 M ( β 2 α 2 s 2 ) , F i = ( Ω k i 2 α 2 s 2 ) ,   F i + 2 = ( 2 α 2 k i k i 2 α 2 s 2 ) ,   ( i = 1 , 2 ) .
Solving (48), (49), and (50), one gets
B 1 = ( h   Γ θ ¯ ˜ ( α ,   0 ,   s ) + M 1 ) ( 2 ξ   F 4   Ω   F 2 ) + ( 2   ξ M 3 Ω   M 2 ) k 2 ( Ω   F 2 2 ξ   F 4 ) k 1 + ( 2 ξ   F 3 Ω   F 1 ) k 2 ,
B 2 = ( h   Γ θ ¯ ˜ ( α ,   0 ,   s ) + M 1 ) ( Ω   F 1 2 ξ   F 3 ) + ( Ω   M 2 2   ξ M 3 ) k 1 ( Ω   F 2 2 ξ   F 4 ) k 1 + ( 2 ξ   F 3 Ω   F 1 ) k 2 ,
N = ( F 4 k 1 F 3 k 2 ) M 2 + ( h   Γ θ ¯ ˜ ( α ,   0 ,   s ) + M 1 ) ( F 1 F 4 F 2 F 3 ) + ( F 1 k 2 F 2 k 1 ) M 3 ( Ω   F 2 2 ξ   F 4 ) k 1 + ( 2 ξ   F 3 Ω   F 1 ) k 2 .
Substituting Equations (51)–(53) into Equation (31) and setting ( z = 0 ) , this gives the surface temperature θ ¯ ˜ ( α , 0 , s ) as follows:
θ ¯ ˜ ( α , 0 , s ) = M 1   ( ( 2 ξ   F 3 Ω   F 1 ) ( β + k 2 ) + ( Ω   F 2   2 ξ F 4 ) ( β + k 1 ) ) + β ( k 1 k 2 ) ( Ω   M 2 + 2   ξ M 3 ) β ( ( 2 F 4   ξ F 2   Ω ) ( h Γ + k 1 ) + ( F 1   Ω 2 F 3   ξ ) ( h Γ + k 2 ) ) .

4. Small Time Approximation

The above calculations and the roots of the auxiliary Equation (28) foretell that converting the solution to the physical domain is too difficult to be performed using the usual methods. So, a small time approximation will be adapted for the roots, conforming to great s-values; this technique was used in [36] according to [38,39].
The roots of Equation (28) are given as
k 1 = 2 α 2 + s 2 + s Γ   ( ε + 1 ) + s   s 2 + 2 s Γ   ( ε 1 ) + Γ 2 ( ε + 1 ) 2 2 ,
k 2 = 2 α 2 + s 2 + s Γ   ( ε + 1 ) s   s 2 + 2 s Γ   ( ε 1 ) + Γ 2 ( ε + 1 ) 2 2 .
Rewrite Equations (55) and (56) by setting p = s 1 in the form k i = p 1 f i ( p ) , where
f 1 = 2 α 2 p 2 + 1 + p Γ ( ε + 1 ) + 1 + 2 p Γ ( ε 1 ) + p 2 Γ 2 ( ε + 1 ) 2 2 , f 2 = 2 α 2 p 2 + 1 + p Γ ( ε + 1 ) 1 + 2 p Γ ( ε 1 ) + p 2 Γ 2 ( ε + 1 ) 2 2 , }
and Γ becomes Γ = ( p 2 +   p   τ 2 + τ 2 2 2 ) p   ( p + τ 1 ) .
Expanding f i ( p ) ,   ( i = 1 , 2 ) in a Taylors series and considering the first four terms only, one gets
f i ( p ) = d i 0 + d i 1 p + d i 2 p 2 + d i 3 p 3 ,
where the coefficients are listed in the Appendix A.
Again, expand [ f i ( p ) ] 1 2 in a Taylor series, considering only the first three terms; this gives the expressions for k 1 and k 2 as
k i = p 1 ( b i 0 + b i 1 p + b i 2 p 2 ) ,   ( i = 1 , 2 )
where
b i 0 = d i 0 , b i 1 = d i 1 2 d i 0 , b i 2 = ( 4 d i 2   d i 0 d i 1 2 ) 8 d i 0 3 / 2 .
Finally, k i can be written in the form
k i = b i 0 s + b i 1 .   ( i = 1 , 2 )

5. Inversion of the Transformation

Owing to the extreme intricacy of obtaining the studied fields in the physical domain, the inverse formula of the double transformation will be used separately. To do that, let f ¯ ˜ ( α , z , s ) be a function in the transformed domain first; the inverse Hankel transform will be performed by utilizing the following formula:
f ¯ ( r , z , s ) = 0 α   J n ( r   α )   f ¯ ˜ ( α , z , s ) d   α ,
where f ¯ ( r , z , s ) is a function in the Laplace domain.
Second, the function f ¯ ( r , z , s ) will be returned to the time domain by applying the Riemann–Sum approximation method using the formula
f ( t ) = e ζ   t t [ 1 2 f ¯ ( ζ ) + ( 1 ) n Re n = 1 N   ( f ¯ ( ζ   + i n π t ) ) ] ,
where N is an integer that indicates the number of terms and should be taken as adequately large. The numerical experiments show that as ζ 4.7 / t , the relation (62) convergences faster [40].

6. Application and Computation

As an example, the computations were implemented on a copper half-space subjected to a laser pulse, the intensity of which is q 0 = 10 4   W / m 2 , assuming that τ 2 = 75 × 10 5   s and β = 10 6   m 1 . To ensure the validity of the obtained solution, θ ( r , 0 , t ) ,   θ ,   w and σ z z were calculated for some values of h . The radial and temporal shapes considered in the computations are chosen to be Gaussian distributions, and thus
Q ( t ) = e ( t a b ) 2 ,
where a and b denote the times at which Q ( t ) becomes maximum and reduces to e 1 , respectively.
The mechanical constants of the considered medium are taken as (see [22])
T 0 = 293   K ρ = 8954   kg   m 3 α T = 1.78 × 10 5   K 1 c E = 383.1   J ( kg ) 1   K 1 ω = 10 3   m k = 386   W   m 1   K 1 λ = 7.76 × 10 10   kg   m 1   s 2   a = 3 × 10 3   s    b = 1 × 10 3   s    A 0 = 0.01     μ = 3.86 × 10 10   kg   m 1   s 2

7. Results and Discussion

Figure 1 represents the surface temperature θ ( r , 0 , t ) calculated for τ 1 = 0.6   τ 2 in the absence of cooling. The figure shows that the surface temperature increases until it reaches its maximum, which clearly shifted to a greater time than that of the maximum of Q ( t ) ( 3 × 10 3 ) . After θ ( r , 0 , t ) reaches its maximum, it begins to decrease and does not reach the zero-value even after the pulse is switched off. Moreover, in the radial direction, θ ( r , 0 , t ) has a weak gradient near the surface; this gradient becomes stronger as r-increases until the temperature vanishes. The behavior of the temperature depends on two main occurrences, namely, the conductivity of the material and the absorbed energy, where at the beginning of the illumination process and due to the increased pulse, the absorbed energy will be greater than the conductivity of the material, so θ ( r , 0 , t ) increases; this behavior lasts up to when the absorbed energy is matched to the conducted energy; exactly at this time, the surface temperature attains its maximum. After that, and according to the chosen pulse shape, the absorbed energy begins to decrease, and the conductivity of the material gains the upper hand; this leads to a decrease in the surface temperature. After the pulse is switched off, the gradient of the temperature becomes smaller than that during the pulse, which is due to the absence of the laser radiation; this leads to a small gradient and, consequently, small conduction of the heat energy.
Figure 2 represents the time-dependent surface temperature, with h as a parameter, calculated for τ 1 =   0.6   τ 2 and r = 5 × 10   4 . As expected, the surface temperature clearly decreases with increasing values of the cooling coefficient; moreover, as the cooling coefficient takes small values, the maximum value of the surface temperature is shifted towards greater time than that of the maximum of Q ( t ) . Beside the general decrease in surface temperature with increasing h , its behavior after it reaches its maximum, and at the times where the laser is switched off, decreases until the surface temperature approximately takes the shape of Q ( t ) . This behavior is due to the contribution of the cooling effect together with the conductivity of the material, especially after the laser is switched off and the absorbed energy is stopped.
Figure 3 represents the time-dependent laser pulse shape and the surface temperature, with τ 1 as a parameter in the absence of cooling, calculated for r = 5 × 10   4 . The figure shows that compared to the chosen temporal pulse shape, Q ( t ) , for τ 1 > τ 2 , the maximum surface temperature becomes smaller and shifts evidently towards greater time than that for τ 1 < τ 2 . This means that the surface temperature for τ 1 > τ 2 needs more time and consumes more energy to reach its maximum; this behavior is similar to the behavior that appears in [31] for the classical coupled theory.
Figure 4 represents the time-dependent laser pulse shape and the surface temperature, with τ 1 as a parameter in the presence of heat losses, calculated for h = 100 , r = 5 × 10   4 , and t = 4 × 10   3 . As seen, the effect of cooling is very pronounced, where the values of the surface temperature are decreased and its maximum is slightly shifted from the time of the maximum of Q ( t ) . Moreover, after the laser is switched off, the temperature approximately reaches zero, which does not occur in the absence of cooling; these behaviors are due to the cooling and the conductivity of the material. The behavior of this figure is slightly similar to the behavior of the previous figure, where the time needed to reach the peak is slightely shifted towards greater time for τ 1 > τ 2 than for τ 1 < τ 2 .
Figure 5 represents the temperature θ as a function of z and t in the absence of cooling, calculated for τ 1 = 0.6   τ 2 and r = 5 × 10   4 . The finite velocity that distinguishes the DPL model appears evidently through the strong gradient at all locations, especially before the pulse is switched off. After the temperature reaches its maximum, a small temperature gradient begins to appear near the surface. The figure shows that the temperature is fully compatible with the time pulse profile, where the temperature increases until it reaches its maximum and then begins to decline. The temperature shows its highest value at the irradiated surface; this value decreases as z increases until it vanishes.
Figure 6 represents the temperature of the medium as a function of z and r in the absence of cooling, calculated for τ 1 = 0.6   τ 2 and t = 4 × 10   3 . It is evident that the maximum temperature occurs at the illuminated surface and the temperature has a weak gradient in a small region adjacent to the target surface, as z and r increases; the gradient becomes strong until the temperature vanishes. This figure agrees with the previous figure, for which the finite velocity of the employed model evidently appears.
Figure 7 represents the temperature versus z, with h as a parameter, calculated for τ 1 =   0.6   τ 2 and r = 5 × 10   4 . As expected, the temperature decreases as h increases; moreover, a weak gradient near the surface is observed; this gradient moves deeper into the medium as h increases. The maximum temperature is slightly shifted, nearer to the surface; this behavior can be observed from the value h = 20 .
Figure 8 represents the temperature versus z in the absence of cooling, with τ 1 as a parameter, calculated for r = 5 × 10   4 and t = 4 × 10   3 . The figure shows that, near the surface, the temperature gradients in the case of τ 1 > τ 2 are stronger than that of τ 1 < τ 2 , which has lower gradients; this behavior is evident for the curve of τ 1 = 0.3   τ 2 . Additionally, the penetration of temperature into the medium increases as τ 1 increases and becomes greater as τ 1 > τ 2 ; this behavior is clearly shown for τ 1 = 2   τ 2 , demonstrated by comparison with the other values of τ 1 . Again, the behavior of the temperature when τ 1 > τ 2 is similar to the behavior of the classical coupled theory of thermoelasticity in [31].
Figure 9 represents the temperature distribution versus z, with τ 1 as a parameter, calculated for h = 100 , r = 5 × 10   4 and t = 4 × 10   3 . The main behavior of the previous figure is observed here for both the gradient of the temperature and the penetration into the medium. The figure shows a pronounced contribution to the cooling, where for τ 1 τ 2 it is clear that the maximum does not occur at the illuminated surface but is shifted to greater z values and the temperature at the surface decreases as τ 1 decreases.
Figure 10 represents the component of the displacement w of the medium as a function of z   and t in the absence of cooling, calculated for τ 1 = 0.6 τ 2 and r = 5 × 10   4 . The negative sign of the displacement refers to its direction, where it occurs in the direction of the free surface. As seen, the displacement does not show a behavior before the time t = 0.002 ; this behavior is attributed to the small temperature at this time; moreover, the displacement increases with increased time. The maximum displacement appears at the surface and decreases as z takes greater values until it vanishes; this behavior is due to temperature behavior.
Figure 11 represents the component of the displacement w of the medium as a function of z and   r in the absence of cooling, calculated for τ 1 = 0.6 τ 2 and t = 4 × 10   3 . This figure agrees with the previous figure, where the maximum displacement is observed at the surface and decreases by increasing both z and r . The figure also shows a negative sign for the displacement in the r-direction; this behavior is due to a decrease in the distance between the particles in the r-direction.
Figure 12 represents the distribution of the displacement w , with h as a parameter, calculated for τ 1 =   0.6   τ 2 and r = 5 × 10   4 . As seen from the figure, the displacement decreases with an increase in h ; this behavior is owed to the decrease in the corresponding temperature (see Figure 7).
Figure 13 represents the displacement w in the absence of cooling, with τ 1 as a parameter, calculated for r = 5 × 10   4 and t = 4 × 10   3 . The figure shows that as τ 1 takes values greater than τ 2 , the displacement decreases, and vice versa. Observing the behaviors of both the temperature and the displacement, one can note that the smallest and deeper temperatures produce the smallest displacements; this conduct can be explained from Figure 8, where for τ 1 > τ 2 , the gradient of the temperature in the vicinity of the surface and the penetration into the medium is greater than in the case of τ 1 < τ 2 ; this means that the first case consumes more energy in heating the surrounding and to penetrate than the second case and, thus, its temperature and displacement is the smallest.
Figure 14 represents the displacement w with the existence of cooling, with τ 1 as a parameter, calculated for r = 5 × 10   4 and t = 4 × 10   3 . The figure shows a clear influence for the heat losses, where the difference between the three curves becomes approximately slight; this is despite the compatibility of the behavior of this figure with the previous figure.
Figure 15 represents the stress component σ z z as a function of z and t in the absence of cooling, calculated for τ 1 = 0.6   τ 2 and r = 5 × 10   4 . The general behavior of σ z z depends on the behavior of both the gradient of the displacement and the temperature (see equation 7); hence, as seen, before the laser pulse reaches its peak, the temperature and the displacement increase, but the temperature has the superiority; this explains the negative values of σ z z . After the laser pulse reaches its peak, the absorbed energy decreases, and then the temperature decreases while the displacement and its gradient are still increasing; this explains the positive values of σ z z near the surface. One can observe that this behavior is clear after the temperature reaches its maximum. By increasing z-values, the influence of the gradient of the displacement begins to decay to permit the temperature to gain the upper hand again; this gives σ z z negative values until the stress vanishes.
Figure 16 represents the stress component σ z z as a function of z and r in the absence of cooling, calculated for τ 1 = 0.6   τ 2 and t = 4 × 10   3 . It is clear that the stress σ z z at the surface satisfies the boundary condition for all values of r. Near the surface, σ z z begins to increase until it reaches a maximum value, then it decreases with increasing z until the stress vanishes (this will be seen in a later figure). The figure also shows that except at the surface ( z = 0 ) , the stress begins with a maximum value in the r-direction, then decreases until it vanishes.
Figure 17 represents the stress component σ z z versus z, with h as a parameter, calculated for τ 1 =   0.6   τ 2 and r = 5 × 10   4 . As seen for different h values, the stresses in the vicinity of the surface are almost identical in their gradients. As z increases, the performance of the temperature with cooling appears evident, as in Figure 7.
Figure 18 represents the stress component σ z z versus z in the absence of cooling, with τ 1 as a parameter, calculated for r = 5 × 10   4 and t = 4 × 10   3 . The figure shows that the stress in the vicinity of the surface is slightly affected by τ 1 ; as z increases, a pronounced effect for τ 1 is observed, where the behavior of the temperature appears with τ 1 (see Figure 8).
Figure 19 represents the stress component σ z z versus z, with τ 1 as a parameter in the existence of cooling, calculated for r = 5 × 10   4 and t = 4 × 10   3 . As seen, the behavior of σ z z at the surface coincides for the three curves; as z increases, the effect of the temperature is clear (see Figure 9).

8. Conclusions

In this paper, the response generated by the volumetric absorption technique of laser pulse in a 2D generalized thermoelastic semi-infinite medium is discussed in the absence and presence of temperature-dependent heat losses by employing the DPL model. From the above discussion, it can be deduced that
  • The results obtained in the absence and presence of cooling are not in conflict with well-known physical conduct, especially for τ 1 < τ 2 .
  • The translation times show clear effects on the behavior of all the studied fields in the absence of cooling and on some studied fields in the presence of cooling.
  • For τ 1 > τ 2 , the thermoelastic response of the DPL model behaves similar to the thermoelastic response of the classical coupled theory of thermoelasticity.
  • The heat losses show a pronounced effect, especially on the temperature, in which its maximum shifted into the medium.

Funding

This research is supported by the deanship of scientific research in Majmaah University under project 1439-5.

Acknowledgments

The author is grateful for the support of the deanship of scientific research in Majmaah University under project 1439-5.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

d 10 = 1 2 ( 1 + ( ε + 1 ) 2 η 0 + η 1 ) ,   d 20 = 1 2 ( 1 + ( ε + 1 ) 2 η 0 η 1 ) d 11 = 1 2 ( ( ε + 1 ) η 2 + Φ 0 2 η 1 ) ,   d 21 = 1 2 ( ( ε + 1 ) η 2 Φ 0 2 η 1 ) d 12 = 1 2 ( 2 α 2 + ( ε + 1 ) η 3 + 1 2 η 1 ( 2 ( ε 1 ) η 3 + ( ε + 1 ) 2 η 4 η 1 2 Φ 1 ) ) d 22 = 1 2 ( 2 α 2 + ( ε + 1 ) η 3 1 2 η 1 ( 2 ( ε 1 ) η 3 + ( ε + 1 ) 2 η 4 η 1 2 Φ 1 ) ) d 13 = 1 2 ( ( ε + 1 ) η 3 τ 1 + ( 2 ( ε 1 ) η 3 τ 1 + ( ε + 1 ) 2 η 5 2 η 1 Φ 0 Φ 2 4 η 1 ) ) d 23 = 1 2 ( ( ε + 1 ) η 3 τ 1 ( 2 ( ε 1 ) η 3 τ 1 + ( ε + 1 ) 2 η 5 2 η 1 Φ 0 Φ 2 4 η 1 ) )
where
η 0 = τ 2 2 τ 1 ,   η 1 = 1 + ( ε 1 ) η 0 + ( ε + 1 ) 2 4 η 0 2 ,     η 2 = ( τ 2 τ 1 + τ 2 2 2   τ 1 2 ) , η 3 = 2 τ 1 2 2 τ 1 τ 2 + τ 2 2 2 τ 1 3 ,   η 4 = 8 τ 1 2 τ 2 2 8 τ 1 τ 2 3 + 4 τ 2 4 4 τ 1 4 η 5 = 2 τ 1 3 τ 2 4 τ 1 2 τ 2 2 + 3 τ 1 τ 2 3 τ 2 4 4 τ 1 5 ,   Φ 0 = 2 ( ε 1 ) η 2 + ( ε + 1 ) 2 η 0 η 2 Φ 1 = Φ 0 2 4   η 1 4 ,   Φ 2 = 2 ( ε 1 ) η 3 + ( ε + 1 ) 2 η 4 η 1 2 Φ 1

References

  1. Biot, M.A. Thermoelasticity and Irreversible Thermodynamics. J. Appl. Phys. 1956, 27, 240–253. [Google Scholar] [CrossRef]
  2. Al-Huniti, N.S.; Al-Nimr, M. Thermoelastic response of a heated thin composite plate using the hyperbolic heat conduction model: Lumped analysis. Int. J. Therm. Sci. 2004, 43, 959–965. [Google Scholar] [CrossRef]
  3. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  4. Lord, H.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  5. Cattaneo, C. Sulla Conduzione Del Calore. Some Asp. Diffus. Theory 1948, 3, 83–101. [Google Scholar]
  6. Cattaneo, C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation. C. R. Chim. 1958, 247, 431–433. [Google Scholar]
  7. Green, A.E.; Lindsay, K.A. Thermoelasticity. J. Elast. 1972, 2, 1–7. [Google Scholar] [CrossRef]
  8. Green, A.E.; Naghdi, P.M. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  9. Hetnarski, R.B.; Ignaczak, J. Soliton-like waves in a low temperature nonlinear thermoelastic solid. Int. J. Eng. Sci. 1996, 34, 1767–1787. [Google Scholar] [CrossRef]
  10. Tzou, D.Y. A Unified Field Approach for Heat Conduction from Macro- to Micro-Scales. J. Heat Transf. 1995, 117, 8–16. [Google Scholar] [CrossRef]
  11. Tzou, D.Y. The generalized lagging response in small-scale and high-rate heating. Int. J. Heat Mass Transf. 1995, 38, 3231–3240. [Google Scholar] [CrossRef]
  12. Zenkour, A.M.; Mashat, D.; Abouelregal, A. The effect of dual-phase-lag model on reflection of thermoelastic waves in a solid half space with variable material properties. Acta Mech. Solida Sin. 2013, 26, 659–670. [Google Scholar] [CrossRef]
  13. Othman, M.I.A.; Abd-Elaziz, E.M. The Effect of Thermal Loading Due to Laser Pulse in Generalized Thermoelastic Medium with Voids in Dual Phase Lag Model. J. Therm. Stress. 2015, 38, 1068–1082. [Google Scholar] [CrossRef]
  14. Lotfy, K. The elastic wave motions for a photothermal medium of a dual-phase-lag model with an internal heat source and gravitational field. Can. J. Phys. 2016, 94, 400–409. [Google Scholar] [CrossRef] [Green Version]
  15. Abbas, I.A.; Mohamed, E.A.E. Dual-Phase-Lag Model on Generalized Magneto-Thermoelastic Interaction in a Functionally Graded Material. Int. J. Acoust. Vib. 2017, 22, 369–376. [Google Scholar] [CrossRef]
  16. Abouelregal, A.; Abo-Dahab, S.M. A two-dimensional problem of a mode-I crack in a rotating fibre-reinforced isotropic thermoelastic medium under dual-phase-lag model. Sadhana 2018, 43, 13. [Google Scholar] [CrossRef] [Green Version]
  17. Kalkal, K.K.; Sheokand, S.K.; Deswal, S. Reflection and transmission between thermoelastic and initially stressed fiber-reinforced thermoelastic half-spaces under dual-phase-lag model. Acta Mech. 2019, 230, 87–104. [Google Scholar] [CrossRef]
  18. Othman, M.I.A.; Abd-Elaziz, E.M.; Hilal, M.I.M. State-space approach to a 2-D generalized thermoelastic medium under the effect of inclined load and gravity using a dual-phase-lag model. Mech. Based Des. Struct. Mach. 2020, 1–17. [Google Scholar] [CrossRef]
  19. Allam, M.N.; Tayel, I.M. Thermal effect on transverse vibrations of a nonhomogeneous rectangular thin plate subjected to a known temperature distribution. Trans. Can. Soc. Mech. Eng. 2020, 44, 452–460. [Google Scholar] [CrossRef]
  20. Hassan, A.; El-Nicklawy, M.; El-Adawi, M.; Nasr, E.; Hemida, A.; El-Ghaffar, O. Heating effects induced by a pulsed laser in a semi-infinite target in view of the theory of linear systems. Opt. Laser Technol. 1996, 28, 337–343. [Google Scholar] [CrossRef]
  21. Ezzat, M.A.; El Karamany, A.S.; Fayik, M.A. Fractional Ultrafast Laser–Induced Thermo-Elastic Behavior in Metal Films. J. Therm. Stress. 2012, 35, 637–651. [Google Scholar] [CrossRef]
  22. Henain, E.F.; Hassan, A.F.; Megahed, F.; Tayel, I.M. Thermo-Elastic Half Space Under Illumination of a Laser Beam by Using Lord and Shulman Theory. J. Therm. Stress. 2014, 37, 51–72. [Google Scholar] [CrossRef]
  23. Youssef, H.M.; El-Bary, A.A. Thermoelastic Material Response Due to Laser Pulse Heating in Context of Four Theorems of Thermoelasticity. J. Therm. Stress. 2014, 37, 1379–1389. [Google Scholar] [CrossRef]
  24. Elhagary, M.A. A Two-Dimensional Generalized Thermoelastic Diffusion Problem for a Thick Plate Subjected to Thermal Loading Due to Laser Pulse. J. Therm. Stress. 2014, 37, 1416–1432. [Google Scholar] [CrossRef]
  25. Zenkour, A.M.; Aboelregal, A.E. The fractional effects of a two-temperature generalized thermoelastic semi-infinite solid induced by pulsed laser heating. Arch. Mech. 2015, 67, 53–73. [Google Scholar]
  26. Bassiouny, E.; Abouelnaga, Z.; Youssef, H.M. One-dimensional thermoelastic problem of a laser pulse under fractional order equation of motion. Can. J. Phys. 2017, 95, 464–471. [Google Scholar] [CrossRef]
  27. Allam, M.N.M.; Tayel, I.M. Generalized thermoelastic functionally graded half space under surface absorption of a laser radiation. J. Theor. Appl. Mech. 2017, 55, 155–165. [Google Scholar] [CrossRef]
  28. Abbas, I.A.; Marin, M. Analytical solution of thermoelastic interaction in a half-space by pulsed laser heating. Physica E Low Dimens. Syst. Nanostruct. 2017, 87, 254–260. [Google Scholar] [CrossRef]
  29. Abouelregal, A.E.; Zenkour, A.M. A generalized thermoelastic medium subjected to pulsed laser heating via a two-temperature model. J. Theor. Appl. Mech. 2019, 57, 631–639. [Google Scholar] [CrossRef]
  30. Tayel, I.M.; Hassan, A.F. Heating a thermoelastic half space with a surface absorption pulsed laser using fractional order theory of thermoelasticity. J. Theor. Appl. Mech. 2019, 57, 489–500. [Google Scholar] [CrossRef]
  31. Tayel, I.M. Thermoelastic Response Induced by Volumetric Absorption of Uniform Laser Radiation in a Half-Space. Coatings 2020, 10, 228. [Google Scholar] [CrossRef] [Green Version]
  32. Hetnarski, R. Coupled One-Dimensional Thermal Shock Problem for Small Times. Arch. Mech. 1961, 13, 295–306. [Google Scholar]
  33. Hetnarski, R. The Fundamental Solution of the Coupled Thermoelastic Problem for Small Times. Arch. Mech. 1964, 1, 16. [Google Scholar]
  34. Sherief, H.H.; Darwish, A.A. A short time solution for a problem in thermoelasticity of an infinite medium with a spherical cavity. J. Therm. Stress. 1998, 21, 811–828. [Google Scholar] [CrossRef]
  35. Sherief, H.H.; Elmisiery, A.E.M.; Elhagary, M.A. Generalized thermoelastic problem for an infinitely long hollow cylinder for short times. J. Therm. Stress. 2004, 27, 885–902. [Google Scholar] [CrossRef]
  36. Sherief, H.H.; Hamza, F.A.; El-Latief, A.M.A. Wave propagation study for axi-symmetric 2D problems of a generalized thermo-visco-elastic half space. J. Therm. Stress. 2019, 42, 835–848. [Google Scholar] [CrossRef]
  37. Debnath, L.; Bhatta, D. Integral Transforms and Their Applications; Taylor and Francis: New York, NY, USA, 2015. [Google Scholar]
  38. Boley, B.A. Discontinuities in integral-transform solutions. Q. Appl. Math. 1962, 19, 273–284. [Google Scholar] [CrossRef] [Green Version]
  39. Boley, B.A.; Hetnarski, R.B. Propagation of Discontinuities in Coupled Thermoelastic Problems. J. Appl. Mech. 1968, 35, 489–494. [Google Scholar] [CrossRef]
  40. Tzou, D.Y. Experimental support for the lagging behavior in heat propagation. J. Thermophys. Heat Transf. 1995, 9, 686–693. [Google Scholar] [CrossRef]
Figure 1. The surface temperature θ ( r , 0 , t ) in the absence of cooling.
Figure 1. The surface temperature θ ( r , 0 , t ) in the absence of cooling.
Mathematics 08 02094 g001
Figure 2. The time-dependent surface temperature vs. t with h as a parameter.
Figure 2. The time-dependent surface temperature vs. t with h as a parameter.
Mathematics 08 02094 g002
Figure 3. The time-dependent surface temperature vs. t with τ 1 as a parameter for h = 0 .
Figure 3. The time-dependent surface temperature vs. t with τ 1 as a parameter for h = 0 .
Mathematics 08 02094 g003
Figure 4. The time-dependent surface temperature vs. t with τ 1 as a parameter for h = 0 .
Figure 4. The time-dependent surface temperature vs. t with τ 1 as a parameter for h = 0 .
Mathematics 08 02094 g004
Figure 5. The temperature θ as a function of z and t in the absence of cooling.
Figure 5. The temperature θ as a function of z and t in the absence of cooling.
Mathematics 08 02094 g005
Figure 6. The temperature θ as a function of z and r in the absence of cooling.
Figure 6. The temperature θ as a function of z and r in the absence of cooling.
Mathematics 08 02094 g006
Figure 7. The temperature θ versus z , with h as a parameter for τ 1 = 0.6   τ 2 .
Figure 7. The temperature θ versus z , with h as a parameter for τ 1 = 0.6   τ 2 .
Mathematics 08 02094 g007
Figure 8. The temperature θ versus z , with τ 1 as a parameter in the absence of cooling.
Figure 8. The temperature θ versus z , with τ 1 as a parameter in the absence of cooling.
Mathematics 08 02094 g008
Figure 9. The temperature θ versus z , with τ 1 as a parameter for h 0 .
Figure 9. The temperature θ versus z , with τ 1 as a parameter for h 0 .
Mathematics 08 02094 g009
Figure 10. The displacement w as a function of z and t in the absence of cooling.
Figure 10. The displacement w as a function of z and t in the absence of cooling.
Mathematics 08 02094 g010
Figure 11. The displacement w as a function of z and r in the absence of cooling.
Figure 11. The displacement w as a function of z and r in the absence of cooling.
Mathematics 08 02094 g011
Figure 12. The displacement w versus z , with h as a parameter for τ 1 = 0.6 τ 2 .
Figure 12. The displacement w versus z , with h as a parameter for τ 1 = 0.6 τ 2 .
Mathematics 08 02094 g012
Figure 13. The displacement w versus z , with τ 1 as a parameter in the absence of cooling.
Figure 13. The displacement w versus z , with τ 1 as a parameter in the absence of cooling.
Mathematics 08 02094 g013
Figure 14. The displacement w versus z , with τ 1 as a parameter in the presence of cooling.
Figure 14. The displacement w versus z , with τ 1 as a parameter in the presence of cooling.
Mathematics 08 02094 g014
Figure 15. The stress σ z z as a function of z and t in the absence of cooling.
Figure 15. The stress σ z z as a function of z and t in the absence of cooling.
Mathematics 08 02094 g015
Figure 16. The stress σ z z as a function of z and r in the absence of cooling.
Figure 16. The stress σ z z as a function of z and r in the absence of cooling.
Mathematics 08 02094 g016
Figure 17. The stress σ z z versus z , with h as a parameter for τ 1 = 0.6 τ 2 .
Figure 17. The stress σ z z versus z , with h as a parameter for τ 1 = 0.6 τ 2 .
Mathematics 08 02094 g017
Figure 18. The stress σ z z versus z with τ 1 as a parameter in the absence of cooling.
Figure 18. The stress σ z z versus z with τ 1 as a parameter in the absence of cooling.
Mathematics 08 02094 g018
Figure 19. The stress σ z z versus z , with τ 1 as a parameter in the presence of cooling.
Figure 19. The stress σ z z versus z , with τ 1 as a parameter in the presence of cooling.
Mathematics 08 02094 g019
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tayel, I.M. Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation. Mathematics 2020, 8, 2094. https://doi.org/10.3390/math8112094

AMA Style

Tayel IM. Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation. Mathematics. 2020; 8(11):2094. https://doi.org/10.3390/math8112094

Chicago/Turabian Style

Tayel, Ismail M. 2020. "Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation" Mathematics 8, no. 11: 2094. https://doi.org/10.3390/math8112094

APA Style

Tayel, I. M. (2020). Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation. Mathematics, 8(11), 2094. https://doi.org/10.3390/math8112094

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