# Surface Wave Propagation in a Microstretch Thermoelastic Diffusion Material under an Inviscid Liquid Layer

Authors: Rajneesh Kumar; Sanjeev Ahuja; S. K. Garg
Journal: Advances in Acoustics and Vibration
 (2014)

---

## Abstract

The present investigation deals with the propagation of Rayleigh type surface waves in an isotropic microstretch thermoelastic diffusion solid half space under a layer of inviscid liquid. The secular equation for surface waves in compact form is derived after developing the mathematical model. The dispersion curves giving the phase velocity and attenuation coefficients with wave number are plotted graphically to depict the effect of an imperfect boundary alongwith the relaxation times in a microstretch thermoelastic diffusion solid half space under a homogeneous inviscid liquid layer for thermally insulated, impermeable boundaries and isothermal, isoconcentrated boundaries, respectively. In addition, normal velocity component is also plotted in the liquid layer. Several cases of interest under different conditions are also deduced and discussed.


---

## Body

## 1. Introduction

A micropolar continuum is a collection of interconnected particles in the form of small rigid bodies undergoing both translational and rotational motions. Typical examples of such materials are granular media and multimolecular bodies whose microstructures act as an evident part in their macroscopic responses. Rigid chopped fibres, elastic solids with rigid granular inclusions, and other industrial materials such as liquid crystals are examples of such materials.The theory of micropolar elasticity introduced and developed by Eringen [1] aroused much interest because of its possible utility in investigating the deformation properties of solids for which the classical theory is inadequate. The micropolar theory is believed to be particularly useful in investigating materials consisting of bar-like molecules which exhibit microrotation effects and which can support body and surface couples. Eringen [2] formulated the theory of micropolar fluids which display the effects of local rotary inertia and couple stresses. This theory can be used to explain the flow of colloidal fluids, liquid crystals, animal blood, and so forth. Eringen [3] and Nowacki [4] extended the theory of micropolar elasticity to heat conducting elastic solids by including thermal effects. Chandrasekharaiah [5] developed a heat-flux dependent generalized theory of micropolar thermoelasticity.Eringen [6] developed the theory of micropolar elastic solid with stretch. Eringen [7] also developed the theory of microstretch thermoelastic solid. A microstretch elastic solid possesses seven degrees of freedom: three for translation, three for rotation, and one for stretch. The material points of microstretch bodies can stretch and contract independently of their rotation and translation. Composite materials reinforced with microcracks, animal bones, foams, chopped elastic fibres, porous media whose pores are filled with gas or inviscid liquid, asphalt or other elastic inclusions, “solid-liquid” crystals, and so forth fall outside the domain of micropolar elasticity; the microstretch continua can be used as mathematical models for these types of solids. A comprehensive review on the micropolar continuum theory has been given in the book by Eringen [8].The thermodiffusion in elastic solids is due to coupling of fields of temperature, mass diffusion, and that of strain in addition to heat and mass exchange with the environment. Nowacki [9–12] developed the theory of thermoelastic diffusion by using coupled thermoelastic model. Dudziak and Kowalski [13] and Olesiak and Pyryev [14], respectively, discussed the theory of thermodiffusion and coupled quasistationary problems of thermal diffusion for an elastic layer. Uniqueness and reciprocity theorems for the equations of generalized thermoelastic diffusion problem, in isotropic media were proved by Sherief et al. [15]. Kumar and Kansal [16] developed the basic equation of anisotropic thermoelastic diffusion based on Green-Lindsay model.Surface waves carry a lot of information about the Earth’s crust and dispersion analysis of surface waves is concerning with the phase velocity and wave number. Further the velocity of seismic waves depends on the elastic parameters of the medium through which they travel. Therefore, many investigators have studied the surface wave propagation by considering different kinds of models. Wu and Zhu [17] studied the propagation of Lamb waves in a plate bordered with inviscid liquid layer on both sides. Zhu and Wu [18] derived the dispersion equations of the Lamb waves of a plate bordered with viscous liquid layer or half space viscous liquid on both sides.Several authors discussed different problems in microstretch thermoelastic solid, notable among them are the spatial decay for the dynamic problems of thermomicrostretch elastic solid developed by Quintanilla [19]. Kumar and Partap discussed [20] the free vibration analysis of waves in a homogeneous isotropic microstretch elastic plate subjected to stress free conditions. Abbas and Othman [21] studied the propagation of plane waves in generalized thermomicrostretch elastic solid with thermal relaxation using finite element method. Othman et al. [22] studied the effect of magnetic field on the general model of equations of rotating generalized thermomicrostretch for a homogeneous isotropic elastic half-space solid whose surface is subjected to a Mode-I Crack in the context of the Green and Naghdi theory.Propagation of Rayleigh surface waves in microstretch thermoelastic continua under inviscid fluid loadings has been investigated by Sharma et al. [23]. The propagation of free vibrations in microstretch thermoelastic homogeneous isotropic, thermally conducting plate bordered with layers of inviscid liquid on both sides subjected to stress free thermally insulated and isothermal conditions has been investigated by Kumar and Partap [24]. Recently, Kumar et al. studied the propagation of Rayleigh waves in an isotropic microstretch thermoelastic diffusion solid half-space [25].However, no attempt has been made to discuss the effect of an imperfect boundary along with the relaxation times in a microstretch thermoelastic diffusion solid half space under a homogeneous inviscid liquid layer for thermally insulated, impermeable boundaries and isothermal, isoconcentrated boundaries, respectively.The exact nature of the layers beneath the earth’s surface is unknown. Therefore, for the purpose of theoretical investigations, one has to consider various appropriate theoretical models. The model under consideration is the motivation of the situation when inviscid liquid layer is resting on the microstretch thermoelastic diffusion medium. Since water, oil, chemicals, and so forth are present inside the earth, the propagation of surface waves in such a layer may be relevant in the exploration of oils and other valuable liquids. The problem of surface wave propagation as a part of exploration seismology is also helpful in various economical activities, like tracing of hydrocarbons and other mineral ores. The present study, in fact, is a step to attempt a more realistic model of the earth’s crust and hence may be helpful in the further investigation of exploration seismology, geophysics, earthquake engineering, and soil dynamics; both theoretically and practically, we hope this study will enhance the knowledge in better understanding of the complexities of the earth medium indirectly.In the present paper, the governing equations of isotropic microstretch thermoelastic diffusion solid half space under a layer of inviscid liquid are formulated and solved for two-dimension case. The phase velocity and attenuation coefficients of the waves are studied numerically for a particular model. The study of present investigation has applications in oil companies which rely on seismic interpretation for selecting the sites for exploratory oil wells because seismic wave methods have higher precision, higher resolution, and most cost-effective technique, as compared to drilling which is expansive and time consuming.


## 2. Basic Equations

Following Eringen [8], Sherief et al. [15], and Kumar and Kansal [16], the equations of motion and the constitutive relations in a homogeneous isotropic microstretch thermoelastic diffusion solid in the absence of body forces, body couples, stretch force, and heat sources are given by
(1)
(
λ
+
2
μ
+
K
)
∇
(
∇
·
u
→
)
-
(
μ
+
K
)
∇
×
∇
×
u
→
+
K
∇
×
φ
→
+
λ
o
∇
φ
*
-
β
1
(
1
+
τ
1
∂
∂
t
)
∇
T
-
β
2
(
1
+
τ
1
∂
∂
t
)
∇
C
=
ρ
∂
2
u
→
∂
t
2
,
(
α
+
β
+
γ
)
∇
(
∇
·
φ
→
)
-
γ
∇
×
(
∇
×
φ
→
)
+
K
∇
×
u
→
-
2
K
φ
→
=
ρ
j
∂
2
φ
→
∂
t
2
,
α
0
∇
2
φ
*
+
ν
1
(
T
+
τ
1
T
˙
)
+
ν
2
(
C
+
τ
1
C
˙
)
-
λ
1
φ
*
-
λ
0
∇
·
u
→
=
ρ
j
0
2
∂
2
φ
*
∂
t
2
,
K
*
∇
2
T
=
β
1
T
0
(
1
+
ε
τ
0
∂
∂
t
)
∇
·
u
→
˙
+
ν
1
T
0
(
1
+
ε
τ
0
∂
∂
t
)
φ
˙
*
+
ρ
C
*
(
1
+
τ
0
∂
∂
t
)
T
˙
+
a
T
0
(
C
˙
+
γ
1
C
¨
)
,
D
β
2
∇
2
(
∇
·
u
→
)
+
D
ν
2
∇
2
φ
*
+
D
a
∇
2
(
T
+
τ
1
T
˙
)
+
(
C
˙
+
ε
τ
0
C
¨
)
-
D
b
∇
2
(
C
+
τ
1
C
˙
)
=
0
,
and constitutive relations are
(2)
t
i
j
=
λ
u
r
,
r
δ
i
j
+
μ
(
u
i
,
j
+
u
j
,
i
)
+
K
(
u
j
,
i
-
ε
i
j
r
φ
r
)
+
λ
o
δ
i
j
φ
*
-
β
1
(
1
+
τ
1
∂
∂
t
)
T
δ
i
j
-
β
2
(
1
+
τ
1
∂
∂
t
)
C
δ
i
j
,
m
i
j
=
α
φ
r
,
r
δ
i
j
+
β
φ
i
,
j
+
γ
φ
j
,
i
+
b
0
ε
m
j
i
φ
,
m
*
,
λ
i
*
=
α
0
φ
,
i
*
+
b
0
ε
i
j
m
φ
j
,
m
,
where λ, μ are Lames elastic constants, α, β, γ, K are micropolar constants, λ
0, λ
1, α
o, b
o are microstretch elastic constants, ρ is the mass density, u
→
=
(
u
1
,
u
2
,
u
3
) is the displacement vector, φ
→
=
(
φ
1
,
φ
2
,
φ
3
) is the microrotation vector, φ
* is the scalar microstretch function, T and T
0 are the small temperature increment and the reference temperature of the body chosen such that |
T
/
T
0
|
≪
1, C is the concentration of the diffusion material in the elastic body, K
* is the coefficient of the thermal conductivity, C
* the specific heat at constant strain, D is the thermoelastic diffusion constant, a, b are, respectively, coefficients describing the measure of thermodiffusion and of mass diffusion effects, β
1
=
(
3
λ
+
2
μ
+
K
)
α
t
1, β
2
=
(
3
λ
+
2
μ
+
K
)
α
c
1, ν
1
=
(
3
λ
+
2
μ
+
K
)
α
t
2, ν
2
=
(
3
λ
+
2
μ
+
K
)
α
c
2, α
t
1, α
t
2 are coefficients of linear thermal expansion, and α
c
1, α
c
2 are the coefficients of linear diffusion expansion. j is the microintertia, j
o is the microinertia of the microelements, t
i
j and m
i
j are components of stress and couple stress tensors, respectively, λ
i
* is the microstress tensor, e
i
j
(
=
(
1
/
2
)
(
u
i
,
j
+
u
j
,
i
)
) are components of infinitesimal strain, e
k
k is the dilatation, δ
i
j is the Kronecker delta, τ
0, τ
1 are diffusion relaxation times with τ
1
≥
τ
0
≥
0, and τ
0, τ
1 are thermal relaxation times with τ
1
≥
τ
0
≥
0. Here τ
0
=
τ
0
=
τ
1
=
τ
1
=
γ
1
=
0 for coupled thermoelastic (CT) model, τ
1
=
τ
1
=
0, ε
=
1, γ
1
=
τ
0 for Lord-Shulman (L-S) model, and ε
=
0, γ
1
=
τ
0 where τ
0
>
0 for Green-Lindsay (G-L) model.In the above equations, a comma followed by a suffix denotes spatial derivative and a superposed dot denotes the derivative with respect to time, respectively.Following Achenbach [26], the field equations can be expressed in terms of velocity potential ψ
l for inviscid fluid as
(3)
∂
2
ψ
l
∂
x
1
2
+
∂
2
ψ
l
∂
x
3
2
=
1
(
α
l
)
2
∂
2
ψ
l
∂
t
2
,
where α
l
=
λ
l
/
ρ
l is the velocity of acoustic fluid, λ
l is the bulk modulus, and ρ
l is the density of the fluid.The velocity componentsu
1
l, u
3
l and the acoustic pressure p
l in the medium M
2 are given by
(4)
u
1
l
=
∂
ψ
l
∂
x
1
,
u
3
l
=
∂
ψ
l
∂
x
3
,
p
l
=
-
ρ
l
∂
ψ
l
∂
t
.


## 3. Formulation of the Problem

We consider a homogeneous, isotropic microstretch thermoelastic diffusion medium (M
1) lying under a uniform homogeneous inviscid liquid layer medium (M
2) of thickness H (Figure 1). We take the rectangular Cartesian coordinate system O
x
1
x
2
x
3 at any point on the plane horizontal surface and x
1-axis in the direction of wave propagation and x
3-axis pointing vertically downward into the half-space so that all particles on a line parallel to x
2-axis are equally displaced; therefore, all the field quantities will be independent of x
2-coordinates. Hence the liquid layer medium M
2 occupies the region -
H
<
x
3
<
0 and x
3
>
0 is occupied by the half-space medium M
1. The plane x
3
=
-
H represents the free surface of the liquid and x
3
=
0 is taken as interface between M
1 and M
2 medium. For the two-dimensional problem, we take
(5)
u
→
(
x
1
,
x
3
,
t
)
=
(
u
1
,
0
,
u
3
)
,
φ
→
=
(
0
,
φ
2
,
0
)
,
φ
*
(
x
1
,
x
3
,
t
)
,
T
(
x
1
,
x
3
,
t
)
,
C
(
x
1
,
x
3
,
t
)
.
We define the following dimensionless quantities:
(6)
(
x
1
′
,
x
3
′
)
=
ω
*
c
1
(
x
1
,
x
3
)
,
(
u
1
′
,
u
3
′
)
=
ρ
c
1
ω
*
β
1
T
o
(
u
1
,
u
3
)
,
φ
2
′
=
ρ
c
1
2
β
1
T
o
φ
2
,
φ
*
′
=
ρ
c
1
2
β
1
T
o
φ
*
,
t
i
j
′
=
t
i
j
β
1
T
o
,
m
i
j
′
=
ω
*
c
1
β
1
T
o
m
i
j
,
λ
i
*
′
=
λ
i
*
ω
*
c
1
β
1
T
o
,
T
′
=
T
T
o
,
C
′
=
β
2
C
ρ
c
1
2
,
t
′
=
ω
*
t
,
τ
o
′
=
ω
*
τ
o
,
τ
0
′
=
ω
*
τ
0
,
τ
1
′
=
ω
*
τ
1
,
τ
1
′
=
ω
*
τ
1
,
h
′
=
c
1
ω
*
h
,
and
(7)
(
u
1
'
l
,
u
3
'
l
)
=
1
c
1
(
u
1
l
,
u
3
l
)
,
p
l
′
=
p
l
μ
,
ψ
l
'
=
ω
*
ψ
l
c
1
2
,
K
n
′
=
c
1
μ
K
n
,
where ω
*
=
ρ
C
*
c
1
2
/
K
*, c
1
2
=
(
λ
+
2
μ
+
K
)
/
ρ, ω
* is the characteristic frequency of the medium.Figure 1
Geometry of the problem.Upon introducing the quantities (6) and (7) in (1), with the aid of (5) and after suppressing the primes, we obtain
(8)
δ
2
∂
e
∂
x
1
+
(
1
-
δ
2
)
∇
2
u
1
-
a
1
∂
φ
2
∂
x
3
+
a
2
∂
φ
*
∂
x
1
-
τ
t
1
∂
T
∂
x
1
-
a
3
τ
c
1
∂
C
∂
x
1
=
∂
2
u
1
∂
t
2
,
δ
2
∂
e
∂
x
3
+
(
1
-
δ
2
)
∇
2
u
3
+
a
1
∂
φ
2
∂
x
1
+
a
2
∂
φ
*
∂
x
3
-
τ
t
1
∂
T
∂
x
3
-
a
3
τ
c
1
∂
C
∂
x
3
=
∂
2
u
3
∂
t
2
,
a
4
∇
2
φ
2
+
a
5
(
∂
u
1
∂
x
3
-
∂
u
3
∂
x
1
)
-
a
6
φ
2
=
∂
2
φ
2
∂
t
2
,
(
δ
1
2
∇
2
-
a
7
)
φ
*
-
a
8
e
+
a
9
τ
t
1
T
+
a
10
τ
c
1
C
=
∂
2
φ
*
∂
t
2
,
∇
2
T
=
a
11
τ
e
0
∂
e
∂
t
+
a
12
τ
e
0
∂
φ
*
∂
t
+
τ
t
0
∂
T
∂
t
+
a
13
τ
c
0
∂
C
∂
t
,
a
14
∇
2
e
+
a
21
∇
2
φ
*
+
a
15
τ
t
1
∇
2
T
+
τ
f
0
∂
C
∂
t
-
a
16
τ
c
1
∇
2
C
=
0
,
where
(9)
(
a
1
,
a
2
)
=
1
ρ
c
1
2
(
K
,
λ
0
)
,
a
3
=
ρ
c
1
2
β
1
T
0
,
(
a
4
,
a
5
,
a
6
)
=
1
j
ρ
(
γ
c
1
2
,
K
ω
*
2
,
2
K
ω
*
2
)
,
δ
2
=
λ
+
μ
ρ
c
1
2
,
(
a
11
,
a
12
,
a
13
)
=
1
K
*
ω
*
(
T
0
β
1
2
ρ
,
β
1
T
0
ν
1
ρ
,
ρ
c
1
4
a
β
2
)
,
(
a
14
,
a
15
,
a
16
)
=
D
ω
*
c
1
2
(
β
1
2
ρ
c
1
2
,
β
2
a
β
1
,
b
)
,
(
a
7
,
a
8
,
a
9
,
a
10
)
=
2
j
0
ω
*
2
(
λ
1
ρ
,
λ
0
ρ
,
ν
1
c
1
2
β
1
,
ν
2
ρ
c
1
4
β
1
β
2
T
0
)
,
δ
1
2
=
c
2
2
c
1
2
,
c
2
2
=
2
α
0
ρ
j
0
,
a
21
=
D
ν
β
2
ω
*
ρ
c
1
4
,
τ
t
1
=
1
+
τ
1
∂
∂
t
,
τ
c
1
=
1
+
τ
1
∂
∂
t
,
τ
f
0
=
1
+
ε
τ
0
∂
∂
t
,
τ
t
0
=
1
+
τ
0
∂
∂
t
,
τ
e
0
=
1
+
ε
τ
0
∂
∂
t
,
τ
c
0
=
1
+
γ
1
∂
∂
t
,
e
=
∂
u
1
∂
x
1
+
∂
u
3
∂
x
3
,
∇
2
=
∂
2
∂
x
1
2
+
∂
2
∂
x
3
2
.
We introduce the potential functions ϕ and ψ through the relations
(10)
u
1
=
∂
ϕ
∂
x
1
-
∂
ψ
∂
x
3
,
u
3
=
∂
ϕ
∂
x
3
+
∂
ψ
∂
x
1
.
In (8); we obtain
(11)
∇
2
ϕ
+
a
2
φ
*
-
τ
t
1
T
-
a
3
τ
c
1
C
=
ϕ
¨
,
(12)
(
1
-
δ
2
)
∇
2
ψ
+
a
1
φ
2
=
ψ
¨
,
(13)
(
a
4
∇
2
-
a
6
)
φ
2
-
a
5
∇
2
ψ
=
φ
¨
2
,
(14)
(
δ
1
2
∇
2
-
a
7
)
φ
*
-
a
8
∇
2
ϕ
+
a
9
τ
t
1
T
+
a
10
τ
c
1
C
=
φ
¨
*
,
(15)
∇
2
T
=
τ
e
0
(
a
11
∇
2
ϕ
˙
+
a
12
φ
˙
*
)
+
τ
t
0
T
˙
+
a
13
τ
c
0
C
˙
,
(16)
a
14
∇
4
ϕ
+
a
21
∇
2
φ
*
+
a
15
τ
t
1
∇
2
T
-
a
16
τ
c
1
∇
2
C
+
τ
f
0
C
˙
=
0
.
In order to solve (11)–(16), we assume
(17)
{
ϕ
,
φ
*
,
T
,
C
,
ψ
l
}
(
x
1
,
x
3
,
t
)
=
{
ϕ
-
,
φ
-
*
,
T
-
,
C
-
,
ψ
-
l
}
(
x
3
)
e
i
ξ
(
x
1
-
c
t
)
,
(18)
{
ψ
,
φ
2
}
(
x
1
,
x
3
,
t
)
=
{
ψ
-
,
φ
-
2
}
(
x
3
)
e
i
ξ
(
x
1
-
c
t
)
,
where ξ is the wave number, ω
=
ξ
c is the angular frequency, and c is phase velocity of the wave. Using (17) in (11) and (14)–(16), we obtain a system of four homogeneous equations in four unknowns ϕ
-, φ
-
*, T
-, and C
- which for the nontrivial solution yields
(19)
D
8
+
A
1
*
D
6
+
B
1
*
D
4
+
C
1
*
D
2
+
D
1
*
=
0
and using (18) in (12) and (13), we obtain a system of two homogeneous equations in two unknowns ψ
- and φ
-
2 which for the nontrivial solution yields
(20)
D
4
+
A
2
*
D
2
+
B
2
*
=
0
,
where D
=
d
/
d
x
3, and the coefficients A
1
*, B
1
*, C
1
*, D
1
*, A
2
*, and B
2
* are given in Appendix A.Let the roots of (19) be denoted by m
p
2  (
p
=
1,2
,
3,4
). Four positive values of c in the descending order will be the velocities of propagation of four possible waves, namely, longitudinal displacement wave (LD), thermal wave (T), mass diffusion wave (MD), and longitudinal microstretch wave (LM), respectively. Similarly, for two roots of (20) denoted by m
p
2  (
p
=
5,6
), the positive values of c in the descending order will be the velocities of propagation of two coupled transverse displacement and transverse microrotational waves (CD I, CD II), respectively. Since we are interested in surface waves only, it is essential that the motion is confined to the free surface of x
3
=
0 in the half-space so that the characteristic roots satisfy the radiation conditions Re
(
m
p
)
≥
0,  (
p
=
1,2
,
3,4
,
5,6
). Thus, the solution of field equations takes the form
(21)
{
ϕ
,
φ
*
,
T
,
C
}
=
∑
p
=
1
4
[
A
p
{
1
,
n
1
p
,
n
2
p
,
n
3
p
}
e
-
m
p
x
3
]
e
i
ξ
(
x
1
-
c
t
)
,
{
ψ
,
φ
2
}
=
∑
p
=
5
6
[
A
p
{
1
,
n
4
p
}
e
-
m
p
x
3
]
e
i
ξ
(
x
1
-
c
t
)
,
where A
p
(
p
=
1,2
,
3,4
,
5,6
) are arbitrary constants.The coupling constantsn
1
p, n
2
p, n
3
p, and n
4
p are given in Appendix B.Substituting the values ofψ
l from (17) in (3) and after simplification, we obtain
(22)
ψ
l
=
[
A
7
cos
⁡
⁡
(
m
7
x
3
)
+
A
8
sin
⁡
(
m
7
x
3
)
]
exp
⁡
⁡
(
i
ξ
(
x
1
-
c
t
)
)
,
where
(23)
m
7
=
±
ξ
(
c
1
2
c
2
(
α
l
)
2
-
1
)
1
/
2
and A
7, A
8 are arbitrary constants.The velocity component and acoustic pressure for mediumM
2, with the help of (3), (4), (7), and (22), are given by
(24)
u
3
l
=
m
7
[
-
A
7
sin
⁡
(
m
7
x
3
)
+
A
8
cos
⁡
⁡
(
m
7
x
3
)
]
×
exp
⁡
⁡
(
i
ξ
(
x
1
-
c
t
)
)
,
p
l
=
ε
l
[
A
7
cos
⁡
⁡
(
m
7
x
3
)
+
A
8
sin
⁡
(
m
7
x
3
)
]
×
exp
⁡
⁡
(
i
ξ
(
x
1
-
c
t
)
)
,
where ε
l
=
i
ρ
l
c
1
2
c
ξ
/
μ.


## 4. Boundary Conditions

The appropriate boundary conditions are as follows.(i)
Vanishing of the pressure at the free surface of the liquid layer,x
3
=
-
H; that is,
(25)
p
l
=
0
.
(ii)
Continuity of normal stress components at the interfacex
3
=
0; that is,
(26)
(
t
33
)
M
1
=
-
p
l
.
(iii)
Discontinuity of the pressure at the interfacex
3
=
0; that is,
(27)
p
l
=
K
n
(
(
u
˙
3
)
M
1
-
(
u
3
)
M
2
)
.
(iv)
Vanishing of shear stress for mediumM
1 at the interface x
3
=
0; that is,
(28)
(
t
31
)
M
1
=
0
.
(v)
Vanishing of the tangential couple stress at the interfacex
3
=
0; that is,
(29)
(
m
32
)
M
1
=
0
.
(vi)
Vanishing of the normal microstress component at the interfacex
3
=
0; that is,
(30)
(
λ
3
*
)
M
1
=
0
.
(vii)
Thermal boundary conditions(31)
∂
T
∂
x
3
+
h
1
T
=
0
,
where
{
h
1
⟶
0
corresponds
Insulated
boundary
condition
h
1
⟶
∞
corresponds
Isolated
boundary
condition
.
(viii)
Consider(32)
∂
C
∂
x
3
+
h
2
C
=
0
,
where
{
h
2
⟶
0
corresponds
to
impermeable
boundary
h
2
⟶
∞
refers
to
isoconcentrated
boundary
,where
(33)
t
33
=
∂
u
3
∂
x
3
+
b
1
∂
u
1
∂
x
1
-
τ
t
1
T
-
τ
c
1
C
+
a
2
φ
*
,
t
31
=
b
2
∂
u
1
∂
x
3
+
b
3
∂
u
3
∂
x
1
-
a
1
φ
2
,
m
32
=
b
4
∂
φ
2
∂
x
3
+
b
5
∂
φ
*
∂
x
1
,
λ
3
*
=
b
6
∂
φ
*
∂
x
3
-
b
5
∂
φ
2
∂
x
1
,
(
b
1
,
b
2
,
b
3
)
=
1
ρ
c
1
2
(
λ
,
μ
+
K
,
μ
)
,
(
b
4
,
b
5
,
b
6
)
=
ω
*
2
ρ
c
1
4
(
γ
,
b
o
,
α
o
)
.


## 5. Derivations of the Secular Equations

Making use of (21) and (24) in the (25)–(32), we obtain a system of eight simultaneous linear equations:
(34)
∑
p
=
1
8
k
1
p
A
p
=
0
,
∑
p
=
1
8
k
2
p
A
p
=
0
,
∑
p
=
1
8
k
3
p
A
p
=
0
,
∑
p
=
1
8
k
4
p
A
p
=
0
,
∑
p
=
1
8
k
5
p
A
p
=
0
,
∑
p
=
1
8
k
6
p
A
p
=
0
,
∑
p
=
1
8
k
7
p
A
p
=
0
,
∑
p
=
1
8
k
8
p
A
p
=
0
,
where
(35)
k
1
p
=
{
0
,
for
(
p
=
1,2
,
3,4
,
5,6
)
,
cos
⁡
⁡
(
m
7
H
)
,
for
(
p
=
7
)
,
-
sin
⁡
(
m
7
H
)
,
for
(
p
=
8
)
,
k
2
p
=
{
m
p
2
-
b
1
ξ
2
+
a
2
n
1
p
+
n
2
p
(
i
ξ
c
τ
1
-
1
)
+
n
3
p
(
i
ξ
c
τ
1
-
1
)
,
for
(
p
=
1,2
,
3,4
)
,
i
ξ
m
p
(
b
1
-
1
)
,
for
(
p
=
5,6
)
,
ε
l
,
for
(
p
=
7
)
,
0
,
for
(
p
=
8
)
,
k
3
p
=
{
i
ξ
c
K
n
m
p
a
3
,
for
(
p
=
1,2
,
3,4
)
,
K
n
ξ
2
c
a
3
,
for
(
p
=
5,6
)
,
ε
l
,
for
(
p
=
7
)
,
-
K
n
m
7
,
for
(
p
=
8
)
,
k
4
p
=
{
-
i
ξ
m
p
(
b
2
+
b
3
)
,
for
(
p
=
1,2
,
3,4
)
,
-
(
b
2
m
p
2
+
b
3
ξ
2
+
a
1
n
4
p
)
,
for
(
p
=
5,6
)
,
0
,
for
(
p
=
7,8
)
,
k
5
p
=
{
i
ξ
b
5
n
1
p
,
for
(
p
=
1,2
,
3,4
)
,
-
b
4
n
4
p
m
p
,
for
(
p
=
5,6
)
,
0
,
for
(
p
=
7,8
)
,
k
6
p
=
{
n
1
p
m
p
b
6
,
for
(
p
=
1,2
,
3,4
)
,
i
ξ
b
5
n
4
p
,
for
(
p
=
5,6
)
,
0
,
for
(
p
=
7,8
)
,
k
7
p
=
{
(
h
1
-
m
p
)
n
2
p
,
for
(
p
=
1,2
,
3,4
)
,
0
,
for
(
p
=
5,6
,
7,8
)
,
k
8
p
=
{
(
h
2
-
m
p
)
n
3
p
,
for
(
p
=
1,2
,
3,4
)
,
0
,
for
(
p
=
5,6
,
7,8
)
.
The system of (34) has a nontrivial solution if the determinant of amplitudes A
P,   (
p
=
1,2
,
3,4
,
5,6
,
7,8
) vanishes which leads to the secular equation
(36)
|
k
i
j
|
8
×
8
=
0
,
(
i
,
j
=
1,2
,
3,4
,
5,6
,
7,8
)
which can further be simplified as
(37)
c
2
=
(
α
l
c
1
)
2
(
[
1
ξ
H
tan
-
1
⁡
(
Δ
12
)
]
2
+
1
)
,
where Δ
12 is given in Appendix B.Equation (37) is Rayleigh type surface wave secular equation of an isotropic microstretch thermoelastic diffusion solid half-space under a uniform homogeneous inviscid liquid layer. The secular equation has complete information about the phase velocity, attenuation coefficient, and wave number of surface waves in such medium.


## 6. Particular Cases

(i)
If we take,K
n
→
∞ in (37), we obtained the frequency equation at the welded contact between an inviscid liquid layer and an isotropic microstretch thermoelastic diffusion solid half space with the changed values of k
i
j under
(38)
k
3
p
=
i
ξ
c
m
p
a
3
,
for
(
p
=
1,2
,
3,4
)
,
k
3
p
=
ξ
2
c
a
3
,
for
(
p
=
5,6
)
,
k
3
p
=
0
,
for
(
p
=
7
)
,
k
3
p
=
-
m
7
,
for
(
p
=
8
)
.
(ii)
TakingH
→
0 and K
n
→
∞, in secular equation (37), we obtained the frequency equation for Rayleigh waves propagating along the free surface of a microstretch thermoelastic diffusion solid half space as
(39)
|
k
21
k
22
k
23
k
24
k
25
k
26
k
41
k
42
k
43
k
44
k
45
k
46
k
51
k
52
k
53
k
54
k
55
k
56
k
61
k
62
k
63
k
64
k
65
k
66
k
71
k
72
k
73
k
74
k
75
k
76
k
81
k
82
k
83
k
84
k
85
k
86
|
6
×
6
,
with the changed values of k
i
j in (37), that is, k
3
p
=
0, for p
=
1,2
,
…
,
6; k
q
7
=
k
q
8
=
0, for q
=
2,3
,
…
,
8. The frequency equation (39) in reduced form is similar as that obtained by Kumar et al. [25].
(iii)
In the absence of microstretch effect, that is, takingj
0
=
α
0
=
b
0
=
λ
0
=
λ
1
=
ν
1
=
ν
2
=
0 in (37), we obtained the frequency equation at the imperfect boundary between inviscid liquid layer and a micropolar thermoelastic diffusion solid half space as
(40)
tan
⁡
(
m
7
H
)
=
-
Δ
3
Δ
4
,
where the values of Δ
3 and Δ
4 are given in Appendix B and with the changed values of k
i
j in (37) as
(41)
k
i
4
=
0
,
for
i
=
1,2
,
…
,
8
,
k
5
j
=
0
,
for
j
=
1,2
,
…
,
4
,
k
6
j
=
0
,
for
j
=
1,2
,
…
,
6
.
(iv)
In the absence of microstretch and micropolarity effects, that is, takingj
0
=
α
0
=
b
0
=
λ
0
=
λ
1
=
ν
1
=
ν
2
=
0 and K
=
γ
=
0 in (37), the frequency equation at an imperfect boundary between inviscid liquid layer and thermoelastic diffusion solid half space is obtained as
(42)
tan
⁡
(
m
7
H
)
=
-
Δ
5
Δ
6
,
where the values of Δ
5 and Δ
6 are given in Appendix B and with the changed values of k
i
j in (37) as
(43)
k
i
4
=
0
,
for
i
=
1,2
,
…
,
8
,
k
5
j
=
0
,
for
j
=
1,2
,
…
,
6
,
k
6
j
=
0
,
for
j
=
1,2
,
…
,
6
.
(v)
In the absence of diffusion effect andK
n
→
∞, the frequency equation (37) reduced to the microstretch thermoelastic solid half space lying under inviscid liquid. The resulting frequency equation in reduced form is similar as that obtained by Sharma et al. [23].
(vi)
If we takeh
1
→
0 and h
2
→
0 in (37), we obtained the frequency equation for the thermally insulated and impermeable boundary with imperfect interface between inviscid liquid layer and a microstretch thermoelastic diffusion solid half space.
(vii)
If we takeh
1
→
∞ and h
2
→
∞ in (37), we obtained the frequency equation for the thermally insolated and isoconcentrated boundary with imperfect interface between inviscid liquid layer and a microstretch thermoelastic diffusion solid half space.


## 7. Numerical Results and Discussion

The analysis is conducted for a magnesium crystal-like material. Following Eringen [27], the values of micropolar parameters are
(44)
λ
=
9.4
×
10
10
N
m
-
2
,
μ
=
4.0
×
10
10
N
m
-
2
,
K
=
1.0
×
10
10
N
m
-
2
,
ρ
=
1.74
×
10
3
Kg
m
-
3
,
j
=
0.2
×
10
-
19
m
-
2
,
γ
=
0.779
×
10
-
9
N
.
Thermal and diffusion parameters are given by
(45)
C
*
=
1.04
×
10
3
JK
g
-
1
K
-
1
,
K
*
=
1.7
×
10
6
J
m
-
1
s
-
1
K
-
1
,
α
t
1
=
2.33
×
10
-
5
K
-
1
,
α
t
2
=
2.48
×
10
-
5
K
-
1
,
T
0
=
0.298
×
10
3
K
,
τ
1
=
0.01
,
τ
0
=
0.02
,
α
c
1
=
2.65
×
10
-
4
m
3
Kg
-
1
,
α
c
2
=
2.83
×
10
-
4
m
3
Kg
-
1
,
a
=
2.9
×
10
4
m
2
s
-
2
K
-
1
,
b
=
32
×
10
5
Kg
-
1
m
5
s
-
2
,
τ
1
=
0.04
,
τ
0
=
0.03
,
D
=
0.85
×
10
-
8
Kgm
-
3
s
,
H
=
1
.
The microstretch parameters are taken as
(46)
j
o
=
0.19
×
10
-
19
m
2
,
α
o
=
0.779
×
10
-
9
N
,
b
o
=
0.5
×
10
-
9
N
,
λ
o
=
0.5
×
10
10
Nm
-
2
,
λ
o
=
0.5
×
10
10
Nm
-
2
,
λ
1
=
0.5
×
10
10
Nm
-
2
.
The physical constants for water are given by ρ
l
=
1.0
×
10
3
Kg
/
m
3, λ
l
=
2.1904
×
10
3
N
/
m
2, and K
n
=
1.4
Ns
m
-
3.MATLAB software 7.04 has been used for numerical computation of the resulting quantities. The given cost for solving the dispersion relation for a given wavenumber is equivalent to solving a system of two independent quadratic equations; in MATLAB the computational efforts is negligible; the CPU cost is also independent of the input parameter set. The values of the phase velocity and attenuation coefficients with wave number are plotted graphically to depict the effect of an imperfect boundary along with the relaxation times in a microstretch thermoelastic diffusion solid half space under a homogeneous inviscid liquid layer for thermally insulated, impermeable boundaries and isothermal, isoconcentrated boundaries, respectively. Normal velocity component has also been determined in the liquid medium withH
=
1 for the above cases and is shown in Figures 2, 3, 4, 5, 6, and 7, respectively. For the thermally insulated and impermeable boundaries, the words LSH10IB and GLH10IB symbolize the graphs for L-S and G-L theories at an imperfect boundary between an inviscid liquid layer and a microstretch thermoelastic diffusion half space, whereas the words LSH10WC and GLH10WC symbolize the graphs for L-S and G-L theories at a welded contact between an inviscid liquid layer and a microstretch thermoelastic diffusion half space, respectively. Similarly, for the isothermal, isoconcentrated boundaries, the words LSH1NIB and GLH1NIB symbolize the graphs for L-S and G-L theories at an imperfect boundary between an inviscid liquid layer and a microstretch thermoelastic diffusion half space, whereas the words LSH1NWC and GLH1NWC symbolize the graphs for L-S and G-L theories at a welded contact between an inviscid liquid layer and a microstretch thermoelastic diffusion half space, respectively.Figure 2
Variation of phase velocity with respect to wave number for thermally insulated and impermeable boundaries.Figure 3
Variation of attenuation with respect to wave number for thermally insulated and impermeable boundaries.Figure 4
Variation of normal velocity component “u
3
l” in liquid half-space with respect to wave number for thermally insulated and impermeable boundaries.Figure 5
Variation of phase velocity with respect to wave number for isothermal and isoconcentrated boundaries.Figure 6
Variation of attenuation with respect to wave number for isothermal and isoconcentrated boundaries.Figure 7
Variation of normal velocity component “u
3
l” in liquid half-space with respect to wave number for isothermal and isoconcentrated boundaries.


### 7.1. Phase Velocity, Attenuation Coefficient, and the Normal Velocity Component at the Thermally Insulated and Impermeable Boundaries

Figure2 depicts the variation of phase velocity with wave number ξ. It is noticed that, for smaller values of ξ, the values of phase velocity for LSH10IB and LSH10WC decrease monotonically, whereas for higher values of ξ, the values of phase velocity decrease smoothly and finally become dispersionless; nevertheless, the trend of variation and behavior of the phase velocity for GLH10IB and GLH10WC is similar to LSH10IB and LSH10WC, respectively; however, the corresponding values are different in magnitude.Figure3 exhibits the variation of attenuation with wave number ξ. The values of attenuation for smaller values of wave number increase sharply for all the cases and become oscillatory for LSH10IB and LSH10WC in 4
≤
ξ
≤
6 and instead increase gradually for GLH10IB and GLH10WC, which becomes stationary for higher values of ξ for all the cases. A notable effect of relaxation times is also observed on the propagation of Rayleigh type surface waves at the imperfect boundary for ξ
=
4. Figure 4 exhibits the variation of normal velocity component u
3
l with wave number ξ. The trend of variation and behavior of u
3
l for all the cases is similar in the range 1
≤
ξ
≤
15; nevertheless, due to imperfect boundary, a monotonic decrease in the values of u
3
l is noticed for GLH10IB while compared with GLH10WC within the range 16
≤
ξ
≤
24.


### 7.2. Phase Velocity, Attenuation Coefficient, and the Normal Velocity Component at Isothermal and Isoconcentrated Boundaries

Figure5 depicts the variation of phase velocity with wave number ξ. It is noticed that the trend of variation and behavior of the phase velocity at the isothermal and isoconcentrated boundaries is similar to Figure 2, respectively, but their corresponding values are different in magnitude. Figure 6 exhibits the variation of attenuation with wave number ξ. The values of attenuation coefficients for smaller values of wave number increase monotonically and become oscillatory for all the cases in 4
≤
ξ
≤
6, which come to be dispersionless for higher values of ξ. Figure 7 exhibits the variation of the normal velocity component with wave number ξ. It is noticed that due to imperfect boundary the normal velocity component for LSH10WC shows a good difference in the magnitude value while compared with LSH10IB, though in both the cases the values get a sudden fall for 18
≤
ξ
≤
30. On the other hand, a worthy difference in the behavior of curves is noticed when compared with the curves for thermally insulated and impermeable boundaries, as shown in Figure 4.


## 7.1. Phase Velocity, Attenuation Coefficient, and the Normal Velocity Component at the Thermally Insulated and Impermeable Boundaries

Figure2 depicts the variation of phase velocity with wave number ξ. It is noticed that, for smaller values of ξ, the values of phase velocity for LSH10IB and LSH10WC decrease monotonically, whereas for higher values of ξ, the values of phase velocity decrease smoothly and finally become dispersionless; nevertheless, the trend of variation and behavior of the phase velocity for GLH10IB and GLH10WC is similar to LSH10IB and LSH10WC, respectively; however, the corresponding values are different in magnitude.Figure3 exhibits the variation of attenuation with wave number ξ. The values of attenuation for smaller values of wave number increase sharply for all the cases and become oscillatory for LSH10IB and LSH10WC in 4
≤
ξ
≤
6 and instead increase gradually for GLH10IB and GLH10WC, which becomes stationary for higher values of ξ for all the cases. A notable effect of relaxation times is also observed on the propagation of Rayleigh type surface waves at the imperfect boundary for ξ
=
4. Figure 4 exhibits the variation of normal velocity component u
3
l with wave number ξ. The trend of variation and behavior of u
3
l for all the cases is similar in the range 1
≤
ξ
≤
15; nevertheless, due to imperfect boundary, a monotonic decrease in the values of u
3
l is noticed for GLH10IB while compared with GLH10WC within the range 16
≤
ξ
≤
24.


## 7.2. Phase Velocity, Attenuation Coefficient, and the Normal Velocity Component at Isothermal and Isoconcentrated Boundaries

Figure5 depicts the variation of phase velocity with wave number ξ. It is noticed that the trend of variation and behavior of the phase velocity at the isothermal and isoconcentrated boundaries is similar to Figure 2, respectively, but their corresponding values are different in magnitude. Figure 6 exhibits the variation of attenuation with wave number ξ. The values of attenuation coefficients for smaller values of wave number increase monotonically and become oscillatory for all the cases in 4
≤
ξ
≤
6, which come to be dispersionless for higher values of ξ. Figure 7 exhibits the variation of the normal velocity component with wave number ξ. It is noticed that due to imperfect boundary the normal velocity component for LSH10WC shows a good difference in the magnitude value while compared with LSH10IB, though in both the cases the values get a sudden fall for 18
≤
ξ
≤
30. On the other hand, a worthy difference in the behavior of curves is noticed when compared with the curves for thermally insulated and impermeable boundaries, as shown in Figure 4.


## 8. Conclusion

The propagation of Rayleigh type surface waves at an imperfect boundary between an inviscid liquid layer and a homogeneous isotropic microstretch thermoelastic diffusion solid half-space subjected to thermally insulated/isothermal and impermeable/isoconcentrated boundaries has been investigated. A secular equation for surface wave propagation in the considered medium is derived. A significant effect of an imperfect boundary along with the relaxation times on the phase velocity, attenuation coefficient, and the normal velocity component for the considered mathematical model has been observed. Effect due to imperfect boundary along with relaxation times for LSH10IB and LSH10WC is not much significant as compared with the case of thermally insulated and impermeable boundaries. In case of GLHINIB and GLH1NWC, the attenuation coefficients have very good impact on imperfect boundary when compared to GLH10IB and GLH10WC, respectively, for the smaller value of wave number. It is perceived that the trend of variations and behavior of the all resulting quantities try to converge toward a particular value with the increase of wave number, except foru
3
l in case of GLH10IB at the insulated/isothermal boundary, which shows a significant difference due to imperfect boundary. Among all the graphs, the normal velocity components for both thermally insulated/isothermal and impermeable/isoconcentrated boundaries show respectable variance due to imperfect boundary.



---

*Source: 518384-2014-08-04.xml*