Applications of bioconvection for tiny particles due to two concentric cylinders when role of Lorentz force is significant.

The bioconvection flow of tiny fluid conveying the nanoparticles has been investigated between two concentric cylinders. The contribution of Lorenz force is also focused to inspect the bioconvection thermal transport of tiny particles. The tiny particles are assumed to flow between two concentric cylinders of different radii. The first cylinder remains at rest while flow is induced due to second cylinder which rotates with uniform velocity. Furthermore, the movement of tiny particles follows the principle of thermophoresis and Brownian motion as a part of thermal and mass gradient. Similarly, the gyro-tactic microorganisms swim in the nanofluid as a response to the density gradient and constitute bio-convection. The problem is modeled by using the certain laws. The numerical outcomes are computed by using RKF -45 method. The graphical simulations are performed for flow parameters with specific range like 1≤Re≤5, 1≤Ha≤5, 0.5≤Nt≤2.5, 1≤Nb≤3, 0.2≤Sc≤1.8, 0.2≤Pe≤1.0 and 0.2≤Ω≤1.0. It is observed that the flow velocity decreases with the increase in the Hartmann number that signifies the magnetic field. This outcome indicates that the flow velocity can be controlled externally through the magnetic field. Also, the increase in the Schmidt numbers increases the nanoparticle concentration and the motile density.

1. Introduction

A nanofluid is a novel class of fluids in which metallic or nonmetallic nanoparticles are scattered over the base fluid. These nanofluids are known for their special heat transfer properties like the presence of nanoparticles pronounced the reflective thermal aspect of base materials [1]. This special property makes these fluids applicable for several engineering applications such as coolants in automobiles, nuclear reactors, solar heaters etc. Some of the special properties of nanoparticles are high specific surface area, higher dispersion stability, reduced particle clogging, and many adjustable properties including thermal conductivity and surface wettability [2]. In this regard, Sheikholeslami [3] studied the impact of porosity and Lorentz force on the heat transfer of nanofluid using Darcy law. Bhamani et al. [4] made contributions for growing heating aspect for turbulent flow of nanofluid inside a pipe. Alrashed et al. [5] modelled a system that describes the flow and the thermal performance of nanofluid comprising of water and MWCNT. Abbas et al. [6] scrutinized the entropy generation in the fully developed flow of nanofluid subjected to velocity slip. Nadeem et al. [7] analyzed the heat transfer characteristics of hybrid nanofluid flowing over a curved surface. The nanomaterials various improved migrated phenomenon with subclass of base fluids was directed in communications Puneeth et al. [8, 9]. These continuations on the analysis of heat transfer of single-phase and double phase nanofluid were extended to ternary nanofluid by Manjunatha et al. [10] which comprises of a base fluid and three different classes of nanoparticles. Song et al. [11] reported the thermal distribution of alumina and copper nanomaterials with ethylene glycol and water base fluid. Oke et al. [12] expressed the role of Coriolis force to express the thermal experience of alumina nanoparticles with diameter 47nm. Animasaun et al. [13] presented the thermo-haphazard prospective of nanoparticles with diverse thermal properties. The inclusion of tiny nanoparticles in different kind of base liquids with meta investigations was directed by Wakif et al. [14].

Pattern-forming convection movements established in suspensions of paddling microorganisms is known as bio-convection. In liquid suspensions of floating microorganisms, the cellular streaming trend was found wherein fluid flow motions proceed downwards in places where elevated levels of microorganisms develop and swim upwards in regions of low concentration. This sort of pattern is determined by factors including the depth of the suspension, as well as the quantity and mobility of the organisms. Khan et al. [15] used the Cattaneo-Christov model to analyze the double diffusion in the nanofluid flow containing motile cells. This was further followed by Yahya et al. [16] to explore the bio-convection in Williamson nanofluid. Similarly, Puneeth et al. [17] considered the bio-convective flow of Williamson nanofluid past a Riga plate. Shi et al. [18] analyzed the impact of activation energy and gyro-tactic microorganisms in enhancing the thermal performance of magneto-cross nanofluid. Waqas et al. [19] performed simulations using numerical methods to analyze the impact of microorganisms swimming in a non-Newtonian nanofluid subject to the magnetic field. Further, Koriko et al. [20] analyzed the process of bio-convection comprising the implosion of microorganisms in a thixotropic nanofluid. Puneeth et al. [21] discussed the homogeneous and heterogeneous chemical reactions with the quartic autocatalysis for the flow of micropolar nanofluid flowing in a channel subjected to thermal radiation for microorganism interference. Azam et al. [22] evaluated the impact of nonlinear radiation in the flow of nanofluid under the influence of motile cells. Balla et al. [23] encountered the chemical reaction aspect while reporting the bioconvection phenomenon for oxytactic microorganisms. Makinde and Animasaun [24] performed the determination of bioconvection pattern subject to autocatalysis reactive species in the upper regime of paraboloid. The fluctuation in bioconvective properties of nanofluids with influence of nonlinear radiative phenomenon in paraboloid revolution was visualized in the work of Makinde and Animasaun [25]. Khan et al. [26] rolled out the impact of Navier slip for nanofluids flow with microorganisms. Khan and Makinde [27] conveyed the heat transfer improvement for nanofluids flow with bioconvection enrolment.

The motion of suspended matter across a fluid in the presence of a temperature gradient is known as thermophoresis. A study to interpret the role of thermophoresis on particle migration and concentration distribution discovered that the concentration distribution gets more non-uniform as the particle size increases. Meanwhile, thermophoresis accentuates non-uniformity in the concentration distribution, with a stronger effect at higher mean concentrations [28]. Many researchers have explored the flow with thermophoresis and Brownian motion as they play an important role in the analysis of heat and mass transfer. For instance, Sheikholeslami et al. [29] studied the impact of magnetic force for tiny particles imposed two cylinders having circular orientation. Mirzaeyan and Toghraie [30] investigated the laminar flow of nanofluid between porous cylinders. Arif et al. [31] justified importance of GO nanoparticle in the MoS2−H2O nanofluid for enhancing its thermal performance. Abbas et al. [32] analyzed the impact of magnetic field on the velocity of nanofluid flowing past a non-linear stretching sheet. Reddy et al. [33] employed the Cattaneo-Christov heat flux model to examine the thermal characteristics of nanofluid flowing past a swirling cylinder. Biswal et al. [34] deterimned numerical simulation to investigate the flow of a nanofluid in a semi-porous channel subjected to the magnetic field. Khan et al. [35] discussed the impact of Joule heating referring to interaction of nano-compounds past a swirling cylinder under the influence of Lorentz force. Aghamiri et al. [36] designed a mathematical model that describes the impact of forced convection on the flow of Ferro-nanofluid flowing in a microchannel consisting of rotating cylinders. Bouzerzour et al. [37] discussed natural convection in a nanofluid flowing in an annular space formed due to the separation between confocal elliptic cylinders at different geometric positions. Ch et al. [38] claimed thermophoresis inspection of Walter’s B nanofluid subject to interaction of buoyancy forces. Khan and Ali [39] worked out the thermophoresis model based on nanofluid properties by using Eyring-Powell fluid model. The thermal statement and role of thixotropic nanomaterials incorporating the thermo-diffusion phenomenon for Riga configuration was addressed by Khan et al. [40].

After presenting a comprehensive literature survey, it has been noticed that bioconvection aspect of nanofluid with various flow configurations have been available. However, the bioconvection applications of tiny fluid conveying the nanoparticles between two concentric cylinders different radii is not focused yet. Moreover, the contribution of Lorentz force for bio-convective model is another important task which is addressed in this model. The flow though moving cylinder is interesting topic and some continuations are performed by researchers [41-50]. This investigation presents the answer of following thermal flow questions:

Which mathematical model is used to inspect the bioconvection of tiny particles moving between concentric cylinders having different radii?

How heat and mass transfer process fluctuated with interaction of tiny fluid conveying the nanoparticles?

What is contribution of Lorentz force to improve the heating phenomenon?

How thermophoresis and Brownian motion parameters pays role to enhance the thermal process?

2. Mathematical model

A laminar flow of a Newtonian fluid containing the nanoparticles in the presence of gyro-tactic microorganisms is assumed to flow between two concentric cylinders. Each of these cylinders are of radii R1 and R2 such that R1>R2. The cylinder with radius R2 is assumed to be stationary and is enclosed in the cylinder of radius R1, whose angular momentum is ω1. The presence of microorganisms in the system stabilizes the dilute nanoparticle suspension and prevent sedimentation in the system. The dilute nanoparticle suspension is assumed which no fluctuation of movement of microorganisms within the system. This enables the microorganisms to move easily and constitute the macroscopic phenomena termed bio-convection. Thermophoresis and Brownian motion prospective of tiny particles is examined in view of Buongiorno thermal model. Hence, these two slip mechanisms are included in the mathematical model so that the results that are obtained will be close to practicality. Furthermore, the cylinders of radii r1 and r2 are maintained at a temperature T1 and T2 respectively along with the concentration C1 and C2 and the motile density N1 and N2 respectively. The flow configuration is shown in Fig 1 using cylindrical coordinates.

Fig 1

Schematic flow diagram.

The exploration of thermal model for all constraints is presented via following equations [36, 37, 46]: with

The non-dimensionless form of model is [36, 37, 46]:

The corresponding boundary conditions are

The dimensionless parameters involved in this study are defined as

3. Solution methodology

The transformed Eqs (6)–(9) along with the boundary conditions (10) are remodeled to initial value problem (IVP). This is further simulated with interpretation of well-known RKF -45 method in acquaintance with the shooting method. For the computation purpose, the infinite boundary conditions are considered at η = 10 and the accuracy of the solution is set to the order of 10−5. The proper step size is determined in this method for ensuring the validity. Further, these two approximations are compared and if they hold a close agreement with each other than the approximation is considered valid. The whole process is repeated if the approximations obtained do not match each other and the computation is repeated till the desired accuracy is obtained.

4. Results and discussion

RKF -45 is used to examine the flow of nanofluid between two concentric cylinders in the presence of self-propelled microorganisms. Using appropriate relationships, the equations are non-dimensionalised. The solution was achieved using RKF -45. The outcomes of this study are interpreted through graphs ((2a)-(5d)) and Tables 1 and 2.

Table 1

Variations of quantities of physical interest for the changes in Re and Ha.

Parameter	Range	−(dudr)_r=0	−(dθdr)_r=0	−(dϕdr)_r=0	−(dXdr)_r=0
Re	1	6.019791941	0.114981271	2.455937464	3.294764449
	2	5.599051794	0.033347737	1.570841211	1.789029468
	3	5.188232292	0.008195391	0.928568857	0.800785493
	4	4.789171985	0.001681899	0.503667618	0.252470081
	5	4.404083642	0.000282309	0.248382388	0.013075776
Ha	1	5.260641448	0.098887375	2.350184677	3.120977258
	2	6.019791941	0.114981271	2.455937464	3.294764449
	3	7.010034918	0.134403776	2.572311573	3.487258489
	4	8.098824067	0.153031691	2.675925633	3.660256028
	5	9.234128219	0.169444086	2.762766988	3.806829769

Table 2

Variations of quantities of physical interest for the changes in Sc, Nb and Nt.

Parameter	Range	−(dθdr)_r=0	−(dϕdr)_r=0	−(dXdr)_r=0
Sc	1	0.106610988	3.185041421	3.742403689
	2	0.110749611	2.801801917	3.507260711
	3	0.114981271	2.455937464	3.294764449
	4	0.119287931	2.145470405	3.103745486
	5	0.123651544	1.868213341	2.932904232
Nb	1	0.482348184	3.112881171	3.700386708
	2	0.241850059	2.630081554	3.402457446
	3	0.114981271	2.455937464	3.294764449
	4	0.052410192	2.364185944	3.237914973
	5	0.023120414	2.307734571	3.202886897
Nt	0.5	0.114981271	2.455937464	3.294764449
	1.0	0.077230578	2.824710126	3.522508617
	1.5	0.051498855	3.191753058	3.748817773
	2.0	0.034117987	3.555669601	3.972797117
	2.5	0.022471965	3.916123096	4.194239191

Fig 2(A)–2(D) shows the impact of Reynolds number (Re) on the nanofluid flow profiles. Following the physical dynamic of Reynolds number, inertial forces grow up when Reynolds number is higher. Such forces show their major impact at the boundary region. As the value of Re goes higher, the viscous force becomes less significant, and the fluid will thus be less viscous and result in faster flow. The up-raise change in velocity due to Re is noted in Fig 2(A) Moreover, the temperature rate of tiny particles increases due to the friction created within the fluid due to its faster flowing rate as shown in Fig 2(B). Similarly, the Fig 2(C) and 2(D) indicated that improvement assigning to Reynolds number enhances the mass concentration and the motile density respectively.

Fig 2

The impact of Reynolds number on the fluid profiles.

The Hartmann number (Ha) is the ratio of electromagnetic forces to the viscous forces which measures the significance of drag forces resulting from electromagnetic induction and viscous forces. The impact of this parameter is shown in Fig 3(A)–3(D). The velocity of the nanofluid flow is seen to be reduced in Fig 3(A) because of the strong Lorentz force produced. This force acts against the flow and opposes the motion of the fluid by creating friction. Also, the friction thus created will generate additional heat within the nanofluid as a result more impressive temperature field is noted (Fig 3(B)) for higher values of Ha. Due to the slow movement of nanofluid, the nanoparticles and microorganisms accumulate at the boundary layer. The improved change in concentration and microorganism profile is noticed in Fig 3(C) and 3(D).

Fig 3

The impact of Hartmann number on the fluid profiles.

The impact of thermophoresis (Nt) on θ(r) and ϕ(r) is shown in Fig 4(A) and 4(B) respectively. The increase in the Nt parameter causes the nanoparticles to move from a hotter region to a colder and the nanoparticles dissipate heat into the fluid. As a consequence, the temperature of the nanofluid increases as shown in Fig 4(A). Meanwhile, the movement of nanoparticles becomes faster with the increase in Nt as a result, ϕ(r) at the boundary layer gowsup as reflected in Fig 4(B). Further, the temperature of the nanofluid increases due to the heat generated because of the collision of nanoparticles. Thus increment with increasing trend in referred to Fig 4(C) the zigzag motion of nanoparticles increases with the increase in Nb. During this zigzag motion, the nanoparticles colloid each other and move away from the boundary region as a result the nanoparticle concentration decreases as depicted in Fig 4(D).

Fig 4

The impact of Nt and Nb on the thermal and concentration profiles.

The impact of the Schmidt numbers on the Φ(r) and X(r) is depicted in Fig 5(A) and 5(B) respectively. The Schmidt numbers are inversely proportional to the diffusivities of their corresponding profiles. As a result, as the concentration Schmidt number (Sc) increases, the diffusivity of the nanoparticles reduces, and Φ(r) at the boundary layer falls, as illustrated in Fig 5(A). As Sb increases, the motile density diffusivity falls, and the motile density at the border layer drops, as seen in Fig 5(B). Furthermore, as illustrated in Fig 5(C), increasing the Peclet number increases motile density at the border layer. When Ω upgrade, the microorganism profile declined as shown in Fig 5(D).

Fig 5

The impact of Sc, Pe and Ω on Φ(r) and X(r).

The variations of wall shear surface force, local Nusselt number, Sherwood number and motile density number is tabulated in Tables 1 and 2 for changes in the fluid parameters. Table 1 displays the impact of Reynolds number and Hartmann number whereas, Table 2 shows the impact of Schmidt number and slip mechanisms. The increase in the Reynolds number resulted in a decrease in Cf Nu, Sh and Nn and the same is tabulated in the first row of Table 1. Whereas, it was noticed that the higher values of Ha increased Cf, Nu, Sh and Nn. Meanwhile, the increase in the Schmidt number increased the Nusselt number whereas it decreased the Sherwood number and motile density number. Furthermore, the increase in the values of Nb, decreased the values of Nu, Sh and Nn, but the higher values of Nt decreased the Nusselt number and enhanced Sherwood number and motile density number.

5. Conclusions

The applications of Lorentz force for bioconvection transport of tiny fluid conveying tiny particles due to two concentric cylinders is presented. For nanofluid flow, the concentric cylinders attained same center but different radius. The governing equations are made dimensionless using appropriate relationships for which simulations have been performed via RKF-45 method. Some of the major outcomes of the study are:

➢ The increment in Reynolds number enhanced the nanoparticle concentration and motile density profiles.

➢ The Hartmann number is dominant over the fluid flow and the higher values of Ha reduce the fluid flow.

➢ The enhanced flow velocity and heating phenomenon is noted for increasing the Reynolds number.

➢ The increase in the thermophoresis parameter enhances the thermal and mass profiles of the nanofluid.

➢ The greater values of Schmidt numbers enhance their corresponding profiles whereas the motile density decreases for higher Peclet numbers.