Effects of Viscous Dissipation on Unsteady MHD Thermo Bioconvection Boundary Layer Flow of a Nanofluid Containing Gyrotactic Microorganisms along a Stretching Sheet

This paper presents a numerical study of the problem of unsteady thermo bioconvection boundary layer flow of a nanofluid containing gyrotactic microorganisms along a stretching sheet under the influence of magnetic field and viscous dissipation. With the help of usual transformation, the governing equations are transformed into unsteady nonlinear coupled partial differential equations. The numerical solution is obtained by using an explicit finite difference scheme. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. From the results it is found that both magnetic parameter and bioconvection Rayleigh number have positive effect on the dimensionless Nusselt number and density number of the motile microorganisms while the opposite behavior became clear in the case of Grashof number and Eckert number. The rescaled velocity, temperature, concentration and the density of motile microorganisms depend strongly on the governing parameters.


Introduction
polymer solution and Bio-fluids), the thermal conductivity of these fluids plays a crucial role on the heat transfer coefficient between the heat transfer medium and the heat transfer surface.Nanoparticles are the particles (that have the random movement in the nanofluid called Brownian motion) with diameter less than 100 nm.By pointing out that the use of solid particles as an additive suspended into the base fluids is a manner for increasing the heat transfer rate.Thus, nanofluid is a highly effective method of enhancing heat transfer.Gorla and Hossain [1], Makinde [2] and Mustafa et al. [3] published papers on nanofluids.Kuznetsov and Nield [4] studied analytically the free convective boundary-layer flow of a nanofluid past a vertical plate.They found that the reduced Nusselt number is a decreasing function of each of bouncy force parameter, Brownian motion parameter and thermophoresis parameter.The problem of boundary-layer flow of a nanofluid past a stretching sheet is investigated numerically [5].
They observed that the reduced Nusselt number is a decreasing function of each dimensionless number, while the reduced Sherwood number is an increasing function of higher r P and a decreasing function of lower r P number for each Le , Nb and Nt numbers.Bachok et al. [6] examined boundary layer flow of nanofluids over a moving surface in a flowing fluid.Their results indicate that dual solutions exist when the plate and the free stream move in the opposite directions.Applications of nanofluid: electronics cooling (Nanofluids have been considered as the working fluid for heat pipes in electronic cooling applications), transportation (The addition of nanoparticles to the standard engine coolant has the potential to improve automotive and heavy-duty engine cooling rates), surface coating, biomedical, etc.Sun and Pop [7] studied numerically the steady-state natural convection heat transfer behavior of nanofluid inside a right-angle triangular enclosure saturated by a porous media.It is found that the heat transfer in the cavity is improved with the increasing of solid volume fraction parameter of nanofluids at low Rayleigh number, but opposite effects appear when the Rayleigh number is high.Ferdows et al. [8] presented the problem of transient mixed convective laminar boundary layer flow of an incompressible, viscous, dissipative, electrically conducting nanofluid from a continuously stretching permeable surface in the presence of magnetic field and thermal radiation flux.The steady boundary layer flow of nanofluid over an exponential stretching surface is investigated analytically by Sohail and Changhoon [9].Na and Pop [10] analyzed an unsteady flow due to a stretching sheet.
MHD (Magneto-hydrodynamics) is the science of the motion of electrically conducting fluids under the influence of applied magnetic forces [11].MHD boundary layer flow problem of a nanofluid through a porous medium over an exponentially stretching sheet was studied by Ferdows et al. [12].Heat and mass transfer in the boundary-layer flow of unsteady viscous nanofluid along a vertical stretching sheet in the presence of magnetic field, thermal radiation, heat generation, and chemical reaction are presented by Eshetu and Shankar [13].They found that the velocity, temperature, and concentration profiles of the unsteady flow are less than the corresponding parts of the steady state flow scenario.Recently, Magnetic field effect in three-dimen-sional flow of an Oldroyd-B nanofluid over a radiative surface is presented by Shehzad et al. [14].
Bioconvection is a phenomenon in which physical laws that govern smaller scales lead to a phenomenon visible on a larger scale [15].The problem of natural convection boundary layer flow about a vertical cone in porous media saturated by a nanofluid due to gyrotactic microorganisms is presented by Mahdy [16].A numerical study of a natural convection about a vertical cone embedded in a non Darcian nanofluid containing gyrotactic microorganisms saturated porous medium is studied by Hady et al. [17].
Bioconvection concerns with suspensions of self-propelled microorganisms [18].Bioconvection in a suspension of gyrotactic motile microorganisms is investigated by Kuznetsov and Avramenko [19].The effect of small particles (that are heavier than water) on the stability of a suspension of motile gyrotactic microorganisms in a horizontal fluid layer of finite depth is investigated by Kuznetsov and Avramenko [16].The stability of thermo-bioconvection of oxytactic bacteria in a porous medium is investigated numerically using a Galerkin method by Kuznetsov [20].Nanofluid bioconvection is generated by the combined effects of buoyancy forces and magnetic field on the interaction of motile microorganisms and nanoparticles [21].Khan and Makinde [22] used Oberbeck-Boussinesq approximation and similarity transformations to investigate MHD laminar boundary layer flow with heat and mass transfer of an electrically conducting water-based nanofluid containing gyrotactic microorganisms along a convectively heated stretching sheet.Recently, Computational investigation of Stefan blowing and multiple-slip effects on buoyancy-driven bioconvection nanofluid flow with microorganisms was studied by Jashim Uddin et al. [23].The unsteady flow of liquid containing nanoparticles and motile gyrotactic microorganisms between two parallel plates while keeping one moving and other fixed is described by Ammarah et al. [24].Finally, the main objective of the present paper is to study the effect of viscous dissipation on unsteady thermo bioconvection boundary layer flow of water based nanofluid containing gyrotactic microorganisms along a stretching surface with the influence of magnetic field.

Mathematical Analysis
Consider a tow-dimensional unsteady thermo bioconvection laminar boundary layer flow of a viscous incompressible, electrically conducting nanofluid containing gyrotactic microorganisms past a vertical stretching sheet under the influence of a uniform transversely applied magnetic felid and viscous dissipation.Bioconvection induced flow only take place in a dilute suspension of nanoparticles.The presence of nanoparticles is assumed to have no effect on the direction in which microorganisms swim and on their swimming velocity.We choose the coordinate system ( ) , x y with the x -axis meas- ured along the stretching sheet in the upward direction and the y -axis measured in the normal direction to the stretching sheet.The origin 0 of the coordinate system is placed fixed.The physical model and coordinate system is shown in Figure 1.From the figure, the sheet is stretched due to two equal and opposite forces are introduced along Figure 1.Flow model and physical coordinate system.the x -axis.Also, the figure shows that a strong magnetic field of strength 0 B is ap- plied in the y direction.w T and w n are the temperature and the density of the mo- tile microorganisms at the wall, respectively, which are kept constant thereafter, and the nanofluid particle fraction on the boundary is passively rather than actively controlled.
T ∞ , ∅ ∞ and n ∞ are the ambient values of the temperature, nanoparticle volume fraction, and density of motile microorganisms, respectively, far away from the plate.
With knowing that, the plate temperature and the density of motile microorganisms are Using the Oberbeck-Boussinesq approximation, the boundary-layer approximations of the continuity, momentum, energy, nanoparticle concentration and conservation for microorganisms equations are: ) where t is time, u and v are the velocity components, x and y are the Carte- sian coordinates, T is the temperature, ∅ is the nanoparticle volume fraction, n is the density of motile microorganisms, υ is the kinematic viscosity, f ρ is the density of the fluid, P ρ is the density of the nanoparticles, m ρ ∞ is the microorganism densi- ty, δ is the fluid electrical conductivity, 0 B is the strength of magnetic field, g is the acceleration due to gravity, β is the volumetric expansion coefficient, γ is the average volume of a micro-organism, k is the thermal conductivity, P C is the spe- cific heat at constant pressure, T D is the thermophoresis diffusion coefficient, B D is the Brownian diffusion coefficient, n D is the diffusivity of microorganisms, b is the chemotaxis constant, and c W is maximum cell swimming speed.
The initial and boundary conditions are: , , , everywhere 0, , 0, 0, 0, 0, , , , at 0 By substituting the following dimensionless quantities: into Equations ( 1)-( 5) gives the following dimensionless equations: ( ) ( ) where the Grashof number r G is defined as ( ) ( ) bouncy ratio parameter Nr is defined as )( ) , the magnetic parameter M is defined as , the Prandtl number r P is defined as , the bioconvection constant σ is defined as The dimensionless form of the initial and boundary conditions are: 0, 0, 0, 0, 0, χ 0 0, 0, 0, 0, χ 0 eve at 0 rywhere 0, The dimensionless form of the local Nusselt number, density number of the motile microorganisms and Skin-friction coefficient, respectively is defined as: where, the local Reynolds number Re x is defined as 0 Re x U x υ = .

Numerical Method
The explicit finite difference method is used to solve the unsteady nonlinear coupled partial differential Equations ( 8)-( 12) with initial and boundary conditions ( 13) and (14).The dimensionless equations are solved for the dependent variables U , V , θ , E and χ as functions of X , Y and t .Here, max 1  8)-( 12) are: subject to initial and boundary conditions , 0, 0, 0, 0 ,0 ,0 ,0 , 0, 0, 0, 0 where where ( ) where ( ) , i j represent the grid points with X and Y coordinates, respectively, and n represents the value of time.The value of time t which represent the steady state condition in the figures is 1 t = .In order to check the accuracy of the present method, the obtained numerical results are compared with previously published works by Khan et al. [25].
Table 1 shows that the present results are in excellent agreement with the results presented by Khan et al. [25].

Results and Discussion
In the following section, the numerical solutions of our problem are discussed and displayed by graphic.The effects of the governing physical parameters on the dimensionlesss velocity, temperature, nanoparticle volume fraction and density of motile microorganisms are illustrated in Figures 2-13.on the velocity distribution.It is found that an increase in M and Nt leads to reduction of the velocity profiles and the velocity profiles increase with increasing r P .The physical meaning of the behavior in Figure 2 is due to the Lorentz force which results from the presence of a magnetic filed in the nanofluid and it works to slow down the velocity.Figure 5 and Nanofluids are produced by suspending nanoparticles made up of metal carbides, nitrides and carbon nanotubes in heat transferring fluids (e.g.water, oil, ethylene glycol, How to cite this paper: Hady, F.M., Mahdy, A., Mohamed, R.A. and Zaid, O.A.A. (2016) Effects of Viscous Dissipation on Unsteady MHD Thermo Bioconvection Boundary Layer Flow of a Nanofluid Containing Gyrotactic Microorganisms along a Stretching Sheet.World Journal of Mechanics, 6, 505-526.
Lewis number Le is defined as B Le D = υ , the bioconvection Lewis number Lb is defined as the height of the plate and max 15 Y = is the boundary layer thickness.It is assumed that 0sizes along the X and Y directions, respectively and 0.0001 t ∆ = , these values gives accurate results.In the case of U t the unsteady state problem, the steady state condition is obtained.Let U ′ , V ′ , θ ′ , E′ and χ ′ denote the values of U , V , θ , E and χ at the end of a time step.The finite difference equations corresponding to Equations ( Re x Nu for different values of Nt and Nb when 0 Gr Nr Rb M Ec Pe Lb σ

Figure 2 .
Figure 2. Effect of magnetic parameter M on velocity profiles.

Figure 3 .
Figure 3.Effect of thermophoresis parameter Nt on velocity profiles.

Figure 6
Figure 6 display the temperature distribution for different values of Eckert number ( ) 0.01 0.9 c E ≤ ≤