MAGNETOHYDRODYNAMICS OF BLASIUS VISCOELASTIC FLUID WITH VISCOUS DISSIPATION AND SUCTION/INJECTION EFFECTS NURATIKAH HANI BINTI ABD RANI A thesis submitted of the requirements for the award of the degree of Master of Science (Engineering Mathematics) Faculty of Science Universiti Teknologi Malaysia JANUARY 2013
22
Embed
MAGNETOHYDRODYNAMICS OF BLASIUS VISCOELASTIC FLUID … · lebih tinggi satu peringkat berbanding dengan persamaan Navier-Stokes atau persamaan lapisan sempadan, maka syarat sempadan
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
MAGNETOHYDRODYNAMICS OF BLASIUS VISCOELASTIC FLUID WITH
VISCOUS DISSIPATION AND SUCTION/INJECTION EFFECTS
NURATIKAH HANI BINTI ABD RANI
A thesis submitted of the
requirements for the award of the degree of
Master of Science (Engineering Mathematics)
Faculty of Science
Universiti Teknologi Malaysia
JANUARY 2013
iii
To my beloved father and mother
Thank you for your support.
iv
ACKNOWLEDGEMENT
Alhamdulillah, I am grateful to Allah S.W.T for His grace and mercy for
giving me time, strength and patience in completing this study. Without His
willingness, I would not be able to complete this study.
I would like to take this opportunity to express my gratitude to my supervisor,
Dr. Sharidan Shafie for all the guidance, knowledge and advice throughout the
completion of my study. All his efforts are very much appreciated. I also would like
to extend my sincere gratitude to Abdul Rahman Mohd Kasim who kindly offered
precious and helpful comments in the preparation of this study.
My appreciation would also go to my beloved family who either directly or
indirectly involved in bringing this successful. Their views and opinions are helpful
indeed.
Last but not least, I would like to express my heartiest appreciation to all
those involved in helping me to complete this study. Unfortunately it is impossible to
list all of them in this limited space. Thank you very much.
v
ABSTRACT
Due to the existence of viscoelastic fluid in technological applications, the
research in the viscoelastic fluids has increase rapidly. In this study, the
magnetohydrodynamics flow for the Blasius viscoelastic fluid along with the effects
of viscous dissipation and suction or injection is considered. Since the equation of
motion in viscoelastic fluid is one order higher than the Navier-Stokes or boundary
layer equations, an extra boundary condition is imposed by augmenting the boundary
condition at infinity. The governing equations are transformed into a non
dimensional boundary layer equation by using non dimensional variables. The
equations are solved numerically by using Keller-box method. Numerical results
consist of the velocity and temperature profiles are presented graphically for assorted
values of magnetic parameter, M, viscoelastic parameter, K, suction or injection
parameter, f w, Prandtl number, Pr and the ratio moving parameter, X . It is found
that, as the values of all parameters increased, the velocity profiles are also increased
but opposite situation occurred in temperature profiles.
vi
ABSTRAK
Oleh kerana wujudnya bendalir likat kenyal di dalam beberapa kegunaan
teknologi, penyelidikan berkenaan bendalir likat kenyal telah meningkat dengan
mendadak. Dalam kajian ini, aliran hidrodinamik magnet untuk bendalir likat kenyal
Blasius berserta dengan kesan pelepasan likat dan sedutan atau suntikan
dipertimbangkan. Disebabkan oleh persamaan gerakan dalam bendalir likat kenyal
lebih tinggi satu peringkat berbanding dengan persamaan Navier-Stokes atau
persamaan lapisan sempadan, maka syarat sempadan tambahan diperlukan dengan
menambah syarat sempadan di infiniti. Persamaan menakluk diubah ke bentuk
persamaan lapisan sempadan tak bermatra, dengan menggunakan pemboleh ubah tak
bermatra. Persamaan diselesaikan secara berangka dengan menggunakan kaedah
Kotak-Keller. Keputusan berangka yang terdiri daripada profil halaju dan suhu
dipersembahkan secara grafik bagi pelbagai nilai parameter magnetik, M, parameter
bendalir likat kenyal, K, parameter sedutan atau suntikan, f w, nombor Prandtl, P r ,
dan parameter nisbah gerakan, X. Didapati bahawa apabila semua nilai parameter
meningkat, profil halaju juga meningkat namun keadaan sebaliknya berlaku bagi
profil suhu.
vii
CHAPTER TITLE PAGE
THESIS STATUS VALIDATION FORM
SUPERVISOR’S DECLARATION
TITLE i
DECLARATION ii
DEDICATION iii
ACKNOWLEDGEMENT iv
ABSTRACT v
ABSTRAK vi
TABLE OF CONTENTS vii
LIST OF TABLE ix
LIST OF FIGURES x
LIST OF SYMBOLS/NOTATIONS xii
1 INTRODUCTION
1.1 Background of the Study 1
1.2 Statement of Problem 3
1.3 Objectives of the Study 3
1.4 Scope of the Study 3
1.5 Methodology of the Study 4
1.6 Significant of the Study 5
TABLE OF CONTENTS
2 LITERATURE REVIEW
2.1 Introduction 7
2.2 Blasius Boundary Layer Equation 7
2.3 Viscoelastic Fluid 9
2.4 Viscous Dissipation 12
2.5 Magnetohydrodynamics (MHD) 13
3 THE DERIVATION OF GOVERNING EQUATION
3.1 Introduction 16
3.2 The Continuity Equation 16
3.3 The Momentum Equation 18
3.4 The Energy Equation 23
3.5 Dimensionless Group 28
4 MATHEMATICAL FORMULATION
4.1 Introduction 29
4.2 The Continuity Equation and Momentum Equation 29
4.3 The Energy Equation 38
5 RESULTS AND DISCUSSION
5.1 Introduction 43
5.2 Results and Discussion 43
6 CONCLUSION
6.1 Introduction 53
6.2 Summary of the Study 53
6.3 Suggestions for Future Study 55
viii
REFERENCES 56
ix
TABLE NO.
5.1
LIST OF TABLE
TITLE
The velocity gradient at the surface f " (0) for X = 0 when K = 0, M = 0, Pr = 1, Ec = 0 and f w = 0
PAGE
45
x
FIGURE NO.
5.1
5.2
5.3
5.4
5.5
5.6
5.7
5.8
LIST OF FIGURES
TITLE
Velocity profiles f '( r ) for various values of M when K = -1, Pr = 1, Ec = 0.5, X = 0, and f w =1 (suction)
Temperature profiles 0 (r ) for various values of M when K = -1, Pr = 1, Ec = 0.5, X=0, and f w = 1 (suction)
Velocity profiles f '( r ) for various values of M when K = -1, Pr = 1, Ec = 0.5, X = 0, and f w = -1 (injection)
Temperature profiles 0 (r ) for various values of M when K = -1, Pr = 1, Ec = 0.5, X=0, and f w = -1 (injection)
Velocity profiles f '( r ) for various values of K when M = 2, Pr = 1, Ec = 0.5, X = 0, and f w = 0 (impermeable plate)
Temperature profiles 0 (r ) for various values of K when M = 2, Pr = 1, Ec = 0.5, X=0, and f w = 0 (impermeable plate)
Velocity profiles f '( r ) for various values of X when M = 2, Pr = 1, Ec = 0.5, K = -1, and f w = 1 (suction)
Temperature profiles 0 (r ) for various values of X when M = 2, Pr = 1, Ec = 0.5, K = -1, and f w = 1 (suction)
PAGE
46
46
47
47
48
48
49
49
5.9 Velocity profiles f '( r ) for various values of X when M = 2, Pr = 1, Ec = 0.5, K = -1, and f w = -1 (injection)
50
xi
510 Temperature profiles 0 (r ) for various values of X when 50 M = 2, Pr = 1, Ec = 0.5, K = -1, and f w = -1 (injection)
511 Temperature profiles 0 (r ) for various values of Pr when 51 M = 2, X= 1, Ec = 0.5, K = -1, and f w = 1 (suction)
512 Temperature profiles 0 (r ) for various values of Pr when 51 M = 2, X= 1, Ec = 0.5, K = -1, f w = -1 (injection)
513 Velocity profiles f '( r ) for various values of f w when 52 M = 2, Pr = 1, Ec = 0.5, K = -1, and X= 0.5
514 Temperature profiles 0 ( r )for various values of f w when 52 M = 2, Pr = 1, Ec = 0.5, K = -1, and X= 0.5
xii
LIST OF SYMBOL/NOTATIONS
A - kinematics tensor
A - kinematics tensor
B0 - magnetic field strength
C - momentum of the mass
Cf - skin friction coefficient
Et - total energy
Ec - Eckert number
F - sum of all forces acting
f - dimensionless stream function
f w - dimensionless injection or suction parameter
I - identify matrix
K - viscoelastic parameter
ko - vortex viscosity
M - magnetic parameter
- local Nusselt number
P - pressure scalar
Pr - Prandtl number
0 - heat flow
Re - Reynolds number
T - temperature
w - velocity vector
xiii
Uw - constant velocity
Ux - dimensional constant
u - velocity component in x-direction
v - velocity component in y-direction
v - fluid velocity
W - work done
x - coordinate in direction of surface motion
y - coordinate in direction normal to surface motion
Greek symbols
a
a 1
a 2
a
a 0
V
0
A
M
v
¥
P
thermal diffusivity
viscoelasticity of fluid
cross-viscosity of fluid
stress tensor
electrical conductivity
dimensionless similarity variable
dimensionless temperature
ratio moving parameter
viscosity coefficient
kinematic viscosity
stream function
fluid density
Subscripts
w - condition at the surface
<x> - condition at ambient medium
xiv
Superscripts
differentiation with respect to V
CHAPTER 1
INTRODUCTION
1.1 Background of the Study
There are two types of fluid which are Newtonian fluid and non-Newtonian
fluid. Newtonian fluid is a type of fluid where the relation between shear stress and
shear rate is proportionally linear. On the other hand, non-Newtonian fluid is a type
of fluid where the relation between shear stress and shear rate is not a linear function
of the spatial variation of velocity at a given temperature and pressure. Non-
Newtonian fluids fail to obey Newton’s law of friction since non Newtonian fluids
form very wide class of different materials, whose only common features are fluidity.
It is now generally recognized that non-Newtonian fluids are more appropriate than
Newtonian fluids in industrial and technological applications and numerous models
were suggested for non-Newtonian fluids with their constitutive equations varying
greatly in complexity (Othman, 2010).
Viscoelastic fluid is a fluid that returns to its original shape either fully or
partially after the applied stress is released. It is also the subject that explains both
the elastic and viscous behaviour of materials. Viscosity describes how powerfully
fluids resist rapid changes in shape. For example, honey. Since it tends to deform
slowly, thus, it is very viscous. When honey is poured from a jar, it moves very
slowly. This is because the internal stresses between molecules increase with relative
velocity between molecules. As the movement of the molecules increase, the
2
resistance to that movement become greater. Viscous materials display a time-
dependent reaction to deformation.
Elastic deformation can be observed in the way a rubber band acts. Besides
rubber bands, elastic deformation can also be observed in springs. The harder a
spring is compressed, the more force is needed to hold it there. Springs are called
linear because if doubling the amount of compression it will double the amount of
force. Though it may not be obviously to the eye, when forces are applied to metals,
they get linearly compressed or stretched out. Elastic materials quickly return to their
original sizes when all forces are removed. Fundamentally, when changing shape, all
materials show some elastic and viscous effects.
Magnetohydrodynamics (MHD) is the study of the interaction of electrically
conducting fluids and electromagnetic forces. The word magnetohydrodynamics is
from the words magneto which means magnetic field, hydro means liquid, and
dynamics means movement. Synonyms of MHD that are less commonly used are the
terms hydromagnetics and magnetofluiddynamics MHD problems arise in a broad
range of situations. The description of MHD flows connects both the equations of
fluid dynamics which is the Navier-Stokes equations and the equations of
electrodynamics which is Maxwell's equations. These equations are mutually
coupled through the Lorentz force and Ohm's law for moving electrical conductors
and these differential equations have to be solved simultaneously, either analytically
or numerically.
Since vicoelastic fluid grabs our attention, thus in this study, problem in
viscoelastic fluid will be solved. To make it more interesting, MHD flow will be
considered in this study.
3
1.2 Statement of Problem
Due to the existence of viscoelastic fluid in many of technological
applications over the past decades, the interest in the viscoelastic fluids has increase
largely. For that reason, this study will investigate the problem on the viscoelastic
model concentrating on classical Blasius problem over flat plate. Besides, the study
also will explore how the boundary layer flow and heat transfer nature execute with
viscous dissipation included into the energy equation. Further, the effects of
viscoelastic parameter and the ration of moving parameter on the skin friction and
heat transfer coefficient with the presence of magnetohydrodynamics (MHD) flow
will be studied completely.
1.3 Objectives of the Study
The objectives of this study are:
1. To carry out the mathematical formulation of governing equations of the
boundary layer flow of a viscoelastic fluid with the effects of
magnetohydrodynamics (MHD) flow.
2. To transform the governing equations into non-dimension equations.
3. To observe the effects of magnetohydrodynamics (MHD) flow to the fluid
flow characteristic with the present of viscous dissipation and suction or
injection effects.
1.4 Scope of the Study
This study will consider two-dimensional viscoelastic fluid model in
Cartesian coordinate and the fluid is assumed to be incompressible. This study will
consider the MHD boundary layer flow with suction or injection, and viscous
dissipation.
4
1.5 Methodology of the Study
1.5.1 Mathematical Modelling - Problem Formulation
The governing boundary layer flow equations will be acquired for the new