Next Article in Journal / Special Issue
Performance Rating and Flow Analysis of an Experimental Airborne Drag-Type VAWT Employing Rotating Mesh
Previous Article in Journal
Why Stable Finite-Difference Schemes Can Converge to Different Solutions: Analysis for the Generalized Hopf Equation
Previous Article in Special Issue
Boundary Layer Stagnation Point Flow and Heat Transfer over a Nonlinear Stretching/Shrinking Sheet in Hybrid Carbon Nanotubes: Numerical Analysis and Response Surface Methodology under the Influence of Magnetohydrodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overlapping Grid-Based Spectral Collocation Technique for Bioconvective Flow of MHD Williamson Nanofluid over a Radiative Circular Cylindrical Body with Activation Energy

by
Musawenkosi Patson Mkhatshwa
Department of Mathematics and Applied Mathematics, University of Limpopo, Private Bag X1106, Sovenga 0727, South Africa
Computation 2024, 12(4), 75; https://doi.org/10.3390/computation12040075
Submission received: 7 March 2024 / Revised: 1 April 2024 / Accepted: 3 April 2024 / Published: 5 April 2024

Abstract

:
The amalgamation of motile microbes in nanofluid (NF) is important in upsurging the thermal conductivity of various systems, including micro-fluid devices, chip-shaped micro-devices, and enzyme biosensors. The current scrutiny focuses on the bioconvective flow of magneto-Williamson NFs containing motile microbes through a horizontal circular cylinder placed in a porous medium with nonlinear mixed convection and thermal radiation, heat sink/source, variable fluid properties, activation energy with chemical and microbial reactions, and Brownian motion for both nanoparticles and microbes. The flow analysis has also been considered subject to velocity slips, suction/injection, and heat convective and zero mass flux constraints at the boundary. The governing equations have been converted to a non-dimensional form using similarity variables, and the overlapping grid-based spectral collocation technique has been executed to procure solutions numerically. The graphical interpretation of various pertinent variables in the flow profiles and physical quantities of engineering attentiveness is provided and discussed. The results reveal that NF flow is accelerated by nonlinear thermal convection, velocity slip, magnetic fields, and variable viscosity parameters but decelerated by the Williamson fluid and suction parameters. The inclusion of nonlinear thermal radiation and variable thermal conductivity helps to enhance the fluid temperature and heat transfer rate. The concentration of both nanoparticles and motile microbes is promoted by the incorporation of activation energy in the flow system. The contribution of microbial Brownian motion along with microbial reactions on flow quantities justifies the importance of these features in the dynamics of motile microbes.

1. Introduction

Investigation of non-Newtonian fluidic models is a subject of curiosity for many researchers because of their nature, engineering, and industrial solicitations, including polymeric melt, crystal growing, dilute polymer solutions, drilling muds, cosmetic products, foods, coated sheets, and glass blowing. From a previous analysis, it is clear that many constitutional correlations have been put forward to flaunt the knotty traits of non-Newtonian fluids. Such entrenchment emerges from the view that the indispensable features of non-Newtonian liquids cannot be scrutinized using the classical Navier–Stokes equations, which are solely applicable to determine Newtonian fluid traits. Numerous researchers have investigated the behaviour of non-Newtonian fluids through various geometrical configurations. Among these geometries we have cylindrical surfaces, which play a significant role in polymer processing systems. Nazar et al. [1] investigated mixed convective flow of micropolar fluids through a horizontal circular cylinder for the case of constant surface temperature. Madhavi et al. [2,3,4] and Gaffar et al. [5] scrutinized third-grade fluid flow through a horizontal cylinder surface. Gaffar et al. [6] also explored boundary layer flow and the heat relocation of tangent hyperbolic fluid across a horizontal circular cylindrical body with first-order thermal and velocity slip. Nagaraj et al. [7] investigated flow and heat relocation in non-Newtonian Eyring–Powell fluids over a horizontal circular cylindrical surface by taking into account suction/blowing and heat convective boundary constraints. Zokri et al. [8] studied the flow of non-Newtonian Jeffrey fluids over a horizontal circular cylinder in the company of mixed convection and viscous heating effects.
The non-Newtonian fluid model for Williamson fluids was introduced by Williamson [9], and it has been confirmed to be among the predominantly crucial non-Newtonian fluids due to its low viscosity and higher shear rate. In particular, the Williamson fluid model predicts that the effective viscosity will decline monotonically with an elevation in the shear rate and vice versa. This pseudo-plastic fluid deportment is one instance of the various ways in which this non-Newtonian liquid might be usable in modern technology and industry. The applicability of Williamson fluid in technologies and industries has generated a major effect on manufacturing processes. Thus, many scholars have been involved in the exploration of features of Williamson fluids in the improvement of the thermal features of this fluid. However, just like any other non-Newtonian fluid, the Williamson fluid is still incapable of fulfilling the demand for superior extremity heat transport due to its dissatisfying thermal conductivity. To upsurge thermal conductivity as well as enhance heat transport performance, the theory of adding nanoparticles (NPs) into the non-Newtonian fluid has proven to be one of the more contemporaneously innovative concepts in the recent past years. The mixture of NPs and the host fluid is called a nanofluid (NF) [10]. Various models have been put forward to scrutinize NFs, of which the Buongiorno NF model [11] is one of them. The Buongiorno [11] NF model is a two-phase model that focuses on Brownian motion and thermophoresis impacts, and this model plays a significant role when the movement of NPs remarkably impacts the fluid flow. Another notable model is the Tiwari–Das NF model [12], which is a single-phase model that is concerned with studying the volume fraction of NPs. Rashad et al. [13] used the Buongiorno NF model to analyze the combined convection flow of a Newtonian NF across a horizontal circular cylinder immersed in a penetrable medium with heat-convective boundary conditions. Using the Tiwari–Das NF model, Zokri [14] also studied the mixed convection flow of Jeffrey NFs past a horizontal cylinder by accounting for suction/injection and heat convective boundary constraints. Merkin [15] initiated the study of fluid flow moving through a horizontal circular cylinder with mixed convection, and many studies were reported thereafter. Mixed convection plays a major role in enhancing the thermal features of heat transport. It also accounts for the general instance of convection that takes place in various industrial and technological solicitations including electronic device cooling, solar energy storage, drying technology, food processing, and float glass manufacturing.
In many engineering, technological, and industrial processes, the disparity between the surface and ambient (free stream) temperatures need not be small. In the instance of huge temperature variations, the fluid and heat transfer features are remarkably influenced. In such instances, accounting for linear Boussinesq approximation and classical linear Rosseland approximation in the fluid and heat transport equations, respectively, is inadequate, and that can lead to a reduction in the quality of the results. Consequently, it is mandatory to account for nonlinear Boussinesq approximations and nonlinear Rosseland thermal radiation in such instances so that accurate results are achieved. Practical examples of nonlinear thermal convection include engineering, geophysical, and astrophysical flows. On the other hand, nonlinear Rosseland heat radiation plays a fundamental role in manufacturing processes with larger temperature differences. EL-Zahar et al. [16] employed the modified Tiwari–Das NF to investigate the magnetized mixed convective flow of hybrid NFs across a horizontal circular cylindrical surface by considering nonlinear thermal radiation. Basha et al. [17] used the Buongiorno NF model to analyze the mixed convection flow of tangent hyperbolic NFs through a horizontal circular cylinder by accounting for nonlinear Boussinesq estimation. Activation energy (AE) is the minimal quantity of energy needed by chemical reactants to go through a chemical reaction. Mass transport procedures together with chemical reactions with Arrhenius AE emerge in areas of chemical engineering and food processing. Huang [18] studied free convective flow via a horizontal circular cylindrical surface immersed in a penetrable medium with AE and suction/injection.
Bioconvection is a phenomenon whereby the flow and thermal features of the fluid are impacted by the existence of microbes. Bioconvection flows emerge in various fascinating solicitations in environmental systems, fuel cells, and biological polymer synthesis [19]. The upward swimming microorganisms give rise to unstable density stratification, and because of that, bioconvection takes place in the fluid. In other words, bioconvection is compelled by instability that emerges due to self-propelled microorganisms swimming to and interacting with NPs and buoyancy [20,21,22,23]. The amalgamation of bioconvection and NF has yielded some exceptional results in micro-fluidic devices such as micro-channels and micro-rectors. Kuznestov and Avramenko [24] initiated the use of NPs to scrutinize the bioconvection of gyrotactic microorganisms. Hill and Pedly [25] provided a review study of the hydrodynamics of bioconvection. Through the use of Navier–Stokes equations, Allouit et al. [26] investigated the bioconvection of gyrotactic microorganisms via a vertical cylinder. Mallikarjuna [27] scrutinized the mixed bioconvective flow of a NF containing gyrotatic microorganisms from a vertical slender cylindrical surface. Rashad et al. [28] investigated mixed the bioconvective flow of a Newtonian NF comprising gyrotactic microbes through a horizontal circular cylinder in the company of heat-convective and zero mass flux boundary conditions.
The literature review, particularly investigations by Rashad et al. [13,28], reveals that there is no study on the bioconvective flow of non-Newtonian Williamson NFs in thermal and mass relocation mechanisms over a cylindrical surface with the following research gaps that result in the novelties of the present analysis.
  • The complex rheology features of Williamson materials are examined using the modified Buongiorno NF model in the presence of an applied magnetic field, porous media, and variable fluid viscosity.
  • Nonlinear Boussinesq approximation and nonlinear Rosseland thermal radiation are incorporated in the momentum and energy equations, respectively, to cater for large temperature differences.
  • Heat generation/assimilation and variable thermal conductivity are incorporated in the heat equation for the characterization of energy transfer mechanisms.
  • The modified Arrhenius AE, chemical and microbial reactions, NP and microbial Brownian diffusions are considered in the NP concentration and motile microorganism equations for characterization of mass and motile microbes transfer mechanisms.
The main objective of the current analysis is to scrutinize the bioconvective flow, heat, mass, and motile microbes transfer features of the magnetized Williamson NF model containing motile gyrotactic microbes through a horizontal circular cylindrical surface encased in a porous medium. The consequences of nonlinear thermal convection, nonlinear heat radiation, heat generation/assimilation, varying fluid attributes, AE with chemical and microbial reactions, and Brownian motion for both NPs and motile microbes are further elaborated as novelties. There has been limited research conducted on the influence of microbial Brownian motion in conjunction with microbial reactions on fluid flow properties [29,30]. Moreover, the entire flow analysis has been subjected to first-order velocity slips, suction/injection, and heat-convective and zero mass flux constraints at the boundary. The consideration of mixed convective flow of non-Newtonian NF through a heated horizontal cylinder plays a fundamental role in geological, technological, and industrial processes including the thermal recuperation of oil, geothermal reservoirs, insulating wires, drag reduction, thrust bearing, assessment of aquifers, and underground nuclear waste storage. An efficient and accurate overlapping grid-based multi-domain bivariate spectral local linearization method (OMD-BSLLM) [31,32,33] is executed to procure approximate solutions to the complex fluid flow model. The impact of various pertinent variables on the flow profiles and physical quantities of engineering attentiveness is scrutinized through graphical representations.

2. Mathematical Modeling

We contemplate the two-dimensional, laminar, steady, incompressible, and viscous nonlinear mixed convective boundary layer flow of Williamson NFs comprising gyrotactic microbes that swim from a horizontal circular cylindrical surface immersed in a porosity medium. The pictorial representation of the flow model and the physical coordinate system is depicted in Figure 1, where a is the radius of the cylinder. The x ¯ and y ¯ coordinates are, respectively, measured alongside the surface of the cylinder, resuming from the lower stagnation point, and normal to the cylinder. The uniform magnetic field intensity is enforced parallel to the fluid movement. Fluctuations in density for the buoyant expression are established by utilizing non-linear Boussinesq estimation, and the basic fluid and NPs are in a thermal equilibrium state. The surface of the cylinder is heated via convection from the heated fluid at the temperature T f , which produces a heat transport coefficient h f . The cylinder’s surface is maintained at a constant temperature T f , which is assumed to be greater than the free stream (ambient) temperature ( T f > T ) for the warmed surface (aiding flow) and lower than the ambient temperature ( T f < T ) for the cooled surface (reversing flow). The velocity of the external flow is represented by u ¯ e ( x ¯ ) = U sin x ¯ a , where U is the ambient velocity. The surface of the cylinder is also maintained at a uniform density of motile microbes n w , whereas the uniform NP volume fraction and density of motile microbes far-off from the cylindrical surface are signified by φ and n , respectively. Due to the absence of agglomeration and accumulation of NPs, the NF suspension is dilute (i.e homogeneous dispersion is attained). Furthermore, it is presumed that NPs have no impact on the direction and velocity of the gyrotactic microbes’ swimming. It is worth noting that the velocity for motile microbes, NPs, and the basic fluid is the same. For the characterization of heat transportation, features of nonlinear thermal radiation, heat sink/source, and variable thermal conductivity are involved in the heat (energy) equation. However, for the characterization of mass and motile microbes transportation, features of AE, chemical and microbial reactions, and NP and microbial Brownian motions are invoked in the equations of NP species concentration and motile microbe density conservation. In view of the aforementioned presuppositions, the entire fundamental equations of the flow problem take the form [13,28].
u ¯ x ¯ + v ¯ y ¯ = 0 , ρ f u ¯ u ¯ x ¯ + v ¯ u ¯ y ¯ = ρ f u ¯ e d u ¯ e d x ¯ + y ¯ μ ( T ) u ¯ y ¯ + 2 Γ y ¯ μ ( T ) u ¯ y ¯ u ¯ y ¯ + σ B 0 2 + μ ( T ) k p ( u ¯ e u ¯ ) + g ( 1 φ ) ρ f β 0 ( T T ) + β 1 ( T T ) 2 sin x ¯ a
g ( ρ p ρ f ) ( φ φ ) sin x ¯ a g ( ρ m ρ f ) γ ( n n ) sin x ¯ a , u ¯ T x ¯ + v ¯ T y ¯ = 1 ρ C p y ¯ κ ( T ) T y ¯ + Q 0 ( T T ) ρ C p 1 ρ C p q r y ¯
+ τ D B φ y ¯ T y ¯ + D n n y ¯ T y ¯ + D T T T y ¯ 2 ,
u ¯ φ x ¯ + v ¯ φ y ¯ = D B 2 φ y ¯ 2 k r 2 T T p exp E b k B T ( φ φ ) + D T T 2 T y ¯ 2 ,
u ¯ n x ¯ + v ¯ n y ¯ + b W c φ y ¯ n φ y = D n 2 n y ¯ 2 k m 2 T T p exp E b k B T ( n n ) + D T T 2 T y ¯ 2 ,
where u ¯ and v ¯ are velocity segments along the x ¯ and y ¯ axes, ρ f , ρ p , ρ f , and ρ m represent the density of the fluid, NPs, base fluid, and microbes, Γ is the material constant, σ signifies electrical conductivity, B 0 stands for the intensity of the magnetic field, μ ( T ) stands for the dynamic viscosity that varies with temperature, κ ( T ) indicates thermal conductivity that varies with temperature, g represents acceleration due to gravity, k p is the penetrability of the porous medium, β 0 and β 1 signify linear and nonlinear volumetric thermal expansion coefficients, γ is the average volume of the motile microbes, Q 0 stands for the heat sink/source coefficient, C p represents the specific heat capacity, τ signifies the NF heat capacity ratio, D T is the thermophoretic diffusion coefficient, D B is the Brownian motion factor, D n is the coefficient of motile microbes diffusion, k r 2 and k m 2 are the chemical and microbial reaction rates, k B is the Boltzmann constant, p is the fitted rate constant, E b is the coefficient of AE, b is the chemotaxis constant, and W c is the speed of the swimming cell. It is notable that the third term of the momentum boundary layer equation introduces the non-Newtonian fluid’s behaviour, while the fourth term represents contributions from the porous media and magnetic field. The remaining three terms on the right-hand side of the momentum equation account for buoyancy forces due to temperature variations, NP concentration, and motile microbes. In the heat equation, the second and third terms on the right-hand side capture the contributions from heat generation/absorption and thermal radiation flux based on the Rosseland approximation. However, the last term addresses the impact of nanoparticles on the thermal behavior of the fluid, encompassing phenomena such as thermophoresis and NP and microbial Brownian motions. The second term of the NP concentration and motile microbes equations incorporates the effects of AE, chemical reactions, and microbial reactions. With regard to the Rosseland diffusion approximation, the radiative flux is furnished as [34,35]
q r = 4 σ * 3 k * T 4 y ¯ = 16 σ * T 3 3 k * T y ¯ ,
where k * represents the Rosseland extinction coefficient, σ * signifies the Stefan–Boltzmann constant, and the term 16 σ * T 3 3 k * is called radiative conductivity. Since features of first-order velocity slip, uniform suction/injection velocity, and thermal convective and zero mass flux conditions are considered at the boundary, the appropriate physical boundary conditions for the flow problem under consideration are given by
u ¯ = N 0 u ¯ y ¯ , v ¯ = V w , κ T y ¯ = h f T f T , D B φ y ¯ + D T T T y ¯ = 0 , n = n w , at y ¯ = 0 , u ¯ u ¯ e ( x ¯ ) , T T , φ φ , n n , as y ¯ ,
where V w is the mass suction/injection velocity and N 0 is the first-order velocity slip factor.
In view of Rashad et al. [13,28], the external flow velocity u ¯ e ( x ¯ ) for the boundary layer equations takes the form u ¯ e ( x ¯ ) = U sin x ¯ a , where U is the free stream velocity. The next dimensionless variables are adopted in order to aid in the attainment of the numerical solutions [28]
ξ = x ¯ a , η = R e y ¯ a , u = u ¯ U , v = R e v ¯ U , u e ( ξ ) = u ¯ e ( x ¯ ) U , θ ( ξ , η ) = T T T f T , ϕ ( ξ , η ) = φ φ φ , χ ( ξ , η ) = n n n w n ,
where R e = a U ν signifies the Reynolds number, T = T [ 1 + ( θ f 1 ) θ ( ξ , η ) ] with θ f = T f T being the surface temperature excess ratio. The correctness of the NF flow and the rate of heat transmission can be achieved by considering the temperature-variant fluid viscosity and heat conductance. Consequently, the fluid velocity follows an exponential variation with temperature and the thermal conductivity fluctuates linearly with temperature, respectively, as follows [36,37,38,39]
μ ( T ) = μ e α θ , κ ( T ) = κ ( 1 + β θ ) ,
where α is the viscosity variation parameter and β ( β > 0 for fluids such as air and water, β < 0 for fluids such as lubricant oil) is the temperature-variant thermal conductivity parameter. Applying the non-dimensional variables into Equations (1)–(5), we obtain the next non-similar PDEs
u ξ + v η = 0 , u u ξ + v u η = u e d u e d ξ + e α θ 2 u η 2 α θ η u η 1 + Λ u η
+ M + ε e α θ ( u e u ) + λ θ + Ω θ 2 N r ϕ R b χ sin ξ , u θ ξ + v θ η = 1 P r 2 θ η 2 + β θ 2 θ η 2 + β θ η 2 + 4 R d 3 P r θ η ( θ f 1 ) θ + 1 3
+ N b θ η ϕ η + N n θ η χ η + N t θ η 2 + Q θ ,
u ϕ ξ + v ϕ η = 1 S c 2 ϕ η 2 K c 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ ϕ + 1 S c N t N b 2 θ η 2 , S b u χ ξ + v χ η + P e χ η ϕ η + ( χ + ω ) 2 ϕ η 2 = 2 χ η 2
K m S b 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ χ + N t N n 2 θ η 2 ,
and the corresponding boundary conditions change to
u = δ u η , v = V w R e 1 / 2 U , θ η = B i [ 1 θ ] , N b ϕ η + N t θ η , χ = 1 at η = 0 , u u e , θ 0 , ϕ 0 , χ 0 as η ,
The flow parameters appearing in Equations (10)–(14) are explicated as magnetic field parameter M = σ B 0 2 a U ρ f , Williamson fluid parameter Λ = 2 Γ R e 1 / 2 U a , Prandtl number P r = μ C p κ , thermal radiation parameter R d = 4 σ * T 3 k * κ , heat source/sink parameter Q = Q 0 a U ρ f C p , porosity parameter ε = a ν k p U , velocity slip parameter δ = R e 1 / 2 N 0 a , thermal Biot number B i = h f a R e 1 / 2 κ , buoyancy ratio parameter N r = ( ρ p ρ f ) φ ( 1 φ ) ρ f β 0 ( T f T ) , bioconvection Rayleigh number R b = ( ρ m ρ f ) γ ( n w n ) ( 1 φ ) ρ f β 0 ( T f T ) , nonlinear convection parameter Ω = β 1 ( T f T ) β 0 , NP Brownian motion parameter N b = τ D B φ ν , microbial Brownian motion parameter N n = τ D n ( n w n ) ν , thermophoresis parameter N t = τ D T ( T f T ) ν T , Schmidt number S c = ν D B , bioconvection Schmidt number S b = ν D n , AE parameter E = E b k B T , chemical reaction rate parameter K c = a k r 2 U , bioconvection Peclet number P e = b W c D n , bioconvection constant ω = n ( n w n ) , microbial reaction rate parameter K m = a k m 2 U , and mixed convection parameter λ = ( 1 φ ) g β 0 ( T f T ) a U 2 = G r R e 2 with G r = ( 1 φ ) g β 0 ( T f T ) a 3 ν 2 being the Grashof number.
Equations (9)–(14), along with boundary conditions (15), can be expressed in dimensionless form by assuming that Ψ = ξ f ( ξ , η ) , θ = θ ( ξ , η ) , ϕ = ϕ ( ξ , η ) , and χ = χ ( ξ , η ) [13,28], where Ψ signifies the dimensionless stream function given by u = Ψ η and v = Ψ ξ . Thus, the continuity equation is automatically satisfied and the momentum, energy, NP concentration, and microorganism equations convert into
e α θ f α θ f 1 + Λ ξ f f 2 + f f + M + ε e α θ sin ξ ξ f
+ sin ξ cos ξ ξ + λ θ + Ω θ 2 N r ϕ R b χ sin ξ ξ = ξ f f ξ f f ξ , 1 P r θ + β θ θ + β θ 2 + 4 R d 3 P r θ ( θ f 1 ) θ + 1 3 + f θ
+ N b θ ϕ + N n θ χ + N t θ 2 + Q θ = ξ f θ ξ θ f ξ ,
1 S c ϕ + f ϕ K c 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ ϕ + N t S c N b θ = ξ f ϕ ξ ϕ f ξ , χ + S b f χ K m S b 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ χ
P e ϕ χ + ( χ + ω ) ϕ + N t N n θ = S b ξ f χ ξ χ f ξ ,
with accompanying boundary conditions
f = δ f , f = f w , θ = B i [ 1 θ ] , N b ϕ + N t θ = 0 , χ = 1 at η = 0 , f sin ξ ξ , θ 0 , ϕ 0 , χ 0 , as η ,
where f w = V w R e 1 / 2 U is the suction parameter when f w > 0 or the injection (blowing) parameter when f w < 0 . It is noteworthy that near the lower stagnation point of the cylinder (i.e., ξ 0 ), sin ξ ξ 1 , sin ξ cos ξ ξ 1 , and Equations (16)–(19) reduce to the subsequent system of ordinary differential equations (ODEs):
e α θ f α θ f 1 + Λ ξ f f 2 + f f + M + ε e α θ 1 f + 1 + λ θ + Ω θ 2 N r ϕ R b χ = 0 ,
1 P r θ + β θ θ + β θ 2 + 4 R d 3 P r θ ( θ f 1 ) θ + 1 3 + f θ + N b θ ϕ + N n θ χ + N t θ 2 + Q θ = 0 ,
ϕ + S c f ϕ K c S c 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ ϕ + N t N b θ = 0 ,
χ + S b f χ P e ϕ χ + ( χ + ω ) ϕ + N t N n θ K m S b 1 + ( θ f 1 ) θ p exp E 1 + ( θ f 1 ) θ χ = 0 ,
and the boundary conditions diminish to
f = δ f , f = f w , θ = B i [ 1 θ ] , N b ϕ + N t θ = 0 , χ = 1 at η = 0 , f 1 , θ 0 , ϕ 0 , χ 0 , as η .
The significant physical quantities of engineering attentiveness are the local wall friction factor C f ( ξ ) , local Nusselt number N u ( ξ ) , and local density number of motile microorganisms N m ( ξ ) that can be derived using the transformations described above, and expressed in dimensionless form as follows:
C f ( ξ ) = R e 1 / 2 τ w ρ f U 2 = e α θ ( ξ , 0 ) ξ f ( ξ , 0 ) + Λ 2 ξ 2 ( f ( ξ , 0 ) ) 2 , N u ( ξ ) = R e 1 / 2 a q w k ( T f T ) = 1 + β θ ( ξ , 0 ) + 4 3 R d ( 1 + ( θ f 1 ) θ ( ξ , 0 ) ) 3 θ ( ξ , 0 ) , N m ( ξ ) = R e 1 / 2 a m w D n ( n w n ) = χ ( ξ , 0 ) ,
where τ w = μ ( T ) u ¯ y ¯ + Γ 2 u ¯ y ¯ 2 y ¯ = 0 (surface shear stress), q w = κ ( T ) + 16 σ * 3 k * T 3 T y ¯ y ¯ = 0 (surface heat flux) and m w = n y ¯ y ¯ = 0 (surface mass flux).

3. Solution Method

Numerical solutions for the non-dimensional PDEs (16)–(19) along with boundary conditions (20) have been determined by applying OMD-BSLLM. This numerical method has been successfully implemented in solving dimensional PDEs for some boundary layer flow problems with different degrees of non-linearity and complexity [31,32,33]. The choice of the numerical scheme for the present flow problem is justified by the following merits:
  • The numerical approach is remarkably accurate, convergent, stable, and resource-efficient when solving problems in both smaller and larger computational domains.
  • The numerical method uses fewer number of grid nodes and iterations to achieve results with good spectral accuracy.
  • This numerical method is well-structured, straightforward, and versatile to program in various computer software such as Matlab and Mathematica.
  • The scheme produces a less dense matrix system due to the use of the overlapping grid approach.
The numerical procedure for obtaining the numerical solutions to the present problem using the OMD-BSLLM is discussed in detail in the following section. The method is applied to linear equations, therefore it is necessary to begin by linearising the nonlinear equations using one of the linearisation techniques available in the literature. To this end, the algorithm called local linearisation method (LLM) is applied in the linearisation of the nonlinear terms in the dimensionless PDEs (16)–(19). The LLM algorithm was identified by Motsa et al. [40,41,42] as an efficient method for solving coupled systems of nonlinear ODEs and PDEs that model boundary layer flow problems. The LLM algorithm used in the simplification and decoupling of Equations (16)–(19) is summarized in Appendix A. The resulting linear iterative scheme is as follows:
a 0 , r f r + 1 + a 1 , r f r + 1 + a 2 , r f r + 1 + a 3 , r f r + 1 + a 4 , r f r + 1 ξ + a 5 , r f r + 1 ξ = a 6 , r ,
b 0 , r θ r + 1 + b 1 , r θ r + 1 + b 2 , r θ r + 1 + b 3 , r θ r + 1 ξ = b 4 , r ,
ϕ r + 1 + c 1 , r ϕ r + 1 + c 2 , r ϕ r + 1 + c 3 , r ϕ r + 1 ξ = c 4 , r ,
χ r + 1 + e 1 , r χ r + 1 + e 2 , r χ r + 1 + e 3 , r χ r + 1 ξ = e 4 , r ,
where r + 1 , r represent current and previous iterations, and the variable coefficients have been obtained using the LLM approach (see Appendix B). The next step involves discussing the decomposition of the time domain ξ [ 0 , ξ F ] , which is the solution domain for the dimensionless PDEs. The time interval 0 ξ ξ F is then segmented into ϖ non-overlapping sub-domains indicated as follows:
Δ ι = [ ξ ι , ξ ι 1 ] , 0 = ξ 0 < ξ 1 < . . < ξ ϖ 1 < ξ ϖ = ξ F , ι = 1 , 2 , 3 , , ϖ ,
where each sub-domain is further partitioned into N ξ + 1 Chebyshev collocation points (see Figure 2). The PDEs are solved in the ϖ non-overlapping sub-domains [ ξ ι , ξ ι 1 ] by starting with the easily calculated initial solutions in the first sub-interval ( Δ 1 ) . These initial solutions, which are initial conditions for the first interval, are obtained by solving the system of ODEs (21)–(24). For the remaining sub-domains, we use the solutions produced at the last grid point of each sub-interval as initial conditions for the next sub-interval. This means that the solution at the last grid point of Δ ι is used as an initial condition in Δ ι + 1 . In accordance with the validity of the spectral collocation method (SCM), which is in the interval [ 1 , 1 ] , the decomposed domain [ ξ ι , ξ ι 1 ] has to be transformed to [ 1 , 1 ] using the following linear transformation:
ξ ¯ = 2 ξ ι 1 ξ ι ξ ξ ι + ξ ι 1 ξ ι ξ ι 1 , ξ [ ξ ι 1 , ξ ι ] , ξ ^ [ 1 , 1 ] , ι = 1 , 2 , 3 , , ϖ .
Because of this transformation, we now have a new collocation variable ξ ¯ [ 1 , 1 ] , which is defined as [43]
ξ ¯ j j = 0 N ξ = cos j π N ξ .
The complete grid in the time variable ( ξ ) takes the following order:
0 = ξ N ξ ( 1 ) , , ξ 0 ( 1 ) = ξ N ξ ( 2 ) , , ξ 0 ( ι 1 ) = ξ N ξ ( ι ) , , ξ 0 ( ϖ ) = ξ F , , 2 ι ϖ .
Next, we discuss domain decomposition for the spatial variable η . Unlike the variable ξ , which is defined in the finite domain [ 0 , ξ F ] , the space variable η is specified in the semi-finite domain [ 0 , ) . In order to apply the method, this semi-finite domain has to be truncated to a finite domain [ 0 , K ] , where K is a finite value that is handpicked with the intention that it conforms to the far field boundary constraints. The resulting truncated domain is then divided into ς overlapping sub-domains indicated as
J ϑ = [ η 0 ( ϑ ) , η N η ( ϑ ) ] , ϑ = 1 , 2 , 3 , , ς ,
where each sub-domain is subsequently partitioned into N η + 1 Chebyshev collocation points (see Figure 3). Yang et al. [44] reported that using a different number of collocation points in the various sub-intervals can lead to difficulties in the implementation of the method. In addition to that condition, the length for each sub-domain must be kept the same for all sub-domains, especially if a linear transformation is utilized. This length is given by L = K ς + ( 1 ς ) ( 1 cos π N η ) / 2 (see Appendix C). Again, the decomposed domain [ η 0 ( ϑ ) , η N η ( ϑ ) ] has to be transformed to [ 1 , 1 ] as well, using the following linear mapping:
η ¯ = 2 η N η ( ϑ ) η 0 ( ϑ ) η η N η ( ϑ ) + η 0 ( ϑ ) η N η ( ϑ ) η 0 ( ϑ ) ,
and because of this transformation, we have a new collocation variable η ¯ [ 1 , 1 ] , which is defined as
η ¯ i i = 0 N η = cos i π N η .
The complete grid in the space variable ( η ) takes the following order:
0 = η N η ( 1 ) , , η 1 ( 1 ) = η N η ( 2 ) , η 0 ( 1 ) = η N η 1 ( 2 ) , , η 1 ( ϑ 1 ) = η N η ( ϑ ) , η 0 ( ϑ 1 ) = η N η 1 ( ϑ ) , , η 0 ( ς ) = K , 2 ϑ ς .
The numerical solutions at different time intervals can be differentiated using f ( ι ) , θ ( ι ) , ϕ ( ι ) , and χ ( ι ) . We remark that for the time variable, equations are solved independently in each sub-domain, but for the space variable, equations are solved contemporaneously across the whole domain of integration. Thus, the linearised iterative scheme (27)–(30) takes the form
a 0 , r 3 f ( ι ) r + 1 η 3 + a 1 , r 2 f ( ι ) r + 1 η 2 + a 2 , r f ( ι ) r + 1 η + a 3 , r f ( ι ) r + 1 + a 4 , r η f ( ι ) r + 1 ξ + a 5 , r f ( ι ) r + 1 ξ = a 6 , r ,
b 0 , r 2 θ ( ι ) r + 1 η 2 + b 1 , r θ ( ι ) r + 1 η + b 2 , r θ ( ι ) r + 1 + b 3 , r θ ( ι ) r + 1 ξ = b 4 , r ,
2 ϕ ( ι ) r + 1 η 2 + c 1 , r ϕ ( ι ) r + 1 η + c 2 , r ϕ ( ι ) r + 1 + c 3 , r ϕ ( ι ) r + 1 ξ = c 4 , r ,
2 χ ( ι ) r + 1 η 2 + e 1 , r χ ( ι ) r + 1 η + e 2 , r χ ( ι ) r + 1 + e 3 , r χ ( ι ) r + 1 ξ = e 4 , r ,
In the first interval ( Δ 1 ) , the initial conditions f ( 1 ) ( 0 , η ) , θ ( 1 ) ( 0 , η ) , ϕ ( 1 ) ( 0 , η ) , and χ ( 1 ) ( 0 , η ) are used to obtained numerical solutions f ( 1 ) ( ξ , η ) , θ ( 1 ) ( ξ , η ) , ϕ ( 1 ) ( ξ , η ) , and χ ( 1 ) ( ξ , η ) . Then, for the remaining time sub-intervals ( Δ ι , ι = 2 , 3 , 4 , , ϖ ) , the continuity conditions are given as
f ( ι ) r + 1 ( ξ ι 1 , η ) = f ( ι 1 ) r + 1 ( ξ ι 1 , η ) , θ ( ι ) r + 1 ( ξ ι 1 , η ) = θ ( ι 1 ) r + 1 ( ξ ι 1 , η ) ,
ϕ ( ι ) r + 1 ( ξ ι 1 , η ) = ϕ ( ι 1 ) r + 1 ( ξ ι 1 , η ) , χ ( ι ) r + 1 ( ξ ι 1 , η ) = χ ( ι 1 ) r + 1 ( ξ ι 1 , η ) ,
are used to generate numerical solutions f ( ι ) ( ξ , η ) , θ ( ι ) ( ξ , η ) , ϕ ( ι ) ( ξ , η ) , and χ ( ι ) ( ξ , η ) ( ι = 2 , 3 , , ϖ ) in the other sub-intervals. The procedure for generating the numerical solution is carried out by approximating the desired solutions, such as f ( ξ , η ) , using the bivariate Lagrange interpolating polynomial that appears as
f ( ι ) ( ξ , η ) F ( ι ) ( ξ , η ) = h = 0 N η l = 0 N ξ F ( ι ) ( η h , ξ l ) L h ( η ) L l ( ξ ) .
For each space sub-domain J ϑ , the first, second, and n t h order derivative matrices in the spatial direction are evaluated at the collocation points ( ξ ¯ i , η ¯ j ) for j = 0 , 1 , 2 , , N ξ as follows [33]:
f ( ι ) η ( η ¯ i , ξ j ) = D F ( ι ) j , 2 f ( ϵ ) η 2 ( η ¯ i , ξ j ) = D 2 F ( ι ) j , n f ( ϵ ) η n ( η ¯ i , ξ j ) = D n F ( ι ) j ,
where the Chebyshev differential matrix D = 2 L D ¯ , and D ¯ in the standard Chebyshev differential matrix of the first order with dimension ( N η + 1 ) × ( N η + 1 ) . However, the matrix D has dimension ( M η + 1 ) × ( M η + 1 ) , where M η = N η + ( N η 1 ) ( ς 1 ) is the combination of all collocation points used in the whole spatial region. On the other hand, F ( ι ) j = f ( ι ) ( η 0 ( ϑ ) , ξ j ) , f ( ι ) ( η 1 ( ϑ ) , ξ j ) , f ( ι ) ( η 2 ( ϑ ) , ξ j ) , , f ( ι ) ( η N η ( ϑ ) , ξ j ) T is the matrix-vector function with dimension ( N η + 1 ) × 1 . We note that the n t h ( n > 1 ) derivative matrix is obtained via matrix multiplication. For each time sub-interval Δ ι , the first-order derivative matrix in the time direction is evaluated at the collocation points ( ξ ¯ i , η ¯ j ) for i = 0 , 1 , 2 , 3 , , N η as [33]
f ( ι ) ξ ( η ¯ i , ξ j ) = l = 0 N ξ d j , l F ¯ ( ι ) l
where d j , l = 2 ξ ι ξ ι 1 d ¯ j , l ( j , l = 0 , 1 , 2 , 3 , , N ξ ) , and d ¯ j , l in the standard Chebyshev differential matrix of the first order with dimension ( N ξ + 1 ) × ( N ξ + 1 ) . It is noteworthy that the matrix-vector function F ¯ ( ι ) l has dimension ( M η + 1 ) × 1 , since it accounts for the solution in the entire spatial domain. A similar procedure is utilized in the approximation of the other unknown functions θ ( ξ , η ) , ϕ ( ξ , η ) , χ ( ξ , η ) , and their corresponding derivatives at the collocation points ( ξ ¯ i , η ¯ j ) . Because of the execution of the overlapping grid strategy, the Chebyshev differential matrix D contains many zero elements and becomes less dense, and the structure of the matrix is illustrated in [32,33,38,39] (see Appendix D). The crucial idea behind the use of SCM lies in the approximation of continuous derivatives with discrete derivatives. Consequently, by replacing the continuous derivatives in Equations (39)–(42) with the discrete derivatives discussed in the previous paragraphs and making use of the initial conditions, we obtain
a 0 , r D 3 + a 1 , r D 2 + a 2 , r D + a 3 , r F ¯ ( ι ) i , r + 1 + a 4 , r j = 0 N ξ 1 d i , j D F ¯ ( ι ) j , r + 1 + a 5 , r j = 0 N ξ 1 d i , j F ¯ ( ι ) j , r + 1 = R ( ι ) 1 , i , r ,
b 0 , r D 2 + b 1 , r D + b 2 , r Θ ¯ ( ι ) i , r + 1 + b 3 , r j = 0 N ξ 1 d i , j Θ ¯ ( ι ) j , r + 1 = R ( ι ) 2 , i , r ,
D 2 + c 1 , r D + c 2 , r Φ ¯ ( ι ) i , r + 1 + c 3 , r j = 0 N ξ 1 d i , j Φ ¯ ( ι ) j , r + 1 = R ( ι ) 3 , i , r ,
D 2 + e 1 , r D + e 2 , r χ ¯ ( ι ) i , r + 1 + e 3 , r j = 0 N ξ 1 d i , j χ ¯ ( ι ) j , r + 1 = R ( ι ) 4 , i , r ,
where the right-hand side of the equations are explicated as R ( ι ) 1 , i , r = a 6 , r a 4 , r d i , N ξ D F ¯ ( ι ) N ξ a 5 , r d i , N ξ F ¯ ( ι ) N ξ , R ( ι ) 2 , i , r = b 4 , r b 3 , r d i , N ξ Θ ¯ ( ι ) N ξ , R ( ι ) 3 , i , r = c 4 , r c 3 , r d i , N ξ Φ ¯ ( ι ) N ξ , and R ( ι ) 4 , i , r = e 4 , r e 3 , r d i , N ξ χ ¯ ( ι ) N ξ . The matching boundary constraints are also evaluated at the collocation points, and they take the form
f ( ι ) r + 1 ( η M η , ξ i ) = f w , h = 0 M η D M η , h f ( ι ) r + 1 ( η h , ξ i ) δ h = 0 M η D M η , h 2 f ( ι ) r + 1 ( η h , ξ i ) = 0 , χ ( ι ) r + 1 ( η M η , ξ i ) = 1 , h = 0 M η D M η , h f ( ι ) r + 1 ( η h , ξ i ) B i θ ( ι ) r + 1 ( η M η , ξ i ) = B i , h = 0 M η D M η , h ϕ ( ι ) r + 1 ( η h , ξ i ) + N t N b h = 0 M η D M η , h 2 θ ( ι ) r + 1 ( η h , ξ i ) = 0 , h = 0 M η D 0 , h f ( ι ) r + 1 ( η h , ξ i ) = sin ξ i ξ i , θ ( ι ) r + 1 ( η 0 , ξ i ) = 0 , ϕ ( ι ) r + 1 ( η 0 , ξ i ) = 0 , χ ( ι ) r + 1 ( η 0 , ξ i ) = 0 .
In the form of an N ξ ( M η + 1 ) × N ξ ( M η + 1 ) matrix arrangement, Equation (48) takes the form
A 0 , 0 A 0 , 1 A 0 , 2 A 0 , N ξ 1 A 1 , 0 A 1 , 1 A 1 , 2 A 1 , N ξ 1 A N ξ 1 , 0 A N ξ 1 , 1 A N ξ 1 , 2 A N ξ 1 , N ξ 1 F ¯ ( ι ) 0 F ¯ ( ι ) 1 F ¯ ( ι ) N ξ 1 = R ( ι ) 1 , 0 R ( ι ) 1 , 1 R ( ι ) 1 , N ξ 1 ,
where
A i , i = a 0 , r D 3 + a 1 , r D 2 + a 2 , r D + a 3 , r I + a 4 , r d i , i D + a 5 , r d i , i I , i = j , A i , j ( 1 ) = a 4 , r d i , j D + a 5 , r d i , j I , i j ,
a ζ , r ( η 0 , ξ i ) a ζ , r ( η 1 , ξ i ) a ζ , r ( η M η , ξ i ) , ζ = 0 , 1 , 2 , 3 , 4 , 5 ,
and I is an ( M η + 1 ) × ( M η + 1 ) unit matrix. The additional equations will also give a similar matrix arrangement of the same dimension. The resulting discrete boundary conditions are incorporated into the main diagonal sub-blocks of each matrix system to yield a new system of linear algebraic equations. The later matrix system is then solved iteratively, resuming with initial approximations that are picked in such a way that they conform to the boundary conditions. The compact form of matrix Equations (48)–(51) is given by
A ( 1 ) F = R 1 , A ( 2 ) Θ = R 2 , A ( 3 ) Φ = R 3 , A ( 4 ) χ = R 4 .
The functions used as initial guesses are given by
f 0 ( η ) = f w + 1 1 + δ e η 1 + η , θ 0 ( η ) = B i 1 + B i e η , ϕ 0 ( η ) = N t B i N b ( 1 + B i ) e η , χ 0 ( η ) = e η .

4. Results and Discussion

A comprehensive parametric study is carried out and numerical results are reported in this section through graphs. The numerical computation is performed for each parameter by assigning fixed values to other related parameters as ξ = π 2 [13,28], ω = 0.2 [28], R b = 0.1 [28], N t = 0.1 [13,28], N b = 0.3 [13,28], S c = 10 [13,28], P e = 1 [28], S b = 1 [28], N r = 0.5 [13,28], P r = 6.8 [28], δ = 0.5 [4], K c = 2 [18], B i = 5 [13,28], K m = 2 , M = 0.5 [4,17], f w [18] = 0.5 , R d = 0.5 [17], θ f = 1.1 [16], Λ = 1 [17], ε = 1 [28], β = 0.5 [37], Q = 0.2 [37], α = 0.3 [39], E = 2 [18], p = 0.2 [18], Ω = 0.1 [17], λ = 5 [13,28], and N n = 0.2 [29]. To corroborate the correctness of the present numerical method, the results for the heat transfer rate θ ( 0 ) were contrasted against those disclosed by Merkin [15], Nazar et al. [1], Rashad et al. [13], Rashad et al. [28], and Zahar et al. [16] in Table 1. Table 1 demonstrates that the numerical values in all considered cases were in substantial coherence. Hence, the selection of the present overlapping grid-based numerical scheme can be rationalized.
The convergence of the OMD-BSLLM is tracked by monitoring the solution errors, as defined in [32]. These errors represent solution-based discrepancies and quantify the number of accurate digits in the approximate solutions at the r-th iteration level. If the numerical scheme is converging, it is anticipated that the error norms will decrease as the number of iterations increases. Figure 4 illustrates the evolution of the solution errors for the approximate numerical solutions of f ( η , ξ ) , θ ( η , ξ ) , ϕ ( η , ξ ) , and χ ( η , ξ ) over iterations in both MD-BSLLM and OMD-BSLLM. The consistent reduction in all solution errors suggests that the numerical methods converge. Full convergence is considered attained when the convergence plots start to level off. It is evident that complete convergence is reached after approximately five iterations for all solutions in the MD-BSLLM approach and six iterations for all solutions in the OMD-BSLLM technique, with solution errors approaching 10 8 and 10 12 in the MD-BSLLM scheme and OMD-BSLLM technique, respectively. The minimal errors observed with the OMD-BSLLM approach validate its superior accuracy compared to the MD-BSLLM scheme. The precision of the overlapping grid-based scheme (OMD-BSLLM) can be assessed by evaluating the residual errors, as defined in [32]. These errors quantify how closely the numerical solutions approximate the true solution of the PDEs (16)–(19). Figure 5 is plotted to demonstrate residual error graphs at different time levels for both OMD-BSLLM and MD-BSLLM (no overlapping) schemes. In all the graphs, it is clear that each graph displays sufficiently small residual errors, and that is an indication of good accuracy in both numerical schemes. The accuracy is also noted to be uniform throughout the time domain, and that is one of the benefits of using the multi-domain approach in the implementation of SCM. Nevertheless, the OMD-BSLLM scheme proves to be more accurate than the MD-BSLLM approach in all the graphs. Furthermore, using more overlapping sub-domains with a smaller number of spatial collocation points yields more precision in the OMD-BSLLM approach.
The stability of the OMD-BSLLM scheme is assessed by analyzing the condition numbers of the coefficient matrices associated with the system of linear algebraic equations that need to be solved. The condition number of the matrix gauges the degree of sensitivity of the solution to variations in the input data and to the round-off errors incurred during the solution process. The condition numbers for the coefficient matrices ( A ( 1 ) , A ( 2 ) , A ( 3 ) , A ( 4 ) ) are presented in Table 2 when the number of collocation points and overlapping sub-intervals are varied. The overlapping grid-based numerical scheme is noted to yield small condition numbers compared to the non-overlapping grid-based scheme. The use of the overlapping grid approach results in coefficient matrices that are better conditioned compared to those generated by the non-overlapping grid-based calculations. The small condition numbers observed with the overlapping grid-based scheme indicate that the problem is well-conditioned, suggesting that the matrix representing the problem is not close to being singular. This characteristic is desirable as it ensures the stability, accuracy, and reliability of the numerical solution to the problem. As the number of overlapping sub-intervals increases and the number of space collocation points decreases, the condition numbers also decrease. This trend suggests that maximizing the number of overlapping sub-intervals while reducing the number of space collocation points further enhances the stability and accuracy of the numerical solutions. These findings affirm that the OMD-BSLLM scheme is the preferred numerical method for solving fluid flow problems akin to the one addressed in the current study.
For brevity, the consequences of other flow parameters such as ε , B i , N r , Λ , R b , N b , N t , P e , P r , S c , and S b on the fluid flow attributes and design quantities of engineering attentiveness have been omitted in the current study. This is because the contributions of these flow parameters are similar to the ones that were adequately discussed in Rashid et al. [13,28] in the scenario of Newtonian fluids and the absence of variable fluid properties, nonlinear heat convection, magnetic field, heat sink/source, nonlinear thermal radiation, bioconvection Brownian diffusion, activation energy, chemical and microbial reactions, suction/injection, and first-order velocity slips. The current study seeks to extend the works of Rashid et al. [13,28] by analyzing, from a mathematical point of view, the fluid flow properties, heat, mass, and motile microbe transfer phenomena when the above features are incorporated in the present bioconvection flow model. Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20 illustrate the impact of various parameters stemming from the aforementioned novel features, including the Williamson fluid ( Λ ) , nonlinear convection ( Ω ) , suction ( f w > 0 ) , velocity slip ( δ ) , variable viscosity ( α ) , magnetic field ( M ) , thermal radiation ( R d ) , temperature ratio ( θ f ) , variable thermal conductivity ( β ) , heat source ( Q > 0 ) , chemical reaction ( K c ) , AE ( E ) , microbial reaction ( K m ) , microbial Brownian motion ( N n ) , and motile microbe ( ω ) parameters, on dimensionless velocity ( f ) , temperature ( θ ) , NP concentration ( ϕ ) , density of motile microbes ( χ ) , wall friction factor ( C f ) , Nusselt number ( N u ) , and density number of motile microbes ( N m ) . These pertinent parameters are chosen in the ranges: Λ = 0, 1, 2, 3; Ω = 0.1, 0.5, 1, 1.5; f w = 0.1, 0.3, 0.5, 0.8; δ = 0.1, 0.3, 0.5, 0.8; α = 0.1, 0.5, 1, 1.5; M = 0.1, 1, 5, 10; R d = 0.3, 0.6, 1, 1.5; θ f = 1.1, 1.3, 1.5, 1.8; 0.1, 0.5, 1, 1.5; Q = 0.1, 0.3, 0.5, 0.8; K c = 0.1, 1, 2, 3, 5; E = 0, 1, 2, 3; K m = 0.1, 1, 3, 5; N n = 0.2, 0.3, 0.5, 0.8; and ω = 0.1, 0.3, 0.5, 0.8. The selected range of values for the present model align with those commonly found in a typical nanofluid. For parameters exhibiting a broad spectrum, we opted for values that ensured the stability of the numerical scheme. Figure 6 is plotted to demonstrate the impact of the Williamson fluid parameter ( Λ ) and nonlinear thermal convection parameter ( Ω ) on the dimensionless velocity. From Figure 6a, it is clear that velocity dissemination is diminished with increasing values of Λ . The variable Λ signifies the ratio of Williamson fluid relaxation time to retardation time. When the Williamson fluid parameter gets larger, the relaxation time of the fluid is noted to be lifted. As a result, resistance in the fluid flow direction is produced, which in turn reduces the NF velocity. Basically, because of the high relaxation time, fluid particles take more time to migrate to their original position; thus, fluid viscosity is improved and the velocity of the fluid particles is diminished. Conversely, velocity dissemination is accelerated with mounting values of Ω . For larger Ω , the disparity between the surface and free stream (ambient) temperatures ( T f T ) is increased, instigating a rise in the fluid velocity.
The consequence of the suction parameter ( f w > 0 ) and first-order velocity slip parameter ( δ ) on the velocity dissemination is portrayed in Figure 7. Figure 7a demonstrates that the velocity boundary layer’s thickness is diminished with the mounting values of the suction parameter. This deportment is owed to the fact that suction upshot plays a fundamental role in the removal of the fluid from the flow system; thus, momentum boundary layer thickness is remarkably decelerated. Figure 7b shows that dimensionless velocity is maximized with the growing velocity slip variable. This is on account of the acceleration in the fluid flow and the corresponding momentum boundary layer thickness as the velocity slip effect is increased. The acceleration of the fluid velocity component with improved velocity slip impact was also disclosed by Gaffar et al. [6] and Prasad et al. [45] in the absence of microbes.
The variation of dimensionless velocity against the transverse coordinate ( η ) when the variable viscosity parameter ( α ) and the applied magnetic field parameter ( M ) are varied is illustrated in Figure 8. It is clear that the velocity field is improved with rising values of the variable viscosity parameter (see Figure 8a). Physically, increasing values of α tend to augment the temperature variation ( T f T ) , consequently, the Williamson fluid bond is weakened and the intensity of the dynamic viscosity of the non-Newtonian material is minimized. As a result, the momentum boundary layer thickness is improved, and the fluid velocity is also accelerated. On the other hand, the dimensionless velocity of our nonlinear mixed convective flow is seen to be increased with mounting values of applied magnetic field parameter. This behaviour is opposite to that of the free convective flows (see Gaffar et al. [6]), where the velocity is declined by rising values of M . From a physics perspective, in an electrically conducting fluid, the magnetic field operating in a transverse direction to the geometry generates a resistive force known as the Lorentz force. The resulting Lorentz force inhibits the velocity of the fluid, as is evident in the free convective flows. In the present problem, it is noteworthy that the velocity component portrays the mounting behaviour because of the existence of nonlinear mixed convection. The similar behaviour is adequately explained by Basha et al. [17] in the absence of gyrotactic microorganisms.
Figure 9 presents the effect of the thermal radiation parameter ( R d ) and surface temperature excess ratio ( θ f ) on the dimensionless temperature curves. From the first figure, it is noted that the temperature curve is enhanced by rising values of the thermal radiation parameter. This is because increments in thermal radiation energize the boundary layer and upsurge the heat radiation energy; consequently, both the fluid temperature and thermal boundary layer width are also improved. The Rosseland radiation absorption is enhanced by the rising values of the heat radiation parameter R d ; thus, radiative heat transmission rate to the fluid is upsurged, leading to an improvement in the temperature dispersion. On the other hand, the non-dimensional temperature is also increased by the escalating values of the temperature ratio parameter since the parameter θ f improves the fluid’s thermal state. From a physics point of view, the greater values of the surface temperature excess ratio imply an upsurge in the temperature difference ( T f T ) , consequently, the temperature of the fluid is elevated. Figure 10 illustrates the ramifications of the variable thermal conductance parameter ( β ) and heat generation parameter ( Q > 0 ) on the temperature curves. Figure 10a reveals that the temperature curve improves with the escalation in β , which is owed to the excess amount of heat transmitted from the surface of the cylinder to the fluid. The parameter β is directly proportional to the temperature difference. As a result, an elevation in the variable thermal conductance parameter causes a remarkable potential for thermal energy to convey more heat energy to the NPs, instigating a rise in the fluid temperature. Figure 10b shows that temperature distribution mounts with the rising values of the heat source parameter. It is worth pointing out that Q > 0 represents the amount of heat produced within the system during the fluid flow. Thus, an increment in Q means that more heat is absorbed by the fluid, leading to an enhancement in the thermal boundary layer’s thickness and fluid temperature.
Figure 11 reveals the impact of the non-dimensional AE parameter ( E ) and chemical reaction rate parameter ( K c ) on the non-dimensional NP species concentration distributions. From Figure 11a, it is clear that increasing the non-dimensional chemical reaction rate lessens the thickness of the non-dimensional concentration dissemination. The thinning in the non-dimensional species concentration dissemination is continuously occurring in conjunction with an enlarged non-dimensional concentration gradient. On the other hand, the rise in the non-dimensional AE impact decreases the non-dimensional wall concentration gradient. The physical justification for such a behaviour is that high AE slows down the chemical reaction process. As a result, both the non-dimensional species concentration profiles (see Figure 11b) and the thickness of the concentration boundary layer improve, leading to a reduction in the non-dimensional wall concentration gradient.
The impact of the AE parameter ( E ) and microbial reaction rate parameter ( K m ) on the non-dimensional density of the motile microbes is portrayed in Figure 12. It is perceived from Figure 12a that the density of the motile microbes is reinforced by the rising values of the AE parameter. Since the microbial reaction process is improved by using AE, the motile microbe concentration profile is augmented to a greater extent. On the other hand, Figure 12b illustrates that the density of motile microbes declines with the augmenting values of the microbial reaction rate parameter. With an elevation in K m , the microorganisms’ diffusivity shrinks, for which less motile microorganism transfer occurs. In this phenomenon, the concentration of motile microorganisms drops as the values of K m increase.
Figure 13 reveals the non-dimensional microorganism concentration dissemination for the different values of the motile microbe parameter ( ω ) and microbial Brownian motion parameter ( N n ) . It is worth pointing out that the microbial Brownian motion is implied by the unrestricted mobility and unpredictable behaviour of motile microbes. The decline in the gyrotactic microorganism concentration profile for high N n indicates fewer motile microbes for the reaction (see Figure 13a). From a physics point of view, when N n gets larger, the random collision between motile microbes increases, which leads to a reduction in motile microorganism transfers. As a result, an elevation in N n resembles a decrement in the motile microbe concentration characteristics. In the other figure, there is an enhancement in the magnitude of motile microbes, in particular near the wall of the cylinder, when the bioconvection constant is raised. From formula ω = n n w n , we note that the motile microbe constant correlates the ambient density of the motile microbes to the density difference across the boundary layer. When ω upsurges, there is a larger density gradient across the boundary layer regime, which promotes the propulsion of motile microbes from the surface to the bulk flow. This results in an enhancement in the density of the motile microbes, in particular near the wall of the cylinder.
The ramification of the Williamson fluid parameter ( Λ ) and nonlinear thermal convection coefficient ( Ω ) on the wall friction factor is disclosed in Figure 14. Figure 14a demonstrates that the surface drag coefficient is a decreasing function of the changing Williamson fluid parameter. Larger values of Λ imply that a long relaxation time is needed to yield excess fluid movement resistance, thus increasing the skin friction factor. Figure 14b shows that wall shear stress is also improved with a rise in the nonlinear temperature parameter since an upsurge in Ω raises the density of the fluid flow, which in turn augments resistance to the flow. Figure 15 illustrates the upshots of the suction parameter ( f w > 0 ) and velocity slip parameter ( δ ) on the surface drag coefficient. It is perceived from Figure 15a that a greater suction parameter leads to a decrement in the skin friction factor. This demonstrates that an improvement in the cylinder surface’s porosity causes an increment in the fluid flow resistance. In the other figure, surface shear stress is also diminished by increasing the values of δ , since the fluid flow is decelerated along the cylinder when δ is raised. As noted in Prasad et al. [45], there is progressive movement in the peak shear stress locations far away from the lower stagnation when δ is elevated. Thus, the repercussion of the wall slip is noted to be significant on the boundary layer features of the Williamson fluid flow from the horizontal cylinder.
Figure 16 highlights the consequences of the variable viscosity parameter ( α ) and applied magnetic field parameter ( M ) on the wall friction factor. From the first figure, it is clear that the skin friction factor is reduced by the rise in α because of a lower resistance to flow on the cylinder surface. When there is less resistance to flow, the fluid moves more freely along the surface, causing a reduction in the shear stress and skin friction coefficient. On the other hand, a rise in the magnetic field parameter augments the surface drag coefficient. This is because the Lorentz force monitors the rate at which the fluid particles tend to move; thus, the drag coefficient increases at the level of the cylinder’s surface.
The upshots of the heat radiation parameter ( R d ) and temperature ratio parameter ( θ f ) on the dimensionless Nusselt number are depicted in Figure 17. From Figure 17a, it is clear that the elevation in R d suggests a prominent increment in the local Nusselt number. This deportment is in harmony with the fact that heat transportation becomes superior with the radiation effect. Figure 17b illustrates that the heat transfer rate escalates with increment in the surface temperature excess ratio since the presence of the temperature ratio parameter θ f = T f T contributes towards enhancement of the temperature of the fluid.
Figure 18 illustrates the distribution of the local Nusselt number for various variable thermal conductance ( β ) and heat generation parameters ( Q > 0 ) . The local Nusselt number is enhanced by rising values of β as seen in Figure 18a. This behaviour is due to the fact that the cylinder distributes additional heat by intensifying β , which in turn accelerate the heat distribution in the system as well as the rate of heat transport. Figure 18b demonstrates that escalation in Q gives rise to higher Nusselt numbers. This is because the mechanisms of heat source produce a layer of warmed fluid near the cylindrical surface; thus, the temperature of the fluid surpasses the temperature of the cylinder’s surface.
The variation in the dimensionless local density number of the motile microbes with the AE parameter ( E ) and microbial reaction parameter ( K m ) is depicted in Figure 19. It is evident that the density number of motile microbes drops when the dimensionless AE parameter is increased. This signifies that as the Reynolds number escalates, the diffusion rate prevails over the motile microbe transfer rate. The lower the density number of motile microbes, the more AE there is. Higher AE implies that more energy is needed for the motile microbes to have a fruitful collision. On the other hand, the density number of motile microorganisms is noted to be remarkably enhanced by strong motile microbial reactions. Figure 20 elucidates the impact of the bioconvection constant ( ω ) and microbial Brownian motion parameter ( N n ) on the density number of motile microbes. Figure 20a reveals that as N n augments, the density number of motile microbes is enhanced. This might be attributed to the fact that motile microbes in the Williamson fluid are moving haphazardly. The other figure demonstrates that the motile microbe transfer rate is reduced with the growing values of the motile microbes parameter.

5. Concluding Remarks

A mathematical model has been described for the bioconvective flow of magneto-Williamson NF comprising motile microbes through a radiative horizontal circular cylinder in a porous medium with nonlinear Boussinesq approximation, heat generation/assimilation, variable fluid properties, AE, chemical and microbial reactions, and Brownian motion for both NPs and motile microorganisms. The modified Buongiorno NF model is employed, and the problem is modeled subject to velocity slip, suction/injection, and heat-convective and zero mass flux boundary conditions. Using appropriate similarity transformations, a non-dimensional partial differential boundary value problem has been derived for the transport features. The resulting coupled partial differential equation system has been solved numerically with appropriate boundary conditions at the wall and in the free stream, using an overlapping grid-based SCM that employs the local linearisation technique to deal with non-linearity in the flow problem. The findings confirmed that the overlapping grid-based numerical scheme should be the most preferred numerical method for solving flow problems similar to the one considered in the current work. The remaining prominent numerical results of the present investigation are recapitulated as follows:
  • Nonlinear thermal convection, velocity slip, magnetic field, and variable viscosity parameters accelerate NF motion, while the Williamson fluid and suction parameters decelerate the fluid flow.
  • The surface drag coefficient is lifted for the Williamson fluid, nonlinear convection, and magnetic field parameters, but decays for suction, velocity slip, and variable viscosity.
  • The NF temperature and heat relocation rate are enhanced by rising thermal radiation, temperature ratio, and fluctuating thermal conductivity parameters.
  • The inclusion of AE assists in raising concentrations of NPs and motile microbes.
  • The incorporation of microbial reactions and Brownian motion contributes to a rise in motile microbes transfer rate and a reduction in the density of motile microbes.
In the near future, the current problem can be extended to the bioconvection flow of other non-Newtonian fluid models by considering the Tiwari–Das nanofluid model with other features such as electro-magnetohydrodynamics, viscous dissipation, Joule heating, and microbial growth. To optimize thermal processes, analysis of entropy generation will also be studied.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The author declares no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Latin SymbolsLatin Symbols
u ¯ , v ¯ Dimensional velocity segments along x ¯ and y ¯ axis ( m s 1 ) k p permeability of the porous medium ( m 2 )
u ¯ e velocity of the external flow ( m s 1 ) u , v Non-dimensional velocity components ( m s 1 )
U Free stream (ambient) velocity ( m s 1 ) G r Grashof number
V w Mass suction/injection velocity ( m s 1 ) B i Thermal Biot number
k r 2 , k m 2 Chemical and microbial reaction rates ( s 1 ) N u Local Nusselt number
aRadius of the cylinder ( m ) C f Skin friction vcoefficient
k B Boltzmann constant ( eV K 1 ) N m Local density number of motile microbes
Q 0 Heat source/sink coefficient ( W m 3 K 1 ) pfitted rate constant
N 0 First order velocity slip factor ( m s 1 ) N b , N n NP and microbes Brownian diffusion parameters
TFluid temperature ( K ) MMagnetic field parameter
T f , T Surface and free stream temperature ( K ) S c Schmidt number
bChemotaxis constant (m) K c , K n Chemical and microbes reaction rate parameters
k * Mean absorption variable ( m 1 ) P e bioconvection Peclet number
m w surface motile microbes flux ( kg m 2 s 1 ) h f Heat transfer coefficient
nMotile microbes concentration ( kg m 3 ) f w Suction/blowing parameter
n w , n Surface and free stream concentrations of
motile microbes ( kg m 3 )
R e Local Reynolds number
D B Brownian motion factor ( m 2 s 1 ) fDimensionless stream function
C p Specific heat capacity ( J kg 1 K 1 ) QHeat generation/absorption parameter
D T Thermophoretic diffusion coefficient ( m 2 s 1 ) R b Bioconvection Rayleigh number
q w Surface heat flux ( W m 2 ) P r Prandtl number
W c Maximum cell swimming speed ( m s 1 ) N r Buoyancy ratio parameter
gAcceleration due to gravity ( m s 2 ) Enon-dimensional activation energy parameter
D n Diffusivity of motile microbes ( m 2 s 1 ) R d Thermal radiation parameter
E b Coefficient of activation energy (eV) N t Thermophoresis force constant
B 0 Magnetic field strength ( kg s 2 A 1 ) S b Bioconvection Schmidt number
Greek SymbolsGreek Symbols
γ Average volume of motile microbes ( kg m 3 ) θ f Temperature ratio parameter
β 0 , β 1 Linear and nonlinear volumetric
thermal expansion coefficients ( K 1 )
η Pseudo-similarity variable
κ Thermal conductivity ( W m 1 K 1 ) Λ Williamson fluid parameter
ν Kinematic viscosity ( m 2 s 1 ) Γ Material constant
μ Dynamic viscosity ( kg m 1 s 1 ) Ω nonlinear convection parameter
θ , ϕ , χ Dimensionless temperature, NP species and
microbes concentration
ε Porosity parameter
φ NP volume fraction ( kg m 3 ) ω Difference variable of microbes concentration
τ w Wall shear stress ( kg m 1 s 2 ) α Viscosity variation parameter
σ * Stefan-Boltzmann constant ( 1.3807 × 10 23 J K 1 ) β Variable thermal conductivity parameter
σ Electrical conductivity ( S m 1 ) τ Nanofluid heat capacity ratio
ρ Density ( kg m 3 ) λ Mixed convection parameter
φ Ambient NP volume fraction ( kg m 3 ) δ Velocity slip parameter
Ψ Stream function ( m 2 s 1 ) ξ Stream-wise (transverse) coordinate
SubscriptsSuperscripts
w , Wall and ambient conditions f , p , m Fluid, NP and motile microbes

Appendix A. LLM Algorithm

  • Solve for f in the first equation assuming that θ , ϕ and χ are known from previous iteration. f r + 1 .
  • Considering that f r + 1 is now known, solve for θ in the second equation, assuming that ϕ and χ are known from previous iteration. θ r + 1 .
  • Considering that f r + 1 and θ r + 1 are now known, solve for θ in the third equation, assuming that χ is known from previous iteration. ϕ r + 1 .
  • Lastly, considering that f r + 1 , θ r + 1 and ϕ r + 1 are now known, solve for θ in the last equation. χ r + 1 .

Appendix B. Variable Coefficients Obtained Using the LLM Approach

a 0 , r = e α θ r + Λ ξ e α θ r f r , a 1 , r = f r α θ r e α θ r + Λ ξ e α θ r f r 2 α Λ ξ θ r f r e α θ r + ξ f r ξ , a 2 , r = 2 f r ( M + ε e α θ r ) ξ f r ξ , a 3 , r = f r , a 4 , r = ξ f r , a 5 , r = ξ f r , a 6 , r = f r f r f r 2 M sin ξ ξ ε e α θ r sin ξ ξ sin ξ cos ξ ξ + Λ ξ f r f r e α θ r α Λ ξ θ r f r 2 e α θ r ξ f r f r ξ + ξ f r f r ξ λ ( θ r + Ω θ r 2 N r ϕ r R b χ r ) sin ξ ξ , b 0 , r = 1 + β θ r + 4 3 R d + 4 ( θ f 1 ) R d θ r + 4 ( θ f 1 ) 2 R d θ r 2
+ 4 3 R d ( θ f 1 ) 3 θ r 3 , b 1 , r = P r f r + 2 β θ r + 8 R d ( θ f 1 ) θ r + 16 R d ( θ f 1 ) 2 θ r θ r + 8 R d ( θ f 1 ) 3 θ r 2 θ r + N b P r ϕ r + N n P r χ r + 2 N t P r θ r + P r ξ f r ξ , b 2 , r = β θ r + 4 R d ( θ f 1 ) θ r + 8 R d ( θ f 1 ) 2 θ r θ r + 4 R d ( θ f 1 ) 3 θ r θ r 2 + 8 R d ( θ f 1 ) 2 θ r 2 + 8 R d ( θ f 1 ) 3 θ r θ r 2 + P r Q , b 3 , r = P r ξ f r , b 4 , r = 4 R d ( θ f 1 ) θ r θ r + 2 ( θ f 1 ) θ r θ r 2 + ( θ f 1 ) 2 θ r θ r 3 + θ r 2 + 4 ( θ f 1 ) θ r θ r 2 + 3 ( θ f 1 ) 2 θ r 2 θ r 2 + β θ r θ r + β θ r 2 + N t P r θ r 2 , c 1 , r = S c f r + S c ξ f r ξ , c 3 , r = S c ξ f r , c 4 , r = N t N b θ r , c 2 , r = S c K c ( 1 + ( θ f 1 ) θ r ) p exp E ( 1 + ( θ f 1 ) θ r ) , e 1 , r = S b f r P e ϕ r + S b ξ f r ξ , e 2 , r = P e ϕ r S b K m ( 1 + ( θ f 1 ) θ r ) p exp E ( 1 + ( θ f 1 ) θ r ) , e 3 , r = S b ξ f r , e 4 , r = P e ω ϕ r N t N n θ r .

Appendix C. Derivation of the Length of Each Sub-Interval (L)

To ascertain the length of individual sub-intervals, we observe that the entire length of the spatial domain is specified by
K = 2 L ϱ + ( 2 L 2 ϱ ) ς 2 1 = 2 L ϱ + ( L ϱ ) ( ς 2 ) = ϱ ( 1 ς ) + ς L ,
where ϱ represents the distance of overlap between two adjacent sub-intervals. It is important to highlight that ϱ equals the difference between η 0 and η 1 ( i.e. , ϱ = η 0 η 1 ) . When examining the first interval J 1 , where η [ 0 , η N η 1 ] , we can establish the length L = η N η 1 . Utilizing the linear transformation
η = L 2 η ¯ + η N η 1 2 ,
facilitates the conversion of the interval [ 0 , η N η 1 ] into [ 1 , 1 ] . Therefore, by employing the Gauss–Lobatto collocation points η ¯ i = cos π i N η , where i = 0 , 1 , 2 , , N η , we acquire
η 0 η 1 = L 2 ( η ¯ 0 η ¯ 1 ) = L 2 1 cos π N η .
As a result, Equation (A2) is transformed into
K = L 2 1 cos π N η ( 1 ς ) + ς L ,
and rearranging the formula given in Equation (A5) to solve for L , we acquire
L = K ς + ( 1 ς ) ( 1 cos π N η ) / 2 .

Appendix D. The Structure of Chebyshev Differentiation Matrix (D)

D = D ¯ ( ς ) 0 , 0 D ¯ ( ς ) 0 , 1 D ¯ ( ς ) 0 , N η 1 D ¯ ( ς ) 0 , N η D ¯ ( ς ) 1 , 0 D ¯ ( ς ) 1 , 1 D ¯ ( ς ) 1 , N η 1 D ¯ ( ς ) 1 , N η D ¯ ( ς ) N η 1 , 0 D ¯ ( ς ) N η 1 , 1 D ¯ ( ς ) N η 1 , N η 1 D ¯ ( ς ) N η 1 , N η D ¯ ( ς 1 ) 1 , 0 D ¯ ( ς 1 ) 1 , 1 D ¯ ( ς 1 ) 1 , N η 1 D ¯ ( ς 1 ) 1 , N η D ¯ ( ς 1 ) 2 , 0 D ¯ ( ς 1 ) 2 , 1 D ¯ ( ς 1 ) 2 , N η 1 D ¯ ( ς 1 ) 2 , N η D ¯ ( ς 1 ) N η 1 , 0 D ¯ ( ς 1 ) N η 1 , 1 D ¯ ( ς 1 ) N η 1 , N η 1 D ¯ ( ς 1 ) N η 1 , N η D ¯ ( 1 ) 1 , 0 D ¯ ( 1 ) 1 , 1 D ¯ ( 1 ) 1 , N η 1 D ¯ ( 1 ) 1 , N η D ¯ ( 1 ) 2 , 0 D ¯ ( 1 ) 2 , 1 D ¯ ( 1 ) 2 , N η 1 D ¯ ( 1 ) 2 , N η D ¯ ( 1 ) N η , 0 D ¯ ( 1 ) N η , 1 D ¯ ( 1 ) N η , N η 1 D ¯ ( 1 ) N η , N η ,

References

  1. Nazar, R.; Amin, N.; Pop, I. Mixed convection boundary-layer flow from a horizontal circular cylinder in micropolar fluids: Case of constant wall temperature. Int. J. Numer. Methods Heat Fluid Flow 2003, 13, 86–109. [Google Scholar] [CrossRef]
  2. Madhavi, K.; Nagendra, N.; Raju, G.S.S.; Prasad, V.R. Non-similarity transformations of MHD flow of third-grade fluid over a horizontal circular cylinder in the presence of radiation and suction. Int. J. Innov. Res. Sci. Eng. Technol. 2017, 6, 73–79. [Google Scholar]
  3. Madhavi, K.; Nagendra, N.; Raju, G.S.S.; Prasad, V.R. Heat transfer analysis of MHD non-Newtonian fluid over a horizontal circular cylinder with Biot number effect. Res. J. Sci. Technol. 2017, 9, 395–399. [Google Scholar] [CrossRef]
  4. Madhavi, K.; Prasad, V.R.; Gaffar, S.S.; Venkatadri, K. Entropy analysis of third-grade MHD convection flows from a horizontal cylinder with slip. Arch. Mech. Eng. 2018, 65, 417–440. [Google Scholar]
  5. Gaffar, S.A.; Prasad, V.R.; Reddy, P.R.; Khan, B.M.H. Radiative flow of third grade non-Newtonian fluid from a horizontal circular cylinder. Nonlinear Eng. 2019, 8, 673–687. [Google Scholar] [CrossRef]
  6. Gaffar, S.A.; Prasad, V.R.; Beg, O.A. Computational analysis of magnetohydrodynamic free convection flow and heat transfer on non-Newtonian tangent hyperbolic fluid from a horizontal circular cylinder with partial slip. Int. J. Appl. Comput. Math. 2015, 1, 651–675. [Google Scholar] [CrossRef]
  7. Nagaraja, L.; Reddy, M.S.; Reddy, M.S.N. Heat transfer in non-Newtonian fluid past horizontal circular cylinder in non-Darcy porous medium with suction/injection effects. Indian J. Sci. Technol. 2017, 10, 1–7. [Google Scholar] [CrossRef]
  8. Zokri, S.M.; Arifin, N.S.; Mohamed, M.K.A.; Kasim, A.R.M.; Mohammad, N.F.; Salleh, M.Z. Mathematical model of mixed convection boundary layer flow over a horizontal circular cylinder filled in a Jeffrey fluid with viscous dissipation effect. Sains Malays. 2018, 47, 1607–1615. [Google Scholar] [CrossRef]
  9. Williamson, R.V. The flow of psedoplastic materials. Ind. Eng. Chem. 1929, 21, 1108–1111. [Google Scholar] [CrossRef]
  10. Choi, S.U.S.; Eastman, J.A. Enhancing thermal conductivity of fluids with nanoparticles. ASME Fluids Eng. Div. 1995, 23, 99–106. [Google Scholar]
  11. Buongiorno, J. Convective transport in nanofluids. J. Heat Transfer 2006, 128, 240–250. [Google Scholar] [CrossRef]
  12. Tiwari, R.J.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  13. Rashad, A.M.; Chamkha, A.J.; Modather, M. Mixed convection boundary-layer flow past a horizontal circular cylinder embedded in a porous medium filled with a nanofluid under convective boundary condition. Comput. Fluids 2013, 86, 380–388. [Google Scholar] [CrossRef]
  14. Zokri, S.M.; Salleh, M.Z.; Arifin, N.S.; Kasim, A.R.M. Lower stagnation point flow of convectively heated horizontal circular cylinder in Jeffrey nanofluid with suction/injection. J. Adv. Res. Fluid Mech. Therm. Sci. 2020, 76, 135–144. [Google Scholar] [CrossRef]
  15. Merkin, J.H. Mixed convection from a horizontal circular cylinder. Int. J. Heat Mass Transf. 1977, 20, 73–77. [Google Scholar] [CrossRef]
  16. Zahar, E.R.; Rashad, A.M.; Saad, W.; Seddek, L.F. Magneto-hybrid nanofluids flow via mixed convection past radiative circular cylinder. Sci. Rep. 2020, 10, 10494. [Google Scholar] [CrossRef] [PubMed]
  17. Basha, H.T.; Sivaraj, R.; Prasad, V.R.; Beg, O.A. Entropy generation of tangent hyperbolic nanofluid flow over a circular cylinder in the presence of nonlinear Boussinesq approximation: A non-similar solution. J. Therm. Anal. Calorim. 2021, 143, 2273–2289. [Google Scholar] [CrossRef]
  18. Huang, C.-J. Arrhenius activation energy effect on free convection about a permeable horizontal cylinder in porous media. Transp. Porous Media 2019, 128, 723–740. [Google Scholar] [CrossRef]
  19. Das, K.; Duari, P.R.; Kundu, P.K. Nanofluid bioconvection in presence of gyrotactic microorganisms and chemical reaction in a porous medium. J. Mech. Sci. Technol. 2015, 29, 4841–4849. [Google Scholar] [CrossRef]
  20. Bees, M.; Hill, N.A. Wavelengths of bioconvection patterns. J. Exp. Biol. 1997, 200, 1515–1526. [Google Scholar] [CrossRef]
  21. Kuznetsov, A. Nanofluid bioconvection: Interaction of microorganisms oxytactic up-swimming, nanoparticle distribution, and heating/cooling from below. Theor. Comput. Fluid Dyn. 2012, 26, 291–310. [Google Scholar] [CrossRef]
  22. Pedley, T.J. Gyrotaxis in uniform vorticity. J. Fluid Mech. 2015, 762, R6. [Google Scholar] [CrossRef]
  23. Khan, S.U.; Al-Khaled, K.; Aldabesh, A.; Awais, M.; Tlili, I. Bioconvection ow in accelerated couple stress nanoparticles with activation energy: Biofuel applications. Sci. Rep. 2021, 11, 3331. [Google Scholar] [CrossRef]
  24. Kuznetsov, A.V.; Avramenko, A.A. Effect of small particles on this 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]
  25. Hill, N.A.; Pedley, T.J. Bioconvection. Fluid Dyn. Res. 2005, 37, 1–20. [Google Scholar] [CrossRef]
  26. Alloui, Z.; Nguyen, T.H.; Bilgen, E. Bioconvection of gravitactic microorganisms in a vertical cylinder. Int. Commun. Heat Mass Transf. 2005, 32, 739–747. [Google Scholar] [CrossRef]
  27. Mallikarjuna, B.; Rashad, A.M.; Chamkha, A.; Abdou, M.M.M. Mixed bioconvection flow of a nanofluid containing gyrotactic microorganisms past a vertical slender cylinder. Front. Heat Mass Transf. 2018, 10, 1–10. [Google Scholar] [CrossRef]
  28. Rashad, A.M.; Nabwey, H.A. Gyrotactic mixed bioconvection flow of nanofluid past a circular cylinder with convective boundary condition. J. Taiwan Inst. Chem. Eng. 2019, 99, 9–17. [Google Scholar] [CrossRef]
  29. Dhlamini, M.; Mondal, H.; Sibanda, P.; Motsa, S.; Sachin, S. A mathematical model for bioconvection flow with activation energy for chemical reaction and microbial activity. Pramana-J. Phys. 2022, 96, 112. [Google Scholar] [CrossRef]
  30. Mkhatshwa, M.P. Bioconvective flow of fourth-grade nanofluid via a stretchable porous surface with microbial activity and activation energy. Case Stud. Therm. Eng. 2023, 51, 103526. [Google Scholar] [CrossRef]
  31. Mkhatshwa, M.P.; Motsa, S.S.; Sibanda, P. Overlapping multi-domain bivariate spectral method for non-Darcian mixed convection chemically reacting flow of micropolar fluid over a flat surface in porous media. J. Appl. Anal. Comput. 2021, 11, 1–25. [Google Scholar]
  32. Mkhatshwa, M.P. Overlapping grid spectral collocation approach for electrical MHD bioconvection Darcy–Forchheimer flow of a Carreau–Yasuda nanoliquid over a periodically accelerating surface. Heat Transf. 2022, 51, 1468–1500. [Google Scholar] [CrossRef]
  33. Mkhatshwa, M.P.; Khumalo, M. Double diffusion and Hall effects on MHD sinusoidal natural convection flow of silver water-based nanofluid from a porous vertical plate. Partial Differ Equ. Appl. Math. 2023, 7, 100516. [Google Scholar] [CrossRef]
  34. Raptis, A. Radiation and free convection flow through a porous medium. Int. Comm. Heat Transf. 1998, 25, 289–295. [Google Scholar] [CrossRef]
  35. Hakiem, M.A.; Rashad, A.M. Effect of radiation on non-Darcy free convection from a vertical cylinder embedded in a fluid saturated porous medium with a temperature-dependent viscosity. J. Porous Media 2007, 10, 209–218. [Google Scholar]
  36. Oyelakin, I.S.; Sibanda, P. Analysis of exponentially varying viscosity and thermal conductivity on a tangent hyperbolic fluid. SeMA J. 2020, 77, 257–273. [Google Scholar] [CrossRef]
  37. Mkhatshwa, M.P.; Khumalo, M. Irreversibility scrutinization onEMHD Darcy–Forchheimer slip flow of Carreau hybrid nanofluid through a stretchablesurface in porous medium with temperature-variant properties. Heat Transf. 2023, 52, 395–429. [Google Scholar] [CrossRef]
  38. Mkhatshwa, M.P.; Motsa, S.S.; Sibanda, P. MHD mixed convective radiative flow of Eyring-Powell fluid over an oscillatory stretching sheet using bivariate spectral method on overlapping grids. Heat Transf. 2021, 50, 655–687. [Google Scholar] [CrossRef]
  39. Mkhatshwa, M.P.; Motsa, S.S.; Sibanda, P. MHD bioconvective radiative flow of chemically reactive Casson nanofluid from a vertical surface with variable transport properties. Int. J. Ambient Energy 2020, 43, 3170–3188. [Google Scholar] [CrossRef]
  40. Motsa, S.S. A new spectral local linearization method for nonlinear boundary layer flow problems. J. Applied Math. 2013, 2013, 423628. [Google Scholar] [CrossRef]
  41. Motsa, S.S. On the practical use of the spectral homotopy analysis method and local linearisation method for unsteady boundary-layer flows caused by an impulsively stretching plate. Numer. Algorithms 2014, 66, 865–883. [Google Scholar] [CrossRef]
  42. Motsa, S. On the new bivariate local linearisation method for solving coupled partial differential equations in some applications of unsteady fluid flows with heat and mass transfer. In Mass Transfer—Advancement in Process Modelling; InTech: Vienna, Austria, 2015. [Google Scholar]
  43. Trefethen, L.N. Spectral Methods in MATLAB; SIAM: Philadelphia, PA, USA, 2000. [Google Scholar]
  44. Yang, H.H.; Seymour, B.R.; Shizgal, B.D. A Chebyshev pseudospectral multi-domain method for steady flow past a cylinder up to Re = 150. Comput. Fluids 1994, 23, 829–851. [Google Scholar] [CrossRef]
  45. Prasad, V.R.; Rao, A.S.; Reddy, N.B.; Beg, O.A. Modelling laminar transport phenomena in a Casson rheological fluid from an isothermal sphere with partial slip in a non-Darcy medium. Thoret. Appl. Mech. 2013, 40, 469–510. [Google Scholar] [CrossRef]
Figure 1. Flow model and physical coordinate system.
Figure 1. Flow model and physical coordinate system.
Computation 12 00075 g001
Figure 2. Dividing the time solution domain into ϖ non-overlapping sub-intervals.
Figure 2. Dividing the time solution domain into ϖ non-overlapping sub-intervals.
Computation 12 00075 g002
Figure 3. Dividing the spatial solution domain into ς overlapping sub-intervals.
Figure 3. Dividing the spatial solution domain into ς overlapping sub-intervals.
Computation 12 00075 g003
Figure 4. Comparison of solution errors ( E f , E θ , E ϕ , and E χ ) when N ξ = 5 and ϖ = 20 . (a) Solution errors against iterations for the MD-BSLLM ( ς = 1 , N η = 100 ) . (b) Solution errors against iterations for the OMD-BSLLM ( ς > 1 , N η = 10 ) .
Figure 4. Comparison of solution errors ( E f , E θ , E ϕ , and E χ ) when N ξ = 5 and ϖ = 20 . (a) Solution errors against iterations for the MD-BSLLM ( ς = 1 , N η = 100 ) . (b) Solution errors against iterations for the OMD-BSLLM ( ς > 1 , N η = 10 ) .
Computation 12 00075 g004
Figure 5. Comparison of residual error approximations for the case of MD-BSLLM ( ς = 1 ) and OMD-BSLLM ( ς > 1 ) when N ξ = 5 and ϖ = 20 for different nodes ( [ N η , ς ] = [ 100 , 1 ] , [ 50 , 2 ] , [ 20 , 5 ] , [ 10 , 10 ] ).
Figure 5. Comparison of residual error approximations for the case of MD-BSLLM ( ς = 1 ) and OMD-BSLLM ( ς > 1 ) when N ξ = 5 and ϖ = 20 for different nodes ( [ N η , ς ] = [ 100 , 1 ] , [ 50 , 2 ] , [ 20 , 5 ] , [ 10 , 10 ] ).
Computation 12 00075 g005
Figure 6. Dimensionless velocity profiles.
Figure 6. Dimensionless velocity profiles.
Computation 12 00075 g006
Figure 7. Dimensionless velocity profiles.
Figure 7. Dimensionless velocity profiles.
Computation 12 00075 g007
Figure 8. Dimensionless velocity profiles.
Figure 8. Dimensionless velocity profiles.
Computation 12 00075 g008
Figure 9. Dimensionless temperature profiles.
Figure 9. Dimensionless temperature profiles.
Computation 12 00075 g009
Figure 10. Dimensionless temperature profiles.
Figure 10. Dimensionless temperature profiles.
Computation 12 00075 g010
Figure 11. Dimensionless nanoparticle concentration profiles.
Figure 11. Dimensionless nanoparticle concentration profiles.
Computation 12 00075 g011
Figure 12. Dimensionless density of the motile microbe profiles.
Figure 12. Dimensionless density of the motile microbe profiles.
Computation 12 00075 g012
Figure 13. Dimensionless density of the motile microbe profiles.
Figure 13. Dimensionless density of the motile microbe profiles.
Computation 12 00075 g013
Figure 14. Dimensionless surface drag coefficient.
Figure 14. Dimensionless surface drag coefficient.
Computation 12 00075 g014
Figure 15. Dimensionless surface drag coefficient.
Figure 15. Dimensionless surface drag coefficient.
Computation 12 00075 g015
Figure 16. Dimensionless surface drag coefficient.
Figure 16. Dimensionless surface drag coefficient.
Computation 12 00075 g016
Figure 17. Dimensionless Nusselt number.
Figure 17. Dimensionless Nusselt number.
Computation 12 00075 g017
Figure 18. Dimensionless Nusselt number.
Figure 18. Dimensionless Nusselt number.
Computation 12 00075 g018
Figure 19. Dimensionless density number of the motile microbes.
Figure 19. Dimensionless density number of the motile microbes.
Computation 12 00075 g019
Figure 20. Dimensionless density number of the motile microbes.
Figure 20. Dimensionless density number of the motile microbes.
Computation 12 00075 g020
Table 1. Comparison of θ ( 0 ) when P r = 1 , θ f = 1 , B i , ξ = β = α = f w = M = Λ = R d =   Q = ε = δ = N r = R b = Ω = N b = N n = N t = S c = S b = E = K c = P e = K m = ω = p = 0 .
Table 1. Comparison of θ ( 0 ) when P r = 1 , θ f = 1 , B i , ξ = β = α = f w = M = Λ = R d =   Q = ε = δ = N r = R b = Ω = N b = N n = N t = S c = S b = E = K c = P e = K m = ω = p = 0 .
λ Merkin [15]Nazar et al. [1]Rashad et al.  [13]Rashad et al. [28]Zahar et al. [16]Present Results
−1.750.41990.42050.42020.42020.41980.4199516
−1.500.45760.46010.45790.45790.45730.4573551
−1.000.50670.50800.50680.50680.50670.5066785
−0.500.54200.54300.54210.542010.54210.5420654
0.000.57050.57100.57060.57060.57050.5704700
0.500.59430.59490.59470.59470.59450.5945339
0.880.60960.61120.61110.61110.61080.6107615
0.890.61100.61160.61140.61140.61120.6111687
1.000.61580.61600.61600.61600.61540.6155872
2.000.64970.65180.65180.65180.65150.6514920
5.000.73150.73200.73190.73190.73150.7315104
Table 2. Condition numbers for varying nodes [ N η , ς ] .
Table 2. Condition numbers for varying nodes [ N η , ς ] .
ς N η M η Cond ( A ( 1 ) ) Cond ( A ( 2 ) ) Cond ( A ( 3 ) ) Cond ( A ( 4 ) )
1 (MD-BSLLM)1001001.4518640e+106.0399799e+062.5107303e+063.1173003e+06
2 (OMD-BSLLM)50991.7962336e+091.2883089e+065.4063508e+056.6439084e+05
4 (OMD-BSLLM)25972.4075014e+082.7726519e+051.1971839e+051.4291087e+05
5 (OMD-BSLLM)20961.2646992e+081.6967207e+057.4692771e+048.7440012e+04
10 (OMD-BSLLM)10911.6621778e+073.7796142e+041.9268267e+041.9579885e+04
20 (OMD-BSLLM)5811.7880098e+069.4107769e+037.4254531e+034.6950726e+03
25 (OMD-BSLLM)4768.0984196e+056.3770832e+036.1905704e+032.9795539e+03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mkhatshwa, M.P. Overlapping Grid-Based Spectral Collocation Technique for Bioconvective Flow of MHD Williamson Nanofluid over a Radiative Circular Cylindrical Body with Activation Energy. Computation 2024, 12, 75. https://doi.org/10.3390/computation12040075

AMA Style

Mkhatshwa MP. Overlapping Grid-Based Spectral Collocation Technique for Bioconvective Flow of MHD Williamson Nanofluid over a Radiative Circular Cylindrical Body with Activation Energy. Computation. 2024; 12(4):75. https://doi.org/10.3390/computation12040075

Chicago/Turabian Style

Mkhatshwa, Musawenkosi Patson. 2024. "Overlapping Grid-Based Spectral Collocation Technique for Bioconvective Flow of MHD Williamson Nanofluid over a Radiative Circular Cylindrical Body with Activation Energy" Computation 12, no. 4: 75. https://doi.org/10.3390/computation12040075

APA Style

Mkhatshwa, M. P. (2024). Overlapping Grid-Based Spectral Collocation Technique for Bioconvective Flow of MHD Williamson Nanofluid over a Radiative Circular Cylindrical Body with Activation Energy. Computation, 12(4), 75. https://doi.org/10.3390/computation12040075

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