Next Article in Journal
Magnetron Sputtering Thin Films as Tool to Detect Triclosan in Infant Formula Powder: Electronic Tongue Approach
Previous Article in Journal
Spectrally Selective Solar Absorber Coating of W/WAlSiN/SiON/SiO2 with Enhanced Absorption through Gradation of Optical Constants: Validation by Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Radiative Bioconvective Flow of a Sisko Nanofluid with Motile Microorganisms

1
School of Science, Huzhou University, Huzhou 313000, China
2
Department of Mathematics, Layyah Campus, Government College University Faisalabad, Layyah 31200, Pakistan
3
Department of Mathematics, COMSATS University Islamabad, Sahiwal 57000, Pakistan
4
Department of Mathematics and Statistics, Riphah International University I-14, Islamabad 44000, Pakistan
*
Author to whom correspondence should be addressed.
Coatings 2021, 11(3), 335; https://doi.org/10.3390/coatings11030335
Submission received: 4 February 2021 / Revised: 27 February 2021 / Accepted: 5 March 2021 / Published: 15 March 2021

Abstract

:
The progressive and enhanced thermal mechanisms of nanoparticles has motivated researchers to give attention to this topic in recent years. The synthesizing and versatile applications of such materials include cooling and heating controlling processes, solar systems, energy production, nanoelectronics, hybrid-powered motors, cancer treatments, and renewable energy systems. Moreover, the bioconvection of nanofluids allows for some motivating applications in this era of bioengineering and biotechnology, such as biofuels, biosensors, and enzymes. With these interesting motivations and applications, this study elucidated upon the three-dimensional bioconvection flow of a Sisko fluid (base fluid) in the presence of a nanofluid over a stretched surface. The additional thermal features of radiation were also incorporated to modify the analysis. The rheological features of shear thinning and shear thickening that are associated with the Sisko nanofluid were comprehensively studied. The problem was formulated using highly nonlinear and coupled differential equations, which were numerically simulated via a shooting scheme. The salient physical applications of flow parameters were graphically underlined in view of shear-thinning and shear-thickening scenarios. The results showed that a decrease in velocity in the presence of buoyancy ratio forces was more conducive to the shear-thinning phenomenon. The increase in temperature profile due the thermal Biot number and surface heating source parameter seemed to be more inflated in the shear-thinning scenario. A lower motile microorganism profile was noted for the bioconvection Lewis number.

1. Introduction

The enhancement in the thermal efficiency of innumerable materials has drawn scientists’ attention due to these materials’ applications in engineering and industrial development. The improvement in materials with a low heat transportation rate (ethylene glycol, water) is necessary for engineering and industrial processes, where base liquids gain importance as an energy source. Nanofluids are used to enhance the proficiency of ordinary convectional heat transfer fluids, which are consistent mixtures of a stable fluid. In vehicles, nanofluids are also used as a coolant to improve the efficiency of the heat exchanger/aeration and cooling systems. Nanofluids are widely used in automobile and industrial cooling, cancer treatments, identification, new sources of petrol, microelectronics cooling, cross-breeding motors, home heating/cooling equipment, and so on. Choi [1] used the term nanomaterials for the first time, bearing in mind the overall aim of increasing the heat conductivity of nanofluids. Buongiorno [2] established a theoretical non-homogeneous equilibrium theory for the molecular diffusion of nanofluids. Thermophoresis and Brownian motion are considered important tools for evaluating the slip of nanoparticles. Iqbal et al. [3] examined the impact of convection on a Sisko fluid with peristaltic movement in an asymmetrical tube. Ali et al. [4] studied the physical importance of chemical reactions and the Lorentzian intensity aspects of Sisko fluid flows in a curved structure. Chaudhuri and Rathore [5] reviewed an observational analysis of the pressure-driven flow and the temperature distribution of a Sisko fluid flowing across parallel plates with viscous dissipation. Nabil and Mostapha [6] investigated the Hall current and Joule heating effects of a peristaltic flow of a Sisko fluid with moderate stenosis via a porous medium in the tapered artery with slip and convective process parameters. Toghraie et al. [7] analyzed blood flow within various arteries using a non-Newtonian Sisko model for use in bioengineering. Agoor et al. [8] investigated peristaltic flow with heat transfer to a Sisko fluid in ciliated arteries. Ferdows et al. [9] examined the interaction of two-dimensional microorganisms with the movement of mixed convection nanoparticles past elongated sheets. Gangadhar et al. [10] studied aspects of an unsteady free convective nanofluid boundary layer flow over a stretch sheet using a spectral relaxation technique. Jawad et al. [11] analyzed the entropy production and mass generation in the unsteady spinning magnetohydrodynamics (MHD) of aqueous nanotube structures with thermal radiation and entropy generation impacts. Ali et al. [12] investigated the core features of the mechanical field of a nanofluid expansion cylinder. Nisar et al. [13] examined the physical elaboration of an Eyring–Powell nanofluid under the effect of activation energy. Muhammad et al. [14] researched the significance of mixed convection in a 3D Eyring–Powell nanoliquid flow using activation energy. Muhammad et al. [15] studied the importance of melting temperature variations in a pressurized base fluid (water), non-liquid (carbon nanotubes (CNTs)+water), and hybrid nanofluid (CNTs + CuO + water) flows.
The analysis of non-Newtonian materials is another captivating research area because of their multidisciplinary characteristics and complicated nature. Many biological liquids (including blood), polymer solutions, liquid crystals, solutions, and lubricating oils provide examples of non-Newtonian behavior. The rheological features of such materials cannot be characterized by employing the Newtonian model. A literature survey indicated that various types of non-Newtonian models have been suggested by researchers in order to specify the complicated nature of such materials. The Sisko fluid model is the most modified model, which involves a combination of viscous and power-law models. The Sisko fluid model successfully predicts the interesting features of shear thickening and shear thinning. Various physical properties of viscous and viscoelastic materials are effectively predicted by adjusting the various parameters that are used in this model [16,17,18,19,20].
The bioconvection phenomenon is based on the convection of various materials at the microscopic level. The process of bioconvection is dynamically associated with the fluctuation in the density gradient of the swimming motile microorganisms. It is commonly noticed that microorganisms swim near the upper portion of the surface. The phenomenon of bioconvection is considered an impulsive motion of unicellular microbes in liquids, which results in a nonuniform distribution because of the associated densities of the fluid molecules. The occurrence of these microorganisms increases the primary density of materials and subsequently instigates a gradient with a floating density, which leads to bioconvection. The major cause of bioconvection is the gradient density of motile microorganisms. This theory has significant potential for applicability. For example, it has been used in biology and medicine for cancer treatments or in microsystems for enzyme biomaterials. In addition, it has recently entered the oil and gas industry, in particular, as a technique for the production of microbial fuel. The bioconvection flow of nanoparticles is important for biomolecules, building structures, nanofluid processing, transportation processes, etc. Wager [21] and Platt [22] first introduced the idea of bioconvection. Gyrotactic motile microorganisms often swim in the opposite direction to gravity, where the essential fluid is water. The principle of nanoliquid bioconvection was first investigated by Kuznetsov [23]. Zohra et al. [24] analyzed the slip flow of a nanofluid in the presence of motile microorganisms that were confined by a rotating disk. Zhang et al. [25] assessed the motion of axially symmetrical stagnation point motions of magnetization rate nanomaterials that were designed with a porous extended rotating disc in the absence of motile gyrotactic microorganisms. Kotha et al. [26] inspected the two-dimensional hydromagnetic flow and heat transport processes of water-based nanofluids, such as gyrotactic motile microorganisms across a vertical surface through heat generation and absorption. An inviscid, electrically conductive, bioconvective peristaltic fluid flow between two expandable discs was documented by Siddiq et al. [27]. Kumaraswamy et al. [28] measured the effects of heat radiation and partial slipping on the magnetohydrodynamic movement of Jeffrey nanoparticles consisting of motile microorganisms through a vertically spread substrate. Abdelmalek et al. [29] investigated various Arrhenius slip flow relations of third-grade nanofluid over a stretching configuration. Waqas et al. [30] recently demonstrated the value of the bioconvection method through the use of nanomaterials and produced significant scientific and industrial applications. Waqas et al. [31] investigated Darcy–Forchheimer nanofluid bioconvection flowing through the expansion of a cylinder/surface with modified thermal and mass transfer, activation energy, and gyrotactic microorganism characteristics. Wang et al. [32] investigated the active Prandtl features of a bioconvective heat-formed magnetizing tangent hyperbolic nanofluid with gyrotactic motile microorganisms, as well as second-order velocity slipping. Several researchers [33,34,35,36] have recently introduced a bioconvective swimming microorganism model to investigate the characteristics of heat transformation.
This paper provides a continuation of the aforementioned research and presents the three-dimensional analysis of the steady and incompressible flow assessment of a Sisko nanofluid accommodating motile microorganisms. The nonlinear thermal radiation relations were imposed to get improve the thermal characteristics of nanofluids. The convective boundary conditions were used to perform the analysis. The novel features of a non-Newtonian Sisko fluid, namely, shear thinning and shear thickening, were observed for all flow parameters. The benchmark conversion of the problem into a dimensionless form by using nonlinear equations, which were numerically solved by using a shooting technique. It is remarked that some theoretical developments already exist in the literature that provide simulations for 3D Sisko nanofluid flows. However, the applications of the bioconvection phenomenon and utilization of nonlinear thermal radiation features have not been reported yet. Therefore, this communication aimed to fill this research gap. A comprehensive physical quantitative analysis was performed for each flow parameter.

2. Mathematical Modeling

We assumed a steady flow assessment of a Sisko nanofluid confined by a bidirectional stretched surface. The nanofluid contained motile microorganisms with a uniform density. In an electrically conducting Sisko nanofluid, the magnetic force induction was incorporated along the z-direction. The magnitude of the magnetic force was not specified due to the small magnetic Reynolds number assumption. The Brownian and thermophoretic properties of the nanofluid were examined via the Buongiorno model. The mathematical model was developed by using a Cartesian coordinate system, where the velocity components u and v were assumed along the stretching surface, while the z axis was assumed to be perpendicular to this surface. The surface stretching velocity was denoted with U w ( x ) , while the V w ( y ) was the stretching velocity along the normal axis, as shown in Figure 1.
The salient features of the thermal radiation were also considered. The governing flow equations of the current flow problem are as follows [16,17]:
u x + v y + w z = 0 ,
u u x + v v y + w w z = a ρ f 2 u z 2 b ρ f z ( u z ) n σ B 0 2 ρ f u + 1 ρ f [ ( 1 C ) ρ f β * * g * ( T T ) ( ρ p ρ f ) g * ( C C ) ( N N ) g * γ ( ρ m ρ f ) ] ,
u v x + v v y + w v z = a ρ f 2 v z 2 b ρ f z ( u z ) n 1 v z σ B 0 2 ρ f v ,
u T x + v T y + w T z = α 1 2 T z 2 + τ ( D B T z C z + D T T ( T z ) 2 ) 16 σ * 3 ( ρ c ) f k * z ( T 3 T z ) ,
u C x + v C y + w C z = D B ( 2 C z 2 ) + D T T ( T z ) 2 ,
u N x + v N y + w N z + [ z ( N C z ) ] b W c ( C w C ) = D m z ( N z ) ,
with the following boundary assumptions [16,17]:
u = U w ( x ) = c x , v = V w ( y ) = d y , w = 0 , k T z = h f ( T w T ) , D B C z + D T T T z = 0 , N = N w a t z = 0 ,
u 0 , v 0 , T T , C C , N N a s z .
where ( a , b ) are the Sisko fluid materials constants, ρ f is the fluid density, n is the behavior index, σ 1 is the electrical conductivity, B 0 is the magnetic field strength, β * * is the volume expansion coefficient, g is the acceleration due to gravity, ν is the kinematic viscosity, T is the nanofluid temperature, ρ p is the density of nanoparticles, ρ m is the density of microorganisms, τ = ( ρ c ) p / ( ρ c ) f is the heat capacity ratio, C is the nanofluid concentration, D T is the thermophoresis diffusion coefficient, D B is the diffusion constant, D m is the microorganism diffusion constant, σ is the Stefan–Boltzmann constant, k denotes the mean absorption coefficient, and b ^ is the chemotaxis constant.
Suitable variables were used to transform the governing equations into dimensionless form, as follows [16,17]:
u = c x f ( ζ ) , v = d y g ( ζ ) , w = c ( c n 2 ρ f / b ) 1 n + 1 [ 2 n n + 1 f + 1 n 1 + n ζ f + g ] x n 1 n + 1 , θ ( ζ ) = T T T f T , φ ( ζ ) = C C C , χ ( ζ ) = N N N w N , ζ = z ( c 2 n b / ρ f ) x 1 n 1 + n . }
Applying these variables to Equations (1)–(8) yields [16,17]:
A f + 2 n n + 1 f f + n ( f ) n 1 f f 2 M 2 f + g f + λ ( θ + N r ϕ + N c χ ) = 0 ,
A g + g g + n ( f ) n 1 g + ( 2 n n + 1 ) f g M 2 g ( n 1 ) g f ( f ) n 2 g 2 = 0 ,
[ { 1 + R d ( 1 + ( θ w 1 ) θ ) 3 } θ ] + P r [ N t θ 2 + 2 n n + 1 f θ + N b θ φ + g θ ] = 0 ,
φ + L e P r [ ( 2 n n + 1 ) f φ + g φ ] + N t N b θ = 0 ,
χ + L b [ ( 2 n n + 1 ) f χ + g χ ] P e [ χ φ + ( Ω + χ ) φ ] = 0 ,
with the following related boundary conditions [16,17]:
f ( 0 ) = 0 , g ( 0 ) = 0 , f ( 0 ) = 1 , g ( 0 ) = α , θ ( 0 ) = B i ( 1 θ ( 0 ) ) , N b ϕ ( 0 ) + N t θ ( 0 ) = 0 , χ ( 0 ) = 1 ,
f ( ) 0 , g ( ) 0 , θ ( ) 0 , φ ( ) 0 , χ ( ) 0 ,
with the material parameter of the Sisko fluid A , the Prandtl number P r , the buoyancy ratio parameter N r , the mixed convection parameter λ , the bioconvection Rayleigh number N c , the magnetic parameter M , the velocity ratio parameter R d , the radiation parameter R d , the temperature ratio parameter θ w , the Brownian motion parameter N b , the thermophoresis parameter N t , the Lewis number L e , the Peclet number P e , the bioconvection Lewis number L b , and the Biot number B i , which are defined as follows:
A = R e b 2 n + 1 R e a , R e a = U w x ρ f a , R e b = U w n 2 x n ρ f b , P r = x U α 1 R e b 2 n + 1 , N r ( = ( ρ p ρ f ) ( C w C ) ρ f ( 1 C ) ( T w T ) β * * ) , λ ( = β * * g * ( 1 C ) ( T w T ) a U w ) , N c ( = γ * ( ρ m ρ f ) ( N w N ) ρ f ( 1 C ) ( T w T ) β * * ) , M = σ B 0 2 ( ρ c ) f , α = d c , R d = 16 σ * T 3 3 k * k , θ w = T f T , B i = h f k x Re b 1 n + 1 , N b = τ D B C ν , N t = τ D B ( T f T ) T ν , L e = α 1 D B , P e ( = b W c D m ) , L b ( = ν D m ) , B i ( = h f k ν a ) .
The relations for the skin friction coefficient, the local Nusselt number, the local Sherwood number, and the motile density number are presented as follows [16,17]:
R e b 1 n + 1 C f x = A f ( 0 ) ( f ( 0 ) ) n , R e b 1 n + 1 C f y = V w U w ( A g ( 0 ) + ( f ( 0 ) ) n 1 g ( 0 ) ) . }
R e b 1 n + 1 N u x = [ 1 + R d { 1 + ( θ w 1 ) θ ( 0 ) } 3 ] θ ( 0 ) .
R e b 1 n + 1 S h x = f ( 0 ) . R e b 1 n + 1 N n x = c ( 0 ) .
where R e a = U w x ρ f / a and R e b = U w 2 n x ρ f / b signify the local Reynolds numbers.

3. Numerical Solution

Owing to the highly nonlinear and coupled nature of the formulated equations, finding an exact equation was not possible. Therefore, the computational shooting algorithm was used in MATLAB (USA). The higher-order differential structure was discretized into first-order differential equations by exploiting appropriate variables. To this end, let us emphasize the following linear operators that were used:
f = s 1 , f = s 2 , f = s 3 , f = s 3 , g = s 4 , g = s 5 , g = s 6 , g = s 6 , θ = s 7 , θ = s 8 , θ = s 8 , φ = s 9 , φ = s 10 , φ = s 10 , χ = s 11 , χ = s 12 , χ = s 12 ,
s 3 = s 2 2 ( 2 n n + 1 s 1 s 3 ) s 4 s 3 + M 2 s 2 λ ( s 7 + N r s 9 + N c s 11 ) ( A + n ( s 3 ) n 1 ) ,
s 6 = ( n 1 ) s 6 s 3 ( s 3 ) n 2 ( 2 n n + 1 ) s 1 s 6 + s 5 2 s 4 s 6 + M 2 s 5 ( A + n ( s 3 ) n 1 ) ,
s 8 = Pr [ ( 2 n n + 1 ) s 1 s 8 + s 4 s 8 + N b s 8 s 10 + N t s 10 2 ] [ 1 + R d ( 1 + ( θ w 1 ) ) 3 ] s 8 3 ,
s 10 = L e Pr [ ( 2 n n + 1 ) s 1 s 10 + s 4 s 10 ] N t N b s 8 ,
s 12 = L b [ ( 2 n n + 1 ) s 1 s 12 + s 4 s 12 ] + P e [ s 12 s 10 + ( Ω + s 11 ) s 10 ] .
Similarly the boundary conditions were transformed as follows:
s 1 ( 0 ) = 0 , s 4 ( 0 ) = 0 , s 2 ( 0 ) = 1 , s 5 ( 0 ) = α , s 8 ( 0 ) = B i ( 1 s 7 ( 0 ) ) , N b s 10 + N t s 8 = 0 , χ ( 0 ) = 1 ,
s 2 0 , s 5 0 , s 7 0 , s 9 0 , s 11 0 a s ζ .
With an iterative process, the numerical simulations were repeated until an accuracy of 10−5 was achieved.

4. Solution Verification

The obtained numerical results were verified to confirm the solution’s accuracy. To this end, the present results were compared as a limiting case with the work of Hayat et al. [17] in Table 1. An excellent solution accuracy between both analyses was noted for n = 0.5 and n = 1.5.

5. Results and Discussion

The physical outcomes for the bioconvection flow of a Sisko nanofluid for the modeled parameters are conveyed in this section. Since it was already communicated that a Sisko fluid model captures the dynamics of the shear-thinning ( n > 1 ) and shear-thickening ( n < 1 ) features, all the graphical results were prepared with such consequences in mind. Figure 2 describes the nature of the mixed convection parameter λ for the velocity profile f when n = 0.5 and n = 1.5 . An increased velocity profile was observed for higher λ values. However, the increase in the velocity was comparatively higher for the shear-thinning case with n = 1.5 . The physical assessment of this flow dynamic was attributed to the fact that the mixed convection constant influenced the buoyancy forces. An increase in λ increased the buoyancy force, and as such, the velocity profile had an increased range. Figure 3 shows the evaluation of f for the bioconvection N c . It was found that an increasing behavior of the bioconvection Rayleigh number decreased the velocity distribution for both cases ( n = 0.5 and 1.5 ) . The effects of the buoyancy ratio parameter N r on the velocity profile f are depicted in Figure 4. From these curves, it was found that the buoyancy ratio parameter reduced the velocity concentration f for both shear thinning and thickening ( n = 0.5 and 1.5 ) . Figure 5 reveals the outcomes of the magnetic parameter M regarding the prediction of shear thickening and thinning aspects for the velocity f . The velocity f reduced with a higher magnetic parameter. Physically, the magnetic parameter influenced the Lorentz force, which controlled the fluid motion parameters in the flow system. The significances of the magnetic parameter M for velocity g is revealed in Figure 6. It was perceived that the velocity of fluid g reduced with the magnetic parameter. For the shear-thickening behavior, the reduction in the velocity was relatively slow. The influence of θ w on the temperature profile θ for both cases is presented in Figure 7. Again, the results were prepared using n = 0.5 and n = 1.5 . It was found that increasing values of the temperature ratio parameter boosted the nanofluid temperature. The thermal radiation was more effective at enhancing the nanofluid temperature. The variation of the Prandtl number P r on the temperature profile θ is presented in Figure 8. It was found that higher values of the Prandtl number decreased the temperature of the fluid θ for both values of the behavior index n . Since the Prandtl number is the ratio of the momentum diffusivity to the thermal diffusivity, an increase in the Prandtl number enhanced the momentum diffusivity but slowed the thermal diffusivity, which caused a reduction in the nanofluid temperature. Figure 9 presents the results for the Biot number B i on the temperature θ . The temperature distribution θ increased the Biot number, which was more comprehensive for shear-thinning case. The physical aspects of such an enhanced graphical trend were due to the fact that the Biot number is related to the heat transfer coefficient. The outcome of the Brownian motion constant N b on ϕ , subject to the shear thinning n = 1.5 and thickening n = 0.5 perspectives, is reflected in Figure 10. Increasing the values of the Brownian motion parameter decreased the concentration of the nanoparticles. The decreasing association of ϕ with N b was expected because it replicated the opposite relation with the nondimensional ϕ . To scrutinize the domination of the thermophoresis parameter N t on the volumetric concentration field ϕ , Figure 11 is presented. It was found that increasing the magnitude of the thermophoresis parameter increased the volumetric concentration of the nanoparticles ϕ . Physically, the thermophoresis phenomenon describes the heated nanoparticles that wandered to the low-temperature space. Figure 12 presents the outcomes of the Lewis number L e on ϕ in view of the shear thinning and shear thickening parameters. Increasing the Lewis number decreased the concentration distribution ϕ . Physically, the decreasing trend in ϕ was deemed to be the reverse relation with mass diffusivity. To disclose the increase in the Peclet number P e on the microorganisms field χ , Figure 13 is presented. The summarized graphical trend shows that χ decreased as P e increased. The reverse effect of the Peclet number on the diffusivity of the motile density caused such a trend. In order to analyze the behavior of χ due to L b , Figure 14 was plotted. An abated profile of χ was noted due to the increased values of L b .
Table 2 presents the impact of the local skin friction coefficients along x-axis f ( 0 ) with respect to various flow parameters when n = 0.5 and n = 1.5. It was noted that an increase in the wall shear force acting along the x-direction produced increases in N r , M and N c , which were more prominent when n = 1.5 . The change in g ( 0 ) due to the flow parameters are presented in Table 3. An increasing change in g ( 0 ) was observed for increases in N r and N c , while a contrasting trend was observed for an increase in λ . Table 4 depicts the salient features of the Nusselt number θ ( 0 ) for different flow parameters. This analysis was performed with the presence and absence of Brownian motion N b and the thermophoresis parameter N t . It is clearly highlighted that the increasing variation in θ ( 0 ) due to increases in λ , P r , and R d was more pronounced in the presence of N b and N t . Table 5 presents the numerical data for the local Sherwood number ϕ ( 0 ) . With increases in P r and N b , the local Sherwood number increased. Table 6 highlights the physical significance of the local microorganism number in terms of various parameters. It was clearly noticed that the motile density number was at a maximum when P e and L b were at their maximum values.

6. Conclusions

The three-dimensional thermally radiative flow of a Sisko nanofluid fluid containing microorganisms over a stretched surface was examined. Quantitative analyses for all physical flow parameters in view of the shear-thinning and shear-thickening aspects was comprehensively discussed. Numerical simulations via the shooting technique were performed. The Brownian motion and thermophoresis features were investigated by incorporating the Buongiorno model. The main outcomes are listed as follows:
An increase in the velocity due to the mixed convection parameter was observed, which was enhanced in the shear-thinning case.
The hindering effects on the velocity due to the bioconvection Rayleigh number was relatively slower when shear-thickening effects were dominant.
The nanofluid temperature was enhanced due to the Biot number and surface-heating source parameter. The enhancement in temperature in the case of shear thickening was qualitatively slower as compared to the shear-thinning case.
An increase in the nanofluid concentration due to the thermophoretic parameter was more progressive when shear-thinning features were considered.
The microorganism profile decreased with the bioconvection Lewis number and Peclet number.

Author Contributions

H.G.-J., presents the results analysis and physical explanation of problem; H.W., presented the flow model; S.U.K., computed the solution; M.I.K., made the literature survey; S.F., analyzed the results; S.H., write the discussion of results. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u , v velocity components (m·s−1) x , y coordinate axes (m)
T temperature(K) C nanoparticles concentration (kg·m−3)
ν kinematic viscosity (m2·s−1) T w wall temperature (K)
T free stream temperature (K) C ambient concentration (kg·m−3)
ρ f density (kg·m−3) D B Brownian diffusion coefficient (m2·s−1)
D T thermophoresis diffusion coefficient (m2·s−1) μ dynamic viscosity (N·s·m−2)
g gravity τ heat capacity ratio
k mean absorption coefficient D m microorganism diffusion constant
σ Stefan–Boltzmann constant q r radiative flux
b ^ chemotaxis constant γ volume suspension coefficient
( a , b ) Sisko fluid materials constants σ 1 electrical conductivity
B 0 magnetic field strength β * * volume expansion coefficient
ρ p density of nanoparticles ρ m density of microorganisms
A Sisko fluid material parameter P r Prandtl number
N r buoyancy ratio parameter λ mixed convection parameter
N c bioconvection Rayleigh number M magnetic parameter
R d radiation parameter θ w temperature ratio parameter
N t thermophoresis parameter N b Brownian motion parameter
L e Lewis number P e Peclet number
L b bioconvection Lewis number B i Biot number
R e local Reynolds number N u x local Nusselt number
S h x local Sherwood number N n x motile density number

References

  1. Choi, S.U.S.; Eastman, J. Enhancing thermal conductivity of fluids with nanoparticles. ASME 2001, 231, 718–720. [Google Scholar]
  2. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  3. Iqbal, N.; Yasmin, H.; Kometa, B.K.; Attiya, A.A. Effects of convection on Sisko fluid with peristalsis in an asymmetric channel. Math. Comput. Appl. 2020, 25, 52. [Google Scholar] [CrossRef]
  4. Ali, M.; Irfan, M.; Khan, W.A.; Sultan, F.; Shahzad, M.; Khan, M. Physical significance of chemical processes and Lorentz’s forces aspects on Sisko fluid flow in curved configuration. Soft Comput. 2020, 24, 16213–16223. [Google Scholar] [CrossRef]
  5. Chaudhuri, S.; Rathore, S.K. An analytical investigation of pressure-driven flow and heat transfer of a Sisko fluid flowing through parallel plates with viscous dissipation. Sadhana 2020, 45, 1–17. [Google Scholar] [CrossRef]
  6. Nabil, T.; Mostapha, D.R. Hall current and joule heating effects on peristaltic flow of a Sisko fluid with mild stenosis through a porous medium in a tapered artery with slip and convective boundary conditions. Sains Malays. 2020, 49, 1175–1190. [Google Scholar]
  7. Toghraie, D.; Esfahani, N.N.; Zarringhalam, M.; Shirani, N.; Rostami, S. Blood flow analysis inside different arteries using non-Newtonian Sisko model for application in biomedical engineering. Comput. Methods Progr. Biomed. 2020, 190, 105338. [Google Scholar] [CrossRef] [PubMed]
  8. Agoor, B.M.; Sayed-Ahmed, M.E.; Alam, H. Peristaltic flow with heat transfer on Sisko fluid in a ciliated arteries. Int. J. Fluid Mech. Therm. Sci. 2020, 6, 70–78. [Google Scholar] [CrossRef]
  9. Ferdows, M.; Reddy, M.G.; Sun, S.; Alzahrani, F. Two-dimensional gyrotactic microorganisms flow of hydromagnetic power-law nanofluid past an elongated sheet. Adv. Mech. Eng. 2019, 11. [Google Scholar] [CrossRef]
  10. Gangadhar, K.; Kannan, T.; Sakthivel, G.; Dasaradha Ramaiah, K. Unsteady free convective boundary layer flow of a nanofluid past a stretching surface using a spectral relaxation method. Int. J. Ambient Energy 2020, 41, 609–616. [Google Scholar] [CrossRef]
  11. Jawad, M.; Shah, Z.; Khan, A.; Khan, W.; Kumam, P.; Islam, S. Entropy generation and heat transfer analysis in MHD unsteady rotating flow for aqueous suspensions of carbon nanotubes with nonlinear thermal radiation and viscous dissipation effect. Entropy 2019, 21, 492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ali, M.; Sultan, F.; Khan, W.A.; Shahzad, M.; Arif, H. Important features of expanding/contracting cylinder for cross magneto-nanofluid flow. Chaos Solitons Fractals 2020, 133, 109656. [Google Scholar] [CrossRef]
  13. Nisar, Z.; Hayat, T.; Alsaedi, A.; Ahmad, B. Significance of activation energy in radiative peristaltic transport of Eyring-Powell nanofluid. Int. Commun. Heat Mass Transf. 2020, 116, 104655. [Google Scholar] [CrossRef]
  14. Muhammad, T.; Waqas, H.; Khan, S.A.; Ellahi, R.; Sait, S.M. Significance of nonlinear thermal radiation in 3D Eyring–Powell nanofluid flow with Arrhenius activation energy. J. Therm. Anal. Calorim. 2020, 143, 1–16. [Google Scholar] [CrossRef]
  15. Muhammad, K.; Hayat, T.; Alsaedi, A.; Ahmad, B. Melting heat transfer in squeezing flow of base fluid (water), nanofluid (CNTs + water), and hybrid nanofluid (CNTs + CuO + water). J. Therm. Anal. Calorim. 2020, 1–18. [Google Scholar] [CrossRef]
  16. Hayat, T.; Taseer Muhammad, B.; Ahmad, S.; Shehzad, A. Impact of magnetic field in three-dimensional flow of Sisko nanofluid with convective condition. J. Magn. Magn. Mater. 2016, 413, 1–8. [Google Scholar] [CrossRef]
  17. Hayat, T.; Taseer Muhammad, B.; Ahmad, S.; Shehzad, A.; Alsaedi, A. On three-dimensional boundary layer flow of Sisko nanofluid with magnetic field effects. Adv. Powder Technol. 2016, 27, 504–512. [Google Scholar] [CrossRef]
  18. Pal, D.; Mandal, G. Magnetohydrodynamic stagnation-point flow of Sisko nanofluid over a stretching sheet with suction. Propuls. Power Res. 2020, 9, 408–422. [Google Scholar] [CrossRef]
  19. Khan, M.I. Transportation of hybrid nanoparticles in forced convective Darcy-Forchheimer flow by a rotating disk. Int. Commun. Heat Mass Transf. 2021, 122, 105177. [Google Scholar] [CrossRef]
  20. Chu, Y.M.; Khan, M.I.; Khan, N.B.; Kadry, S.; Khan, S.U.; Tlili, I.; Nayak, M.K. Significance of activation energy, bio-convection and magnetohydrodynamic in flow of third grade fluid (non-Newtonian) towards stretched surface: A Buongiorno model analysis. Int. Commun. Heat Mass Transf. 2020, 118, 104893. [Google Scholar] [CrossRef]
  21. Wager, H. On the effect of gravity upon the movements and aggregation of Euglena Viridis, Ehrb., and other micro-organisms. Philos. Trans. R. Soc. Lond. B 1911, 201, 333–390. [Google Scholar]
  22. Platt, J.R. “Bioconvection pattern” in cultures of the free-swimming organism. Science 1961, 133, 1766–1767. [Google Scholar] [CrossRef] [PubMed]
  23. Kuznetsov, A. Bio-thermal convection induced by two different species of microorganisms. Int. Commun. Heat Mass Transf. 2011, 38, 548–553. [Google Scholar] [CrossRef]
  24. Zohra, F.T.; Uddin, M.J.; Basir, F.; Ismail, A.I.M. Magnetohydrodynamic bio-nano-convective slip flow with Stefan blowing effects over a rotating disc. Proc. Inst. Mech. Eng. Part N J. Nanomater. Nanoeng. Nanosyst. 2019, 234, 83–97. [Google Scholar] [CrossRef]
  25. Zhang, T.; Khan, S.U.; Imran, M.; Tlili, I.; Waqas, H.; Ali, N. Activation energy and thermal radiation aspects in bioconvection flow of rate-type nanoparticles configured by a stretching/shrinking disk. J. Energy Resour. Technol. 2020, 142, 1–19. [Google Scholar] [CrossRef]
  26. Kotha, G.; Kolipaula, V.R.; Rao, M.V.S.; Penki, S.; Chamkha, A.J. Internal heat generation on bioconvection of an MHD nanofluid flow due to gyrotactic microorganisms. Eur. Phys. J. Plus 2020, 135, 1–19. [Google Scholar] [CrossRef]
  27. Siddiq, M.K.; Ashraf, M. Bioconvection of micropolar nanofluid with modified Cattaneo–Christov theories. Adv. Mech. Eng. 2020, 12. [Google Scholar] [CrossRef]
  28. Kumaraswamy Naidu, K.; Harish Babu, D.; Harinath Reddy, S.; Satya Narayana, P.V. Radiation and partial slip effects on MHD jeffrey nanofluid containing gyrotactic microorganisms over a stretching surface. J. Therm. Sci. Eng. Appl. 2020, 1–28. [Google Scholar] [CrossRef]
  29. Abdelmalek, Z.; Ullah Khan, S.; Waqas, H.; ANabwey, H.; Tlili, I. Utilization of second-order slip, activation energy, and viscous dissipation consequences in the thermally developed flow of third-grade nanofluid with gyrotactic microorganisms. Symmetry 2020, 12, 309. [Google Scholar] [CrossRef] [Green Version]
  30. Abdelmalek, Z.; Khan, S.U.; Waqas, H.; Riaz, A.; Khan, I.A.; Tlili, I. A mathematical model for bioconvection flow of Williamson nanofluid over a stretching cylinder featuring variable thermal conductivity, activation energy, and second-order slip. J. Therm. Anal. Calorim. 2020, 1–13. [Google Scholar] [CrossRef]
  31. Waqas, H.; Imran, M.; Muhammad, T.; Sait, S.M.; Ellahi, R. On bio-convection thermal radiation in Dar-cy–Forchheimer flow of nanofluid with gyrotactic motile microorganism under Wu’s slip over-stretching cylinder/plate. In-ternational Journal of Numerical Methods for Heat & Fluid Flow.tion and Brownian Motion. Mach. Learn. Res. 2020, 4, 51. [Google Scholar]
  32. Wang, Y.; Waqas, H.; Tahir, M.; Imran, M.; Jung, C.Y. Effective prandtl aspects on bio-convective thermally developed magnetized tangent hyperbolic nanoliquid with gyrotactic microorganisms and second-order velocity slip. IEEE Access 2019, 7, 130008–130023. [Google Scholar] [CrossRef]
  33. Li, Y.; Waqas, H.; Imran, M.; Farooq, U.; Mallawi, F.; Tlili, I. A Numerical exploration of modified second-grade nanofluid with motile microorganisms, thermal radiation, and Wu’s slip. Symmetry 2020, 12, 393. [Google Scholar] [CrossRef] [Green Version]
  34. Mahdy, A.; Nabwey, H.A. Microorganisms time-mixed convection nanofluid flow by the stagnation domain of an impulsively rotating sphere due to Newtonian heating. Results Phys. 2020, 103347. [Google Scholar] [CrossRef]
  35. Khan, S.U.; Waqas, H.; Muhammad, T.; Imran, M.; Ullah, M.Z. Significance of activation energy and Wu’s slip features in Cross nanofluid with motile microorganisms. Commun. Theor. Phys. 2020, 72, 105001. [Google Scholar] [CrossRef]
  36. Ansari, M.S.; Otegbeye, O.; Trivedi, M.; Goqo, S.P. Magnetohydrodynamic bio-convective casson nanofluid flow: A numerical simulation by paired quasilinearisation. J. Appl. Comput. Mech. 2020. [Google Scholar] [CrossRef]
Figure 1. Flow geometry.
Figure 1. Flow geometry.
Coatings 11 00335 g001
Figure 2. Change in f ( ζ ) with λ .
Figure 2. Change in f ( ζ ) with λ .
Coatings 11 00335 g002
Figure 3. Change in f ( ζ ) with N c .
Figure 3. Change in f ( ζ ) with N c .
Coatings 11 00335 g003
Figure 4. Change in f ( ζ ) with N r .
Figure 4. Change in f ( ζ ) with N r .
Coatings 11 00335 g004
Figure 5. Change in f ( ζ ) with M .
Figure 5. Change in f ( ζ ) with M .
Coatings 11 00335 g005
Figure 6. Change in g ( ζ ) with M .
Figure 6. Change in g ( ζ ) with M .
Coatings 11 00335 g006
Figure 7. Change in θ ( ζ ) with θ w .
Figure 7. Change in θ ( ζ ) with θ w .
Coatings 11 00335 g007
Figure 8. Change in θ ( ζ ) with P r .
Figure 8. Change in θ ( ζ ) with P r .
Coatings 11 00335 g008
Figure 9. Change in θ ( ζ ) with B i .
Figure 9. Change in θ ( ζ ) with B i .
Coatings 11 00335 g009
Figure 10. Change in ϕ ( ζ ) with N b .
Figure 10. Change in ϕ ( ζ ) with N b .
Coatings 11 00335 g010
Figure 11. Change in ϕ ( ζ ) with N t .
Figure 11. Change in ϕ ( ζ ) with N t .
Coatings 11 00335 g011
Figure 12. Change in ϕ ( ζ ) with L e .
Figure 12. Change in ϕ ( ζ ) with L e .
Coatings 11 00335 g012
Figure 13. Change in χ ( ζ ) with P e .
Figure 13. Change in χ ( ζ ) with P e .
Coatings 11 00335 g013
Figure 14. Change in χ ( ζ ) with L b .
Figure 14. Change in χ ( ζ ) with L b .
Coatings 11 00335 g014
Table 1. Comparative analysis for the obtained numerical solution with Hayat et al. [17] when M = 0.2 , α = 0.1 , λ = N r = N c = 0 .
Table 1. Comparative analysis for the obtained numerical solution with Hayat et al. [17] when M = 0.2 , α = 0.1 , λ = N r = N c = 0 .
A Hayat et al. [17]Present Results
R e b 1 n + 1 C f x R e b 1 n + 1 C f x
n = 1 n = 2 n = 1 n = 2
0.01.03941.00501.03951.0048
0.51.27301.20761.27311.2075
1.01.47001.47141.47001.4714
Table 2. Numerical amplitudes of the local skin friction coefficient with respect to various parameters.
Table 2. Numerical amplitudes of the local skin friction coefficient with respect to various parameters.
Parameters f ( 0 )
λ N r N c M n = 0.5 n = 1.5
0.2
0.4
0.8
0.10.10.41.0726
1.0709
1.0675
1.4402
1.4392
1.4370
0.10.8
1.6
2.2
1.2735
1.2746
1.2748
1.6531
1.6554
1.6573
0.8
1.6
2.2
1.2999
1.3308
1.3542
1.6796
1.7127
1.7380
0.2
0.4
0.8
1.1952
1.2731
1.4160
1.5934
1.6512
1.7614
Table 3. Numerical amplitudes of the local skin coefficient with respect to various parameters.
Table 3. Numerical amplitudes of the local skin coefficient with respect to various parameters.
Parameters g ( 0 )
λ N r N c M n = 0.5 n = 1.5
0.2
0.4
0.8
0.10.10.40.0911
0.0910
0.0908
0.1477
0.1446
0.1444
0.10.8
1.6
2.2
0.2983
0.2990
0.2997
0.4405
0.4499
0.4504
0.8
1.6
2.2
0.2978
0.2988
0.2996
0.4367
0.4435
0.4521
0.2
0.4
0.8
0.2816
0.2986
0.3633
0.4307
0.4410
0.4856
Table 4. Numerical amplitudes of the local Nusselt number with respect to various parameters.
Table 4. Numerical amplitudes of the local Nusselt number with respect to various parameters.
Parameters θ ( 0 ) θ ( 0 )
( N t = N b = 0.3 ) ( N t = N b = 0.0 )
λ N r N c M P r R d n = 0.5 n = 1.5 n = 0.5 n = 1.5
0.2
0.4
0.8
0.10.10.42.00.60.3885
0.3891
0.3932
0.5128
0.5234
0.5297
0.3371
0.3154
0.2904
0.4231
0.4346
0.4454
0.10.8
1.6
2.2
0.5360
0.5352
0.5346
0.4700
0.4686
0.4676
0.4267
0.4054
0.3909
0.4134
0.3889
0.3667
0.8
1.6
2.2
0.5347
0.5325
0.5307
0.4686
0.4655
0.4631
0.4378
0.4205
0.4056
0.4102
0.4067
0.3956
0.2
0.4
0.8
0.5443
0.5367
0.5199
0.4753
0.4711
0.4617
0.4479
0.4334
0.4056
0.4156
0.3924
0.3854
2.5
3.0
3.5
0.5921
0.6383
0.6775
0.5259
0.5733
0.6145
0.5303
0.5554
0.5934
0.4865
0.5276
0.5467
0.8
1.4
1.8
0.4917
0.4271
0.3963
0.4288
0.3713
0.3453
0.4165
0.3704
0.32565
0.3776
0.3467
0.3156
Table 5. Numerical amplitudes of the local Sherwood number with respect to various parameters.
Table 5. Numerical amplitudes of the local Sherwood number with respect to various parameters.
Parameters ϕ ( 0 )
λ N r N c M P r N t L e N b n = 0.5 n = 1.5
0.2
0.4
0.8
0.10.10.42.00.32.00.8145
0.8256
0.9310
0.7954
0.8167
0.8966
0.10.8
1.6
2.2
0.8040
0.8028
0.8019
0.7050
0.7030
0.7014
0.8
1.6
2.2
0.8021
0.7987
0.7960
0.7029
0.6983
0.6947
0.2
0.4
0.8
0.8165
0.8050
0.7799
0.7129
0.7067
0.6925
2.5
3.0
3.5
0.8882
0.9575
1.0162
0.7889
0.8600
0.9218
0.1
0.4
0.8
0.8180
0.8045
0.7956
0.7870
0.7756
0.7434
3.0
4.0
5.0
0.7991
0.7952
0.7922
0.7012
0.6976
0.6950
0.1
0.4
0.8
0.8990
1.0234
1.1014
0.7854
0.7901
0.9912
Table 6. Numerical amplitudes of the local microorganism number with respect to various parameters.
Table 6. Numerical amplitudes of the local microorganism number with respect to various parameters.
Parameters χ ( 0 )
λ N r N c M P e L b n = 0.5 n = 1.5
0.2
0.4
0.8
0.10.10.40.12.01.0853
1.1012
1.1145
0.9534
0.9611
0.9917
0.10.8
1.6
2.2
1.0751
1.0738
1.0728
0.9429
0.9405
0.9386
0.8
1.6
2.2
1.0724
1.0678
1.0642
0.9378
0.9336
0.9287
0.2
0.4
0.8
1.0916
1.0762
1.0422
0.9539
0.9450
0.9243
0.5
1.0
1.5
1.3580
1.7185
2.0874
1.2022
1.5295
1.8625
3.0
4.0
5.0
0.2375
0.2382
0.2388
0.2285
0.2290
0.2293
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ge-JiLe, H.; Waqas, H.; Khan, S.U.; Khan, M.I.; Farooq, S.; Hussain, S. Three-Dimensional Radiative Bioconvective Flow of a Sisko Nanofluid with Motile Microorganisms. Coatings 2021, 11, 335. https://doi.org/10.3390/coatings11030335

AMA Style

Ge-JiLe H, Waqas H, Khan SU, Khan MI, Farooq S, Hussain S. Three-Dimensional Radiative Bioconvective Flow of a Sisko Nanofluid with Motile Microorganisms. Coatings. 2021; 11(3):335. https://doi.org/10.3390/coatings11030335

Chicago/Turabian Style

Ge-JiLe, Hu, Hassan Waqas, Sami Ullah Khan, Muhammad Ijaz Khan, Shahid Farooq, and Sajjad Hussain. 2021. "Three-Dimensional Radiative Bioconvective Flow of a Sisko Nanofluid with Motile Microorganisms" Coatings 11, no. 3: 335. https://doi.org/10.3390/coatings11030335

APA Style

Ge-JiLe, H., Waqas, H., Khan, S. U., Khan, M. I., Farooq, S., & Hussain, S. (2021). Three-Dimensional Radiative Bioconvective Flow of a Sisko Nanofluid with Motile Microorganisms. Coatings, 11(3), 335. https://doi.org/10.3390/coatings11030335

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop