Next Article in Journal
Some Asymptotic Properties of a Kernel Bispectum Estimate with Different Multitapers
Next Article in Special Issue
Prediction and Optimization of Pile Bearing Capacity Considering Effects of Time
Previous Article in Journal
On Some Improved Class of Estimators by Using Stratified Ranked Set Sampling
Previous Article in Special Issue
Activation Energy Performance through Magnetized Hybrid Fe3O4PP Nanofluids Flow with Impact of the Cluster Interfacial Nanolayer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Power-Law Nanofluid Flow over a Stretchable Surface Due to Gyrotactic Microorganisms

1
Department of Mathematics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
2
Department of Basic Engineering Science, Faculty of Engineering, Menoufia University, Shebin El-Kom 32511, Egypt
3
Department of Mechanical Engineering, College of Engineering, Prince Mohammad Bin Fahd University, Al Khobar 31952, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Aswan University, Aswan 81528, Egypt
5
Department of Information Technology, Fanshawe College London, London, ON N5Y 5R6, Canada
6
Basic and Applied Sciences Department, College of Engineering and Technology, Arab Academy for Science & Technology and Maritime Transport (AASTMT), Aswan Branch 81528, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(18), 3285; https://doi.org/10.3390/math10183285
Submission received: 24 July 2022 / Revised: 30 August 2022 / Accepted: 7 September 2022 / Published: 9 September 2022
(This article belongs to the Special Issue Modeling and Simulation in Engineering, 2nd Edition)

Abstract

:
This study aims to learn more about how the flow of a power-law nanofluid’s mixed bio-convective stagnation point flow approaching a stretchable surface behaves with the presence of a passively controlled boundary condition. The governing equations incorporate the motile bacterium and nanoparticles, and the current model includes Brownian motion and thermophoresis effects. The governing equations are transformed into ordinary differential equations, which are then numerically solved using the Runge–KuttaFehlberg (RKF) with the shooting technique. The controlling parameters are chosen as follows: the velocity ratio parameter, ε, is taken between 0.1 and 1.5; the mixed convection parameter, λ, is considered in the range 0–3; the buoyancy ratio parameter is considered in the range between 0.1 and 4; the bio-convection parameter, Rb, is taken in the range 0–1; nanofluid parameters are taken in the range 0.1–0.7; the bioconvection Schmidt number is considered in the range 0.1–3; the Prandtl number is taken between 1–4; and the Schmidt number is taken between 1 and 3. The Nusselt number, skin friction, and nanoparticle volume fraction profiles are shown graphically to observe the impact of several parameters under consideration. Both the Schmidt number and the Brownian motion parameter are shown to significantly increase the Sherwood number. Thermophoresis, however, has been proven to lower the Sherwood number. Furthermore, the bioconvection constant and Peclet number both help to slow down the rate of mass transfer. The presented theoretical investigation has a considerable role in engineering, where nanofluid flow is applied to organize a bioconvection process to develop power generation and mechanical energy. One of the more essential features of bioconvection is the aggregation of nanoparticles with motile microorganisms requested to augment the stability, heat, and mass transmission.
MSC:
76D05; 76D10; 80A10; 80A19

1. Introduction

The importance of nanofluids has been increasing with time, and investigators intend to study the behavior of nanofluids subjected to heat transfer systems. Nanofluids and their implications in the industrial sector have grown more because of their homogeneous nature in thermal conductivity and rudimentary heat transfer. Typical fluids such as propylene glycol, water, and ethylene glycol, among others, have poor heat transfer properties. Nano-fluid, a homogenous solid–liquid mixture, is applied to enhance the thermal conductivity of the base fluid. Choi [1] and Das et al. [2] exposed the most remarkable characteristic of nanofluids: thermal conductivity is dependent on temperature. Nanofluids are extensively used in energy and biomedical applications (nano-drug delivery, cancer therapeutics, and nano cryosurgery) [3,4,5]. Buongiorno [6] examined the effect of convective transport in nanofluid by observing the heat transfer properties of the Brownian motion and thermophoresis. Numerous studies are presented in the literature to demonstrate the uses of nanofluid in various fields [7,8,9,10].
The study of heat relocation in power-law (non-Newtonian) fluids has received significant attention because of its application in different branches of modern technologies, industries, and engineering like polymer-thickened oils, liquid crystals, polymeric suspensions, and physiological fluid mechanics. Furthermore, instances showing a diversity of non-Newtonian characteristics contain biological fluids, pharmaceutical formulations, toiletries, synthetic lubricants, cosmetics, paints, and foodstuffs (see Irvine and Karni [11]). Various models have been reported to analyze the non-Newtonian attitude toward fluids. Among these patterns, which are famous for following the empirical Ostwald–de Waele model, in which the shear stress varies according to a power function of the strain rate, gained much acceptance. The theoretical analysis of power-law fluid was first scrutinized by Schowalter [12] and Acrivos et al. [13]. Rashad et al. [14] examined the power-law nanofluid flow across a vertical cone in a porous medium. Later on, several studies were analyzed by many investigators [5,15,16,17].
Bioconvection occurs as the natural microbe swims upwards. Thus, the microorganism is denser than the base fluids. Because there are so many microbes on the upper surface of the foundation, it becomes weak. As a result, the bacteria decrease and promote bioconvection, and the microbes’ return to swimming strengthens the process. This bacterial emigration into the water raises the temperature and mass transfer in the environment. Because of medical advancements, microscopic kinds have significantly contributed to the improvement in human life. Microorganisms are essential for life and it cannot exist without them. Decreasing the length of the cavities and the cell resistance enables the construction of continuum numerical patterns. It is often approved that the nanoparticles of concentration distribution are massive relative to the cell pivot. Bioconvection occurs as combined nanofluids are addressed using heat and mass conversion. Platt [18] further proposed the concept of bioconvection, which characterizes the micro-structural flow brought on by gradation density, in addition to motile gyro-tacticmicro-organisms. The first study of the bioconvection of gyro-tactic motile bacteria, including nanoparticles, was conducted by Kuznetsov and Avramenko [19]. Kuznetsov [20] presented key discoveries of nanofluidbioconvection in a horizontal porous layer, including both nanoparticles and gyro-tactic motile bacteria. Khan and coworkers [7,21,22,23,24] investigated the free convection of non-Newtonian nanofluid numerically down a vertical plate in a porous media. They considered that the medium contains gyrotactic microorganisms and that the temperature, nanoparticle concentration, and motile microbe density are all controlled in the plate.
The local Nusselt, Sherwood, and density numbers are found to be strongly influenced by nanofluid and bioconvection parameters. Additionally, they provided a mathematical model to examine the flow of a water-based nanofluid containing gyrotactic micro-organisms around a truncated cone with a convective boundary condition at the surface. It has been discovered that, as a surface grows rougher, the densities of the mobile microorganisms, Sherwood number, Nusselt number, and skin friction all increase. Nabwey and collaborators [8,9,10] investigated the flow of a nanofluid containing gyrotactic bacteria over an isothermal cone surface in the presence of viscous dissipation and Joule heating. Consideration was given to the combined effects of a transverse magnetic field and Navier slip in the flow. They also considered mixed bioconvective flow on a vertical wedge in a Darcy porous media filled with a nanofluid that contains both nanoparticles and gyrotactic bacteria. To combine energy and concentration equations with passively controlled boundary conditions, the effects of thermophoresis and Brownian motion are considered. They found that the buoyancy ratio and magnetic field parameters increase the local skin friction coefficient, Nusselt number, Sherwood number, and local density of the motile microorganism’s number. Ishak et al. [25] analyzed the flow of a two-dimensional stagnation point of an incompressible viscous fluid near a steady mixed convection boundary layer flow across a stretching vertical sheet in its plane. For an aiding flow, they demonstrated that both the skin friction coefficient and the local Nusselt number increase with the buoyancy ratio; however, when the Prandtl number increases, only the local Nusselt number increases and the skin friction coefficient decreases.
However, the current study intends to close the knowledge gap regarding the mixed bioconvective stagnation point flow of a power-law nanofluid towards a stretchable surface. The nanofluid flow will be mathematically modeled using Buongiorno’s two-component, including the Brownian movement and thermophoresis aspects. Moreover, modeling of the motile microorganism and nanoparticles is addressed in the governing equations.

2. Governing Equations

The present study considers the mixed bioconvection flow of a non-Newtonian power-law nanofluid having motile microorganisms and obeying the Ostwald–de Waele model [18] near the stagnation point at a heated stretchable vertical surface coinciding with the plane y = 0. The flow is confined to the region y > 0, where x and y are the cartesian coordinates. The model suggested by Buongiorno is utilized for the nanofluid attitude in which the influence of thermophoresis and Brownian movement is taken into consideration. The stretchable surface is maintained at a constant temperature Tw along with the constant density of motile microorganisms Nw. The ambient temperature, concentration, and motile microorganisms are symbolized as T∞, C∞, and N∞, respectively. It is considered that the stretching velocity is given by Uw(x) = cx and the velocity of external flow in the neighborhood of the stagnation point at x = y = 0 is given by U∞(x) = ax, where a and c are positive constants. Although, in the past, this reality is eliminated, in a practical situation, there is no nanoparticle flux on the boundaries, and the values of the nanofluid fraction C adapt to the concentration distribution. This means that we consider passively controlled boundary conditions as proposed by Kuznetsov and Nield [26].
Under the above-mentioned assumptions, the physical description of the problem under consideration in Figure 1 with Boussinesq approximations can be expressed as follows [8,24]:
u x + v y = 0
u u x + v u y = U d U d x + 1 ρ f τ x y y + 1 ρ f ( 1 C ) g β ( T T ) ( ρ p ρ f ) g ( C C ) ( ρ m ρ f ) g γ ( N N )
u T x + v T y = α 2 T y 2 + τ D B C y T y + D T T T y 2
u C x + v C y = D B 2 C y 2 + D T T 2 T y 2
u N x + v N y + b W c C y N C y = D m 2 N y 2
The boundary conditions can be set in the following form [14]:
u = U w ( x ) , v = 0 , T = T w , D B C y + D T T T y = 0 ,   N = N w at   y = 0
u = U ( x ) , T = T ,   C = C ,   N = N as   y
In the present problem, we have ∂u/∂y > 0 when a/c > 1 (the ratio of free stream velocity and stretching velocity), which gives the shear stress where τ x y = K u / y n . Here, K is the consistency coefficient and n is the power-law fluid. It is noted that, when n = 1, the fluid is Newtonian; when n < 1, the fluid is called pseudoplastic power-law fluid; and when n > 1, it is called dilatants power-law fluid.
Here ( u , v ) are the velocity components, along with the ( x , y ) directions. Here, T, C, and N are the fluid temperature, concentration, and density of motile microorganisms, respectively, and g is the gravitational acceleration. Wc denotes the maximum cell swimming speed, α stands for thermal diffusivity, β is the coefficient of thermal expansion, γ represents the average volume of a microorganism, ρf represents the density of the fluid, ρp denotes the density of the particles, ρm is the density of the microorganism, τ = ( ρ c ) p / ( ρ c ) f is the nanofluid heat capacity ratio, ρ c f is the heat capacity of the base fluid and ρ c p is the effective heat capacity of the nanoparticle material, Dn stands for diffusivity of the microorganisms, DB represents the Brownian diffusion coefficient, and DT stands for the thermophoretic diffusion coefficient of the microorganisms.
It is observed that the continuity equation is automatically determined by specifying the upgraded stream function, such that u = ψ / y , v = ψ / x , and exhibits the following non-dimensional variables:
θ = T T T w T , ϕ = C C C , χ = N N N w N ,   ψ = K / ρ f c 1 2 n 1 / ( n + 1 ) x 2 n / ( n + 1 ) F ( η ) ,   η = y c 2 n K / ρ f 1 / ( n + 1 ) x ( 1 n ) / ( 1 + n )
Equations (1)–(5) take the following non-dimensional form [8,9]:
n F ( n 1 ) F + 2 n 1 + n F F F 2 + ε 2 + λ θ N r ϕ R b   χ = 0
θ + Pr 2 n 1 + n F θ + N b θ ϕ + N t θ 2 = 0
ϕ + S c 2 n n + 1 F ϕ + N t N b θ = 0
χ + S b 2 n 1 + n F χ P e ϕ χ + ( χ + σ ) ϕ = 0
F ( 0 ) = 1 , F ( 0 ) = 0 , θ ( 0 ) = 1 , N b ϕ ( 0 ) + N t θ ( 0 ) = 0 , χ ( 0 ) = 1
F ( ) = ε , θ ( ) = 0 , ϕ ( ) = 0 , χ ( ) = 0
where the prime denotes differentiation with respect to η and λ = G r x Re x 2 is the dimensionless mixed convection parameter, N r = ( ρ p ρ f ) C ρ f β ( T w T ) ( 1 C ) stands for the buoyancy ratio parameter, G r x = ( 1 C ) ρ f g β ( T w T ) ( K / ρ f ) 2 is the local Grashof number, Re x = ( c x ) 2 n x n K / ρ f is the local Reynolds number, N b = τ D B C c x 2 Re x 2 / ( n + 1 ) is the Brownian motion parameter, N t = τ D T ( T w T ) T c x 2 Re x 2 / ( n + 1 ) is the thermophoresis parameter, Pr = c x 2 α Re x 2 / ( n + 1 ) is the Prandtl number and S c = c x 2 D B Re x 2 / ( n + 1 ) is the Schmidt number, R b = ( N w N ) γ ( ρ m ρ f ) ρ f β ( T w T ) ( 1 C ) is the bioconvection Rayleigh number, and ε = a / c stands for the ratio of the velocity parameter. S b = c x 2 Re x 2 1 + n D m denotes the bioconvection Schmidt number, σ = N N w N is the bioconvection constant, and P e = b W c D m represents the bioconvection Peclet number.
Various engineering quantities of interest like local skin friction Cf, the local Nusselt number Nux, and the local density of the motile microorganisms’ number Nnx are explored for the present nanofluid flow model. These quantities are defined as follows:
C f = 2 τ w ( c x ) 2 ρ f ;   N u x = q w x k f T w T ;   N n x = q n x D n T w T
For more clarification, the expressions of the shear stress τ w , wall flux qw (i.e., heat flux), and qn (i.e., motile microorganisms density flux) are given as follows:
τ w = K u y y = 0 ;   q w = k f T y y = 0 ;   q n = D n N y y = 0
By invoking the transformations of Equation (8), Equations (15) and (16) are reduced as follows:
1 2 C f Re x 1 / ( 1 + n ) = F ( 0 ) n ;   N u x Re x 1 / ( 1 + n ) = θ ( 0 ) ;   N n x Re x 1 / ( 1 + n ) = χ ( 0 )

3. Numerical Solution and Code Validation

Using MAPLE 22, Equations (8)–(11) subject to boundary conditions (12) were solved numerically. In order to solve boundary value issues numerically, this software, by default, employs a four-fifths order Runge–Kutta–Fehlberg approach. Its reliability and precision have been repeatedly demonstrated in numerous heat transfer articles. The unity coefficient of the term was changed to a continuation ( 101 100 λ ) , or ( 10 9 λ ) was used in the dsolve command in order to speed convergence for all values of the governing parameters selected for this investigation. Without making this adjustment, MAPLE produces results that do not conform to the asymptotic values, but produces results that have a sharp angle at which the axis intersects. The numerical solution to challenging ODE boundary value problems in Maple’s help section contains more details on resolving the convergence challenges. Using a value of 8 for the similarity variable η max , the asymptotic boundary conditions from Equation (12a,b) were substituted as follows.
F ( 8 ) = ε , θ ( 8 ) = 0 , ϕ ( 8 ) = 0 , χ ( 8 ) = 0
The selection of η max = 8 guaranteed that all numerical solutions appropriately approximated the asymptotic values. This is a crucial feature that is frequently missed in the literature on boundary layer fluxes.
To authenticate the model’s validity, we have compared the skin friction values for several values ε in Table 1, Table 2 and Table 3, while Table 4 compares the heat transfer values with the existing literature. It exhibits good agreement for various parameters, indicating that our numerical solution is valid.

4. Results and Discussion

In this study, the mixed bioconvective stagnation-point flow of a power-law nanofluid over a stretchy sheet was computationally studied using the Runge–Kutta–Fehlberg method of the seventh order (RKF7) in conjunction with the shooting method. The influence of gyrotactic microorganisms at the surface is taken into account. Figure 2a shows the effects of the bioconvection Rayleigh number and buoyancy parameter on the dimensionless velocity, with all other parameters held constant. We notice that the dimensionless velocity increases significantly in the vicinity of the surface and then drops to the boundary layer edge with the Rayleigh number Rb. The dimensionless velocity overshoots at the region of the surface owing to the existence of buoyant forces. The Rayleigh number increases the buoyancy forces because of bioconvection, increasing the dimensionless velocity. The effects of mixed convection and velocity ratio parameters on the dimensionless velocity are presented in Figure 2b. The convergence rate depends upon the velocity ratio parameter. As the velocity ratio increases, the convergence rate increases. Within the hydrodynamic boundary layer, the mixed convection parameter also plays an important role. When the mixed convection parameter is increased, the dimensionless velocity increases.
The impacts of nanofluid parameters Nb and Nt are explained in Figure 3a. The Brownian motion parameter Nb keeps particles moving in a fluid. This keeps particles from settling, resulting in colloidal solutions that are more stable. It helps in enhancing the dimensionless velocity within the boundary layer. On the other side, the thermophoresis parameter generates a force due to temperature difference. Nanoparticles are transported towards the lower temperature zone by this force.
Consequently, the dimensionless velocity increases inside the boundary layer. In heat transfer, the bioconvection Schmidt number is equivalent to the Prandtl number. With a Schmidt number of one, momentum and mass transfer by diffusion are alike, and the velocity and concentration boundary layers are almost identical. An increasing bioconvection Schmidt number reduces the dimensionless velocity and hence the hydrodynamic boundary layer thickness, as shown in Figure 3b.
The effects of the thermophoresis parameter on the dimensionless temperature are presented in Figure 4 for several fluids. Thermophoresis is more important in a mixed convection process, where the flow is generated by the buoyancy force caused by a temperature differential. The nanoparticles move in the direction of a temperature drop, and decreasing the bulk density improves the heat transfer process. It is worth noting that nanoparticles transport thermal energy from high-temperature areas to lower-temperature areas. As a result, the thickness of the thermal boundary layer thickens as the thermophoresis parameter Nt increases. This fact is explained in Figure 4. The Prandtl number determines the thickness of the thermal boundary layer. As the Prandtl number increases, the momentum diffusivity dominates the behavior and, as a result, the thickness of the thermal boundary layer decreases.
Figure 5a demonstrates the influence of nanofluid parameters on the dimensionless concentration at different thermophoresis parameters Nt = 0.4 and Nt = 0.5. The thermophoresis and Brownian processes cause the nanoparticle distribution to become non-uniform throughout the domain, as shown in Figure 5a. When the particle concentration is low and the Rayleigh number is low, the distribution of nanoparticles is more uniform. Nanoparticle concentration rises when Nt increases because the thermophoresis parameter represents the movement of particles due to the temperature gradient.
Figure 5b depicts the influence of the Schmidt number Sc on the dimensionless concentration for several fluids. The dimensionless concentration overshoots near the surface and drops to zero, as observed. The reason is that the dimensionless velocity rises towards the surface before falling to zero. It is worth noting that, in the surface region, the dimensionless concentration rises with the Schmidt number. This is because the mass diffusivity diminishes, resulting in a lower concentration of nanoparticles. As a result, the concentration boundary layer thickness decreases as the Schmidt number Sc distance from the cone surface increases and increases as the Schmidt number Nt decreases. However, in the passively controlled model, it is assumed that there is no nanoparticle flux at the plate and that its particle fraction value adjusts accordingly. Thus, the passively controlled nanofluid model [26] can be used in practical applications. A numerical survey is then performed for all four profiles embodying the velocity, temperature, nanoparticle volumetric fraction, and density of motile microorganisms.
The effects of the bioconvection Schmidt number for different values of Peclet number are presented in Figure 6a. The large values of the Peclet number suggest an advectively dominated distribution, while a smaller value indicates a dispersed flow. As the Peclet number increases, the boundary layer thickness of motile microorganisms increases. Conversely, the bioconvection Schmidt number tends to suppress the boundary layer thickness. It is observed that the dimensionless motile microorganism density decreases with the bioconvection Schmidt number, but increases with the Peclet number. This can be attributed to the substantial decrease in mass diffusivity, which generates a lower concentration. Figure 6b illustrates the influence of the bioconvection constant for different values of the Brownian motion parameter. Nanoparticles are not self-propelled in a nanofluid. Brownian motion and the thermophoresis effect cause them to move. Motile microorganisms are mixed with a dilute suspension of nanoparticles to increase mass transfer and microscale mixing, as well as nanofluid stability in the flow. For this reason, the Brownian motion parameter tends to decrease the dimensionless microorganisms while the bioconvection constant tends to increase the microorganism’s density.
For several values of the bioconvection Rayleigh number and dimensionless mixed convection parameter, the variation in skin friction with the buoyancy ratio parameter is shown in Figure 7a. The skin friction is observed to rise with the bioconvection Rayleigh number and dimensionless mixed convection parameter, but decreases slightly with the buoyancy parameter. The skin friction is exceptionally high in the absence of the buoyancy effect and subsequently tends to decrease along with the buoyancy parameter. However, compared with buoyancy force convection, increasing the bioconvection Rayleigh number and dimensionless mixed convection parameter enhanced the convection heat. These results are anticipated; in the interim, heat is produced as a result of increased skin friction, resulting in the formation of a layer of hot fluid close to the surface. It is interesting to note that, as long as skin friction values are positive, we can see that drag force is imparted to the surface via the fluid. Figure 7b explains the effects of nanofluid parameters on the dimensionless heat transfer for different fluids. It is important to note that the Nusselt number decreases with both nanofluid parameters, but increases with the Prandtl number. It is commonly known that the Nusselt number represents the ratio of convective to conductive heat transfer; as a result, with more significant Nusselt numbers, heat convection predominates and the Nusselt number decreases as Nb and Nt increase.
The ratio of momentum diffusivity to mass diffusivity is known as the Schmidt number. In heat transmission, this is like the Prandtl number. When there is simultaneous momentum and mass transmission, this is used to characterize flows. The Sherwood number measures the efficacy of mass transfer at the surface. The variation in the Sherwood number with the Schmidt number is depicted in Figure 8a for different values of the thermophoresis parameter and Brownian number. It is perceived that the Sherwood number increases significantly with the Schmidt number and Brownian motion parameter. However, the thermophoresis parameter tends to reduce the Sherwood number.
The variation in the dimensionless local density number of the motile microorganisms with bioconvection parameters is depicted in Figure 8b. The motile microorganism mass transfer rate is significantly increased when the bioconvection Schmidt number increases. The bioconvection Peclet number and bioconvection constant reduce the rate of motile microorganism mass transfer. This is because Pe is directly related to W c (maximum cell swimming speed) and inversely proportional to D m (the diffusivity of microorganisms). As a result, larger Pe values diminish microorganism speed and reduce microorganism diffusivity. This will result in lower microorganism concentrations in the border layer and a higher rate of motile micro-organism mass transfer.

5. Conclusions

The mixed bioconvective flow of water-based nanofluid containing micro-organisms is investigated numerically using a Runge-Kutta–Fehlberg method of the seventh order (RKF7) coupled with a shooting method. The effects of bioconvection and nanofluid parameters on the dimensionless variables and quantities of interest are investigated numerically. The results are presented graphically and are compared with the existing data. In the RKF7 method, the following are the primary conclusions that can be deduced from the study:
  • Besides the bioconvection Schmidt number, all other parameters enhance the dimensionless velocity inside the boundary layer.
  • Thermal boundary layer thickness decreases with the increasing Prandtl number.
  • The nanofluid parameters generate non-uniform nanoparticle distribution throughout the domain.
  • Dimensionless motile microbe density decreases as the bio-convection Schmidt number rises, but rises when the Peclet number falls.
  • Skin friction is slightly reduced by the buoyancy parameter, but is slightly increased by the dimensionless mixed convection parameter and the bioconvection Rayleigh number.
  • Nanofluid parameters reduce the Nusselt number, while the Prandtl number enhances it.
  • Schmidt number and Brownian motion parameter both significantly boost the Sherwood number. Thermophoresis, on the other hand, tends to lower the Sherwood number.
  • The bioconvection Peclet number and bioconvection constant contribute to slowing down the rate of mass transfer in motile microorganisms.

Author Contributions

All authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a and cpositive constants
Cconcentration
Cflocal skin-friction coefficient
Dndiffusivity of the microorganisms
DBBrownian diffusion coefficient
DTthermophoretic diffusion coefficient of the microorganisms
F′dimensionless velocity
ggravitational acceleration
Grxlocal Grashof number
Sbbioconvection Schmidt number
ScSchmidt number
Kconsistency coefficient
kfthermal conductivity
Ndensity of motile microorganisms
npower-law fluid
NbBrownian motion number
Nrbuoyancy ratio parameter
Ntthermophoresis number
NuxNusselt number
Nnxdensity number
PebioconvectionPéclet number
PrPrandtl number
qwall heat flux
Rbbioconvection Rayleigh number
Rexlocal Reynolds number
Ttemperature
u and vdimensionless velocity component in the x-direction
Uw(x)stretching velocity
U(x)free stream velocity
wcmaximum cell swimming speed
xstreamwise coordinate
ytransverse coordinate
Greek Symbols
αthermal diffusivity
βcoefficient of thermal expansion
γaverage volume of a microorganism
σbioconvection constant
ηpseudo-similarity variable
θdimensionless temperature
φnanoparticle volume fraction
ψnon-dimensional stream function
λmixed convection parameter
ε ratio of velocity parameter
χdimensionless density of motile microorganisms
μdynamic viscosity
νkinematic viscosity
ρf density of the fluid
ρfdensity of the base fluid
ρpdensity of the particles
ρmdensity of the microorganism
(ρc)f heat capacity of the fluid
(ρc)peffective heat capacity of the nanoparticle material
ρdensity
Subscripts
wcondition at the wall
condition at infinity

References

  1. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. In Proceedings of the 1995 ASME International Mechanical Engineering Congress and Exposition, San Francisco, CA, USA, 12–17 November 1995; pp. 99–105. [Google Scholar]
  2. Das, K.; Putra, N.; Thiesen, P.; Roetzel, W. Temperature dependence of thermal conductivity enhancement for nanofluid. J. Heat Transf. 2003, 125, 567–574. [Google Scholar] [CrossRef]
  3. Gorla, R.S.R.; El-Kabeir, S.M.M.; Rashad, A.M. Heat transfer in the boundary layer on a stretching circular cylinder in a nanofluid. J. Thermophys. Heat Transf. 2011, 25, 183–186. [Google Scholar] [CrossRef]
  4. Chamkha, A.J.; Rashad, A.M. Natural convection from a vertical permeable cone in nanofluid saturated porous media for uniform heat and nanoparticles volume fraction fluxes. Int. J. Numer. Methods Heat Fluid Flow 2012, 22, 1073–1085. [Google Scholar] [CrossRef]
  5. Chamkha, A.J.; Rashad, M.; Gorla, R.S.R. Non-Similar Solutions for Mixed Convection along with a Wedge Embedded in a Porous Medium Saturated by a Non-Newtonian Nanofluid: Natural Convection Dominated Regime. Int. J. Numer. Methods Heat Fluid Flow 2014, 24, 1471–1786. [Google Scholar] [CrossRef]
  6. Buongiorno, J. Convective transport in nanofluids. ASME J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  7. Khan, W.A.; Rashad, A.M.; Abdou, M.M.M.; Tlili, I. Natural Bioconvection Flow of A Nanofluid Containing Gyrotactic Microorganisms About A Truncated Cone. Eur. J. Mech./B Fluids 2019, 75, 133–142. [Google Scholar] [CrossRef]
  8. Chamkha, A.J.; Nabwey, H.A.; Abdelrahman, Z.M.A.; Rashad, A.M. Mixed Bioconvective Flow Over a Wedge in Porous Media Drenched with a Nanofluid. J. Nanofluids 2020, 9, 24–35. [Google Scholar] [CrossRef]
  9. Nabwey, H.A.; El-Kabeir, S.M.M.; Rashad, A.M.; Abdou, M.M.M. Viscous Dissipation and Joule Heating Effects on MHD Bioconvection Flow of a Nanofluid Containing Gyrotactic Microorganisms over a Vertical Isothermal Cone. J. Nanofluids 2020, 9, 242–255. [Google Scholar] [CrossRef]
  10. Nabwey, H.A.; El-Kabeir, S.M.M.; Rashad, A.M.; Abdou, M.M.M. Gyrotactic Microorganisms Mixed Convection Flow of Nanofluid over a Vertically Surfaced Saturated Porous Media. Alex. Eng. J. 2021, 61, 1804–1822. [Google Scholar] [CrossRef]
  11. Irvine, T.F.; Karni, J. Non-Newtonian fluid flow and heat transfer. In Handbook of Single-Phase Convective Heat Transfer; Kakac, S., Shah, R.K., Aung, W., Eds.; Wiley: New York, NY, USA, 1987; pp. 20.1–20.57. [Google Scholar]
  12. Schowalter, W.R. The application of boundary-layer theory to power-law pseudoplastic fluids: Similar solutions. AICHE J. 1960, 6, 24–28. [Google Scholar] [CrossRef]
  13. Acrivos, A.; Shah, M.J.; Petersen, E.E. Momentum and heat transfer in laminar boundary layer flows of non-Newtonian fluids past external surfaces. AICHE J. 1960, 6, 312–318. [Google Scholar] [CrossRef]
  14. Rashad, A.M.; EL-Hakiem, M.A.; Abdou, M.M.M. Natural convection boundary layer of a non-Newtonian fluid about a permeable vertical cone embedded in a porous medium saturated with a nanofluid. Comput. Math. Appl. 2011, 62, 3140–3151. [Google Scholar] [CrossRef]
  15. Rashad, A.M.; Chamkha, A.J.; Abdou, M.M.M. Mixed Convection Flow Of Non-Newtonian Fluid From Vertical Surface Saturated In A Porous Medium Filled with A Nanofluid. J. Appl. Fluid Mech. 2013, 6, 301–309. [Google Scholar]
  16. Chamkha, A.J.; Abbasbandy, S.; Rashad, A.M. Non-Darcy Natural Convection Flow of a Non-Newtonian Nanofluid over a Cone Saturated in a Porous Medium with Uniform Heat and Volume Fraction Fluxes. Int. J. Numer. Methods Heat Fluid Flow 2015, 25, 422–437. [Google Scholar] [CrossRef]
  17. Nabwey, H.A.; Boumazgour, M.; Rashad, A.M. Group Method Analysis of Mixed Convection Stagnation-Point Flow Of Non-Newtonian Nanofluid over a Vertical Stretching Surface. Indian J. Phys. 2017, 91, 731–742. [Google Scholar] [CrossRef]
  18. Platt, J.R. Bioconvection Patterns in Cultures of Free-Swimming Organisms. Science 1961, 133, 1766–1767. [Google Scholar] [CrossRef]
  19. Kuznetsov, A.V.; Avramenko, A.A. Effect of small particles on the stability of bioconvection in a suspension of gyrotactic microorganisms in a layer of finite depth. Int. Commun. Heat Mass Transf. 2004, 31, 1–10. [Google Scholar] [CrossRef]
  20. Kuznetsov, A.V. The onset of nanofluidbioconvection in a suspension containing both nanoparticles and gyrotactic microorganisms. Int. Commun. Heat Mass Transf. 2010, 37, 1421–1425. [Google Scholar] [CrossRef]
  21. Waqas, H.; Khan, S.U.; Hassan, M.; Bhatti, M.M.; Imran, M. Analysis of the bioconvection flow of modified second-grade nanofluid containing gyrotactic microorganisms and nanoparticles. J. Mol. Liq. 2019, 291, 111231. [Google Scholar] [CrossRef]
  22. Khan, W.A.; Uddin, M.J.; Ismail, A.I. Free convection of non-Newtonian nanofluids in porous media with gyrotactic microorganisms. Transp. Porous Media 2013, 97, 241–252. [Google Scholar] [CrossRef]
  23. Khan, W.A.; Rashad, A.M. Combined effects of radiation and chemical reaction on heat and mass transfer by MHD stagnation-point flow of a Micropolar fluid towards a stretching surface. J. Niger. Math. Soc. 2017, 36, 219–238. [Google Scholar]
  24. Beg, O.A.; Uddin, M.J.; Khan, W.A. Bioconvective Non-Newtonian Nanofluid Transport in Porous Media Containing Micro-Organisms in a Moving Free Stream. J. Mech. Med. Biol. 2015, 15, 5. [Google Scholar] [CrossRef]
  25. Ishak, A.; Nazar, R.; Pop, I. Mixed convection boundary layers in the stagnation point flow toward a stretching vertical sheet. Meccanica 2006, 41, 509–518. [Google Scholar] [CrossRef]
  26. Kuznetsov, A.V.; Nield, D.A. Natural convective boundary-layer flow of a nanofluid past a vertical plate: A revised model. Int. J. Therm. Sci. 2014, 77, 126–129. [Google Scholar] [CrossRef]
  27. Mahapatra, T.R.; Nandy, S.K.; Gupta, A.S. Magnetohydrodynamic stagnation-point flow of a power-law fluid towards a stretching surface. Int. J. Non-Linear Mech. 2009, 44, 123–128. [Google Scholar] [CrossRef]
  28. Khan, M.; Malik, R.; Munir, A.; Khan, W.A. Flow and heat transfer to Siskonanofluid over a nonlinear stretching sheet. PLoS ONE 2015, 10, e0125683. [Google Scholar]
Figure 1. Schematic functionality of the flow.
Figure 1. Schematic functionality of the flow.
Mathematics 10 03285 g001
Figure 2. Variation of dimensionless velocity with (a) the buoyancy ratio parameter and bioconvection Rayleigh number and (b) the dimensionless mixed convection parameter and velocity ratio parameter.
Figure 2. Variation of dimensionless velocity with (a) the buoyancy ratio parameter and bioconvection Rayleigh number and (b) the dimensionless mixed convection parameter and velocity ratio parameter.
Mathematics 10 03285 g002
Figure 3. Variation of dimensionless velocity with (a) nanofluid parameters and (b) bioconvection Schmidt number.
Figure 3. Variation of dimensionless velocity with (a) nanofluid parameters and (b) bioconvection Schmidt number.
Mathematics 10 03285 g003
Figure 4. Variation in dimensionless temperature with the thermophoresis parameter for different fluids.
Figure 4. Variation in dimensionless temperature with the thermophoresis parameter for different fluids.
Mathematics 10 03285 g004
Figure 5. Variation in the dimensionless nanofluid concentration with (a) nanofluid parameters and (b) Schmidt number for different fluids.
Figure 5. Variation in the dimensionless nanofluid concentration with (a) nanofluid parameters and (b) Schmidt number for different fluids.
Mathematics 10 03285 g005
Figure 6. Variation of dimensionless motile micro-organism density with (a) the bioconvection Schmidt number with different Peclet numbers and (b) the bioconvection constant and Brownian motion parameter.
Figure 6. Variation of dimensionless motile micro-organism density with (a) the bioconvection Schmidt number with different Peclet numbers and (b) the bioconvection constant and Brownian motion parameter.
Mathematics 10 03285 g006
Figure 7. Variation in skin friction with (a) buoyancy ratio parameters for different values of bioconvection Rayleigh number and dimensionless mixed convection parameter and (b) variation in the Nusselt number with the Prandtl number for different values of nanofluid parameters.
Figure 7. Variation in skin friction with (a) buoyancy ratio parameters for different values of bioconvection Rayleigh number and dimensionless mixed convection parameter and (b) variation in the Nusselt number with the Prandtl number for different values of nanofluid parameters.
Mathematics 10 03285 g007
Figure 8. (a) Variation in the Sherwood number with the Schmidt number for different nanofluid parameters and (b) variation in the density number of the motile microorganisms.
Figure 8. (a) Variation in the Sherwood number with the Schmidt number for different nanofluid parameters and (b) variation in the density number of the motile microorganisms.
Mathematics 10 03285 g008
Table 1. Comparison of numerical values of F 0 at λ = 0 ,   n = 1 .
Table 1. Comparison of numerical values of F 0 at λ = 0 ,   n = 1 .
ε Ishak et al. [10]Khan and Rashad [25]Present Results
0.1−0.9694−0.96939−0.969386154
0.2−0.9181−0.91811−0.918107089
0.5−0.6673−0.66726−0.667263673
22.01752.0175032.0175028007
34.72944.72928211.751990603
Table 2. Comparison of F ( 0 ) for various values of Pr at n = 1 and λ = Nr = 0.
Table 2. Comparison of F ( 0 ) for various values of Pr at n = 1 and λ = Nr = 0.
PrWang [23]Gorla and Sidawi [24]Khan and Pop [8]Present Results
0.070.06560.06560.06630.0659
0.20.16910.16910.16910.1690
0.70.45390.53490.45390.4539
2.00.91130.91130.91130.9113
7.01.89541.89051.89541.8954
20.03.35393.35393.35393.3539
70.06.46226.46226.46216.4622
Table 3. Comparison of F ( 0 ) for various values of n and a/c at n = 1 and λ = Nr = 0.
Table 3. Comparison of F ( 0 ) for various values of n and a/c at n = 1 and λ = Nr = 0.
nMahapatra et al. [27]Present ResultsMahapatra et al. [27]Present Results
ε = 1.1ε = 1.5
0.40.10350.10431.20191.2020
0.60.11930.12380.16911.0170
0.80.14070.142100.94340.9431
1.00.16430.164120.90950.9089
1.20.18880.188710.89370.8932
1.50.22570.225040.88530.8840
Table 4. Comparison of numerical values of θ 0 at λ = ε = N b = N t = 0 ,   n = 1 .
Table 4. Comparison of numerical values of θ 0 at λ = ε = N b = N t = 0 ,   n = 1 .
Pr Khan et al. [28]Present Results
0.70.45390.45445
20.9113−0.91135
71.8954−1.89540
203.3539−3.35391
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nabwey, H.A.; A. Khan, W.; Rashad, A.M.; Mabood, F.; Salah, T. Power-Law Nanofluid Flow over a Stretchable Surface Due to Gyrotactic Microorganisms. Mathematics 2022, 10, 3285. https://doi.org/10.3390/math10183285

AMA Style

Nabwey HA, A. Khan W, Rashad AM, Mabood F, Salah T. Power-Law Nanofluid Flow over a Stretchable Surface Due to Gyrotactic Microorganisms. Mathematics. 2022; 10(18):3285. https://doi.org/10.3390/math10183285

Chicago/Turabian Style

Nabwey, Hossam A., Waqar A. Khan, A. M. Rashad, Fazal Mabood, and Taha Salah. 2022. "Power-Law Nanofluid Flow over a Stretchable Surface Due to Gyrotactic Microorganisms" Mathematics 10, no. 18: 3285. https://doi.org/10.3390/math10183285

APA Style

Nabwey, H. A., A. Khan, W., Rashad, A. M., Mabood, F., & Salah, T. (2022). Power-Law Nanofluid Flow over a Stretchable Surface Due to Gyrotactic Microorganisms. Mathematics, 10(18), 3285. https://doi.org/10.3390/math10183285

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