In this paper, heat and mass transportation flow of swimming gyrotactic microorganisms (microbes) and solid nanoparticles under the viscous dissipation effect is investigated. The flow model PDEs are renovated with ordinary ones using suitable boundary layer approximations. The system governing the flow model dimensionless equations as well as boundary conditions is numerically treated with the SOR (successive over relaxation) technique. The flow, heat, and mass transport characteristics are examined against the prime parameters. A comparison is examined to be in a good agreement with the earlier results. It is found here that flow and thermal characteristics of the problem are substantially affected by the porous medium. The outcomes evidently point out that porous medium causes an enhancement in the skin friction and density of the motile microbes. Further, the rate of heat transport is devaluated at the surface of sheet due to viscous dissipation and elevated due to suction.

1. Introduction

An enhancement in heat transfer due to nanofluids is essentially required in various thermal systems. Nanofluids possess better thermal conductivity as equated to usual fluids. Several researchers interpreted the nanofluid flows under different assumptions. The flow phenomenon and thermal variations of unsteady radiative nanofluid flow induced by the magnetized spinning disk was investigated by Acharya et al. [1]. The presence of partial slips at the disk surface has been executed in this analysis. Salahuddin et al. [2] examined 3D stagnation point flow of hybrid nanofluid passed along a stretchable heated wavy cylinder under the impact of variable thickness and slip conditions. Shafiq et al. [3] determined the heat, mass, and motile microorganisms transfer rates in the magnetohydrodynamic (MHD) convective flow of a tangent hyperbolic material. A boundary layer flow of an incompressible Williamson fluid over a permeable radiative stretched surface was investigated by Shafiq and Sindhu [4]. Both electric and magnetic fields were taken into account. Imran et al. [5] interpreted the effects of melting phenomena and nonlinear thermal radiation in Cross-nanofluid bioconvection flow with motile microorganisms with a convective boundary over a cylinder. Brownian motion and thermophoresis diffusion were also taken into account in this mathematical model. Bilal et al. [6] presented a numerical simulation of microorganisms, carbon nanotubes (CNTs), and ferric oxide water-based hybrid nanofluid flow induced by a wavy fluctuating spinning disc with energy propagation. The nanofluid was synthesized in the presence of CNTs and magnetic nanoparticulates.

Recently, the intention of research community has shifted towards the study of gyrotactic microorganisms’ in nanofluid flows. A microorganism is a living being that can evolve, replicate, respond to stimuli, and have a structured order. In this way, any protists, animals, fungi, and bacteria will be an entity. These species can be categorized in various ways. Unicellular and multicellular animals are two essential groups. The microorganism is a much tiny living being that cannot be observed without an unaided eye. Algae, archaea, dust mite, and bacteria are all the examples of microorganisms. Taxes (locomotive behavior of microbes) may describe the microorganism’s movement. Some microbes are electrotaxis, rheotaxis, aerotaxis, chemotaxis, phototaxis, and gyrotaxis. The swimming microbs show fixed movement because of the torque. The torque produced from viscous drag and gravity in a flow causes the motion of gyrotactic microorganisms. The necessary factors, oxygen for breathing, light for photosynthetic, and search of nutrients (food) force microorganisms to move.

Enhanced nanofluid constancy, microscale fraternization, increased mass transmission, and application in microvolumes are the benefits that are obtained by incorporating the mobile microbes into the nanoparticle suspension [7]. Nanofluids with microbes can be used in many biomicrosystems such as nanochips to determine the toxicity and cellulose optimization. It is also beneficial to use the nanofluids with microbes in enzyme biosensors and microfluidics such as micro volumes and micromixers powdered. Our findings may also increase microbial fuel cell performance. Furthermore, gyrotcatic microbes can enhance the stability of nanofluids in a flow [8]. A microbial-enhanced oil restoration often utilizes bioconvection mechanisms in order to retain the shift in permeability, as nutrients and organisms are poured in the oil-carrying layer [9].

In nanofluid flows, the influences of gyrotactic organisms have been deeply studied. The transmission of heat and mass in the flow of gyrotactic organisms and nanoparticles was scrutinized by Ramzan et al. [10]. This research contained viscous desperation with nonlinear thermal radiation, slip effects, and convective boundary values. As the values of Peclet and Lewis number raised, a decreasing trend in the density of motile microbes was observed. In an investigation on the flow of nanofluids consisting of gyrotactic microorganisms, Chakraborty et al. [11] noticed a decrease in the thermal transport volume, boundary layers thickness, and flux of motile microorganisms through the fluid medium by increasing the magnetic field factor whereas the absorption of nanoparticles improved. Tausif et al. [12] explained the flow of nanofluid including motile microbes and nanoparticles with numerous slip influences. To gain the essential results in polymer industry, the values of numerous slip parameters were escalated, due to which the properties of fluid like microorganisms’ mass and flux rate or heat transference rate can be declined/improvised. Iqbal et al. [13] investigated the two-dimensional steady flow of nanofluid containing nanoparticles and gyrotactic microorganisms. The stagnation point flow of discursively striking nanofluid was included in this research. Atif et al. [14] explored that the density of motile microorganisms which was increased by the buoyancy ratio parameters but depreciated by the micropolar parameters. They integrated the influences of Joule heating and thermal radiation. The increasing effect of the density of microorganisms due to Peclet and Lewis numbers were investigated by Acharya et al. [15]. The numerical solutions were examined by Keller box method which marked the decreasing effect of the concentration of microbes due to elevation in values of the Peclet number. Waqas et al. [16] discussed some important properties such as Arrhenius energy activation, nonlinear radiation, and gyrotactic microorganisms in the flow of electrically conducted Maxwell nanofluid. In this study, bioconvection across the porous media due to microbes has also been discussed in detail.

Zuhra et al. [17] developed the mathematical model in which two nanoliquids (Williamson and Casson) were used in the flow and heat transfer within porous media upon its vertical solid surface against autocatalysis chemical reaction and gyrotactic microbes. In order to simplify the governing equations, nondimensional variables were used and the coupled ODEs were numerically solved using the homotopic analytic method (HAM). This study clarified how microorganism’s density was reduced by the parameters (e.g., heterogeneous chemical reaction and porous media). Aziz et al. [18] have performed numerical work on the boundary layer flow of the bioconvective nanofluid over a horizontal plate embedded in porous media and claimed that the bioconvection parameters had a marginal effect on the movement of mobile microorganisms. Shaw et al. [19] employed a spectral relaxation technique to solve nonlinear differential equations describing the flow of gyrotactic motile organisms and nanofluid. This bioconvection flow was considered over a sphere plunged within a permeable medium (under the combined influence of MHD and viscous dissipation). Kuznetsov [20] described the uncertain behavior of nonoscillatory and oscillatory conditions for nanofluid flow containing oxytactic organisms, which led to density stratification because of oxytactic motile organisms or the distribution of nanoparticles. Sarkar et al. [21] erected a bioconvective fluid in a permeable medium which holds gyrotactic microorganisms over an expanding surface. Aurangzaib et al. [22] oriented, for Darcy-Brinkman porous medium, a detailed study of bioconvection nanofluid together with gyrotactic microorganisms over a convectively heated wedge. Kessler [23] described, almost for the first time, the concept of bioconvection through a porous medium.

One of the aims to present this study is to determine how much bioconvection due to gyrotactic microorganisms affect the flow through a permeable media in the coexistence of viscous dissipation. Even though literature is available on the work under consideration but still bioconvective flow due to gyrotactic microbes and nanofluids under the influence of viscous dissipation through permeable media, specifically, over a nonlinear expanding/shrinking surface is not explored. The novel results of the present work, acquired by varying the preeminent parameters, are elaborated and illustrated with the help of graphs and tables.

2. Description of Physical Model

A steady flow of gyrotactic microorganisms and nanoparticles is considered within a permeable media. Viscous dissipation effect has also been taken in the energy equation. The velocity and the swimming path way of organisms are not influenced by the nanoparticles. However, volume fraction of nanoparticles, which is larger than 1%, can interrupt the motion of microorganisms. To attain the desired stability of bioconvection, a mixture of solid nanoparticles and microbes is added in the base liquid. An ambient concentration of microbes as well as ambient concentration and temperature of the fluid far away from the surface of sheet is represented by , , and (see Figure 1).

The relevant equations of the problem, in case of porous media and viscous dissipation, may be composed as [2427]

Here, -axis is taken normal to the sheet having velocity component , and -axis is considered to be in the direction of flow having velocity component . Moreover, the terms involved in the flow model equations such as , , , , , , , , , , , , , , and represent the thermophoretic coefficient, thermal conductivity, diffusivity of microorganisms, ambient concentration, cell swimming speed, temperature of fluid, nanofluid heat capacity ratio, Brownian diffusion coefficient, Darcy permeability, specific heat constant, kinematic viscosity, microorganism concentration, concentration of fluid, ambient temperature, and density, respectively. The BCs (boundary conditions) at sheet surface and far away from sheet can be written as

where is the mass flux which describes suction (if ) and injection (if ). Whereas epitomizes ambient velocity and signifies stretching velocity, here, is fixed () which expresses the nonlinearity of shrinking/stretching sheet. Nondimensional variables are utilized (see Aini Mat et al. [28]) to convert partial differential equations (1)–(5) into ordinary differential equations:

By using the coordinates of equation (8) in equations (2)–(5), the following nondimensional equations are achieved:

The analogous boundary conditions (6) and (7) take the form:

where denotes the injection/suction, is the stretching (if ) or contracting (if ) of surface parameter. The other prevailing parameters involved in equations (9)–(12) are as follows:

where is the motile microbes’ parameter, is the Lewis number, is the Eckert number, is the thermophoresis parameter, is the Peclet number, is the porosity parameter, is the bioconvection Schmidt number, is the Brownian motion, and is the Prandtl number

Further, dimensionless coordinates of engineering interest like skin friction, local density number of motile microbes, and Nusselt and Sherwood numbers may be defined as where local Reynolds number is denoted as .

3. Numerical Solution

Now, we employ SOR method to construct an algorithm [2933]. The order of equation (9) is reduced by one after modifying

Now, equations (9) (12) have been transformed as with boundary conditions:

Using finite differences, equations (17)–(20) take the following form:


The succeeding approximation process is realistic to attain the numerical solution of the problem: (a)To satisfy the BCs given in equation (21), the initial guesses are chosen for , and (b)The Simpson’s rule is applied to obtain from equation (16) which is used later in equation (22) to obtain (c)Finite difference discretization is utilized to simplify the system of equations (17)–(20), and then, the SOR method is employed to solve the resultant system of equations (22)–(25).(d)After the approximations , and are obtained, the process is repetitive until the functions , and converge(e)When the following inequality is fulfilled, then, the procedure of calculations is stopped:

We have taken value of as in all the calculations performed here.

3.1. Comparisons and Code Validation

The simulation results are associated with the previously accomplished results (Table 1). The present results are in good contact with the previous ones. This evaluation approves the correctness of our code. Table 1 also designates that the local Sherwood number increases by escalating suction parameter s.

The numerical data is selected from the available literature to associate our graphical outcomes with the existing ones (see Figure 2). Results are compared in the absence of porosity as well as bioconvection parameters. The comparison is in good contact with the literature results. The impact of Brownian motion parameter is to promote the temperature in the flow regime.

4. Results and Discussions

In this section, numerical results are discussed with the help of physical interpretations. The system of equations (17)–(20) is a numerically tackled subject to boundary conditions (21) for distinct parameters. The SOR method is utilized to attain the numerical solutions. The motile microbe’s density distribution , temperature , velocity , and concentration exhibit asymptotic expression for the estimated step size . Mass and heat transport attributes can be analyzed by selecting several estimation values to the prime parameters, instead of utilizing unique field measurements and fluid characteristics. In this way, the parametric values should be modified for the particular applications of the work (Ahmad et al. [3638]). Thus, we choose several values of the preeminent parameters such as Le, and otherwise specified.

Table 2 depicts that the density of local motile microbes and the skin friction are enhanced as the values of constant m increased positively; also, the same behavior is noticed for and .

Table 3 elucidates the effect of dimensionless porosity parameter on the physical quantities. The rates of heat and mass transfer are not significantly affected by the porosity of the medium which is due to the fact that the concentration and energy equation do not involve the porosity parameter. Moreover, only substantial change is noticed in shear stress. It may be due to the presence of (porosity parameter) in dimensionless momentum equation. As the pore size of the medium is larger than size of the microorganisms, so microbes’ movement is not influenced by the flow through porous media. The structure of the gyrotactic microorganisms is not influenced by the permeable media having much smaller permeability (or high porosity). The dilute suspension, for which bioconvection takes place, may disturb the desired concentration in the presence of porous media having low porosity (or very small pores). It is noticed here that an increment in the porosity of the medium provides better solutions for bioconvection parameters as compared to smaller porosity. It is worth mentioning here that the density of microorganisms is more stable when porosity is large.

Now, if we increase the values of Eckert number (e.g., ) then no change in skin friction is observed. It is noticed here that an increment in leads to a significant decrease in heat transport rate, and it tends to increase the mass transport rate as well as density of the motile microbes (see Table 4). However, viscous dissipation does not affect the bioconvection through porous medium, and the density enhances with viscous dissipation effect.

The rates of heat transportation for various values of , , and are provided in Table 5. This table demonstrates the fact that the thermophoresis and viscous dissipation effects tend to decrease the rate of heat transportation whereas the Prandtl number increases it. The suction of fluids and nanoparticles leads the fluid towards the stretching surface which causes an increase in the friction. However, with an increase in and , the rates of mass transport increase as predicted in Table 6.

It is noticed that the Brownian motion raises the rate of mass transfer whereas the thermophoretic parameter is a decreasing function for the rate of heat transfer. The bulk amount of heat is generated inside the fluid due to the collision of the nanoparticles. Because of this fluid and the presence of base fluid, the particles start motion in a crisscross manner. The thermophoretic parameter generates large amount of heat due to particles collision which enlarges the temperature but declares an opposite behavior for heat transportation rate.

The density of the microorganisms on the surface of the sheet increases because of bioconvection which collect the fluid alongside the sheet. Table 7 elucidates the improvement in the density of the microbes on the surface of the sheet when the values of , , and increase. The transportation of mass toward fluid from sheet occurs if the fluid concentration at the surface is less than the concentration of nanoparticles. The low concentration on the sheet causes due to the improvement in . It is observed that the shear stress increases due to the porosity parameter whereas the motile microbes density raises due to the bioconvection parameters (e.g., , , and ) correspondingly.

The different properties of describe the various cases such as shrinking and stretching of sheet is denoted by and , respectively, and there is no boundary layer when . The flow of planar stagnation point is associated with . The behavior of temperature, velocity, and microbes density profiles for distinct and is depicted in Figures 35. In both the cases, shrinking and stretching of sheet, respectively, velocity of fluid increases. On the other hand, and are depressed with the effect of positive constant . Figure 6 elucidates the consequences of injection/suction parameter and porosity parameter on streamwise velocity . The velocity elevates with an increase in and .

Figure 7 evidently reveals that the Prandtl number reduces the temperature in the flow regime. The profile of temperature is shown in Figure 8 for different values of . As we raise the values of both these parameters the temperature profile displays an increasing trend, the more heat transports in the fluid as the values of the Eckert number increase. This phenomenon affects the boundary layer thickness and causes an elevation in the temperature.

Thermophoretic and Brownian motions are accountable to the transport of nanoparticles. On the opposite, movement of mobile microbes powered by bioconvection is self-propelled. Therefore, the motion of tiny microbes is autonomous. The single-celled plants (phytoplankton) and single-celled green algae (chlamydomonas) can be analyzed as exemplary gyrotactic microbes. These are microscopic photosynthetic and green species that are useful to investigate the movement of microorganisms in fluid flows. These species can typically be contained in water or marine environment. The chlamydomonas and phytoplankton move upward due to gravitational force.

The effect of Lewis and Schmidt numbers for the bioconvection of motile organisms can be seen in Figure 9. A declining behavior in concentration profile is associated with the enlarging values of Le and Sc. The density of motile microbes depreciates while the viscous diffusion rate escalates due to increase in the bioconvection Lewis number.

Figure 10 clarifies that if we go on increasing the estimations of the microorganism’s parameter and the bioconvective Peclet number; then, the profiles are demoted. The bioconvection exists due to the up swimming of microorganisms on the upper surface of the fluid. In this problem, the viscous drag and the gravitational torques cause the motion of microorganisms. An increment in the values of the Peclet number upsurges the rate of advective transport associated to the dispersion and, on the other side, the flux of self-swimming motile microbes increases rapidly. A deceleration in the thickness of the motile organisms near the surface occurs due to the production of the rate of swimming motile microbes by the bioconvection Peclet number in the fluid. The bioconvection parameters substantially affect the density of the motile microbes. Thermal or temperature inversion occurs when a layer of hot fluid settles over a layer of cooler fluid that lies near the boundary. Viscous dissipation effect causes thermal inversion (reversal) due to which more heat is generated near boundary. This is the reason that some graphs deviate (little bit) from the imposed boundary conditions. Figures 11 and 12 are portrayed against concentration to analyze the effects of Suction parameter and Lewis number. It is found here that the the effect of both these parameters, together with Brownian motion parameter and positive constant, is to reduce the concentration.

5. Conclusions

This paper discloses the novel study regarding the bioconvective flow of gyrotactic microorganisms and nanofluids through a porous media. The flow is taken over a nonlinear expanding/shrinking surface. Bioconvection flows of microbes together with nanoparticles have much useful employments in several fields of science and biotechnology. The system of flow model equations is initially simplified by the similarity transformations, and then, successive over relaxation scheme is used to solve the resultant system numerically. From the outcomes of this examination, it is interpreted that the usage of nanoparticles together with gyrotactic microorganisms drives the flow towards consistency and stability. The major outcomes are listed below: (i)An increment in the motile microbes’ parameter and bioconvection Schmidt number leads towards an enhancement in the diffusion rate of microorganisms. On the other hand, a decreasing trend in the density profiles of mobile microorganisms is observed due to an increase in bioconvection parameters(ii)The viscous dissipation tends to decelerate the rate of heat transportation on the surface of sheet and elevates the temperature (iii)The porosity parameter negligibly affects the density of motile microbes in the absence of viscous dissipation (), while it substantially affects the density by increasing viscous dissipation(iv)The dimensionless stream wise velocity profiles seem to be magnifying with an increase in the values of porosity, injection/suction, and stretching/shrinking parameters


:Ambient temperature
:Fluid temperature
:Darcy permeability
:Microorganism concentration
:Cell swimming speed
:Stream function
:Brownian motion diffusivity
:Diffusivity of microorganisms
:Heat capacity of nanofluid
:Sherwood number
:Fluid density
:Dynamic viscosity
:Thermophoretic coefficient
:Rate constant of chemical reaction
:Similarity variable
:Dimensionless temperature
:Density of motile microorganisms
:Stretching velocity
:Nusselt number
:Ambient concentration
:Positive constant for nonlinear stretching sheet
:Thermal conductivity
:Ambient velocity
:Ambient microorganisms concentration
:Fuid concentration
:Skin friction
:Kinematics viscosity

Data Availability

The data used to obtain and finding this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors carried out the proofs and conceived of the study. All authors read and approved final form of the manuscript.