# Rotation and Magnetic Field Effect on Surface Waves Propagation in an Elastic Layer Lying over a Generalized Thermoelastic Diffusive Half-Space with Imperfect Boundary

Authors: S. M. Abo-Dahab; Kh. Lotfy; A. Gohaly
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The aim of the present investigation is to study the effects of magnetic field, relaxation times, and rotation on the propagation of surface waves with imperfect boundary. The propagation between an isotropic elastic layer of finite thickness and a homogenous isotropic thermodiffusive elastic half-space with rotation in the context of Green-Lindsay (GL) model is studied. The secular equation for surface waves in compact form is derived after developing the mathematical model. The phase velocity and attenuation coefficient are obtained for stiffness, and then deduced for normal stiffness, tangential stiffness and welded contact. The amplitudes of displacements, temperature, and concentration are computed analytically at the free plane boundary. Some special cases are illustrated and compared with previous results obtained by other authors. The effects of rotation, magnetic field, and relaxation times on the speed, attenuation coefficient, and the amplitudes of displacements, temperature, and concentration are displayed graphically.


---

## Body

## 1. Introduction

The foundations of magnetoelasticity were presented by Knopoff [1] and Chadwick [2] and developed by Kaliski and Petykiewicz [3]. An increasing attention is devoted to the interaction between magnetic field and strain field in a thermoelastic solid due to its many applications in the fields of geophysics, plasma physics, and related topics. All papers quoted above assumed that the interactions between the two fields take place by means of the Lorentz forces appearing in the equations of motion and by means of a term entering Ohm’s law and describing the electric field produced by the velocity of a material particle, moving in a magnetic field. The most ideal interface model, as it is known, is called perfect bond interface where the displacement and traction are continuous across the interface. However, interfaces are seldom perfect. Therefore, various imperfect models such as three-phase and linear models like spring models have been introduced by Yu et al. [4], Yu [5], and Benveniste [6]. Perhaps the most frequently studied imperfect interface model is the smooth bond interface, where the normal components of the displacements and traction are continuous across the interface, while the shear traction components are zero on the interface. Lord and Shulman [7] formulated a generalized theory of thermoelasticity with one thermal relaxation time, who obtained a wave equation by postulating a new law of heat conduction instead of classical Fourier’s law. Green and Lindsay [8] developed a temperature rate-dependent thermoelasticity that includes two thermal relaxation times and does not violate the classical Fourier’s law of heat conduction, when the body under consideration has a center of symmetry. Hetnarski and Ignaczak [9] introduced a review and presentation of generalized theories of thermoelasticity. Diffusion can be defined as the random walk of an assemble of particles from regions of high concentration to that of low concentration. Nowadays, there is a great deal of interest in the study of phenomena due to its applications in geophysics and electronic industry. In integrated circuit fabrication, diffusion is used to introduce “depants” in controlled amounts into semiconductor substance. In particular, diffusion is used to form the base and emitter in bipolar transistors, integrated resistors, and the source/drain in metal oxide semiconductor (MOS) transistors and polysilicon gates in MOS transistors. In most of the applications, the concentration is calculated using Fick’s law. This is simple law which does not take into consideration the mutual interaction between the introduced substance and the medium into which introduced. Study of the diffusion phenomenon is used to improve the conditions of oil extractions. These days’ oil companies are interested in the process of thermoelastic diffusion for more efficient extraction of oil from oil deposits. Until recently, thermodiffusion in solids, especially in metals, was considered as a quantity that is independent of body deformation. Practice, however, indicates that the process of thermodiffusion could have a very considerable influence on the deformation of the body. Thermodiffusion in elastic solid is due to the coupling of temperature, mass diffusion, and strain in addition to the exchange of heat and mass with the environment. Nowacki [10–13] developed the theory of thermoelastic diffusion by using coupled thermoelastic model. This implies infinite speed of propagation of thermoelastic waves. Olesiak and Pyryev [14] investigated the theory of thermoelastic diffusion and coupled quasistationary problems of thermal diffusion for an elastic layer. They studied the influence of cross effects arising from the coupling 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. Sherief et al. [15] developed the generalized theory of thermoelastic diffusion with one relaxation time which allows finite speeds of propagation of waves.Recently, Sherief and Saleh [16] investigated the problem of a thermoelastic half-space in the context of the theory of generalized thermoelastic diffusion with one relaxation time. Singh [17] discussed the reflection phenomena of waves from free surface of a thermoelastic diffusion with one relaxation time and with two relaxation times in [18]. Aouadi [19–21] investigated different problems in thermoelastic diffusion. Sharma and Walia [22, 23] discussed the effect of rotation on Rayleigh waves in the piezothermoelastic half-space. Kumar and Kansal [24] discussed the propagation of Rayleigh waves on free surface in transversely isotropic thermoelastic diffusion. Kumar and Kansal [25] derived the basic equations for generalized thermoelastic diffusion and discussed the Lamb waves. Dawn and Chakraborty [26] studied Rayleigh waves in Green-Lindsay’s model of generalized thermoelastic media. Kumar and Chawla [27] investigated the effect of rotation and stiffness on surface waves propagation in an elastic layer lying over a generalized thermodiffusive elastic half-space with imperfect boundary. New contributions on waves propagation in thermoelastic media have been discussed [28–32].In this paper, linear model is adopted to represent the imperfectly bonded interface conditions. The linear model is simplified and idealized situation of imperfectly bonded interface, where the discontinuities in displacements at interfaces have a linear relationship with the interface stresses. Taking these applications into account, the surface waves propagation at imperfect boundary between an isotropic elastic layer and isotropic thermodiffusive elastic half-space with rotation in the context of Green-Lindsay theory is investigated. The phase velocity and attenuation coefficients of wave propagation have been computed from the secular equations. The amplitudes of displacements, temperature, concentration, and specific loss are computed and depicted graphically to make clear the influence of magnetic field, rotation, stiffness, relaxation times, and diffusion on the phenomena and compare with the practical results.


## 2. Basic Equations

The basic governing equations for homogenous generalized thermodiffusive solid in the absence of heat and mass diffusion sources are as follows (Singh [17]).(i)
Constitutive relations are(1)σji=2μεji+λεkk-β1T+τ1T˙-β2C+τ′C˙δji.(ii)
Equation of motion in the rotating frame of reference is(2)μui,jj+λ+μuj,ij-β1T+τ1T˙,i-β2C+τ′C˙,i+Fi=ρu→¨+Ω→×Ω→×u→+2Ω→×u→.i,where Ω→×Ω→×u→ is the centripetal acceleration due to the time varying motion only and 2Ω→×Ω→. is the Coriolis acceleration:(3)Ω→×Ω→×u→=-Ω2u,0,-Ω2w=-Ω2u,0,w2Ω→×u→.=2Ωw˙,0,-Ωu˙=2Ωw˙,0,-u˙,where(4)F→=J→×B→.Consider that the medium is a perfect electric conductor; we take the linearized Maxwell equations governing the electromagnetic field, taking into account absence of the displacement current (SI): (5)curlh→=J→,curlE→=-μe∂h→∂t,divh→=0,divE→=0,where (6)h→=curlu→×H0→,where we have used(7)H→=H0→+h→x,z,t,H0→=0,H,0(8)B→=μeH→=μeHJ→;then(9)Fx=BJz=μeHH∂2u∂x12+∂2w∂x1∂x3=μeH2∂2u∂x12+∂2w∂x1∂x3,Fz=BJx=μeHH∂2u∂x3∂x1+∂2w∂x32=μeH2∂2u∂x3∂x1+∂2w∂x32.(iii)
Equation of heat conduction is(10)ρcET˙+τ0T¨+β1T0ε˙kk+aT0C˙+τ0C¨=KT,ii.(iv)
Equation of mass diffusion is(11)Dβ2εkk,ii+DaT+τ1∂T∂t,ii+C˙-DbC+τ′∂C∂t,ii=0.Here, the medium is rotating with angular velocityΩ=Ωn_, where n_ is the unit vector along the axis of rotation and this equation of motion includes two additional terms, namely,(i)
the centripetal accelerationΩ→×Ω→×u→ due to time-varying motion,(ii)
the Carioles acceleration2Ω→×u→.,whereβ1=3λ+2μαt and β2=3λ+2μαc, λ and μ are Lame’s constants, αt is the coefficient of linear thermal expansion, ρ and CE are, respectively, the density and specific heat at constant strain, a, b are, respectively, coefficient describing the measure of thermoelastic diffusion effects and of diffusion effects, T0 is the reference temperature assumed to be such that T/T0≪1, τ0,τ1 are thermal relaxation times with τ1≥τ0>0 and τ0, τ′ are diffusion relaxation times with τ′≥τ0>0, and ui are components of displacement vector. T(x,y,z) is the temperature change and C is the concentration; σij=σji, εij=ui,j+uj,i/2 are, respectively, the components of stress and strain tensor.The symbols correspond to partial derivative and time derivative, respectively.Following Bullen [33], the equations of motion and constitutive relations in isotropic elastic medium are given by(12)λe+μeuj,ije+μeui,jje+Fie=ρeu¨ie,(13)σjie=λeθeδij+2μeεije,i=j=1,2,3,where(14)Θe=uk,ke,(15)εije=ui,je+uj,ie2,i=j=1,2,3,and ue=u1e,u2e,u3e is the displacement vector, ρe is the density of the isotropic medium and λe,μe are Lame’s constants, σjie=σije are components of stress tensor, and δij is the Kronecker delta.


## 3. Formulation of the Problem

As shown in Figure1, we consider an isotropic elastic layer (Medium M1) of thickness H overlaying a homogeneous, isotropic, generalized thermodiffusive elastic half-space in rotating frame of reference (Medium M2). The origin of the coordinate system (x,y,z) is taken at any point on the horizontal surface and x2-axis in the direction of wave propagation and x3-axis taking vertically downward into half-space, so that all particles on a line parallel to x2-axis are equally displaced. Therefore, all the field quantities will be independent of x2-axis coordinate. The interface between isotropic elastic layer and thermodiffusive elastic half-space with rotation has been taken at an imperfect boundary. The displacement vector u, temperature T, concentration C, and rotation for medium M2 are taken as(16)u→=u1,0,u3,Tx1,x3,t,Cx1,x3,t,Ω=0,Ω,0and displacement vector u→e for the layer (Medium M1) is taken as(17)u→e=u1e,0,u3e.We define the dimensionless quantities(18)xi′=ω1∗xiv1,ui′=ω1∗uiv1,t′=ω1∗t,τ1′=ω1∗τ,T′=β1Tρv12,C′=β2Cρv12,τ0′=ω1∗τ0,τ0′=ω1∗τ0,τ′′=ω1∗τ′,σij′=σija1T0,Ω′=Ωω1∗,Kn′=v1knβ1T0ω1∗,K1′=v1k1β1T0ω1∗,v12=λ+2μ+μeH2ρ,ω1∗=ρcEv12k1.Figure 1
Schematic of the problem.Upon introducing the quantities in (1)–(3), (8), and (14)–(16), after suppressing the primes, with the aid of (17) and (18), we obtain(19)u1,11+δ2u3,13+δ1u1,33-τΘT.1-τc′C,1=u¨1-Ω2u1+2Ωu˙3,δ1u3,11+δ2u1.31+u3,33-τθT.3-τc′C,3=u¨3-Ω2u3-2Ωu˙1,∇2T=τθ0T˙+χ1τf0C˙+χ2e˙q1∗∇2e+q2∗τθ′∇2T-q3∗τc′∇2+C˙=0,(20)(u1,11e+u3,13e)δ32+∇2u1eδ42=u¨1e,(u1,13e+u3,33e)δ32+∇2u3eδ42=u¨3e,where(21)δ1=μ(λ+2μ+μeH),δ2=λ+μ+μeH(λ+2μ+μeH),χ1=aT0v12β1ω1∗kβ2,χ2=T0β12ρkω1∗,q1∗=Dω1∗β22ρv12,q2∗=Dω1∗bv12,q3∗=Dω1∗β2aβ1v12,∇2=∂2∂x2+∂2∂Z2,e=∂u∂x+∂w∂z,δ42=ρev12μe,δ32=ρev12μeH2+λe+μe,δ52=ρev12λe+2μeμeH2,τθ′=1+τ1∂∂t,τθ0=1+τ0∂∂t,τc′=1+τ′∂∂t,τf0=1+τ0∂∂t,RH2=μeHλ+2μ.For an isotropic elastic layer, we introduce potential functionsΦ and Ψ through the relations(22)u1e=∂Φ∂x1-∂Ψ∂x3,u3e=∂Φ∂x3+∂Ψ∂x1.Substituting from (22) into (20), we have(23)Φ,11+Φ,33-δ52Φ¨=0,(24)Ψ,11+Ψ,33-δ42Ψ¨=0.


## 4. Solution of the Problem

To solve (19), (22), and (23), we assume the solution in the form(25)u1,u2,T,C,Φ,Ψ=1,W,S,R,P1,P2Uexpiξx1+mx3-ct,where c=ω/ξ is the nondimensional phase velocity, ω is the frequency, m is still parameter 1, and W, S, and R are, respectively, the amplitude ratio of u1, u3, T, and C with respect to u1. Substituting the values of u1, u3, T, and C from (25) in (19), we obtain(26)1+δ1m2-c21+Λ2+δ2m-2iΛc2W+icω-1τ11S+τ21R=0,δ2m-2Λic2+δ1+m2-c21+Λ2W+icω-1mτ11S+τ21R=0,χ2C+χ2CmW-iω-1c2χ1τ20R-iω-1iω1+m2+c2τ10S=0,q1∗1+m2+q1∗1+m2mw-iq2∗1+m2τ11cω-1S+iq3∗1+m2τ21cω-1+c3ω-2R=0,where(27)Λ=Ωω-1,τ11=1-iωτ1,τ21=1-iωτ′,τ10=1-iωτ0,τ20=1-iωτ0.The system of (26) has a nontrivial solution if the determinant of the coefficient [1,W,S,R]T vanishes, which yield to the following polynomial characteristic equation:(28)m8+A∗m6+B∗m4+C∗m2+D∗=0.The constantsA∗, B∗, C∗, and D∗ are given in Appendix A. The characteristic equation (28) is biquadratic in m2 and hence possesses four roots, mp2; p=1,2,3,4 corresponding to four roots; there exist three types of quasilongitudinal waves and one quasitransverse wave. The formal expression for displacement, temperature, and concentration satisfying the radiation condition Remp≥0 can be written as(29)u1,u3,T,C=∑P=141,n1p,n2p,n3pexpiξx+impz-ct.Substituting the values ofΦ and Ψ from (25) in (23) and (24) and with the aid of (22), we obtain(30)u1e=iξB1c5+B2s5+ξm6D1s6-D2c6expiξx-ct,u3e=ξm5B2c5-B1s5+iξD1c6+D2s6expiξx-ct,where(31)m5=c2δ52-1,m6=c2δ42-1,c5=cosξm5x3,c6=cosξm5x3,s5=sinξm5x3,s6=sinξm6x3,where AP  P=1,2,3,4, B1,B2,D1,D2 are arbitrary constants. The coupling constants n1p,n2p,n3p  p=1,2,3,4 are given in Appendix B.


## 5. Boundary Conditions

In this paper, linear model is adopted to represent the imperfectly bonded interface conditions. The boundary conditions are the vanishing of the normal stress; Maxwell’s electromagnetic stress tensorτij is given by(32)τij=μeHihj+Hjhi-Hkhkδijand tangential stress at free surface. The discontinuities in displacements have linear relations with stresses, continuity of normal, Maxwell’s electromagnetic stress tensor and tangential stress, vanishing of the gradient of temperature, and concentration at the interface between the isotropic elastic layer and isotropic thermodiffusive elastic half-space. Mathematically, these can be written as follows.(i)
Mechanical conditions:(33)σ33eM1+τ33e=0σ31eM1+τ31e=0,ppppppppppppppppppppppppx3=-H,σ33eM1+τ33e=σ33M2+τ33σ31eM1+τ31e=σ31M2+τ31σ33M1+τ33e=knu3M2-u3M1σ31M1+τ31e=ktu3M2-u3M1;ppppppppppppppipipppppppx3=0.(ii)
Temperature condition(34)∂T∂x3=0,x3=0.(iii)
Concentration condition(35)∂C∂x3=0,x3=0,wherekn and kt are the normal and transverse stiffness of layer which have dimension Nm-3.


## 6. Derivation of the Secular Equations

Substituting the value ofu1,u3,T,C,u1e,u3e from (29)-(30) into (33)–(35), with the aid of (1), (10) and (13)-(14), after simplification we obtain(36)b5tanξm5HΔ1+b5Δ2+b6s6c5Δ3-b6c6c5Δ4=0,where Δ1=Rij7×7 the entries Rij of the determinant are given in Appendix C and Δ2 is obtained by replacing the first column of Δ1 by [R11∗00R41∗R51∗R61∗R71∗]T,  Δ3 is obtained by replacing the second column of Δ2 by [R12∗00R42∗R52∗R62∗R72∗]T, and Δ4 is obtained by replacing the third column of Δ3 by [R13∗00R43∗R53∗R63∗R73∗]T. The entries of ΔP(P=1,2,3,4) are given in Appendix C.If we write(37)c-1=v-1+iω-1G,then ξ=F+iG where F=ω/v and G are real numbers. Also the roots of characteristic equations are in general complex.Hence, assume thatmp=pp+iqp so that exponent in the plane wave solutions in (25) becomes(38)iFx1-mpAx3-vt-FGFx1+mpRx3,where(39)mpA=GFpp+qp,mpR=pp-iGFqp.This shows thatv is the propagation velocity and G is the attenuation coefficient of the wave. Upon using the representation (37) in secular equation (36), the values of propagation speed v and attenuation coefficient G of wave propagation can be obtained.


## 7. Particular Cases

(i) Normal Stiffness. In this case, kn≠0,  kt→∞, and the secular equation (36) remains the same. But the following will be replaced in the values of Δpp=1,2,3,4:(40)R71=-iξβ1T0,R72=0,R73=ξm6β1T0,R74=R75=R76=R77=1,R71∗=0,R72∗=-iξ,R73∗=0.(ii) Tangential Stiffness. In this case kt≠0, kn→∞, and the secular equation (36) remains the same with the change of values of Δp(p=1,2,3,4) by taking(41)R61=-iξβ1T0,R62=0,R63=ξm61β1T0,R64=n12,R65=n11,R66=n13,R67=n14,R61∗=0,R62∗=-ξm51β1T0,R63∗=-iξ1β1T0.(iii) Welded Contact. In this case, kn→∞, kt→∞, and the secular equation (36) remains the same, but the value of Δp(p=1,2,3,4) is given by replacing(42)R61=-iξ1β1T0,R62=0,R63=ξm61β1T0,R64=n12,R65=n11,R66=n13,R67=n14,R61∗=0,R72∗=-ξm51β1T0,R63∗=-iξ1β1T0,R61=-iξ1β1T0,R62=0,R63=ξm61β1T0R64=n12,R65=n11,R66=n13,R67=n14,R61∗=0,R72∗=-ξm51β1T0,R63∗=-iξ1β1T0.


## 8. Special Cases

Case (i). If we take Ω=0, that is, in the absence of rotation effect, the frequency equation (36) will reduce to the frequency equation for an isotropic elastic layer and a homogenous isotropic thermodiffusive elastic half-space without rotation.Case (ii). If we take H=0,  Ω=0, τ′>0, and τ1>0 and in the absence of diffusion effect, that is, b1=b3=a=b=0, (36) will reduce to the frequency equation for Rayleigh wave(43)2-c2δ122m-12+m-22+m-1m-2-1+c2-4m-1m-2m-3m-1+m-2=0.Here, m3=1-c2/δ12 and m-pp=1,2 are the roots of (28), obtained by taking b1=b3=a=b=0 and τ′>0, τ1>0. The resulting equation (43) is similar to (22) as given by Dawn and Chakraborty [26].Case (iii). In the absence of isotropic elastic layer and thermal and diffusion effects, we obtain the frequency equation corresponding to isotropic elastic half-space by changing the dimensionless quantities into physical quantities as(44)2-c2c222=41-c2c121/21-c2c221/2,where c12=λ+μ+μe/ρ, c22=μ/ρ.The frequency equation (36) is the same as derived in Ewing et al. [34].


## 9. Surface Displacements, Temperature Change, and Concentration

The amplitude of surface displacements, temperature change, and concentration at the surfacez=0 during Rayleigh wave propagation in the cases of stress free, vanishing of the gradient of temperature, and concentration of the half-space are(45)u1=E-AexpiFx1-vt,u3=R-AexpiFx1-vt,T=W-AexpiFx1-vt,C=L-AexpiFx1-vt,where(46)A=A1exp-Fx1,E¯=F1∗-F2∗+F3∗-F4∗F1∗,R¯=n11F1∗-n12F2∗+n13F3∗-n14F4∗F1∗,W¯=n21F1∗-n22F2∗+n23F3∗-n24F4∗F1∗,L¯=n31F1∗-n32F2∗+n33F3∗-n34F4∗F1∗and FP∗ (P=1,2,3,4) are given in Appendix D.


## 10. Specific Loss

The specific loss is the ratio of energy(Δw) dissipated in taking a specimen through a stress cycle, to the elastic energy (w) stored in the specimen when the strain is maximum. The specific loss is the most direct method of defining internal friction for a material. For a sinusoidal plane wave of small amplitude, Kolsky [35] shows that the specific loss Δw/w equals 4π times the absolute value of imaginary part of ξ to the real part of ξ; that is,(47)Δww=4πImξReξ=4πiGF=4πiGω/v=4πivGω=4π1FiG.


## 11. Numerical Results and Discussion

Following Sherief and Saleh [16], we take the following values of relevant parameter for the copper material:(48)λ=7.76×1010Kgm-1s-2,μ=3.86×1010Kgm-1s-2,T0=0.293×103Kρ=8.954×103kgm-3,a=1.2×104m2s-2k-1,b=9×105m5s-2gk-1cE=3831×103JKg-1K-1,αt=1.78×10-5k-1,αc=1.98×10-4k-1m3gK=0.383×103Wm-1k-1,D=0.85×10-8Sm-3g-1,τ0=0.07s,τ1=0.03,τ0=0.04,τ1=0.07s.The elastic constants for granite are given by Bullen [33]:(49)λe=μe=2.238×103J·Kgm-1K-1,ρe=2.65×103J·kgm-3.The concentration change, phase velocity and attenuation coefficient of wave propagation, displacement, stresses, temperature in the context of Green Lindsay (GL) theory of thermoelastic diffusion with variations of magnetic field, and rotation in 2D and 3D have been computed for various values of nondimensional wave number and calculated numerically and represented graphically in Figures2–22.Figure 2
Variation of concentration with respect to the wave number with variation of rotation.Figure 3
Variation of phase velocity with respect to the wave number with variation of rotation.Figure 4
Variation of attenuation coefficient with respect to the wave number with variation of rotation.Figure 5
Variation of the displacementu with respect to the wave number with variation of rotation.Figure 6
Variation of normal stressσxx with respect to the wave number with variation of rotation.Figure 7
Variation of shear stressσxy with respect to the wave number with variation of rotation.Figure 8
Variation of the temperatureθ with respect to the wave number with variation of rotation.Figure 9
Variation of the concentration charge with respect to the wave number with variation of the magnetic field.Figure 10
Variation of phase velocity with respect to the wave number with variation of the magnetic field.Figure 11
Variation of the attenuation coefficient with respect to the wave number with variation of the magnetic field.Figure 12
Variation of the displacementu with respect to the wave number with variation of magnetic field.Figure 13
Variation of shear stressσxx with respect to the wave number with variation of the magnetic field.Figure 14
Variation of shear stressσxy with respect to the wave number with variation of the magnetic field.Figure 15
Variation of the temperatureθ with respect to the wave number with variation of magnetic field.Figure 16
Variation of the concentration charge with respect to the wave number andx-axis.Figure 17
Variation of the phase velocity with respect to the wave number andx-axis.Figure 18
Variation of the attenuation coefficient with respect to the wave number andx-axis.Figure 19
Variation of displacementu with respect to the wave number and x-axis.Figure 20
Variation of the stressσxx with respect to the wave number and x-axis.Figure 21
Variation of the stressσxy with respect to the wave number and x-axis.Figure 22
Variation of the temperature with respect to the wave number andx-axis.Figure2 displays the variation of concentration change with different values of rotation with respect to the wave number; it appears that the concentration change increases with an increasing of the wave number and if there is no rotation and small wave number takes larger values than in the presence of rotation and large values of the wave number. Figures 3 and 4 clear the phase velocity and attenuation coefficients variations with respect to the wave number with variation of the rotation; it is concluded that the phase velocity begins from 1.2 but the attenuation coefficient begins from zero if the wave number equals zero and tends to zero if the wave number tends to infinity; also, it is seen that they increase and decrease periodically with an increasing of the wave number and rotation that indicate the interruption of the phase velocity and attenuation coefficient with the largest values of the wave number and rotation that agree with the practical results. Figure 5 shows the variation of the displacement component with respect to the wave number with variation of the rotation; it is concluded that it takes various values for zeros value of the wave number and tends to zero if the wave number tends to infinity; also, it seems that it increases and decreases periodically with an increasing of the wave number and rotation.Figures6 and 7 make clear the variation of the stresses components with respect to the wave number with variation of the rotation; it appears that it begins from zero for zeros value of the wave number and tends to zero if the wave number tends to infinity; also, it is shown that it increases and decreases periodically with an increasing of the wave number and rotation. Figure 8 displays the variation of the temperature with various values of the wave number and rotation; one can see that it begins from zero for zeros wave number and tends to zero, decreases and increases tends to zero if the wave number tends to infinity; it is clear that the temperature decreases with increasing rotation values that physically indicates the negative influence of rotation on the temperature that takes into consideration engineering and structures. Also, it is seen that, in the absence of rotation, all values of the concentration change, phase velocity, attenuation coefficient, displacement, stresses components, and temperature take a smooth behavior comparing with the corresponding values in the presence of rotation that take harmonic behavior.Figure9 shows the variation of concentration change with various values of magnetic field with respect to the wave number; it appears that the concentration change increases with an increasing of the wave number and decreases and arrives to the unity if the wave number tends to infinity; if the magnetic field is absent, the concentration change smallest than the corresponding values for small values of the wave number and inverses to the largest values for the large values of the wave number. Figures 10 and 11 make obvious the phase velocity and attenuation coefficients variations with respect to the wave number with variation of the magnetic field; it is concluded that the phase velocity begins from 1.2 but the attenuation coefficient begins from zero if the wave number equals zero and tends to zero if the wave number tends to infinity; also, it appears that they increase and decrease periodically with an increasing of the wave number and magnetic field that indicate the interruption of the phase velocity and attenuation coefficient with the largest values of the wave number and magnetic field that agree with the practical results. From Figure 12, it is concluded that the displacement component takes various values for zeros value of the wave number and tends to zero if the wave number tends to infinity; also, it is shown that it increases and decreases with an increasing of the wave number and magnetic field.Figures13 and 14 display the variation of the normal and shear stresses with respect to the wave number with variation of the magnetic field; it appears that it begins from zero for zeros value of the wave number and tends to zero if the wave number tends to infinity; also, it is shown that the normal stress decreases and then increases arriving to zero with the large values of the wave number and magnetic field. It is seen that the shear stress increases and periodically decreases with the variation of the wave number and magnetic field and begins from zero at zeros value of the wave number and tends again to zero if the wave number tends to infinity. Figure 15 dispatches the variation of the temperature with various values of the wave number and magnetic field; we can show that it begins from zero for zeros wave number and tends to zero for infinity wave number and increases and decreases with variation of the wave number and magnetic field.Totally, it is clear that the temperature decreases with increasing rotation values that physically indicates the negative influence of rotation on the temperature that takes into consideration engineering and structures.Figures16–22 display the 3D variation on the concentration change, phase velocity, attenuation coefficient, displacement, stresses components, and temperature, respectively, with respect to the wave number and x-axis if Ho=105 and Ω=0.1.It is obvious from Figure16 that the concentration change decreases to arrive to zero with the increased values of the wave number; it increases and decreases to its minimum value and after that increases with an increasing of x-axis. From Figures 17 and 18, it appears that the phase velocity begins from 1.2 but the attenuation coefficient begins from zero if the wave number equals zero and tends to zero if the wave number tends to infinity; also, it appears that it increases and decreases periodically and arrives to zero with an increasing of x-axis and the wave number. Figure 19 shows that the displacement decreases, increases periodically with the variation of the wave number, and tends to zero as the wave number tends to infinity. It appears that the displacement decreases and increases with the increased values of x-axis.From Figure20, we concluded that the normal stress decreases and then increases arriving to zero as the wave number tends to zero and increases with the small values of x-axis and after that decreases to tend to zero. Figure 21 shows that the shear stress component increases and decreases periodically with the variation of the wave number but increases and decreases with an increasing of x-axis.Finally, it appears that the temperature increases and decreases arrive to zero as the wave number tends to infinity but decreases with an increasing of the smallx-axis and after that increases.


## 12. Concluding Remarks

Surface waves at imperfect boundary between isotropic elastic layer of finite thickness and isotropic thermodiffusive elastic half-space with magnetic field, stiffness, and rotation with two thermal relaxation times (GL) model are illustrated. The secular equation in compact form has been derived. The concentration change, phase velocity, attenuation coefficient, displacement, stresses, and temperature are displayed graphically. The amplitudes of displacements, temperature, and concentration are computed at the free plane boundary and presented graphically. Specific loss of energy is obtained and depicted graphically.The analysis to be carried will be useful in the design and construction of rotating sensors, engineering, structures, and surface acoustic waves devices and the following remarks have been concluded.(1)
If there is no rotation, small wave number takes larger values than in the presence of rotation and large values of the wave number.(2)
The phase velocity, attenuation coefficient, displacement, stresses components, and temperature begin from zero for zeros value of the wave number and tend to zero if the wave number tends to infinity; also, it is seen that they increase and decrease periodically with an increasing of the wave number and rotation and magnetic field.(3)
The temperature decreases with increasing rotation values that physically indicates the negative influence of rotation on the temperature that takes into consideration engineering and structures and acoustic and rotating sensors.



---

*Source: 671783-2015-08-12.xml*