Finite Element Simulation of Multi-Slip Effects on Unsteady MHD Bioconvective Micropolar Nanofluid Flow Over a Sheet with Solutal and Thermal Convective Boundary Conditions

In this article, the intention is to explore the flow of a magneto-hydrodynamic (MHD) bioconvective micro-polar Nanofluid restraining microorganism. The numerical solution of 2-D laminar bioconvective boundary layer flow of micro-polar nanofluids are presented. The phenomena of multi-slip, convective thermal and Solutal boundary conditions have been integrated. A system of non-linear partial differential equations are transformed into the system of coupled nonlinear ordinary differential equations by applying appropriate transformations, the transformed equations are then solved by applying the variational finite element method (FEM). The fascinating features of assorted velocity parameter, microrotation, temperature, microorganism compactness, solutal and nanoparticles concentration have been inspected. The rate of heat transfer, the skin friction coefficient, couple stress and Sherwood number for microorganisms have also been discussed graphically and numerically. The investigations illustrated that increase in material parameters causes a reduction in microorganism compactness, concentration and temperature. As a result of enhancement in the unsteadiness parameter, the fluid velocity, concentration of microorganisms and the temperature are observed to be declines. Energy and microorganism compactness profile affected by the improvement in the buoyancy ratio parameter. As the improvement in results of buoyancy ratio parameter effects on improvement in the energy and the microorganism compactness profile while the velocity profile is condensed. In the end, rationalized convergence of the finite element solution has been inspected; the computations are found out via depreciating the mesh size.


Introduction
The transfer of heat is the basic feature of the gigantic appliance through applications. The heat transfer depends on the thermal conductivity of the functioning fluids as measures up to the aptitude of thermal appliances and the structure. Nanofluids are deferments of nanoparticles in fluids that demonstrate the momentous enrichment of their properties at retiring nanoparticle concentrations. There are numerous publications on nanofluids that concerned to understand their behavior so that it can be employed wherever the straight heat transfer enrichment is dominant as in various industrial applications, nuclear reactors, transportation, electronics as well as biomedicine and food. nanofluid has been testified, where the heat transfer can be condensed or increased. Wen [1] represents about the inadequate considerate of the construction and mechanism of nanofluids and their applications. Exploration of nanofluids has proceeded quickly since its enhanced thermal conductivity was first reported about a decade ago, though much controversy and inconsistency have been reported. Yang et al. [2] examined the viability of the new application, the stable TiO 2 nanofluids constructed to examine the formaldehyde and total volatile organic compounds (TVOC) distillation as well as sterilization. They also investigated the heat recovery characteristics of nanofluids by testing the gas-liquid aerated heat transfer behavior. The learning of magnetohydrodynamics plays an imperative role as it has importance in industry and much another flow phenomenon. From the last few decades, the endeavor to explore fluid physics at the nano and micro scales has improved. The fluid flow in the micro-scale is different from that on the macroscale. Areas of biological, microcirculation, molecular machinery, atherogenesis and micro-fluidics have motivated researchers to work on micro-scale phenomena in the field of fluid mechanics. It is desirable to investigate the large scale applications of non-Newtonian fluids in different industries, such as molten plastics, imitation fibers, nuclear desecrate disposal, slurries, transpiration cooling and petroleum tanks, and their characteristic performance. The heat transfer, resolute by the thermal radiation, has huge applications in diverse scientific processes, nuclear power plants, space vehicles, gas turbines, satellites, comprising missiles and in many impulsive devices for jets. Bio-convection arrangement which is a common phenomenon, they generally take place due to up-swimming of germs (tiny living organism like as microbes, microorganisms that are observable in microscopes) that is slightly impenetrable. By up-swimming, the top portion of the liquid layers of postponement becomes too dense. Due to the assembly of microorganisms, it happens to wobbly. In the geophysical phenomenon, Thermo bio-convection plays an imperative role, resembling in the phenomenon of thermopiles in which motile microorganisms roam warm springs. Moreover in the field of oil convalescence enhances by microbial. In oil-bearing layers are further microorganisms and nutrients to systematize permeability discrepancy. Bio-convective flows with magnetohydrodynamics (MHD) have quite a lot of applications in different fields such as hyperthermia, Treatment of several arterial syndromes, the illustration of copper wires and diminution of blood in surgeries and many other fields. A Casson fluid has an inestimable thickness at zero rates of shear and has zero thickness at an inestimable rate of shear and relent stress in which no flow happens. Therefore, Casson-fluid has to be described as shear contraction fluid/pseudoplastic fluid. The flow features of blood can precisely explain by Casson fluid at small shear rates. Pradhan et al. [3] Investigation are conceded out for free convective flow of an electrically executed micropolar fluid over an absorbent able stretching sheet in the existence of porous medium. Sharada and Shankar [4] numerically investigated the consequences of Soret/Dufour, chemical reaction, thermal radiation on the fluid flow, mass and heat transfers of Casson fluid through the exponentially elongating surface. Moradikazerouni [5] scrutinized the numerical heat transfer and thermal deformation of a micro-channel heat sink with various ingress channels. Pushpalatha et al. [6] investigate that the thermos-diffusion and thermos-diffusion impact on Casson fluid over an unsteady stretching surface with thermal radiation and magnetic field. Also, it described that Casson-parameter and unsteadiness parameter contains the proclivity to disparage the velocity distribution. Vo et al. [7] represents the influences of numerous shapes of γ-AlOOH nanoadditives (i.e., cylindrical, brick, blade and platelet) on the fluid flow and the heat transfer aspects of a sinusoidal grooved channel. Also, scrutinized the impact of nanoadditive fraction and Reynolds number on the performance of nanofluid in the grooved. Ullah et al. [8] analyzed the buoyancy-driven unsteady flow of Casson nanofluid. Also, see the consequences of Brownian motion and thermophoresis on the flow fields that are deliberated in the existence of a magnetic field. Uddin et al. [9] deliberated about free convective boundary layer flow of nanofluid over a penetrable vertical plate tinplate porous medium along with thermal convective boundary conditions. Alsarraf et al. [10] represents that a two-phase mixture model is used to discuss the turbulent boehmite alumina nanofluid inside a double-pipe heat exchanger and appraise the effect of nanoparticle concentration and Reynolds number on the heat exchanger performance. Moradikazerouni [11] investigated the influence of heat sink thermal demonstration and of the environmental parameters on the heat sink's stability against thermal stress, which are important for manufacturers' targets. They also inspected the effects of parameters on heat the transfer rate from the surface of the heat sink. Ranjbarzadeh et al. [12] scrutinized experimentally that the improvement of thermal properties of deionized water (base fluid), stability and the thermal conductivity of the silica-water nanofluid. The thermal conductivity of the organized nanofluid has been experimentally inspected in the wide range of temperatures and solid volume fractions. Naseem et al. [13] investigated the study of MHD bio-convection non-Newtonian nanofluid flow over a stretching sheet and conclude that, the motile concentration profile compact for greater Péclet and bio-convection Lewis number. It is acknowledged that in space, the gravity effect is decreased as a result of both the thermal buoyancy consequence and solutal buoyancy consequence are also decreased. The g-jitter/lingering accelerations originate commencing the diversity of resources like crew movements, perfunctory vibrations like as propels, motors, spacecraft exercises, distinctive drag and the gradient of earth's weightiness [14]. The movement of macroscopic convective fluid flow due to the compactness is greater than water, which is entitled bio-convection [15,16]. It is valuable in fuel cell machinery, biodiesel fuels, bio-reactors, bio-micro organisms and gas-bearing, and so forth. The advantages of summing up the motile microorganisms to nanoparticles are enhancing the mass transporting and perk up the constancy of nanofluid. Dissimilar bio-convection nanofluid sculpts were studied based on various kinds of microorganisms and in the motion system of the fluid. Kuznetsov [17] examined that a horizontal layer of nanofluid that restrain gyro tactic microorganisms, then the particles may either lift up or sadden critical Rayleigh numbers, while the microorganisms have always subverting effect. In view of oxytactic motile microorganisms, bio-convection in nanofluid drenched in porous media, he also notorious that the nanofluid constancy is stated by nanoparticles allocation, compactness stratification stimulated by a perpendicular temperature gradient, by up-swimming of oxytactic microorganisms. Reference [18] revealed that the existence of oxytactic microorganisms reduces the nanofluid postponement consistency and eradicates the oscillatory inconsistency as compare to non-oscillatory inconsistency. Recently, Rahman et al. [19] investigated that the hydro-magnetic multi slip flow of water-based nanofluid through a wedge along with convective surface in the subsistence of heat production or assimilation. Ibrahim and Shankar [20] evaluated the magneto-hydrodynamic boundary layer flow and the heat transfer of nanofluid past through penetrable stretching sheet using thermal, velocity and solutes slip boundary conditions. Ma et al. [21] scrutinized that the heat transference of alumina/water nanofluid in a square inclusion through an inclined magnetic field and also investigate the impact of Ra, Hartman number, magnetic field angle, enclosure angle, the volume fraction of nanoadditives and aspect ratio on the heat transfer. Das [22] obtained a very important numerical inquiry on the convective heat transfer recital of nanofluids over a penetrable stretching surface in the existence of thermal buoyancy, partial slip and internal heat production/fascination. Abbas et al. [23] calculated the impacts of radiation in the existence of the same magnetic field for the nanofluid on a curled stretching surface by assimilating the consequences of slip. In recent times, Makinde et al. [24] have seen numerically the radiation consequences on chemically reacting MHD nanofluid manipulated by the heat source/sink and also illustrated the radiation effects in collective heat and mass relocating investigation for mixed convection flow in excess of the perpendicular plane with radiation and chemical reaction. Aziz et al. [25] using the Buongiorno model [26] deliberated the boundary layer flow through a porous medium crammed with nanofluid and gyrotactic microorganism. They give a correspondence solution and illustrated that the influences of the dimensionless parameters on the behavior of motile microbes. Several experimental consequences have given the verifications to support the slip condition [27]. In recent publications [28][29][30][31][32][33] investigated the significant significance of nanofluids in distinctive types of microsystems. The fluid motion of macroscopic convective flow caused by density, which is vaguely higher than water called bio-convection. It is valuable in fuel cell equipment, bio-micro systems, bio-reactors, bio-diesel fuels and gas-bearing sedimentary, and so forth. The reimbursement of accumulation motile microorganisms to nanoparticles is to improve the mass transfer and get better the stability of nanofluid. Different bio-convection nanofluid models were deliberated on the basis of different categories of microorganisms and the directional motion system of the fluid.
Inspired by the above literature, and in the applications of numerous areas that have been discussed, an investigation of the impact of multi-slip and solutal boundary conditions on MHD unsteady bioconvective micropolar nanofluid restraining gyrotactic microorganism, heat and mass transfer effect over a stretching/shrinking sheet (which have not been discussed before) was carried out. The main intent of contemporary study is the analysis of the radiative MHD Micropolar nanofluid having micro-organisms. Furthermore, the article is made more fascinating by the usage of solutal and thermal boundary conditions with radiative heat flux in the unsteady Micropolar nanofluid fluid flow over the stretching sheet. The deportment of existing parameters is demonstrated graphically through an appropriate discussion. After that, suitable similarities have been used for transformation; the governing non-linear partial differential equations are composed in a non-linear system of ordinary differential equations (ODEs). The resulting system of non-linear ODEs has been solved numerically with a proficient and authenticated variational finite element method (FEM) along with the boundary conditions. The influences of various parameters are studied graphically. Furthermore, the graphical narration of Nusselt number and microorganism flux is accessible and the skin friction behavior and also the impact of different parameters of the flow is numerically inspected. After that, the numerical comparison of the existing results has been presented and discussed with graphs. In view of this study, transient flow with slip effects with the existence of mixed convection and chemical reaction on the sheet/disk can be observed.

Problem Description
Consider the 2-D unsteady boundary layer MHD bioconvective flow of an incompressible fluid and micropolar nanofluid viscous flow over an electrically conducted stretching/shrinking sheet has been considering for exploration. The solutal and slip impacts have also been integrated. We made a decision with the Cartesian system with the purpose that the measurement of the sheet would be taken with x-axis, y-axis, and the x-axis would be chosen beside the stretching sheet, the vertical axis would be perpendicular to the surface as given in Figure 1. U(x, t) = ax/(1 − λt) is the non-uniform velocity of sheet, since a is the stretching rate along x-axis but λt is positively constant as by condition λt < 1. The uniform magnetic field is supposed to be a function of distance from the origin, where B 0 is the potency of the magnetic field and is applied along the direction of the positive y-axis. The magnetic Reynolds numeral is very small as compared to the induced magnetic field. The free stream temperature is supposed to be T ∞ , the free mass concentration is as C ∞ and the concentration of microorganisms is n ∞ . The governing equations for flow can be expressed (Ullah et al. [8]).

∂u ∂x
where x and y are the co-ordinates with the x and y-axis, u is the velocity components along the x-axis and v with the y-axis and α, ν, σ, are the thermal diffusivity, kinetic viscosity, and electrical conductivity, respectively. And symbolizes that ρ f is nanofluid density, ρ p is nanoparticle density, ρ m is microorganism density and g is the gravity. The temperature is T, C is the solutal concentration; n is the concentration of microorganisms and T 0 , C 0 and n 0 represents the reference temperature, reference concentration of nanoparticles and reference concentration of the microorganisms, respectively. T w , C w , n w and T ∞ , C ∞ , n ∞ symbolizes the surface temperature, the concentration of nanoparticles, the concentration of the microorganisms, free stream temperature, compactness of nanoparticles and the concentration of microorganisms, respectively, and N represents the angular velocity. D M is the molecular diffusivity; D T is the thermal diffusivity, D B is the Brownian diffusivity, in Equation (3), j where κ f is the eddy viscosity and Rosseland eradicative heat flux is q r , which is delineated by q r = −( 4σ * 3κ * ∂T 4 ∂y )j wherever σ * is Stefan-Boltzmann numeral and κ * is the mean assimilation coefficients. And in Equation (6), n is the compactness of microorganisms, the chemotaxis invariable is b, W c is velocity cell swimming. The boundary conditions for the mathematical model are followed as (see [8]): wherever the injection/suction velocity is is the temperature of the sheet and C w (x, t), E w (x, t) are concentrations at surface of the below form [34]: where T 0 , C 0 and n 0 are the reference temperature, solutal concentration and bio-convection concentration respectively, such that 0 Generally, the stream function Ψ is defined as u = ∂Ψ ∂y and v = − ∂Ψ ∂x which satisfies the equation of continuity automatically. Following are the similarity transformations to solve the above Equations (see [34][35][36]): In view of above similarity transformations Equation (9), PDEs from (1)-(6) transform into the following system of nonlinear ordinary differential equations (ODEs): for the above problem, the transformed boundary conditions are: The primes show the differentiation with respect to η. The using parameters in Equations (10)-(14) are defined as: where unsteadiness parameter σ = λ/a, M is the magnetic field parameter, Pr is the Prandtl number, Nb is the Brownian motion parameter, Nt is the thermo-phoresis parameter, β is the Deborah number, Sc is the Schmidt numeral, also P e , λ, G r , R ex , K, E c , L b , R d , R b are bioconvection Peclet numeral, buoyancy parameter, Grashof numeral, Reynolds numeral, micro-polar constant, Eckert numeral, Lewis parameter, radiation parameter, Rayleigh number and Ω the microorganism concentration difference. The corporal measurements of the most significant are the skin friction coefficient, Sherwood number, Nusselt number and the compactness of the motile microorganisms. The explanation of these dimensionless corporal measurements is given as: when Equation (7) substituted in Equations (17)-(19) as a result final dimensionless form is obtained; , where the reduced skin friction is C f r, reduced nusselt number is Nur, density number of the motile microorganisms is Nn x and the reduced sherwood number is Shr.

Implementation of Method
The finite element method has been implemented to acquire the numerical solution of the system of non-linear boundary value problem that is given in (10)-(14) subjected with the boundary conditions (15)- (16). The finite element method is more proficient and more consistent than other numerical methods such as ADM, HPM, FDM and so on. FEM is very capable and has been applied to learning the different problems in fluid mechanics and computational fluid dynamics, solid mechanics, mass transfer and heat transfer in several other fields [35,37]. The general detail of the finite element method (FEM) is explained in Reddy [38], providing the general aspects of the variational finite element method and furthermore found that the finite element method (FEM) is employed in commercial software like MATLAB, ANSYS, ADINA and ABAQUS. Swapna et al. and Rana et al. [39][40][41] elucidate that the variational finite element method solves the boundary value problem very proficiently and correctly. To solve the non-linear system of Equations (10)-(14) by using the finite element method (FEM), along with boundary conditions (15)-(16), first we have to consider: The Equations (10)-(14) takes the form The corresponding boundary conditions reduce to following form:

Variational Formulations
The variational structure of the Equations (21)-(26) over a distinctive two-nodded domain (η c , η c+1 ) is given as where w 1 , w 2 , w 3 , w 4 , w 5 and w 6 are weight functions and are viewed as variations in f , h, g, θ, φ, ξ respectively.

Finite Element Formulations
The equation of finite element model that is attained from the Equations (29) to (34) by exchanging the FEM approximation as the following: with where the shape function ψ e i are the functions for the element (η c , η c + 1) that are taken as where [W mn ] and b m (m, n = 1,2,3,4,5,6) are the matrices and are given by: j=1h j ψ e j ,θ = ∑ 3 j=1θ j ψ e j andφ = ∑ 3 j=1φ j ψ e j are assumed to be known. After the grouping of element equations, a resulting system of non-linear equations is acquired; after this, it necessitates an iterative method to solve it for a proficient solution. The computation of f , g, θ, φ and ξ are then approved out for a higher level and continued it until the required 0.00005 is not obtained. From the outcomes that in Table 1 illustrates the convergences, as we calculated for the escalating number of elements, n = 100, 180, 280, 450, 600, 850. From the outputs, it is shown that as the number of elements ascending more, than there is no significant difference in the values of f , g, θ, φ and ξ can be seen as the number of elements increases beyond 850, so the results at n = 850 elements are reported.

Results and Discussion
The numerical computations were executed for the temperature, velocity, solutal and micropolar nanofluid volume fraction functions for a different evaluations of physical parameters, as Prandtl number Pr, magnetic parameter M, Unsteadiness σ, Nr the buoyancy ratio parameter, Brownian motion Nb, Rayleigh number Rb, thermophoresis Nt, Schmidt number Sc, Pe Peclet number, λ are the buoyancy parameter, G r the Grashof number, Reynolds number R ex , K the micropolar constant, E c the Eckert number, L b the bioconvection Lewis parameter, the heat generation coefficient, R b the Rayleigh numeral, Ω the microorganism concentration difference parameter, buoyancy parameter λ, f w is Suction/Injection. The detail of the present results and the appraisal of flow velocity is made with the exact solution which consents the validity of the finite element method. It can be seen clearly in Tables 2-4. To ensure the accuracy of presented numerical results, the results obtained by finite element method for skin friction co-efficient for the steady and unsteady flow have been compared with the numerical results that have been previously reported in existing studies in Table 2, and regarding our results, there is admirable conformity among our outcomes and previously available research articles that approve the cogency and the accuracy of the current results that are obtained by the finite element method (FEM). Table 2 describes the results of the heat transfer rate that is acquired by the finite element method which is compared with the results of earlier studies and with the accurate solution of references [34,37,42]. We found that our results are in absolute accord and the grid invariance test was executed to sustain accuracy for 5-decimal points.   Table 3 represents the results of the heat transfer rate acquired by FEM that are estimated with the results of earlier studies and the exact result of Ishak et al. [44] in the particular case (Nt = Nb = 0). It has been observed that the attained numerical results are in complete accord and a grid invariance test has been performed to preserve four decimal points of correctness.Of all the calculations according to this critique, ξ(η) have been elected as 10 −6 . For more consistency and justifications of the results, the consequences of − f (0) and g (0) that have been detailed by Hsiao [48] and Eldabe [47] are effectively replicated and shown in Table 4. Figure 2 is evidence for the effect of M on the velocity function with and without hydro-dynamic slip. The results show that in each case the velocity component decreases with increasing value of M. Actually, the magnetic parameter M generated Lorentz force reminiscent of a drag force; that is why there is deceleration of motion of the fluid. Figure 2 also shows that in both conditions the suction f w deceases the momentum boundary layer thickness that can be controlled with f w . To observe the consequences of the buoyancy parameter λ on the velocity profile with non-existence and the existence of hydrodynamic slip Figure 3 is presented. In this case, we have to perceive that the momentum boundary layer declines with the escalating values of the buoyancy parameter. As to see the effects of the buoyancy ratio parameters on the velocity profile with and without hydro-dynamic slip, as shown in Figure 4. The dimensionless velocity profile f (η) is decreasing while the opposite is the case for Nr in the absence of slip and the velocity slip thickness is decreasing with the existence of slip. In Figure 5 it is noted that, due to the disparity in the velocity profile with and without hydro-dynamic slip there is an increase in the bio-convection Rayleigh number Rb. In order to increase Rb, the buoyancy force is increased, due to which the velocity of the fluid is decelerated. A similar flow prototype was observed by reference [49] in their work.  The impact of σ on velocity profile for S f = 0 and S f = 0.5 is revealed in Figure 6. It is observed that escalating values of σ decrease the fluid velocity in the locality of the stretching sheet. Actually, this illustrates that fluid velocity contiguous to the sheet is less than the velocity of the stretching sheet while slither takes place. Escalating σ adequately means more fluid slides over the sheet and the flow slump is close to the sheet more rapidly. Furthermore, it shows that the fluid velocity close to the sheet will not remain the same as the sheet stretching velocity as the slide transpire, for that reason; there is indeed velocity. An analogous tendency was inspected by reference [36] in their work. With the increasing values of Nc and Rd, Figure 7 illustrates that both the temperature and gyrotactic microorganism concentration increases with an increase in the value of M. The main rationale of this exploit is that temperature and concentration gradients decline as existing transitory through moving fluids and as a result, ascend in thermal and concentration boundary layer thickness. Figure 8 shows the effect of the Nb on temperature θ(η) which reveals that an enrichment in Nb causes an increase in θ(η) at different values of Eckert number Ec. A greater thermal boundary layer width is created with high Nb values while higher concentration boundary layer thickness is achieved with lesser Nb values. Both the thermal profile and the associated boundary layer thickness are elevated for the increasing assessments of Nb. Figure 9 is shown to envisage the inconsistency in temperature θ(η) due to the thermophoresis parameter Nt. Also, it was proved that temperature was enhanced as the thermal boundary layer grows to be thicker for both the steady and unsteady cases.  Figure 10 demonstrates the influence of Nb on φ(η). The decline in value of Nb is related to an increase in the value of φ(η). Actually, Nb exacerbates the particles of the fluid regimes; for that reason, an increment is observed in the concentration profile. Figure 11 demonstrates the influence of the magnetic parameter M on φ(η) for both Sc = 1 and Sc = 2. An acclivity in M causes an increase in the value of φ(η). These curves show an increase in the values of M and enrichment in the concentration profile. Figure 12 shows the motile compactness profile ξ(η) for the escalating values of the magnetic numeral M with and without hydro-dynamic slip. These curves revealed that an increase in values of M gives an enrichment in the motile compactness profile φ(η). The density profile ξ(η) is a declining function of Lewis number Lb for P e = 2 and P e = 5 as exposed in Figure 13. The elevated values of the Lb mean a decrease in the diffusivity of microorganisms, so as a result, ξ(η) is reduced.  Figure 14 illustrates the impacts of Pe on the microorganism profile. It illustrates that the microorganism profile is enhanced while escalating the bio-convection Peclet number Pe. While Pe contains interrelation by means of maximum cell swimming movement, chemotaxis is constant and the opposite relation in the diffusivity of the microorganism. As a result, for higher values of Pe, the microorganisms profile will be enlarged. In Figure 15 it is indomitable that the motile concentration profile ξ(η) for unsteady parameter σ is enhanced. Furthermore, a decrease in associated boundary layer thicknesses with and without hydro-dynamic slip is also eminent. The discrepancy in motile concentration profile ξ(η) due to attachment of assorted parameters and the effect of the material parameter K on the motile concentration profile ξ(η), is demonstrated in Figure 16. An enhancement in the value K causes the motile concentration profile ξ(η) to decline.   Figure 17 that the skin friction coefficient declines with the ascending value of the slip parameter, magnetic and unsteadiness parameter and Figure 18 shows the influence of M, σ and S f on g (0). Figure 19 allows the investigation of the effect of M and σ on Nusselt number with S f = 0 and f w = 0.5. It can be observed that the mass transfer ratio enhances as M and σ increases. It is also resolute that in the case of nanofluid the mass transfer ratio is larger in disparity by viscous fluid. Lastly, Figure 20 divulged the discrepancy of local Sherwood number for escalating values of magnetic parameter M the bio-convection Lewis parameter Lb and bio-convection Peclet number Pe. Intriguingly, the mass transfer rate is seen as an increasing function of Lb and Pe. While escalating values of M decreases the mass transfer rate. Thus, it can exemplified that the mass transfer rate decreases as M enhance. All these interpretations are consistent with the consequences that are shown in Table 5. In Table 5 we analyze the variation of physical parameters M, λ, N r , Rb, σ, P r , L e on skin friction co-efficient velocity g (0), − f (0), Nusselt number −θ (0), Sherwood number −ξ (0). The following results are concluded from Table 5: (1) The Skin-friction coefficient and velocity are increasing while reducing the local Nusselt and density number through improvement in the Magnetic parameter.

Conclusions
An unsteady the flow of magneto-hydrodynamic (MHD) bioconvective micropolar nanofluid containing gyrotactic microorganism has been inspected numerically. The study was performed to analyze the solution of a 2-D laminar bioconvective boundary layer flow of water-based micropolar nanofluids over a stretching sheet. The mathematical sculpt was rehabilitated into the set of coupled ordinary differential equations by applying suitable transformations; the resulting system of ordinary differential equations was evaluated numerically with a proficient and validated finite element method (FEM). FEM is a numerical scheme for the estimation of a parameter-dependent problem that is used to find tht the accuracy of the approximations depreciates as the parameter tends to a limiting value. For the given problem, the present results in Table 3 shows the accuracy, up to four decimal places, of the fast convergence for all values of the parameter. Furthermore, for a particular case, the comparison of the current results with previous studies of the presented model has been discussed. A parametric study has been performed to investigate the influence of various parameters on the velocity, temperature, micro-rotation, gyrotactic microorganism volume fraction and the motile microorganism's profile. Moreover, local skin friction, heat transfer rate and the motile microorganism's thickness rate are scrutinized graphically. A few motivating elucidations from the present work are affirmed as,

•
The fluid velocity, temperature and gyrotactic microorganism concentration are found to have declined with enhancement in the unsteadiness parameter.

•
The Skin-friction coefficient decreases with the increasing value of slip parameters, magnetic and unsteadiness parameter but the effect is the opposite for escalating values of thermal buoyancy, suction parameter and solutal buoyancy.

•
The temperature declines, whereas the gyrotactic microorganism concentration decreases with enhancement in values of magnetic parameter M, Brownian motion parameter Nb, thermophoresis parameter Nt.

•
Due to an enhancement in the material parameter K, unsteady parameter as σ, bio-convection parameter λ and bioconvection Lewis parameter Lb, the motile concentration profile ξ(η) is declined.

•
Mass transfer rate increases as Brownian motion parameter Nb increases whereas decreases with enhancement in thermophoresis parameter Nt.

Conflicts of Interest:
The authors declare no conflicts of interest.

Microorganism concentration Rd
Radiation parameter Rb Rayleigh number T w Wall temperature of the fluid T ∞ Temperature of the fluid far away from the sheet U(x, t) Velocity of sheet D T