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.

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 [1013] 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 [1921] 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 [2832].

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(ii)Equation of motion in the rotating frame of reference is where is the centripetal acceleration due to the time varying motion only and is the Coriolis acceleration:whereConsider 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): where where we have usedthen(iii)Equation of heat conduction is(iv)Equation of mass diffusion is

Here, the medium is rotating with angular velocity , where is the unit vector along the axis of rotation and this equation of motion includes two additional terms, namely,(i)the centripetal acceleration due to time-varying motion,(ii)the Carioles acceleration ,

where and , and are Lame’s constants, is the coefficient of linear thermal expansion, and are, respectively, the density and specific heat at constant strain, , are, respectively, coefficient describing the measure of thermoelastic diffusion effects and of diffusion effects, is the reference temperature assumed to be such that , , are thermal relaxation times with and , are diffusion relaxation times with , and are components of displacement vector. is the temperature change and is the concentration; , 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 bywhereand is the displacement vector, is the density of the isotropic medium and are Lame’s constants, are components of stress tensor, and is the Kronecker delta.

3. Formulation of the Problem

As shown in Figure 1, we consider an isotropic elastic layer (Medium ) of thickness overlaying a homogeneous, isotropic, generalized thermodiffusive elastic half-space in rotating frame of reference (Medium ). The origin of the coordinate system is taken at any point on the horizontal surface and -axis in the direction of wave propagation and -axis taking vertically downward into half-space, so that all particles on a line parallel to -axis are equally displaced. Therefore, all the field quantities will be independent of -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 , temperature , concentration , and rotation for medium are taken asand displacement vector for the layer (Medium ) is taken asWe define the dimensionless quantities

Upon introducing the quantities in (1)–(3), (8), and (14)–(16), after suppressing the primes, with the aid of (17) and (18), we obtainwhere

For an isotropic elastic layer, we introduce potential functions and through the relationsSubstituting from (22) into (20), we have

4. Solution of the Problem

To solve (19), (22), and (23), we assume the solution in the formwhere is the nondimensional phase velocity, is the frequency, is still parameter , and , , and are, respectively, the amplitude ratio of , , , and with respect to . Substituting the values of , , , and  from (25) in (19), we obtainwhereThe system of (26) has a nontrivial solution if the determinant of the coefficient vanishes, which yield to the following polynomial characteristic equation:

The constants , , , and are given in Appendix A. The characteristic equation (28) is biquadratic in and hence possesses four roots, ; 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 can be written as

Substituting the values of and from (25) in (23) and (24) and with the aid of (22), we obtainwherewhere   , are arbitrary constants. The coupling constants    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 is given byand 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:(ii)Temperature condition(iii)Concentration condition

where and are the normal and transverse stiffness of layer which have dimension

6. Derivation of the Secular Equations

Substituting the value of from (29)-(30) into (33)–(35), with the aid of (1), (10) and (13)-(14), after simplification we obtainwhere the entries of the determinant are given in Appendix C and is obtained by replacing the first column of by ,   is obtained by replacing the second column of by , and is obtained by replacing the third column of by . The entries of are given in Appendix C.

If we writethen where and are real numbers. Also the roots of characteristic equations are in general complex.

Hence, assume that so that exponent in the plane wave solutions in (25) becomeswhere

This shows that is the propagation velocity and is the attenuation coefficient of the wave. Upon using the representation (37) in secular equation (36), the values of propagation speed and attenuation coefficient of wave propagation can be obtained.

7. Particular Cases

(i) Normal Stiffness. In this case, ,  , and the secular equation (36) remains the same. But the following will be replaced in the values of :

(ii) Tangential Stiffness. In this case , , and the secular equation (36) remains the same with the change of values of by taking

(iii) Welded Contact. In this case, , , and the secular equation (36) remains the same, but the value of is given by replacing

8. Special Cases

Case (i). If we take , 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 ,  , , and and in the absence of diffusion effect, that is, , (36) will reduce to the frequency equation for Rayleigh waveHere, and are the roots of (28), obtained by taking and , . 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 aswhere , .

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 surface during Rayleigh wave propagation in the cases of stress free, vanishing of the gradient of temperature, and concentration of the half-space arewhereand () are given in Appendix D.

10. Specific Loss

The specific loss is the ratio of energy dissipated in taking a specimen through a stress cycle, to the elastic energy 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 equals times the absolute value of imaginary part of to the real part of ; that is,

11. Numerical Results and Discussion

Following Sherief and Saleh [16], we take the following values of relevant parameter for the copper material:The elastic constants for granite are given by Bullen [33]:

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 Figures 222.

Figure 2 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.

Figures 6 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.

Figure 9 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.

Figures 13 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.

Figures 1622 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 -axis if and

It is obvious from Figure 16 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 -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 -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 -axis.

From Figure 20, 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 -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 -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 small -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.

Appendices

A. Constant of (28)

Consider the following:

B. Constant of (29)

Consider the following:where

C. Determinant Elements of (36)

Consider the following:

D. Given in (46)

Consider the following:

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.