Thermal Radiation , Joule Heating , and Viscous Dissipation Effects on MHD Forced Convection Flow with Uniform Surface Temperature

In this paper, we studied the effects of thermal radiation, Joule heating and viscous dissipation on forced convection flow in a magnetohydrodynamics (namely MHD) pump in rectangular channel with uniform surface temperature. Numerical results were obtained by solving the nonlinear governing momentum and energy equations with steady state fully developed assumptions by finite difference method. The Lorentz force in momentum and Joule heating, and viscous dissipation in energy equation with the Rossel and approximation are assumed to increase the knowledge of the details of the temperature and flow field in order to design a MHD pump. The purpose of this study is the parametric study of a Newtonian fluid in a MHD pump. The values of maximum velocity, fully developed Nusselt number for different values of magnetic density flux, Brinkman number, viscous heating and radiation number are obtained. However, the maximum temperature stays almost constant with magnetic field, as current increases, the velocity and the temperature increase too. Besides, the increase of thermal radiation number causes the increase in effective thermal conductivity and decrease in thermal boundary layer and the Nusselt number at wall.


Introduction
Laminar heat transfer to incompressible viscous Newtonian and non-Newtonian fluids in tubes is a problem of considerable industrial significance, due to its application especially in heat exchanger designs, packed-bed combustors, catalytic converters, materials-processing applications, such as extrusion, metal forming, continuous casting, as well as wire and glass fiber drawing.This problem has been therefore considered in many papers.A detailed review of literature being performs here.The work of Brinkman [1] appears to be the first theoretical work dealing with heat generated by viscous dissipation.The temperature distribution in the entrance region of a circular pipe at the wall of which was maintained at either the constant temperature of the entering fluid or constant heat-flux was examined.The highest temperatures were, not surprisingly, discovered to be localized in the wall region.Lin [2] studied laminar heat-transfer to a non-Newtonian Couette flow with pressure gradient using the power-law model.The effects of pressure gradient and viscous dissipation on the heat transfer were discussed.Using a functional analysis method, Lahjomri [3] analytically studied thermally-developing laminar Hartman flow through a parallel-plate channel, with a prescribed transversal uniform magnetic field, including both viscous dissipation, Joule heating and axial heat-conduction with uniform heat-flux.In a recent study, Nield et al. [4] investigated the thermal development of forced convection in a parallel plate channel filled by a saturated porous medium, with walls held at a uniform temperature, and with the effects of axial conduction and viscous dissipation included.Davaa et al. [5] numerically studied fully-developed laminar heattransfer to non-Newtonian fluids flowing between parallel plates with the axial movement of one of the plates with an emphasis on the viscous-dissipation effect.Increasing the Brinkman number increased the heat-transfer rates at the heated wall when the movement direction of the upper plate was the same as the direction of the main flow, while the opposite is true for the movement of the upper plate in the opposite direction.Hashemabadi et al. [6] obtained an analytical solution to predict the fully-developed, steady and laminar heat-transfer of viscoelastic fluids between parallel plates.Their results emphasized the significant effect of viscous heating on the Nusselt number.Furthermore recently Hamid et al. [7] studied the effects of radiation on Marangoni convection over a flat surface, and Chamkha et al. [8] studied that effect on a non-isothermal wedge not in a channel.
The aim of the present study is to investigate analytically the effect of viscous dissipation on steady-state laminar heat transfer in a Poiseuille flow between plane parallel plates in a rectangular channel with a simultaneous pressure gradient, Lorentz force and joule heating effect.The effect of Brinkman number on the temperature profile and the Nusselt number is obtained for the constant temperature boundary-condition from radiate surfaces.

Governing Equations and Numerical Method
Consider steady, hydrodynamically and thermally fully-developed, laminar flow of an incompressible fluid between two parallel plates (Figure 1).The thermal conductivity and the thermal diffusivity of the fluid are considered to be independent of temperature.Fluid viscosity is modeled by the power law.Furthermore magnetic and electric properties are constant.The both plate is assumed to be in a stationary.
The axial heat-conductions in the fluid and in the wall are neglected.The momentum equation in the x-direction is described as final stage In which Lorentz force can be written as: ( ) where E is electric field intensity in z-direction, B is magnetic density in y-direction and σ is electrical conductivity of the fluids.The conservation of energy, including the effect of the viscous dissipation, can be written as follows: where the second term on the right-hand side is the viscous-dissipation term and third term is joule heating.Also the Energy equation by considering the Rossel and approximation for the radiative flux is used which assumes that the radiative heat flux vector in a gray optically thick medium corresponds to the blackbody emissive power ( ( )( ) ) and suggest the effective thermal conductivity ( ( ) ) [9].Under the following boundary conditions: When the temperature profile is fully developed one can write 0 x θ ∂ = ∂ (Bejan [10]).The thermal boundary conditions of Equation ( 5) are )

Results and Discussion
To solve the governing nonlinear differential equations, an in-house finite-volume code is developed and utilized.The code is based on structured grid.It generates a uniform grid consisting of a number of nodes in the yand z-directions.This method is utilized to assure high accuracy.A minimum of 2000 nodes are assumed.The residuals for velocity and temperature are set at 1e − 6 .The program uses the input data and solve the momentum equations to obtain velocity field along the MHD pump.Then, temperature distribution is obtained by solving the energy equation.For verification purposes, Figure 3 shows the current as it is compared with a computational model known case (flow inside a rectangular duct) by Wang et al. [11].As shown, the volumetric flow rates for both studies show a similar trend.
The effect of magnetic flux density on maximum velocity and volumetric flow rate of the MHD pump is shown by Figure 4.The effect of electric current on velocity profile and volumetric flow rate of the MHD pump is shown by Figure 5.As shown, the maximum velocity and then volumetric flow rate increase with increasing magnetic flux density and electric current.The effect of Brinkman number and MHD heat generation on temperature profiles are presented of Figure 6 and Figure 7 respectively.The Brinkman number increase the temperature profile near the wall while the joule heating increase the temperature profile near the center of the channel.As demonstrated the maximum temperature increase with increasing electric current but is not varied considerably with magnetic flux density variations.Also Figure 8 and Figure 9 present the linear increasing effect of Brinkmann number and magnetic flux on Nusselt number, respectively.Reverse effect for radiative heating on Nusselt number presented on Figure 10.

Conclusion
In this study, a finite-difference code is used for solving momentum and energy equations and the velocity and the temperature field has been derived for the laminar, steady, convective heat-transfer problem in a Poiseuille flow between plane parallel plates with a simultaneous pressure-gradient and the Lorentz force.Interest has been focused on the influence of the viscous dissipation and internal heating ratio.The effect of the Brinkman number on the Nusselt number has been discussed in terms of the energy balance.The maximum velocity and volumetric flow rate increase with increasing magnetic flux density and electric current.The maximum temperature increases with increasing electric current but is not varied considerably with magnetic flux density variations.As illustrated the Brinkman number increases the temperature profile near the wall while the joule heating increases the temperature profile near the center of the channel.The results show that for the thermally fully-developed, a strong influence of viscous dissipation on the heat transfer for higher values of Br (Br > 1) has been observed, while this influence is found to be negligible for lower values of Br.Finally by increase of Brinkmann number and magnetic flux, and decrease of radiative heating has linear increasing effect on Nusselt number.

Figure 3 .
Figure 3.Comparison of the flow velocity as a function of channel width by [11] (channel depth = 2 mm, I = 0.5 A , B = 18 mT).

Figure 5 .
Figure 5. Maximum velocity as a function of electric current, where B is kept constant at 0.4.