Soret-Dufour Effects on the MHD Flow and Heat Transfer of Microrotation Fluid over a Nonlinear Stretching Plate in the Presence of Suction

In this work, the Micropolar fluid flow and heat and mass transfer past a horizontal nonlinear stretching sheet through porous medium is studied including the Soret-Dufour effect in the presence of suction. A uniform magnetic field is applied transversely to the direction of the flow. The governing differential equations of the problem have been transformed into a system of non-dimensional differential equations which are solved numerically by Nachtsheim-Swigert iteration technique along with the sixth order Runge-Kutta integration scheme. The velocity, microrotation, temperature and concentration profiles are presented for different parameters. The present problem finds significant applications in hydromagnetic control of conducting polymeric sheets, magnetic materials processing, etc.


Introduction
The natural convection processes involving the combined mechanism of heat and mass transfer are encountered in many natural and industrial transport processes such as hot rolling, wire drawing, spinning of filaments, metal extrusion, crystal growing, continuous casting, glass fiber production, paper production, and polymer processing, etc. 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.Goody [2] considered a neutral fluid.Sakiadis [3] 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 Blasiuss problem of boundary flow over a semi-infinite flat plate due to entrainment of ambient fluid.
Micropolar fluids, distinctly non-Newtonian in nature, are referred to those that contain micro-constituents belonging to a class of complex fluids with nonsymmetrical stress tensor.These fluids respond to micro-rotational motions and spin inertia, and therefore can support couple stress and distributed body torque which are not achievable with the classical Navier-Stokes equations or the viscoelastic flow models.The Micropolar fluid models are developed to make an analysis of the flow characteristics of physiological fluids (blood containing corpuscles), colloidal suspensions, paints, liquid crystal suspensions, concentrated silica particle suspensions, oils containing very fine suspensions, industrial contaminants containing toxic chemicals, lubricants, organic/inorganic hybrid nano-composites and clay which are fabricated by melt intercalation etc. Eringen [4] first designed the study on micropolar fluid making an analysis on the theory of micropolar fluids which provided a mathematical model for non-Newtonian behavior.
Crane [5] noted that usually the sheet is assumed to be inextensible, but situations may arise in the polymer industry in which it is necessary to deal with a stretching plastic sheet.For examples, materials manufactured by aerodynamic extrusion processes and heat-treated materials traveling between a feed roll and a wind-up roll or on a conveyor belt possess the characteristics of a moving continuous stretching surface.Moreover lots of metallurgical processes occupy the system of cooling of continuous strips or filaments by drawing them through a quiescent fluid and that in the process of drawing, these strips are sometimes stretched.
An important matter is that the final product depends to a great extent on the rate of cooling.By drawing such strips in an electrically conducting fluid subjected to a magnetic field, the rate of cooling can be controlled and a final product of desired characteristics can be achieved.The study of heat and mass transfer is necessary for determining the quality of the final product.Sparrow [6] explained a parameter named Rosseland approximation to describe the radiation heat flux in the energy equation in his book.
Moreover, the thermal-diffusion (Soret) effect, for instance, has been utilized for isotope separation, and in mixtures between gases with very light molecular weight (Hz, He) and of medium molecular weight (Nz, air) the diffusion-thermo (Dufour) effect was found to be of a considerable magnitude such that it cannot be ignored, described by Eckert and Drake [17] in their book.Recently plenty of investigators [18][19][20][21] are getting interest work on Soret-Dufour effects.
From the above referenced work and the numerous possible industrial applications of the problem, it is of paramount interest in this study in order to clarify the parametric behavior of magneto-hydrodynamic flow of free convection of a micropolar fluid over a nonlinear stretching sheet in the presence of dynamic effects of suction, thermal-diffusion and diffusion-thermo.

Mathematical Model
We consider the isothermal, steady, laminar, hydromagnetic free convection flow of an incompressible micropolar fluid flowing past a nonlinear stretching sheet coinciding with the plane , the flow being confined in the region .The flow under consideration is also subjected to a strong transverse magnetic field with a constant intensity along the y-axis.
Two equal and opposite forces are introduced along the x-axis so that the surface is stretched keeping the origin fixed.The flow configurations and the coordinate system are shown in Figure 1.We assume that the velocity of a point on a sheet is proportional to its distance from the slit.We assume that all the fluid properties are isotropic and constant.Under the usual boundary layer and Boussinesq approximations, the governing equations for the problem under consideration can be written as follows:   In Equation (2) the Darcian porous drag force term is defined by the term,  [15]).We note that in the viscous shear diffusion term, , the Newtonian kinematic viscosity is now supplemented by the Eringen micropolar vortex viscosity, S. In the present work, we assume that the micro-inertia per unit mass j is a constant.Also, positive or negative n indicate the acceleration and deceleration of the sheet from the extruded slit respectively.Here   , u v are the fluid velocity com- ponents in the x-, y-directions respectively, N is the microrotation, T is the temperature,  the kinematic viscosity, is  is the electric conductivity, g is the acceleration due to gravity,  the volumetric coefficient of thermal expansion, 0 B the uniform magnetic field strength, p k is the Darcy permeability of porous medium, j is the Microinertia per unit mass,  the thermal conductivity of the fluid, c p the specific heat at constant pressure, m D is the chemical molecular diffusivity, T is the Mean fluid temperature and c s is the Concentration susceptibility.the The appropriate boundary conditions suggested by the physical conditions are: 1) on the plate surface at 0 y  : 2) matching with the quiescent free stream as : where the subscripts w and  refer to the wall and boundary layer edge, respectively.The relationship between the microrotation function N and the surface shear u y   is chosen for investigating the effect of different surface conditions for the microrotation of the micropolar fluid elements.The conditions are generally of importance in micropolar boundary layer analysis.When microrotation parameter , we obtain which represents no-spin condition i.e. the microelements in a concentrated particle flow-close to the wall are not able to rotate (Rahman [15]).Finally and  are the constants.

Similarity Solutions
The partial differential Equations (1) to ( 5) are transformed into non-dimensional form by mean of following dimensionless variables Implementing Equation ( 8) into Equations ( 1) to (5) produces the following ordinary differential equations:     and corresponding boundary conditions are reduce to: where the primes denote differentiation with respect to   is the suction parameter.

Skin Friction, Rate of Heat and Mass Transfer
The where is the Reynolds number.And hence the values proportional to the skin-friction coefficient, couple stress, Nusselt number and Sherwood number are and Re 0 ,

Numerical Computation
The numerical solutions of the non-linear differential Equations ( 9) to (12) under the boundary conditions (13) have been performed by applying a shooting method namely Nachtsheim and Swigert [22] iteration technique (guessing the missing values) along with sixth order Runge-Kutta iteration scheme.We have chosen a

Results and Discussions
For the purpose of discussing the results of the flow field represented in the Figure 1, the numerical calculations are presented in the form of non-dimensional velocity, microrotation, temperature and concentration profiles.The values of buoyancy parameter Gr is taken to be both positive to represent cooling of the plate.The parameters are chosen arbitrarily where Pr = 0.71 corresponds physically to air at 20˚C, Pr = 1.0 corresponds to the electrolyte solution such as salt water and Pr = 7.0 corresponds to water, and and 1.0 corresponds to hydrogen, water vapor and methanol respectively at 25˚C and 1 atmosphere.The values of Dufour and Soret numbers are chosen in such a way their production is constant provided that the meat temperature is kept constant as well.f .The numerical re- sults for the velocity, microrotation, temperature and concentration profiles are displayed in Figures 6-33.
Figure 6 shows the effect of vortex viscosity parameter  ( 0.2, 0.5, 1, 1.5   are chosen) on the velocity profiles.From here we see that velocity profiles decrease with the increase of  .Figure 7 demonstrates the effect of  on the microrotation profiles.From this figure it is seen that microrotation increases very evidently with the increase of the vortex viscosity parameter .It is also understood that as the vortex viscosity increases the rotation of the micropolar constituents gets induced in most part of the boundary layer where kinematic viscosity dominates the flow.From Figure 8 it is found that the temperature profiles increase for the increase of   .The effect of vortex viscosity parameter on the concentration profile is not so noteworthy displayed in Figure 9.


The effects of the surface nonlinearity constant n are characterized in the Figures 10-13.At the beginning the velocity profiles decrease with the increase of the value of   1, 2, 3, 4 n n  but far away from the plate they increase after 1.76   displayed in Figure 10. Figure 11 expresses that the microrotation profiles at the beginning increase extensively but at a distance from the plate they overlap and start to decrease very slowly.Figures 12 and 13 enlighten the temperature and the concentration profiles for the increasing influence of the parameter respectively.n Figures 14-17 exhibit the velocity, microrotation, temperature and concentration profiles for the different values of the Eckert number Ec (0.03, 0.1, 0.5 and 1.0).Figure 14 demonstrates that the effect of the Ec on velocity profiles very significant.We observe that velocity increases rapidly with increasing the value of Ec.From the Figure 15 we notice that microrotation profiles decreases with the increase of the value of Ec. Figure 16 presents the increasing effect of Ec on the temperature profiles.The concentration profiles decrease with the increase of the value of Ec illustrated in Figure 17.number Sr influences the concentration profiles to a great extent.Quantitatively when 1.0

 
and Du increases from 0.5 to 1, there is 23.08% decrease in the concentration value, whereas the corresponding decrease is 23.07%when Du increases from 2 to 3.
Figure 30 displays that the suction parameter f w has strong effect on the velocity profiles.With the increase of the value of fw the velocity profiles decrease.Elaborately when 1.0

 
and Du decreases from 0 to 0.5, there is 176.79%decrease in the concentration value, whereas the corresponding decrease is 21.47% when Du increases from 1 to 3. It is observed that, when suction f w increases, the microrotation increase monotonically seen in Figure 31.These Figures 32 and 33 indicate that temperature as well as concentration profiles decrease with the increase of suction velocity or mass transfer parameter frequently.

Conclusions
In the present paper, Soret Dufour effect on the boundary layer flow and heat transfer of microrotation fluid over a nonlinear stretching plate in the presence of suction has been studied.The governing momentum and energy equations were transformed to a set of non linear ordinary differential equations by employing the appropriate similarity transformations and solve numerically for various combinations of problem parameters.From the present study the following conclusions are made: 1) The effect of vortex viscosity parameter on velocity and microrotation is prominent; 7) The effect of suction parameter   fw is dominat- ing on the velocity, microrotation, temperature and concentration profiles.So, using suction boundary layer growth can be stabilized.

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

(
non-dimensional y-coordinate) and S parameters of engineering interest for the present problem are the skin friction coefficient   the wall shear stress, couple stress, the rate of heat transfer and the local surface mass flux respectively.The dimensionless skin-friction coefficient, Couple stress, Nusselt number and Sherwood number for impulsively started plate are given by satisfy the convergence criterion of 10 −6 in all cases.The value of   has been found to each iteration loop by  determined when the values of the unknown boundary conditions at 0   not change to successful loop with error less than 10 −6 .In order to verify the effects of the step size   , we have run the code for our model with three different step sizes as Δη = 0.01, Δη = 0.005 and Δη = 0.001, and in each case we have found excellent agreement among them shown in Figures 2-5 .

Figures 18 -. 5 Figure 6 .
Figures 14-17 exhibit the velocity, microrotation, temperature and concentration profiles for the different values of the Eckert number Ec (0.03, 0.1, 0.5 and 1.0).Figure14demonstrates that the effect of the Ec on velocity profiles very significant.We observe that velocity increases rapidly with increasing the value of Ec.From the Figure15we notice that microrotation profiles decreases with the increase of the value of Ec.Figure16presents the increasing effect of Ec on the temperature profiles.The concentration profiles decrease with the increase of the value of Ec illustrated in Figure17.Figures 18-21 represent the influence of the constant parameter  for the values 0, 1, 2, 5   .All the profiles except microrotation profiles decrease with the increase of .The effects of  are very significant smooth on the distributions.The microrotation profiles increase with the increase of the value of .It is observed from the Figure 22 that with the increase

Figure 9 .
Figure 9. Distribution of concentration profiles for .

Figure 12 .
Figure 12.Distribution of temperature profiles for n.

2 )n 3 ) 5 )
Nonlinearity of the stretching surface is effective on the boundary layer flow;   Effect of Eckert number   Ec is uniform; 4) Boundary layer growth can be controlled by using constant parameter    ; The Dufour   Du effect is significant; 6) The Soret number   Sr plays a role on concentra- tion;