^{1}

^{*}

^{2}

A numerical analysis has been carried out to study the problem of plane stagnation double-diffusive MHD convective flow with convective boundary condition in a porous media. The governing nonlinear partial differential equations have been reduced to systems of nonlinear ordinary differential equations by the similarity transformations. The transformed equations are solved numerically by using the classical fourth order Runge-Kutta method together with the shooting technique implemented on a computer program. The effects of the physical parameters are examined on the velocity, temperature and concentration profiles. Numerical data for the skin-friction coefficients, Nusselt and Sherwood numbers have been tabulated for various parametric conditions and are also shown graphically and discussed.

The effect of free convection on the flow of a viscous incompressible fluid past a plane stagnation point on an infinite vertical plate has many important scientific and technological applications in the astrophysical, geophysical and engineering problems. The heating of rooms and buildings by the use of radiators is a familiar example of heat transfer by convection. Heat and mass transfer play an important role in manufacturing industries for the design of fins, steel rolling, nuclear power plants, gas turbines and various propulsion device for aircraft, combustion and furnace design, materials processing, energy utilization, temperature measurements. Most studies have been concerned with constant surface velocity and temperature, but for many practical applications the surface undergoes convective cooling or heating.

The growing interest in the field of flow through porous media stems from the fact that heat and mass transfer through porous media occurs in many engineering, geophysical and biological applications. Convection flows in porous media are particularly important in environmental studies involving air and water pollution. In order to assess the impact of an emission source, it is necessary to predict both the trajectory and the diffusion of the plum or jet at several distances away from the source. In practice, most emissions of a buoyant effluent are accompanied by a discharge momentum flux, hence the need to study the combined effect of heat and mass transfer in a buoyant driven flow.

The subject of Magneto hydrodynamics has attracted the attention of a large number of scientists due to its diverse application. In astrophysics and geophysics, it is applied to study the stellar and solar structures, interstellar matter and radio propagation through the ionosphere. In engineering, it finds its application in MHD pumps and MHD bearing. The study of effects of magnetic field on free convection flow is important in liquid metals, electrolytes, and ionized gases. At the high temperature attained in some engineering devices, gas, for example, can become ionized and so becomes an electrical conductor. The ionized gas or plasma can be made to interact with the magnetic field and alter heat transfer and friction characteristic. Recently, it is of great interest to study the effect of magnetic field on the temperature distribution and heat transfer when the fluid is not only an electrical conductor but also when it is capable of emitting and absorbing thermal radiation.

Heat transfer characteristics in stagnation point flow has also been studied by several authors, for example Hiemenz [

Heat and mass transfer play an important role in manufacturing industries for the design of fins, steel rolling, nuclear power plants, gas turbines and various propulsion device for aircraft, combustion and furnace design, materials processing, energy utilization, temperature measurements. Most studies have been concerned with constant surface velocity and temperature, but for many practical applications the surface undergoes stretching and cooling or heating that cause surface velocity and temperature variations. Crane [

The problem of plane stagnation double diffusive convective flow with convective boundary conditions in a porous media has receive little or no attention by researchers hence the motivation for this study.

A two-dimensional body is placed in a stream of quiescent fluid. We consider heat transfer near the upstream stagnation line, where the flow is assumed to be laminar. The problem is restricted to the case of a plane plate perpendicular to the stream. It is assumed that the fluid properties are constant except for the fluid viscosity which vary as an inverse linear function of temperature as follows.

where u and v are the velocity components in the x and y directions respectively, T is the temperature, is the free stream temperature, is the reference temperature, is the reference concentration and is the free stream concentration. g is the acceleration due to gravity, is the coefficient of expansivity, D is the molecular diffusion coefficient, C is the mass concentration, is the free stream temperature, is the free stream concentration, U(x) is the free stream velocity and S(x) is the body force, K is the permeability, p is the pressure, B is the magnetic field, is the kinematic viscosity, is the specific heat capacity, is the fluid density, is the electrical conductivity, is the thermal diffusivity.

In order to solve Equations (1)-(5) above we introduce the following dimensionless quantities.

Substituting the above variables in Equation (4a) into Equations (1)-(4), we obtain

Subject to

where G_{T} is the Grashof number, G_{C} is the solutal Grashof number, Sc is the Schmidt number, Pr is the Prandtl number, H is the magnetic parameter, Da is the Darcy number, and Bi is the Biot number. And prime denote differentiation with respect to η.

The formulation of the problem that account for, the plane stagnation double diffusive convective flow with convective boundary conditions in a porous media was accomplished. The nonlinear Equations (5)-(8) are solved numerically using the classical fourth order Runge-Kutta method together with the shooting technique implemented on a computer program written in Maple (14). A convenient step size was chosen to obtain the desired accuracy. The values for the velocity, temperature, concentration profiles and the skin-friction coefficients, Nusselt and Sherwood numbers have been obtain and tabulated for various parametric conditions, as presented in

In

Different values of the Solutal Grashof number, velocity profiles in the boundary layer are shown in

The effect of Hartmann number on the velocity profiles is shown in

The temperature profile for different values of the Prandtl number is shown in

The Temperature profiles for various values of the Biot number is depicted in