Effect of Brinkman Number and Magnetic Field on Laminar Convection in a Vertical Plate Channel
Total Page:16
File Type:pdf, Size:1020Kb
Science World Journal Vol 12(No 4) 2017 www.scienceworldjournal.org ISSN 1597-6343 EFFECT OF BRINKMAN NUMBER AND MAGNETIC FIELD ON LAMINAR CONVECTION IN A VERTICAL PLATE CHANNEL 1* 2 3 Joseph Kpop Moses ; Peter Ayuba and Abubakar S. Magaji Full Length Research Article 1, 2, 3 Department of Mathematical Sciences, Kaduna State University – Nigeria. *Corresponding Author’s e-mail address: [email protected] ABSTRACT Gav et al (1992) showed experimentally for mixed convection that The effect of Brinkman number and magnetic field on laminar the increasing buoyancy parameter will make the reversal flow convection in a vertical plate channel with uniform and wider and deeper taking into consideration the effects of viscous asymmetric temperatures has been studied. The dimensionless dissipation. Barletta (1998) demonstrated that in the case of form of momentum and energy balanced equations has been upward flow reversal decreases with viscous dissipation. solved using one term perturbation series solution. The solution of However, in the case of downward flow, the flow reversal the dimensionless velocity u and temperature θ were obtained. increases with viscous dissipation. He also calculated criteria for The numerical values of the skin friction and the Nusselt number the onset of flow reversal for both directions of flow, but without were obtained and tabulated. It was found that the effect of taking viscous dissipation into account. magnetic parameter enhances the effect of flow velocity by Barletta (1999) investigated the fully developed and laminar suppressing the turbulence while Brinkman number has convection in a parallel – plate channel by taking into account accelerating effect on the temperature profile. both viscous dissipation and buoyancy. Uniform and asymmetric temperatures are prescribed at the channel walls. The velocity Keywords: MHD, Laminar, Convection, vertical channel, field is considered as parallel. A perturbation method is employed Brinkman number to solve the momentum balance equation and energy balance equation. A comparison between the velocity and temperature NOMENCLATURES profiles in the case of laminar forced convection with viscous ν - Kinematic viscosity dissipation is performed in order to point out the effect of 훼 - Thermal conductivity buoyancy. The case of convective boundary conditions is also 휇- Dynamic viscosity discussed. 푇푏 - Bulk temperature Daniel et.al. (2010) investigated the fully developed and laminar 푈푚 − Mean velocity convection in a vertical parallel plate channel by taking into ∆푇 - Temperature difference account both viscous dissipation and buoyancy. Uniform and 푀 표푟 퐻푎 − magnetic parameter asymmetric temperatures are prescribed at the channel walls. 퐺푟 -Grashof number The velocity field is considered as parallel. A perturbation method 푅푒 -Reynolds number is employed to solve the momentum balance equation and the - Acceleration due to gravity, energy balance equation. A comparison with the velocity and 퐵푟 - Brinkman number temperature profiles in the case of laminar forced convection with various dissipation is performed in order to point out the effect the INTRODUCTION effect of buoyancy. The case of convective boundary condition is The study of magnetohydrodynamic (MHD) flows with viscous also discussed. dissipation and buoyancy effects has many importance in Joseph et.al. (2013) studied the MHD forced convection in a industrial, technological and geothermal applications such as horizontal double – passage with uniform wall heat flux by taking temperature plasmas, cooling of nuclear reactors, liquid metals into account the effect of magnetic parameter. They assumed the etc. A comprehensive review of the study of MHD flows in flow of the fluid to be laminar, two – dimensional, steady and fully relations to the applications to the above areas has been made by developed. The fluid is incompressible and the physical properties several authors. are constants. Dileep and Rashin (2012) investigated the fully Mixed convection in a vertical parallel – plate channel with developed magntohydrodynamic mixed convection in a vertical uniform wall temperatures has been studied by Tao (1960). channel, partially filled with clear fluid and partially, filled with a Javerci (1975) dealt with the effects of viscous dissipation and fluid – saturated porous medium is by taking into account viscous Joule heating on the fully developed MHD flow with heat transfer and ohmic dissipations. Alim, et al (2007) investigated the in a channel. The exact solution of the energy equation was pressure work and viscous dissipation effects on MHD natural derived for constant heat flux with small magnetic Reynolds convection along a sphere. The laminar natural convection flow number. Sparrow et al (1984) observed experimentally flow from a sphere immersed in a viscous incompressible fluid in the reversal in pure convection in a one – sided heated vertical presence of magnetic field has been considered in this channel. The mixed convection in a vertical parallel – plate investigation. channel with asymmetric heating, where one plate is heated and Barik (2013) analyzed the free convection heat and mass transfer the other is adiabatic has been studied by Habchi and Acharya MHD flow in a vertical channel in the presence of chemical (1986). reaction. Ranuka et al (2009) studied the MHD effects of unsteady heat convective mass transfer flow past an infinite vertical porous plate with variable suction, where the plate’s Effect of Brinkman Number and Magnetic Field an Laminar Convection in a 58 Vertical Plate Channel Science World Journal Vol 12(No 4) 2017 www.scienceworldjournal.org ISSN 1597-6343 2 2 2 2 temperature oscillates with the same frequency as that of variable 푑 푇 ν 푑푈 푒퐵0 푈 훼 2 + ( ) + = 0 (2) suction velocity with soret effects. The governing equations are 푑푌 푐푝 푑푌 0푐푝 solved numerically by using implicit finite difference method. Where the acceleration due to gravity is g, ν is the kinematic 휇 Saleh and Hashim (2009) analysed flow reversal phenomena of viscosity which is defined by ν= , 푃 = 푝 + 휌푏푋 is the the fully – developed laminar combined free and forced MHD 푏 difference between the pressure p and the hydrostatic convection in a vertical plate – channel. The effect of viscous pressure P. 훼 is thermal conductivity and 휇 = 휌 휈 is the dissipation is taken into account. Tamia and Samad (2010) 푏 dynamic viscosity studied and analyzed the radiation and viscous dissipation effects Equations (2) and the boundary condition on푇, shows that 푃 on a steady two – dimensional magneto hydrodynamics free 휕푃 depends on 푋; 푇 depends only on 푌; 푇 is a constant; is convection flow along a stretching sheet with heat generation. 푏 휕푋 The non – linear partial differential equations governing the flow a constant. field under consideration have been transformed by a similarity The bulk temperature 푇푏 and the mean velocity 푈푚 (Barletta, transformation into a system of non – linear ordinary differential 1999) are given by. equations and then solved numerically by applying Nachtsheim – 1 퐿 1 퐿 Swigert shooting iteration technique together with sixth order and (3) 푇푏 = ∫0 푈푇푑푌 푈푚 = ∫0 푈푑푌 Runge – Kutta integration scheme. Resulting in non – 퐿푈푚 퐿 dimensional velocity, temperature and concentration profiles are The boundary conditions are: then presented graphically for different values of the parameters 푈(0) = 0 , 푇(0) = 푇0, of physical engineering interest. 푈(퐿) = 0 , 푇(퐿) = 푇0 (4) Idowu et.al (2013) investigated viscous dissipation and buoyancy Due to the symmetry of the problem, 푈 푎푛푑 푇 are symmetric effects on laminar convection with transpiration. Uniform and functions of 푌, so that they must fulfill the condition asymmetric temperatures are prescribed at the channel walls. 푈′(0) = 0, 푇′(0) = 0 (5) The velocity field is considered as parallel. A perturbation method Introducing the dimensionless quantities below 푈 푇−푇 푌 −퐿2 푑푃 is employed to solve the momentum balance equation and the 푢 = , 휃 = 푏 , 푦 = , 휆 = , 퐺푟 = energy balance equation. A comparison with the velocity and 푈푚 ∆푇 퐿 휇푈푚 푑푋 64퐿3훽∆푇 , temperature profiles in the case of laminar forced convection with 휈2 various dissipation is performed in order to point out the effect of 4퐿푈 퐺푟 휇푈2 푅푒 = 푚 , 퐸 = , 퐵푟 = 푚 buoyancy. The case of convective boundary condition is also 휈 푅푒 훼∆푇 discussed. Where the temperature difference ∆푇 is given by 2 Das et.al. (2015) studied the fully developed mixed convection 휇푈푚 ∆푇 = flow in a vertical channel filled with nanofluid in the presences of 푘 uniform transverse magnetic field. Anjali Devi and Ganga (2010) While the Grashof number 퐺푟 is always positive, the Reynolds studied the Viscous and Joule dissipation effects on MHD non – number 푅푒 and the parameter 퐸 can be either positive or linear flow and heat transfer past a stretching porous surface negative. In particular, in the case of upward flow (푈푚 > embedded in a porous medium under a transverse magnetic field. 0), both 푅푒 and 퐸 are positive, while, for downward flow(푈푚 < Patrulescu, et al (2010) investigated the steady mixed convection 0), these dimensionless parameters are negative. flow in a vertical channel for laminar and fully developed flow regime. In the modeling of heat transfer the viscous dissipation Therefore, as a consequence of the dimensional quantities term was also considered. Temperature on the right wall is together with the momentum balance equation (1), the energy assumed constant while a mixed boundary condition (Robin balance equation (2) and equations (3) – (5) yields; boundary condition) is considered on the left wall. 푑2푢 퐸휃 The mixed convection of Newtonian fluid between vertical parallel − 푀푢 = − + 휆 plates channel with MHD effect and variation in Brinkman number 푑푦2 16 2 2 has been studied by Alizadeh Rasul et.al. (2014). The effect of 푑 휃 2 푑푢 2 + 퐵푟푀푢 = −퐵푟 ( ) (6) MHD and Brinkman number on laminar mixed convection of 푑푦 푑푦 훽2퐿2 Newtonian fluid between vertical parallel plates channel was Where 푀 = 0 0 which is the magnetic parameter and 0휈 investigated by Salehi et.al (2013).