Entropy Generation in MHD Mixed Convection Non-Newtonian Second-Grade Nanoliquid Thin Film Flow through a Porous Medium with Chemical Reaction and Stratification

Chemical reaction in mixed convection magnetohydrodynamic second grade nanoliquid thin film flow through a porous medium containing nanoparticles and gyrotactic microorganisms is considered with entropy generation. The stratification phenomena, heat and mass transfer simultaneously take place within system. Microorganisms are utilized to stabilize the suspended nanoparticles through bioconvection. For the chemical reaction of species, the mass transfer increases. The governing equations of the problem are transformed to nonlinear differential equations through similarity variables, which are solved through a well known scheme called homotopy analysis method. The solution is expressed through graphs and illustrations which show the influences of all the parameters. The residual error graphs elucidate the authentication of the present work.


Introduction
Entropy is a thermodynamic quantity representing the unavailability of a system's thermal energy for conversion into mechanical work, often interpreted as the degree of disorder or randomness in the system. Entropy generation is an important part of the mechanical system. In the context of the present situation, irreversibility (entropy generation) is reliant to heat transmission, viscous dissipation and effects of magnetic field as well as diffusion. Irreversibility is studied with a sharp rate. The present study also discusses bioconvection and chemical reaction. Bio-convection is associated with a macroscopic flow of convecting fluid under the effect of density gradients due to hydrodynamic propulsion such as swimming of motile microorganisms. Mixed convection has several applications including heat dismissal into the atmosphere like lakes, rivers and seas. It is also utilized in the storage of thermal energy like solar ponds and heat transfer thermal sources such as the condensers of power plants. The applications of biosynthesized nanoparticles in a wide spectrum of potential areas are targeted drug delivery, cancer treatment gene therapy and DNA analysis, antibacterial agents, biosensors, enhancing reaction rates, separation science and magnetic resonance imaging (MRI). Entropy generation and bioconvection is a hot subject for scientists and researchers. Some published works and their few findings are given in Table 1. Table 1. Some recent developments related to the topic under investigation.

Author Names Author Works Some Outcomes
Khan et al. [1] Entropy generation Entropy generation increases with Reynolds number Ishaq et al. [2] Entropy generation Entropy generation decreases with Eyring-Powell parameter Hayat et al. [3] Entropy generation Entropy generation increases with Reynolds number Khan et al. [4] Bioconvection Bioconvection decreases with reduced heat transfer parameter Zuhra et al. [5] Bioconvection Gyrotactic microorganisms depreciates with magnetic field parameter Khan [6] Bioconvection Stratification increases with second-grade fluid parameter Palwasha et al. [7] Bioconvection Stratification increases with gravitational forces Raees et al. [8] Bioconvection Bioconvection depends on upper plate Zuhra et al. [9] Bioconvection Microorganisms decrease with increasing Brownian motion parameter Pedley [10] Instability Instability of the system is due to microorganisms Xu et al. [11] Mixed convection Nanoparticles favor the mixed convection Aziz et al. [12] Free convection Instability of the system increases due to microorganisms Bhatti et al. [13] Chemical reaction Mass transfer increases with chemical reaction Raees et al. [14] Bioconvection' Passively controlled nanofluid model provides better results Ramzan et al. [15] Chemical reaction Mass transfer increases with chemical reaction Anjalidevi [16] Chemical reaction Mass transfer increases with chemical reaction The study of applied magnetic field plays a significant role in heat and mass transfer, controlling momentum in boundary layer nanofluid flow through a shrinking/stretching surface. Moreover, magnetohydrodynamics is utilized in diverse engineering applications like cooling of nuclear reactor, electromagnetic casting and plasma confinement. In microfluidics complexities, magnetohydrodynamics is helpful to create non-pulsating and continuous movement in the structure of complex micro-channels. In biomedical engineering, magnetohydrodynamics is applicable for magnetic drug targeting in cancer diseases. On behalf of these applications, Khan et al. [17] explored that due to the strong applied magnetic field, flow becomes slow and temperature becomes high. They also found the residual errors for the authentication of the solution. Palwasha et al. [18] discussed the intensified magnetic field which causes the Hall current effect in two dimensions case. They found that magnetic field effect is a highly dominant factor to the flow. Khan et al. [19] presented the study to show the effects of applied magnetic field on graphene nanoparticles in a thin film Eyring Powell nanofluid involving velocity slip condition. Their results reveal that the thermal conductivity increases with increasing the concentration of graphene nanoparticles. Zuhra et al. [20] also provided the study on graphene nanoparticles in two dimensional unsteady flow and heat transfer in thin liquid film non-Newtonian (Casson and Williamson) nanofluids suspension past a stretching sheet in the presence of transverse magnetic field and non-uniform heat source/sink. Their study shows that the temperature increases with temperature dependent heat source sink parameter. Khan et al. [21] analyzed the magnetohydrodynamic study in cylindrical coordinates for nanoliquid thin film using spraying phenomenon in thermal system. Apart from various effects of different parameters on different profiles, they also discussed the pressure role in the system. Focusing on the MHD flow, Khan et al. [22] presented Brownian motion and thermophoresis effects on MHD mixed convective thin film second-grade nanofluid flow with Hall effect and heat transfer past a stretching sheet, manifesting that the flow becomes three dimensional due to Hall effect.
The study of motion through a porous space is also very important in engineering and industrial processes. A porous media is a material having pores, i.e., oxides and metals, where these pores are covered with a fluid. Such type of motion has two main benefits, as its area of dissipation is higher as compared to conventional fins that improve the heat convection. Further, the regular flow of a fluid along the whole individual beads blends the fluid effectively. The nanofluid flow through a porous media reveals enhanced thermal properties, i.e., convective heat transfer coefficients as compared with a base material and higher thermal conductivity. Consequently, a porous media having nanoparticles plays a vital role in the enhancement of heat transfer and this mechanism is very important from an application and theoretical point of view. Such types of flow can be seen in filtration, oil flow, heat exchangers, drying process and geothermal systems. Moreover, nanofluid flow through a porous media has many applications in thermal engineering which are favorable to improve the thermal properties of lubricants and oil, electronic cooling equipment and vehicle cooling. Applying Navier-Stokes equations to the physical problems, Khan et al. [23] examined the important role of flow and transfer in a second-grade fluid through a porous medium past a stretching sheet in which the porous medium carries a leading role. Their results show a close agreement with the published work. It is very important to explain how nanofluid flow is possible in a porous medium. Thus, without special precautions, nanoparticles will be simply absorbed by the porous matrix. Basically, the porous matrix will work as a filter for nanoparticles. This situation has been described, explained and modeled in the paper by Wu et al. [24]. The physical situation described in this paper reveals that the work on porous media filled by nanofluids is not just a mathematical challenge, but is based on a deep physical understanding of nanofluid flows. This demonstrates that we are simulating here a real physics problem of mixed convection flow and heat transfer in a porous medium filled by a nanofluid. Kuznetsove [25] has compiled the studies on the developments of bioconvection in porous media caused by either gyrotactic or oxytactic microorganisms, with the limitation of average pore size is larger than the size of a microorganism, which allows the smooth movement of the microorganisms. Tham et al. [26] numerically studied the steady mixed convection boundary layer flow on a horizontal circular cylinder with a constant surface temperature embedded in a porous medium saturated by a nanofluid containing both nanoparticles and gyrotactic microorganisms in a stream flowing vertically upwards for both cases of a heated and cooled cylinder. The objective of that paper is to understand the behavior of nanofluid suspensions containing nanoparticles, gyrotactic microorganisms and entropy generation at the fundamental level.
Motivated from the above discussion, the present study discusses the entropy generation and mixed convection in gravity-driven MHD second grade nanofluid flow containing both nanoparticles and gyrotactic microorganisms through a porous medium with chemical reaction along a convectively heated vertical surface. Employing suitable transformations, the basic governing equations of the problem are converted to dimensionless form which have been solved using a powerful analytic tool called Homotopy Analysis Method (HAM) [27]. The influences of all the pertinents parameters on velocity, temperature, nanoparticle concentration, density motile microorganisms and entropy generation profiles have been shown graphically and illustrated.

Basic Equations
A motion of a magnetohydrodynamic two dimensional, time independent, laminar and an incompressible second grade nanofluid falling downwards through a porous medium along a vertical solid surface due to gravity is considered. The uniform incoming flow on the right side of the plate has a constant temperature T ∞ at x = 0 and the nanofluid on the left side of the plate has another constant temperature T f . To avoid the bioconvection instability, it is assumed that the nanofluid is dilute. In addition, the assumption is taken for the stability of the nanoparticles suspended in the base fluid so that the nanoparticles do not agglomerate in the fluid. The passive control of nanoparticles volume fraction at the boundary (on the solid wall) is taken. It is assumed that the microorganisms have constant distributions on the wall. It is most important to note that the base fluid is water so that the microorganisms can survive. The assumption is also maintained that the existence of nanoparticles has few effects on the motion of the microorganisms. A transverse magnetic field of strength B 0 is applied in the positive y-direction, normal to the plate surface. There is no applied voltage and the magnetic Reynolds number is small, so the induced magnetic field and Hall effects are negligible. The physical model and coordinate system are shown in Figure 1. Applying the above-mentioned assumptions, the equations for the nano-bioconvection model are written in the following form as in [4,[13][14][15] ∂u ∂x Here dp where U is the free stream velocity, a is constant, u and v are the velocity components in the x-and y-directions. α 1 is the normal stress moduli, ∂y is the average swimming velocity vector of the oxytactic microorganisms in which b is the chemotaxis constant and W c is the maximum cell swimming speed. The subscripts p, f and f ∞ denote respectively the solid particles, the nanofluid and the base fluid at far field. ∆ρ = ρ cell − ρ f ∞ is the density difference between a cell and base fluid density at far field, µ f is the dynamic viscosity, γ av is the average volume of microorganisms, σ is the electrical conductivity and ρ f is the density of the nanoliquid. ν f = µ f ρ f is the kinematic viscosity, β is the coefficient of volumetric volume expansion of a second grade nanofluid, g is the acceleration due to gravity, C is the nanopartical volume fraction, k = k 0 x is the permeability of porous medium and k F = x is the Forchheimer resistance factor where k 0 is the initial permeability. N is the number density of motile microorganisms, C ∞ is the ambient nanofluid volume fraction, λ = k 1 ρ f is the thermal diffusivity of the nanofluid in which k 1 is the thermal conductivity, k r is the rate (constant) of chemical reaction, τ = (ρc) p (ρc) f is the ratio of nanoparticle heat capacity and the base fluid heat capacity, D B is the Brownian diffusion coefficient, D n is the diffusivity of microorganisms, T ∞ is the ambient temperature, D T is the thermophoretic diffusion coefficient and T is the temperature inside the boundary layer.
The boundary conditions for the nano-bioconvection flow are where h f (x) denotes the heat transfer coefficient due to T f , N w is the wall concentration of microorganisms and N ∞ is the ambient concentration of microorganisms. Note that, in order to satisfy the boundary conditions at infinity, we have to set N ∞ = 0. Introducing the transformations for nondimensional variables f, θ, φ, Ω and similarity variable ζ as where ψ is the stream function such that u = ∂ψ ∂y and v = − ∂ψ ∂x , x and y are the Cartesian coordinates along surface and normal to it. Equation (9) automatically satisfies mass conservation Equation (1). With the help of Equation (9), the Equations (2)-(5), (7) and (8) yield the following six Equations (10)- (15).

Parameter Names Symbols/Notations Defined Values
Dimensionless second-grade nanofluid parameter γ 1 Buoyancy ratio parameter Nr

Entropy Generation
The subject of the present analysis has a close relation with engineering and industrial processes so it is interesting to investigate the rate of entropy generation. The volumetric entropy generation of the fluid can be written as where R denotes the ideal gas constant and D represents the diffusivity. In Equation (16), on the right hand side, the first term represents the irreversibility due to heat transfer, the second term is entropy generation due to viscous dissipation, the third term is second grade fluid friction irreversibility, the fourth term is the local entropy generation due to the effect of magnetic field and the fifth to eighth terms are the irreversibilities due to diffusion effect. The characteristic entropy generation rate is where L is the characteristic length. The entropy generation E G is the ratio between the entropy generation rate E gen and the characteristic entropy generation rate E 0 i.e., Applying Equation (9), the non-dimensional form of Equation (18) for entropy generation number where is the diffusive constant parameter due to gyrotactic microorganisms concentration, are respectively the temperature, nanoparticles concentration and microorganisms concentration differences. Due to entropy generation rate, the second-grade nanofluid parameter and magnetic field parameter, here become respectively

Analytical Solution of the Problem by Homotopy Analysis Method
The initial approximations and the auxiliary linear operators for velocity, temperature, nanoparticles concentration and motile microorganism concentration are presented as The following properties are related to the linear operators where C i (i = 1-9) are the arbitrary constants.

Zeroth-Order Deformation Problems
Introducing the nonlinear operator ℵ as ℵ where p is an embedding parameter such that p ∈ [0, 1]. The zeroth-order deformation problems are expressed as where h represents the auxiliary non-zero parameter.
Equations (27)-(30) have the boundary conditions in Equations (31)-(34) respectively Assuming p = 0 and p = 1, the following results hold Similarly Using Taylor series expansion and Equations (27)-(30), gives the following results The convergence of the series is closely related to h. The value of h is taken in such a way that the series in Equations (39)-(42) converge at p = 1, then Equations (39)-(42) result in

Velocity Profile
This section depicts the influence of numerous dimensionless parameters on the velocity f (ζ). Figure 6 shows that the velocity profile decreases significantly with an enhancement in the values of second grade nanofluid parameter γ 1 . Velocity is minimum for non-Newtonian nanofluids in general. Figure 7 demonstrates that the non-dimensional velocity profile f (ζ) increases for the greater values of reduced heat transfer parameter γ 2 due to slip condition. Figure 8 reveals that the non-dimensional velocity profile f (ζ) depreciates for the high values of porosity parameter γ 3 . The holes of the medium creates resistance to the flow thereby decreasing the motion. Figure 9 is prepared for the inertial parameter γ 4 . An increment in this parameter decreases the velocity f (ζ). In Figure 10, it is seen that the non-dimensional velocity f (ζ) amplifies with increasing magnitude of chemical reaction parameter γ 5 which shows that velocity enhances with the strengthening of chemical reaction. Figure 11 shows that the non-dimensional velocity f (ζ) increases with the increasing magnitude of buoyancy parameter Gr. Gravitational force favors the second grade nanofluid flow along the vertical solid surface. The buoyancy ratio parameter Nr is showing its effect in Figure 12. An enhancement in the parameter Nr makes a sharp advancement in the nanofluid flow so velocity f (ζ) becomes high. In Figure 13, the velocity f (ζ) decreases for increasing quantities of bioconvection Rayleigh number Rb which characterizes the effect of up-swimming of microorganisms. It should be noted that there is no effect of Rb on the boundary layer thickness. For the case of pure natural convection, Rb = 0 which implies that microorganisms have the same density as water and, therefore, they do not contribute to density stratification (or which corresponds to the situation with zero concentration of microorganisms). Since the microorganisms swim on average in the upward direction, their presence produces a destabilizing effect and reduces the critical value of Rb. By increasing the average concentration of corresponding microorganisms in the suspension, Rb can be increased. Figure 14 depicts the role of Lewis number Le. Velocity f (ζ) is found to be maximized with greater values of Le. Nanoparticles concentrations are consistently depressed with greater lewis number; nevertheless the over-riding influence is imposed by the lewis number on the diffusion of nanoparticles. The role of Schmidt number Sc is projected in Figure 15 which shows that velocity f (ζ) amplifies with the effective diffusion of microorganisms and the acceleration of nanoparticles. Figure 16 demonstrates that the velocity f (ζ) decelerates with the rising values of magnetic field parameter M. It is due to the fact that the magnetic field exerts a force, known as Lorentz force which suppresses the velocity. Figure 17 illustrates that the velocity f (ζ) enhances by the high values of the Prandtl number Pr. It is evident that high quantities of Pr make enhancement in the mass transfer rate so f (ζ) elevates due to gravity.

Temperature Profile
Temperature profile exhibits the interesting behavior for various parameters. Figure 18 shows the influence of second grade nanofluid parameter γ 1 . It demonstrates that the temperature θ(ζ) increases with the non-Newtonian effect of the nanofluid. The thermal slip effect is led by reduced heat transfer parameter γ 2 and is exhibited in Figure 19. The graph projects that with the elevation of reduced heat transfer parameter γ 2 , temperature θ(ζ) diminishes consequently the decrement comes into play in the thermal boundary layer thickness. Figure 20 reveals that as the porosity parameter γ 3 increases, the surface area of the porous media increases which drops down the temperature. Figure 21 depicts that as the values of inertial parameter γ 4 increases, the temperature diminishes in the boundary region. Figure 22 clearly shows that temperature increases with increasing chemical reaction parameter γ 5 . It is seen that temperature θ(ζ) increases with the sharpness of chemical reaction. Indeed, some chemical reactions like fissions and fusions make the temperature high. The thermophoresis parameter Nt gives rise in temperature with larger values of Nt as demonstrated by Figure 23. The thermophoresis parameter Nt is directly proportional to the heat transfer associated with the hot nanofluid. Thus an increment in the thermophoresis parameter Nt causes reduction in the heat transfer at the surface, consequently the temperature at the surface of the plate increases. Physically, it is due to the fact that with an enhancement of the thermophoresis parameter Nt, the random motion of the particles increases, resulting in an increase in the temperature profile. Due to the Brownian motion, the particles are much closer and so enhance the heat flow among the particles. Figure 24 depicts that the nondimensional temperature θ(ζ) decreases for the increasing values of magnetic field parameter M. It is due to the gravity because temperature rises at height and it falls down in tending towards ground. The influence of Prandtl number Pr has been projected in Figure 25 showing that the temperature θ(ζ) diminishes with the increasing values of Pr. The reason is that the increasing values of Prandtl number increase the viscosity of the fluid which absorbs the heat as a result, cooling phenomenon is generated.

Nanoparticle Concentration Profile
Nanoparticles having one or more dimensions of the order 100 nm or less have attracted great attention due to their unusual and fascinating properties and applications. Figure 26 displays that the concentration decreases with the non-Newtonian second grade nanofluid parameter γ 1 . The non-Newtonian second grade nanofluid is saturated so there is less opportunity for the diffusion of the nanoparticles, therefore, the concentration decreases. Figure 27 shows the effect of reduced heat parameter γ 2 on concentration φ(ζ). It is a common phenomenon that concentration depreciates with decreasing temperature so in the present case, the concentration φ(ζ) falls down with decreasing the reduced heat parameter γ 2 . Figure 28 focuses on the role of porosity parameter γ 3 . It provides that the concentration is enhanced with the rising values of porosity parameter γ 3 . The performance of inertial parameter γ 4 is shown in Figure 29 against the non-dimensional concentration profile φ(ζ). Concentration makes development with inertia since saturation is easily received. The highlights of the impact of chemical reaction parameter γ 5 is revealed in Figure 30. It is reported that for constructive (or generation) chemical reaction, the mass transfer shows an increasing behavior i.e., it is seen that the concentration profile increases significantly for the larger values of the chemical reaction parameter γ 5 . Figure 31 projects the role of Brownian motion parameter Nb. Due to Brownian motion, the random motion of the particles increases, which favors the non-dimensional nanoparticle concentration profile φ(ζ). Figure 32 illustrates that by enhancing the thermophoresis parameter Nt, the concentration profile φ(ζ) depreciates. The nanopartical concentration boundary layer thickness is an increasing function of the thermophoresis parameter Nt. This indicates that an increment in the thermophoresis parameter Nt induces resistance to the diffusion of the solute and this assists in the reduction of the concentration gradient at the surface. It is interesting to note that the effect of the thermophoresis parameter is more pronounced for the nanoparticle volume fraction than on the fluid velocity and temperature profiles. Figure 33 is specified for the effect of Lewis number Le on the concentration field φ(ζ) which reveals that the concentration becomes rich for the different values of the Lewis number Le. The fact is that the performance of the nanoparticles movement is developed, which enhances the nanoparticles concentration. An increment in the Lewis number Le causes the reduction in the mass diffusivity of the nanofluid which in turn reduces the solute concentration. Figure 34 shows that the nondimensional concentration field φ(ζ) increases as the magnetic field parameter M increases. It is due to the existence of Lorentz force which shows resistivity to the over all motion, consequently, nanoparticle concentration increases. The magnetic field parameter increases the friction between nanofluid layers so the concentration boundary layer thickness gets enhanced. The influence of Prandtl number Pr on the nanoparticle concentration profile φ(ζ) is depicted in Figure 35. Prandtl number is the fluid viscosity parameter so by increasing Prandtl number, the internal viscous forces of the fluid is strong and the concentration goes to peak which causes no special capability for further nanoparticles concentration. Therefore, concentration φ(ζ) becomes weak with larger values of Pr.

Gyrotactic Microorganism Concentration
There have been tremendous developments in the field of microorganism produced nanoparticles and their applications. The characteristics of the density motile microorganisms with respect to different parameters are discussed in the following. Figure 36 elucidates that the microorganism concentration profile Ω(ζ) grows up with increasing magnitude of second grade nanofluid parameter γ 1 . The reason is that increasing the quantity γ 1 increases the microorganism propulsion by providing the best platform for the growth of microorganism and hence explains the reason why the microorganism concentration boundary layer of second grade nanofluid enhances. Figure 37 shows the influence of reduced heat parameter γ 2 on the microorganism concentration Ω(ζ) profile. Intensity of temperature is not good for the survival of microorganism concentration. Therefore, in the present case, the microorganism concentration Ω(ζ) increases with moderate temperature. Figure 38 projects that the the microorganism concentration Ω(ζ) decreases with amplifying the porosity parameter γ 3 . In porous medium there is a resistance of holes to the flow of microorganisms so their progress is not significant in that case. Figure 39 illustrates that the non-dimensional microorganism concentration Ω(ζ) profile depreciates with the rising values of inertial parameter γ 4 . This is due to the fact that the inertial parameter affects the motion of the microorganisms since the environment is not clear to them. Figure 40 reveals that microorganism concentration field Ω(ζ) amplifies with the chemical reaction parameter γ 5 . Figure 41 shows that the microorganism concentration Ω(ζ) decreases against the Brownian motion parameter Nb because the propulsion of microorganisms is weak. Figure 42 demonstrates that by increasing thermophoresis parameter Nt, the microorganism concentration profile Ω(ζ) increases. This is due to the close relation of concentration with moderate temperature and it is expected that a higher thermophoresis parameter allows a deeper penetration of the microorganism concentration. This achieves enhancement in species in the nanofluid and elevates the motile microorganism concentration profile Ω(ζ). This indicates that a stronger thermophoresis effect induces a resistance in the mass diffusion. In Figure 43, the microorganism concentration field Ω(ζ) shows an upshot for the greater values of Lewis number Le. This happens because as the Lewis number Le increases, the viscous diffusion rate increases which in turn reduces the velocity at the surface and so elevates the density of the microorganisms. The effect of Schmidt number Sc on the microorganism concentration profile Ω(ζ) is projected in Figure 44. When the magnitude of Schmidt number Sc increases, the microorganism concentration profile Ω(ζ) decreases. This is attributed to the low propulsion of microorganisms. Figure 45 shows the effect of bioconvection Peclet number Pe on the microorganism concentration profile Ω(ζ). It indicates that the microorganism concentration profile Ω(ζ) grows large for different values of bioconvection Peclet number Pe. This is due to the fact that microorganisms diffuse more and more with larger values of bioconvection Peclet number Pe. The bioconvection Peclet number Pe helps to increase the speed of the microorganisms in respect of the nanofluid and so the density of the microorganisms is reduced near the surface. The behavior is highly influenced with an increase in the Lewis number Le which again helps in the reduction of the mass diffusivity. Figure 46 shows that by increasing magnetic field parameter M, the nondimensional motile density function profile Ω(ζ) decreases because the magnetic field opposes the flow due to Lorentz force.

Entropy Generation Analysis
It is necessary to study the irreversibility of the present system in terms of second-grade nanofluid, nanoparticles and microorganisms. Figure 47 shows that the irreversibility N G (ζ) of the system increases for high values of Reynolds number Re. Similarly, the values of Br make enhancement in entropy generation rate N G (ζ) which has the same increasing behavior in Figure 48. Different values of magnetic field parameter M in Figure 49 presents an increasing effect on the entropy generation rate N G (ζ). In Figure 50, the nondimensional diffusive constant parameter γ 6 assumes various high values (ranging from 0.50 to 1.20) with deep effects observations. γ 6 is related to the diffusivity. It shows that entropy generation rate N G (ζ) increases with increasing γ 6 where at the same time the temperature difference parameter θ w reveals a unique activity, including the depreciation of entropy generation rate N G (ζ) in Figure 51. The formation of curves of the entropy generation rate N G (ζ) in Figure 52 is due to the various values of the concentration difference parameter φ w , a reactive species making increment in N G (ζ). In Figure 53, it is shown that microorganisms concentration difference parameter Ω w favors the entropy generation rate.

Residual Errors
Residual errors are shown in the graphs to show the authentication of the problem. By fixing all the others parameters, different values are assigned to convergence control parameter h to show the effects in

Conclusions
This article discusses the analytical solution of entropy generation and mixed convection in gravity-driven non-Newtonian second grade nanofluid flow containing both nanoparticles and gyrotactic microorganisms through a porous medium with a chemical reaction along a convectively heated vertical solid surface. The thin film considered in this work contains the mixture of nanoparticles and motile gyrotactic microorganisms [14]. The solution of the problem has been obtained by using analytical technique called Homotopy Analysis Method (HAM) for the velocity, temperature, concentration and microorganism concentration fields, which has been displayed in the diagrams. From (i) The velocity f (ζ) depreciates for the porosity parameter γ 3 , inertial parameter γ 4 , bioconvection Rayleigh number Rb and magnetic field parameter M while it elevates for the second grade nanofluid parameter γ 1 , reduced heat transfer parameter γ 2 , chemical reaction parameter γ 5 , buoyancy parameter Gr, buoyancy ratio parameter Nr, Lewis number Le, Schmidt number Sc and Prandtl number Pr.
The temperature θ(ζ) diminishes for the reduced heat transfer parameter γ 2 , porosity parameter γ 3 , inertial parameter γ 4 , magnetic field parameter M and Prandtl number Pr while it elevates for the second grade nanofluid parameter γ 1 , chemical reaction parameter γ 5 and thermophoresis parameter Nt. (iii) The nanoparticles concentration φ(ζ) diminishes for second grade nanofluid parameter γ 1 , reduced heat transfer parameter γ 2 , thermophoresis parameter Nt and Prandtl number Pr while it elevates for the porosity parameter γ 3 , inertial parameter γ 4 , chemical reaction parameter γ 5 , Brownian motion parameter Nb, Lewis number Le and magnetic field parameter M. (iv) The microorganism concentration Ω(ζ) diminishes for the porosity parameter γ 3 , inertial parameter γ 4 , Brownian motion parameter Nb, Schmidt number Sc and magnetic field parameter M while it elevates for the second grade nanofluid parameter γ 1 , reduced heat parameter γ 2 , chemical reaction parameter γ 5 , Lewis number Le, thermophoresis parameter Nt and bioconvection Peclet number Pe.
(v) Entropy generation rate N G (ζ) diminishes with temperature difference parameter θ w while it elevates for Reynolds number Re, Brinkman number Br, magnetic field parameter M, diffusive constant parameter γ 6 , nanoparticles concentration difference parameter φ w and microorganism concentration difference parameter Ω w . (vi) Residual errors graphs are self explanatory for the efficiency of HAM solution. Acknowledgments: The authors are thankful to the Higher Education Commission (HEC) Pakistan for providing the technical and financial support. The authors also express their deep and hearty thanks to the honorable reviewers for useful comments and information on the literature.

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

Abbreviations
The following abbreviations are used in this manuscript: Chemotaxis