# Plane Waves and Fundamental Solutions in Heat Conducting Micropolar Fluid

Authors: Rajneesh Kumar; Mandeep Kaur
Journal: Journal of Fluids
 (2016)

---

## Abstract

In the present investigation, we study the propagation of plane waves in heat conducting micropolar fluid. The phase velocity, attenuation coefficient, specific loss, and penetration depth are computed numerically and depicted graphically. In addition, the fundamental solutions of the system of differential equations in case of steady oscillations are constructed. Some basic properties of the fundamental solution and special cases are also discussed.


---

## Body

## 1. Introduction

Eringen [1] developed the theory of microfluids, in which microfluids possess three gyration vector fields in addition to its classical translatory degrees of freedom represented by velocity field. Eringen introduced the micropolar fluids [2] which are subclass of these fluids, in which the local fluid elements possess rigid rotations without stretch. Micropolar fluids can support couple stress, the body couples, and asymmetric stress tensor and possess a rotational field, which is independent of the velocity of fluid. Anisotropic fluids, liquid crystals with rigid molecules, magnetic fluids, cloud with dust, muddy fluids, biologicaltropic fluids, and dirty fluids (dusty air and snow) over airfoil can be modeled more realistically as micropolar fluids. Ariman et al. [3, 4] studied microcontinuum fluid mechanics. Říha [5] discussed the theory of heat conducting micropolar fluid with microtemperature. Eringen and Kafadar [6] developed polar field theories. Brulin [7] discussed linear micropolar media. Flow and heat transfer in a micropolar fluid past with suction and heat sources were discussed by Agarwal and Dhanapal [8]. Payne and Straughan [9] investigated critical Rayleigh numbers for oscillatory and nonlinear convection in an isotropic thermomicropolar fluid. Gorla [10] studied combined forced and free convection in the boundary layer flow of a micropolar fluid on a continuous moving vertical cylinder. Eringen [11] investigated the theory of microstretch and bubbly liquids. Aydemir and Venart [12] investigated the flow of a thermomicropolar fluid with stretch. Yerofeyev and Soldatov [13] discussed a shear surface wave at the interface of an elastic body and a micropolar liquid. The theory of elastic and viscoelastic micropolar liquids was studied by Yeremeyev and Zubov [14]. Hsia and Cheng [15] discussed longitudinal plane waves propagation in elastic micropolar porous media. Hsia et al. [16] studied propagation of transverse waves in elastic micropolar porous semispaces.Construction of fundamental solution of systems of partial differential equations is necessary to investigate the boundary value problems of the theory of elasticity and thermoelasticity. The fundamental solutions in the classical theory of coupled thermoelasticity were firstly studied by Hetnarski [17, 18]. Hetnarski and Ignaczak [19] studied generalized thermoelasticity. Svanadze [20–25] constructed the fundamental solutions in the microcontinuum field theories. Kumar and Kansal [26] investigated the fundamental solution in the theory of thermomicrostretch elastic diffusive solids. Fundamental solution in the theory of micropolar thermoelastic diffusion with voids was studied by Kumar and Kansal [27]. Recently, Kumar and Kansal [28] discussed plane waves and fundamental solution in the generalized theories of thermoelastic diffusion. Kumar and Kansal [29] studied propagation of plane waves and fundamental solution in the theories of thermoelastic diffusive materials with voids. The information related to fundamental solutions of differential equations is contained in the books of Hörmander [30, 31].The main objective of the present paper is to study the propagation of plane waves in heat conducting micropolar fluid. Several qualitative characterizations of the wave field, such as phase velocity, attenuation coefficient, specific loss, and penetration depth, are computed and depicted graphically for different values of frequency. The representation of fundamental solution of system of equations in the case of steady oscillations is obtained in terms of elementary functions. Some particular cases have also been deduced.


## 2. Basic Equations

In three-dimensional spaceE
3, let x
=
(
x
1
,
x
2
,
x
3
) be the points of the Euclidean space, t represents the time variable, and D
x
=
(
∂
/
∂
x
1
,
∂
/
∂
x
2
,
∂
/
∂
x
3
).Following Ciarletta [32], the basic equations for homogeneous, isotropic heat conducting micropolar fluids without body forces, body couples, and heat sources are given by (1)
D
1
v
+
λ
+
μ
∇
∇
·
v
+
K
∇
×
Ψ
-
b
∇
T
-
c
0
∇
ϕ
∗
=
0
,
D
2
Ψ
+
α
+
β
∇
∇
·
Ψ
+
K
∇
×
v
=
0
,
K
∗
∇
2
T
-
b
T
0
∇
·
v
=
ρ
0
a
T
0
∂
T
∂
t
,
ρ
0
∂
ϕ
∗
∂
t
=
∇
·
v
,where(2)
D
1
=
μ
+
K
Δ
-
ρ
0
∂
∂
t
,
D
2
=
γ
Δ
-
I
∂
∂
t
-
2
K
,where λ, μ, K, α, β, γ, and c
0 are material constants of the fluid. v and Ψ are the velocity vector and microrotation velocity vector, ρ
0 is the density, I is a scalar constant with the dimension of moment of inertia of unit mass, K
∗ is the thermal conductivity, a
T
0 is the specific heat at constant strain, T
0 is the absolute temperature, T is the temperature change, ϕ
∗ is the variation in specific volume, b
=
3
λ
+
2
μ
+
K
α
T, where α
T is the coefficient of linear thermal expansion, and Δ is the Laplacian operator.For convenience, the following nondimensional quantities are introduced:(3)
x
i
′
=
ω
∗
c
1
2
x
i
,
v
i
′
=
1
ω
∗
c
1
2
v
i
,
ψ
2
′
=
ρ
0
c
1
2
b
T
0
ψ
2
,
t
′
=
ω
∗
t
,
ϕ
∗
=
ρ
0
ϕ
∗
,
T
′
=
T
T
0
,where ω
∗
=
λ
/
I
,
c
1
2
=
λ
+
2
μ
+
K
/
ρ
0
,
a
n
d
ω
∗ is the characteristic frequency of the medium.Making use of (2) and (3) in basic equations (1) and after suppressing the primes, we obtain(4)
δ
1
Δ
v
+
g
r
a
d
d
i
v
v
+
δ
2
c
u
r
l
Ψ
-
δ
3
∇
T
-
δ
4
∇
ϕ
∗
=
δ
5
v
.
,
(5)
δ
7
g
r
a
d
d
i
v
Ψ
+
Δ
Ψ
-
2
δ
6
Ψ
+
δ
8
c
u
r
l
v
=
δ
9
Ψ
.
,
(6)
Δ
T
=
δ
11
T
+
δ
10
∇
·
v
.
,
(7)
∂
∂
t
ϕ
∗
=
∇
·
v
.Making use of (7) in (4), we obtain(8)
δ
1
Δ
v
.
+
g
r
a
d
d
i
v
v
.
+
δ
2
c
u
r
l
Ψ
.
-
δ
3
∇
T
.
-
δ
4
∇
∇
·
v
=
δ
5
v
¨
,where(9)
δ
1
=
μ
+
K
λ
+
μ
,
δ
2
=
K
b
T
0
λ
+
μ
ω
∗
ρ
0
c
1
2
,
δ
3
=
b
T
0
λ
+
μ
ω
∗
,
δ
4
=
c
0
λ
+
μ
ρ
0
ω
∗
,
δ
5
=
ρ
0
c
1
2
λ
+
μ
,
δ
6
=
K
c
1
2
γ
ω
∗
,
δ
7
=
α
+
β
γ
,
δ
8
=
K
ρ
0
c
1
4
γ
b
T
0
,
δ
9
=
I
c
1
2
γ
,
δ
10
=
b
c
1
2
K
∗
,
δ
11
=
ρ
0
a
T
0
c
1
2
K
∗
.For two-dimensional problem, we take(10)
v
=
v
1
x
1
,
x
3
,
0
,
v
3
x
1
,
x
3
,
Ψ
=
0
,
Ψ
2
x
1
,
x
3
,
0
.The relation between dimensionless velocity components v
1 and v
3 and nondimensional velocity potential functions ϕ and ψ is expressed as(11)
u
1
=
∂
ϕ
∂
x
1
-
∂
ψ
∂
x
3
,
u
3
=
∂
ϕ
∂
x
3
+
∂
ψ
∂
x
1
.Making use of (10)-(11) in (5), (6), and (8), we obtain(12)
δ
1
+
1
∂
∂
t
-
δ
4
∇
2
ϕ
-
δ
3
∂
∂
t
T
-
δ
5
∂
2
ϕ
∂
t
2
=
0
,
δ
1
∂
∂
t
∇
2
-
δ
5
∂
2
∂
t
2
ψ
+
δ
2
∂
Ψ
2
∂
t
=
0
,
∇
2
Ψ
2
-
δ
8
∇
2
ψ
-
2
δ
6
Ψ
2
-
δ
9
∂
Ψ
2
∂
t
=
0
,
∇
2
T
-
δ
11
∂
∂
t
T
-
δ
10
∇
2
ϕ
=
0
,where ∇
2
=
∂
2
/
∂
x
1
2
+
∂
2
/
∂
x
3
2.


## 3. Solution of Plane Waves

For plane harmonic waves, we assume the solution of the form(13)
ϕ
,
ψ
,
Ψ
2
,
T
=
ϕ
¯
,
ψ
¯
,
Ψ
2
¯
,
T
¯
e
x
p
ι
ξ
x
1
l
1
+
x
3
l
3
-
ω
t
,where ω is the circular frequency and ξ is the complex wave number. ϕ
¯
,
ψ
¯
,
Ψ
2
¯
,
T
¯ are undetermined amplitude vectors that are independent of time t and coordinates x
m
(
m
=
1,3
). l
1 and l
3 are the direction cosines of the wave normal onto x
1
x
3-plane with the property l
1
2
+
l
3
2
=
1.Using (13) in (12), we obtain(14)
ξ
2
δ
1
∗
ι
ω
+
δ
4
+
δ
5
ω
2
ϕ
¯
+
δ
3
ι
ω
T
¯
=
0
,
ι
ω
δ
1
ξ
2
+
δ
5
ω
2
ψ
¯
+
δ
2
ι
ω
Ψ
2
¯
=
0
,
-
ξ
2
-
2
δ
6
+
δ
9
ι
ω
Ψ
2
¯
+
δ
8
ξ
2
ψ
¯
=
0
,
-
ξ
2
+
δ
11
ι
ω
T
¯
+
δ
10
ξ
2
ϕ
¯
=
0
.The system of (14) will have nontrivial solution, if the determinant of the coefficients ϕ
¯
,
ψ
¯
,
Ψ
2
¯
,
T
¯ vanishes which on expansion yield(15)
G
1
ξ
4
+
G
2
ξ
2
+
G
3
=
0
,
G
4
ξ
4
+
G
5
ξ
2
+
G
5
=
0
,where(16)
G
1
=
δ
1
∗
ι
ω
+
δ
4
,
G
2
=
-
δ
11
ι
ω
δ
1
∗
ι
ω
+
δ
4
+
δ
5
ω
2
+
δ
3
δ
10
ι
ω
,
G
3
=
-
δ
5
δ
11
ι
ω
3
,
G
4
=
-
δ
1
ι
ω
,
G
5
=
-
δ
5
ω
2
+
δ
1
ι
ω
δ
9
ι
ω
-
2
δ
6
-
δ
2
δ
8
ι
ω
,
G
6
=
δ
5
ω
2
δ
9
ι
ω
-
2
δ
6
.Solving (15), we obtain eight roots of ξ, in which four roots of ξ, that is, ξ
1
,
ξ
2
,
ξ
3, and ξ
4, correspond to positive x
3-direction and other four roots of ξ, that is, -
ξ
1
,
-
ξ
2
,
-
ξ
3, and -
ξ
4, correspond to negative x
3-direction. Now and after, we will restrict our work to positive x
3-direction. Corresponding to roots ξ
1
,
ξ
2
,
ξ
3, and ξ
4 there exist four waves in descending order of their velocities, that is, two coupled longitudinal waves and two coupled transverse waves.The expressions for phase velocity, attenuation coefficient, specific loss, and penetration depth of above waves are derived as follows.(i) Phase Velocity. The phase velocities are given by(17)
V
i
=
ω
R
e
ξ
i
;
i
=
1,2
,
3,4
,where V
1
,
V
2
,
V
3
,
V
4 are the phase velocities of two coupled longitudinal waves and two coupled transverse waves, respectively.(ii) Attenuation Coefficient. The attenuation coefficients are defined as(18)
Q
i
=
I
m
ξ
i
;
i
=
1,2
,
3,4
,where Q
1
,
Q
2
,
Q
3
,
Q
4 are the attenuation coefficients of two coupled longitudinal waves and two coupled transverse waves, respectively.(iii) Specific Loss. The specific loss is the ratio of energy (W
¯) dissipated in taking a specimen through a rate of stress cycle, to the elastic energy (W) stored in the specimen when the rate of 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 [33] shows that the specific loss (
W
¯
/
W
) equals 4π times the absolute value of the imaginary part of ξ to the real part of ξ; that is,(19)
S
i
=
W
¯
W
i
=
4
π
I
m
ξ
i
R
e
ξ
i
;
i
=
1,2
,
3,4
,where S
1
,
S
2
,
S
3
,
S
4 are the specific loss of two coupled longitudinal waves and two coupled transverse waves, respectively.(iv) Penetration Depth. The penetration depths are defined by(20)
P
i
=
1
I
m
ξ
i
;
i
=
1,2
,
3,4
,where P
1
,
P
2
,
P
3
,
P
4 are the penetration depths of two coupled longitudinal waves and two coupled transverse waves, respectively.


## 4. Steady Oscillations

Let us assume the solution of the form(21)
v
¯
x
,
t
,
Ψ
¯
x
,
t
,
T
f
¯
x
,
t
=
R
e
v
,
Ψ
,
T
f
e
-
ι
ω
t
,where ω
>
0 is the frequency of oscillation.Making use of (21) into (5), (6), and (8), the system of equations of steady oscillations is obtained as(22)
δ
1
Δ
+
ι
δ
5
ω
v
+
δ
∗
g
r
a
d
d
i
v
v
+
δ
2
c
u
r
l
Ψ
+
δ
3
∇
T
f
=
0
,
Δ
+
μ
∗
Ψ
+
δ
7
g
r
a
d
d
i
v
Ψ
+
δ
8
c
u
r
l
v
=
0
,
Δ
+
δ
11
ι
ω
T
f
-
δ
10
d
i
v
v
=
0
,where μ
∗
=
δ
9
ι
ω
-
2
δ
6
,
δ
∗
=
1
-
ι
δ
4
/
ω.The matrix differential operator is taken as(23)
F
D
x
=
F
g
h
D
x
7
×
7
,where(24)
f
m
n
D
x
=
δ
1
Δ
+
ι
δ
5
ω
δ
m
n
+
δ
∗
∂
2
∂
x
m
∂
x
n
,
F
m
+
3
,
n
D
x
=
δ
2
∑
r
=
1
3
ε
m
r
n
∂
∂
x
r
,
F
m
7
D
x
=
δ
3
∂
∂
x
m
,
F
m
+
3
,
n
+
3
D
x
=
Δ
+
μ
∗
δ
m
n
+
δ
7
∂
2
∂
x
m
∂
x
n
,
F
m
+
3
,
n
D
x
=
-
δ
8
∑
r
=
1
3
ε
m
r
n
∂
∂
x
r
,
F
m
+
3,7
D
x
=
F
7
,
n
+
3
=
0
,
F
77
D
x
=
Δ
+
δ
11
ι
ω
,
F
7
n
D
x
=
-
δ
10
∂
∂
x
n
,
m
,
n
=
1,2
,
3
.The above system of (22) can be represented in the following form:(25)
F
D
x
V
x
=
0
,where V
=
(
v
,
Ψ
,
T
f
) is a seven component vector function on E
3.Let us assume that(26)
δ
1
≠
0
.
F is an elliptic differential operator Hörmander [30], if condition (26) is fulfilled.Definition 1.
The fundamental solution of the system of (21)-(22) (the fundamental matrix of operator F) is the matrix G
(
x
)
=
G
g
h
x
7
×
7, satisfying condition [30](27)
F
D
x
G
x
=
δ
x
I
x
.Here δ is the Dirac delta, I
=
δ
g
h
7
×
7 is the unit matrix, and x
∈
E
3.Now furtherG
(
x
) in terms of elementary functions is constructed.


## 5. Fundamental Solution of System of Equation of Steady Oscillations

We consider the following system of equations:(28)
δ
1
Δ
v
+
δ
∗
g
r
a
d
d
i
v
v
+
δ
8
c
u
r
l
Ψ
-
δ
10
∇
T
f
+
ι
δ
5
ω
v
=
H
′
,
(29)
Δ
+
μ
∗
Ψ
+
δ
7
g
r
a
d
d
i
v
Ψ
+
δ
2
c
u
r
l
v
=
H
′
′
,
(30)
δ
3
d
i
v
v
+
Δ
+
δ
11
ι
ω
T
f
=
L
,where H
′
,
H
′
′ are three-component vector function on E
3 and L is scalar function on E
3.The system of (28)–(30) can be written as(31)
F
t
r
D
x
V
x
=
Q
x
,where F
t
r is the transpose of matrix F, Q
=
(
H
′
,
H
′
′
,
L
), and x
∈
E
3.The following equations are obtained by applying the operatord
i
v to (28) and (29):(32a)
Δ
δ
1
∗
+
δ
5
ι
ω
d
i
v
v
-
δ
10
Δ
T
f
=
d
i
v
H
′
,
(32b)
Δ
δ
7
∗
+
μ
∗
d
i
v
Ψ
=
d
i
v
H
′
′
,
(32c)
δ
3
d
i
v
v
+
Δ
+
δ
11
ι
ω
T
f
=
L
,where δ
1
+
1
=
δ
1
∗
,
δ
7
+
1
=
δ
7
∗.Equations (32a) and (32c) can be written as(33)
Γ
1
Δ
d
i
v
v
=
Φ
1
,
(34)
Γ
1
Δ
T
f
=
Φ
2
,where(35)
Γ
1
Δ
=
e
∗
d
e
t
Δ
δ
1
∗
+
δ
5
ι
ω
-
δ
10
Δ
δ
3
Δ
+
δ
11
ι
ω
,
e
∗
=
1
δ
1
∗
,
(36)
Φ
1
=
e
∗
Δ
+
δ
11
ι
ω
d
i
v
H
′
+
δ
10
Δ
L
,
(37)
Φ
2
=
e
∗
-
δ
3
d
i
v
H
′
+
Δ
δ
1
∗
+
ι
ω
δ
5
L
.It can be seen that(38)
Γ
1
Δ
=
∏
m
=
1
2
Δ
+
λ
m
2
,and λ
m
2
,
m
=
1,2
,
3 are the roots of the equation Γ
1
(
-
k
)
=
0 (with respect to k).From (32b), it is seen that(39)
Δ
+
λ
5
2
d
i
v
Ψ
=
1
δ
7
∗
d
i
v
H
′
′
,where λ
5
2
=
μ
∗
/
δ
7
∗.Applying the operatorsΔ
+
μ
∗ and δ
8
c
u
r
l to (28) and (30), respectively, we obtain(40)
Δ
+
μ
∗
δ
1
Δ
v
+
δ
∗
g
r
a
d
d
i
v
v
+
δ
5
ι
ω
2
v
+
δ
8
Δ
+
μ
∗
c
u
r
l
Ψ
=
Δ
+
μ
∗
H
′
+
δ
10
g
r
a
d
T
f
,
(41)
δ
8
Δ
+
μ
∗
c
u
r
l
Ψ
=
-
δ
2
δ
8
c
u
r
l
c
u
r
l
v
+
δ
8
c
u
r
l
H
′
′
.Now(42)
c
u
r
l
c
u
r
l
v
=
g
r
a
d
d
i
v
v
-
Δ
v
.Using (41) and (42) in (40), we obtain(43)
Δ
+
μ
∗
δ
1
+
δ
8
Δ
+
δ
5
ι
ω
Δ
+
μ
∗
v
=
-
Δ
+
μ
∗
δ
∗
-
δ
2
δ
8
g
r
a
d
d
i
v
v
+
Δ
+
μ
∗
H
′
+
δ
10
g
r
a
d
T
f
-
δ
8
c
u
r
l
H
′
′
.Applying the operator Γ
1
(
Δ
) to (43) and using (33)-(34), we obtain(44)
Γ
1
Δ
Δ
+
μ
∗
δ
1
+
δ
8
Δ
+
δ
5
ι
ω
Δ
+
μ
∗
v
=
-
Δ
+
μ
∗
δ
∗
-
δ
2
δ
8
g
r
a
d
Φ
1
+
Δ
+
μ
∗
Γ
1
Δ
H
′
+
δ
10
g
r
a
d
T
f
-
Γ
1
Δ
δ
8
c
u
r
l
H
′
′
.The above equation can be rewritten as(45)
Γ
1
Δ
Γ
2
Δ
v
=
Φ
′
,where(46)
Γ
2
Δ
=
f
∗
d
e
t
δ
1
Δ
+
δ
5
ι
ω
δ
8
Δ
-
δ
2
Δ
+
μ
∗
,
f
∗
=
1
δ
1
,
(47)
Φ
′
=
f
∗
-
Δ
+
μ
∗
δ
∗
-
δ
8
g
r
a
d
Φ
1
+
Δ
+
μ
∗
Γ
1
Δ
H
′
+
δ
10
g
r
a
d
Φ
2
-
δ
8
Γ
1
Δ
c
u
r
l
H
′
′
.It can be seen that(48)
Γ
2
Δ
=
Δ
+
λ
3
2
Δ
+
λ
4
2
,where λ
3
2
,
λ
4
2 are the roots of the equation Γ
2
(
-
k
)
=
0 (with respect to k).Applying the operatorsδ
2
c
u
r
l and (
δ
1
Δ
+
δ
5
ι
ω
) to (29) and (30), respectively, we obtain(49)
δ
1
Δ
+
δ
5
ι
ω
δ
2
c
u
r
l
v
=
δ
2
c
u
r
l
H
′
-
δ
2
δ
8
c
u
r
l
c
u
r
l
Ψ
,
(50)
Δ
+
μ
∗
δ
1
Δ
+
δ
5
ι
ω
Ψ
+
δ
7
δ
1
Δ
+
δ
5
ι
ω
g
r
a
d
d
i
v
Ψ
+
δ
2
δ
1
Δ
+
δ
5
ι
ω
c
u
r
l
v
=
δ
1
Δ
+
δ
5
ι
ω
H
′
′
.Now(51)
c
u
r
l
c
u
r
l
Ψ
=
g
r
a
d
d
i
v
Ψ
-
Δ
Ψ
,Using (47) and (51) in (50), we obtain(52)
Δ
+
μ
∗
δ
1
Δ
+
δ
5
ι
ω
Ψ
+
δ
7
δ
1
Δ
+
δ
5
ι
ω
g
r
a
d
d
i
v
Ψ
+
δ
2
c
u
r
l
H
′
-
δ
2
δ
8
g
r
a
d
d
i
v
Ψ
-
Δ
Ψ
=
δ
1
Δ
+
δ
5
ι
ω
H
′
′
.The above equation may also be written as(53)
Δ
+
μ
∗
δ
1
Δ
+
δ
5
ι
ω
+
δ
2
δ
8
Δ
Ψ
=
-
δ
7
δ
1
Δ
+
δ
5
ι
ω
-
δ
2
δ
8
g
r
a
d
d
i
v
Ψ
-
δ
2
c
u
r
l
H
′
+
δ
1
Δ
+
δ
5
ι
ω
H
′
′
.Applying operator (
Δ
+
λ
5
2
) to (53) and using (39), we obtain(54)
Δ
+
λ
5
2
δ
1
Δ
2
+
Δ
δ
1
μ
∗
+
δ
5
ι
ω
+
δ
2
δ
8
+
δ
5
ι
ω
μ
∗
Ψ
=
-
δ
7
δ
1
Δ
+
δ
5
ι
ω
-
δ
2
δ
8
g
r
a
d
1
δ
7
∗
d
i
v
H
′
′
-
δ
2
Δ
+
λ
5
2
c
u
r
l
H
′
+
Δ
+
λ
5
2
δ
1
Δ
+
δ
5
ι
ω
2
H
′
′
.The above equation can also be written as(55)
Γ
2
Δ
Δ
+
λ
5
2
Ψ
=
Φ
′
′
,where(56)
Φ
′
′
=
f
∗
-
δ
7
δ
1
Δ
+
δ
5
ι
ω
-
δ
2
δ
8
1
δ
7
∗
g
r
a
d
d
i
v
H
′
′
-
δ
2
Δ
+
λ
5
2
c
u
r
l
H
′
+
Δ
+
λ
5
2
δ
1
Δ
+
δ
5
ι
ω
H
′
′
.From (34), (45), and (55), we obtain(57)
Θ
Δ
V
x
=
Φ
^
x
,where(58)
Φ
^
=
Φ
′
,
Φ
′
′
,
Φ
2
,
Θ
Δ
=
Θ
g
h
Δ
7
×
7
,
Θ
m
m
Δ
=
Γ
1
Δ
Γ
2
Δ
=
∏
q
=
1
4
Δ
+
λ
q
2
,
Θ
m
+
3
,
n
+
3
Δ
=
Γ
2
Δ
Δ
+
λ
5
2
,
Θ
g
h
Δ
=
0
,
Θ
77
Δ
=
Γ
1
Δ
,
m
,
n
=
1,2
,
3
,
g
,
h
=
1,2
,
3
,
…
,
7
,
g
≠
h
.Equations (37), (47), and (56) can be rewritten in the form(59)
Φ
′
=
f
∗
Δ
+
μ
∗
Γ
1
Δ
J
+
q
11
Δ
g
r
a
d
d
i
v
H
′
+
q
21
Δ
c
u
r
l
H
′
′
+
q
31
Δ
g
r
a
d
L
,
(60)
Φ
′
′
=
q
12
Δ
c
u
r
l
H
′
+
f
∗
Δ
+
λ
5
2
δ
1
Δ
+
δ
5
ι
ω
J
+
q
22
Δ
g
r
a
d
d
i
v
H
′
′
,
(61)
Φ
2
=
q
13
Δ
d
i
v
H
′
+
q
33
Δ
L
,where J
=
δ
g
h
3
×
3 is the unit matrix.In (59)–(61), we have used the following notations:(62)
q
11
Δ
=
e
∗
f
∗
-
Δ
+
μ
∗
δ
∗
-
δ
8
Δ
+
δ
11
ι
ω
-
Δ
+
μ
∗
δ
10
δ
3
,
q
21
Δ
=
-
f
∗
δ
8
Γ
1
Δ
,
q
31
Δ
=
e
∗
f
∗
δ
8
-
Δ
+
μ
∗
δ
∗
δ
10
+
δ
10
Δ
+
μ
∗
Δ
δ
1
∗
+
δ
5
ι
ω
,
q
12
Δ
=
-
f
∗
δ
2
Δ
+
λ
5
2
,
q
22
Δ
=
-
f
∗
δ
7
∗
δ
7
δ
1
Δ
+
δ
5
ι
ω
-
δ
2
δ
8
,
q
13
Δ
=
-
e
∗
δ
3
,
q
33
Δ
=
e
∗
Δ
δ
1
∗
+
δ
5
ι
ω
.Now from (59)–(61), we have(63)
Φ
^
x
=
R
t
r
D
x
Q
x
,where(64)
R
=
R
m
n
7
×
7
,
R
m
n
D
x
=
f
∗
Δ
+
μ
∗
Γ
1
Δ
δ
m
n
+
q
11
Δ
∂
2
∂
x
m
∂
x
n
,
R
m
,
n
+
3
D
x
=
q
12
Δ
∑
r
=
1
3
ε
m
r
n
∂
∂
x
r
,
R
m
7
D
x
=
q
13
Δ
∂
∂
x
m
,
R
m
+
3
,
n
D
x
=
q
21
Δ
∑
r
=
1
3
ε
m
r
n
∂
∂
x
r
,
R
m
+
3
,
n
+
3
D
x
=
f
∗
Δ
+
λ
5
2
δ
1
Δ
+
δ
5
ι
ω
δ
m
n
+
q
22
Δ
∂
2
∂
x
m
∂
x
n
,
R
m
+
3,7
D
x
=
R
7
,
m
+
3
D
x
=
0
,
R
7
n
D
x
=
q
31
Δ
∂
∂
x
n
,
R
77
D
x
=
q
33
Δ
,
m
,
n
=
1,2
,
3
.The following relation is obtained from (31), (57), and (61):(65)
Θ
V
=
R
t
r
F
t
r
V
.The above relation implies that(66)
R
t
r
F
t
r
=
Θ
,
F
D
x
R
D
x
=
Θ
Δ
.Let us assume that(67)
λ
m
2
≠
λ
n
2
≠
0
,
m
,
n
=
1,2
,
3,4
,
5
,
m
≠
n
.Let(68)
Y
x
=
Y
r
s
x
7
×
7
,
Y
m
m
x
=
∑
n
=
1
4
r
1
n
ς
n
x
,
Y
m
+
3
,
m
+
3
x
=
∑
n
=
3
5
r
2
n
ς
n
x
,
Y
77
x
=
∑
n
=
1
2
r
3
n
ς
n
x
,
Y
u
v
x
=
0
,
m
=
1,2
,
3
,
v
,
w
=
1,2
,
…
,
7
,
v
≠
w
,where(69)
ς
n
x
=
-
1
4
π
x
e
x
p
ι
λ
n
x
,
n
=
1,2
,
…
,
5
,
r
1
l
=
∏
m
=
1
,
m
≠
l
4
λ
m
2
-
λ
l
2
-
1
,
l
=
1,2
,
3,4
,
r
2
v
=
∏
m
=
3
,
m
≠
v
5
λ
m
2
-
λ
v
2
-
1
,
v
=
3,4
,
5
,
r
3
g
=
∏
m
=
1
,
m
≠
g
2
λ
m
2
-
λ
g
2
-
1
,
g
=
1,2
.Now the following lemma will be proved.Lemma 2.
The matrixY defined above is the fundamental matrix of operator Θ
(
Δ
); that is,(70)
Θ
Δ
Y
x
=
δ
x
I
x
.Proof.
To prove the lemma, it is proved that(71a)
Γ
1
Δ
Γ
2
Δ
Y
11
x
=
δ
x
,
(71b)
Γ
2
Δ
Δ
+
λ
5
2
Y
44
x
=
δ
x
,
(71c)
Γ
1
Δ
Y
77
x
=
δ
x
.We find that(72)
r
11
+
r
12
+
r
13
+
r
14
=
0
,
∑
j
=
2
4
r
1
j
λ
1
2
-
λ
j
2
=
0
,
∑
j
=
1
2
r
1
j
∏
m
=
3
4
λ
m
2
-
λ
j
2
=
0
,
∑
j
=
3
4
r
1
j
∏
m
=
1
3
λ
m
2
-
λ
j
2
=
0
,
r
14
λ
1
2
-
λ
4
2
λ
2
2
-
λ
4
2
λ
3
2
-
λ
4
2
=
1
,
Δ
+
λ
m
2
ς
n
x
=
δ
x
+
λ
m
2
-
λ
n
2
ς
n
x
,
m
,
n
=
1,2
,
3,4
.Now consider(73)
Γ
1
Δ
Γ
2
Δ
Y
11
x
=
Δ
+
λ
2
2
Δ
+
λ
3
2
Δ
+
λ
4
2
∑
n
=
1
4
r
1
n
δ
x
+
λ
1
2
-
λ
n
2
ς
n
x
=
Δ
+
λ
2
2
Δ
+
λ
3
2
Δ
+
λ
4
2
∑
n
=
2
4
r
1
n
λ
1
2
-
λ
n
2
ς
n
x
=
Δ
+
λ
3
2
Δ
+
λ
4
2
∑
n
=
2
4
r
1
n
λ
1
2
-
λ
n
2
δ
x
+
λ
2
2
-
λ
n
2
ς
n
x
=
Δ
+
λ
3
2
Δ
+
λ
4
2
∑
n
=
3
4
r
1
n
λ
1
2
-
λ
n
2
λ
2
2
-
λ
n
2
ς
n
x
=
Δ
+
λ
4
2
∑
n
=
2
4
r
1
n
λ
1
2
-
λ
n
2
λ
2
2
-
λ
n
2
δ
x
+
λ
3
2
-
λ
n
2
ς
n
x
=
Δ
+
λ
4
2
ς
4
x
=
δ
x
.In the similar way, (71b) and (71c) can be proved.
The following matrix is now introduced:(74)
G
x
=
R
D
x
Y
x
.From (66), (70), and (74), it is obtained that(75)
F
D
x
G
x
=
F
D
x
R
D
x
Y
x
=
δ
x
I
x
.Hence, G
(
x
) is a solution to (27).
Hence, the following theorem has been proved.Theorem 3.
The matrixG
(
x
) defined by (60) is the fundamental solution of system of (22).


## 6. Basic Properties of the Matrix G(x)

Property 1.
Each column of the matrixG
(
x
) is the solution of the system of (22) at every point x
∈
E
3 except the origin.Property 2.
The matrixG
(
x
) can be written as(76)
G
=
G
m
n
7
×
7
,
G
m
n
x
=
R
m
n
D
x
Y
11
x
,
G
m
,
n
+
3
x
=
R
m
,
n
+
3
D
x
Y
44
x
,
G
m
p
x
=
R
m
p
D
x
Y
77
x
,
m
=
1,2
,
…
,
7
,
n
=
1,2
,
3
,
p
=
7
.


## 7. Numerical Results and Discussion

The following values of relevant parameters for numerical computations are taken.Following Singh and Tomar [34], the values of micropolar constants are taken as (77)
λ
=
0.15
×
1
0
8
N
s
e
c
m
-
2
,
μ
=
0.03
×
1
0
8
N
s
e
c
m
-
2
,
K
=
0.2
×
1
0
5
N
s
e
c
m
-
2
,
γ
=
0.0222
×
1
0
5
N
s
e
c
,
ρ
0
=
0.8
×
1
0
3
k
g
m
-
3
,
I
=
0.00400
×
1
0
-
16
N
s
e
c
2
m
-
2
.Thermal parameters are taken as of comparable magnitude: (78)
T
0
=
0.196
K
,
K
∗
=
0.89
×
1
0
2
N
s
e
c
-
1
K
-
1
,
c
0
=
0.005
×
1
0
11
N
2
s
e
c
2
m
-
6
,
a
=
1.5
×
1
0
5
m
2
s
e
c
-
2
K
-
2
.The variations of phase velocities, attenuation coefficients, with respect to frequency have been shown in Figures 1–4 and 5–8, respectively. In Figures 1–8, solid line corresponds to heat conducting micropolar fluid for K
=
0.2 (MT1) and dash line corresponds to heat conducting micropolar fluid for K
=
0.4 (MT2).Figure 1
Variation of phase velocity with frequencyω.Figure 2
Variation of phase velocity with frequencyω.Figure 3
Variation of phase velocity with frequencyω.Figure 4
Variation of phase velocity with frequencyω.Figure 5
Variation of attenuation coefficient with frequencyω.Figure 6
Variation of attenuation coefficient with frequencyω.Figure 7
Variation of attenuation coefficient with frequencyω.Figure 8
Variation of attenuation coefficient with frequencyω.


### 7.1. Phase Velocity

It is noticed from Figures1–4 that the magnitudes of the phase velocities V
i
,
i
=
1,2
,
3,4, for MT1 and MT2 increase with increase in frequency. The phase velocity V
1 for MT1 is greater than the phase velocity for MT2, while for V
2
,
V
3
,
V
4 the behavior is reversed. This shows that as the value of micropolar constant increases, the phase velocity V
1 decreases, while other phase velocities increase. There is slight difference in the phase velocity V
2 for MT1 and MT2.


### 7.2. Attenuation Coefficient

Figures5–8 depict that the magnitudes of attenuation coefficient Q
i
,
i
=
1,2
,
3,4, get increased with increase in wave number ω. It is noticed that, with increase in the value of K, the magnitude of attenuation coefficient increases; that is, the values of attenuation coefficients for MT2 are greater than the values for MT1 that shows the effect of micropolarity.


## 7.1. Phase Velocity

It is noticed from Figures1–4 that the magnitudes of the phase velocities V
i
,
i
=
1,2
,
3,4, for MT1 and MT2 increase with increase in frequency. The phase velocity V
1 for MT1 is greater than the phase velocity for MT2, while for V
2
,
V
3
,
V
4 the behavior is reversed. This shows that as the value of micropolar constant increases, the phase velocity V
1 decreases, while other phase velocities increase. There is slight difference in the phase velocity V
2 for MT1 and MT2.


## 7.2. Attenuation Coefficient

Figures5–8 depict that the magnitudes of attenuation coefficient Q
i
,
i
=
1,2
,
3,4, get increased with increase in wave number ω. It is noticed that, with increase in the value of K, the magnitude of attenuation coefficient increases; that is, the values of attenuation coefficients for MT2 are greater than the values for MT1 that shows the effect of micropolarity.


## 8. Conclusion

In the present paper, we have studied the propagation of plane waves in heat conducting micropolar fluid. The magnitudes of phase velocities and attenuation coefficients are depicted numerically and presented graphically with respect to frequency. Appreciable micropolarity effect is observed on these amplitudes. It is noticed that the magnitudes of phase velocityV
i
,
i
=
2,3
,
4, and attenuation coefficient Q
i
,
i
=
1,2
,
3,4, for K
=
0.4 remain more than the magnitude for K
=
0.2. This reveals that the more the value of micropolar constant, the more the magnitude of phase velocity and attenuation coefficient. The fundamental solution in heat conducting micropolar fluid in case of steady oscillations in terms of elementary functions is also constructed in the present study.



---

*Source: 1453613-2016-07-13.xml*