Elastodynamic Response of Thermal Laser Pulse in Micropolar Thermoelastic Mass Diffusion Medium

The present investigation deals with the deformation in micropolar thermoelastic diffusion medium due to inclined load subjected to thermal laser pulse. Normal mode analysis technique is used to solve the problem. The inclined load is assumed to be a linear combination of a normal load and a tangential load. The closed form expressions of normal stress, tangential stress, couple stress, temperature distribution, and mass concentration are obtained. A computer program has been developed to derive the physical quantities numerically.The variation of normal stress, tangential stress, coupled stress, temperature change, andmass concentration is depicted graphically to show the effect of relaxation times and mass concentration. Some particular cases of interest are deduced from the present investigation.


Introduction
Modern engineering structures are often made up of materials possessing the internal structure.Polycrystalline materials and materials with fibrous or coarse grain structure come in this category.The classical theory of elasticity is not sufficient to explain the phenomenon of high frequency short wave length and ultrasonic waves.So micropolar theory was developed to overcome the shortcomings of the classical theory of elasticity by considering the granular structure of the material of the medium.The micropolar theory of elasticity is applied to materials for problems where the classical theory of elasticity fails owing to the microstructure of the material.The linear theory of micropolar elasticity was developed by Eringen [1].Under this theory, solids can undergo macrodeformations and microrotations.Also they can support couple stresses in addition to force stresses.Nowacki [2][3][4] extended the micropolar theory of elasticity to include the thermal effects.
Diffusion is defined as the spontaneous movement of the particles from a high-concentration region to the lowconcentration region, and it occurs in response to a concentration gradient expressed as the change in the concentration due to change in position.Thermal diffusion utilizes the transfer of heat across a thin liquid or gas to accomplish isotope separation.Today, thermal diffusion remains a practical process to separate isotopes of noble gases, for example, Xenon, and other light isotopes, for example, Carbon, for research purposes.In most of the applications, the concentration is calculated using Fick's law.This is a simple law which does not take into consideration the mutual interaction between the introduced substance and the medium into which it is introduced or the effect of temperature of this interaction.However, there is a certain degree of coupling with temperature and temperature gradients as temperature speeds up the diffusion process.The thermodiffusion in elastic solids is due to coupling of fields of temperature and mass diffusion process.The thermodiffusion in elastic solids is due to coupling of fields of temperature, mass diffusion, and of strain in addition to heat and mass exchange with the environment.Nowacki [2][3][4][5] developed the theory of thermoelastic diffusion by using coupled thermoelastic model.Dudziak and Kowalski [6], respectively, discussed the theory of thermodiffusion and coupled quasi-stationary problems of thermal diffusion for an elastic layer.They studied the influence of cross-effects arising from the coupling 2 Journal of Thermodynamics of the fields of temperature, mass diffusion, and strain due to which the thermal excitation results in additional mass concentration and generates additional fields of temperature.
Laser technology has a vital application in nondestructive materials testing and evaluation.When a solid is heated with a laser pulse, it absorbs some energy which results in an increase in localized temperature.This causes thermal expansion and generation of the ultrasonic waves in the material.The irradiation of the surface of a solid by pulsed laser light generates wave motion in the solid material.There are generally two mechanisms for such wave generation, depending on the energy density deposited by the laser pulse.At high energy density, a thin surface layer of the solid material melts, followed by an ablation process whereby particles fly off the surface, thus giving rise to forces that generate ultrasonic waves.At low energy density, the surface material does not melt, but it expands at a high rate and wave and wave motion is generated due to thermoelastic processes.
Very rapid thermal processes (e.g., the thermal shock due to exposure to an ultrashort laser pulse) are interesting from the stand point of thermoelasticity, since they require a coupled analysis of the temperature and deformation fields.A thermal shock induces very rapid movement in the structural elements, giving rise to very significant inertial forces and thereby an increase in vibration.Rapidly oscillating contraction and expansion generate temperature changes in materials susceptible to diffusion of heat by conduction [7].This mechanism has attracted considerable attention due to the extensive use of pulsed laser technologies in material processing and nondestructive testing and characterization [8,9].The so-called ultrashort lasers are those with pulse durations ranging from nanoseconds to femtoseconds.In the case of ultrashort-pulse laser heating, the high intensity energy flux and ultrashort duration lead to very large thermal gradients or ultrahigh heating may exist at the boundaries.In such cases, as pointed out by many investigators, the classical Fourier model, which leads to an infinite propagation speed of the thermal energy, is no longer valid [10].Researchers have proposed several models to describe the mechanism of heat conduction during short-pulse laser heating, such as the parabolic one-step model [11], the hyperbolic one-step model [12], and the parabolic two-step and hyperbolic two-step models [13,14].It has been found that usually the microscopic two-step models, that is, parabolic and hyperbolic two-step models, are useful for thin films.
Scruby and Drain [15] considered the point source model to study the ultrasonic generation by lasers.They studied the heated surface by laser pulse irradiation in the thermoelastic system as a surface center of expansion (SCOE).They also discussed the applications of laser technology in flaw detection and acoustic microscopy.Rose [16] later presented a more exact mathematical basis.Point source model explains main features of laser-generated ultrasound waves but this model fails to explain precursor in epicenter waves.Later introducing the thermal diffusion McDonald [17] and Spicer et al. [18] proposed a new model known as laser-generated ultrasound model.This model reported excellent agreement between theory and experiment for metal materials.But due to the optical penetration effect, this model cannot be applied to the study of laser-generated ultrasound in nonmetallic material directly.The optical absorption occurs at the surface layer in metallic materials, and the heat penetration is resulting due to heat diffusion.In nonmetallic materials, the laser beam can penetrate the specimen to some finite depth and induced a buried bulk-thermal source, so the features of the laser-generated ultrasound will be significantly different from that in metallic materials.
Dubois et al. [19] experimentally demonstrated that penetration depth plays a very important role in the laserultrasound generation process.Ezzat et al. [20] discussed the thermoelastic behavior in metal films by fractional ultrafast laser.Al-Huniti and Al-Nimr [21] investigated the thermoelastic behavior of a composite slab under a rapid dual-phase lag heating.The comparison of one-dimensional and two-dimensional axisymmetric approaches to the thermomechanical response caused by ultrashort laser heating was studied by Chen et al. [22].Kim et al. [23] studied thermoelastic stresses in a bonded layer due to pulsed laser radiation.Thermoelastic material response due to laser pulse heating in context of four theorems of thermoelasticity was discussed by Youssef and El-Bary [24].Theoretical study of the effect of enamel parameters on laser induced surface acoustic waves in human incisor was studied by Yuan et al. [25].A two-dimensional generalized thermoelastic diffusion problem for a thick plate under the effect of laser pulse thermal heating was studied by Elhagary [26].
In this research, taking into account the mass concentration effect and radiation of ultrashort laser, we have established a model for micropolar thermoelastic medium with mass diffusion.The disturbance due to inclined loads has been studied in the proposed problem.The normal stress, tangential stress, coupled tangential stress, temperature distribution, and mass concentration are obtained numerically.

Basic Equations
Following Eringen [27], Sherief et al. [28], Kumar and Kansal [29], and Al-Qahtani and Datta [30], the basic equations for homogeneous, isotropic micropolar generalized thermoelastic solid with mass diffusion in the absence of body forces and body couples are given by The plate surface is illuminated by laser pulse given by the heat input: where  0 is the energy absorbed.The temporal profile () is represented as follows: Here  0 is the pulse rise time.The pulse is also assumed to have a Gaussian spatial profile in  1 : where  is the beam radius and as a function of the depth,  3 , the heat deposition due to the laser pulse is assumed to decay exponentially within the solid: Equation ( 7) with the aid of ( 8), (9), and (10) takes the following form: Here , , , , , and  are material constants,  is mass density, u = ( 1 ,  2 ,  3 ) is the displacement vector,  = ( 1 ,  2 ,  3 ) is the microrotation vector,  is temperature,  0 is the reference temperature of the body chosen,  is the concentration of the diffusion material in the elastic body,  * is the coefficient of the thermal conductivity,  * is the specific heat at constant strain,  is the thermoelastic diffusion constant,  is the coefficient describing the measure of thermodiffusion,  is the coefficient describing the measure of mass diffusion effects,  is the microinertia,  1 = (3+2+ ) 1 ,  2 = (3 + 2 + ) 1 ,  1 is the coefficient of linear thermal expansion,  1 is the coefficient of linear diffusion expansion,   are components of stress,   are components of couple stress,   are components of strain,   is the dilatation,   is Kronecker delta function,  0 and  1 are the diffusion relaxation times, and  0 and  1 are thermal relaxation times with  0 ≥  1 ≥ 0.Here  0 =  1 =  0 =  1 =  1 = 0 for coupled thermoelastic (CT) theory model.Consider  1 =  1 = 0,  = 1,  1 =  0 for Lord-Shulman (LS) model and  = 0 and  1 =  0 where  0 > 0 for Green-Lindsay (GL) model.In the above equations symbol "," followed by a suffix denotes differentiation with respect to spatial coordinates and a superposed dot ("⋅") denotes the derivative with respect to time, respectively.

Formulation of the Problem
We consider a micropolar generalized thermoelastic solid with mass diffusion medium with rectangular Cartesian coordinate system  1  2  3 with  3 -axis pointing vertically downward the medium.
Suppose that an inclined line load  0 per unit length is acting on -axis and its inclination to -direction is  (Figure 1).
For two-dimensional problems, we take the displacement vector and microrotation vector as For further consideration it is convenient to introduce in (1)-( 5) the dimensionless quantities defined by The displacement components  1 and  3 are related to the nondimensional potential functions  and  as follows: Substituting the values of  1 and  3 from ( 14) in ( 1)-( 4) with the aid of ( 12), we obtain Here,  1 = ( + )/

Solution of the Problem
The solution of the considered physical variables can be decomposed in terms of the normal modes as in the following form: Here  is the angular velocity and  is wave number.Making use of ( 21) and ( 15)-( 19) after some simplifications, we obtain where D = / The solution of the above system of ( 22)-( 25) satisfying the radiation conditions that (, , ,  2 , ) → 0 as  3 → ∞ is given as follows: where  2  ( = 1, 2, 3) are the roots of ( 22) and  2  ( = 4, 5) are the roots of characteristic equation of (25) and Substituting the values of , , ,  2 , and  from ( 26) in ( 5)-( 6) and using ( 12)-( 13) and ( 21) and solving the resulting equations, we obtain Consider   =     ,  = 1, 2, . . ., 5. Here,

Boundary Conditions
We consider normal and tangential force acting at the surface  3 = 0 along with vanishing of couple stress in addition to thermal and mass concentration boundaries considered at  3 = 0 and for  0 = 0. Mathematically, this can be written as where  1 and  2 are the magnitude of the applied force.
Substituting the expression of the variables considered into these boundary conditions, we can obtain the following equations: Equations ( 41) are solved by using the matrix method as follows: 5.1.Applications.Inclined Line Load.For an inclined line load  0 we have Making use of (42) in (41), we obtain the corresponding expressions of normal stress, tangential stress, couple stress, temperature distribution, and mass concentration due to inclined load.

Particular Cases
(i) If we take  1 =  1 = 0 and  = 1, in (41) with the aid of (43), we obtain the corresponding expressions of stresses, displacements, and temperature distribution for micropolar mass diffusion thermoelastic halfspace with one relaxation time.
(ii) If we take  = 0 in (41) with the aid of (43), the corresponding expressions of stresses, displacements, and temperature distribution are obtained for micropolar mass diffusion thermoelastic half-space with two relaxation times.
(iv) This is the case where  = 0 corresponds to normal load.
(v) This is the case where  = /2 corresponds to tangential load.

Special Case
Micropolar Thermoelastic Medium.If we neglect the diffusion effect in (41) with the aid of (43), we obtain the corresponding expressions of stresses, displacements, and temperature for micropolar generalized thermoelastic half-space.

Numerical Results and Discussions
The analysis is conducted for a magnesium crystal-like material.For numerical computations, following Eringen [31], the values of physical constants are as follows: Following Dhaliwal and Singh [32], thermal and diffusion parameters are given by A comparison of the dimensionless form of the field variables for the cases of micropolar thermoelastic mass diffusion medium with a laser pulse (MPL) and micropolar thermoelastic mass diffusion medium without a laser pulse (MP), for two different values of inclined angle  ( = 0 and  = /2), is presented in Figures 2-6.The values of all physical quantities for all cases are shown in the range 0 ≤  3 ≤ 2.
Solid lines and dash lines correspond to micropolar thermoelastic mass diffusion medium with a laser pulse  (MPL) and micropolar thermoelastic mass diffusion medium without a laser pulse (MP), respectively, for  = 0.
Solid lines with central symbol and dash lines with central symbol correspond to micropolar thermoelastic mass diffusion medium with a laser pulse (MPL) and micropolar thermoelastic mass diffusion medium without a laser pulse (MP), respectively, for  = /2.The computations were carried out in the absence and presence of laser pulse ( 0 = 10 5 , 0) and, on the surface of plane  1 = 1,  = 0.1.source.Initially the value of normal stress  33 decreases as  3 increases.The angle of inclination effect on the normal stress can be seen in this figure.For  = 0 the normal stress initially decreases and then increases and thereafter approaches the boundary surface away from the application of normal force.

Insulated Impermeable Boundary.
Figure 3 shows the variation of tangential stress  31 with the distance  3 .We notice that initially  31 shows similar behavior for  = 0 and  = /2 for MP.The initial behaviors of  31 in case of MPL for  = 0 and  = /2 is opposite to each other due to the difference between the angles of inclination of load applied.Away from the load applied the value of  31 approaches boundary surface.
Figure 4 displays the variation of coupled tangential stress  32 with the distance  3 .It is noticed that initially the behavior of  32 for all values of angle of inclination is similar.Also it is observed that the variation in coupled tangential stress for  = 0 is comparable to the variation for  = /2.The values of  32 for both values of  remain oscillatory and approach the boundary surface away from the applied force.
Figure 5 displays the variation of temperature  with distance  3 .The values of temperature  attain the maximum value near the application of source and away from it decreases and approaches the boundary surface for all the angles of inclination.In case of micropolar thermoelastic medium without laser pulse for  = 0 the temperature initially shows increasing trend but later on this trend is similar to other curves.
Figure 6 displays the variation of mass concentration  with distance  3 .The value of  decreases near the application of the source and then approaches the boundary surface far away from the load applied.

Conclusion
The problem of laser irradiation on micropolar thermoelastic mass diffusion medium is a significant problem in continuum mechanics.It is observed that the physical quantities are also affected by the different nonclassical theories of thermoelasticity with mass diffusion.It is observed from Figures 2-6 that angle of inclination effects on stress components, temperature change, and mass concentration depend upon the distance  3 .
The present problem has a significant application in geophysics and electronics engineering.The effect of diffusion is used to improve the conditions of oil extractions (seeking ways of more efficiently recovering oil from oil deposits).Also the study of thermal and diffusion effects plays an important role in understanding many seismological processes.Nowadays, there is a great deal of interest in the study of this phenomenon due to its application in geophysics and electronic industry.
The results obtained here are useful in engineering problems particularly in the determination state of stresses in a micropolar thermoelastic mass diffusion medium.Also, any particular case of special interest can be derived by assigning suitable values to the parameters and functions in the problem.

Figure 3 :
Figure 3: Variation of tangential stress with respect to distance  3 .

Figure 4 :
Figure 4: Variation of coupled tangential stress with respect to  3 .

Figure 5 :
Figure 5: Variation of temperature distribution with respect to  3 .