BUOYANCY EFFECTS OF THERMAL AND MASS DIFFUSION ON MHD NATURAL CONVECTION PAST FINITE VERTICAL, POROUS FLAT PLATE

The flow of an electrically conducting incompressible fluid due to buoyancy effects of thermal and mass diffusion past a finite vertical porous plate with constant suction was investigated in the presence of uniform transverse magnetic field. The problem has been solved for velocity, temperature and concentration profiles. The equations governing the flow are solved numerically using finite difference method for various values of Grashof parameter ranging from 0 to – 1. The results obtained are then presented using tables and graphs. It was noted that a decrease in Grashof parameter leads to an increase in primary, secondary, temperature and concentration profile. AMS Subject Classification: 76E06

The phenomenon of natural convection arises in fluids when temperature change causes density variation leading to buoyancy forces acting on the fluid particles.These flows which are driven by temperature differences abound in nature and have been studied extensively because of their application such as hot rolling and crystal growing glass fibre.In nature there exist flows that are caused not only by temperature differences but also by concentration differences.In industries many transport processes exists in which heat and mass transfer takes place simultaneously as a result of combined buoyancy effect of thermal diffusion and mass diffusion through chemical species.Free convection flows involving coupled heat and mass transfer occur frequently in nature.For this flows, the driving forces arises due to temperature and concentration variation of the fluid.For example in atmospheric flows, thermal convection resulting from heating of the earth by sunlight is affected by difference in water vapour concentration.A comprehensive literature on various aspects of natural convectional flow and its application coupled heat and mass transfer (Double diffusion) driven by buoyancy due to temperature and concentration variations in a saturated porous medium , has several important application in geothermal and geophysical engineering such as the migration of moisture through the air contained in fibrous insulation, the extraction of geothermal energy, underground disposal of nuclear waste, and the spreading of chemical contaminants through water saturated soil.The problem of combined heat and mass transfer over electrically conducting fluid in MHD free convection adjacent to a vertical surface are analysed by taking into account the effects of ohmic heating and viscous dissipation.Buoyancy forces resulting from thermal and mass diffusion is of considerable interest in nature and many industrial applications such as Geo-physics, solidification of binary alloys, drying process and chemical engineering.The analytical 1-D model compared favourable with the experiment results when modified to account for the hall parameter, the parameter was evaluated and correlation was developed.The hall parameter had a strong dependence disturbed the power output from the electrodes of high temperatures.

Literature Review
[3] presented free convection with mass transfer flow for a micro-polar fluid bounded by vertical infinite surface under the action of transverse magnetic field, the differential equations governing the problem were solved numerically by means of order Runge-kutta method.It was observed that the velocity decreases as prandtl number increases.Finally the concentration distribution is affected only by Schmidt number and decreases as this parameter increases.[4] analysed the effect of thermal radiation and viscous dissipation on MHD free convection flow over a semi-infinite vertical porous plate.The network simulation method was used to solve the boundary layer equation based on the finite difference formulation.It was found that an increase in viscous dissipation leads to an increase of both velocity and temperature profiles and a decrease in the velocity profiles.Finally an increase in the suction parameter leads to an increase in the local skin friction and Nusselt number.Also, when the Schmidt number was increased, the concentration level was decreased resulting in decreased fluid velocity.In addition, it was found that the skin friction coefficient increased due to increase in concentration buoyancy effect which it decreased due to increase in either heat absorption coefficient or Schmidt number.
[1] considered unsteady 2-D laminar, boundary layer flow over viscous, incompressible, electrically conducting and heat absorbing fluid along a semiinfinite vertical permeable moving plate in the presence of a uniform transverse magnetic field.Thermal and concentration buoyancy were discussed.The dimensionless governing equations for this investigation are solved analytically using two term harmonic and non-harmonic function, it was found that when the solution Grashof number increased, the concentration buoyancy effects were enhanced and thus, fluid velocity increased .however, the presence of heat absorption effect caused a reduction in the fluid temperature which resulted in a decrease in fluid velocity.Also, when the Schmidt number was increased, the concentration level was decreased resulting in decreased fluid velocity.In addition, it was found that the skin friction coefficient increased due to increase in concentration buoyancy effect which it decreased due to increase in either heat absorption coefficient or Schmidt number.
[2] considered heat transfer in viscous flow along a plane wall with periodic suction and heat source.The effects of various parameters on heat transfers in a 3-D laminar boundary layer past a flat plate in the presence of a heat source when a sinusoidal transverse suction velocity was applied to the walls were studied.
[7] discussed their study on natural convection heat transfer in vertical open channel flows with discrete heating.Significant influences on the local heat transfer characteristics were noted, while only slight effects on the average characteristics were observed.[6] Investigated unsteady MHD flow through an inclined closed rectangular channel with upper and lower surfaces of varying permeability.In the same year [7] considered natural convection flow over a vertical frustum of a cone with constant heat flux.The numerical values of the velocity and temperature functions required to calculate the surface skin friction and heat transfer rates were discussed for various values of prandtl number.
Quadri and [1] considered simultaneous heat and mass transfer by natural convection for a vertical semi infinite plate embedded in a fluid saturated porous medium in the presence of wall suction and injection, heat generation or absorption effects, porous medium inertia and thermal dispersion effect.In general, the porous medium thermal dispersion effects increase the temperature of the fluid causing higher flow rate along the surface.However this seems not to be the case in their study as the peak values of the temperature and velocity profile were lowered as porous medium thermal dispersion parameter increases.
[5] investigated coupled heat and mass transfer by natural convection from a horizontal line course in a saturated porous medium.An analysis of boundary layer coupled heat and mass transfer from a line source in saturated porous medium was reported.The equations were numerically integrated for a range of Lewis number and buoyancy ratio, for aiding and opposing buoyancy effects.
The influence of viscous heating dissipation effect in Natural convection flows, shows that the heat transfer rates are reduced by an increase in the dissipation parameter.However the interaction of the radiation with mass transfer of hydromagnetic dynamics dissipative fluid in the presence of heat source has received little attention hence present study is attempted.[4] solved the problem of MHD flow with wavy porous boundary and the influences of heat source, suction and waviness of the boundaries of the flow were numerically analysed.The K-e' model is more reliable and likely to provide a superior performance for flow boundary layers under strong adverse pressure gradients.This was confirmed by comparing their findings with the experimental results.The research work cited above is the motivation of my current work.This

Mathematical Formulations
Taking x and y as the coordinates parallel and normal to a flat plate respectively and u, v and w are taken as the velocities in the x, y, z direction respectively.Since the plate is finite in length and for a two dimensional free convective fluid flow, the physical variables are functions of x, y and t.
The equation of momentum along the x -axis is The equation of momentum along the y -axis is The equation of energy is The equation of concentration is Final set of non-dimensional governing equations Initial conditions are At y = 0, u 0 (0,i) = 1, w 0 (0,i) = 0, θ 0 (0,i) = 1, c 0 (0,i) = 1.y > 0, u 0 (k,i) = 0, w 0 (k,i) = 0, θ 0 = 1, c 0 (k,i) = 0. k > 0 and all i the boundary conditions takes the form At y = 0, u n (0 of the finite difference method, the program is run using smaller values of △y and △x such as 0.0001, 0.00005, 0.000015.

Primary Velocity Profiles
From Figure 2, it is observed that as the Grashof parameter decreases primary velocity profiles increases and as the distance from the plates increase, the curve exhibits an increase in primary velocity profile from curves Gr = 0, Gr = −0.2,Gr = −0.4,Gr = −0.6,Gr = −0.8,Gr = −1 respectively (heating of the plate by free convection current).

Secondary Velocity Profiles
From Figure 3, it is observed that as the Grashof parameter decreases from 0 to -1 the secondary velocities profiles increases but as distance from the plate increases the curve exhibits a decrease in the velocity profile from the curves Gr = 0, Gr = −0.2,Gr = −0.4,Gr = −0.6,Gr = −0.8,Gr = −1 respectively (heating of the plate by free convection current).

Temperatures Profile
From Figure 4, it is observed that as the Grashof parameter decreases from 0 to -1 the temperature profiles increases, also as the distance from the plate increases the curve exhibits an increase in temperature profiles from curves

Concentration Profile
From the Figure 5, it is observed that as the Grashof parameter decreases from 0 to -1 the concentration profiles increases but as distance from the plate increases the concentration profile decreases exponentially as exhibited by the curves Gr = 0, Gr = −0.2,Gr = −0.4,Gr = −0.6,Gr = −0.8,Gr = −1 respectively (heating of the plate by free convection current).

Conclusion
An analysis of the buoyancy effects of thermal and mass diffusion on velocities, temperature and concentration profiles on MHD natural convection past a finite vertical plate has been carried out.In all the cases considered, the velocity was resolved into two components and the work restricted to laminar boundary layer.
In the study the results obtained for Gr < 0 when the temperature of the plate is less than that of the fluid in the free stream region, which implies that the heat was being transferred from the fluid to the plate which refers to the heating of the plate by free convection currents.
The results obtained for various values of Grashof parameter were presented

∞
Nomenclature and List of Symbols C Dimensionless concentration of species kg/m3 C + Dimensional concentration of the species in free stream kg/m3 D Diffusion coefficient m 2 /s Gr Grashof number Re Reynolds number Sc Schmidt number D Gradient operator i ∂ ∂x + j ∂ ∂y − k ∂ ∂z 1. Introduction