Effects of Thermophoresis on Unsteady MHD Free Convective Heat and Mass Transfer along an Inclined Porous Plate with Heat Generation in Presence of Magnetic Field

An analysis of Thermophoresis effect on unsteady magneto-hydrodynamic free convection flow over an inclined porous plate with time dependent suction in presence of magnetic field with heat generation has been considered by employing Nachtsheim-Swigert shooting iteration technique along with sixth order Runge-Kutta integration scheme. Resulting non-dimensional velocity, temperature and concentration profiles are then presented graphically for different values of the parameters entering into the problem. Finally, the effects of the pertinent parameters on the skin-friction coefficient, the rate of heat transfer (Nusselt number) and wall deposition flux (Stanton number), which are of physical interest, are exhibited in tabular form.


Introduction
The free convection processes involving the combined mechanism of heat and mass transfer are encountered in many natural processes, in many industrial applications and in many chemical processing systems.The study of free convective mass transfer flow has become the object of extensive research as the effects of heat transfer along with mass transfer effects are dominant features in many engineering applications such as rocket nozzles, cooling of nuclear reactors, high sinks in turbine blades, high speed aircrafts and their atmospheric re-entry, chemical devices and process equipments.S. Ostrach [1], the initiator of the study of convection flow, made a technical note on the similarity solution of transient free convection flow past a semi infinite vertical plate by an integral method.B. C. Sakiadis [2] analyzed the boundary layer flow over a solid surface moving with a constant velocity.This boundary layer flow situation is quite different from the classical Blasius problem of boundary flow over a semi-infinite flat plate due to entrainment of ambient fluid.L. E. Erickson, L. T. Fan and V. G. Fox [3] extended the work of Sakiadis for suction or injection of a smooth surface.S. L. Goren [4] studied thermophoresis in laminar flow over a horizontal flat plate.He found the deposition of particles on cold plate and particles free layer thickness in hot plate case.E. M. Sparrow [5] explained a parameter named Rosseland approximation to describe the radiation heat flux in the energy equation in his book.R. S. R. Gorla [6] depicted the application of linearly stretched surface in electrochemistry.G. M. Homsy, F. T. Geyling and K. L. Walker [7] solved Blasius series solution.A. Raptis and C. Perdikis [8] studied numerically free convection flow through a porous medium bounded by a semi-infinite vertical porous plate.M. Epstein, G. M. Hauser and R. E. Henry [9] analyzed the thermophoresis in natural convection for a cold vertical surface.M. A. Alabraba, A. R. Bestman and A. Ogulu [10] studied the interaction of mixed convection with thermal radiation in laminar boundary layer flow taking into account the binary chemical reaction and Soret-Dufour effects.A. Sattar and M. Hossain [11] investigated the unsteady free convective flow, with Hall currents and mass transfer, past an accelerated vertical porous plate in the presence of a transverse magnetic field while assuming the plate temperature and concentration to be functions of time.M. A. Hossain and H. S. Takhar [12] analyzed the effect of radiation using the Rosseland diffusion approximation that leads to non similar boundary layer equation gov-erning the mixed convection flow of an optically dense viscous incompressible fluid past heated vertical plate with a free uniform stream velocity and surface temperature.
In recent years, it is found that the thermophoresis phenomenon has many practical applications in removing small particles from gas streams, in determining exhaust gas particle trajectories from combustion devices, and in studying the particulate material deposition on turbine blades.It has been found that thermophoresis is the dominant mass transfer mechanism in the modified chemical vapor deposition (MCVD) process as currently used in the fabrication of optical fiber performs.Thermophoretic deposition of radioactive particles is considered to be one of the important factors causing accidents in nuclear reactors.J.-S.Lin, C.-J. Tsai and C.-P. Chang [13] investigated the suppression of particle deposition from flow through a tube with circular cross-section when the wall temperature exceeds that of the gas.M. A. Seddek [14] investigated finite element method for the effects of thermophoresis and chemical reaction on a boundary layer hydro magnetic flow with heat and mass transfer over a heat surface and found that the velocity, temperature and concentration profiles reduce with the increase of thermophoretic parameter.M. S. Alam, M. M. Rahman and M. A. Sattar [15] considered the effects of heat generation and thermophoresis on steady, laminar, hydromagnetic, two-dimensional flow with heat and mass transfer along a semi-infinite, permeable inclined flat surface.M. A. Samad and M. E. Karim [16] considered an unsteady two dimensional MHD flow of a viscous incompressible and electrically conducting fluid along a vertical plate with time dependent suction under the influence of a uniform magnetic field.M. A. A. Mahmoud [17] made a report on the effects of thermal radiation and variable viscosity on the unsteady hydromagnetic flow of an electrically conducting fluid past an infinite vertical porous plate in the presence of viscous dissipation and time dependent suction.M. A. Samad, M. E. Karim and D. Mohammad [18] calculated numerically the effect of thermal radiation on steady MHD free convectoin flow taking into account the Rosseland diffusion approximaion.P. Loganathan and P. P. Arasu [19] analyzed the effects of thermophoresis particle deposition on non-Darcy MHD mixed convective heat and mass transfer past a porous wedge in the presence of suction or injection.M. Ferdows, Nazmul and M. OTA [20] descried that in the presence of uniform magnetic field with viscous dissipation at the wall, the thermophoretic parameter is one of the most useful parameter to control the boundary layer of the fluid.O. D. Makinde and P. O. Olanrewaju [21] investigated unsteady mixed convection with the thermo-diffusion and diffusion-thermal effects past a vertical porous plate moving through a binary mixture in the presence of radiative heat transfer and n-th order Arrhenius-type chemical reaction.N. Ghara, S. L. Maji, S. Das, R. Jana and S. K. Ghosh [22] analyzed the unsteady MHD Couette flow of a viscous fluid between two infinite non-conducting horizontal porous plates with the consideration of both Hall currents and ion-slip.
The present paper is the investigation of the thermophoretic effects of an electrically conducting viscous incompressible fluid interaction with heat generation on the flow over an inclined porous plate in the presence of heat and mass transfer permitted by a transversely applied uniform magnetic field taking into account the Rosseland diffusion approximation.The investigation is based on known similarity analysis and the local similarity solutions are obtained numerically.

Mathematical Analysis
In this work we considered an unsteady two-dimensional MHD flow of a viscous incompressible and electrically conducting fluid past an inclined porous plate with an angle  to the vertical embedded in a porous medium under the influence of a uniform magnetic field.Initially the flow is assumed to be in the x-direction, which is chosen along the plate in the upward direction and y-axis normal to it.The plate and the fluid are at a constant temperature T  in a stationary condition with concentration level C  at all points.At time t > 0 the plate is assumed to be moving in the upward direction with the velocity   U t and there is a suction velocity   0 v t taken to be a function of time, the temperature of the plate raised to   T t and the concentration level at the plate is raised to   C t where and The plate is considered to be of infinite length, all derivatives with respect to x vanish and so the physical variables are functions of y and t only.The flow configuration and the coordinate system are shown in the Figure 1.
The fluid is assumed to have constant properties except that the influence of the density variations with temperature and concentration, which are considered only in the body force term, and is considered to be gray, absorbing emitting radiation but non-scattering medium, and the Rosseland approximation is used to describe the radioactive heat flux in the energy equation.A uniform magnetic field of strength B 0 is applied normal to the plate parallel to y-direction.
Under the usual boundary layer and Boussinesq approximation and using the Darcy-Forchhemier model, the flow and heat transfer in the presence of radiation are governed by the following equations.
Continuity Equation

Momentum Equation
Energy Equation where u and v are the velocity components along x-and y-directions respectively, υ is the kinematic viscosity, ρ is the density of the fluid, g 0 is the acceleration due to gravity, β is the coefficient of volume expansion, β * is the volumetric coefficient of expansion with concentration, α is the inclination of the plate, σ is the electric conductivity, B 0 is the uniform magnetic field strength (magnetic induction), k is the Darcy permeability, b is the empirical constant, c p is the specific heat at constant pressure, T and T ∞ are the fluid temperature within the boundary layer and in the free-stream respectively, while C and C ∞ are the corresponding concentrations.Also Q 0 is the heat generation, D m is the coefficient of mass diffusivity and V T is the thermophoretic velocity.
The effects of thermophoresis are being taken into account to help in the understanding of the mass deposition variation on the surface.We further assume that 1) the mass flux of particles is sufficiently small so that the main stream velocity and temperature fields are not affected by the thermo-physical processes experienced by the relatively small number of particles; 2) due to the boundary layer behavior the temperature gradient in the y-direction is much larger than that in the x-direction and hence only the thermophoretic velocity component which is normal to the surface is of importance; 3) the fluid has constant kinematic viscosity and thermal diffusivity, and that the Boussinesq approximation may be adopted for unsteady laminar flow; 4) the particle diffusivity is assumed to be constant, and the concentration of particles is sufficiently dilute to assume that particle coagulation in the boundary layer is negligible; and 5) the magnetic Reynolds number is assumed to be small so that the induced magnetic field is negligible in comparison to the applied magnetic field.
Initially (t = 0) the fluid and the plate are at rest.Thus the no slip boundary conditions at the surface of the plate for the above problem for t > 0 are: The effect of thermophoresis is usually prescribed by means of an average velocity that a particle will acquire when exposed to a temperature gradient.For boundary layer analysis it is found that the temperature gradient along the plate is much lower than the temperature gradient normal to the surface, i.e., T y x T .So the component of thermophoretic velocity along the plate is negligible compared to the component of its normal to the surface.As a result, the thermophoretic velocity V T which appears in Equation ( 4) can be written as where K is the thermophoretic coefficient which ranges in value from 0.2 to 1.2 as indicated by G. K. Batchelor and C. Shen [23] and is defined from the theory of L. Talbot, R. K. Cheng, A. W. Schefer and D. R. Wills [24] by: where  and p  are the thermal conductivities of the fluid and dif- fused particles, respectively and Kn is the Knudsen number.
A thermophoretic parameter  can be defined (see A.
With this similarity parameter , a similarity variable is then introduced as δ y   (10) In terms of this length scale, a convenient solution of the Equation ( 1) can be taken as where v 0 is the mass transfer parameter, which is positive for suction and negative for injection.Following M. M. A. Samad and M. M. Rahman [27], we see that

 
U t , and are now considered to have the following form: where n is a non-negative integer and, U 0 , T 0 and C 0 are respectively the free stream velocity, mean temperature and concentration.Here , is the value of  at t = t 0 .0 δ Now to make the Equations ( 2)-( 4) dimensionless, we introduce the following transformations: Now introducing all the above similarity variables in Equations ( 2) and (3) we have the following dimensionless ordinary non-linear differential equations where is the modified Grashof number, As a result the corresponding boundary conditions for t > 0 take the form 1, 1, 1 at 0 0, 0, 0 as

Skin-Friction Coefficient, Nusselt Number and Stanton Number
The parameters of engineering interest for the present problem are the skin-friction coefficient, local Nusselt number and the local Stanton number which indicate physically wall shear stress, rate of heat transfer and wall deposition flux respectively.The skin-friction coefficient is given by   the local Nusselt number may be written as and the Stanton number can be written as Thus the values proportional to the skin-friction coefficient, Nusselt number and the Stanton number are   0 f  , and respectively.

Numerical Computation
The numerical solutions of the nonlinear differential Equations ( 14)-( 16) under the boundary condition (17) have been performed by applying a shooting method namely P. R. Nachtsheim and P. Swigert [28] iteration technique (guessing the missing value) along with sixth order Runge-Kutta integration scheme.We have chosen a step size of  = 0.01 to satisfy the convergence criterion of 10 −6 in all cases.The value of   was found to each iteration loop by   =  + .The maximum value of   to each group of parameters v 0 , Gr, Gm, α, M, Da, Fs 1 , n, Pr, Q, Sc and τ determined when the value of the unknown boundary conditions at  = 0 not change to successful loop with error less than 10 −6 .

Results and Discussion
For the purpose of discussing the results, the numerical calculations are presented in the form of non-dimensional velocity and temperature profiles.Numerical computations have been carried out for different values of the suction parameter (v 0 ), Grashof number (Gr), modified Grashof number (Gm), angle of the inclination    , magnetic field parameter (M), Darcy number (Da), modified Forchhemier number (Fs 1 ), constant parameter (n), Prandtl number (Pr), heat generation (Q), Schmidt number (Sc) and the thermophoresis parameter (τ).The values of Grashof number (Gr) are taken be large from  the physical point of view.The large Grashof number values correspond to free convection problem.
Figures 5-7 display the effects of the suction parameter v 0 on the velocity, temperature and concentration profiles respectively.It is observed that, when suction v 0 increases, all the profiles decrease.The effects on velocity and temperature profiles are significant but on concentration are small.The decreasing suction restrains the heat transfer coefficient.
We take α = 0˚, 30˚, 40˚ and 60˚ as the values of the angle of the inclination of the plate to explain the effect of α.As α increases, the effect of the buoyancy force decreases because of the multiplication factor cosα and hence the velocities decrease with the increase of α swiftly shown in Figure 8. Figure 9 explains that the temperatures as well as the concentration profiles are completely beyond the influence of the inclination of the plate.M = 0, 0.5, 1.0, 2.0 are taken to examine the effect of the magnetic field parameter M on the velocity, temperature and concentration fields shown in Figures 10   17 and 18.The positive value of Q represents source i.e. heat generation in the fluid.From Figure 17, it is observed that due to the generation of heat the buoyancy force increases which in turn gives higher velocity in the boundary layer.This is corroborated by Figure 18 where it is seen that the temperatures indicated by solid lines do indeed rapidly increase as Q increases and there are very tiny fluctuations in the concentration profiles pointed out by dashed lines.
In Figures 19 and 20 we have plotted th menonless velocity, temperature and concentration profiles showing the effect of thermophoretic parameter τ.The dropping effect of τ on velocity is seen in Figure 19.There is no effect on temperature profiles marked by solid line for the variation of the values of τ in Figure 20.We also observe from Figure 20 specified by the dashed lines that the thermophoretic parameter τ affects the concentration profiles very considerably.The concentration profiles decrease with the increase of τ.

Conclusions
In this paper we ha action with unsteady MHD boundary layer flow past a continuously moving vertical plate with suction.From the present study we can make the following conclusions:  Using trolled.Suction stabilizes the hydrodynamic, thermal

Figure 1 .
Figure 1.Flow configuration and coordinate system.

Figure 2 .
Figure 2. Velocity profiles for different values of .

Figure 3 .Figure 4 .
Figure 3. Temperature profiles for different values of .

Figure 10 .Figure 11 .Figure 12 .Figure 13 . 5 Figure 14 . 5 Figure 15 . 5 Figure 16 .
Figure 10.Velocity profiles for different values of M.of the Darcy number Da on the velocity, temperature and concentration fields shown in Figures12 and 13expressing that the velocity increases slightly with the increase of Da but no effect is shown on the temperature

Finally, the effe
iction C f , local Nusselt number Nu and local Sherwood number Sh are shown in Tables 1-5.

f , Nu an for differen es of Q.
Q2. C