 Open Access
 Total Downloads : 1978
 Authors : S. Saha, S. Chakrabarti
 Paper ID : IJERTV2IS70010
 Volume & Issue : Volume 02, Issue 07 (July 2013)
 Published (First Online): 01072013
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Impact Of Magnetic Field Strength On Magnetic Fluid Flow Through A Channel
S. Saha1, S. Chakrabarti2
1Dept. of Mechanical Engineering, Dr. Sudhir Chandra Sur Degree Engineering College, 540 Dum Dum Road,
Kolkata – 700074, West Bengal, India
2Department of Mechanical Engineering, Bengal Engineering and Science University, Shibpur, Howrah 711103,
West Bengal, India
ABSTRACT
In this paper an attempt has been made to study the flow characteristics of a magnetic fluid flowing through a channel under the action of magnetic field keeping the magnet along axial direction and changing its length considering FHD phenomenon with respect to the formation of recirculating bubbles, velocity contour and average stagnation pressure. The continuity and momentum equations have been discretized by a control volume based finite difference method. Power law scheme is used to discretize the convective terms. The discretized equations have been solved iteratively by SIMPLE algorithm, using linebyline ADI method. The distribution of grid nodes is nonuniform and staggered in both coordinating directions. The shape and size of the bubble and average stagnation pressure have a greater importance in case of various practical applications, such as magnetic hyperthermia, targeted transport of drugs, magnetic wound or cancer tumour treatment, opening the blockage in the artery etc.

INTRODUCTION
Numerical study of magnetic fluid flow through a channel has become a very demanding research area for researchers. Over the past few decades, eminent researchers have extensively worked on dynamics of magnetic fluid in the presence of magnetic field. They have made an endeavour to implicate it in the industry, medical technology and bioengineering. Magnetic fluid or ferrofluid is well known for its wide use in the industry and in the field of medicine, e.g. magnetic drug targeting, ferrofluid sealing, highgradient magnetic separation, magnetic devices for cell separation, targeted transport of drugs, hyperthermia and reduction of bleeding, cancer treatment, etc. The role of formation of recirculating bubble, velocity contour and average stagnation pressure plays an important role in the stage of above applications.
Literature has been enriched with experimental, numerical and theoretical activities on magnetic fluid flow under the action of magnetic field. From review of literature, it is felt that first work in that field has been carried out by Loukopoulos and E.
E. Tzirzilakis [1]. They have numerically presented two dimensional, laminar, incompressible biomagnetic fluid flow in a rectangular channel under the action of a magnetic field. E. E. Tzirzilakis [2] has presented numerical simulation on laminar, incompressible, three dimensional, fully developed viscous flow of a Newtonian biometric fluid in a straight or curved rectangular duct under the action of spatially varying or uniform magnetic field. E. E. Tzirzilakis et al. [3] have presented two dimensional, incompressible Newtonian, electrically conducting turbulent flow between two parallel plates under the action of localized magnetic field. Hayat et al. [4] have studied the effect of magnetic field on
a turbulent biomagnetic, electrically conducting fluid flow in a rectangular channel. Ikbal et al. [5] have theoretically investigated the atherosclerotic arteries dealing with mathematical models that represent nonNewtonian flow of blood through a stenosed artery in the presence of a transverse magnetic field of magnetic number ranging from 0 to 4 and Reynolds number from 300 to 1000. E. E. Tzirzilakis

has studied the fundamental problem of the biomagnetic fluid flow in a channel with stenosis under the influence of a steady localized magnetic field. Papadopolos et al.

have numerically studied the laminar, incompressible, fully developed viscous flow of a biomagnetic fluid in a rectangular duct under the influence of an applied
which flows through a channel, under the action of variable axial magnetic field. The origin of the Cartesian coordinate system is located at the leading edge of the lower plate and the flow is subjected to a magnetic source which is placed much closed to the lower plate and below it. The flow at the entrance is assumed to be uniform and fully develop at the exit. A schematic diagram of the computational domain is shown in fig.1. In this study, the dimensional velocity components and the pressure are governed by the mass and momentum conservation equations. Thus the dimensional governing equation along the x, y direction with considering FHD are written as follows,
u v 0
magnetic field of magnetic number range of 3Ã—106 to 5Ã—106 and Reynolds number range from 200 to 2000.
x y
————– (1)
u u p
H 2 u
2 u
From the abovementioned review of
u x v y x 0 M
x x 2
y 2
literature, it has been observed that in all the
cases the length of the magnet has not been
———— (2)
varied and magnet has been kept along the transverse direction and also they have not
u v v v p M
0
0
x y y
H 2 v
2
2
y x
2 v
y 2
studied the flow characteristics with respect
——— (3)
to shape and size of the recirculating bubble, velocity contour and average stagnation pressure, so far. Hence in the present work, an attempt has been taken to study the flow characteristics of a biomagnetic fluid flowing through a channel under the action of magnetic field keeping the magnet along axial direction and changing its length considering FHD phenomenon with respect to the formation of recirculating bubbles and velocity contour generated by MATLAB Software, and average stagnation pressure respectively.


MATHEMATICAL FORMULATION

Governing equations
In the above dimensional equations,
is the fluid density, is the viscosity, H is the magnetic field strength intensity, M is the magnetization and 0 is the magnetic permeability under vacuum.
The magnetization linear equation for isothermal cases is given by,
M= H ———(4)
Where is constant called magnetic susceptibility.
In our study, we have considered the following dimensionless variables,
Length:
The keyline of this work is steady, laminar, incompressible biomagnetic fluid
x*
x , y*
W
y , L* L W W
———– (5)
Velocity:
u*
u , v* v U U
—————– (6)
Y*
D (0, W)
C (L, W)
Pressure:
p*
p
U 2
—————— (7)
Inflow
Outflow
Magnetic field strength:
A (0, 0)
E (L/2, 0) B (L, 0) X*
H * H
H 0
————– (8)
Magnetic source
Fig.1. Schematic diagram of the computational domainand contours of the magnetic
field strength H.
Thus with the help of the above non dimensional variables, the governing equation (1) (3) are transformed to the following equations.
0
0
u* v*

Boundary conditions

Three types of boundary conditions have been considered to the present problem. They are as follows,
x*
y*
————– (9)

At the walls: No slip condition is
used, i.e. u* 0, v* 0 .
u* u *
*
*
*
*
p*
* H *
1 1 u* 1 u*

At the inlet: Uniform flow condition
u v

M n H
x*
y*
x*
x*
Re x* x* y* y*
is specified and the transverse
———— (10)
velocity is set to zero, i.e., u* 1, v* 0 .


At the exit: Fully developed
v* v *
*
*
*
*
p*
* H 1 1 v 1 v
*
*
* *
condition is assumed and hence
u v
Mn H
x*
y*
y*
y* Re x* x* y* y*
gradients is set to zero,
————– (11)
i.e., u*
x* 0 , v*
x* 0 .
Where Mn is the Magnetic number arising from FHD which is given by,
n 0 0
n 0 0
M H 2 /U 2
————– (12)
Another nondimensional parameter which will affect the biomagnetic fluid flow under consideration is Reynolds number, which is given by,

NUMERICAL SIMULATION
The nondimensional partial differential continuity and momentum equations (9)(11) have been solved according to the SIMPLE method in the finite volume formalism by use of a non uniform and staggered grid in both coordinating directions allowing higher grid node concentrations in the region close to the walls and close to the higher rate of change of magnetic field. The convection and diffusion terms have been discretized with the help of Power law scheme [8]. The
Re UW
————– (13)
discretized equations have been solved iteratively by using linebyline ADI (Alternating directional implicit) method. For all calculations, the length of the
channel is considered to be 220. During computation, the numerical mesh is considered to be comprising of 160 Ã— 60 grid nodes in x and y direction respectively.

RESULT AND DISCUSSION
During our study, a series of numerical simulation is performed and the parameters those affect the flow characteristics, are identified as,

Reynolds number, Re = 100.

Magnetic number arising from fero hydrodynamics (FHD),
2Ã—104 Mn 5Ã—104.

Uniform velocity at inlet and fully develop at exit.

The length of the magnet ranging from 4 nondimensional units to 22 non dimensional
units and magnet is placed at the middle of computational domain keeping equal unit of nondimensional length in each side of the middle point of the channel.

Variation of streamline and Velocity contours
In all the cases, we have shown the variation of streamline contours and velocity contours of the flow subjected to axial magnetic field. Fig.2 and Fig.3 show the streamline and velocity contour of the channel flow without considering magnetic number for a typical value of Re of 100. Here all the streamlines are straight, parallel and perpendicular to the cross section of the channel as there is no external force acting on the fluid particle during its motion along the streamline. The velocity profile is uniform at the inlet of the channel and fully develop after covering a short distance and also remain fully develop up to the exit of the channel. Fig.4.1. and Fig.4.2.show the effect of length of magnet on streamline and
velocity contour for a typical value of Re of 100 and Magnetic number of 2Ã—104. Fig.5.1 and Fig.5.2 show the effect of magnetic
number on streamline and velocity contour for a typical value of Re of 100 and magnet of 20 unit of nondimensional length. In all the cases, separation phenomenon occurs just before the entrance of the magnetic field and reattaches just after the magnetic field. It is also observed that the size of the recirculation zone increases along the axial direction and decreases along transverse direction with the increase in the length of the Magnet and Magnetic number as the imaginary external magnetic lines of force increases with increasing length of the magnet and magnetic number. This magnetic force, along the axial direction, acts on the magnetic fluid particle which is flowing through the channel. From fig. 4.1, 4.2, 5.1 and 5.2, it has been also observed that the complete recirculating zone has formed at that zone where the velocity profile developed a small velocity loop, which may be a positive loop or a negative loop i.e. the complete recirculating zone has formed at that zone where the velocity suddenly changes from negative to positive and again from positive to negative and vice versa. Moreover, with the increase in recirculating zone, it can be expected that the rate of heat transfer as well as temperature of the working fluid will increase mainly due to conversion of kinetic energy of the working fluid in to heat energy. Due to the increasing recirculating zone, retention of blood or magnetic fluid at a particular area will increase as a result of which, toxide medicine, if injected at that location, will get time to be sedimented on the required area, without affecting the important part of the body like heart, kidney, liver, lung etc. Recirculating bubble may create turbulence in the fluid flow which is used for proper mixing in combustor, carburetor etc.
1
0.9
0.8
0.7
0.6
0.5
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.2. Streamline contour without Considering Magnetic number at Re of 100
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0.4
0.3
0.2 0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.4.1.5. Mn20000,
Re100, Magnet Length22unit
Fig.4.1. Effect of Length of Magnet on streamline contour for
fixed Magnetic Number and Reynold Number
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.3. Velocity profile without Considering Magnetic number at Re of 100
1
0.9
0.8
0.7
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.4.2.1. Mn20000,
Re100, Magnet Length4unit
Fig.4.1.1. Mn20000,
Re100, Magnet Length4unit
1
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 22
Fig.4.1.2. Mn20000,
Re100, Magnet Length6unit
1
0.9 0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.4.1.3. Mn20000,
Re100, Magnet Length16unit
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.4
0.3
0.2
0.1
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.4.2.2. Mn20000,
Re100, Magnet Length6unit
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.4.2.3. Mn20000,
Re100, Magnet Length16unit
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.4.1.4. Mn20000,
Re100, Magnet Length20unit
0
20.24 2 50 95 110 125 165 195 215
X*
Fig.4.2.4. Mn20000,
Re100, Magnet Length20unit
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.4.2.5. Mn20000,
Re100, Magnet Length22unit
Fig.4.2. Effect of Length
of Magnet on Velocity Contour for Fixed Magnetic and Reynold Number
1
0.9
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.1
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.5.2.1. Considering
Magnetic Number Mn2Ã—10
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.8
0.7
0.2
0.1
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.5.1.1. Considering Magnetic Number Mn 2Ã—104
1
0.9
0.8
0
20.24 20 50 95 110 125 165 195 215
X*
Fig.5.2.2. Considering Magnetic Number Mn3Ã—104
1
0.9
0.8
0.7
0.6
Y*
Y*
0.5
0.4
0.3
0.2
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
0.1
0
20.24 50 95 110 125 165 195 215
X*
Fig.5.2.3. Considering Magnetic Number Mn4Ã—104
Fig.5.1.2. Considering Magnetic Number Mn3Ã—104
1
1
0.9
0.8
0.7
0.6
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.5
Y*
Y*
0.4
0.3
0.2
0.1
0
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.5.1.3. Considering Magnetic Number Mn4Ã—104
1
0.9
0.8
20.24 20 50 95 110 125 165 195 215
X*
Fig.5.2.4. Considering Magnetic Number Mn5Ã—104
Fig.5.2. Effect of Magnetic number on Velocity profile at Re of 100 and magnet length 20 unit
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 20 40 60 80 100 120 140 160 180 200 220
Fig.5.1.4. Considering Magnetic Number Mn5Ã—104
Fig.5.1. Effect of Magnetic number on the streamline contour for a fixed Reynold number

Variation of average stagnation pressure along the length of the channel
Fig. 6.1 show the variation average stagnation pressure without considering magnetic number. From the Figure it is cleared that average stagnation pressure gradually decrease along the length of the
channel due to friction. Fig.6.2 show the
700000
effect of length of magnet on average
600000
stagnation pressure for a typical value of Re
500000
of 100 and Magnetic number of 2Ã—104.
400000
Fig.6.3 show the effect of magnetic number on average stagnation pressure for a typical value of Re of 100 and magnet of 20 unit of nondimensional length. From the Fig.6.2 and 6.3 it is clear that average stagnation pressure increases with increase in magnetic
number as well as with the increase in length
Mn60000 Mn40000
Mn30000 Mn20000
Mn60000 Mn40000
Mn30000 Mn20000
300000
p
p
sav
sav
200000
100000
0
90 95 100 105 110 115 120 125
90 95 100 105 110 115 120 125
700
650
600
550
500
450
400
350
300
700
650
600
550
500
450
400
350
300
0 50 100 150 200 250
x*
Fig.6.3. considering
of the magnet due to formation of recirculating bubble from just after the inlet to just before the outlet of the channel. But there is a sharp rise of average stagnation pressure at the magnetic zone. It has also been observed that as we increase the magnetic number and length of the magnet, its impact has been started earlier on the magnetic fluid flow. The generation of high pressure at that magnetic zone may compress the stenosis, formed inside the artery which will lead to opening of the blockage in the artery.
5
0
5
av
av
10
p
p
15
20
25
30
0 50 100 150 200 250
*
*
Fig.6.1. AveraXge stagnation pressure without considering
Magnetic Number and considering Re100
22unit 20unit
18unit 16unit
22unit 20unit
18unit 16unit
700000
700
650
600
550
500
450
400
350
300
700
650
600
550
500
450
400
350
300
600000
500000
400000
p
p
sav
sav
300000
90 95 100 105 110 115 120 125
90 95 100 105 110 115 120 125
200000
100000
0
0 50 100 150 200 250
x*
Fig.6.2. considering Re100 and Mn20000
Re100 and Magnet length 20 unit
Fig.6. Effect of Length of Magnet and Magnetic number on average stagnation pressure



CONCLUSION
The shape and size of the bubble and average stagnation pressure having greater importance in case of various practical applications, such as magnetic hyperthermia and reduction of bleeding, the development of magnetic devices for cell separation, targeted transport of drugs using magnetic particles as drug carriers, magnetic wound or cancer tumour treatment etc, as magnet can control the magnetic particle. It is expected that the rate of heat transfer as well as temperature of the flowing fluid will increase as the recirculating zone increases. It occurs mainly due to conversion of kinetic energy of the flowing fluid in to heat energy. Due to the increasing recirculating zone, retention of blood or magnetic fluid at a particular area will increase as a result of which, toxide medicine, if injected at that location, will get time to be sedimented on the required area, without affecting the important part of the body like heart, kidney, liver, lung etc. Recirculating bubble may create turbulence in the fluid flow which is used for proper mixing in combustor, carburetor etc. The generation of high pressure at that magnetic zone may compress the stenosis, formed inside the artery which will lead to opening of the blockage in the artery.
NOMENCLATURE
H Magnetic field intensity, Amp.m1 H0 Magnetic field intensity, Amp.m1 at the middle of the channel
L Length of the channel, m M Magnetization, Amp. m1
Mn Magnetic Number arising from FHD P Static pressure, Nm2
Re Reynolds Number

Velocity in xdirection, ms1 U Average velocity, ms1

Velocity in ydirection, ms1
W Height of the channel, m x, y Cartesian coordinates
Density, kg m3
Dynamic viscosity, kg m1s1

Md. A. Ikbal, S. Chakraborty, K. K. L.Wong,
J. Mazumder, P. K. Mandal, 2009, Unsteady Response of NonNewtonian Blood Flow Through a Stenosed Artery in Magnetic Field,
J. Computational and Applied Mathematics, 230, 245259.

E.E. Tzirtzilakis, 2008, Biomagnetic Fluid Flow in a Channel with Stenosis, J. Physica, 237(D), 6681.

P. K. Papadopoulos, E. E. Tzirtzilakis, 2004, Biomagnetic flow in a curved square duct under the influence of an applied magnetic field
J. Fluid Physics, 16, 29522962.

S. V. Patankar, 1980, Numerical heat transfer and fluid flow, Hemisphere Publication.
Âµ0 Magnetic permeability of vacuum, Tesla. m. Amp.1
Magnetic susceptibility
Superscripts
* Dimensionless terms
REFERENCES:

V. C. Loukopoulos, E. E. Tzirtzilakis, 2004, Biomagnetic Channel Flow in Spatially Varying Magnetic Field, Int. J. Engineering Science, 42, 571590.

E. E. Tzirtzilakis, 2006, Mathematical Modeling and Simulation of Blood Flow in Magnetic Field, 2nd Int. Conference From Scientific Computing to Computational Engineering.

E. E. Tzirtzilakis, M. Xenos, V. C. Loukopoulos, N. G. Kafpussias, 2006, Turbulent Biomagnetic Fluid Flow in a Rectangular Channel under the action of a Localized Magnetic Field, Int. J. Engineering Science, 44, 12051224.

T. Hayat, A. Afsar, M. Khan, S. Asghar, 2007, Peristaltic Transport of a Third Order Fluid under the Effect of a Magnetic Field, Computers and Mathematics with Applications, 53, 10741087.