Abstract

This work addresses 3D bioconvective viscoelastic nanofluid flow across a heated Riga surface with nonlinear radiation, swimming microorganisms, and nanoparticles. The nanoparticles are tested with zero (passive) and nonzero (active) mass flux states along with the effect of thermophoresis and Brownian motion. The physical system is visualized via high linearity PDE systems and nondimensionalized to high linearity ordinary differential systems. The converted ordinary differential systems are solved with the aid of the homotopy analytic method (HAM). Several valuable and appropriate characteristics of related profiles are presented graphically and discussed in detail. Results of interest such as the modified Hartmann number, mixed convection parameter, bioconvection Rayleigh number, and Brownian motion parameter are discussed in terms of various profiles. The numerical coding is validated with earlier reports, and excellent agreement is observed. The microorganisms are utilized to improve the thermal conductivity of nanofluid, and this mechanism has more utilization in the oil refinery process.

1. Introduction

“Bioconvection” is known to be the convective movement within sight of swimming microorganisms. In this convective mode, the cells with bottom-heavy cells tend to swim at an angle to vertical, and this process is known as gyrotactic [1]. Therefore, the gyrotactic microorganisms are stable in the upper layer of the fluid, and consequently, stratification of the top-heavy fluid layer will become unbalanced. Thus, the system, which consists of a gyrotactic microorganism, induces one of the exciting characters in heat transfer that is “stability.” The reason is that nanofluids that have higher stability tend to improve the thermal efficiency of the heat exchanger (any energy systems). Hosseinzadeh et al. [2] examined the gyrotactic microorganism influence over a cylindrical surface with cross fluid flow. Mogharrebi et al. [3] present the MHD nanofluid flow towards a rotating cone with motile oxytactic microorganisms. Nowadays, the research on nanofluid through a Riga plate becomes an exciting area of research. For instance, mixed convective nanofluid flowed a Riga plate is studied numerically and analytically in [4]. It is shown that the desired size of the nanoparticle influences the skin friction coefficient. Ahmad et al. [5] studied the role of nanofluid past a heated vertical Riga plate numerically. Influence on viscous dissipation and thermal radiation of nanofluid flow between the Riga plate is explored numerically [6]. They have used carbon nanotubes as the nanoparticle, and it is shown that by varying the radiation parameter, the local heat transfer rate elevates. The application of the Cattaneo-Christov approach heat generation and absorption for the second-grade fluid that passed through the Riga plate is presented numerically [7]. In line with the application of the Riga plate, the stagnation flow on the vacillating Riga plate is studied numerically [8]. The variable thicked Riga plate for melting heat transfer application is explored numerically, and it is reported that for higher values of modified Hartmann number, the velocity profile distribution increases [9]. Recently, several research papers have been devoted to the study of nanofluid and their applications in practical situations [1013].

Further, the study on heat transfer characteristics of nanofluid is widely studied. For instance, viscoelastic fluid with Newtonian heating is numerically studied [1417]. It is shown that Al2O3/water nanofluid exhibits higher performance evaluation criteria. Besides, the study of viscoelastic nanofluid through the Riga plate is extensively studied. A nonuniform heat flux unsteady viscoelastic fluid that is unsteady is numerically investigated [18]. The viscoelastic nanofluid over an unsteady surface that is stretchable is evaluated numerically [19]. Also, the magnetic field effect on the Maxwell viscoelastic nanofluid over a plate that is moving at a uniform velocity is numerically assessed [20]. The temperature and velocity relaxation time influencing the heat transfer rate of the nanofluid are concluded by them. Besides that, the nonlinear effects on the MHD stagnation flow of viscoelastic nanofluid are explored numerically [21]. The Homotopy Analysis Method (HAM) is employed by Hayat et al. [22] to solve the 3D flow of a viscoelastic nanofluid over a stretching surface. The same research group extended their work towards the viscoelastic model for various applications [2325]. A slew of researchers studies the gyrotactic microorganism impact. For example, mixed convection of nanofluid containing third-grade nanomaterial containing gyrotactic microorganisms is figured out numerically [26]. Walters B nanofluid with the incorporation of gyrotactic microorganisms is evaluated numerically [27]. Acharya et al. [28] reported the effects of solar radiation bioconvection nanofluid with gyrotactic microorganisms. The suspension of microorganisms induced with the effect of the magnetic field is reported [29, 30]. Also, numerous researches on gyrotactic microorganisms are studied and explored [3035]. Analysis of active and passive controls with the chemical reaction of nanofluid is analyzed [36]. Using the homotopy perturbation method, unsteady nanofluid with active and passive controls is developed by Acharya et al. [37]. Over a bent surface, the active-passive control of dihydrogen monoxide nanofluid is explored numerically [38]. These controls are studied in various nanofluids with diverse applications [3944].

By considering the earlier reports, it is concluded that no studies were found to analyze the bioconvection of viscoelastic nanofluid on a 3D Riga surface with a comparison of active and passive control. The present study explores the 3D viscoelastic nanofluid flow across the heated Riga surface with nonlinear radiation and heat generation/absorption effects. The effects of several parameters are compared with the active and passive control model of nanoparticles. The homotopy analytic method [4551] is employed to study the present nonlinear ODE systems. The results are discussed in terms of various profiles. The numerical coding of the present study is validated with earlier reports. The relevant research is applied to multiple engineering streams like bioengineering, chemical, nuclear, thermal, and mechanical.

2. Problem Development

We have considered bioconvective viscoelastic nanofluid flow with in the direction in the direction over a 3D Riga surface with gyrotactic microorganisms swimming. The surface is expanding in all three directions, namely, , and It is assumed that nanoparticles do not influence the swimming microorganism; also, the nanoparticles are assumed to be stable in the fluid layer. Also, the nanoparticles have no effect on the velocity and temperature of the swimming microorganisms. The geometrical configuration is shown in Figure 1. By considering the above assumptions, the governing equations are described below [42]:

With boundary conditions

Now, we state the nondimensional similarity variables below:

Using the transformations in Equation (1), we get

Boundary condition (2) and (3) in expressions of is developed:

The dimensionless variables are

The nondimensional structure of surface drag force and heat transfer rate , mass transfer rate , and microorganism density number is stated as

3. Solution Approach: Homotopy Analysis Method (HAM)

The primary assumptions of the homotopy analytic method are stated as follows:

The auxiliary linear operators , , , , and are derived as

The above linear operators satisfying

The appropriate solutions are where denote the arbitrary conditions.

4. Convergence Analysis

To analyze the convergence of the numerical code, the convergence analysis is carried out for the value of The variable has control over the convergence and divergence of the numerical code, and it is shown in Figure 2. In the case of active and passive control, the values of in the range of . For the convergence study in the direction of and axes, the value of is in the range of -2.0 to 0.5, respectively. The order of approximation is shown in Tables 1 and 2. Table 3 reports that the value of is compared with earlier results with Qayyum et al. [14].

5. Results and Discussions

In this work, viscoelastic nanofluid past a stretching 3D Riga surface with gyrotactic microorganisms is explored numerically. The governing parameters such as velocity () parameters, temperature () parameter, , modified Hartmann number (), viscoelastic parameter (), mixed convection parameter (), buoyancy ratio parameter ), bioconvection Rayleigh number bioconvection Peclet number ), heat generation/absorption parameter (), thermophoresis parameter (), and microorganism concentration difference parameter () are studied.

Figure 3(a) describes the impact of (modified Hartmann’s number) on the velocity profile. It is observed that with an increase in , velocity profiles decrease in both directions (). This is due to the increase in leading to an increase in Lorentz force and consequently velocity profile decreases. In Figure 3(b), it is noticed that as the viscoelastic parameter () increases, the velocity profile decreases. In general, tensile stress is generated by viscoelasticity of the fluid. This stress opposes the fluid motion, and finally, the velocities of and directions are decays when the values of are improved. The increase in velocity profile is noted in Figure 3(c) with an enhancement in the mixed convection parameter (). By definition, is the ratio between the buoyancy force and the viscous force. Figure 3(d) reveals that an increase in buoyancy ratio parameter velocity profile decreases. In this work, the thermal and concentration forces are considered, and those forces provide smaller resistance consequently due to this reason that the velocity profile decreases. Also, it is observed in Figure 3(e) that the extending values of bioconvection Rayleigh number have a tendency to diminish the velocity profile. In Figure 3(f), it is noted that an increase in the bioconvection Peclet number produces the diminishing performance in the swimming speed of microorganisms that cause a decreasing trend in the velocity profile.

From Figure 4(a), it is evident that the increase in temperature profile increases with an increase in bioconvection Rayleigh number (). In fact, for larger values of , boundary layer thickness increases which increases the buoyancy forces; subsequently, temperature profile increases for both the cases of active and passive controls. The impact of temperature profile with regard to the thermophoresis parameter is discussed in Figure 4(b). With an increase in the thermophoresis parameter, temperature profile is increased. With an increase in , more nanoparticles shifted to the colder place from the hotter one, so the temperature profile increases. An increase in temperature ratio parameter increases the temperature profile as observed in Figure 4(c). The ratio of temperature at the surface () to the temperature at free stream () is mathematically defined as . For nonlinear radiation, the value of must be higher than 1. Moreover, an increasing trend found in the temperature along the surface is noted for larger values. Consequently, thermal boundary and associated layer thickness improve. The effect of the heat generation/absorption parameter, on the temperature profile, is proclaimed in Figure 4(d). As we increase , more movement between the nanoparticles is reported; therefore, the temperature profile increases. The slowing moment of the liquid flow due to the modified Hartmann number leads to producing a heat energy, and it is highly pronounced in the lesser value of modified Hartmann number , as seen in Figure 4(e).

Figure 5 illustrates the role of nanoparticle volume fraction on (a) thermophoresis parameter () and (b) Brownian motion parameter (), for active and passive control cases. It is noted that for both behaviours of active and passive control, the value increases. As we increase the value, more resistive force is applied to the surface of the nanoparticle; thus, to increase the heat transfer, the nanoparticle volume fraction has to be increased. This increasing trend is reported in Figure 5(a) for both the cases of active and passive controls. As overturn to the above discussion, it is noted in Figure 5(b) that the volume fraction of nanoparticles decreases with an increase in the Brownian motion parameter, . The fact is that the larger values of the Brownian motion parameter increase the length of the mean free path of the nanoparticle, which incredibly decreases the volume fraction of the nanoparticles.

The influence of the bioconvection Rayleigh number on the concentration of microorganisms is shown in Figure 6(a). With the increase in the Rayleigh number, the movement of microorganisms is increased; consequently, the concentration of microorganisms also increases. The role of the microorganism’s concentration difference parameter in concentration is explored in Figure 6(b). For both the cases of active and passive controls, the trend is reported to be the same. In Figure 6(c), it is seen that the concentration of microorganisms increases as the bioconvection Lewis number decreases. It is due to the variation in the temperature difference between the nanoparticles. It is noted in Figure 6(d) that the trend of active and passive control is in an opposite direction for the variation in the bioconvection Peclet number, . As we increase the value of Pe, the diffusivity of microorganisms is decreased. Therefore, the concentration of microorganisms is decreased for passive control. Whereas the scenario is in active control, the trend is increased. The thermophoresis parameter effect on concentration of microorganisms is plotted in Figure 6(e). With the increase of the low heated nanoparticle is moved into the high heated surface which is induced by the increase in . Hence, the concentration on microorganism is increased. In final, the variation in mixed convection parameter is analyzed in Figure 6(f). The increase in and the concentration on microorganisms is decreased. While for larger values of , temperature variation between the nanoparticles increases; therefore, this strongly affects the concentration of microorganisms.

Figures 7 and 8 depict the heat transfer rate () for different combinations of pertinent parameters. From Figure 7, the heat transfer rate enhances for higher values of modified Hartmann number (), and it decays for radiation parameter ). Figure 8 elucidates that heat transfer rates reduce for larger values of and . Figure 9 explores the combination of and in the mass transfer rate (). It is noted from this figure that the mass transfer rate enhances for the Brownian motion parameter for active control while an opposite behaviour is noticed for the thermophoresis parameter for both cases. Figures 10 and 11 portray the motile microorganism’s density number rate () for different related parameters. From these figures, we found that the combined parameters and show the inverse effects on the microorganism’s density number rate for active and passive controls.

6. Conclusion

In this work, the active and passive control of viscoelastic nanofluid flow over a 3D Riga plate is investigated analytically by the homotopy technique. Effects are studied with the influence of the gyrotactic microorganism in bioconvective heat transfer. The following outcomes are observed: (i)An increase in the modified Hartmann number produces the diminishing impact on velocity profile for both and directions(ii)Fluid velocity decays in both and directions for the augmentation in bioconvection Rayleigh number(iii)The microorganism profile is an increasing function of bioconvection Peclet number for active control, and it is decreasing for passive control(iv)Mass transfer rate enhances for Brownian motion parameter while an opposite behaviour is noticed for thermophoresis parameter for both cases(v)An increase in bioconvection Peclet number produces the inverse phenomena in microorganism profile for active and passive controls

Nomenclature

:Positive constants (-)
:Concentration
:Ambient concentration
:Surface concentration of nanoparticles
:Specific heat
:Skin friction in direction
:Skin friction in direction
:Chemotaxis constant (m)
:Brownian diffusion coefficient (m2 s-1)
:Thermophoretic diffusion coefficient (m2 s-1)
:Microorganism’s diffusion coefficient (m2 s-1)
:Velocity similarity function in direction (—)
:Velocity similarity function in direction
:Heat generation/absorption parameter
:Thermal conductivity
:Lewis number
:Bioconvection Lewis number
:Brownian motion parameter
:Thermophoresis parameter
:Surface concentration of microorganisms
:Ambient concentration of microorganisms
:Buoyancy ratio parameter
:Nusselt number
:Microorganisms’ density number
:Prandtl number
:Bioconvection Peclet number
:Modified Hartmann number
:Dimensional heat generation/absorption coefficient
:Bioconvection Rayleigh number
Rd:Radiation parameter
:Sherwood number
:Temperature
:Ambient temperature
:Velocity of the sheet
:Velocity components
:Maximum cell swimming speed
:Cartesian coordinates (m).
Greeks
:Material parameter of fluid
:Volume expansion coefficient
:Viscoelastic parameter
:Concentration similarity function
:Stretching ratio
:Similarity parameter
:Microorganisms’ similarity function
:Mixed convection parameter
:Kinematic viscosity
:Ratio of the effective heat capacity
:Temperature similarity function
:Density
:Density of nanofluid
:Density of nanoparticles
:Density of microorganism’s particles
:Electrical conductivity
:Stream function
:Microorganisms’ concentration difference parameter .

Data Availability

The raw data supporting the conclusions of this article will be made available by the corresponding author without undue reservation.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.