## 1 Introduction

In fluid dynamics, bioconvection [1,2,3] occurs when microorganisms, which are denser than water, swim upwards. The upper surface of the fluid becomes thicker due to the assemblage of microorganisms. As a result, the upper surface becomes unstable and microorganisms fall down, which creates bioconvection. Bioconvection continues to be explored widely because of its enormous applications in the field of pharmaceutical industry, purification of cultures, microfluidic devices, mass transport enhancement and mixing, microbial enhanced oil recovery and enzyme biosensors. Bioconvection systems could be categorized based on the directional motion of different species of microorganisms. In particular, gyrotactic microorganisms are the ones whose swimming direction is dependent on a balance between gravitational and viscous torques [4,5]. Oyelakin et al. [6] pondered the impact of bioconvection and motile gyrotactic microorganisms on the Casson nanofluid past a stretching sheet and observed that the microorganism profile decreases as a result of an increment in the Peclet number. Saini and Sharma [7] explored the effects of bioconvection and gyrotactic microorganisms on the nanofluid flow over a porous stretching sheet. It is noted that the Lewis number escalates the bioconvection process. Dhanai et al. [8] explored the impact of bioconvection on the fluid flow over an inclined stretching sheet and assessed that the microorganism density profile is enhanced with an improvement in the bioconvection Schmidt number. Mahdy [9] pondered the effects of motile microorganisms on the fluid past a stretching wedge and noted that a positive variation in the Peclet number leads to an augmentation in the microorganism profile. Avinash et al. [10] pondered the impact of bioconvection and aligned magnetic field on the nanofluid flow over a vertical plate and concluded that the heat transfer rate increases with an improvement in the Lewis number. Makinde and Animasaun [11] studied the effects of magnetohydrodynamics (MHD), bioconvection, nonlinear thermal radiation and nanoparticles on fluid past an upper horizontal surface of a paraboloid of revolution and found that the Brownian motion boosts the concentration profile. Khan et al. [12] studied the impact of MHD, gyrotactic microorganisms, slip condition and nanoparticles on the fluid flow over a vertical stretching plate; it was observed that the magnetic field suppresses the dimensionless velocity inside the boundary layer. Later, the effects of different features of the gyrotactic microorganisms on the fluid flow are analyzed in various investigations [13,14,15].

Nanotechnology has been considered the most substantial and fascinating forefront area in physics, engineering, chemistry and biology. The thermal conductivity of a nanofluid is greater than that of the base fluid. The thermal conductivity of the fluid is considered to be enhanced by the nanoparticles present in the fluid. Buongiorno [16] established a model to examine the thermal conductivity of nanofluids. Baby and Ramaprabhu [17] analyzed the heat transport of fluids using graphene nanoparticles. They reported that the thermal conductivity of hydrogen-exfoliated graphene is enhanced with an increment in the volume fraction of the nanoparticles. Khan and Gorla [18] pondered the mass transfer of the nanofluid flow over a convective sheet using the Keller box scheme and noted that the heat transfer rate is high in the dilatant fluids compared with that in the pseudoplastic fluids. Das [19] discussed the rotating flow of a nanofluid with respect to the constant heat source. A boost in the volume fraction of nanoparticles was observed to cause an increment in the thermal boundary layer thickness. Gireesha et al. [20] considered the Hall impact on a dusty nanofluid and concluded that the skin friction coefficient decreases due to an improvement in the Hall current.

The experimental and the theoretical scientific studies of the non-Newtonian liquids together with MHD have achieved a considerable attention of researchers because of their adequate applications in the field of aeronautics, chemical, mechanical, civil and bio-engineering. The fluid becomes electrically conducting under the effect of MHD like ionized gases, plasmas and liquid metals such as mercury. The impact of MHD and nonlinear thermal radiation on the Sisko nanofluid flow over a nonlinear stretching surface is premeditated by Prasannakumara et al. [21]. Rashidi et al. [22] pondered the MHD viscoelastic fluid together with the Soret and Dufour effects and observed that the velocity profile decreases with an improvement in the magnetic parameter. Kothandapani and Prakash [23] studied the effect of magnetic field on peristaltic tangent hyperbolic nanofluid past a asymmetric channel. Gaffar et al. [24] showed the tangent hyperbolic fluid flow over a cylinder together with the MHD and partial slip effects. Nagendramma et al. [25] analyzed the tangent hyperbolic fluid flow over a stretching sheet together with the MHD effect. Das et al. [26] investigated the impact of magnetic field, chemical reaction and double-diffusive convection on the Casson fluid flow past a stretching plate and noted that the skin friction coefficient decreases as a result of an augmentation in the Grashof number. Sravanthi and Gorla [27] examined the effect of the Maxwell nanofluid flow over an exponentially stretching sheet together with MHD, chemical reaction and heat source/sink.

Double-diffusion phenomena describe a form of convection driven by two different density gradients, holding distinct rate of diffusion. Double-diffusive convection occurs in a variety of scientific disciplines such as oceanography, biology, astrophysics, geology, crystal growth and chemical reactions [28]. Nield and Kuznetsov [29] scrutinized the nanofluid past a porous medium along with the double-diffusive convection effect. The impact of double-diffusive convection on the fluid flow over a square cavity is analyzed by Mahapatra et al. [30]. Gireesha et al. [31] discussed the Casson nanofluid past a stretching sheet along with the MHD and double-diffusive convection. Rana and Chand [32] explored the effect of double-diffusive convection on viscoelastic fluid and deduced that a Rayleigh number increases with an improvement in the Soret parameter. Gaikwad et al. [33] have monitored the fluid flow above a stretching sheet together with double-diffusive convection and found that an augmentation in the Nusselt number takes place with an improvement in the Dufour parameter. Kumar et al. [34] inspected the influence of nanoparticles and double diffusion on viscoelastic fluid and monitored that an increase in the velocity field occurs with an increment in the Dufour Lewis number.

Convection is a process common to particles, gases and vapours. Convection occurs when a fluid is in motion and that motion carries with it a material of interest such as the particles or the droplets of an aerosol. There are two types of convection: free convection and forced convection. In free convection or natural convection, the fluid motion cannot led by external sources such as fans, pumps, and suction devices etc. Gravity is the main driving force in the case of free convection. Free convection has various environmental and industrial applications such as plate tectonics, oceanic currents, formation of microstructures during the cooling of molten metals, fluid flows around shrouded heat dissipation fins, solar ponds and free air cooling without the aid of fans. In forced convection, the fluid motion is generated externally with the help of pumps, fans, suction devices, etc. This mechanism has enormous applications in our daily life such as heat exchangers, central heating system, steam turbines and air conditioning. Mixed convection is the situation in which both free convection and forced convection are of comparable order. Mixed convection is of great interest to researchers due to its enormous applications in the industrial and engineering sectors. Ibrahim and Gamachu [35] found the numerical solution of the mixed convective Williamson nanofluid past a stretching sheet by the Galerkin finite element method. Shateyi and Marewo [36] adopted the spectral quasi-linearization method to achieve the numerical solution of the mixed convective magneto Jeffrey fluid flow over an exponentially stretching sheet together with the thermal radiation and observed that the fluid velocity improves with an augmentation in the buoyancy parameter. Nalinakshi et al. [37] found the numerical solution of the mixed convective fluid past a vertical stretched plate using a nonlinear shooting method. El-Aziz and Tamer Nabil [38] gave the numerical solution for the problem of the MHD and Hall current effect on mixed convective fluid past a stretching sheet using the homotopy analysis method (HAM) and noted that a positive variation in the Hall current parameter leads to an increase in the velocity field. Beg et al. [39] employed an explicit finite difference scheme to yield the solution of the magneto mixed convection nanofluid flow over a stretchable surface under the effect of MHD and viscous dissipation. The numerical solution of the gravity-driven Navier–Stokes equation has been reported by Zhang et al. using a finite difference method [40]. Pal and Chatterjee [41] studied the impact of the Soret and Dufour effects along with nonlinear thermal radiation on the double-diffusive convective fluid past a stretchable surface and achieved the numerical solution for problem using the Runge–Kutta–Fehlberg method along with the shooting scheme. They noted that the velocity field increases with an enhancement in the Grashof number.

The aim of this study was to construct a mathematical model that describes a form of convection driven by two different density gradients, which have different rates of diffusion (double-diffusive convection). So far, no reviews have been reported on the non-Newtonian fluid past a stretching sheet embedded with nanoparticles, double-diffusive convection and motile gyrotactic microorganisms.

## 2 Mathematical formulation

Figure 1 displays the effect of tangent hyperbolic nanofluid past a stretching sheet with stretching velocity *u*_{w} = *ax* along the *x*-axis. When the Reynolds number is assumed to be small, the induced magnetic field can be neglected compared with the applied magnetic field *B*_{0}, which is applied transversely to the surface. *T*_{w}, *γ*_{w}, *C*_{w} and *N*_{w} denote the temperature, solute concentration, concentration of nanoparticles and density of the motile gyrotactic microorganisms at the wall, respectively, whereas *T*_{∞}, *γ*_{∞}, *C*_{∞} and *N*_{∞} denote the ambient temperature, solute concentration, concentration of nanoparticles and density of the motile gyrotactic microorganisms, respectively. The fluid has further been assumed to contain the gyrotactic microorganisms. The microorganisms present in the fluid move towards light. The “bottom heavy” mass of the microorganisms orients its body and enables them to move against the gravity *g*, which is called as gyrotactic phenomena. The presence of microorganisms is considered to be beneficial for the suspension of the nanofluid. To maintain the stability of convection, the motion of microorganisms has been taken, independent of that of the nanoparticles. The double-diffusive fluid flow over a stretching sheet embedded with gyrotactic microorganisms has not been explored yet, and we want to rectify this problem in this study.

The governing equations include some important effects that have eminent involvement in the industries and engineering fields. The momentum equation includes bioconvection and MHD. MHD has been used in many engineering processes such as nuclear reactor, MHD power generation, in which heat energy is directly converted into electrical energy, Yamato-1 boat incorporating a superconductor cooled by liquid helium and microfluidics. A microorganism or microbe is an organism that is so small that it can be seen only through a microscope (invisible to the naked eye). The presence of microorganisms in the fluid becomes the core area of the research during the past decade. The presence of microorganisms in the base fluid causes a “stabilization” or “destabilization” in the motion of nanoparticles. The microorganisms have various applications in genetic engineering, wastewater engineering, agricultural engineering and chemical engineering. The temperature equation and concentration equations are embedded with nanoparticles and double-diffusive convection. Nanoparticles are used to enhance the thermal conductivity of the fluid and used in tissue engineering, mechanical engineering, nanomedicine, environmental engineering, etc. Double diffusion portrays the form of convection conducted by two different density gradients. There are various examples in environmental engineering such as Arctic Ocean study and Lake Kivu, in which magma, sand and materials of different densities are diffused with water. The same situation is applicable in our modelled problem, in which microorganisms and nanoparticles of different densities are diffused together in the fluid. The last governing equation tells us about the impact of gyrotactic microorganisms present in the fluid. Various types of microorganisms such as algae, fungi, protozoa and bacteria are suspended in the fluid. These microorganisms swim in the fluid under the combination of gravitational and viscous torques (gyrotactic) in fluid flow. The gyrotactic microorganisms have enormous contribution to genetic engineering, microbial engineering and soil engineering. Under the usual boundary layer approximations, the equations of conservation of mass, momentum, thermal energy, solute, concentration of nanoparticles and gyrotactic microorganisms take the following forms [11,12,13,31,32,34]:

*ρ*

_{f}” depicts the fluid density and “

*ρ*

_{p}” represents the density of nanoparticles. The similarity transformations [38] are as follows:

Invoking equation (8), equation (1) is automatically satisfied and equations (2)–(6) become:

Distinct physical parameters arising after the conversion of PDEs into ODEs are as follows:

The important quantities of interest like rate of shear stress *C*_{f} and heat as well as mass transfer rates Nu_{x} and Sh_{x} and Sh_{x,n} and Nn_{x} are as follows:

*τ*

_{w},

*q*

_{w},

*q*

_{m},

*q*

_{mn}and

*q*

_{n}are as follows:

By substituting equation (17) into equation (16) and using the similarity transformation, the quantities defined in equation (17) are nondimensionalized as follows:

## 3 Numerical scheme

The dimensionless system of equations (9)–(13) along with the boundary condition (14) should be handled with the help of the numerical scheme called the implicit finite difference method (Keller box technique) [42,43] for distinguished parameters that emerged during numerical simulation of the problem. Such type of differential equations in this article can usually be solved with the help of other numerical techniques such as shooting method, HAM and bvp4c [27,31,32,33,34,44,45,46,47,48,49]. In this study, the standard Keller box method has been used. This numerical technique is quite effective and flexible to solve the parabolic-type boundary value problems of any order, is unconditionally stable and attains remarkable accuracy. The Keller box scheme is numerically more stable and converges using less iterations compared with other numerical techniques. Figure 2 shows the flow chart procedure of the Keller box method. By adopting the new variables *z*_{1}, *z*_{2}, *z*_{3}, *z*_{4}, *z*_{5} and *z*_{6},

The dimensionless equations (9)–(13) are transformed into first-order differential equations (ODEs) as follows:

The transformed boundary conditions are as follows:

Figure 3 portrays the mesh structure for central difference approximations. The stepping procedure for the selection of the nodes in the case of domain discretization is as follows:

The derivatives of equations (20)–(24) are approximated by employing the central difference at the midpoint

After linearization of the above-mentioned system of equations, the subsequent block-tridiagonal block structure:

*A*is the

*j*×

*j*tridiagonal matrix of block size 11 × 11, and

*δ*and

*R*are the column matrices of

*j*rows. Now equation (38) has been tackled using the LU factorization method with lower triangular matrix

*L*and upper triangular matrix

*U*enumerated as follows:

To solve the problem numerically, the domain of the problem has been considered [0,*η*_{max}] instead of [0,∞), where *η*_{max} = 16 and the step size is *h*_{j} = 0.01. All the numerical results achieved in this problem are subjected to an error tolerance of 10^{−5}.

Table 1 displays the comparison analysis of the given numerical scheme results with Ibrahim [40].

## 4 Results and discussion

To discuss the outcomes, the behaviour of various pertinent parameters against the Nusselt number, the Sherwood number, motile density profile, velocity field, temperature field, mass fraction field and solute profile is monitored. Table 2 exhibits the behaviour of distinguished parameters on heat transfer at the boundary, mass fraction field and the motile microorganisms density profile for thermophoresis parameter (Nt) = 0.1, Prandtl number (Pr) = 6.2, Lewis number (Le) = 0.5, Dufour Lewis number (Ld) = 0.1 and mixed convection parameter *Λ* = 0.1. The heat transfer rate diminishes in the case of magnetic parameter *M*, Weissenberg number (We), modified Dufour parameter (Nd), power law index *n*, nanofluid Lewis number (Ln) and buoyancy ratio parameter (Nr), whereas an embellishment in the Nusselt number is seen for the Brownian motion parameter (Nb) and the bioconvection Rayleigh number (Nc). The Nusselt number has shown no variation in the case of microorganism concentration difference parameter *σ*, Peclet number (Pe) and bioconvection Lewis number (Lb). The mass fraction field depreciates in the case of *M*, Nc, nanofluid Lewis number (Ln) and buoyancy ratio parameter (Nr), but a positive variation is observed for Nb, We, Nd and *n*, whereas static behaviour is seen for *σ*, Pe and Lb. Furthermore, the number of motile microorganisms has been seen to increase in the case of positive variation in *M*, *σ*, Pe, Lb and Ln, but the situation is opposite in the case of Nr, Nc, *n*, We, Nd and Nb.

Variation in *Λ* = 0.1 are fixed

M | Nc | Nr | Nb | We | n | σ | Pe | Lb | Nd | Ln | −θ′(0) | −ξ′(0) | −χ′(0) |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|

0.1 | 0.5 | 0.5 | 0.1 | 0.3 | 0.2 | 0.5 | 1 | 1 | 0.1 | 2 | 0.93786 | 1.48950 | 1.30837 |

0.2 | 0.93815 | 1.50424 | 1.32242 | ||||||||||

0.3 | 0.93841 | 1.51812 | 1.33562 | ||||||||||

0.1 | 2.03841 | 2.71812 | 3.45016 | ||||||||||

0.3 | 2.03853 | 2.72591 | 2.63562 | ||||||||||

0.5 | 2.03865 | 2.73339 | 2.64400 | ||||||||||

0.1 | 2.05487 | 3.12893 | 2.93505 | ||||||||||

0.2 | 2.05483 | 3.12731 | 2.93360 | ||||||||||

0.3 | 2.05480 | 3.12569 | 2.93214 | ||||||||||

0.4 | 0.33911 | 7.72999 | 11.6733 | ||||||||||

0.5 | 0.55922 | 7.70119 | 11.6306 | ||||||||||

0.6 | 0.68789 | 7.69032 | 11.6146 | ||||||||||

0.1 | 0.82446 | 4.07381 | 6.23334 | ||||||||||

0.2 | 0.82438 | 4.06941 | 6.22659 | ||||||||||

0.3 | 0.82431 | 4.06476 | 6.21943 | ||||||||||

0.3 | 0.82367 | 4.03450 | 6.17277 | ||||||||||

0.4 | 0.82271 | 3.98592 | 6.09800 | ||||||||||

0.5 | 0.82134 | 3.90893 | 5.97975 | ||||||||||

0.1 | 0.82438 | 4.06941 | 4.69593 | ||||||||||

0.2 | 0.82438 | 4.06941 | 5.07859 | ||||||||||

0.3 | 0.82438 | 4.06941 | 5.46126 | ||||||||||

0.1 | 0.82438 | 4.06941 | 1.03179 | ||||||||||

0.5 | 0.82438 | 4.06941 | 3.31479 | ||||||||||

1 | 0.82438 | 4.06941 | 6.22659 | ||||||||||

0.5 | 0.82438 | 4.06941 | 3.31479 | ||||||||||

1 | 0.82438 | 4.06941 | 6.22659 | ||||||||||

1.5 | 0.82438 | 4.06941 | 9.18276 | ||||||||||

0.1 | 0.89710 | 2.21688 | 3.50681 | ||||||||||

0.2 | 0.81434 | 2.20547 | 3.48796 | ||||||||||

0.3 | 0.78569 | 2.18375 | 3.45366 | ||||||||||

1 | 0.89710 | 2.21688 | 3.50681 | ||||||||||

2 | 0.85497 | 3.26596 | 5.04071 | ||||||||||

3 | 0.82438 | 4.06941 | 6.22659 |

Figure 4 exhibits the effect of the magnetic parameter *M* on the velocity profile *f*′(*η*). It has been found that an increase in *M* decreases the velocity profile. Actually, the resistive force called the Lorentz force is generated due to the application of the magnetic field to the electrically conducting fluid. As a result, the velocity of the fluid reduces. Figure 5 indicates the effect of *n* on the velocity field *f*′(*η*). The parameter decides the viscosity of the fluid or how much viscous the fluid is. The fluid behaves like shear thinning for the case of *n* < 1, shear thickening for the larger values of *n* > 1 and Newtonian in the case of *n* = 1. The velocity of the fluid decreases in the case of *n* > 1, and as a result, the velocity field diminishes. Figure 6 depicts the effect of *f*′(*η*) on We. The Weissenberg number is defined as the ratio of viscous forces to the inertial forces. This parameter is important to study the fluid flow behaviour. The Weissenberg number actually depicts the elastic nature of the fluid. It is noted that the higher values of the Weissenberg number indicate the solid nature of the fluid, while lower values of the Weissenberg number depict the liquid nature of the fluid. It is clear that an augmentation in the Weissenberg number leads to a reduction in the velocity of the fluid. Figure 7 highlights the variation in the temperature profile *θ*(*η*) against the various values of *M* and observed that an electric current in the presence of magnetic field generates a Lorentz force. This force resists the motion of the fluid; hence, additional heat is produced, which enhances the fluid temperature. Figure 8 highlights the behaviour of temperature field *θ*(*η*) against Pr. The Prandtl number is a dimensionless quantity, which is defined as the ratio of momentum diffusivity to thermal diffusivity and has important application in the study of boundary layer concept. The thermal diffusivity dominates in the case of Pr ≪ 1, whereas momentum diffusivity dominates in the case of Pr ≫ 1. It is observed that the fluids with small Prandtl number are free flowing liquids with high thermal conductivity and favourable choice for heat conducting fluids. The thermal conductivity of the fluid decreases with an augmentation in the value of Pr, and the heat transfer decelerates, which decreases the temperature of the flow field, and as a result, a decrease in the temperature is observed.

Figure 9 portrays the effect of Brownian diffusion parameter (Nb) on the temperature distribution *θ*(*η*). Brownian motion is actually the random motion of the particles suspended in the fluid. The temperature of the fluid increases as a result of the random collision of particles suspended in the liquid, which further leads to an expected improvement in the temperature profile *θ*(*η*). Figure 10 explores the effect of the thermophoresis parameter (Nt) on the temperature distribution *θ*(*η*). In the thermophoresis process, smaller particles migrate from the region having high temperature to the region having low temperature, which ultimately causes an improvement in the fluid temperature.

Figure 11 shows the behaviour of temperature profile *θ*(*η*) against the different values of Nd. The situation in which heat and mass transfer happens simultaneously in a moving fluid affecting each other causes a cross-diffusion. The mass transfer caused by temperature gradient is called the Soret effect, whereas the heat transfer caused by concentration is called the Dufour effect. The Dufour number implies the effect of the concentration on the thermal energy flux in the flow. It is found that a variation in the modified Dufour number leads to a monotonic enhancement in the temperature field *θ*(*η*). Figure 12 highlights the effect of Nb on the mass fraction field. Brownian diffusion and thermophoresis parameters emerge as a result of an inclusion of nanoparticles into the fluid. Brownian diffusion and thermophoresis parameters help to understand the motion of the nanoparticles in the fluid. It is verified that the higher values of Nb are the root cause to boost the random motion among the nanoparticles present in the fluid. This results in the decrease in the concentration of the fluid.

Figure 13 describes the effect of Nt on the mass fraction field. It is observed that increasing values of Nt push nanoparticles away from the warm surface. The density of the concentration boundary layer upsurges due to an augmentation in the value of Nt, which leads to an embellishment in the mass fraction field. Figure 14 portrays the effect of the nanofluid Lewis parameter (Ln) on the mass fraction field. The Lewis number is defined as the ratio of thermal diffusivity to the mass diffusivity, and it is the prominent factor to study the heat and mass transfer. It is observed that the concentration profile decreases due to the dependence of the Lewis number on the Brownian diffusion coefficient, which means that an augmentation in the Brownian diffusion coefficient brings about a decrease in the concentration profile and the nanofluid Lewis number.

Figure 15 shows the effect of Peclet number (Pe) on the microrotation distribution *χ*(*η*). The Peclet number is the prominent factor to study the microorganisms swimming in the fluid. The Peclet number is defined as the ratio of maximum cell swimming speed to diffusion of microorganisms. Diffusion is the process in which a substance moves from an area of high concentration to an area of low concentration. It explains the movement of the substances in the fluid. It is found that diffusivity of microorganisms is decreased in the case of an augmentation in Pe. As a result, the microrotation distribution declines. Figure 16 depicts the effect of the bioconvection Lewis number (Lb) on the microrotation distribution. Similar to Figure 14, an augmentation in Lb results in a decrease in the diffusivity of microorganisms, which results in the reduction of the motile density profile.

Figure 17 portrays the effect of microorganism concentration difference parameter *σ* on the motile density profile. It is observed that by increasing the value of *σ*, the concentration of microorganisms in ambient fluid is decreased. Figure 18 delineates the effect of the regular Lewis number (Le) on the solute profile *γ*(*η*). The Lewis number is defined as the ratio of thermal diffusivity to mass diffusivity. As seen in Figure 13, the Lewis number is related to the Brownian diffusion coefficient. It is observed that a positive variation in Brownian diffusion leads to a decrease in the concentration of particles. Thus, a positive variation in the Lewis number (Le) leads to a decrease in the solute profile. Figure 19 portrays the relationship between the Dufour Lewis number (Ld) and the solute profile *γ*(*η*). The Dufour Lewis number depicts the influence of temperature gradient on the concentration field. It is perceived that the concentration gradient excites the flow with an enhancement in the thermal energy, which results in an increase in the solute profile. Figure 20 depicts the effect of mass fraction field on Nb for the distinguished values of the nanofluid Lewis number (Ln). It is also observed that due to the random collision of molecules, the heat transfer process escalates and nanoparticle diffusion reduces, which results in an increment in the Sherwood number.

Figure 21 elucidates the performance of Nt on the mass fraction field for various values of Ln. It is found that in the presence of the thermophoretic force, the nanoparticles present close to the hot boundary have been shifted towards the cold fluid, which decreases the thermal boundary layer and heightens the nanofluid Lewis number. An upsurge in Nt escalates nanofluid Lewis number (Ln) and further leads to an augmentation in the mass fraction field. Figure 22 presents the effect of microorganism concentration difference parameter *σ* on the density number of microrotation distribution for different values of Peclet number (Pe). A positive variation in *σ* lessens the thickness of the boundary layer and leads to an increment in the concentration of the motile gyrotactic microorganisms. Figure 23 elucidates the conduct of the Dufour Lewis number (Ld) on the solutal Sherwood number for different values of the Prandtl number. The Lewis number is defined as the ratio of thermal diffusivity to momentum diffusivity. It is observed that an enhancement in Lewis number drives more heat within the fluid, which brings about an augmentation in the Prandtl number. It is noteworthy that a positive variation in the Dufour Lewis parameter leads to an augmentation in the solutal Sherwood number. Figure 24 elaborates the effect of Lewis number (Le) on the solutal Sherwood number. It has been observed that the solutal Sherwood number increases with an augmentation in the Lewis number.

## 5 Concluding remarks

This article elaborates the effects of nanoparticles and double-diffusive convection along with motile gyrotactic microorganisms on the non-Newtonian fluid past a stretching sheet. To our knowledge, no model has been developed so far to see the impact of gyrotactic microorganisms and double-diffusive convection simultaneously on the non-Newtonian hyperbolic tangent nanofluid, and furthermore, a numerical technique (Keller box) has been used to achieve the numerical solution of the problem. A comparison with the previous literature was made to check the reliability of our proposed numerical scheme. The results are quite promising. Some of the key findings of the present investigation are as follows:

- An improvement in the Weissenberg number (We) leads to a decrease in the velocity profile.
- The mass fraction field shows an opposite behaviour as a result of variation in the nanofluid Lewis number (Ln).
- A positive variation in the Peclet number (Pe) leads to a decrease in the solute profile.
- The microrotation distribution profile declines with an improvement in the bioconvection Lewis number (Lb) and microorganism concentration difference parameter
*σ*. - The solute profile is decreased with an enhancement in the regular Lewis number (Le).

## Nomenclature

a | stretching rate |

B_{0} | magnetic field strength |

C_{∞} | ambient concentration |

C_{∞} | ambient solute concentration at the wall |

C_{f} | skin friction coefficient |

C_{p} | specific heat |

C_{w} | solute concentration at the wall |

D_{B} | Brownian diffusion |

D_{m} | diffusivity of the microorganisms |

D_{T} | thermophoresis diffusion |

g | gravity |

G_{T} | Grashof number |

k | thermal conductivity |

Lb | bioconvection Lewis number |

Ld | Dufour Lewis number |

Le | Lewis number |

Ln | nanofluid Lewis number |

M | magnetic parameter |

n | power law index |

N_{∞} | ambient density of the motile microorganisms |

Nb | Brownian motion |

Nd | modified Dufour parameter |

Nt | thermophoresis parameter |

N_{w} | density of the motile microorganisms at the wall |

Pe | Peclet number |

Pr | Prandtl number |

q_{m} | wall mass flux |

q_{w} | wall heat flux |

Re | Reynolds number |

T_{∞} | ambient temperature |

T_{w} | wall temperature |

u,v | velocity components |

U_{w} | stretching velocity along the |

W_{c} | cell swimming speed |

We | Weissenberg number |

α | thermal diffusivity |

γ | solute profile |

μ | dynamic viscosity |

ν | kinematic viscosity |

Λ | mixed convection parameter |

ρ | density of the fluid |

ρ_{m} | density of the microorganisms |

ρ_{p} | density of the nanoparticles |

σ | microorganism concentration difference |

σ_{1} | electric conductivity |

τ_{w} | wall shear stress |

ϕ | volume fraction of the nanoparticles |

## References

- [2]↑
Nield DA, Kuznetsov AV. The onset of bio-thermal convection in a suspension of gyrotactic microorganisms in a fluid layer: oscillatory convection. Int J Therm Sci. 2006;45:990–7.

- [3]↑
Alloui Z, Nguyen TH, Bilgen E. Numerical investigation of thermo- bioconvection in a suspension of gravitactic microorganisms. Int J Heat Mass Transf. 2007;50:1435–41.

- [4]↑
Sokolov A, Goldstein RE, Feldchtein FI, Aranson IS. Enhanced mixing and spatial instability in concentrated bacterial suspensions. Phys Rev E. 2009;80(3):031903.

- [5]↑
Tsai T, Liou D, Kuo L, Chen P. Rapid mixing between ferro-nano fluid and water in a semi-active Y-type micromixer. Sens Actuators, A. 2009;153(2):267–73.

- [6]↑
Oyelakin IS, Mondal S, Sibanda P, Sibanda D. Bioconvection in Casson nanofluid flow with gyrotactic microorganisms and variable surface heat flux. Int J Biomath. 2019;12:1950041.

- [7]↑
Saini S, Sharma YD, Double-diffusive bioconvection in a suspension of gyrotactic microorganisms saturated by nanofluid. Int J Appl Fluid Mech. 2019;12:271–80.

- [8]↑
Dhanai R, Rana P, Kumar L. Lie group analysis for bioconvection MHD slip flow and heat transfer of nano fluid over an inclined sheet: multiple solutions. J Taiwan Inst Chem Eng. 2016;66:283–91.

- [9]↑
Mahdy A. Gyrotactic microorganisms mixed convection nanofluid flow along an isothermal vertical wedge in porous media. Int J Mech Aero Indus Mechatron Manuf Eng. 2017;11(4):840–50.

- [10]↑
Avinash K, Sandeep N, Makinde OD, Animasaun IL. Aligned magnetic field effect on radiative bioconvection flow past a vertical plate with thermophoresis and Brownian motion. Defect Diffus Forum. 2017;377:127–40.

- [11]↑
Makinde OD, Animasaun IL. Thermophoresis and Brownian motion effects on MHD bioconvection of nanofluid with nonlinear thermal radiation and quartic chemical reaction past an upper horizontal surface of a paraboloid of revolution. J Mol Liq. 2016;21:733–43.

- [12]↑
Khan WA, Makinde OD, Khan ZH. MHD boundary layer flow of a nanofluid containing gyrotactic microorganisms past a vertical plate with Navier slip. Int J Heat Mass Transf. 2014;74:285–91.

- [13]↑
Mutuku WN, Makinde OD. Hydromagnetic bioconvection of nanofluid over a permeable vertical plate due to gyrotactic microorganisms. Comput Fluids. 2014;95:88–97.

- [14]↑
Khan WA, Makinde OD. MHD nanofluid bioconvection due to gyrotactic microorganisms over a convectively heat stretching sheet. Int J Therm Sci. 2014;81:118–24.

- [15]↑
Makinde OD, Animasaun IL. Bioconvection in MHD nanofluid flow with nonlinear thermal radiation and quartic autocatalysis chemical reaction past an upper surface of a paraboloid of revolution. Int J Therm Sci. 2016;109:159–71.

- [17]↑
Baby TT, Ramaprabhu S. Enhanced convective heat transfer using graphene dispersed nanofluids. Nanoscale Res Lett. 2011;6:289.

- [18]↑
Khan WA, Gorla RSR. Heat and mass transfer in power-law nanofluids over a nonisothermal stretching wall with convective boundary condition. Int J Heat Mass Trans. 2012;134:112001.

- [19]↑
Das K. Flow and heat transfer characteristics of nanofluids in a rotating frame. Alexandria Eng J. 2014;53:757–66.

- [20]↑
Gireesha BJ, Mahanthesh B, Thammanna GT, Sampathkumar PB. Hall effects on dusty nano fluid two-phase transient flow past a stretching sheet using KVL model. J Molec Liq. 2018;256:139–47.

- [21]↑
Prasannakumara BC, Gireesha BJ, Krishnamurthy MR, Ganesh KK. MHD flow and nonlinear radiative heat transfer of Sisko nanofluid over a nonlinear stretching sheet. Info Med Unlocked. 2017;99:123–32.

- [22]↑
Rashidi MM, Ali M, Rostami B, Rostami P, Xie GN. Heat and mass transfer for MHD viscoelastic fluid flow over a vertical stretching sheet with considering Soret and Dufour effects. Math Probl Eng. 2015;2015:1–12.

- [23]↑
Kothandapani K, Prakash J. Influence of heat source, thermal radiation, and inclined magnetic field on peristaltic flow of a hyperbolic tangent nanofluid in a tapered asymmetric channel. IEEE Trans Nanobiosci. 2015;14:385–92.

- [24]↑
Gaffar SA, Prasad VR, Baeg OA. Computational analysis of magnetohydrodynamic MHD free convection flow and heat transfer of non-Newtonian tangent hyperbolic fluid from a horizontal circular cylinder with partial slip. Int J Appl Comput Math. 2015;1:651–75.

- [25]↑
Nagendramma V, Leelarathnam A, Raju K, Shehzad A, Hussain T. Doubly stratified MHD tangent hyperbolic nanofluid flow due to permeable stretched cylinder. Results Phys. 2018;9:23–32.

- [26]↑
Das M, Mahanta G, Shaw S, Parida SB. Unsteady MHD chemically reactive double-diffusive Casson fluid past a flat plate in porous medium with heat and mass transfer. Heat transfer – Asian Res. 2019;48:1761–77.

- [27]↑
Sravanthi CS, Gorla RSR. Effects of heat source/sink and chemical reaction on MHD Maxwell nanofluid flow over a convectively heated exponentially stretching sheet using homotopy analysis method. Int J Appl Mech Eng. 2018;23:137–59.

- [28]↑
Hussain S, Mehmood K, Sagheer M, Yamin M. Numerical simulation of double diffusive mixed convective nanofluid flow and entropy generation in a square porous enclosure. Int J Heat Mass Trans. 2018;122:1283–97.

- [29]↑
Nield DA, Kuznetsov AV. The Cheng-Minkowycz problem for the double-diffusive natural convective boundary layer flow in a porous medium saturated by a nanofluid. Int J Heat Mass Trans. 2011;54:374–8.

- [30]↑
Mahapatra R, Pal D, Mondal S. Effects of buoyancy ratio on double-diffusive natural convection in a lid-driven cavity. Int J Heat Mass Trans. 2013;57:771–85.

- [31]↑
Gireesha BJ, Archana M, Prasannakumara BC, Gorla RSR, Makinde OD. MHD three dimensional double diffusive flow of Casson nanofluid with buoyancy forces and nonlinear thermal radiation over a stretching surface. Int J Num Meth Heat Fluid Flow. 2018;125:290–309.

- [32]↑
Rana GC, Chand R. Stability analysis of double-diffusive convection of Rivlin-Ericksen elastico-viscous nanofluid saturating a porous medium: a revised model. Forsch Ingenieurwes. 2015;79:1–2.

- [33]↑
Gaikwad SN, Malashetty MS, Prasad KR. An analytical study of linear and non-linear double diffusive convection with Soret and Dufour effects in couple stress fluid. Int J Non-Linear Mech. 2007;42:903–13.

- [34]↑
Kumar GK, Gireesha BJ, Manjunatha S, Rudraswamy NG. Effect of nonlinear thermal radiation on double-diffusive mixed convection boundary layer flow of viscoelastic nanofluid over a stretching sheet. Int J Mech Mater Eng. 2017;12:18.

- [35]↑
Ibrahim W, Gamachu D. Finite element method solution of mixed convection flow of Williamson nanofluid past a radially stretching sheet. Heat Transfer-Asian Res. 2019;12:1–23.

- [36]↑
Shateyi S, Marewo GT. Numerical solution of mixed convection flow of an MHD Jeffery fluid over an exponentially stretching sheet in the presence of thermal radiation and chemical reaction. Open Phys. 2018;16:249–59.

- [37]↑
Nalinakshi N, Dinesh PA, Chandrashekar DV. Shooting method to study mixed convection past a vertical heated plate with variable fluid properties and internal heat generation, Mapana. J Sci. 2017;13:31–50.

- [38]↑
El-Aziz MA, Tamer Nabil T. Homotopy analysis solution of hydromagnetic mixed convection flow past an exponentially stretching sheet with Hall Current. Math Prob Engg. 2012;2012:26.

- [39]↑
Beg OA, Khan MS, Ifsana Karim I. Explicit numerical study of unsteady hydromagnetic mixed convective nanofluid flow from an exponentially stretching sheet in porous media. Appl Nanosci. 2014;4:943–57.

- [40]↑
Zhang T, Salama A, Sun S, Zhong H. A compact numerical implementation for solving Stokes equations using matrix-vector operations. Procedia Comp Sci. 2015;51:1208–18.

- [41]↑
Pal D, Chatterjee S. Convective-radiative double-diffusion heat transfer in power-law fluid due to a stretching sheet embedded in non-Darcy porous media with Soret–Dufour effects. Int J Comput Methods Eng Sci Mech 2019;20:269–82.

- [42]↑
Keller HB. Numerical methods for two-point boundary value problems. New York: Dover publications; 1992.

- [43]↑
Cebeci T, Bradshaw P. Physical and computational aspects of convective heat transfer. New York: Springer Verlag; 1988.

- [44]↑
Pal D, Mondal SK. Magneto-bioconvection of Powell Eyring nanofluid over a permeable vertical stretching sheet due to gyrotactic microorganisms in the presence of nonlinear thermal radiation and Joule heating. Int J Ambient Energy. 2019;7:3723–31.

- [45]↑
Khan NS. Bioconvection in second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms, Brazilian. J Phys. 2019;48:227–41.

- [46]↑
Zaman S, Gul M. Magnetohydrodynamic bioconvective flow of Williamson nanofluid containing gyrotactic microorganisms subjected to thermal radiation and Newtonian conditions. Int J Ambient Energy. 2019;479:22–8.

- [47]↑
Wubshet I. MHD flow of a tangent hyperbolic fluid with nanoparticles past a stretching sheet with second order slip and convective boundary condition. Results Phys. 2017;7:3723–31.

- [48]↑
Xua H, Pop I. Mixed convection flow of a nanofluid over a stretching surface with uniform free stream in the presence of both nanoparticles and gyrotactic microorganisms. Int J Heat Mass Trans. 2014;75:610–63.

- [49]↑
Reddy BSK, Krishna MV, Surya KV, Rao N, Vijaya B. HAM Solutions on MHD flow of Nano-fluid through Saturated Porous medium with Hall effects. Mater Today. 2018;5:120–31.