Bioconvection is the phenomenon of pattern creation found in the aqueous solution of motile microorganisms when they respond to specific stimuli by swimming in certain directions named taxes Platt [1]. According to Ghorai and Hill [2], gyrotaxis is swimming driven by the balance of torque due to gravity acting on a bottom-heavy cell and torque owing to viscous forces resulting from local shear flows. According to the study of Raees et al. [3], bioconvection in nanofluids holds great potential in the Colibri micro-volumes spectrometer as well as in improving the stability of nanofluids. They also showed the application of bioconvection in the field of microbial enhanced oil recovery, which involves injecting selected microorganisms into the reservoir and multiplying them in situ to reduce the residual oil left in the reservoir once secondary recovery is exhausted. Waqas et al. [4] discussed the effects of activation energy and thermal radiation in the bioconvection flow of third-grade nanofluid. Ramzan et al. [5] presented bioconvective Reiner-Rivlin nanofuid flow over a rotating disk. It was reported that motile density distribution diminishes for large values of bioconvective Lewis numbers. Manjunatha et al. [6] explicated quartic autocatalysis of homogeneous and heterogeneous reactions in the bioconvective flow of radiating micropolar nanofluid. Koriko et al. [7] exploited active and passive controls of bioconvection flow of thixotropic fluid with Catteno-Christov phenomenon. They reported that density of motile microorganisms profile diminish for stronger values of bioconvective parameter for both active and passive controls of nanoparticles. MHD bioconvection flow of Casson nanofluid over a rotating disk was explored by Jawad et al. [8]. It was reported that density of motile gyrotactic microorganisms distribution declines with larger values of bioconvection Peclet number.

For many years, the dynamics of fluid with small particles has been a major topic of research in industry. It is a well-known fact that each particle in this type of fluid typically rotates independently of the fluid’s velocity. Micropolar fluids are viscous fluids composed of rigid, randomly oriented particles suspended in a viscous liquid. Eringen [9] proposed the concept of micropolar fluid in an attempt to characterize the behavior of a specific fluid including polymeric additives and naturally occurring fluids such as colloidal fluid flow, liquid crystals, and animal blood. Heat and mass transfer of MHD micropolar fluid flow along a permeable channel was examined by Agarwal et al. [10]. Shamshuddin et al. [11] studied characteristics of thermophoresis and Brownian motion on radiative reactive micropolar fluid flow towards continuously moving flat plate. Numerical simulation of stagnation point flow in magneto micropolar fluid over a stretchable surface under influence of activation energy and bilaterial reaction was explored by Shamshuddin et al. [12]. Influence of viscous dissipation on MHD flow of micropolar fluid over a slendering stretching surface with modified heat flux model was presented by Kumar et al. [13]. Kumar et al. [14] considered physical aspects on unsteady MHD free convective stagnation point flow of micropolar fluid over a stretching surface. Koriko et al. [15] analysed boundary layer flow of micropolar fluid with combined nonlinear thermal radiation and thermal stratification. They observed that when there is a high concentration of microstructures, the rate of rise in vertical velocity is faster than when there is a low concentration of microstructures. Homotopy analysis of MHD free convection flow of micropolar fluid subjected stratified environment was examined by Koriko et al. [16].

Heat transfer has become an essential component of all industrial activities. As a result, it must be added, removed, or moved from one phase to the next. Heat transfer efficiency may be improved by increasing the thermal conductivity of the working fluid. Choi and Eastman [17] pioneered the concept of heat transfer enhancement when they developed an efficient approach to adjust heat transfer rate utilizing nanoparticles in their experiment. Nanofluids are employed in a wide range of engineering applications because of their enhanced heat transfer efficiency in diverse thermal systems. Nanofluids may find applications in engine coolant, automatic transmission fluid, brake fluid, gear lubrication, engine oil, and greases Senthilraja et al. [18]. Shah et al. [19] scrutinized analytic solution of Brownian motion and thermophoretic diffusion effects on Maxwell fluid containing tiny particles. It was reported that both the velocity and temperature distributes are enhanced for larger Brownian motion and thermophoretic parameter. Amar and Kishan [20] examined radiation effects on Mhd boundary layer flow of a nanofluid. Mhd stagnation point flow of a nanofluid with solar radiation effects was studied by Ghasemi and Hatami [21]. Rasool and Zhang [22] considered characteristics of chemical reaction and convective boundary conditions in Powell-Eyring nanoflid flow along a radiative Riga plate. It was reported that Brownian motion factor and thermophoresis improve the thermal boundary layer. Fayyadh et al. [23] discussed the influence of Biot number on the convective heat transfer of Darcy-Forchheimer nanofluid flow over stretched zero mass flux surface in the presence of magnetic field. Nayak et al. [24] examined partial slip and viscous dissipation effects on the radiative tangent hyperbolic nanofluid flow via a vertical permeable riga plate with internal heating. It was reported that the concentration of nanoparticles increases for larger thermophoretic force in the flow domain owing to the dragging process.

Motivated by all of the research on non-Newtonian fluid boundary layer flow past a vertical surface. The aim of this work is to investigate the bioconvection flow of micropolar fluid containing nanoparticles and gyrotactic microorganisms, with Catteneo-Christov model at the instances of strong and weak concentrations of microelements.

In accordance with boundary layer theory, a steady, buoyant convective boundary layer equation of micropolar nanofluid is investigated to achieve our goals. The flow is considered to be in the x-direction, which runs vertically along the plate, with y-direction being normal to it. Nanoparticles are considered to have little effect on the direction and speed with which microorganisms swim. Bioconvection flow is thought to occur exclusively in dilute nanoparticle suspensions. It is important to note that the base fluid is water in order for the microorganisms to survive. Under the foregoing assumptions with Boussinesq approximation and following the works of Koriko et al. [15], Koriko et al. [16], Hayat et al. [25] and Saleem et al. [26] the governing equations for mass, momentum, energy, and density of gyrotactic microorganisms in two-dimensional micropolar nanofluid may be expressed as follows;

Continuity Equation

∂ u ∂ x + ∂ v ∂ y = 0 , (1)

Momentum Equation

u ∂ u ∂ x + v ∂ u ∂ y = ( μ + χ ρ f ) ∂ 2 u ∂ y 2 + χ ρ f ∂ H ∂ y − σ B o 2 u ρ f + ( 1 − C ∞ ) ρ f g β ( T − T ∞ ) − g ( ρ p − ρ f ) ( C − C ∞ ) − g γ ( ρ m − ρ f ) ( N − N ∞ ) ρ f , (2)

Angular Momentum Equation

u ∂ H ∂ x + v ∂ H ∂ y = ( μ + χ / 2 ρ j ) ∂ 2 H ∂ y 2 − χ ρ j ( 2 H + ∂ u ∂ y ) , (3)

Energy equation of nanoparticles

u ∂ T ∂ x + v ∂ T ∂ y = κ ρ C p ∂ 2 T ∂ y 2 + τ ( D B ∂ T ∂ y ∂ C ∂ y + D T T ∞ ( ∂ T ∂ y ) 2 ) + λ 1 ( u 2 ∂ 2 T ∂ x 2 + v 2 ∂ 2 T ∂ y 2 + 2 u v ∂ 2 T ∂ y ∂ x + u ∂ u ∂ x ∂ T ∂ x + u ∂ v ∂ x ∂ T ∂ y + v ∂ u ∂ y ∂ T ∂ x + v ∂ v ∂ y ∂ T ∂ y ) , (4)

Concentration equation of nanoparticles

u ∂ C ∂ x + v ∂ C ∂ y = D B ∂ 2 C ∂ y 2 + D T T ∞ ∂ 2 T ∂ y 2 + λ 2 ( u 2 ∂ 2 C ∂ x 2 + v 2 ∂ 2 C ∂ y 2 + 2 u v ∂ 2 C ∂ y ∂ x + u ∂ u ∂ x ∂ C ∂ x + u ∂ v ∂ x ∂ C ∂ y + v ∂ u ∂ y ∂ C ∂ x + v ∂ v ∂ y ∂ C ∂ y ) , (5)

Density of gyrotactic microorganisms equation

u ∂ N ∂ x + v ∂ N ∂ y + b W c C w − C ∞ [ ∂ ∂ y ( N ∂ C ∂ y ) ] = D m ∂ 2 N ∂ y 2 (6)

The equations above are subjected to the following boundary conditions:

u = a x ,     v = 0 ,     H = − n ∂ u ∂ y ,     T = T w ,     C = C w ,     N = N w     at       y = 0 , (7)

u → 0 ,     H → 0 ,     T → T ∞ ,     C → C ∞ ,     N → N ∞     as       y → ∞ (8)

where u and v are velocity components in x and y directions respectively, c h is the chemotaxis constant, W c is the maximum cell swimming speed, m is the velocity power index, T f is the local fluid, T is the temperature of the fluid, C is the nanoparticle, N is the density of motile micro-organisms, p is the pressure, ρ f , ρ p , ρ m are the density of nanofluid, nanoparticles and micro-organisms respectively, D B , D T , D m are the Brownian diffusion coefficient, thermophoresis diffusion coefficient and diffusivity of micro-organisms respectively, κ , σ are the thermal and electrical conductivity of the fluid respectively, α is the thermal diffusivities, τ is the ratio of the effective heat capacitance of the nanoparticle to that of the base fluid, j is the micro-inertial density, χ is the vortex viscosity, n is the surface condition parameter and varies from 0 to 1. The study considers two cases when n = 0 and n = 1 / 2 . When n = 0 (called strong concentration) indicates concentrated particle flows in which the microelements close to the wall are unable to spin, H = 0 near the wall and H is the micro-rotation or angular velocity whose direction of rotation is in the xy plane. When n = 1 / 2 , the anti-symmetric component of the stress tensor vanishes, indicating a weak concentration of microelements, when n = 1 , the situation is utilized for simulating turbulent boundary layer flows. The micropolar parameter, often known as the material

parameter, is K = χ μ , K ≠ 0 for micropolar fluid and K = 0 for classical Newtonian fluid.

For the sake of non-dimensionalization and parameterization of Equations (2)-(5) and (6) subject to boundary conditions (7) and (8), the following similarity variables are considered corresponding to the following models:

η = a 1 2 ϑ 1 2 y ,         ψ ( x , y ) x a 1 2 ϑ 1 2 = f ( η ) ,         H = a x ( a ϑ ) 1 / 2 , θ ( η ) = T − T ∞ T w − T ∞ ,       ϕ = C − C ∞ C w − C ∞ ,       w ( η ) = N − N ∞ N w − N ∞ , (9)

when the stream function ψ ( x , y ) is introduced the continuity Equation (1) is automatically satisfied then Equations (2)-(6) become

( 1 + K ) d 3 f d η 3 − ( d f d η ) 2 + f d 2 f d η 2 + K d p d η − M d f d η + G r θ − N r ϕ + R b w (10)

( 1 + K 2 ) d 2 p d η 2 + f d p d η − p d f d η − K ( 2 p + d 2 f d η 2 ) = 0 (11)

d 2 θ d η 2 + P r f d θ d η − P r δ t ( f 2 d 2 θ d η 2 + f d f d η d θ d η ) + P r N b d ϕ d η d θ d η + P r N t ( d θ d η ) 2 = 0 (12)

d 2 ϕ d η 2 + L e f d ϕ d η − L e δ n ( f 2 d 2 ϕ d η 2 + f d f d η d ϕ d η ) + N t N b d 2 θ d η 2 = 0 (13)

d 2 w d η 2 + S c m f d w d η − P e ς d 2 ϕ d η 2 − P e w d 2 ϕ d η 2 − P e d ϕ d η d w d η = 0 (14)

subject to the boundary conditions (7) and (8) which becomes;

f ( 0 ) = 0 ,     d f ( 0 ) d η = 1 ,     p ( 0 ) = − n d 2 f d η 2 ,     θ ( 0 ) = 1 ,     ϕ ( 0 ) = 1 ,     w ( 0 ) = 1 , (15)

d f ( ∞ ) d η = 0 ,     p ( ∞ ) = 0 ,     θ ( ∞ ) = 0 ,     ϕ ( ∞ ) = 0 ,     w ( ∞ ) = 0 (16)

In the above equations primes denote differentiation with respect to η . The dimensionless velocity, temperature, concentration and density of microorganisms

are represented as f ( η ) , p ( η ) , θ ( η ) , ϕ ( η ) and w ( η ) . P r = ϑ α is the Prandtl number, M = σ B o 2 ρ f a is the magnetic parameter, G r = g ρ f ∞ β ( 1 − C ∞ ) ( T f − T ∞ ) ρ f a 2 x is the modified local Grashof number, N r = g ( C w − C ∞ ) ( p ρ − p f ) ρ f a 2 x is the buoyancy-ratio parameter, R b = g γ ( C w − C ∞ ) ( N − N ∞ ) ρ f a 2 x is the bioconvection Rayleigh number, N b = ( ρ C ) p ϑ ( ρ C ) f D B ( C w − C ∞ ) is the Brownian motion parameter, N t = ( ρ C ) p ϑ ( ρ C ) f D T T ∞ ( T w − T ∞ ) is the thermophoresis parameter, δ t = λ 1 a is the relaxation time parameter of temperature, L e = ϑ D B is the Lewis number, δ n = λ 2 a is the relaxation time parameter of nanoparticle volume fraction, ς = N ∞ N w − N ∞ is gyrotactic microorganisms concentration difference parameter, S c m = ϑ D B is the Schmidth number for diffusing motile microorganisms and P e = b W c D n is the Peclet number.

The dimensionless form of skin friction coefficient,

R e 1 / 2 C f x = f ″ ( 0 ) (17)

The dimensionless form of wall motile microorganisms number,

R e 1 / 2 Q n x = − w ′ ( 0 ) (18)

where R e = a x 2 ϑ is the local Reynolds number.

HAM provides a great freedom to express solutions of a given non-linear problem by means of different base functions. Non-linear problems can be approximated more efficiently by choosing a proper set of base functions, mainly because, the convergence region and rate of a series are determined by the base functions used to express the solution [27] [28] [29] [30] [31]. In view of the boundary conditions (15) & (16), f ( η ) , p ( η ) , θ ( η ) , ϕ ( η ) and w ( η ) can be expressed by the set of base functions in the form

〈 η j exp ( − n j ) | j ≥ 0 , n ≥ 0 〉 (19)

The solutions f ( η ) and θ ( η ) can be represented in a series form as

f ( η ) = a 0 , 0 0 + ∑ n = 0 ∞     ∑ k = 0 ∞     a n , k k η k exp ( − n j ) (20)

p ( η ) = ∑ n = 0 ∞     ∑ k = 0 ∞     b n , k k η k exp ( − n j ) (21)

θ ( η ) = ∑ n = 0 ∞     ∑ k = 0 ∞     c n , k k η k exp ( − n j ) (22)

ϕ ( η ) = ∑ n = 0 ∞     ∑ k = 0 ∞     d n , k k η k exp ( − n j ) (23)

w ( η ) = ∑ n = 0 ∞     ∑ k = 0 ∞     e n , k k η k exp ( − n j ) (24)

in which a n , k k , b n , k k , c n , k k , d n , k k and e n , k k are the coefficients. As long as such a set of base functions are determined, the auxiliary function H ( η ) , the initial approximation f o ( η ) , p o ( η ) , θ o ( η ) , ϕ o ( η ) and w o ( η ) and the auxiliary linear operator L f , L p , L θ , L ϕ and L w must be chosen in such a way that all solutions exist and can be expressed by these sets of base functions. Therefore, in framing the Homotopy Analysis Method (HAM), we apply the rule of solution expressions in choosing the auxilliary function H ( η ) , initial approximation f o ( η ) , p o ( η ) , θ o ( η ) , ϕ o ( η ) and w o ( η ) .

Invoking the rule of solution expressions for f ( η ) , p ( η ) , θ ( η ) , ϕ ( η ) and w ( η ) for Equations (15) & (17) together with (18)-(21) are expressed as

f o ( η ) = 1 − exp ( − η ) ,   p o ( η ) = − n ∗ exp ( − η ) ,   θ o ( η ) = exp ( − η ) , ϕ o ( η ) = exp ( − η ) ,   w o ( η ) = exp ( − η ) (25)

with linear operators L f , L θ , L ϕ and L w as

L f [ f ( η ; q ) ] = ∂ 3 f ( η ; q ) ∂ η 3 − ∂ f ( η ; q ) ∂ η (26)

L p [ p ( η ; q ) ] = ∂ 2 p ( η ; q ) ∂ η 2 − p ( η ; q ) (27)

L θ [ θ ( η ; q ) ] = ∂ 2 θ ( η ; q ) ∂ η 2 − θ ( η ; q ) (28)

L ϕ [ ϕ ( η ; q ) ] = ∂ 2 ϕ ( η ; q ) ∂ η 2 − ϕ ( η ; q ) (29)

L w [ w ( η ; q ) ] = ∂ 2 w ( η ; q ) ∂ η 2 − w ( η ; q ) (30)

The linear operators L f , L p , L θ , L ϕ and L w have the following properties

L f [ C 1 + C 2 exp ( − η ) + C 3 exp ( η ) ] = 0 , L p [ C 4 exp ( − η ) + C 5 ] = 0 , L θ [ C 6 exp ( − η ) + C 7 ] = 0 , L ϕ [ C 8 exp ( − η ) + C 9 ] = 0 , L w [ C 10 exp ( − η ) + C 11 ] = 0 (31)

in which C 1 , C 2 , C 3 , C 4 , C 5 , C 6 , C 7 , C 8 , C 9 , C 10 , and C 11 are constants.

The flow of micropolar fluid in a water-based solution containing both nanoparticles and gyrotactic microorganisms has been analyzed for various values of emerging parameters. The effects of various parameters on velocity f ′ ( η ) , microrotation p ( η ) , temperature θ ( η ) , nanoparticles concentration ϕ ( η ) , and motile microorganism density profiles w ( η ) have been studied and reported. The theoretical values of the governing parameters for both cases of weak and strong concentration of microelements are M = G r = N r = R b = 0.5 , δ t = δ n = N b = N t = ς = 0.1 , S c = P e = 1.0 , P r = 1.2 , K = 0.4 , L e = 2.0 . Table 1 displays the numerical values of the local skin friction coefficient, plate couple stress, local Nusselt number, local Sherwood number, and local density number for different values of the Micropolar parameter K, while the other parameters are held constant at strong microstructure concentrations. The local skin friction coefficient is found to be a decreasing function of K, with the first two values increasing plate couple stress and the final two entries lowering plate couple stress. The first two K entries show an increase in the local Nusselt number, whereas the latter two show a reduction. The local Nusselt number is a decreasing function of K, and the first two entries of K improve the local density number, while the latter two entries of K slightly improve it.

Table 2 displays the numerical values of the local skin friction coefficient, plate couple stress, local Nusselt number, local Sherwood number an local density number for various values of Micropolar parameter K when the other parameters are held constant at low microstructure concentrations. The first two K values indicate a decrease in the local skin friction coefficient, whereas the latter two entries indicate an increase in the local skin friction coefficient. The first and last two entries of K represent an increase and decrease in plate couple stress respectively. It is also worth noting that the local Nusselt number and Sherwood number are both decreasing functions of K, whereas the first two and final two entries of K cause an increase in the local density number. Figures 1(a)-5(b) show the effect of micropolar parameter K on various profiles when n = 0 & 0.5 . Figure 1(a) and Figure 1(b) show that the velocity profiles increase as the magnitude of the micropolar parameter K augment when n = 0 & 0.5 . In Figure 2(a) it is noticed that larger values of micropolar parameter K correspond to an increase in the microrotation profile when n = 0 . It is also demonstrated that the

particle near the wall does not rotate because n = 0 depicts concentrated particle flows in which micro-elements close to the wall are unable to rotate.

However as shown in Figure 2(b) when the micropolar parameter K grows, the microrotation profile diminishes. Figure 3(a) indicates that increasing the value of K yields a reduction in the temperature distribution when n = 0 but Figure 3(b) shows a more massive decrease in the temperature distribution as the value of K increases. Figure 4(a) shows that increasing values of K leads to decrease in the concentration of nanoparticles distribution when n = 0 whereas increasing K results in increase in the concentration of nanoparticles within

the domain 0 ≤ η ≤ 2.56 and thereafter all profiles merged together and later decline towards the freestream. Figure 5(a) and Figure 5(b) depict the influence of K on the density of motile microorganisms profiles when n = 0 and n = 0.5 . Figure 5(a) indicates that increasing K results in a slight decline in the density of motile microorganisms profiles when n = 0 but increases the density of motile microorganisms when n = 0.5 as revealed in Figure 5(b). The observation in Figure 5(b) is completely different from that of Figure 5(a) mainly because microorganisms are free to rotate at the wall allowing the profiles to enhance particularly at the bottom layer of the wall.

Figure 6(a) and Figure 6(b) depict the effects of M on the microrotation profiles when n = 0 and n = 0.5 respectively. In Figure 6(a) it noticed that when n = 0 all profiles fused together at the peak 0.26 and afterwards fall for larger values of M towards the freestream, whereas in Figure 6(b) it is observed that there is a conspicuous augmentation in microrotation profile when n = 0.5 as M increases.

In Figure 7(a) when n = 0 it is revealed that raising the value of N_t produces an increase in the microrotation distribution along the wall. Furthermore, all profiles merged at η = 1.6 , resulting in decrease in the microrotation profile. Figure 7(b) shows that when n = 0.5 , incremental values of N_t cause a slight increase in microrotation profile, but for larger values of N_t, the profiles decline towards the freestream. Figure 8(a) and Figure 8(b) demonstrate impacts of N_t on temperature profiles when n = 0 and n = 0.5 . Figure 8(a) illustrates that increasing the magnitude of N_t improves the temperature distribution when n = 0 . Figure 8(b) reveals that when n = 0.5 , there is an increase in the temperature profile within the region 0 ≤ η ≤ 2.89 as N_t increases.

Figures 9(a)-10(b) depict the effects of the gyrotactic microorganisms concentration difference parameter ς for various profiles when n = 0 & 0.5 . In

Figure 9(a) and Figure 9(b) when n = 0 & 0.5 , incremental values of ς leads to temperature distribution augmentation. Physically, increasing the temperature of the fluid causes the bonds that hold the molecules of particles together to break down, allowing the microorganisms swimming within the fluid to interact and swim back to bottom layer of the fluid, allowing bioconvection to occur again. Figure 10(a) and Figure 10(b) indicate that increasing ς cause a noticeable reduction in the density of motile microorganisms profile for both n = 0 & 0.5 .

In this communication, the influence of pertinent micropolar fluid flow parameters on nanoparticles and gyrotactic microorganisms is investigated. Both instances of strong and weak microelement concentrations are taken into consideration. The controlling flow equations were modeled and converted using suitable similarity variables. The series solution was determined using the Homotopy Analysis Method (HAM). The research produced the following important findings:

1) Incremental values of micropolar parameter K signifies enhancement of velocity profiles when n = 0 & n = 0.5 , while the microrotation profiles increase when n = 0 and declines when n = 0.5 .

2) Increasing values of micropolar parameter K corresponds to diminishing of temperature profiles when n = 0 & n = 0.5 .

3) Incremental values of magnetic parameter M corresponds to decline of microrotation profile when n = 0 while augmentation of microrotation is observed when n = 0.5 .

4) Temperature profiles are increasing functions of gyrotactic microorganisms concentration difference parameter ς for both n = 0 & n = 0.5 , while density of motile microorganisms profiles are decreasing functions of ς for both n = 0 & n = 0.5 .

The corresponding author “Dr. Tosin Oreyeni” would like to express his sincere appreciations to the respected reviewers for their technical and insightful comments that helped to enrich meaningfully the present investigation.

The authors declare no conflicts of interest regarding the publication of this paper.

Koriko, O.K., Omowaye, A.J., Popoola, A.O., Oreyeni, T., Adegbite, A.A., Oni, E.A., Omokhuale, E. and Altine, M.M. (2022) Insight into Dynamics of Hydromagnetic Flow of Micropolar Fluid Containing Nanoparticles and Gyrotactic Microorganisms at Weak and Strong Concentrations of Microelements: Homotopy Analysis Method. American Journal of Computational Mathematics, 12, 267-282. https://doi.org/10.4236/ajcm.2022.122017