American Journal of Applied Mathematics
Volume 4, Issue 1, February 2016, Pages: 53-61

Effects of Temperature Dependent Viscosity on Magnetohydrodynamic Natural Convection Flow Past an Isothermal Sphere

Mwangi Wanjiku Lucy*, Mathew Ngugi Kinyanjui, Surindar Mohan Uppal

Department of Pure and Applied Mathematics, Jomo Kenyatta University of Agriculture and Technology, Nairobi, Kenya

Email address:

(L. W. Mwangi)

To cite this article:

Mwangi Wanjiku Lucy, Mathew Ngugi Kinyanjui, Surindar Mohan Uppal. Effects of Temperature Dependent Viscosity on Magnetohydrodynamic Natural Convection Flow Past an Isothermal Sphere. American Journal of Applied Mathematics. Vol. 4, No. 1, 2016, pp. 53-61. doi: 10.11648/j.ajam.20160401.15


Abstract: In this study, the effects of temperature dependent viscosity on MHD natural convection flow past an isothermal sphere are determined. The uniformly heated sphere is immersed in a viscous and incompressible fluid where viscosity of the fluid is taken as a non-linear function of temperature. The Partial Differential Equations governing the flow are transformed into non dimensional form and solved using the Direct Numerical Scheme and implemented in MATLAB. The numerical results obtained are presented graphically and in tables and are discussed. In this study, it has been observed that increasing the Magnetic parameter M leads to decrease in velocity, temperature, skin friction and the rate of heat transfer. It has also been noted that increase in the Grashof number Gr leads to increase in velocity and temperature whereas increase in the values of eta leads to increase in temperature but there is a decrease in velocity. These results are applicable to engineers in designing electricity plants which have higher life expectancy.

Keywords: Natural Convection Flow, Temperature Dependent Viscosity, Magnetohydrodynamic (MHD), Isothermal Sphere


1. Introduction

Electrically conducting fluids in presence of strong magnetic field is an important phenomenon in our day-to-day lives. The interaction of the current with the magnetic field changes the motion of the fluid and produces an induced magnetic field. An isothermal sphere in this study is made of non-conducting material and is uniformly heated. Molla et. al (2005) studied the magnetohydrodynamic natural convection flow from a sphere with uniform heat flux in the presence of heat generation. They observed that for increased values of heat generation parameter, there was an increase in local skin friction, velocity and temperature profiles but there was a decrease in the local rate of heat transfer. Temperature distribution was observed to increase whereas velocity distribution, local rate of heat transfer and local skin friction coefficient were observed to decrease slightly as the Magnetic parameter  increased. The MHD natural convection flow from an isothermal sphere has many applications in the engineering processes. This has been due to its cost and practical advantages which has led to its usefulness. For example, the fluid flow is considered to be of importance in the geothermal and hydroelectric plants in Kenya.

Temperature is one of the factors that may cause variation in the viscosity of a given fluid. Soares et al (2010) investigated the effects of temperature dependent viscosity on forced convection heat transfer from a cylinder in cross flow of power-law fluids and found out that the variation of viscosity with temperature is shown to have an effect on both the local and the surface averaged values of the Nusselt number.

It has been found out that the velocity and temperature of the flow change more or less with the variation of the flow parameters and that for higher values of Prandtl number Pr, both velocity and temperature decreases such that there exists a local maximum value of velocities (Kabir et. al, 2013). Alam et. al (2007) investigated the effects of viscous dissipation on MHD natural convection flow over a sphere in the presence of heat generation and found out that velocity increases as the values of viscous parameter increases. Also, they found out that velocity distribution increase as the values of heat generation parameter increases.

Shrama et. al(2010) studied the effect of temperature dependent electrical conductivity on steady natural convection flow of viscous incompressible low Prandtl (Pr<<1) electrical conducting fluid along an isothermal vertical plate in the presence of transverse magnetic field and exponentially decaying heat generation. They found out that fluid velocity increases in the presence of heat generation due to increase in the electrical conductivity parameter while it decreases due to increase in the magnetic field intensity. They also deduced that fluid temperature increases in the presence of volumetric rate of heat generation or due to increase in magnetic field intensity while it decreases due to increase in electrical conductivity parameter. They also explained that the skin friction coefficient increases with the increase in electrical conductivity parameter or in the presence of volumetric rate of heat generation while it decreases due to increase in the Prandtl number. Finally, they concluded that the rate of heat transfer increases with the increase in the Prandtl number of electrical conductivity parameter which is not much effective while it decreases due to increase in the magnetic field intensity.

Chaudhary et. al (2007) studied combined heat and mass transfer effects on MHD free convection flow past an oscillating plate embedded in porous medium and found out that the concentration decreases with an increase in Schmidt number. They found out that in case of cooling of the plate (Gr >0), the velocity decreases with an increase in the phase angle, magnetic parameter, Schmidt number and Prandtl number while it increases with an increase in the value of Grashof number and modified Grashof number, permeability parameter and time. In case of heating of the plate (Gr< 0), the velocity increases with an increase in magnetic parameter, Schmidt number and Prandtl number while it decreases with an increase in the value of Grashof number and modified Grashof number, permeability parameter and time. The Skin friction increases with an increase in Schmidt number, Prandtl number and Magnetic parameter while it decreases with an increase in the value of Grashof number, modified Grashof number, permeability parameter and time. They also concluded that Nusselt number increase with an increase in the Prandtl number while temperature decreases with an increase in the value of Prandtl number.

Molla et. al (2005) studied Magnetohydrodynamic natural convection flow from an isothermal sphere with temperature dependent heat generation and they found out that when the values of heat generation parameter Q increases, there is an increase in the local skin friction coefficient Cfx but there is a decrease in the local rate of heat transfer. Both the velocity and temperature profiles increase significantly when the value of heat generation parameter increases. They also found out that the local rate of heat transfer and the local skin friction coefficient decreases when the value of magnetic parameter increases. Finally, they concluded that increased values of magnetic parameter leads to decrease in the velocity distribution whereas there is an increase in the temperature distribution.

Natural convection flow along an isothermal vertical plate with temperature dependent viscosity and heat generation has been studied by Molla et. al (2014) and they deduced that the effect of viscosity variation parameter and Rayleigh number decreases the skin friction coefficient whereas increasing the local average rate of heat transfer. They also found out that the momentum and thermal boundary layer becomes thinner when the values of viscosity -variation parameter increases. They have also found out that both viscosity and velocity distribution increases with the effect of Rayleigh number. This has also led to significant decrease in temperature distributions whereas the thickness of momentum boundary layer is enhanced. They concluded that the decrease in the local and average rate of heat transfer is more significant.

Miraj et. al (2011) studied the effects of viscous dissipation and radiation on MHD free convection flow along a sphere with joule heating and heat generation and they found out that velocity profiles increases with the increasing values of radiation parameter, heat generation parameter, magnetic parameter, joule heating parameter and viscous dissipation parameter whereas the profiles decreases with increasing values of Prandtl number. They also discovered that the skin friction coefficient increases with the increasing values of heat generation parameter whereas decreases for the increasing values of magnetic parameter. Finally, they found out that the rate of heat transfer increases for increasing values of radiation parameter whereas decreases for increasing values of heat generation parameter, joule heating parameter and viscous dissipation parameter.

Molla et al (2012) have studied the MHD natural convection flow from an isothermal cylinder under the consideration of temperature dependent viscosity and they found out that increasing the values of magnetic parameter M and viscosity variation parameter leads to decrease in the local skin friction coefficient, Grashof number and the local Nusselt number. It was also seen that the velocity distribution decreases as well as the temperature distribution increases with the increasing values of the magnetic parameter M and viscosity variation parameter.

Hague et. al (2014) studied the effects of viscous dissipation on MHD natural convection flow over a sphere with temperature dependent thermal conductivity in presence of heat generation and they found out that the velocity and temperature of the fluid within the boundary layer increases with increasing thermal conductivity variation parameter, heat generation parameter and viscous dissipation parameter. They also noted that the skin friction along the surface of the sphere increases with increasing thermal conductivity variation parameter, heat generation parameter and viscous dissipation parameter but decreases for the increasing values of M. The heat transfer rate from the surface was found to decrease with the increasing value of magnetic parameter, thermal conductivity variation parameter, heat generation parameter and viscous dissipation parameter.

This study, presents the findings of the studies made on the effects of temperature dependent viscosity on the MHD natural convection flow past an isothermal sphere taking viscosity as a non-linear function of temperature and analyzing the results using Direct Numerical Scheme to a simple analytical case.

2. Formulation

The configuration of the problem considered in this study is as shown in the diagram below:

Figure 1. Physical Model of Problem.

The diagram above shows a two-dimensional MHD Laminar free convective flow past a uniformly heated sphere Center  and radius   which is immersed in a viscous and incompressible fluid having a temperature dependent viscosity. In this study, viscosity is considered as a non-linear function of temperature.

2.1. Governing Equations

2.1.1. Equation of Conservation of Mass

This equation states that mass can neither be created nor destroyed under normal circumstances. It is obtained from the fact that the mass of fluid entering and leaving a volume in the flow field have the same mass balance.

(2.1)

Where  is the velocity in the X, Y direction and

(2.2)

2.1.2. Equation of Motion

The expression of the equation in tensor form is given as;

(2.3)

Where is the temporal acceleration,is the convective acceleration, is the pressure gradient, is the force due to viscosity and  represents the body forces vector in X and Y directions.

2.1.3. Equation of Induced Magnetic Force

The Induced magnetic force for the magnetic fluid is given as:

(2.4)

Where B is the Magnetic Field Density vector, Mg is magnetization of the magnetic material and represents the Induced free current.

2.1.4. The Energy Equation

The general equation is given as:

(2.5)

Where viscous energy dissipation term  is defined as:

(2.6)

2.2. Non-dimensional Numbers

Non-dimensional numbers are applicable in this study in order to transform the obtained results obtained in the model into any other dynamically similar case. In this work, the following numbers have been used;

2.2.1. Grashof Number, Gr

This number usually occurs in natural convection problem and is usually defined as:

2.2.2. Prandtl Number, Pr

This number gives the ratio of viscous force to the thermal and is defined as:

2.2.3. Nusselt Number, Nu

This number is important in the heat transfer problems and is defined as;

2.2.4. Magnetic Parameter, M

This is defined as;

2.2.5. Skin Friction Coefficient, Cf

This is defined as;

2.2.6. Pressure Coefficient, Cp

This is the ratio of pressure number to the inertia force and it gives the importance of pressure force to the inertia force. It is given as;

2.3. Problem Modelling

In this study, the fluid viscosity  is assumed to vary as an inverse function of temperature T, in the form;

(2.7)

Where is the temperature of the fluid andis the temperature of the ambient fluid.

The non-dimensional variables used in transforming equations (2.1), (2.3),(2.4) and (2.5) into non- dimensional form are:

,,,

,,,

(2.8)

Equations (2.8) are used together with the non-dimensional numbers given above.

Equation of conservation of mass reduces to:

(2.9)

Substituting equations (2.8) in Equation (2.1), we obtain;

(2.10)

The expression of the equation of motion reduces to;

(2.11)

Substituting equations (2.7), (2.8) and the non-dimensional numbers above in (2.11) and simplifying, we obtain;

(2.12)

(2.13)

Thus, the non-dimensionalized equation of motion becomes:

(2.14)

On simplification of the equation of induced magnetic field, we obtain the Induced Force which is taken to be along the x-axis in our study and is given as:

(2.15)

This result is used in the last term of equation (2.3).

In energy equation (2.5), the flow viscosity is not large; hence the Viscous Energy dissipation  is small and is neglected.

 is considered to be along the y-axis and thus the energy equation becomes;

(2.16)

Substituting equations (2.8) in equation (2.16) and simplifying, we obtain;

(2.17)

Taking, therefore equation (2.17) becomes;

(2.18)

Boundary Conditions

u = v = 0, q = 1, at y = 0       (2.19a)

u ® 0, q ® 0, as y ® ¥      (2.19b)

2.4. Method of Solution

Direct Numerical Scheme (DNS) is applied in solving equations (2.10), (2.14) and (2.18) subject to the boundary conditions (2.19). The following new set of transformations is introduced;

(2.20)

(2.21)

Thus, substituting equations (2.20) and (2.21) into equation (2.10), we get;

(2.22)

Substituting transformations (2.20) in equation (2.14), we get;

(2.23)

Substituting transformations (2.20) in equation (2.18), we get;

(2.24)

Boundary Conditions;

,  at  any             (2.25a)

,  at ,             (2.25b)

, as ,             (2.25c)

Writing equations (2.22), (2.23) and (2.24) in terms of central finite differences, we get the following equations;

(2.26)

Making  the subject of the formulae, we get;

(2.27)

Writing equation (2.23) in terms of finite differences, we get;

(2.28)

Consider the term;

;

This can be written as;

(2.29)

Putting this in equation (2.28) and making the subject of the formulae, we get;

(2.30)

Writing equation (2.24) in terms finite differences, we get;

(2.31)

From equation (2.31), consider;

(2.32)

Putting equation (2.32) in (2.31) and making the subject of the formulae, we get;

(2.33)

The physical quantities to be obtained are the shearing stress (the rate of skin friction) and the rate of heat transfer. The equations to be used in obtaining the graphs are;

(2.34)

(2.35)

Writing the equation in terms of finite differences, we have;

(2.36)

(2.37)

Equations (2.27), (2.30), (2.33), (2.36) and (2.37) are the final equations which are solved using MATLAB in order to obtain the results shown below.

3. Results and Discussion

In this study, we have investigated effects of temperature dependent viscosity on MHD natural convection flow past an isothermal sphere. Here, we have considered the viscosity of the fluid as a non-linear function of temperature. The numerical solutions start from the lower stagnation point , round the sphere to the upper stagnation point where. The MHD Parameter values are taken as  and the variable eta and Grashof number is given as.

Figure 2. Velocity varying Magnetic Parameter M.

Figure 3. Temperature varying Magnetic Parameter M.

From figure 2 above, velocity varying MHD parameter is presented as shown. Setting . It can be seen that, as the magnetic parameter M increases the velocity profile decrease. From figure 3, as the magnetic parameter increases, the temperature profile decreases slightly. This is because the interaction of the magnetic field and the moving electric charge carried by the fluid induces a force which tends to oppose the fluid motion. Velocity increases and then decreases slowly and finally approaches to zero near the surface of the sphere. This shows that there is a maximum local velocity within the boundary layer.

Figure 4. Velocity varying Grashof number Gr.

Figure 5. Temperature varying Grashof number Gr.

Figure 4 and figure 5 shows the Velocity and Temperature profiles respectively. , , ,.

It can be seen that increase in Grashof number  leads to increase in the velocity (figure 4) as well increase in  leads to increase in temperature as in figure 5. This is because as temperature of fluid increases, viscosity of the fluid decreases and therefore, the drag force in the fluid is reduced and thus the velocity of the fluid increases.

 

Figure 6. Velocity varying eta.

Figure 7. Temperature varying eta.

Figure 6 and figure 7 shows the velocity and temperature distribution against the variable Y for different values of etawhile, , . It can be observed that the velocity decreases and temperature distribution increases with the increasing values of the eta, . It can be deduced that at each value of eta, the velocity profile has a local maximum value within the boundary layer.

Figure 8. Skin friction varying Magnetic parameter M.

Figure 8 represents skin friction varying Magnetic parameter whereas figure 9 represents Rate of Heat Transfer varying the Magnetic Parameter.. In this case the values are taken at the stagnation point of the sphere, . It can be observed that increasing the value leads to decrease in skin friction and rate of Heat transfer. This is due to the fact increasing the value of  leads to increase in Lorentz force which opposes the motion hence decreasing the velocity and temperature gradients thus causing a decrease in Skin Friction as well as the Rate of heat transfer. In figure 9, it can be observed that there is an inverse behavior of the curves. This can be explained from the fact that viscosity is taken to vary inversely proportional to temperature thus a varying change of the curves.

Figure 9. Heat Transfer varying Magnetic Parameter M.

Table 1. The results of skin friction and Rate of heat Transfer (Nusselt number) varying Magnetic parameter while Pr=0.073, eta=0.015, Gr=4.

M=0 M=1 M=5 M=10
X0 Skin friction Nusselt number X0 Skin friction Nusselt number X0 Skin friction Nusselt number X0 Skin friction Nusselt number
0.02 0.4015 1.9581 0.02 0.3914 1.9498 0.02 0.3558 1.9202 0.02 0.3198 1.8890
0.40 8.0984 1.8704 0.40 7.9050 1.8642 0.40 7.2220 1.8417 0.40 6.5319 1.8181
0.79 15.4196 1.7279 0.79 15.0242 1.7250 0.79 13.6387 1.7143 0.79 12.2564 1.7029
1.17 21.3508 1.5315 1.17 20.7584 1.5336 1.17 18.7002 1.5409 1.17 16.6756 1.5481
1.55 24.7245 1.2949 1.55 23.9855 1.3039 1.55 21.4388 1.3359 1.55 18.9665 1.3679
1.94 24.2875 1.0665 1.94 23.5241 1.0830 1.94 20.9063 1.1409 1.94 18.3863 1.1987
2.32 19.3930 0.9343 2.32 18.7799 0.9546 2.32 16.6764 1.0262 2.32 14.6496 1.0972
2.71 10.9739 0.9644 2.71 10.6430 0.9816 2.71 9.4990 1.0428 2.71 8.3833 1.1045
3.09 1.2630 1.1182 3.09 1.2298 1.1281 3.09 1.1125 1.1642 3.09 0.9942 1.2023

4. Conclusions

In this study, the effects of viscosity which is a non-linear function of temperature on MHD natural convection flow past an isothermal sphere has been determined. Numerical results of the study have been obtained using the Direct Numerical Scheme. From this study, the following conclusions can be made;

Increase in Magnetic parameter leads to decrease in velocity profile as well as the temperature profile.

Increase in the Grashof number leads to increase in velocity and temperature.

Increase in the values of eta  leads to increase in temperature and decrease in velocity.

Increase in the value of leads to a decrease in skin friction and the rate of heat transfer

Acknowledgement

I would like to thank my Uncle Joseph Kimama, Auntie Grace Wambui and Marie Dawood Trust Foundation for their financial support.

Nomenclature

Roman Symbols  

Radius of the sphere, m

Specific heat at constant pressure, J/deg kg

Skin friction coefficient

Dimensionless stream function

Acceleration due to gravity, g

Grashof number

Thermal conductivity of the fluid, W/mK

MHD parameter

Nusselt number

Prandtl number

Heat flux at the surface, W/m2

Temperature of the fluid, K

Temperature of the ambient fluid, K

Temperature at the surface, K

Dimensionless velocity component

Axis direction

Fluid velocity in the x, y direction, m3/s
Body forces vector in x and y directions, N

Magnetic field intensity, T

Local skin Friction coefficient

Magnetization of the magnetic Material

Heat Generation Parameter

Radial distance from the symmetrical axis to the surface of the sphere, m

The center of the sphere
Greek Symbols Quantity

Volumetric coefficient of thermal expansion, (Co)-1

Shearing stress, N/m2

Density of the ambient fluid, Kg/m3

Density of the fluid, Kg/m3

Viscosity of the fluid, Ns/m2

Dimensionless temperature function

Strength of magnetic field, A/m

Electric conduction, S/m

Gradient operator

References

  1. Alam M., Alim M.A, Chodhury M.K.(2007).Viscous Dissipation effects on MHD Natural Convection Flow over a sphere in the presence of Heat Generation. Nonlinear Analysis Modelling and Control. 12(4):447-459.
  2. Miraj M., Alim M.A. and Andallah L.S. (2011).Effects of Viscous Dissipation and Radiation on Magnetohydrodynamic Free Convection Flow along sphere with joule Heating and Heat Generation. Thanmasat International Journal. Of Science and Technology. 16(4):52-64.
  3. Kabir K.H., Alim A. and Andallah L.S (2013).Effects of viscous dissipation on MHD natural convection flow along a vertical wavy surface. Journal of Theoretical and Applied Physics.7(31):4-8.
  4. Soares A.A, Ferreira J.M., Caramelo L., Anacleto J., and Chhabra R.P. (2010).The effects of temperature-dependent Viscosity on Forced Convection heat transfer from a cylinder in crossflow of power-law fluids. International Journal of heat and mass transfer. 53(21-22): 4728-4740.
  5. Morris S. (1982).The effects of a strongly temperature dependent viscosity on a slow flow past a hot sphere. Journal of Fluid Mechanics. 124:1-26.
  6. Molla M., Suvash C. S., and Khan M.A.I. (2012). MHD Natural Convection Flow from an Isothermal Circular Cylinder under the consideration of Temperature Dependent Viscosity. Engineering Computations.29(8):875-887.
  7. Huang M.J., and Chen C.K.(1987).Laminar Free Convection Flow from a Sphere with Blowing and Suction.ASME Journal of Heat transfer.86(4),537-541.
  8. Molla M.M, Hossain M.A. and Taher M.A. (2005). Magnetohydrodynamic natural Convection Flow on a sphere with Uniform Heat flux in presence of Heat Generation. Acta Mechanica. 186:75-86.
  9. Shrama P.R. and Gurminder S. (2010).Steady MHD natural convection flow with variable electrical conductivity and heat generation along an isothermal vertical plate. Tamkang Journal of Science and Engineering.13(3), 235-242.
  10. Chaudhary R.C. and Jain A. (2007).Combined Heat and Mass Transfer Effects on MHD Free Convection Flow Past an Oscillating Plate Embedded in Porous Medium. Rom. Journal. Physics. 52(5-7), 505-524.
  11. Molla M., Taher M.A., Chowdhury M.K. and Hossain A. (2005). Magnetohydrodynamic Natural Convection Flow from an Isothermal Sphere with Temperature Dependent Heat Generation. Nonlinear Analysis: Modelling and Control. 10(4), 349-363.
  12. Molla M.M, Biswas A., Al-Mamun A. and Hossain A. (2014). Natural Convection Flow along an Isothermal Vertical Flat Plate with Temperature Dependent Viscosity and Heat Generation. Journal of Computational Engineering, Volume 2014(2014), Article ID 712147, 13 pages.
  13. Perot J.B. and Subramanian V. (2007). A Discrete Calculus Analysis of the Keller Box Scheme and a Generalization of the Method to Arbitrary Meshes. Journal of Computational Physics. 226(1),494-508.
  14. Haque R., Alam M., Ali M.M and Sheikh N.M.A. (2014).Effects of viscous dissipation on MHD natural convection flow over a sphere with temperature dependent thermal conductivity in presence of heat generation. Journal of Computer and Mathematical Sciences.5(1), (1-122).
  15. Kang'ethe G.(2011). Magnetohydrodynamic Flow in Porous Media over A Stretching Surface in a Rotating System with Heat and Mass Transfer (Doctoral Thesis).Retrieved from http://ir.jkuat.ac.ke/bitstream/handle/123456789/867/Giterere,Kang'ethe_PhD.Applied%20Mathematics-2011.pdf?sequence=1.

Article Tools
  Abstract
  PDF(1167K)
Follow on us
ADDRESS
Science Publishing Group
548 FASHION AVENUE
NEW YORK, NY 10018
U.S.A.
Tel: (001)347-688-8931