** American Journal of Computational Mathematics ** Vol. 3 No. 2 (2013) , Article ID: 33534 , 8 pages DOI:10.4236/ajcm.2013.32015

Effects of Viscous Dissipation on MHD Natural Convection Flow along a Vertical Wavy Surface with Heat Generation

^{1}Department of Mathematics, Mohammadpur Kendriya University College, Dhaka, Bangladesh

^{2}Department of Mathematics, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh

^{3}Department of Mathematics, Jahangirnagar University, Savar, Bangladesh

Email: kfzkabir@gmail.com, a0alim@gmail.com, andallahls@gmail.com

Copyright © 2013 K. H. Kabir et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received February 16, 2013; revised March 20, 2013; accepted April 17, 2013

**Keywords:** Natural Convection; Wavy Surface; Magnetohydrodynamics; Viscous Dissipation; Prandtl Number; Heat Generation

ABSTRACT

The effect of external magnetic field and internal heat generation or absorption on a steady two-dimensional natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface has been investigated. The governing boundary layer equations are first transformed into a non-dimensional form using suitable set of dimensionless variables. The transformed boundary layer equations are solved numerically using the implicit finite difference method, known as Keller-box scheme. Numerical results for velocity, temperature, skin friction, the rate of heat transfer are obtained for different values of the selected parameters, such as viscous dissipation parameter (Vd), heat generation parameter (Q), magnetic parameter (M) and presented graphically and discussed. Streamlines and isotherms are presented for selected values of heat generation parameter and explained.

1. Introduction

The natural convection boundary layer flow about a heated vertical wavy surface has received a great deal of attention due to its relation to practical applications of complex geometries. There is also a model problem for the investigation of heat transfer from roughened surfaces in order to understand heat transfer enhancement. The natural convection along a vertical wavy surface was first studied by Yao [1] and used an extended Prantdl’s transposition theorem and a finite-difference scheme. He proposed a simple transformation to study the natural convection heat transfer from isothermal vertical wavy surfaces, such as sinusoidal surface. Moulic and Yao [2] also investigated mixed convection heat transfer along a vertical wavy surface. Alam et al. [3] have also studied the problem of free convection from a wavy vertical surface in presence of a transverse magnetic field. The combined effects of thermal and mass diffusion on the natural convection flow of a viscous incompressible fluid along a vertical wavy surface have been investigated by Hossain and Rees [4]. Wang and Chen [5] investigated transient force and free convection along a vertical wavy surface in micropolar fluid. Hossain et al. [6] have studied the problem of natural convection of fluid with temperature dependent viscosity along a heated vertical wavy surface. Natural and mixed convection heat and mass transfer along a vertical wavy surface have been investigated by Jang [7,8]. Recently, Molla et al. [9] have studied natural convection flow along a vertical wavy surface with uniform surface temperature in presence of heat generation/absorption. Tashtoush and Al-Odat [10] investigated magnetic field effect on heat and fluid flow over a wavy surface with a variable heat flux. Hossain [11] investigated the natural convection flow past a permeable wedge for the fluid having temperature dependent viscosity and thermal conductivity. Very recently, Parveen and Alim [12] investigated Joule heating effect on Magnetohydrodynamic natural convection flow along a vertical wavy surface with viscosity dependent on temperature. The present study is to incorporate the effects of the viscous dissipation on MHD natural convection flow along a uniformly heated vertical wavy surface with heat generation.

Numerical results have been obtained in terms of local skin friction coefficient and the rate of heat transfer in terms of local Nusselt number, and the velocities as well as the temperature profiles for a selection of relevant physical parameters are shown graphically.

2. Formulation of the Problem

Steady two dimensional laminar free convection boundary layer flow of a viscous incompressible and electrically conducting fluid along a vertical wavy surface in presence of uniform transverse magnetic field is considered. It is assumed that the wavy surface is electrically insulated and is maintained at a uniform temperature T_{w}. The fluid is stationary above the wavy plate and is kept at a temperature T_{¥}. The surface temperature T_{w} is greater than the ambient temperature T_{¥} that is, T_{w} > T_{¥}. The flow configuration of the wavy surface and the two-dimensional cartesian coordinate system are shown in Figure 1.

The boundary layer analysis outlined below allows being arbitrary, but our detailed numerical work assumed that the surface exhibits sinusoidal deformations. The wavy surface may be defined by

(1)

where a is the amplitude and L is the wave length associated with the wavy surface.

The governing equations of such flow of magnetic field in presence of heat generation/absorption with viscosity variation along a vertical wavy surface under the usual Boussinesq approximations can be written in a dimensional form as:

Continuity Equation

(2)

Figure 1. Physical model and coordinate system.

(3)

(4)

(5)

where are the dimensional coordinates along and normal to the tangent of the surface and are the velocity components parallel to, g is the acceleration due to earth gravity, P is the dimensional pressure of the fluid, T is the temperature of the fluid in the boundary layer, C_{P} is the specific heat at constant pressure, μ is the dynamic viscosity of the fluid in the boundary layer region depending on the fluid temperature, ρ is the density, is the kinematic viscosity, where

,_{,} k is the thermal conductivity of the fluid, β is the volumetric coefficient of thermal expansion, b_{0} is the strength of magnetic field,_{ }is the electrical conductivity of the fluid and is the Laplacian operatorwhere

The boundary conditions for the present problem are

(6)

Using Prandtl’s transposition theorem to transform the irregular wavy surface into a flat surface as extended by Yao [1] and boundary layer approximation, the following dimensionless variables are introduced for non-dimensionalizing the governing equations

(7)

where θ is the dimensionless temperature function and are the dimensionless velocity components parallel to and Gr is the Grashof number. Now introducing the dimensionless dependent and independent variables into Equations (2)-(5), and the following dimensionless form of the governing equations is obtained after ignoring terms of smaller orders of magnitude in the Grashof number Gr.

(8)

(9)

(10)

(11)

It is worth noting that the σ_{x} and σ_{xx} indicate the first and second derivatives of σ with respect to x, therefore,

and.

In the above equations Pr, Vd, M and Q are respectively known as the Prandtl number, viscous dissipation parameter, magnetic parameter and heat generation parameter which are defined as:

(12)

For the present problem the pressure gradient is zero. Thus, the elimination of from Equations (9) and (10) leads to

(13)

The corresponding boundary conditions for the present problem then turn into

(14)

Now we introduce the following transformations to reduce the governing equations to a convenient form:

(15)

where is the dimensionless stream function, η is the dimensionless similarity variable and ψ is the stream function that satisfies the continuity Equation (8) and is related to the velocity components in the usual way as

(16)

Introducing the transformations given in Equation (15) and using (16) into Equations (13) and (11) are transformed into the new co-ordinate system. Thus the resulting equations are

(17)

(18)

The boundary conditions (14) now take the following form:

(19)

Here prime denote the differentiation with respect to η.

However, once we know the values of the functions f and q and their derivatives, it is important to calculate the values of the rate of heat transfer in terms of local Nusselt number Nu_{x} and the shearing stress t_{w} in terms of the local skin friction coefficient C_{fx} from the following relations:

(20)

where

(21)

Here is the unit normal to the surface.

Using the transformation (15) and (21) into Equation (20) the local skin friction coefficient C_{fx} and the rate of heat transfer in terms of the local Nusselt number Nu_{x} take the following forms:

(22)

(23)

For the computational purpose the period of oscillations in the waviness of this surface has been considered to be π.

3. Results and Discussion

We have investigated the effects of viscous dissipation on natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface. Although there are five parameters of interest in the present problem, the effects of Prandtl number Pr, viscous dissipation Vd, magnetic parameter M, the heat generation parameter Q and the amplitude of the wavy surface a on the surface shear stress in terms of local skin friction coefficient, the rate of heat transfer in terms of the local Nusselt number, the velocity and temperature profiles, the streamlines and the isotherms. Numerical values of local shearing stress and the rate of heat transfer are calculated from Equations (22) and (23) in terms of the skin-friction coefficients C_{fx} and Nusselt number Nu_{x} respectively for a wide range of the axial distance variable x starting from the leading edge for different values of the parameters Pr, Vd, M, Q and a. Solutions are obtained in terms of velocity profiles, temperature profiles against h and the skin friction coefficients C_{fx}, the rate of heat transfer in terms of the Nusselt number Nu_{x} at any position of x presented graphically for selected values of magnetic parameter M = 0.0, 0.5, 1.5, 2.5, 5.0, viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0 and heat generation parameter Q = 0.0, 0.5, 1.0, 1.5, 2.0.

The effects for different values of magnetic parameter M on the velocity and temperature profiles with, and have been presented graphically in Figures 2(a) and (b). It is seen from the Figures 2(a) that for the values of magnetic parameter M = 0.0, 0.5, 1.5, 2.5, 5.0 the velocity decreasing upto the position of from the wall. At the position of

velocity becomes constant that is velocity profiles meet at a point and then cross the side and increasing with magnetic parameter M. This is because of the velocity profiles having lower peak values for higher values of magnetic parameter M tend to decreases comparatively slower along h-direction than velocity profiles with higher peak values for lower values of magnetic parameter M. The maximum values of velocities are recorded as 0.49091, 0.47285, 0.43933, 0.40912 and 0.34643 for magnetic parameter M = 0.00, 0.50, 1.50, 2.50, 5.00 respectively which occur at the same position. Here, it is observed that at, the maximum velocity decreases by 29.43% as the magnetic parameter M change from 0.0 to 5.0. The values of temperature are recorded as 0.70911, 0.71666, 0.73133, 0.74536 and 0.77742 for magnetic parameter M = 0.00, 0.50, 1.50, 2.50, 5.00 at the same position of 1.23788 and the temperature increases by 9.63%. In Figures 3(a) and (b) the effects for different values of the viscous dissipation parameter Vd on the velocity and temperature profiles with and have been shown graphically. It has been seen from Figure 3(a) that as the viscous dissipation parameter Vd increases, the velocities rising up to the position of for viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0 and from that position of h velocities fall down slowly and finally approaches to zero. It is also observed from Figure 3(b) that as the viscous dissipation parameter Vd increases, the temperature profiles increases. The maximum values of velocities are recorded as 0.48686, 0.48888, 0.49091, 0.49499 and 0.49909 for viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0 respectively which occur at the same position and the maximum velocity increases by 2.51%. Temperatures are recorded as 0.69855, 0.70380, 0.70911, 0.71991 and 0.73095 for viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0 respec-

(a)(b)

Figure 2. (a) Velocity and (b) Temperature profiles against h for different values of M with Pr = 0.72, a = 0.2, Vd = 10.0 and Q = 0.4.

(a)(b)

Figure 3. (a) Velocity and (b) Temperature profiles against h for different values of Vd with Pr = 0.72, a = 0.2, M = 0.1 and Q = 0.4.

tively at the same position of and the temperature profiles increases by 4.64%. Both the velocity and temperature profiles accumulate nearly in the following points where and respectively for viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0. That is, velocity boundary layer thickness and thermal boundary layer thickness are unchanged. The effects for different values of the heat generation parameter Q = 0.0, 0.5, 1.0, 1.5, 2.0 on the velocity and temperature profiles with, , and have been presented graphically in Figures 4(a) and (b) respectively. For the higher values of the heat generation parameter Q both the velocity and the temperature increases.

In Figures 5(a) and (b) effects of magnetic parameter M = 0.0, 0.5, 1.5, 2.5, 5.0 on skin friction and the rate of heat transfer with and v have been presented. From Figure 5(a) it is found that skin friction decreases significantly for greater magnetic field strength. This is physically realizable as the magnetic field retards the velocity field and consequently reduces the frictional force at the wall. However rate of heat transfer opposite pattern due to the higher values of magnetic parameter M which are presented in Figure 5(b). The different values of viscous dissipation parameter Vd with and. of the skin friction coefficients and the rate of heat transfer are shown graphically in Figures 6(a) and (b) respectively. In this case the values of local skin friction coefficient C_{fx} are recorded to be 1.12663, 1.48796, 2.03818, 4.01586 and 7.94542 for Vd = 0.0, 5.0, 10.0, 20.0, 30.0 which occur at same point. From the Figure 6(a), it is observed that at, the skin friction coefficient increases by 605.23% due to the higher value of viscous dissipation parameter Vd. However, the values of rate of heat transfer are found to be −0.07240, −3.63124, −11.82060, −66.89554, and −284.88541 for Vd = 0.0, 5.0, 10.0, 20.0, 30.0 which occur at same point. The rate of heat transfer coefficient decreases by 393.39% due to the increased value of viscous dissipation parameter Vd = 0.0, 5.0, 10.0, 20.0, 30.0. It is seen from the Figure 6(b) that for higher values of viscous dissipation parameter the rate of heat transfer decreases that is heat transfer slows down for higher viscous dissipation parameter Vd. In Figures 7(a) and (b) the skin friction coefficient C_{fx} and local rate of heat transfer Nu_{x} for different values of heat generation parameter Q = 0.0, 0.5, 1.0, 1.5, 2.0 with and Vd = 10.0 have been displayed. It is observed from the Figure 7(a) that for higher values of heat generation parameter Q, skin friction increasing up to the axial position of x = 0.6 and then skin friction becomes constant for all values of Q that is, skin friction coefficient meet together at the position of x = 0.6 and cross the sides that means after the axial position of x = 0.6 skin friction is decreasing for higher values of heat generation parameter but frictional force at the wall always rising towards downstream. It is seen from the Figure 7(b) that for higher values of heat generation parameter the rate of heat transfer decreases down to the axial position of x = 1.3 and then rate of heat transfer becomes constant for all values of Q that is, rate of heat transfer meet together at the position of x = 1.3 and cross the sides that means after the axial position of x = 1.3 rate of heat transfer is increasing up for higher values of heat generation.

In Figure 8(a) and (b) are shown that the effects of streamlines and isotherms for different values of heat generation parameter Q = 0.0 and 2.0 with amplitude of waviness of the surface, Prandtl number Pr = 0.72, magnetic parameter M = 0.1 and viscous dissipation parameter Vd = 10.0 respectively. The value of stream function ψ is 0.00 near the wall and then ψ increases gradually within the boundary layer and away from the wall. On the other hand the isolines of temperature (iso-

(a)(b)

Figure 4. (a) Velocity and (b) Temperature profiles against h for different values of Q with Pr = 0.72, a = 0.2, M = 0.1 and Vd = 10.0.

(a)(b)

Figure 5. (a) Skin friction coefficient and (b) Rate of heat transfer against x for different values of M with Pr = 0.72, a = 0.2, Vd =10.0 and Q = 0.4.

(a)(b)

Figure 6. (a) Skin friction coefficient and (b) Rate of heat transfer against x for different values of Vd with Pr = 0.72, a = 0.2, M = 0.1 and Q = 0.4.

therms) distribution shows that temperature decreases significantly as the values of the heat generation parameter Q increases which have been presented in Figure 8(b). The value of isotherm is 1.0 at the wall and isotherms decreases slowly along the y-direction and finally approach to zero.

(a)(b)

Figure 7. (a) Skin friction coefficient and (b) Rate of heat transfer against x for different values of Q with Pr = 0.72, a = 0.2, M = 0.1 and Vd =10.0.

(a)(b)

Figure 8. (a) Streamlines and (b) Isotherms for Q = 0.0 (Red solid lines), Q = 2.0 (Black dashed lines), with Pr = 0.72, a = 0.2, M = 0.1 and Vd = 10.0.

4. Conclusions

The effects of the Prandtl number Pr, the magnetic parameter M, the viscous dissipation parameter Vd, the heat generation parameter Q and the amplitude of wavy surface a on MHD natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface have been studied. From the present investigation the following conclusions may be drawn:

The temperature and the rate of heat transfer coefficient increase for increasing values of magnetic parameter. The velocity decreases and at the position of becoming constant that is velocity profile meets at the point and then crosses the side and increases with magnetic parameter. The local skin friction coefficient decreases due to the increased value of magnetic parameter.

The velocity and the temperature rise up and the local skin friction coefficient increase due to the higher values of viscous dissipation parameter Vd which cause reduction of the rate of heat transfer.

The velocity, temperature and the skin friction coefficient enhance for higher values of internal heat generation parameter Q but for the same reason the rate of heat transfer reduces.

REFERENCES

- L. S. Yao, “Natural Convection along a Vertical Wavy Surface,” ASME Journal of Heat Transfer, Vol. 105, No. 3, 1983, pp. 465-468. doi:10.1115/1.3245608
- S. G. Moulic and L. S. Yao, “Mixed Convection along Wavy Surface,” ASME Journal of Heat Transfer, Vol. 111, No. 4, 1989, pp. 974-979. doi:10.1115/1.3250813
- M. A. Alim, S. Alam and M. Miraj, “Effects of Temperature Dependent Thermal Conductivity on Natural Convection Flow along a Vertical Wavy Surface with Heat Generation,” International Journal of Engineering & Technology, Vol. 11, No. 6, 2011, pp. 60-69.
- M. A. Hossain and D. A. S. Rees, “Combined Heat and Mass Transfer in Natural Convection Flow from a Vertical Wavy Surface,” Acta Mechanica, Vol. 136, No. 3-4, 1999, pp. 133-141. doi:10.1007/BF01179253
- C. C. Wang and C. K. Chen, “Transient Force and Free Convection along a Vertical Wavy Surface in Micropolar Fluid,” International Journal of Heat and Mass Transfer, Vol. 44, No. 17, 2001, pp. 3241-3251. doi:10.1016/S0017-9310(00)00329-X
- K. C. A. Alam, M. A. Hossain and D. A. S. Rees, “Magnetohydrodynamic Free Convection along a Vertical Wavy Surface,” Journal of Applied Mechanical Engineering, Vol. 1, 1997, pp. 555-566.
- J. H. Jang, W. M. Yan and H. C. Liu, “Natural Convection Heat and Mass Transfer along a Vertical Wavy Surface,” International Journal of Heat Mass Transfer, Vol. 46, No. 6, 2003, pp. 1075-1083. doi:10.1016/S0017-9310(02)00361-7
- J. H. Jang and W. M. Yan, “Mixed Convection Heat and Mass Transfer along a Vertical Wavy Surface,” International Journal of Heat Mass Transfer, Vol. 47, No. 3, 2004, pp. 419-428. doi:10.1016/j.ijheatmasstransfer.2003.07.020
- M. M. Molla, M. A. Hossain and L. S. Yao, “Natural Convection Flow along a Vertical Wavy Surface with Uniform Surface Temperature in Presence of Heat Generation/Absorption,” International Journal of Thermal Sciences, Vol. 43, No. 2, 2004, pp. 157-163. doi:10.1016/j.ijthermalsci.2003.04.001
- B. Tashtoush and M. Al-Odat, “Magnetic Field Effect on Heat and Fluid flow over a Wavy Surface with a Variable Heat Flux,” Journal of Magnetism and Magnetic Materials, Vol. 268, No. 3, 2004, pp. 357-363. doi:10.1016/S0304-8853(03)00547-X
- M. A. Hossain, S. Kabir and D. S. A. Rees, “Natural Convection of Fluid with Temperature Dependent Viscosity from Heated Vertical Wavy Surface,” Zeitschrift für Angewandte Mathematik und Physik, Vol. 53, No. 1, 2002, pp. 48-57. doi:10.1007/s00033-002-8141-z
- N. Parveen and M. A. Alim, “Joule Heating Effect on Magnetohydrodynamic Natural Convection flow along a Vertical Wavy Surface with Viscosity Dependent on Temperature,” International Journal of Energy & Technology, Vol. 3, 2011, pp. 1-10.

Nomenclature

C_{fx} Local skin friction coefficient

C_{p} Specific heat at constant pressure [J·kg^{−1}·K^{−1}]

f Dimensionless stream function

g Acceleration due to gravity [ms^{−2}]^{}

Gr Grashof number

k Thermal conductivity [Wm^{−1}·K^{−1}]

k_{¥} Thermal conductivity of the ambient fluid [Wm^{−1}·K^{−1}]

L Characteristic length associated with the wavy surface [m]

n Unit normal to the surface

Nu_{x} Local Nusselt number

P Pressure of the fluid [Nm^{−2}]

Pr Prandtl number

Q Heat generation parameter

Q_{0} Heat generation constant

q_{w} Heat flux at the surface [Wm^{−2}]

T Temperature of the fluid in the boundary layer [K]

T_{w} Temperature at the surface [K]

T_{¥}_{ }Temperature of the ambient fluid [K]

u, v Dimensionless velocity components along the (x, y) axes [ms^{−1}]

x, y Axis in the direction along and normal to the tangent of the surface

Greek Symbols

α Amplitude of the surface waves

β Volumetric coefficient of thermal expansion [K^{−1}]

η Dimensionless similarity variable

θ Dimensionless temperature function

ψ Stream function [m^{2}·s^{−1}]

μ Viscosity of the fluid [kg·m^{−1}·s^{−1}]

μ_{¥} Viscosity of the ambient fluid

ν Kinematic viscosity [m^{2}·s^{−1}]

ρ Density of the fluid [kg·m^{−3}]

σ_{0 }Electrical conductivity

τ_{w} Shearing stress