Magneto-Bioconvection Flow of Williamson Nanofluid over an Inclined Plate with Gyrotactic Microorganisms and Entropy Generation

Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1 Department of Mathematics, Adeleke University, Ede 240003, Osun State, Nigeria; tundeayusuf04@gmail.com 2 Department of Information Technology, Fanshawe College, London, ON N5Y 5R6, Canada; mabood1971@yahoo.com 3 Department of Studies and Research in Mathematics, Davangere University, Davangere 577126, Karnataka, India; dr.bcprasanna@gmail.com 4 Department of Mechanical Engineering, University of West Attica, 12244 Athens, Greece * Correspondence: sarris@uniwa.gr


Introduction
The industrial cooling system relies on a standard base fluid, which is not sufficient to meet the requirements of industrial processes. To overcome this issue, a modern phase of nanotechnology is proposed to achieve a higher cooling rate in the performance of industrial processes. Nanofluids can be exploited in cars as cooling and shock absorbers and in improving refrigeration/air-conditioners. Inspired by these applications, several researchers examined the flow of nanoliquids past diverse surfaces [1][2][3][4][5]. The Williamson liquid is an important class of pseudo-plastic liquid model. The study of the boundary layer stream of pseudoplastic liquid is of great interest due to its widespread use in industry such as polymer sheet extraction, emulsion sheets and melting of high-density polymers. Inspired by these, Prasannakumara et al. [6] examined the consequence of chemical reaction in a Williamson nanoliquid stream through an extending sheet. Mabood et al. [7] deliberated on the radiative stream of Williamson liquid past an incessantly moving intense surface. Rasool et al. [8] exemplified the consequences of chemical reactions on a Williamson nanoliquid stream through a widening flat surface with entropy production.
Haq et al. [9] deliberated on the behavior of gyrotactic microorganisms on Williamson nanoliquid stream through a stretchy cylinder.
In the last few decades, a number of investigations into energy production in several fluid stream conditions have been carried out by researchers focusing on the generation of entropy. Research studies were motivated by recognizing the significance of entropy production in liquid flow and its importance in some industrial uses such as air conditioners, refrigerators, heat pumps, fire engines and much more. Entropy management efforts were first discussed by Bejan [10]. Later, Yusuf et al. [11] deliberated the entropy production in the magnetohydrodynamic (MHD) Williamson nanoliquid stream past a stretchy sheet with chemical reaction. Azam et al. [12] scrutinized the production of entropy on a Williamson nanoliquid stream. Bhatti et al. [13] explored the entropy production on a Williamson nanoliquid stream. For a better biological process, a better liquid mixture is needed, which is why gyrotactic microorganisms are essential. The main reason for bioconvection refers to the gradient size of gyrotactic microorganisms. This phenomenon has great potential for use. For example, it has been used for bio-microsystems of enzyme biosensors and biomedicine treatment of cancer. Waqas et al. [14] scrutinized the time-dependent MHD stream of Williamson liquid on taking account of motile gyrotactic microorganism with suspended nanoparticles. Khan et al. [15] examined the bioconvective stream of Williamson liquid with the suspension of nanoparticles through stretching sheet with gyrotactic microorganisms. Hayat et al. [16] inspected the influence of radiation effect on an MHD stream of non-Newtonian nanoliquid through a stretching cylinder with gyrotactic microorganisms. Jayadevamurthy et al. [17] explored the performance of a moving rotating disk on stream of hybrid nanoliquid on taking account of gyrotactic microorganisms.
The state of radiation contributes significantly to the formation of high temperatures. Non-linear heat transfer conditions in a stretched surface play a major role in fluid conversion and have many applications in chemical engineering and metallurgy. Over the past decades, several researchers have scrutinized the influence of thermal radiation effect on different fluid flow patterns with the existence of porous medium. Hosseinzadeh et al. [18] scrutinized the influence of chemical reaction and radiation effect on the non-Newtonian liquid stream through a plate on taking account of permeable medium. Dogonchi et al. [19] examined the influence of porous medium and thermal radiation effect on copper-water nanoliquid stream though a cylinder. Yusuf et al. [20] inspected the consequences of the radiation effect on a Darcy-Forchhiemer stream of hybrid nanoliquid through a stretchy sheet with a porous medium. Shehzad et al. [21] examined the impact of radiation effect on a micropolar nanoliquid stream through a stretching sheet on taking account of the porous medium. The diffusion of species on taking account of chemical reaction in the interface layer stream has abundant applications in atmospheric flows, fibrous insulation, water and air pollutions and several chemical engineering issues. Inspired by these applications, several researchers scrutinized the consequences of chemical reactions on different fluid streams on taking account of several influencing factors. Khan et al. [22][23][24] deliberated the influence of chemical reaction on non-Newtonian liquid streams through stretchy surfaces with production of entropy. Alsaadi et al. [25] inspected the impact of chemical reaction on a non-Newtonian nanoliquid stream through a stretching sheet with entropy production. Kotresh et al. [26] examined the impact chemical reaction on a nanofluid stream through a disk.
In view of the above literature survey, it is obvious to the best of the writer's knowledge, that no efforts have so far been initiated with respect to studying the encouragement of chemical reaction on the Williamson nanoliquid along an aligned semi-infinite plate embedded in a porous medium, which constitutes the novelty of the current study. In this investigation, we scrutinized entropy production on bioconvective flow of Williamson nanofluid along an aligned semi-infinite plate embedded in a porous medium under the influence of thermal radiation and gyrotactic microorganisms.

Problem Formulation
Mixed convective flow of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium under the influence of thermal radiation and gyrotactic microorganisms is examined. The flow configuration is depicted in Figure 1. The considered fluid is conducted in an imposed applied magnetic field of strength B 0 parallel to the y-axis. This semi-infinite plate is inclined about the vertical direction making an angle Ω. The stretchy surface velocity is u w = ax, where a is a positive constant. Here, we assume a two-dimensional steady flow and the fluid is assumed to be an incompressible. Also, the magnetic Reynolds number and the induced electric field are assumed to be small and negligible.
In this investigation, we scrutinized entropy production on bioconvective flow of Williamson nanofluid along an aligned semi-infinite plate embedded in a porous medium under the influence of thermal radiation and gyrotactic microorganisms.

Problem Formulation
Mixed convective flow of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium under the influence of thermal radiation and gyrotactic microorganisms is examined. The flow configuration is depicted in Figure 1. The considered fluid is conducted in an imposed applied magnetic field of strength 0 B parallel to the y-axis. This semi-infinite plate is inclined about the vertical direction making an angle Ω . The stretchy surface velocity is w u ax = , where a is a positive constant. Here, we assume a two-dimensional steady flow and the fluid is assumed to be an incompressible. Also, the magnetic Reynolds number and the induced electric field are assumed to be small and negligible. It is adopted that, at the boundary, the temperature, concentration and microorganism concentrations are characterized by a convective heating process and that nanoparticles have no influence on the movement of microorganisms. The influence of thermal radiation along with chemical reaction is accounted for. Williamson fluid is a pseudoplastic type of non-Newtonian fluid for which the stress component of the velocity and the extra shear stress xy S is expressed as (Nadeem et al. [27]) Here, 0 μ is the viscosity, Δ 0 > is the material constant of the fluid and γ is defined We assumed that the fluid properties are invariant. The magnetic Reynolds number, electric field due to polarization of charges, and Hall effects are assumed to be negligible. Under the above assumptions and boundary It is adopted that, at the boundary, the temperature, concentration and microorganism concentrations are characterized by a convective heating process and that nanoparticles have no influence on the movement of microorganisms. The influence of thermal radiation along with chemical reaction is accounted for. Williamson fluid is a pseudoplastic type of non-Newtonian fluid for which the stress component of the velocity and the extra shear stress S xy is expressed as (Nadeem et al. [27]) S xy = µ 0 1 − ∆ . γ ∂u ∂y + ∂v ∂x . Here, µ 0 is the viscosity, ∆ > 0 is the material constant of the fluid and . γ is defined as . We assumed that the fluid properties are invariant. The magnetic Reynolds number, electric field due to polarization of charges, and Hall effects are assumed to be negligible. Under the above assumptions and boundary layer approximations, the governing model can be expressed by exploring the Buongiorno nanofluid model given as follows: ∂u ∂x Subject to the boundary conditions For optically thick media, Rosseland approximation is expressed as where k s is the Rosseland mean spectral absorption coefficient and the Stefan-Boltzmann radiation is denoted as σ * . Therefore, Taylor series expansion about T ∞ Neglecting the higher order term of (T − T ∞ ) to get Combining Equations (8) and (10) yield Also in view of Equation (11), Equation (3) therefore becomes Now the following similarity variables are introduced Equation (1) is satisfied trivially and the dimensionless model are (1 + Rd)θ + Pr f θ + EcPr And the boundary conditions where prime represents the differentiation w.r.t (η) only, and β = ∆x

Entropy Generation
The effective expression in the entropy generation E G is characterized by entropy generation via heat transfer with thermal radiation, Joule dissipation, porous medium, fluid friction and diffusion effect. This is expressed as: The non-dimensional entropy generation number is given by The Bejan number is defined as the fraction part involving the entropy generation via heat transfer and diffusion and total entropy and is given as The irreversibility ratio is also given as 2 , Interesting quantities portraying the skin friction C f and the Nusselt number Nu, Sherwood number Sh and the nanoparticle Sherwood number N n at the vertical walls are expressed below where p w ,q w , j w and j n are defined as: The dimensionless forms are

Numerical Solution
In the present investigation, an effective numerical scheme, Runge-Kutta-4 with a shooting method, has been explored to examine the stream model Equations (14)- (17) in conjunction with boundary constraints (18) for varying values of governing parameters. The choice of step size and convergence criteria is carefully made to be 0.001 and 10 −6 respectively. In the absence of Williamson and mixed convection parameters for the special case of the model, the current results are in accordance with those conveyed by [28], as depicted in Table 1, and this serves as a benchmark for the accuracy of our numerical procedure.

Results and Discussion
The main aspect of this section is to scrutinize the entropy production in a mixed convective stream of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium, taking account of gyrotactic microorganisms and thermal radiation. The described flow pattern is expressed in terms of PDEs. The framed PDEs are reduced to ODEs by opting suitable similarity variables. To understand this, the ODE's with the boundary constraints, they are solved by applying an efficient numerical scheme Runge-Kutta-4 with shooting technique. In this section, the behavior of the magnetic field parameter (M), Williamson parameter (β), porosity parameter (K), buoyancy forces ratios (A), angle inclination parameter (Ω) mixed convection parameter (λ), Biot numbers due to heat transfer (B 2 ), thermal radiation (Rd), bioconvection Schmidt number (Sc 2 ), the Biot numbers due to mass transfer (B 3 ) and the bioconvection Peclet number (Pe) on the dimensionless velocity, thermal, concentration, microorganism, mass transfer rate (Sh), heat transfer rate (Nu), skin friction ..coefficient, and motile microorganisms transfer rate (Nn), on entropy generation rate via friction variation (N G f ), entropy generation rate via heat variation (N Gh ), the irreversibility ratio (Φ), the total entropy generation rate (N G ) and the Bejan number (Be) are revealed in Figures [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]. inclination in M creates a resistive force to the fluid motion which automatically declines the velocity gradient. Further, the rate of declination in the velocity gradient is faster in the presence of a porous medium when compared to the absence of a porous medium. Physically, the presence of a porous medium creates the resistance to the fluid motion due to which fluid velocity decreases. The sway of mixed convection parameter on the velocity gradient in the presence and absence of Williamson parameter is displayed in Figure 2b. It is seen from the plotted figure that boost up values of λ escalate the velocity gradient. Materially, upsurge in λ produces a higher thermal buoyancy force, which is the reason for the inclination in the velocity of the liquid. Further, rate of inclination is slightly faster for the case β = 0 when compared to β = 0.2. Figure 3a displays the sway of buoyancy forces ratios on velocity gradient in the presence and absence of thermal radiation parameter. One can notice from the plotted figure that, boost up values of A declines the velocity of the fluid motion. Further, the rate of declination in the velocity gradient is faster in the presence of thermal radiation parameter. Figure 3b portrays the impact of angle inclination parameter on velocity gradient for the case B 1 = 0 and B 1 = ∞. The result reveals that declination in angle inclination parameter declines the velocity of the fluid motion. Further, the rate of declination in the velocity gradient is little bit faster for the case B 1 = ∞ when compared to B 1 = 0.  Figure 2a portrays the influence of on velocity gradient for dual values of porosity parameter. Here, the escalation in values of decays the velocity gradient. Physically, inclination in M creates a resistive force to the fluid motion which automatically declines the velocity gradient. Further, the rate of declination in the velocity gradient is faster in the presence of a porous medium when compared to the absence of a porous medium. Physically, the presence of a porous medium creates the resistance to the fluid motion due to which fluid velocity decreases. The sway of mixed convection parameter on the velocity gradient in the presence and absence of Williamson parameter is displayed in Figure 2b. It is seen from the plotted figure that boost up values of escalate the velocity gradient. Materially, upsurge in produces a higher thermal buoyancy force, which is the reason for the inclination in the velocity of the liquid. Further, rate of inclination is slightly faster for the case = 0 when compared to = 0.2. Figure 3a displays the sway of buoyancy forces ratios on velocity gradient in the presence and absence of thermal radiation parameter. One can notice from the plotted figure that, boost up values of declines the velocity of the fluid motion. Further, the rate of declination in the velocity gradient is faster in the presence of thermal radiation parameter. Figure 3b portrays the impact of angle inclination parameter on velocity gradient for the case = 0 and = ∞. The result reveals that declination in angle inclination parameter declines the velocity of the fluid motion. Further, the rate of declination in the velocity gradient is little bit faster for the case = ∞ when compared to = 0.

Heat Transfer Characteristics
Figure 4a reveals the sway of +M on the thermal gradient for dual values of porosity parameter. It is observed from the figure that escalation in the values of M improves the thermal gradient. Further, the rate of inclination in the thermal gradient is faster in the presence of a porosity meter. Physically, the presence of a porous medium with inclination in magnetic parameter produces resistance to the fluid motion, which automatically upsurges the thermal gradient. The sway of the thermal radiation parameter on the thermal gradient for the case Pr = 0 and Pr = 1 is illustrated in Figure 4b. Here, the heightening of Rd heightens the thermal gradient. Physically, inclined values of Rd produce inner heat, which is the reason for the inclination of the thermal gradient. Here, the rate of inclination in the thermal gradient of the fluid is faster for a lower Pr value case when compared to a higher Pr value case. The sway of λ on the thermal gradient for the cases B 1 = 0.5 and B 1 = 1 is displayed in Figure 5a. It is detected from the plotted graph that the rise in values of λ improves the thermal gradient. Further, the rate of inclination is faster for case B 1 = 1 when compared to B 1 = 0.5. Figure 5b portrays the dual solutions for the impact of Ec on the thermal gradient. It is noticed in the plotted figure that escalating values of Ec decline the thermal gradient. Moreover, the rate of declination in heat transfer is faster for the case N t = 0.1 when compared to N t = 0.5.

Heat Transfer Characteristics
Figure 4a reveals the sway of + on the thermal gradient for dual values of porosity parameter. It is observed from the figure that escalation in the values of improves the thermal gradient. Further, the rate of inclination in the thermal gradient is faster in the presence of a porosity meter. Physically, the presence of a porous medium with inclination in magnetic parameter produces resistance to the fluid motion, which automatically upsurges the thermal gradient. The sway of the thermal radiation parameter on the thermal gradient for the case = 0 and = 1 is illustrated in Figure 4b. Here, the heightening of heightens the thermal gradient. Physically, inclined values of produce inner heat, which is the reason for the inclination of the thermal gradient. Here, the rate of inclination in the thermal gradient of the fluid is faster for a lower value case when compared to a higher value case. The sway of on the thermal gradient for the cases = 0.5 and = 1 is displayed in Figure 5a. It is detected from the plotted graph that the rise in values of improves the thermal gradient. Further, the rate of inclination is faster for case = 1 when compared to = 0.5. Figure 5b portrays the dual solutions for the impact of on the thermal gradient. It is noticed in the plotted figure that escalating values of decline the thermal gradient. Moreover, the rate of declination in heat transfer is faster for the case = 0.1 when compared to = 0.5.     Figure 6a portrays the multiple solution for the influence of M on mass transfer. One can notice from the plotted figure that escalation in M escalates the concentration gradient. Further, the rate of inclination is faster in the presence of a porous medium. The sway of B 2 on the concentration gradient for dual cases of δ is depicted in Figure 6b. It is detected from graph that a rise in the values of B 2 improves the mass transfer. The further rate of inclination is slower in the presence of the chemical reaction rate parameter. The multiple solutions for the fluctuation in the concentration gradient for diverse values of Schmidt number are depicted in Figure 7a. Here, a rise in the values of the Schmidt number declines the mass transfer characteristics. In the physical fact of sight, an increase in the Sc lessens the molecular diffusivity and it results in the declination of the concentration gradient. Further, the rate of declination in mass transfer is faster for case N b = 0.2 when compared to N b = 0.1. Figure 7b portrays the impact of N t on the concentration gradient in both the presence and absence of A. One can notice from the graph that the upsurge in N t improves the mass transfer. Further, the rate of inclination in the concentration gradient is faster in the presence of A when compared to the absence of A. , on temperature. Figure 6a portrays the multiple solution for the influence of on mass transfer. One can notice from the plotted figure that escalation in escalates the concentration gradient. Further, the rate of inclination is faster in the presence of a porous medium. The sway of on the concentration gradient for dual cases of is depicted in Figure 6b. It is detected from graph that a rise in the values of improves the mass transfer. The further rate of inclination is slower in the presence of the chemical reaction rate parameter. The multiple solutions for the fluctuation in the concentration gradient for diverse values of Schmidt number are depicted in Figure 7a. Here, a rise in the values of the Schmidt number declines the mass transfer characteristics. In the physical fact of sight, an increase in the Sc lessens the molecular diffusivity and it results in the declination of the concentration gradient. Further, the rate of declination in mass transfer is faster for case = 0.2 when compared to = 0.1. Figure 7b portrays the impact of on the concentration gradient in both the presence and absence of . One can notice from the graph that the upsurge in improves the mass transfer. Further, the rate of inclination in the concentration gradient is faster in the presence of when compared to the absence of .

Microorganism Density Characteristics
The sway of the bioconvection Peclet number on the microorganism profile in the absence and the presence of magnetic parameter is portrayed in Figure 8a

Microorganism Density Characteristics
The sway of the bioconvection Peclet number on the microorganism profile in the absence and the presence of magnetic parameter is portrayed in Figure 8a. Here, a rise in the values of Pe declines the Microorganism density characteristics. The rate of declination in the microorganism profile is faster in the absence of the magnetic parameter when compared to the presence of M. The influence of bioconvection Schmidt number on microorganism profile for the cases B 3 = 5 and B 3 = ∞ is illustrated in Figure 8b. Here, the inclined values of the bioconvection Schmidt number declines the microorganism profile. The rate of declination in the microorganism profile is faster for case B 3 = 5 when compared to B 3 = ∞.

Microorganism Density Characteristics
The sway of the bioconvection Peclet number on the microorganism profile in the absence and the presence of magnetic parameter is portrayed in Figure 8a. Here, a rise in the values of declines the Microorganism density characteristics. The rate of declination in the microorganism profile is faster in the absence of the magnetic parameter when compared to the presence of . The influence of bioconvection Schmidt number on microorganism profile for the cases = 5 and = ∞ is illustrated in Figure 8b. Here, the inclined values of the bioconvection Schmidt number declines the microorganism profile. The rate of declination in the microorganism profile is faster for case = 5 when compared to = ∞.

Skin Friction Coefficient, Nusselt Number, Sherwood Number and the Density Number of the Motile Microorganisms
The influence of on / versus in the presence and absence of the Williamson parameter is represented in Figure 9a. Here, the heightening of decays the friction factor. Moreover, the rate of declination in skin friction is faster in the absence of and skin friction acts as a decreasing function of the magnetic parameter. The sway of on / versus for dual values of Ω is illustrated in Figure 9b. Here, the inclined values of improve the skin friction coefficient. Further, the rate of inclination in the friction factor is faster for case Ω = π/6 when compared to Ω = π/3 and the friction factor acts as a decreasing function of . The dual solutions for the influence of K on

Skin Friction Coefficient, Nusselt Number, Sherwood Number and the Density Number of the Motile Microorganisms
The influence of K on Re r 1/2 C f versus M in the presence and absence of the Williamson parameter is represented in Figure 9a. Here, the heightening of K decays the friction factor. Moreover, the rate of declination in skin friction is faster in the absence of β and skin friction acts as a decreasing function of the magnetic parameter. The sway of λ on Re r 1/2 C f versus A for dual values of Ω is illustrated in Figure 9b. Here, the inclined values of λ improve the skin friction coefficient. Further, the rate of inclination in the friction factor is faster for case Ω = π/6 when compared to Ω = π/3 and the friction factor acts as a decreasing function of A. The dual solutions for the influence of K on Re r −1/2 Nu versus M are portrayed in Figure 10a. One can notice from plotted figure that an upsurge in K declines the rate of heat transfer. The rate of declination in the Nusselt number is faster for case β = 0.1 when compared to β = 0.3 and the Nusselt number acts as a decreasing function of the magnetic parameter. The sway of Rd on Re r −1/2 Nu versus A for cases B 2 = 0.1 and B 2 = ∞ is portrayed in Figure 10b. It can be detected from the plotted figure that a rise in values of Rd improves the rate of heat transfer.    Furthermore, the rate of inclination in the heat transfer rate is faster for case B 2 = 0.1 when compared to B 2 = ∞. The twin solutions for the influence of K on Re r −1/2 Sh versus M are portrayed in Figure 11a. Here, an escalation in K increases the mass transfer rate. Further, the rate of inclination in the mass transfer rate is faster for lower values of the Williamson parameter and also the Sherwood number acts as a growing function of M. Figure 11b displays the sway of the Schmidt number on the rate of mass transfer versus N t for cases B 2 = 1 and B 2 = ∞. Here, the boost up values of the Schmidt number decline the mass transfer rate. Furthermore, the rate of declination in the mass transfer rate is faster for case B 2 = ∞ when compared to B 2 = 1 and the Sherwood number acts as a decreasing function of A. The sway of the porosity parameter on the microorganisms transfer rate versus the magnetic parameter in the presence and absence of the Williamson parameter is displayed in Figure 12a. Here, an inclination in K improves the rate of the microorganism's transfer. Further, the rate of inclination in the microorganism's transfer rate is faster in the absence of β and the microorganism transfer rate acts as a growing function of M. The sway of Pe on the microorganism transfer rate versus Sc 2 for cases B 3 = 1 and B 3 = ∞ is depicted in Figure 12b. One can notice from the plotted graph that a rise in the values of Pe improves the microorganism transfer rate and the rate of inclination in the microorganism transfer rate is faster for case B 3 = ∞ when compared to B 3 = 1. Further, the microorganism transfer rate acts as an increasing function of Sc 2 . function of . The sway of on the microorganism transfer rate versus for cases = 1 and = ∞ is depicted in Figure 12b. One can notice from the plotted graph that a rise in the values of improves the microorganism transfer rate and the rate of inclination in the microorganism transfer rate is faster for case = ∞ when compared to = 1. Further, the microorganism transfer rate acts as an increasing function of .

Entropy Generation, Bejan Number and Irreversibility Ratio
The influence of on the entropy generation rate due to heat variation versus in the presence and absence of the porosity parameter is portrayed in Figure 13a. It is seen from the figure that the upsurge in improves the . Further, the rate of inclination in is faster in the absence of , and acts as a decreasing function of . The sway of on the entropy generation rate due to heat variation versus for the cases = 1 and = 2 is portrayed in Figure 13b. It is observed from the plotted graph that the boost up values of improve the . Further, the rate of inclination in is faster for higher values. Here, acts as a growing function of . The influence of on the entropy generation rate due to friction action versus in the presence and absence of porosity parameter is portrayed in Figure 14a. It is seen from the figure that an upsurge in improves the . Further, the rate of inclination in is faster in the presence of and acts as an increasing function of . The sway of on the entropy generation rate due to friction action versus for cases = 1 and = 2 is portrayed in

Entropy Generation, Bejan Number and Irreversibility Ratio
The influence of β on the entropy generation rate due to heat variation versus M in the presence and absence of the porosity parameter is portrayed in Figure 13a. It is seen from the figure that the upsurge in β improves the N Gh . Further, the rate of inclination in N Gh is faster in the absence of K, and N Gh acts as a decreasing function of M. The sway of Rd on the entropy generation rate due to heat variation versus Ec for the cases Pr = 1 and Pr = 2 is portrayed in Figure 13b. It is observed from the plotted graph that the boost up values of Rd improve the N Gh . Further, the rate of inclination in N Gh is faster for higher Pr values. Here, N Gh acts as a growing function of Ec. The influence of β on the entropy generation rate due to friction action versus M in the presence and absence of porosity parameter is portrayed in Figure 14a. It is seen from the figure that an upsurge in β improves the N G f . Further, the rate of inclination in N G f is faster in the presence of K and N G f acts as an increasing function of M. The sway of Rd on the entropy generation rate due to friction action versus Ec for cases Pr = 1 and Pr = 2 is portrayed in Figure 14b. It is observed from the plotted graph that the boost up values of Rd declines the N G f . Further, the rate of declination in N G f is slower for higher Pr values. Here, N Gh acts as a growing function of Ec. The sway of β on the entropy generation profile in the presence and absence of K is portrayed in Figure 15a. Here, inclination in β escalates the entropy generation rate. The rate of inclination in the entropy generation rate is faster in the presence of K and N G acts as an increasing function of M. The control of Rd on the entropy generation rate versus Ec for cases Pr = 1 and Pr = 2 is portrayed in Figure 15b. It is observed from the plotted graph that the boost up values of Rd improve the N G . Moreover, the rate of inclination in N G is faster for higher Pr values and N G acts as a growing function of Ec. . Further, the rate of declination in is slower for higher values. Here, acts as a growing function of . The sway of on the entropy generation profile in the presence and absence of is portrayed in Figure 15a. Here, inclination in escalates the entropy generation rate. The rate of inclination in the entropy generation rate is faster in the presence of and acts as an increasing function of . The control of on the entropy generation rate versus for cases = 1 and = 2 is portrayed in Figure  15b. It is observed from the plotted graph that the boost up values of improve the . Moreover, the rate of inclination in is faster for higher values and acts as a growing function of .    The influence of on versus in the presence and absence of the porosity parameter is portrayed in Figure 16a. it is seen from the figure that an upsurge in declines the Bejan number. Further, the rate of declination in is faster in the presence of and acts as a decreasing function of . The power of on the Bejan number versus for cases = 1 and = 2 is portrayed in Figure 16b. It is observed from the plotted graph that the boost up values of improve the . Further, the rate of inclination in is faster for higher values and acts as a growing function of . The sway of on the irreversibility ratio in the presence and absence of is portrayed in Figure 17a. Here, the inclination in escalates the irreversibility ratio. Further, the rate of inclination in the irreversibility ratio rate is faster in the presence of and Φ acts as a cumulative function of . The sway of on the irreversibility ratio versus for the cases = 1 and = 2 is portrayed in Figure 17b. It is observed from the plotted graph that the boost up values of declines the Φ. Further, the rate of declination in Φ is faster for higher values. Here, Φ acts as a decreasing function of . Figure 18a portrays the streamlines pattern. Generally, streamlines are the path of imaginary particles suspended in the liquid and carried along with it. Figure 18b portrays the isotherms pattern. Physically, an isotherm is the relationship between the fluid and concentrations of a solid particle, used to describe states of no change in the sorption process. The influence of β on Be versus M in the presence and absence of the porosity parameter is portrayed in Figure 16a. it is seen from the figure that an upsurge in β declines the Bejan number. Further, the rate of declination in Be is faster in the presence of K and Be acts as a decreasing function of M. The power of Rd on the Bejan number versus Ec for cases Pr = 1 and Pr = 2 is portrayed in Figure 16b. It is observed from the plotted graph that the boost up values of Rd improve the Be. Further, the rate of inclination in Be is faster for higher Pr values and N G acts as a growing function of Ec. The sway of β on the irreversibility ratio in the presence and absence of K is portrayed in Figure 17a. Here, the inclination in β escalates the irreversibility ratio. Further, the rate of inclination in the irreversibility ratio rate is faster in the presence of K and Φ acts as a cumulative function of M. The sway of Rd on the irreversibility ratio versus Ec for the cases Pr = 1 and Pr = 2 is portrayed in Figure 17b. It is observed from the plotted graph that the boost up values of Rd declines the Φ. Further, the rate of declination in Φ is faster for higher Pr values. Here, Φ acts as a decreasing function of Ec. Figure 18a portrays the streamlines pattern. Generally, streamlines are the path of imaginary particles suspended in the liquid and carried along with it. Figure 18b portrays the isotherms pattern. Physically, an isotherm is the relationship between the fluid and concentrations of a solid particle, used to describe states of no change in the sorption process. and = 2 is portrayed in Figure 17b. It is observed from the plotted graph that the boost up values of declines the Φ. Further, the rate of declination in Φ is faster for higher values. Here, Φ acts as a decreasing function of . Figure 18a portrays the streamlines pattern. Generally, streamlines are the path of imaginary particles suspended in the liquid and carried along with it. Figure 18b portrays the isotherms pattern. Physically, an isotherm is the relationship between the fluid and concentrations of a solid particle, used to describe states of no change in the sorption process.

Conclusions
In this article, a numerical simulation for an incompressible, two-dimensional, steady, mixed convective flow of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium under the influence of gyrotactic microorganisms and thermal radiation is examined. Applications of this scrutiny are comprised of magnetohydrodynamic propulsion thermo-fluid dynamics materials processing and laminar magneto-aerodynamics. The characteristics of velocity, thermal, concentration, entropy generation and microorganisms' profiles have been deliberated under the influence of in-

Conclusions
In this article, a numerical simulation for an incompressible, two-dimensional, steady, mixed convective flow of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium under the influence of gyrotactic microorganisms and thermal radiation is examined. Applications of this scrutiny are comprised of magnetohydrodynamic propulsion thermo-fluid dynamics materials processing and laminar magneto-aerodynamics. The characteristics of velocity, thermal, concentration, entropy generation and microorganisms' profiles have been deliberated under the influence of in-

Conclusions
In this article, a numerical simulation for an incompressible, two-dimensional, steady, mixed convective flow of Williamson nanofluid along an inclined semi-infinite plate embedded in a porous medium under the influence of gyrotactic microorganisms and thermal radiation is examined. Applications of this scrutiny are comprised of magnetohydrodynamic propulsion thermo-fluid dynamics materials processing and laminar magnetoaerodynamics. The characteristics of velocity, thermal, concentration, entropy generation and microorganisms' profiles have been deliberated under the influence of involved flow controlling physical parameters. The main outcomes of the present investigation are as follows: The escalation in values of M declines the velocity gradient and improves the thermal and concentration gradients. The inclination in radiation parameter improves the thermal gradient. The rise in values of the Schmidt number declines the mass transfer characteristics but the converse trend is depicted for the boost up values of N t . The escalating values of Ec decline the thermal gradient. The rise in values of Pe and Sc 2 declines the microorganism density characteristics. The inclination in β and Rd escalates the entropy generation rate. The upsurge in β declines the Bejan number and the irreversibility ratio but the converse trend is depicted in both the profiles for inclined values of Rd. The upsurge in K declines the rate of heat transfer but the converse trend is depicted for inclined values of Rd. The escalation in K increases the mass transfer rate but the converse trend is depicted for inclined values of Sc 1 . The inclination in K and Pe improves the rate of a microorganism's transfer.