Two-Layered Blood Flow Through A Composite Stenosis in The Presence of A Magnetic Field

Download as pdf or txt
Download as pdf or txt
You are on page 1of 12

International Journal of Application or Innovation in Engineering& Management (IJAIEM)

Web Site: www.ijaiem.org Email: [email protected], [email protected]


Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 30


ABSTRACT
Blood flow through a composite stenosis in the presence of a magnetic field was investigated by examining the effects that an
external, uniform, transverse magnetic field has on its velocity, flow rate and three flow characteristics (impedance, wall shear
stress and shear stress at the stenosis throat). A two-layered model, consisting of a cell-free peripheral plasma layer and a particle-
fluid suspension in the core region, was used to represent blood flowing through a composite stenosis in an artery. The effects of
introducing the magnetic field, as well as increasing its intensity were examined. The fluid's velocity and flow rate were reduced
when the magnetic field was introduced as well as when its intensity was increased. All three of the flow characteristics increased
when the magnetic field was introduced and as its intensity increased. These effects make it more difficult for the blood to flow
and may even lead to wall vessel damage and plaque rupture. When hematocrit, stenosis height and tube length were varied the
same trends were observed in the presence of the external magnetic field as those observed in the absence of the magnetic field.
The two-fluid model flow characteristics were lower than those of the one-fluid model. This knowledge can aid in the improvement
of existing diagnostic tools, especially those used in diagnosing cardiovascular diseases.

Keywords: blood flow, stenosis, magnetic field, impedance, shear stress

1. INTRODUCTION
According to the World Health Organisation (WHO), an estimated 17 million deaths that occur annually are due to
cardiovascular related complications. The 23.3 million deaths predicted by 2030 ensure that cardiovascular diseases will
remain the single leading cause of death worldwide and a popular subject of scientific research [29]. The two main causes
of the associated deaths are due to coronary heart disease (disease of the blood vessels supplying the heart muscle) and
cerebrovascular disease (disease of the blood vessels supplying the brain). A heart attack or a stroke usually follows when
blood flow to the heart or the brain respectively is blocked. Experimental and theoretical blood flow studies which include
magnetic effects have been useful in the design of medical devices that are used in diagnosing cardiovascular diseases.
In several situations, it may be necessary to subject the human body to magnetic fields although the full implications of its
effects on biological fluid dynamics are still being researched. Magnetic Resonance Imaging (MRI) is a tool that is
frequently used in the field of medicine to study blood flow. Low intensity, low frequency magnetic energy has been used
to treat chronic pain secondary to tissue ischemia and non healing or slowly healing ulcers. This energy appeared to affect
biological processes due to its electrically induced charges. New types of cellular phones, levitated trains and
superconducting magnetic energy storage have been designed to utilise electro-magnetic frequencies [26]. Other devices
that contain magnetic fields have already been used in cell separation, targeted drug transport carriers, cancer tumour
treatment, the reduction of bleeding in surgeries and injuries at accidents and wars, magnetic band aids and tracers [28].
Further studies on the effects of an external magnetic field on blood flow will allow greater flexibility for future designs of
such devices.
In medical research, the effect of electro-magnetic fields was first investigated by Kolin [10]. Such studies involving the
flow through a channel in the presence of a transverse magnetic field are encountered in a variety of applications such as
magnetohydrodynamics (MHD) generators, pumps, accelerators and flow meters. The application of an external magnetic
field was then seen as a means to aid in the regulation of the movement of blood in the human circulatory system. Since
blood is electrically conducting, MHD principles were easily applied to its flow. By Lenz's law, the Lorentz's force acting
on the blood opposes its motion thereby reducing its velocity. This reduction in velocity proves useful in the treatment of
the cardiovascular diseases which involve accelerated blood flow, for example in hypertension and hemorrhages. It is also
used to provide analgesia, decrease healing time of fractures and to act as a treatment for depression. Myocardial
infarctions could even be triggered by man-made magnetic fields since it was found that magnetic fields exceeding 0.05 T
can be dangerous to the heart and may result in irreversible changes. Therefore the majority of devices that utilise
magnetic fields contain static magnetic fields rather than electo-magnetic ones. Even static magnetic fields can be
dangerous since they are known to cause a change in calcium dynamics and alter skin blood flow [1].
Pioneering blood flow models have treated blood flow as a single phase, homogeneous, Newtonian, viscous fluid [11].
Non Newtonian models for blood flow were also developed since blood, being a suspension of corpuscles, does not always
exhibit the characteristics of a Newtonian fluid especially in tubes with small diameters [3]. Many blood flow studies also
involved the presence of stenosis (an abnormal growth which occurs under diseased conditions in any location of the
cardiovascular system) [30]. Stenosis has been modeled as a single stenosis (symmetric, bell shaped, non-symmetric,
Two-layered Blood Flow through a Composite
Stenosis in the Presence of a Magnetic Field.

Alana R. Sankar
1
, Sreedhara Rao Gunakala
2
and Donna M.G. Comissiong
3
.

1,2,3
Department of Mathematics and Statistics, The University of the West Indies, Trinidad and Tobago.

International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 31

composite) or a multiple stenosis (overlapping, irregularly shaped) and its effect on blood flow characteristics has been
examined [11], [20]-[23]. Magnetic effects on single fluid (Newtonian, Casson, Micropolar or Couple Stress fluid) models
were also examained [5], [14], [27]. One fluid models examining the effects of MHD on blood flow in stenosed arteries
were proposed by Ikbal et al., Singh and Rathee, Misra et al. and Sut [7], [18], [16], [24].
However, based on Haynes' theoretical analysis and Cokelet's experimental observations, a single-phase homogeneous
viscous fluid analysis cannot be used to realistically describe blood flow in blood vessels of diameter m 1000 s [4], [6].
Additionally, the discovery that hematocrit (red cell volume fraction) in human blood is approximately 45% and has a
strong influence on the blood flow properties of blood has led to the treatment of the red blood cells (erythrocytes) as
discrete particles by Skalak [19]. Thus, Srivastava and Srivastava proposed a two-phase macroscopic model (i.e. a
suspension of red cells in plasma) for blood in small vessels (of diameter m 2400 s ) where the individuality of the red
cells (of diameter m 8 ) is significant [20]. Further experimental work provided evidence for the existence of a layer of
plasma containing no cells located layer near the walls and a core region containing all the erythrocytes suspended in
plasma [2]. Therefore two layered models can be used to describe blood flow in vessels of diameter m. 1000 s
To study blood flow in the presence of MHD, two-layered models were proposed [1], [15]. Ponalagusamy and Selvi
studied oscillatory blood flow in the presence of a magnetic field through a mild stenosis using a two-layered model where
both layers were represented as Newtonian fluids in the presence of a uniform magnetic field using a finite difference
method [15]. A more realistic two-layered model for blood flow through small blood vessels consists of the peripheral cell
free layer and a core two-phase macroscopic layer. This model has been used to study the effect of stenosis on blood flow
characteristics [9], [16], [22], [23]. A magnetic study was conducted by Bali and Awasthi where the red blood cells were
treated as magnetic particles to examine the effect this would have on the fluid's velocity, flow rate and effective viscosity
but in the absence of stenosis and an external magnetic field [1].
Thus the objective of the present analysis is to provide a model to examine the effect of an external magnetic field on
blood flow through a composite stenosis which has not yet been examined in previous literature. This will be done by
using a two-layered blood flow model consisting of a cell-free peripheral layer and a core region of erythrocytes in plasma
flowing through a composite stenosis in the presence of an external, transverse magnetic field. This two-layered model for
blood flow provides a more realistic model for flow in small arteries since the existence of the peripheral layer and the red
blood cells in the plasma can no longer be neglected. The study of a composite stenosis is important since the shape of the
stenosis which is manifested in the arteries varies. The effect that the external magnetic field has on the fluid's velocity,
flow rate and three of its blood flow characteristics, resistance to flow, wall shear stress and shear stress at the stenosis
throat will be examined. Knowledge generated in this area can lead to the improvement of existing diagnostic tools for a
more effective treatment of patients suffering from cancer, hypertension, myocardial infarction, stroke and paralysis. It
can also provide an understanding of the effect that a magnetic field may have on different diseased arterial systems.

2. METHODOLOGY

Application of a uniform, transverse magnetic field in the axial direction
Figure 1 The geometry of the stenosed artery

Consider a two layered model for the flow of blood through a cylindrical vessel which consists of a central core layer of
radius
1
R which is a suspension of red blood cells (erythrocytes) in plasma and a peripheral plasma layer of thickness
)
1
( R R modeled as a Newtonian, viscous fluid. The blood flow is assumed to be axisymmetrical, fully developed, steady
and laminar. The artery is assumed to be uniform, circular and rigid. It is also assumed that a mild, axially nonsymmetric
but radially symmetric composite stenosis has manifested in the artery. A uniform, transverse magnetic field is applied in
the axial direction. The wall near to the stenosis development is treated as being solid in nature. It is assumed that the
length of the artery is much greater than its radius so that entrance, end and special wall effects can be neglected. The
stenosis geometry manifested in the arterial segment is described in Figure 1 as [9], [13], [16]
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 32

2
0
); )(
1
, (
0 0
2
) , 1 (
0
)) (
1
), ( ( L
d z d d z
L R R
z R z R
+ s s = o o o
0
2
0
;
2
0
0
2
cos 1
0
2
)
1
, (
) , 1 ( L d z
L
d
L
d z
L R
+ s s + + =
(
(

|
|
.
|

\
| t
o o
o
otherwise. ); , 1 ( o =
where ) (z R R ~ is the radius of the tube with constriction and
0
R are the radius of the tube without any constriction. The
length of the tube is denoted by L and the length of the stenosis by
0
L . The stenosis starts at d z = . In the unobstructed
region, the ratio of the radius of the central core to that of the tube is o . The maximum height of the stenosis and
bulging of the interface at the location where
2
0
L
d z + = in the stenotic region is )
1
, ( o o .
In the core region, )
1
0 ( R r s s , a two-phase macroscopic model is used. The following additional assumptions must
be made to simplify the analysis due to the complicated structure of blood [20]:
(i) the red cell is spherical in shape and rigid,
(ii) the artery wall is rigid and of infinite length,
(iii) the flow has rotational symmetry,
(iv) a modified Stokes drag law governs the interaction between the two phases,
(v) red cell-cell interaction and Brownian motion is neglected.
Using a continuum approach and these assumptions, the following are the governing equations for the two-phase blood
flow in the region
1
0 R r s s in the presence of Lorentz's force.
For the fluid (plasma) phase:
Equation of axial momentum
.
2
0
) (
2
) ( ) 1 ( ) 1 ( ) 1 (
f
u B
f
u
p
u CS
f
u C
s
C
z
p
C
r
f
u
f
v
z
f
u
f
u
t
f
u
f
C o + V +
c
c
=
c
c
+
c
c
+
c
c

|
|
.
|

\
|
(1)
Equation of radial momentum
). (
2
1
2
) ( ) 1 ( ) 1 ( ) 1 (
f
v
p
v CS
f
v
r
C
s
C
r
p
C
r
f
v
f
v
z
f
v
f
u
t
f
v
f
C + V +
c
c
=
c
c
+
c
c
+
c
c

|
|
.
|

\
|
|
|
.
|

\
|
(2)
Equation of continuity
. 0 ) 1 ( ) 1 (
1
=
c
c
+
c
c
(

f
u C
z
f
v C r
r r
(3)

For the particle (erythrocyte) phase:
Equation of axial momentum
). (
p
u u CS
z
p
C
r
u
v
z
u
u
t
u
C
f
p
p
p
p
p
p
+
c
c
=
c
c
+
c
c
+
c
c
|
|
.
|

\
|
(4)
Equation of radial momentum
). (
p
v v CS
r
p
C
r
v
v
z
v
u
t
v
C
f
p
p
p
p
p
p
+
c
c
=
c
c
+
c
c
+
c
c
|
|
.
|

\
|
(5)
Equation of continuity
( ) ( ) . 0
1
=
c
c
+
c
c
p p
Cu rCv
z r r
(6)
Here ) , ( z r are two-dimensional cylindrical polar coordinates with z measured along the axis of the tube and r
measured normal to the tube axis,
2
2
1
2
z
r
r
r r
c
c
+
c
c
c
c
= V
|
.
|

\
|
is a two-dimensional Laplacian operator, ) , (
f
v
f
u and
) , (
p
v
p
u are the (axial, radial) components of the fluid and particle velocities. C denotes the volume fraction density of
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 33

the particles, p is the pressure,
s
C
s
~ ) ( is the mixture viscosity (apparent or effective viscosity), S is the drag
coefficient of interaction for the force exerted by one phase on the other,
f
and
p
are the actual densities of the
material constituting the fluid (plasma) and the particle (erythrocytes) phases respectively;
f
C ) 1 ( is the fluid phase
density and
p
C is the particulate phase density, and the subscripts f and p denote the quantities associated with the
plasma (fluid) and erythrocyte (particle) phases respectively.
The electrical conductivity of the fluid is o and
0
B is the component of the constant uniform magnetic field which was
applied in the axial direction. The suspension viscosity,
s
, has been chosen as the following empirical relation [3], [20]
,
1
0
) (
mC
C
s
=

(7)
where
. ) 69 . 1 exp(
1107
49 . 2 exp )
2
10 7 (
(

|
.
|

\
|
+

= C
T
C m (8)
Here
0
is the fluid viscosity (suspending medium). The temperature of the blood,T , is measured on the absolute scale
(K). Charm & Kurland performed experiments to investigate the suspension viscosity using a cone and plate viscometer.
A good agreement was found within 10% of the tested cases [3]. It was also found to be reasonably accurate for 6 . 0 s C
[3], [21]-[23]. The expression for the drag coefficient of interaction, S , from the classical Stokes drag was modified to
account for the finite particulate fractional volume is [25]
2
) 3 2 (
3
2
1
)
2
3 8 ( 3 4
2
0
0
5 . 4
C
C C C
a
S

+ +
=
|
|
|
.
|

\
|

(9)
with
0
a as the red blood cell radius.
Using Navier-Stokes equations in the presence of the Lorentz's force, the governing equations in the peripheral region,
R r R s s
1
, can be written as
0
2
0 0
2
0
0
0
0
0
0
0
u B u
z
p
r
u
v
z
u
u
t
u
o V +
c
c
=
c
c
+
c
c
+
c
c
|
|
.
|

\
|
(10)
0
2
1
2
0
0
0
0
0
0
0
v
r
r
p
r
v
v
z
v
u
t
v
|
|
.
|

\
|
|
|
.
|

\
|
V +
c
c
=
c
c
+
c
c
+
c
c
(11)
where )
0
,
0
( v u are the (axial, radial) components and
0
and the density of the fluid.
The boundary conditions (standard no slip conditions of velocities and the shear stresses at the tube wall and the
interface) are given as [13], [16], [21-23]
, 0 at , 0 and 0
,
1
at , and
0
, at , 0
0
= = =
c
c
=
c
c
= = =
= =
r
f r
p
u
r
f
u
R r
f p f
u u
R r u
t
t t
where
r
u
p c
c
=
0
0
t and
r
f
u
s
C
f c
c
= t ) 1 ( are the shear stresses of the peripheral and central regions respectively.
It was also assumed that there is no radial flow along the axis of the tube ( 0
0
= v ) and the axial velocity gradient of the
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 34

fluid and particle phases of the streaming blood is equal to zero
|
|
.
|

\
|
=
c
c
=
c
c
0 i.e.
r
p
v
r
f
v
.
3. RESULT
The following non-dimensional variables are used in this analysis

,
0
0 0
,
0
2
0
2
0 2
,
0
0
0 *
,
0
0 *
,
0
) , ,
0
(
)
*
,
*
,
*
0
( ,
0
) , ,
0
(
)
*
,
*
,
*
0
(
,
0
2
0 *
,
0
*
,
0
*
,
0
1
*
1
,
0
0
*
0
,
0
*
,
0
*
,
0
*
,
0
*

o
o
R U
Re
R B
M
U
pR
p
R
t U
t
U
p
v
f
v v
p
v
f
v v
U
p
u
f
u u
p
u
f
u u
R S
S
s
s
R
R
R
R
R
R
R
L
L
R R
d
d
R
z
z
R
r
r
= = = = = =
= = = = = = = = =

where M is the Hartmann number and Re the Reynolds number. After non-dimensionalising all the equations and
dropping the * notation for simplicity, the following system is obtained


(12)
,
1
0 ,
2
) (
) (
) 1 ( ) 1 ( R r
f
u M
f
u
p
u CS
f
u
r
r
r r
C
s
C
dz
dp
C s s +
c
c
c
c
=
|
.
|

\
|

(13)
,
1
0 ), ( R r
p
u
f
u CS
dz
dp
C s s = (14)
subject to
, at , 0
0
R r u = = (15)
,
1
at , ) 1 (
0
and
0
R r
r
f
u
s
C
r
u
f
u u =
c
c
=
c
c
= (16)
, 0 at , 0 = =
c
c
=
c
c
r
r
p
u
r
f
u
(17)
in the presence of a composite stenosis given by
2
0
); )(
1
, (
0
2
) , 1 ( )) (
1
), ( (
L
d z d d z
L
z R z R + s s = o o o (18)
0
2
0
;
2
0
0
2
cos 1
2
)
1
, (
) , 1 ( L d z
L
d
L
d z
L
+ s s + + =
(
(

|
|
.
|

\
| t
o o
o (19)
otherwise. ); , 1 ( o = (20)
The expressions for velocities,
f
u u ,
0
and
p
u are found to be
,
1
,
2
1
) (
0 2
) (
0 1 0
R r R
dz
dp
M
Mr K k Mr I k u s s + =
|
.
|

\
|
(21)
,
1
0 ,
2
1
) (
0
3
R r
dz
dp
M
r M J k u
f
s s =
|
.
|

\
|
(22)
,
1
0 ,
1
2
1
) (
0
3
R r
dz
dp
S
M
r M J k u
p
s s + =
|
.
|

\
|
|
|
.
|

\
|
(23)
where
,
) (
0
) (
0 2
2
1
1
MR I
MR K k
dz
dp
M
k

=
|
.
|

\
|

,
1
,
0
2
0
1
R r R u M u
r
r
r r dz
dp
s s
c
c
c
c
=
|
.
|

\
|
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 35

| |
,
)]
1
(
0
)
1
(
1
)
1
(
1
) 1 )(
1
(
0
)[ (
0
)]
1
(
1
) 1 )(
1
(
0
)
1
(
1
)
1
(
0
)[ (
0
)
1
(
0
)
1
(
1
)
1
(
1
) 1 )(
1
(
0
2
1
2

|
.
|

\
|
+ +

+
=
R M J MR I R M J
s
C MR I MR K
R M J
s
C MR K MR K R M J MR I
R M J MR I R M J
s
C MR I
dz
dp
M
k




,
)
1
(
1
) 1 (
)
1
(
1 1
)
1
(
1 2
3 R M J
s
C
MR I k MR K k
k


=
and
s
C

) 1 (
1

=
where
n
J is the Bessel functions of first kind of order n ,
n
I is the modified Bessel function of the first kind of order n ,
and
n
K is the modified Bessel function of the second kind of order n .
The flow flux, Q , is now calculated as
,
1
1
0
] ) 1 [(
0
2

} } + + =
R
R
R
dr
p
Cu
f
u C r dr ru Q t
which simplifies to

|
|
|
.
|

\
|
+
|
|
.
|

\
|
+
+ =
S
CR
M
R
dz
dp
R M J R
M
k
MR RK MR K R
M
k
MR I R MR RI
M
k
Q
2
2
1
2
2
2
)
1
(
1 1
(
3
)) (
1
)
1
(
1 1
(
2
))
1
(
1 1
) (
1
(
1
2

t
(24)
Since the sum of the fluxes in the two layers is equal to the total flux, it follows that R R o =
1
and oo o =
1
[23]. Thus
the pressure drop ) at , 0 at ( L z p z p p = = = A across the stenosis is calculated as
} = A
|
.
|

\
|
L
dz
dz
dp
p
0

}
+
+ }
+
+
+
}
+
+ } = A
|
.
|

\
|
|
.
|

\
|
|
.
|

\
|
|
.
|

\
|
L
L d
dz
R
dz
dp
L d
L
d
dz
R
dz
dp
L
d
d
dz
R
dz
dp
d
dz
R
dz
dp
p
0
) 20 .( Eq from
0
2
0
) 19 .( Eq from
2
0
) 18 .( Eq from
0
) 20 .( Eq from
(25)
where
|
.
|

\
|

dz
dp
is obtained from Eq. (24).
To obtain the integrals involved in Eq. (25) in closed form poses some difficulty so they will be evaluated numerically.
Non-dimensionless expressions for the impedance (flow resistance), , the wall shear stress,
w
t and the shear stress
at the stenosis throat,
s
t are obtained by using
.
3
0
0
4
0
with
1
2
and
2
where ,
0
) , (
) , (
,
4
0
0
8
0
and where ,
0
R
R
dz
dp R
s
dz
dp R
w
s w
s w
R
L
Q
p
t

t
o
t t
t
t t
t t
t

=
=
= = =
=
A
= =
(

|
.
|

\
|
|
.
|

\
|

International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 36

Here
0
and
0
t are the impedance and shear stress respectively for a Newtonian fluid in the absence of stenosis and
magnetic field (i.e., 0 , 0 = = M C ), and ) , , (
s w
t t are (impedance, wall shear stress and shear stress at the stenosis
throat) in their dimensional forms.

4. DISCUSSION/ANALYSIS
Plots of the analytical solutions are provided for the following parameter values: d =0;
0
L =1; L =1, 2, 5; C =0, 0.1, 0.2,
0.3, 0.4, 0.5, 0.6; o =0, 0.05, 0.10, 0.15, T =25.5
0
C, to observe the effects of the external magnetic field on flow
characteristics [11], [16], [17], [30]. Values of the Hartmann number, M =0, 2, 4, 6, 8, 10 are used in this analysis [18],
[26]. The parameter o is computed using
0
/ 1 R c o = in which ) (C c c ~ represents the peripheral layer of thickness
as a function of cell concentration. In order for the computation to be done, the temperature was chosen as T =25.5
0
C, to
be able to use Haynes' analysis where ) ( m c =6.18, 4.67, 3.60, 3.12, 2.58, 2.18 corresponds to hematocrit (%) =10, 20,
30, 40, 50, 60 respectively. The parameter
s
is computed with the use of Eqns. (7) and (8). Simpson's rule is used to
evaluate the integral involved in Eq. (25) numerically.
The effect of the external magnetic field on the fluid's velocity, flow rate and three flow characteristics (impedance, wall
shear stress, shear stress at the stenosis throat) will be analysed. Knowledge of how the magnetic field affects flow in the
presence of stenosis can aid in the understanding angina while knowledge of how it affects wall shear stress and the stress
at the stenosis throat can help in the theory of plaque rupture. For abnormal hearts, an increase in wall shear stress may
result in paralysis or death. The study of the resistance to flow (impedance) is physiologically important as it is a
reflection of whether the amount of blood being supplied to the vital organs is sufficient or not [15]. Thus the study of
these hydrodynamic factors is valuable to understanding the development and progression of arterial diseases.
From Figure 2, as the Hartmann number (magnetic parameter) increases, the velocity of the fluid decreases in both the
layers of the flow. It can also be seen that in the absence of the external magnetic field ( 0 = M ), the fluid's velocity is
higher than in its presence ( 0 > M ). Therefore, the presence of an external magnetic field reduces blood's velocity and as
the intensity of the magnetic field is increased it is further reduced which is in agreement with the work of Bali and
Awasthi [2]. This occurs due to the opposing Lorentz's force that is introduced when the magnetic field is applied. Hence
magnetic fields can be used to control blood flow.

Figure 2 The effect of varying the Hartmann number, M , on the fluid's velocity

The flow rate is higher in the absence of the magnetic field than in its presence (see Figure 3). As the intensity of the
magnetic field increases (i.e. M increases) the flow rate decreases. Also generally, as the hematocrit, C , increases, the
flow rate decreases both in the absence and presence of the magnetic field [2].

Figure 3 Variation of flow rate with hematocrit, C for various values of the Hartmann number, M
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 37

Figure 4 shows that as the Hartmann number, M , increases the impedance to flow, , increases and that the presence of
the magnetic field causes the impedance to increase. Thus, the presence of the magnetic field increases the resistance to
flow and as its intensity is increased, the resistance is further increased. This means it becomes more difficult for the
blood to flow through the vessels. Ponalagusamy and Selvi also observed this based on their model [15]. It must also be
noted that in the presence of a magnetic field of any intensity, as the hematocrit increases, the impedance increases and as
the stenosis height increases, the impedance increases. This means that as the number of red blood cells increases and the
larger the stenosis gets, the harder it is for the blood to flow in the vessels which was also observed in the absence of the
magnetic field by Sankar et al. [16].

Figure 4 Impedance versus stenosis height for different M and C values

From Figure 5, in the presence of the external magnetic field ( 2 = M ), as the artery length increases impedance to flow,
, decreases. Therefore the longer the artery is without any branching or curving, the easier it is for blood to flow in that
region even in the presence of a magnetic field. Also, the two layered model's flow impedance was slightly lower than in
the one-fluid model.

Figure 5 Impedance versus stenosis height for different artery lengths (with M =2, C =0.4)

The presence of an external magnetic field also affects the wall shear stress. From Figure 6, it is observed that as
Hartmann number, M , increases, the wall shear stress increases and that the wall shear stress is higher in the presence of
the magnetic field than in its absence. Therefore the presence of the external magnetic field causes a greater force to be
exerted on the vessel wall and as the intensity of the field is increased the force also increases which in the long term can
damage the walls [5]. Additionally, in the presence of the magnetic field, as the hematocrit, C , increases, the wall shear
stress increases which can lead to further damage.
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 38


Figure 6 Wall shear stress distribution,
w
t , in the stenotic region for different M and C values

As the stenosis height, o , increases the shear stress on the wall increases in the presence of the magnetic field (Figure
7). Therefore the larger the stenosis gets, the greater the force is exerted on the wall as it increases rapidly in the upstream
of the stenosis and achieves it maximum at the throat of the stenosis (
2
0
L
d z + = ) before it begins to decrease
downstream to the endpoint of the constriction profile (
0
L z = ). This trend was also observed in the absence of the
magnetic field [16]. The two layered model's wall shear stress was slightly lower than in the one-fluid model.

Figure 7 Wall shear stress distribution,
w
t , in the stenotic region for different stenosis heights, o (with M =2, C =0.4)
From Figure 8, as the Hartmann number, M , is increased the shear stress at the stenosis throat increases and the shear
stress is higher in the presence of the magnetic field than in its absence ( M =0). Thus exposure to an external magnetic
field increases the shear stress at the stenosis throat and as the intensity increases the shear stress is further increased
which may lead to plaque rupture [14]. Also in the presence of the magnetic field as the hematocrit, C , increases the
shear stress at the stenosis throat also increases which can further contribute to the possibility of plaque rupture.
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 39


Figure 8 Shear stress at the stenosis throat,
s
t , versus stenosis height, o for different M and C values
In the presence of the external magnetic field, as the stenosis height increases, the shear stress at the stenosis throat also
increases and that the two layered model's shear stress at the stenosis throat was slightly lower than in the one-fluid
model (see Figure 9). Thus there is an increased need for the consideration of the peripheral layer (a two-layered flow)
when examining blood flow through small blood vessels for different diseases that occur such as plasma cell dyscrasias
(hematocrit=28.00%, o =0.816), Hb SS-sickle cell (hematocrit=24.80%, o =0.795), hypertension-controlled
(hematocrit=43.13%, o =0.928), hypertension-uncontrolled (hematocrit=43.25%, o =0.925) and polycythemia
(hematocrit=63.20%, o =0.990) [2].

Figure 9 Shear stress at the stenosis throat,
s
t , versus stenosis height, o (with M =2)

5. CONCLUSION
In this analysis of blood flow in small vessels in the presence of a magnetic field where a composite stenosis is
manifested, a two layered model is used to model blood. The two layers are a Newtonian peripheral plasma layer and a
core region of erythrocytes suspended in plasma (i.e., particle-fluid mixture). The fluid's velocity, flow rate and three of
its flow characteristics are examined.
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 40

Both the fluid's velocity and flow rate are reduced as the magnetic field is introduced and as its intensity is increased.
This explains why magnetic fields can be used as a control mechanism for blood flow. The impedance to flow is found to
increase as the Hartmann number (magnetic parameter), M , is introduced and as it increases. Thus in the presence of the
magnetic field it is more difficult for the blood to flow through the vessels. This effect is accentuated as the intensity of
the field increases. Both the wall shear stress and shear stress at the stenosis throat increases as Hartmann number, M , is
introduced and as the intensity of the field is increased. This has serious implications because it can cause damage to the
vessel walls and even lead to plaque rupture.
Flow impedance is found to increase with hematocrit and stenosis height but it decreases with tube length. Therefore in
the presence of the external magnetic field, it is harder for the blood to flow if the number of red blood cells increases, if
the stenosis increases in size radially or if a very short blood vessel is being considered. The shear stress on the wall
increases rapidly in the upstream of the stenosis throat, achieves its maximum at the throat of the stenosis, then decreases
downstream. It increases as stenosis height and hematocrit increases, that is, in the presence of a magnetic field, if the
number of red blood cells increases or the stenosis increases in size radially, the wall stress gets larger, which can lead to
it becoming damaged. In the presence of the magnetic field, variations in the shear stress at the stenosis throat and
impedance are similar with respect to all parameters. All three flow characteristics are lower in the two fluid model than
those of the one-fluid model. Therefore as a physiological application to the analysis, the peripheral layer becomes
particularly important in the study of diseased vessels since the thickness of the peripheral layer (and hematocrit in the
blood) varies for different diseases. Some of the popular diseases which affect hematocrit and peripheral plasma layer
thickness of the blood when they occur are plasma cell dyscrasias, Hb SS-sickle cell, hypertension and polycythemia.
In conclusion, low intensity magnetic fields can be beneficial because they reduce the blood's velocity. As the intensity
increases however, it can lead to fatal effects such as blood vessel wall damage and plaque rupture. The restrictions
present in this study (including the rigid wall, steady and fully developed flow and constant thickness of the peripheral
layer) do not prevent it from being useful in understanding the flow phenomena. By considering a fully developed flow
with 1
0
/ << R o , closed form solutions are obtained but this can only be applied to early vessel constriction (mild,
composite stenosis). These solutions are used to find the relationships that the flow characteristics have with the variables
which are indeed valid even though the parameter
0
/ R o is limited to values up to 0.15 due to the possible flow
separation that can occur even at small Reynolds numbers [30]. Studies on two layered flows are not only important to
understanding physiological flows (blood flow, protein diffusion, microorganisms movement, particle deposition on the
respiratory tract) but also to other areas of research such as powder technology, petroleum transport, aerosol filtration,
wastewater treatment, fluidization, mining, power plant piping, corrosive particles in engine oil flow, environmental
pollution, lunar ash flows, atmospheric fallout, combustion, agriculture and food technologies [8], [12]. These studies
involving the effect of the magnetic field on blood flow in stenosed arteries can lead to the improvement of existing
diagnostic tools which contain magnetic fields. The results obtained without considering the presence of a peripheral
layer, stenosis and an external magnetic field, can be used to analyse a mixture of particles in a fluid within a circular
cylinder in any physical situation.

References
[1] R. Bali, U. Awasthi, Mathematical Model of Blood Flow in Small Blood Vessel in the Presence of Magnetic Field,
Applied Mathematics, 2(2), pp. 264-269, 2011.
[2] G. Bugliarello, J. Sevilla, Velocity distribution and other characteristics of Steady and Pulsatile Blood Flow in fine
glass tubes, Biorheology, 7, pp. 85-107. 1970.
[3] S. Charm, G. Kurland, Blood flow and Microcirculation, New York, John Wiley, 1974.
[4] G. Cokelet, The Rheology of Human Blood: In Biomechanics, Prentice-Hall, New Jersey, Englewood Cliffs, 1972.
[5] K. Halder, S. Ghosh, Effect of a magnetic field on blood flow through an indented tube in the presence of
erythrocytes, Indian Journal of Pure and Applied Mathematics, 25, pp. 345-352, 1994.
[6] R. Haynes, Physical basis of the dependence of blood viscosity on tube radius, American Journal of Physiology,
198, pp. 1193-1205, 1960.
[7] M. Ikbal, S. Chakravarty, K. Wong, Mazumda, P. Mandal, Unsteady response of Non-Newtonian Blood Flow
through a Stenosed Artery in Magnetic Field, Journal of Computational and Applied Mathematics, 230, pp. 243-
259, 2009.
[8] B. Jha, C. Ampere, Unsteady MHD two-phase Couette flow of fluid-particle suspension, Applied Mathematical
Modelling, 37, pp. 1920-1931, 2013.
[9] P. Joshi, A. Pathak, B. Joshi, Two-layered model of blood flow through composite stenosed artery, Applications
and Applied Mathematics, 4, pp. 343-354, 2009.
[10] A. Kolin, Electromagnetic flow meter: principle of method and its applications to blood flow measurements,
Proceedings of the Society for Experimental Biology and Medicine, 35, pp. 53, 1936.
[11] D. McDonald, On steady flow modelled cellular stenosis, Journal of Biomechanics, 12, pp. 13-50, 1979.
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: [email protected], [email protected]
Volume 2, Issue 12, December 2013 ISSN 2319 - 4847

Volume 2, Issue 12, December 2013 Page 41

[12] A. Medhavi, U. Singh, A two-layered Suspension Flow induced by Peristaltic Waves, International Journal of Fluid
Mechanics, 35, pp. 258-272, 2008.
[13] A. Medhavi, R. Srivastava, Q. Ahmad, V. Srivastava, Two-phase Arterial Blood Flow through a Composite
Stenosis, e-Journal of Science and Technology, 7, pp. 83-94, 2012.
[14] J. Misra, A. Sinha, G. Shit, Mathematical modelling of Blood Flow in a Porous Vessel having Double Stenoses in
the Presence of an External Magnetic Field, International Journal of Biomathematics, 4(2), pp. 207-225, 2011.
[15] R. Ponalagusamy, R. Selvi, Blood flow in stenosed arteries with radially variable viscosity, peripheral plasma layer
thickness and magnetic field, Meccanica, 2013.
[16] A. Sankar, S. Rao Gunakala, D. Comissiong, Two-Layered Suspension Blood Flow Through a Composite Stenosis,
Journal of Mathematics Research, 5(4), pp. 26-38, 2013.
[17] J. Shukla, R. Parihar, B. Rao, Effects of stenosis on non-Newtonian flow of blood in an artery, Bulletin of
Mathematical Biology, 42, pp. 283-294, 1980.
[18] J. Singh, R. Rathee, Analysis of Non-Newtonian Blood Flow through Stenosed Vessel in Porous Medium under the
effect of Magnetic Field, International J ournal of the Physical Sciences, 6(10), pp. 2497-2506, 2011.
[19] R. Skalak, Mechanics of Microcirculation: In Biomechanics, Its Foundation and Objectives, New Jersey, Prentice-
Hall, Englewood Cliffs, 1972.
[20] L. Srivastava, V. Srivastava, On two-phase model of pulsatile blood flow with entrance effects, Biorheology, 20,
pp. 761-777, 1983.
[21] V. Srivastava, Particle-fluid suspension model of blood flow through stenotic vessels, International Journal of
Biomedical Computing, 38, pp. 141-154, 1995.
[22] V. Srivastava, M. Saxena, Two-layered model of Casson fluid flow through stenotic blood vessels: applications to
the cardiovascular system, J ournal of Biomechanics, 27, pp. 921-928, 1994.
[23] V. Srivastava, M. Saxena, Suspension model for blood flow through stenotic arteries with a cell-free plasma layer,
Mathematical Biosciences, 139, pp. 79-102, 1997.
[24] D. Sut, Study of Blood Flow with Effects of Slip in Arterial Stenosis due to Presence of Transverse Magnetic Field,
International Journal of Mathematical Archive, 3(3), pp. 983-993, 2012.
[25] C. Tam, The drag on a cloud of spherical particles in low Reynolds number flows, Journal of Fluid Mechanics, 38,
pp. 537-546, 1969.
[26] G. Varshney, K. Katiyar, S. Kumar, Effect of Magnetic Field on the Blood Flow in Artery having Multiple Stenosis:
a Numerical Study, International Journal of Engineering, Science and Technology, 2(2), pp. 67-82, 2010.
[27] S. Verma, A. Srivastava, Effect of Magnetic Field on Steady Blood Flow through an Inclined Circular Tube,
International Journal of Engineering Research and Applications, 3(4), pp. 428-432, 2013.
[28] P. Voltairas, D. Fotiadis, L. Michalis, Hydrodynamics of magnetic drug targeting, Journal of Biomechanics, 35,
pp. 813-821, 2002.
[29] World Health Organisation, Cardiovascular diseases (CVDs). 2013. [Online]. Available:
http://www.who.int/mediacentre/factsheets/fs317/en/index.html. [Accessed: Oct. 15, 2013].
[30] D. Young, Fluid mechanics of arterial stenosis, Journal of Biomechanical Engineering, 101, pp. 157-175, 1979.


AUTHOR

Alana Sankar received the B.Sc. General (First Class Honours in Mathematics and Computer Science) in 2009
and the M.Sc. (with distinction in Mathematics) in 2011. She is currently pursuing a Ph. D. in Mathematics at
the University of the West Indies, St. Augustine, Trinidad and Tobago.

You might also like