Ann.
Acad.
Rom.
Sci.
Ser.
Math.
Appl.
ISSN
2066-6594
Vol.
18,
No.
4/2026
STRUCTURAL
STABILITY
OF
MICROPOLAR
FLUID
FLOW
IN
A
CURVED
PIPE
WITH
DARCY–LAPWOOD–FORCHHEIMER–
BRINKMAN
MODEL
∗
A.
Prabal
†
M.
Devakar
‡
DOI 10.56082/annalsarscimath.2026.4.385
Abstract
In
this
article,
we
address
the
structural
stability
of
the
Darcy–
Lapwood–Forchheimer–Brinkman
model
for
a
micropolar
fluid
flowing
through
a
curved
pipe.
A
rigorous
analysis
is
carried
out
to
prove
the
continuous
dependence
of
solutions
of
the
model
on
several
physical
and
geometric
parameters
such
as
the
Reynolds
number,
curvature
ratio,
Forchheimer
coefficient,
and
micropolarity
parameter.
Keywords:
Darcy–Lapwood–Forchheimer–Brinkman
model,
micropo-
lar
fluid,
curved
pipe,
structural
stability,
continuous
dependence.
MSC:
35Q35,
76A05,
76S05,
35B30.
1
Introduction
and
mathematical
model
A
central
question
in
mathematical
modelling
is
whether
solutions
depend
continuously
on
variations
in
coefficients,
boundary
conditions,
or
the
gov-
∗
Accepted
for
publication
on
October
20,
2026
†
ds20mth001@students.vnit.ac.in
,
Department
of
Mathematics,
Visvesvaraya
Na-
tional
Institute
of
Technology
Nagpur,
440010
India
‡
m
devakar@mth.vnit.ac.in
,
Department
of
Mathematics,
Visvesvaraya
National
In-
stitute
of
Technology
Nagpur,
440010
India
385
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
386
erning
partial
differential
equations
themselves.
This
continuous
dependence
of
the
solutions
with
respect
to
the
model’s
structural
components
is
known
as
structural
stability.
For
real
life
problems,
variations
in
coefficients
of
the
governing
partial
differential
equations
correspond
to
changes
in
the
con-
stitutive
parameters,
therefore,
establishing
the
continuous
dependence
is
crucial
for
validating
the
mathematical
model’s
real-world
applicability.
A
comprehensive
treatment
of
structural
stability
and
related
continuous
de-
pendence
results
can
be
found
in
the
monographs
by
Ames
and
Straughan
[
1
]
and
Straughan
[
19
].
For
Newtonian
fluids,
structural
stability
of
the
flows
through
porous
media
with
the
Darcy
and
Brinkman
model
has
been
ex-
tensively
investigated
by
Morro
and
Straughan
[
13
],
Ames
and
Payne
[
2
],
and
Payne
and
Straughan
[
14
–
16
].
For
the
Brinkman–Forchheimer
model,
related
stability
results
of
Newtonian
fluid
flows
were
established
by
Payne
and
Straughan
[
17
],
C
¸elebi
et
al.
[
6
,
7
],
and
Liu
[
11
].
Lukaszewicz
[
12
]
estab-
lished
continuous
dependence
of
two-dimensional
micropolar
fluid
flows
on
the
microrotation
viscosity
coefficient.
To
the
best
knowledge
of
the
authors,
structural
stability
of
the
flows
in
curved
pipe
is
not
yet
explored.
There-
fore,
in
this
paper,
we
investigate
the
structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
with
Darcy-Lapwood-Forchheimer-Brinkman
model
on
the
pertinent
physical
and
geometric
parameters.
The
mathematical
model
of
the
flow
problem
under
consideration
is
presented
in
Subsection
1.1
.
A
priori
estimates
that
are
required
for
establishing
structural
stability
are
presented
in
Section
2
.
Main
theorems
establishing
the
structural
stability
of
the
flow
model
are
presented,
along
with
their
proofs,
in
Section
3
.
1.1
Mathematical
model
Figure
1:
Flow
Geometry.
We
consider
the
flow
of
a
micropolar
fluid
within
a
stationary
curved
pipe
of
arbitrary
cross-section
of
diameter
2
a
,
which
is
characterized
by
a
A.
Prabal,
M.
Devakar
387
constant
centerline
radius
R
>
a
.
A
typical
section
of
the
curved
pipe
is
shown
in
Figure
1
.
We
assume
that
the
pipe’s
cross-section,
Ω
⊂
R
2
,
is
a
bounded
domain
with
Lipschitz
boundary
∂
Ω.
The
flow
is
driven
by
a
constant
pressure
gradient
along
the
axis
of
the
pipe
and
it
is
assumed
to
be
axisymmetric.
Under
these
assumptions,
using
toroidal
coordinate
system
[
3
],
the
partial
differential
equations
governing
the
flow
of
micropolar
fluid
through
porous
medium
with
Darcy-Lapwood-Forchheimer-Brinkman
model
are
given
by
[
8
–
10
,
18
]
∇
·
(
H
u
)
=
0
,
(1)
ReH
(
u
·
∇
)
u
−
Reδu
2
3
ˆ
a
2
+
Reδu
2
u
3
ˆ
a
3
=
f
−
ReH
∇
p
+
ReG
ˆ
a
3
+
mH
∇
×
ν
+
mδν
3
ˆ
a
1
+
(1
+
m
)
∇
·
(
H
∇
u
)
−
(1
+
m
)
δ
2
H
u
2
ˆ
a
2
−
(1
+
m
)
δ
2
H
u
3
ˆ
a
3
−
σ
2
H
u
−
ReσC
F
H
|
u
|
u
,
(2)
ReHJ
(
u
·
∇
)
ν
−
ReJδu
3
ν
3
ˆ
a
2
+
ReJδu
3
ν
2
ˆ
a
3
=
l
−
2
mHν
ν
ν
+
mH
∇
×
u
+
mδu
3
ˆ
a
1
+
n
1
∇
·
(
H
∇
ν
ν
ν
)
−
n
1
δ
2
H
(
ν
2
ˆ
a
2
+
ν
3
ˆ
a
3
)
+(
n
2
+
n
3
)
H
∇
(
∇
·
ν
ν
)
+
(
n
2
+
n
3
)
δ
∇
ν
2
−
(
n
2
+
n
3
)
δ
2
H
ν
2
ˆ
a
2
(3)
in
Ω
,
where,
the
velocity
vector
u
=
(
u
1
,
u
2
,
u
3
)
and
the
microrotation
vector
ν
=
(
ν
1
,
ν
2
,
ν
3
)
are
independent
of
x
3
,
Re
=
ρUa/µ
is
the
Reynolds
num-
ber,
δ
=
a/R
∈
[0
,
1)
is
the
curvature
ratio,
σ
=
a/
√
K
is
the
porosity
parameter,
H
=
1
+
δx
2
is
the
scale
factor,
G
is
the
constant
axial
pres-
sure
gradient
so
that
the
total
pressure
is
p
(
x
1
,
x
2
)
−
G
x
3
,
m
=
k
µ
is
the
micropolarity
parameter,
J
is
the
nondimensional
gyration
parameter,
and
n
1
,
n
2
and
n
3
are
nondimensional
gyro-viscosity
parameters.
We
assume
that
min
{
n
1
,
n
2
+
n
3
}
>
0.
The
constants
µ
and
k
represents
the
usual
New-
tonian
viscosity
and
the
microrotation
viscosity.
ρ
is
the
density
of
the
fluid,
K
is
the
permeability
of
the
porous
medium
and
C
F
is
Forchheimer
coeffi-
cient.
The
functions
f
and
l
are
the
prescribed
dimensionless
body
force
and
dimensionless
bode
couple,
respectively,
U
is
the
characteristic
velocity
and
{
ˆ
a
1
,
ˆ
a
2
,
ˆ
a
3
}
is
the
orthonormal
basis
of
the
toroidal
coordinate
system
[
3
,
18
].
The
differential
operators
∇
=
∂
∂x
1
,
∂
∂x
2
,
0
and
∇·
=
∂
∂x
1
+
∂
∂x
2
.
Since
the
curved
pipe
is
stationary,
the
no-slip
and
the
hyperstick
bound-
ary
conditions
to
be
satisfied
are,
respectively,
u
i
=
0
and
ν
i
=
0
on
∂
Ω
for
i
=
1
,
2
,
3
.
(4)
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
388
2
A
priori
estimates
The
authors
have
established
the
existence
and
uniqueness
of
a
solutions
to
the
boundary
value
problem
(BVP)
(
1
)-(
4
)
in
[
18
].
More
precisely,
it
was
proved
that
if
f
,
l
∈
(
L
2
(Ω))
3
then
for
sufficiently
small
values
of
δ
,
there
exists
(
u
,
p,
ν
)
∈
X
×
L
2
0
(Ω)
×
H
1
0
(Ω)
3
such
that
the
triplet
(
u
,
p,
ν
)
satisfy
(
1
)-(
4
)
in
the
sense
of
distributions,
where
X
=
u
=
(
u
1
,
u
2
,
u
3
)
∈
H
1
0
(Ω)
3
:
∂
(
Hu
1
)
∂x
1
+
∂
(
Hu
2
)
∂x
2
=
0
.
It
was
also
proved
that
the
solution
(
u
,
p,
ν
)
is
unique
for
sufficiently
small
Re
and
m
.
A
priori
estimates
for
u
and
ν
,
which
were
presented
in
[
18
],
can
be
simplified
as
u
1
,
2
,
Ω
=
∇
u
2
≤
C
1
(1
+
m
)
H
0
m
(1+
δ
)
ν
ν
1
,
2
,
Ω
+
F
2
+
f
2
≡
E
u
(5)
and
ν
1
,
2
,
Ω
=
∇
ν
ν
2
≤
C
2
m
(1+
δ
)
(1+
m
)
H
0
F
2
+
f
2
+
l
2
n
1
H
0
−
C
3
δ
(
n
2
+
n
3
)
≡
E
ν
,
(6)
where
C
1
,
C
2
and
C
3
are
positive
constants
depending
only
on
Ω,
F
=
(0
,
0
,
Re
G
)
and
H
0
=
inf
Ω
H
.
To
derive
an
estimate
for
p
,
we
express
equation
(
2
)
as
Re
∇
p
=
(1
+
m
)
∇
2
u
+
F
,
(7)
in
which
F
=
f
H
+
F
H
−
Re
(
u
·
∇
)
u
+
Reδ
H
(
u
2
3
ˆ
a
2
−
u
2
u
3
ˆ
a
3
)
+
δ
H
∂
u
∂x
2
−
(1
+
m
)
δ
2
H
2
(
u
2
ˆ
a
2
+
u
3
ˆ
a
3
)
−
σ
2
u
−
ReσC
F
|
u
|
u
+
m
∇
×
ν
+
mδ
H
ν
3
ˆ
a
1
.
(8)
From
the
expression
of
F
,
it
is
clear
that
F
∈
H
−
1
(Ω)
3
;
therefore,
using
Lemma
IV.1.9
of
[
4
],
we
get
Re
p
2
≤
C
4
(1
+
m
)
u
1
,
2
,
Ω
+
F
∗
,
(9)
A.
Prabal,
M.
Devakar
389
where
C
4
is
a
positive
constant
depending
only
on
Ω
and
·
∗
denotes
the
operator
norm
in
H
−
1
(Ω)
3
.
Now,
we
find
F
∗
to
get
appropriate
estimate
for
p
.
For
υ
υ
∈
H
1
0
(Ω)
3
,
F
,υ
υ
H
−
1
,H
1
0
=
Ω
f
H
+
F
H
−
Re
(
u
·
∇
)
u
+
Reδ
H
(
u
2
3
ˆ
a
2
−
u
2
u
3
ˆ
a
3
)
+
δ
H
∂
u
∂x
2
−
(1
+
m
)
δ
2
H
2
(
u
2
ˆ
a
2
+
u
3
ˆ
a
3
)
−
σ
2
u
−
ReσC
F
|
u
|
u
+
m
∇
×
ν
+
mδ
H
ν
3
ˆ
a
1
·
υ
υ
dx
(10)
Using
Sobolev
embedding
[
5
]
and
H¨
older’s
inequality,
we
obtain
F
,υ
υ
υ
H
−
1
,H
1
0
≤
C
5
1
H
0
f
2
+
F
2
+
Re
1
+
σC
F
+
δ
H
0
∇
u
2
2
+
σ
2
+
δ
H
0
+
(1
+
m
)
δ
2
H
2
0
∇
u
2
+
m
1
+
δ
H
0
∇
ν
2
∇
υ
υ
2
,
and
hence
F
∗
≤
C
5
1
H
0
f
2
+
F
2
+
Re
1
+
σC
F
+
δ
H
0
E
2
u
+
σ
2
+
δ
H
0
+
(1
+
m
)
δ
2
H
2
0
E
u
+
m
1
+
δ
H
0
E
ν
,
(11)
where
C
5
is
a
positive
constant
depending
only
on
Ω.
Using
(
11
)
in
(
9
),
we
get
an
estimate
for
p
as
p
2
≤
E
p
,
(12)
E
p
=
C
6
1
H
0
Re
f
2
+
F
2
+
1
+
σC
F
+
δ
H
0
E
2
u
+
1
Re
σ
2
+
δ
H
0
+
(1
+
m
)
1
+
δ
2
H
2
0

E
u
+
m
Re
1
+
δ
H
0
E
ν
,
where
C
6
is
a
positive
constant
depending
only
on
Ω.
Using
a
priori
estimates
(
5
),
(
6
)
and
(
12
),
in
the
next
section,
we
inves-
tigate
the
structural
stability
of
the
flow
problem
(
1
)-(
4
);
particularly,
we
establish
the
continuous
dependence
of
the
solution
on
the
key
physical
and
geometric
parameters,
namely,
the
Reynolds
number
(
Re
),
the
curvature
ra-
tio
(
δ
),
the
Forchheimer
coefficient
(
C
F
),
and
the
micropolarity
parameter
(
m
).
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
390
3
Structural
stability
We
assume
that
f
,
l
∈
L
2
(Ω)
3
and
the
values
of
δ
,
Re
and
m
are
sufficiently
small
to
ensure
the
existence
and
uniqueness
of
the
weak
solutions
of
the
problem
(
1
)-(
4
).
Throughout
this
section,
we
denote
C
as
a
generic
positive
constant
depending
only
on
Ω.
The
main
results
of
the
paper
are
presented
in
the
following
four
sub-
sections;
indeed,
Theorems
1
,
2
,
3
and
4
present
the
continuous
dependence
of
the
solution
on
parameters
Re
,
δ
,
C
F
and
m
,
respectively.
3.1
Continuous
dependence
on
the
Reynolds
number
(
Re
)
Let
(
u
,
p,
ν
)
and
(
v
,
q,
η
)
be
the
solutions
of
the
BVP
(
1
)-(
4
)
corresponding
to
the
Reynolds
numbers
Re
1
and
Re
2
,
respectively.
Then
the
difference
(
w
,
π,
φ
)
:=
(
u
−
v
,
p
−
q,
ν
−
η
)
satisfies
the
BVP
∇
·
(
H
w
)
=
0
,
(13)
Re
1
H
(
u
·
∇
)
w
+
(
w
·
∇
)
v
+
ReH
(
v
·
∇
)
v
−
Re
1
δ
(
u
3
+
v
3
)
w
3
ˆ
a
2
−
Reδv
2
3
ˆ
a
2
+
Re
1
δ
(
w
2
u
3
+
v
2
w
3
)ˆ
a
3
+
Reδv
2
v
3
ˆ
a
3
=
Re
G
ˆ
a
3
−
Re
1
H
∇
π
−
ReH
∇
q
+
mH
∇
×
φ
+
mδφ
3
ˆ
a
1
+(1
+
m
)
∇
·
(
H
∇
w
)
−
(1
+
m
)
δ
2
H
(
w
2
ˆ
a
2
+
w
3
ˆ
a
3
)
−
σ
2
H
w
−
Re
1
σC
F
H
(
|
u
|
u
−
|
v
|
v
)
−
ReσC
F
H
|
v
|
v
,
(14)
Re
1
HJ
(
u
·
∇
)
φ
φ
φ
+
(
w
·
∇
)
η
+
ReHJ
(
v
·
∇
)
η
−
Re
1
Jδu
3
φ
3
ˆ
a
2
−
Re
1
Jδw
3
η
3
ˆ
a
2
−
ReJδv
3
η
3
ˆ
a
2
+
Re
1
Jδu
3
φ
2
ˆ
a
3
+
Re
1
Jδw
3
η
2
ˆ
a
3
+
ReJδv
3
η
2
ˆ
a
3
=
−
2
mHφ
φ
φ
+
mH
∇
×
w
+
mδw
3
ˆ
a
1
+
n
1
∇
·
(
H
∇
φ
)
−
n
1
δ
2
H
(
φ
2
ˆ
a
2
+
φ
3
ˆ
a
3
)
+(
n
2
+
n
3
)
H
∇
(
∇
·
φ
)
+
(
n
2
+
n
3
)
δ
∇
φ
φ
−
(
n
2
+
n
3
)
δ
2
H
φ
2
ˆ
a
2
(15)
in
Ω
,
and
w
=
0
=
φ
on
∂
Ω
,
(16)
A.
Prabal,
M.
Devakar
391
where
w
=
(
w
1
,
w
2
,
w
3
),
φ
=
(
φ
1
,
φ
2
,
φ
3
),
w
i
=
u
i
−
v
i
,
φ
i
=
ν
i
−
η
i
(
i
=
1
,
2
,
3)
and
Re
=
Re
1
−
Re
2
.
The
main
result
regarding
the
continuous
dependence
on
the
Reynolds
number
(
Re
)
is
stated
in
the
following
theorem.
Theorem
1.
If
(
w
,
π,φ
φ
)
is
the
solution
of
(
13
)
-
(
16
)
,
then,
for
sufficiently
small
values
of
Re
1
,
m
and
δ
,
w
→
0
in
X
and
φ
φ
→
0
in
(
H
1
0
(Ω))
3
as
Re
2
→
Re
1
.
Proof.
Taking
the
L
2
(Ω)-inner
product
of
(
14
)
with
w
,
we
get
Re
1
Ω
H
(
u
·
∇
)
w
·
w
dx
+
Re
1
Ω
H
(
w
·
∇
)
v
·
w
dx
+
Re
Ω
H
(
v
·
∇
)
v
·
w
dx
−
Re
1
δ
Ω
v
3
w
3
w
2
dx
−
Reδ
Ω
v
2
3
w
2
dx
+
Re
1
δ
Ω
v
2
w
2
3
dx
+
Reδ
Ω
v
2
v
3
w
3
dx
=
Re
G
Ω
w
3
dx
−
Re
1
Ω
H
∇
π
·
w
dx
−
Re
Ω
H
∇
q
·
w
dx
+
m
Ω
H
(
∇
×
φ
φ
)
·
w
dx
+
mδ
Ω
φ
3
w
1
dx
+(1
+
m
)
Ω
∇
·
(
H
∇
w
)
·
w
dx
−
(1
+
m
)
δ
2
Ω
1
H
(
w
2
2
+
w
2
3
)
dx
−
σ
2
Ω
H
w
·
w
dx
−
Re
1
σC
F
Ω
H
(
|
u
|
u
−
|
v
|
v
)
·
w
dx
−
ReσC
F
Ω
H
|
v
|
v
·
w
dx.
(17)
Since
(
|
u
|
u
i
−
|
v
|
v
i
)
w
i
=
1
2
(
|
u
|
+
|
v
|
)
w
2
i
+
1
2
(
|
u
|
−
|
v
|
)
2
(
|
u
|
+
|
v
|
)
and
H
is
a
positive
function,
we
have
Ω
H
(
|
u
|
u
−
|
v
|
v
)
·
w
dx
≥
0.
Also,
considering
the
nonnegativity
of
Ω
1
H
(
w
2
2
+
w
2
3
)
dx
together
with
(
13
)
and
the
fact
that
Ω
H
(
u
·
∇
)
w
·
w
dx
=
0,
the
equality
(
17
)
becomes
(1
+
m
)
H
0
∇
w
2
2
+
σ
2
H
0
w
2
2
≤
Re
1
Ω
H
(
w
·
∇
)
v
·
w
dx
+
Re
1
δ
Ω
v
3
w
3
w
2
dx
+
Re
1
δ
Ω
v
2
w
2
3
dx
+
Re
Ω
H
(
v
·
∇
)
v
·
w
dx
+
Reδ
Ω
v
2
3
w
2
dx
+
Reδ
Ω
v
2
v
3
w
3
dx
+
Re
G
Ω
w
3
dx
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
392
+
ReσC
F
Ω
H
|
v
|
v
·
w
dx
+
m
Ω
H
(
∇
×
φ
φ
)
·
w
dx
+
mδ
Ω
φ
3
w
1
dx
.
(18)
Using
H¨
older’s
inequality,
Sobolev
inequality
and
the
bound
(
5
),
the
terms
of
the
right
hand
side
of
(
18
)
can
be
estimated
as
follows:
Re
1
Ω
H
(
w
·
∇
)
v
·
w
dx
≤
Re
1
H
1
w
4
∇
v
2
w
4
≤
C
Re
1
E
v
∇
w
2
2
,
Re
1
δ
Ω
v
3
w
3
w
2
dx
≤
Re
1
δ
v
2
w
2
4
≤
C
Re
1
δE
v
∇
w
2
2
,
Re
Ω
H
(
v
·
∇
)
v
·
w
dx
≤
|
Re
|
H
1
v
4
∇
v
2
w
4
≤
C
|
Re
|
E
2
v
∇
w
2
,
Reδ
Ω
v
2
3
w
2
dx
≤
|
Re
|
δ
v
2
4
w
2
≤
C
|
Re
|
δE
2
v
∇
w
2
,
Re
G
Ω
w
3
dx
≤
|
Re
||
G
||
Ω
|
w
2
≤
C
|
Re
|∇
w
2
,
ReσC
F
Ω
H
|
v
|
v
·
w
dx
≤
|
Re
|
σC
F
H
1
|
v
|
v
2
w
2
≤
C
|
Re
|
σC
F
∇
v
2
2
∇
w
2
(see
[
18
])
≤
C
|
Re
|
σC
F
E
2
v
∇
w
2
,
m
Ω
H
(
∇
×
φ
)
·
w
dx
≤
mH
1
∇
×
φ
φ
2
w
2
≤
C
m
∇
φ
φ
φ
2
∇
w
2
,
A.
Prabal,
M.
Devakar
393
and
mδ
Ω
φ
3
w
1
dx
≤
mδ
φ
φ
2
w
2
≤
C
mδ
∇
φ
φ
2
∇
w
2
.
The
preceding
estimates
transform
the
inequality
(
18
)
into
(1
+
m
)
H
0
∇
w
2
2
≤
C
Re
1
(1
+
δ
)
E
v
∇
w
2
2
+
|
Re
|
(1
+
δ
+
σC
F
)
E
2
v
+
1
∇
w
2
+
m
(1
+
δ
)
∇
φ
φ
2
∇
w
2
.
(19)
Taking
the
L
2
(Ω)-inner
product
of
(
15
)
with
φ
φ
,
we
get
Re
1
J
Ω
H
(
u
·
∇
)
φ
·
φ
dx
+
Re
1
J
Ω
H
(
w
·
∇
)
η
·
φ
φ
dx
+
ReJ
Ω
H
(
v
·
∇
)
η
·
φ
dx
−
Re
1
Jδ
Ω
w
3
η
3
φ
2
dx
−
ReJδ
Ω
v
3
η
3
φ
2
dx
+
Re
1
Jδ
Ω
w
3
η
2
φ
3
dx
+
ReJδ
Ω
v
3
η
2
φ
3
dx
≤
−
2
m
Ω
Hφ
φ
φ
·
φ
φ
dx
+
m
Ω
H
(
∇
×
w
)
·
φ
dx
+
mδ
Ω
w
3
φ
1
dx
+
n
1
Ω
∇
·
(
H
∇
φ
φ
)
·
φ
φ
dx
−
n
1
δ
2
Ω
1
H
(
φ
2
2
+
φ
2
3
)
dx
+(
n
2
+
n
3
)
Ω
H
∇
(
∇
·
φ
φ
)
·
φ
φ
dx
+
(
n
2
+
n
3
)
δ
Ω
(
∇
φ
φ
)
·
φ
φ
dx
−
(
n
2
+
n
3
)
δ
2
Ω
1
H
φ
2
2
dx.
(20)
Integration
by
parts
allows
us
to
replace
the
integral
(
n
2
+
n
3
)
Ω
H
∇
(
∇
·
φ
φ
φ
)
·
φ
φ
dx
by
−
(
n
2
+
n
3
)
Ω
H
(
∇
·
φ
φ
φ
)
2
dx
−
(
n
2
+
n
3
)
δ
Ω
φ
2
(
∇
·
φ
φ
φ
)
dx
in
the
inequality
(
20
);
subsequently,
using
the
nonnegativity
of
Ω
H
(
∇
·
φ
φ
φ
)
2
dx
,
Ω
1
H
(
φ
2
2
+
φ
2
3
)
dx
and
Ω
1
H
φ
2
2
dx
and
the
fact
that
Ω
H
(
u
·∇
)
φ
φ
)
·
φ
dx
=
0,
we
obtain
n
1
H
0
∇
φ
2
2
+
2
mH
0
φ
φ
2
2
≤
Re
1
J
Ω
H
(
w
·
∇
)
η
·
φ
dx
+
Re
1
Jδ
Ω
w
3
η
3
φ
2
dx
+
Re
1
Jδ
Ω
w
3
η
2
φ
3
dx
+
(
n
2
+
n
3
)
δ
Ω
φ
2
(
∇
·
φ
φ
φ
)
dx
+
(
n
2
+
n
3
)
δ
Ω
(
∇
φ
φ
φ
)
·
φ
φ
φ
dx
+
ReJ
Ω
H
(
v
·
∇
)
η
·
φ
φ
φ
dx
+
ReJδ
Ω
v
3
η
3
φ
2
dx
+
ReJδ
Ω
v
3
η
2
φ
3
dx
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
394
+
m
Ω
H
(
∇
×
w
)
·
φ
dx
+
mδ
Ω
w
3
φ
1
dx
.
(21)
The
terms
present
in
right
hand
side
of
(
21
)
can
be
estimated
as
Re
1
J
Ω
H
(
w
·
∇
)
η
·
φ
dx
≤
Re
1
JH
1
w
4
∇
η
2
φ
4
≤
C
Re
1
JE
η
∇
w
2
∇
φ
φ
2
,
Re
1
Jδ
Ω
w
3
η
3
φ
2
dx
≤
Re
1
Jδ
w
4
η
2
φ
φ
φ
4
≤
C
Re
1
JδE
η
∇
w
2
∇
φ
φ
2
,
(
n
2
+
n
3
)
δ
Ω
φ
2
(
∇
·
φ
)
dx
≤
(
n
2
+
n
3
)
δ
φ
φ
2
∇
·
φ
φ
2
≤
C
(
n
2
+
n
3
)
δ
∇
φ
φ
2
2
,
(
n
2
+
n
3
)
δ
Ω
(
∇
φ
φ
)
·
φ
φ
φ
dx
≤
(
n
2
+
n
3
)
δ
∇
φ
φ
φ
2
φ
φ
2
≤
C
(
n
2
+
n
3
)
δ
∇
φ
φ
φ
2
2
,
ReJ
Ω
H
(
v
·
∇
)
η
·
φ
φ
dx
≤
|
Re
|
JH
1
v
4
∇
η
2
φ
φ
4
≤
C
|
Re
|
JE
v
E
η
∇
φ
2
,
ReJδ
Ω
v
3
η
3
φ
2
dx
≤
|
Re
|
Jδ
v
4
η
2
φ
φ
4
≤
C
|
Re
|
JδE
v
E
η
∇
φ
2
,
m
Ω
H
(
∇
×
w
)
·
φ
φ
dx
≤
mH
1
∇
×
w
2
φ
2
≤
C
m
∇
w
2
∇
φ
φ
2
,
and
mδ
Ω
w
3
φ
1
dx
≤
mδ
w
2
φ
φ
2
≤
C
mδ
∇
w
2
∇
φ
φ
φ
2
.
A.
Prabal,
M.
Devakar
395
Combining
all
these
bounds
with
(
21
),
we
obtain
n
1
H
0
∇
φ
φ
2
2
≤
C
(1
+
δ
)(
m
+
Re
1
JE
η
)
∇
w
2
∇
φ
2
+(
n
2
+
n
3
)
δ
∇
φ
φ
2
2
+
|
Re
|
J
(1
+
δ
)
E
v
E
η
∇
φ
φ
2
.
(22)
Adding
(
19
)
and
(
22
),
we
get
(1
+
m
)
H
0
∇
w
2
2
+
n
1
H
0
∇
φ
2
2
≤
C
Re
1
(1
+
δ
)
E
v
∇
w
2
2
+
|
Re
|
(1
+
δ
+
σC
F
)
E
2
v
+
1
∇
w
2
+(1
+
δ
)(
m
+
Re
1
JE
η
)
∇
w
2
∇
φ
φ
φ
2
+
(
n
2
+
n
3
)
δ
∇
φ
φ
2
2
+
|
Re
|
J
(1
+
δ
)
E
v
E
η
∇
φ
2
.
(23)
By
Young’s
inequality,
we
obtain
(1
+
m
)
H
0
∇
w
2
2
+
n
1
H
0
∇
φ
2
2
≤
(1
+
m
)
H
0
2
∇
w
2
2
+
n
1
H
0
2
∇
φ
φ
φ
2
2
+
C
(
n
2
+
n
3
)
δ
∇
φ
φ
φ
2
2
+
C
Re
1
(1
+
δ
)
E
v
+
(1
+
δ
)
2
(
m
+
Re
1
JE
η
)
2
n
1
H
0
∇
w
2
2
+
C
(1
+
δ
+
σC
F
)
E
2
v
+
1
2
(1
+
m
)
H
0
+
J
2
(1
+
δ
)
2
E
2
v
E
2
η
n
1
H
0
Re
2
.
(24)
After
simplification,
(
24
)
yields
(1
+
m
)
H
0
2
−
C
Re
1
(1
+
δ
)
E
v
+
(1
+
δ
)
2
(
m
+
Re
1
JE
η
)
2
n
1
H
0

∇
w
2
2
+
n
1
H
0
2
−
C
(
n
2
+
n
3
)
δ
∇
φ
φ
2
2
≤
C
(1
+
δ
+
σC
F
)
E
2
v
+
1
2
(1
+
m
)
H
0
+
J
2
(1
+
δ
)
2
E
2
v
E
2
η
n
1
H
0
Re
2
.
(25)
It
follows
from
inequality
(
25
)
that
the
required
result
holds
whenever
Re
1
≤
1
,
m
≤
2
and
δ
≤
3
for
some
positive
constants
1
,
2
and
3
.
Remark
1.
The
convergence
of
v
to
u
in
X
and
η
to
ν
in
H
1
0
(Ω)
3
as
Re
2
→
0
follows
from
the
inequality
(
25
)
when
Re
1
=
0
.
Therefore,
it
can
be
concluded
that
the
BVP
(
1
)
-
(
4
)
admits
a
regular
perturbation
with
respect
to
the
Reynolds
number
Re
.
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
396
3.2
Continuous
dependence
on
the
curvature
ratio
(
δ
)
Suppose
that
(
u
,
p,
ν
)
and
(
v
,
q,
η
)
are
the
solutions
of
the
BVP
(
1
)-(
4
)
cor-
responding
to
the
curvature
ratios
δ
1
and
δ
2
,
respectively.
Let
(
w
,
π,
φ
)
:=
(
u
−
v
,
p
−
q,
ν
−
η
).
Then
(
w
,
π,
φ
)
satisfies
the
BVP
∇
·
(
H
δ
1
w
)
+
˜
δ
∇
·
(
x
2
v
)
=
0
,
(26)
ReH
δ
1
(
u
·
∇
)
w
+
(
w
·
∇
)
v
+
˜
δRex
2
(
v
·
∇
)
v
−
Reδ
1
w
3
(
u
3
+
v
3
)ˆ
a
2
−
˜
δRev
2
3
ˆ
a
2
+
Reδ
1
(
w
2
u
3
+
v
2
w
3
)ˆ
a
3
+
˜
δRev
2
v
3
ˆ
a
3
=
−
ReH
δ
1
∇
π
−
˜
δRex
2
∇
q
+
mH
δ
1
∇
×
φ
φ
+
˜
δmx
2
∇
×
η
+
mδ
1
φ
3
ˆ
a
1
+
˜
δmη
3
ˆ
a
1
+
(1
+
m
)
∇
·
(
H
δ
1
∇
w
)
+
˜
δ
(1
+
m
)
∇
·
(
x
2
∇
v
)
−
(1
+
m
)
δ
2
1
H
δ
1
(
w
2
ˆ
a
2
+
w
3
ˆ
a
3
)
−
˜
δ
(1
+
m
)
δ
1
+
δ
2
H
δ
1
−
δ
2
2
x
2
H
δ
1
H
δ
2
(
v
2
ˆ
a
2
+
v
3
ˆ
a
3
)
−
σ
2
H
δ
1
w
−
˜
δσ
2
x
2
v
−
ReσC
F
H
δ
1
(
|
u
|
u
−
|
v
|
v
)
−
˜
δReσC
F
x
2
|
v
|
v
,
(27)
ReH
δ
1
J
(
u
·
∇
)
φ
φ
+
(
w
·
∇
)
η
+
˜
δReJx
2
(
v
·
∇
)
η
−
ReJδ
1
(
u
3
φ
3
+
w
3
η
3
)ˆ
a
2
−
˜
δReJv
3
η
3
ˆ
a
2
+
ReJδ
1
(
u
3
φ
2
+
w
3
η
2
)ˆ
a
3
+
˜
δReJv
3
η
2
ˆ
a
3
=
−
2
mH
δ
1
φ
φ
φ
−
2
˜
δmx
2
η
+
mH
δ
1
∇
×
w
+
˜
δmx
2
∇
×
v
+
mδ
1
w
3
ˆ
a
1
+
˜
δmv
3
ˆ
a
1
+
n
1
∇
·
(
H
δ
1
∇
φ
φ
)
+
˜
δn
1
∇
·
(
x
2
∇
η
)
−
n
1
δ
2
1
H
δ
1
(
φ
2
ˆ
a
2
+
φ
3
ˆ
a
3
)
−
˜
δn
1
δ
1
+
δ
2
H
δ
1
−
δ
2
2
x
2
H
δ
1
H
δ
2
(
η
2
ˆ
a
2
+
η
3
ˆ
a
3
)
+
(
n
2
+
n
3
)
H
δ
1
∇
(
∇
·
φ
φ
)
+
˜
δ
(
n
2
+
n
3
)
x
2
∇
(
∇
·
η
)
+
(
n
2
+
n
3
)
δ
1
∇
φ
2
+
˜
δ
(
n
2
+
n
3
)
∇
η
2
−
(
n
2
+
n
3
)
δ
2
1
H
δ
1
φ
2
ˆ
a
2
−
˜
δ
(
n
2
+
n
3
)
δ
1
+
δ
2
H
δ
1
−
δ
2
2
x
2
H
δ
1
H
δ
2
η
2
ˆ
a
2
(28)
in
Ω
,
and
w
=
0
=
φ
on
∂
Ω
,
(29)
where
w
=
(
w
1
,
w
2
,
w
3
),
φ
=
(
φ
1
,
φ
2
,
φ
3
),
w
i
=
u
i
−
v
i
,
φ
i
=
ν
i
−
η
i
(
i
=
1
,
2
,
3),
˜
δ
=
δ
1
−
δ
2
and
H
δ
j
=
1
+
δ
j
x
2
(
j
=
1
,
2).
The
principal
result
concerning
the
continuous
dependence
on
the
cur-
vature
ratio
(
δ
)
is
presented
below:
A.
Prabal,
M.
Devakar
397
Theorem
2.
Let
(
w
,
π,φ
φ
φ
)
be
the
solution
of
the
boundary
value
problem
(
26
)
-
(
29
)
.
Then
w
→
0
in
X
and
φ
φ
φ
→
0
in
(
H
1
0
(Ω))
3
as
δ
2
→
δ
1
for
sufficiently
small
values
of
Re
,
m
and
δ
1
.
Proof.
Taking
the
L
2
(Ω)-inner
product
on
both
sides
of
(
27
)
with
w
,
we
get
Re
Ω
H
δ
1
(
u
·
∇
)
w
·
w
dx
+
Re
Ω
H
δ
1
(
w
·
∇
)
v
·
w
dx
+
˜
δRe
Ω
x
2
v
·
∇
)
v
·
w
dx
−
Reδ
1
Ω
(
u
3
+
v
3
)
w
3
w
2
dx
−
˜
δRe
Ω
v
2
3
w
2
dx
+
Reδ
1
Ω
(
w
2
u
3
+
v
2
w
3
)
w
3
dx
+
˜
δRe
Ω
v
2
v
3
w
3
dx
=
−
Re
Ω
H
δ
1
∇
π
·
w
dx
−
˜
δRe
Ω
x
2
∇
q
·
w
dx
+
m
Ω
H
δ
1
(
∇
×
φ
)
·
w
dx
+
˜
δm
Ω
x
2
(
∇
×
η
)
·
w
dx
+
mδ
1
Ω
φ
3
w
1
dx
+
˜
δm
Ω
η
3
w
1
dx
+(1
+
m
)
Ω
∇
·
(
H
δ
1
∇
w
)
·
w
dx
+
˜
δ
(1
+
m
)
Ω
∇
·
(
x
2
∇
v
)
·
w
dx
−
(1
+
m
)
δ
2
1
Ω
1
H
δ
1
(
w
2
2
+
w
2
3
)
dx
−
˜
δ
(1
+
m
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
v
2
w
2
+
v
3
w
3
)
dx
+
˜
δδ
2
2
(1
+
m
)
Ω
x
2
H
δ
1
H
δ
2
(
v
2
w
2
+
v
3
w
3
)
dx
−
σ
2
Ω
H
δ
1
w
·
w
dx
−
˜
δσ
2
Ω
x
2
v
·
w
dx
−
ReσC
F
H
δ
1
Ω
(
|
u
|
u
−
|
v
|
v
)
·
w
dx
−
˜
δReσC
F
Ω
x
2
|
v
|
v
·
w
dx.
(30)
Using
the
identity
Ω
H
δ
1
(
u
·∇
)
w
·
w
dx
=
0,
and
the
nonnegativity
of
the
integrals
Ω
1
H
δ
1
(
w
2
2
+
w
2
3
)
dx
and
Ω
(
|
u
|
u
−
|
v
|
v
)
·
w
dx
,
the
equality
(
30
)
gives
(1
+
m
)
H
δ
1
∇
w
2
2
+
σ
2
H
δ
1
w
2
2
≤
Re
Ω
H
δ
1
(
w
·
∇
)
v
·
w
dx
+
˜
δRe
Ω
x
2
v
·
∇
)
v
·
w
dx
+
Reδ
1
Ω
v
3
w
3
w
2
dx
+
˜
δRe
Ω
v
2
3
w
2
dx
+
Reδ
1
Ω
v
2
w
2
3
dx
+
˜
δRe
Ω
v
2
v
3
w
3
dx
+
Re
Ω
H
δ
1
∇
π
·
w
dx
+
˜
δRe
Ω
x
2
∇
q
·
w
dx
+
m
Ω
H
δ
1
(
∇
×
φ
φ
φ
)
·
w
dx
+
˜
δm
Ω
x
2
(
∇
×
η
)
·
w
dx
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
398
+
mδ
1
Ω
φ
3
w
1
dx
+
˜
δm
Ω
η
3
w
1
dx
+
˜
δ
(1
+
m
)
Ω
∇
·
(
x
2
∇
v
)
·
w
dx
+
˜
δ
(1
+
m
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
v
2
w
2
+
v
3
w
3
)
dx
+
˜
δδ
2
2
(1
+
m
)
Ω
x
2
H
δ
1
H
δ
2
(
v
2
w
2
+
v
3
w
3
)
dx
+
˜
δσ
2
Ω
x
2
v
·
w
dx
+
˜
δReσC
F
Ω
x
2
|
v
|
v
·
w
dx
.
(31)
By
H¨
older’s
and
Sobolev
inequalities
together
with
(
26
),
the
terms
of
the
right-hand
side
of
(
31
)
can
be
estimated
as
follows:
Re
Ω
H
δ
1
(
w
·
∇
)
v
·
w
dx
≤
C
ReE
v
∇
w
2
2
,
˜
δRe
Ω
x
2
v
·
∇
)
v
·
w
dx
≤
C
|
˜
δ
|
ReE
2
v
∇
w
2
,
Reδ
1
Ω
v
3
w
3
w
2
dx
≤
C
Reδ
1
E
v
∇
w
2
2
,
˜
δRe
Ω
v
2
3
w
2
dx
≤
C
|
˜
δ
|
ReE
2
v
∇
w
2
,
Re
Ω
H
δ
1
∇
π
·
w
dx
=
Re
Ω
π
∇
·
(
H
δ
1
w
)
dx
=
Re
Ω
π
−
˜
δ
∇
·
(
x
2
v
)
dx
(by
(
26
))
≤
|
˜
δ
|
Re
π
2
x
2
∇
·
v
+
v
2
2
≤
C
|
˜
δ
|
Re
(
p
2
+
q
2
)
∇
v
2
≤
C
|
˜
δ
|
Re
(
E
p
+
E
q
)
E
v
,
˜
δRe
Ω
x
2
∇
q
·
w
dx
=
˜
δRe
Ω
q
∇
·
(
x
2
w
)
dx
=
˜
δRe
Ω
q
(
x
2
∇
·
w
+
w
2
)
dx
≤
|
˜
δ
|
Re
q
2
x
2
∞
∇
·
w
2
+
w
2
≤
C
|
˜
δ
|
ReE
q
∇
w
2
,
A.
Prabal,
M.
Devakar
399
m
Ω
H
δ
1
(
∇
×
φ
φ
)
·
w
dx
≤
C
m
∇
φ
φ
2
∇
w
2
,
˜
δm
Ω
x
2
(
∇
×
η
)
·
w
dx
≤
C
|
˜
δ
|
mE
η
∇
w
2
,
mδ
1
Ω
φ
3
w
1
dx
≤
C
mδ
1
∇
φ
2
∇
w
2
,
˜
δm
Ω
η
3
w
1
dx
≤
C
|
˜
δ
|
mE
η
∇
w
2
,
˜
δ
(1
+
m
)
Ω
∇
·
(
x
2
∇
v
)
·
w
dx
=
˜
δ
(1
+
m
)
Ω
x
2
∇
v
:
∇
w
dx
≤
C
|
˜
δ
|
(1
+
m
)
E
v
∇
w
2
,
˜
δ
(1
+
m
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
v
2
w
2
+
v
3
w
3
)
dx
≤
C
|
˜
δ
|
(1
+
m
)
(
δ
1
+
δ
2
)
H
0
δ
1
E
v
∇
w
2
,
˜
δδ
2
2
(1
+
m
)
Ω
x
2
H
δ
1
H
δ
2
(
v
2
w
2
+
v
3
w
3
)
dx
≤
C
|
˜
δ
|
(1
+
m
)
δ
2
2
H
0
δ
1
H
0
δ
2
E
v
∇
w
2
,
˜
δσ
2
Ω
x
2
v
·
w
dx
≤
C
|
˜
δ
|
σ
2
E
v
∇
w
2
,
and
˜
δReσC
F
Ω
x
2
|
v
|
v
·
w
dx
≤
C
|
˜
δ
|
ReσC
F
E
2
v
∇
w
2
,
where
H
0
δ
1
:=
inf
Ω
H
δ
1
and
H
0
δ
2
:=
inf
Ω
H
δ
2
.
In
view
of
these,
the
inequality
(
31
)
becomes
(1
+
m
)
H
0
δ
1
∇
w
2
2
≤
C
Re
(1
+
δ
1
)
E
v
∇
w
2
2
+
|
˜
δ
|∇
w
2
Re
(1
+
σC
F
)
E
2
v
+
ReE
q
+
mE
η
+(1
+
m
)
1
+
δ
1
+
δ
2
H
0
δ
1
+
δ
2
2
H
0
δ
1
H
0
δ
2
E
v
+
σ
2
E
v
+
|
˜
δ
|
Re
(
E
p
+
E
q
)
E
v
+
m
(1
+
δ
1
)
∇
φ
φ
2
∇
w
2
.
(32)
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
400
Taking
the
L
2
(Ω)-inner
product
on
both
sides
of
(
28
)
with
φ
φ
φ
,
we
get
ReJ
Ω
H
δ
1
(
u
·
∇
)
φ
·
φ
dx
+
ReJ
Ω
H
δ
1
(
w
·
∇
)
η
·
φ
φ
dx
+
˜
δReJ
Ω
x
2
(
v
·
∇
)
η
·
φ
φ
dx
−
ReJδ
1
Ω
w
3
η
3
φ
2
dx
−
˜
δReJ
Ω
v
3
η
3
φ
2
dx
+
ReJδ
1
Ω
w
3
η
2
φ
3
dx
+
˜
δReJ
Ω
v
3
η
2
φ
3
dx
=
−
2
m
Ω
H
δ
1
φ
·
φ
φ
dx
−
2
˜
δm
Ω
x
2
η
·
φ
dx
+
m
Ω
H
δ
1
(
∇
×
w
)
·
φ
φ
dx
+
˜
δm
Ω
x
2
(
∇
×
v
)
·
φ
φ
dx
+
mδ
1
Ω
w
3
φ
1
dx
+
˜
δm
Ω
v
3
φ
1
dx
+
n
1
Ω
∇
·
(
H
δ
1
∇
φ
φ
)
·
φ
φ
dx
+
˜
δn
1
Ω
∇
·
(
x
2
∇
η
)
·
φ
φ
dx
−
n
1
δ
2
1
Ω
1
H
δ
1
(
φ
2
2
+
φ
2
3
)
dx
−
˜
δn
1
(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
η
2
φ
2
+
η
3
φ
3
)
dx
+
˜
δn
1
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
(
η
2
φ
2
+
η
3
φ
3
)
dx
+
(
n
2
+
n
3
)
Ω
H
δ
1
∇
(
∇
·
φ
φ
)
·
φ
φ
dx
+
˜
δ
(
n
2
+
n
3
)
Ω
x
2
∇
(
∇
·
η
)
·
φ
φ
dx
+
(
n
2
+
n
3
)
δ
1
Ω
∇
φ
2
·
φ
φ
dx
+
˜
δ
(
n
2
+
n
3
)
Ω
∇
η
2
·
φ
dx
−
(
n
2
+
n
3
)
δ
2
1
Ω
1
H
δ
1
φ
2
2
dx
−
˜
δ
(
n
2
+
n
3
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
η
2
φ
2
dx
+
˜
δ
(
n
2
+
n
3
)
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
η
2
φ
2
dx.
(33)
After
simplification,
(
33
)
turns
out
to
be
n
1
H
0
δ
1
∇
φ
φ
2
2
+
2
mH
0
δ
1
φ
φ
φ
2
2
≤
ReJ
Ω
H
δ
1
(
w
·
∇
)
η
·
φ
φ
dx
+
˜
δReJ
Ω
x
2
(
v
·
∇
)
η
·
φ
dx
+
ReJδ
1
Ω
w
3
η
3
φ
2
dx
+
˜
δReJ
Ω
v
3
η
3
φ
2
dx
+
ReJδ
1
Ω
w
3
η
2
φ
3
dx
+
˜
δReJ
Ω
v
3
η
2
φ
3
dx
+
2
˜
δm
Ω
x
2
η
·
φ
φ
dx
+
m
Ω
H
δ
1
(
∇
×
w
)
·
φ
φ
dx
+
˜
δm
Ω
x
2
(
∇
×
v
)
·
φ
dx
+
mδ
1
Ω
w
3
φ
1
dx
+
˜
δm
Ω
v
3
φ
1
dx
+
˜
δn
1
Ω
∇
·
(
x
2
∇
η
)
·
φ
φ
φ
dx
A.
Prabal,
M.
Devakar
401
+
˜
δn
1
(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
η
2
φ
2
+
η
3
φ
3
)
dx
+
˜
δn
1
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
(
η
2
φ
2
+
η
3
φ
3
)
dx
+
(
n
2
+
n
3
)
δ
1
Ω
φ
2
(
∇
·
φ
φ
)
dx
+
˜
δ
(
n
2
+
n
3
)
Ω
x
2
(
∇
·
η
)(
∇
·
φ
φ
)
dx
+
˜
δ
(
n
2
+
n
3
)
Ω
φ
2
(
∇
·
η
)
dx
+
(
n
2
+
n
3
)
δ
1
Ω
∇
φ
2
·
φ
φ
φ
dx
+
˜
δ
(
n
2
+
n
3
)
Ω
∇
η
2
·
φ
φ
φ
dx
+
˜
δ
(
n
2
+
n
3
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
η
2
φ
2
dx
+
˜
δ
(
n
2
+
n
3
)
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
η
2
φ
2
dx
.
(34)
The
terms
appearing
on
the
right-hand
side
of
(
34
)
admit
the
following
bounds.
ReJ
Ω
H
δ
1
(
w
·
∇
)
η
·
φ
φ
dx
≤
C
ReJE
η
∇
w
2
∇
φ
2
,
˜
δReJ
Ω
x
2
(
v
·
∇
)
η
·
φ
φ
φ
dx
≤
C
|
˜
δ
|
ReJE
v
E
η
∇
φ
φ
φ
2
,
ReJδ
1
Ω
w
3
η
3
φ
2
dx
≤
C
ReJδ
1
E
η
∇
w
2
∇
φ
2
,
˜
δReJ
Ω
v
3
η
3
φ
2
dx
≤
C
|
˜
δ
|
ReJE
v
E
η
∇
φ
φ
2
,
2
˜
δm
Ω
x
2
η
·
φ
φ
φ
dx
≤
C
|
˜
δ
|
mE
η
∇
φ
φ
φ
2
,
m
Ω
H
δ
1
(
∇
×
w
)
·
φ
dx
≤
C
m
∇
w
2
∇
φ
φ
2
,
˜
δm
Ω
x
2
(
∇
×
v
)
·
φ
φ
φ
dx
≤
C
|
˜
δ
|
mE
v
∇
φ
φ
φ
2
,
mδ
1
Ω
w
3
φ
1
dx
≤
C
mδ
1
∇
w
2
∇
φ
2
,
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
402
˜
δm
Ω
v
3
φ
1
dx
≤
C
|
˜
δ
|
mE
v
∇
φ
φ
2
,
˜
δn
1
Ω
∇
·
(
x
2
∇
η
)
·
φ
φ
φ
dx
=
˜
δn
1
Ω
x
2
∇
η
:
∇
φ
φ
dx
≤
C
|
˜
δ
|
n
1
E
η
∇
φ
φ
2
,
˜
δn
1
(
δ
1
+
δ
2
)
Ω
1
H
δ
1
(
η
2
φ
2
+
η
3
φ
3
)
dx
≤
C
|
˜
δ
|
n
1
(
δ
1
+
δ
2
)
H
0
δ
1
E
η
∇
φ
2
,
˜
δn
1
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
(
η
2
φ
2
+
η
3
φ
3
)
dx
≤
C
|
˜
δ
|
n
1
δ
2
2
H
0
δ
1
H
0
δ
2
E
η
∇
φ
φ
2
,
(
n
2
+
n
3
)
δ
1
Ω
φ
2
(
∇
·
φ
φ
)
dx
≤
C
(
n
2
+
n
3
)
δ
1
∇
φ
2
2
,
˜
δ
(
n
2
+
n
3
)
Ω
x
2
(
∇
·
η
)(
∇
·
φ
φ
φ
)
dx
≤
|
˜
δ
|
(
n
2
+
n
3
)
E
η
∇
φ
φ
φ
2
,
˜
δ
(
n
2
+
n
3
)
Ω
φ
2
(
∇
·
η
)
dx
≤
|
˜
δ
|
(
n
2
+
n
3
)
E
η
∇
φ
φ
2
,
(
n
2
+
n
3
)
δ
1
Ω
∇
φ
2
·
φ
φ
dx
≤
C
(
n
2
+
n
3
)
δ
1
∇
φ
2
2
,
˜
δ
(
n
2
+
n
3
)
Ω
∇
η
2
·
φ
φ
dx
≤
C
|
˜
δ
|
(
n
2
+
n
3
)
E
η
∇
φ
2
,
˜
δ
(
n
2
+
n
3
)(
δ
1
+
δ
2
)
Ω
1
H
δ
1
η
2
φ
2
dx
≤
C
|
˜
δ
|
(
n
2
+
n
3
)(
δ
1
+
δ
2
)
H
0
δ
1
E
η
∇
φ
φ
2
,
and
˜
δ
(
n
2
+
n
3
)
δ
2
2
Ω
x
2
H
δ
1
H
δ
2
η
2
φ
2
dx
≤
C
|
˜
δ
|
(
n
2
+
n
3
)
δ
2
2
H
0
δ
1
H
0
δ
2
E
η
∇
φ
φ
φ
2
.
A.
Prabal,
M.
Devakar
403
Taking
into
account
the
preceding
estimates,
(
34
)
gives
n
1
H
0
δ
1
∇
φ
φ
φ
2
2
≤
C
δ
1
(
n
2
+
n
3
)
∇
φ
φ
2
2
+
(1
+
δ
1
)
m
+
ReJE
η
∇
w
2
∇
φ
φ
2
+
|
˜
δ
|∇
φ
φ
2
(
n
1
+
n
2
+
n
3
)
1
+
δ
1
+
δ
2
H
0
δ
1
+
δ
2
2
H
0
δ
1
H
0
δ
2
E
η
+
mE
η
+
m
+
ReJE
η
E
v

.
(35)
Adding
(
32
)
and
(
35
),
we
get
(1
+
m
)
H
0
δ
1
∇
w
2
2
+
n
1
H
0
δ
1
∇
φ
φ
2
2
≤
C
Re
(1
+
δ
1
)
E
v
∇
w
2
2
+
δ
1
(
n
2
+
n
3
)
∇
φ
φ
φ
2
2
+
(1
+
δ
1
)
m
+
ReJE
η
∇
w
2
∇
φ
φ
φ
2
+
|
˜
δ
|
D
1
∇
w
2
+
|
˜
δ
|
D
2
∇
φ
φ
2
+
|
˜
δ
|
Re
(
E
p
+
E
q
)
E
v
,
(36)
where
D
1
=
Re
(1
+
σC
F
)
E
2
v
+
ReE
q
+
mE
η
+
(1
+
m
)
1
+
δ
1
+
δ
2
H
0
δ
1
+
δ
2
2
H
0
δ
1
H
0
δ
2
E
v
+
σ
2
E
v
and
D
2
=
(
n
1
+
n
2
+
n
3
)
1
+
δ
1
+
δ
2
H
0
δ
1
+
δ
2
2
H
0
δ
1
H
0
δ
2
E
η
+
mE
η
+
m
+
ReJE
η
E
v
.
Applying
Young’s
inequality,
(
36
)
yields
(1
+
m
)
H
0
δ
1
2
∇
w
2
2
+
n
1
H
0
δ
1
2
∇
φ
φ
2
2
≤
C
Re
(1
+
δ
1
)
E
v
∇
w
2
2
+
δ
1
(
n
2
+
n
3
)
∇
φ
φ
2
2
+
(1
+
δ
1
)
2
m
+
ReJE
η
2
(1
+
m
)
H
0
δ
1
∇
φ
2
2
+
D
2
1
(1
+
m
)
H
0
δ
1
+
D
2
2
n
1
H
0
δ
1
˜
δ
2
+
Re
(
E
p
+
E
q
)
E
v
|
˜
δ
|
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
404
which
implies
(1
+
m
)
H
0
δ
1
2
−
C
Re
(1
+
δ
1
)
E
v
∇
w
2
2
+
n
1
H
0
δ
1
2
−
C
δ
1
(
n
2
+
n
3
)
−
C
(1
+
δ
1
)
2
m
+
ReJE
η
2
(1
+
m
)
H
0
δ
1
∇
φ
φ
2
2
≤
C
D
2
1
(1
+
m
)
H
0
δ
1
+
D
2
2
n
1
H
0
δ
1
˜
δ
2
+
C
Re
(
E
p
+
E
q
)
E
v
|
˜
δ
|
.
(37)
The
desired
result
follows
from
(
37
)
for
sufficiently
small
values
of
Re
,
m
,
and
δ
1
.
Remark
2.
In
the
particular
case
when
δ
1
=
0
and
δ
2
→
0
,
inequality
(
37
)
implies
that
the
solution
v
converges
to
u
in
X
and
η
converges
to
ν
in
H
1
0
(Ω)
3
.
Hence,
we
conclude
that
the
solution
of
the
curved
pipe
flow
converges
to
that
of
the
straight
pipe
flow
when
δ
→
0
.
3.3
Continuous
dependence
on
the
Forchheimer
coefficient
(
C
F
)
Let
(
u
,
p,
ν
)
and
(
v
,
q,
η
)
be
the
solutions
of
the
BVP
(
1
)-(
4
)
correspond-
ing
to
the
Forchheimer
coefficients
C
F
1
and
C
F
2
,
respectively.
Then
the
differences
(
w
,
π,
φ
)
:=
(
u
−
v
,
p
−
q,
ν
−
η
)
solve
the
problem
∇
·
(
H
w
)
=
0
,
(38)
ReH
(
u
·
∇
)
w
+
(
w
·
∇
)
v
−
Reδu
3
w
3
ˆ
a
2
−
Reδv
3
w
3
ˆ
a
2
+
Reδu
3
w
2
ˆ
a
3
+
Reδw
3
v
2
ˆ
a
3
=
−
ReH
∇
π
+
mH
∇
×
φ
φ
+
mδφ
3
ˆ
a
1
+
(1
+
m
)
∇
·
(
H
∇
w
)
−
(1
+
m
)
δ
2
H
(
w
2
ˆ
a
2
+
w
3
ˆ
a
3
)
−
σ
2
H
w
−
ReσC
F
1
H
(
|
u
|
u
−
|
v
|
v
)
−
C
F
ReσH
|
v
|
v
,
(39)
A.
Prabal,
M.
Devakar
405
ReHJ
(
u
·
∇
)
φ
φ
φ
+
(
w
·
∇
)
η
−
ReJδu
3
φ
3
ˆ
a
2
−
ReJδw
3
η
3
ˆ
a
2
+
ReJδu
3
φ
2
ˆ
a
3
+
ReJδw
3
η
2
ˆ
a
3
=
−
2
mHφ
φ
+
mH
∇
×
w
+
mδw
3
ˆ
a
1
+
n
1
∇
·
(
H
∇
φ
)
−
n
1
δ
2
H
(
φ
2
ˆ
a
2
+
φ
3
ˆ
a
3
)
+
(
n
2
+
n
3
)
H
∇
(
∇
·
φ
)
+
(
n
2
+
n
3
)
δ
∇
φ
2
−
(
n
2
+
n
3
)
δ
2
H
φ
2
ˆ
a
2
(40)
in
Ω
,
and
w
=
0
=
φ
on
∂
Ω
,
(41)
where
w
=
(
w
1
,
w
2
,
w
3
),
φ
=
(
φ
1
,
φ
2
,
φ
3
),
w
i
=
u
i
−
v
i
,
φ
i
=
ν
i
−
η
i
(
i
=
1
,
2
,
3)
and
C
F
=
C
F
1
−
C
F
2
.
The
following
theorem
establishes
the
continuous
dependence
of
the
so-
lution
on
the
Forchheimer
coefficient
(
C
F
).
Theorem
3.
Suppose
that
(
w
,
π,φ
φ
)
solves
the
BVP
(
38
)
-
(
41
)
.
Then,
for
sufficiently
small
values
of
Re
,
m
and
δ
,
w
→
0
in
X
and
φ
φ
→
0
in
(
H
1
0
(Ω))
3
as
C
F
2
→
C
F
1
.
Proof.
Taking
the
L
2
(Ω)-inner
product
of
(
39
)
with
w
yields
the
estimate
(1
+
m
)
H
0
∇
w
2
2
≤
C
Re
(1
+
δ
)
E
v
∇
w
2
2
+
m
(1
+
δ
)
∇
φ
2
∇
w
2
+
C
F
ReσE
2
v
∇
w
2
.
(42)
Taking
the
L
2
(Ω)-inner
product
of
(
40
)
with
φ
,
we
obtain
the
estimate
n
1
H
0
∇
φ
φ
2
2
≤
C
(1
+
δ
)
m
+
ReJE
η
∇
w
2
∇
φ
φ
φ
2
+
δ
(
n
2
+
n
3
)
∇
φ
φ
2
2
.
(43)
Adding
(
42
)
and
(
43
)
and
using
the
Young’s
inequality,
we
get
(1
+
m
)
H
0
2
−
C
Re
(1
+
δ
)
E
v
∇
w
2
2
+
n
1
H
0
−
C
δ
(
n
2
+
n
3
)
−
C
(1
+
δ
)
2
m
+
ReJE
η
2
(1
+
m
)
H
0
∇
φ
φ
2
2
≤
C
Re
2
σ
2
E
4
v
(1
+
m
)
H
0
C
F
2
.
(44)
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
406
The
inequality
(
44
)
shows
that
∇
w
2
→
0
and
∇
φ
φ
2
→
0
as
C
F
→
0
for
sufficiently
small
values
of
Re
,
m
and
δ
.
Hence,
the
required
result
is
proved.
Remark
3.
For
a
particular
case
C
F
1
=
0
,
when
C
F
2
→
0
,
the
inequality
(
44
)
implies
that
the
solution
to
the
Darcy–Lapwood
–Forchheimer–Brinkman
model
for
micropolar
fluid
flow
in
a
curved
pipe
converges
to
that
of
the
Darcy–Lapwood–Brinkman
model.
3.4
Continuous
dependence
on
the
micropolarity
parameter
(
m
)
Let
(
u
,
p,
ν
)
and
(
v
,
q,
η
)
be
the
solutions
of
the
BVP
(
1
)-(
4
)
with
micropo-
larity
parameters
m
1
and
m
2
,
respectively.
The
difference
fields
(
w
,
π,
φ
)
:=
(
u
−
v
,
p
−
q,
ν
−
η
),
then
satisfy
the
problem
∇
·
(
H
w
)
=
0
,
(45)
ReH
(
u
·
∇
)
w
+
(
w
·
∇
)
v
−
Reδu
3
w
3
ˆ
a
2
−
Reδv
3
w
3
ˆ
a
2
+
Reδu
3
w
2
ˆ
a
3
+
Reδw
3
v
2
ˆ
a
3
=
−
ReH
∇
π
+
m
1
H
∇
×
φ
+
mH
∇
×
η
+
m
1
δφ
3
ˆ
a
1
+
mδη
3
ˆ
a
1
+(1
+
m
1
)
∇
·
(
H
∇
w
)
+
m
∇
·
(
H
∇
v
)
−
(1
+
m
1
)
δ
2
H
(
w
2
ˆ
a
2
+
w
3
ˆ
a
3
)
−
m
δ
2
H
(
v
2
ˆ
a
2
+
v
3
ˆ
a
3
)
−
σ
2
H
w
−
ReσC
F
H
(
|
u
|
u
−
|
v
|
v
)
,
(46)
ReHJ
(
u
·
∇
)
φ
φ
+
(
w
·
∇
)
η
−
ReJδu
3
φ
3
ˆ
a
2
−
ReJδw
3
η
3
ˆ
a
2
+
ReJδu
3
φ
2
ˆ
a
3
+
ReJδw
3
η
2
ˆ
a
3
=
−
2
m
1
Hφ
φ
−
2
mH
η
+
m
1
H
∇
×
w
+
mH
∇
×
v
+
m
1
δw
3
ˆ
a
1
+
mδv
3
ˆ
a
1
+
n
1
∇
·
(
H
∇
φ
φ
φ
)
−
n
1
δ
2
H
(
φ
2
ˆ
a
2
+
φ
3
ˆ
a
3
)
+
(
n
2
+
n
3
)
H
∇
(
∇
·
φ
φ
)
+(
n
2
+
n
3
)
δ
∇
φ
2
−
(
n
2
+
n
3
)
δ
2
H
φ
2
ˆ
a
2
(47)
in
Ω
,
and
w
=
0
=
φ
on
∂
Ω
,
(48)
A.
Prabal,
M.
Devakar
407
where
w
=
(
w
1
,
w
2
,
w
3
),
φ
=
(
φ
1
,
φ
2
,
φ
3
),
w
i
=
u
i
−
v
i
,
φ
i
=
ν
i
−
η
i
(
i
=
1
,
2
,
3)
and
m
=
m
1
−
m
2
.
Our
main
result
on
the
continuous
dependence
on
the
micropolarity
pa-
rameter
(
m
)
is
formulated
below.
Theorem
4.
If
(
w
,
π,φ
φ
φ
)
is
the
solution
of
(
45
)
-
(
48
)
,
then
w
→
0
in
X
and
φ
φ
φ
→
0
in
(
H
1
0
(Ω))
3
as
m
2
→
m
1
,
provided
the
values
of
Re
,
m
1
and
δ
are
sufficiently
small.
Proof.
Adopting
the
same
line
of
reasoning
from
the
preceding
results,
we
infer
that
(1
+
m
1
)
H
0
2
−
C
Re
(1
+
δ
)
E
v
∇
w
2
2
+
n
1
H
0
2
−
C
δ
(
n
2
+
n
3
)
−
C
(1
+
δ
)
2
m
1
+
ReJE
η
2
(1
+
m
1
)
H
0
∇
φ
φ
2
2
≤
C
E
η
(1
+
δ
)
+
E
v
1
+
δ
2
H
0

2
(1
+
m
1
)
H
0
+
E
η
+
E
v
(1
+
δ
)
2
n
1
H
0
m
2
.
(49)
It
is
clear
from
(
49
)
that,
for
sufficiently
small
values
of
Re
,
m
1
and
δ
,
the
required
result
holds.
Remark
4.
It
is
well
known
that
in
the
absence
of
microrotational
effects
the
system
of
micropolar
equations
reduced
to
the
system
of
Navier-Stokes
equations;
it
is
also
known
that
when
the
microrotation
viscosity
equals
to
zero,
the
linear
velocity
vector
and
the
microrotation
vector
are
uncoupled
and
the
global
motion
is
unaffected
by
the
microrotation
vector
[
8
].
The
same
can
be
achieved
from
(
49
)
by
taking
m
1
=
0
and
allowing
m
2
→
0
.
References
[1]
K.A.
Ames
and
B.
Straughan,
Non-Standard
and
Improperly
Posed
Problems
,
Academic
Press,
San
Diego,
1997.
[2]
K.A.
Ames
and
L.E.
Payne,
On
stabilizing
against
modeling
errors
in
a
penetrative
convection
problem
for
a
porous
medium,
Math.
Models
Methods
Appl.
Sci.
4
(1994),
733-740.
[3]
N.
Arada,
On
generalized
Newtonian
fluids
in
curved
pipes,
SIAM
J.
Math.
Anal.
48
(2016),
1210-1249.
Structural
stability
of
micropolar
fluid
flow
in
a
curved
pipe
408
[4]
F.
Boyer
and
P.
Fabrie,
Mathematical
Tools
for
the
Study
of
the
Incom-
pressible
Navier-Stokes
Equations
and
Related
Models
,
Springer,
New
York,
2012.
[5]
H.
Brezis,
Functional
Analysis,
Sobolev
Space
and
Partial
Differential
Equations
,
Springer,
New
York,
2011.
[6]
A.O.
C
¸elebi,
V.K.
Kalantarov
and
D.
U˘
gurlu,
Continuous
dependence
for
the
convective
Brinkman–Forchheimer
equations,
Appl.
Anal.
84
(2005),
877-888.
[7]
A.O.
C
¸elebi,
V.K.
Kalantarov
and
D.
U˘
gurlu,
On
continuous
depen-
dence
on
coefficients
of
the
Brinkman–Forchheimer
equations,
Appl.
Math.
Lett.
19
(2006),
801-807.
[8]
A.C.
Eringen,
Theory
of
micropolar
fluids,
J.
Math.
Mech.
16
(1966),
1-18.
[9]
M.H.
Hamdan,
Single-phase
flow
through
porous
channels
a
review
of
flow
models
and
channel
entry
conditions,
Appl.
Math.
Comput.
62
(1994),
203-222.
[10]
C.T.
Hsu
and
P.
Cheng,
Thermal
dispersion
in
a
porous
medium,
Int.
J.
Heat
Mass
Transfer
33
(1990),
1587-1597.
[11]
Y.
Liu,
Convergence
and
continuous
dependence
for
the
Brinkman–
Forchheimer
equations,
Math.
Comput.
Modell.
49
(2009),
1401-1415.
[12]
G.
Lukaszewicz,
Long
time
behavior
of
2D
micropolar
fluid
flows,
Math.
Comput.
Modell.
34
(2001),
487-509.
[13]
A.
Morro
and
B.
Straughan,
Continuous
dependence
on
the
source
pa-
rameters
for
convective
motion
in
porous
media,
Nonlinear
Anal.
18
(1992),
307-315.
[14]
L.E.
Payne
and
B.
Straughan,
Structural
stability
for
the
Darcy
equa-
tions
of
flow
in
porous
media,
Proc.
R.
Soc.
A
454
(1974),
1691-1698.
[15]
L.E.
Payne
and
B.
Straughan,
Stability
in
the
initial-time
geometry
problem
for
the
Brinkman
and
Darcy
equations
of
flow
in
porous
media,
J.
Math.
Pures
Appl.
75
(1996),
225-271.
[16]
L.E.
Payne
and
B.
Straughan,
Analysis
of
the
boundary
condition
at
the
interface
between
a
viscous
fluid
and
a
porous
medium
and
related
modelling
questions,
J.
Math.
Pures
Appl.
77
(1998),
317-354.
A.
Prabal,
M.
Devakar
409
[17]
L.E.
Payne
and
B.
Straughan,
Convergence
and
continuous
dependence
for
the
Brinkman–Forchheimer
equations,
Stud.
Appl.
Math.
102
(1999),
419-439.
[18]
A.
Prabal
and
M.
Devakar,
On
the
existence,
uniqueness
and
regularity
of
solutions
of
micropolar
fluid
flow
through
porous
medium
in
a
curved
pipe,
Nonlinear
Anal.
Real
World
Appl.
91
(2026),
104616.
[19]
B.
Straughan,
Stability
and
Wave
Motion
in
Porous
Media
,
Applied
Mathematical
Sciences,
Springer,
New
York,
2008.