Darcy resistance flow of Sutterby nanofluid with microorganisms with applications of nano-biofuel cells.

The objective of current research is to endorse the thermal aspect of Sutterby nanofluid containing the microorganisms due the stretched cylinder. The features of nonlinear thermal radiation, Darcy resistance and activation energy are also incorporated to inspect the thermal prospective. The problem is further extended with implementation of modified Fourier and Fick's theories. The results are presented for the stretched cylinder and also for stationary plate. The numerical formulation for the problem is presented by following the shooting technique. The comparative numerical is performed to verify the computed simulations. The results convey that the presence of Darcy resistance parameter enhanced the velocity more effectively for stretched cylinder. A reduction in velocity due to Sutterby fluid parameter and buoyancy ratio parameter has been observed. Moreover, the temperature profile enhanced with larger sponginess parameter more effectively for stretching cylinder.

Introduction

In today's growing world of technology, most engineers, scientists and researchers are focused on the analysis of nanoliquid due to their physical applications in the field of applied sciences. Nanofluids are known to be the accumulation of micro-shape solid objects in the convectional fluid. The accumulation of micro-shaped solid objects in the convection fluid is considered to be nanofluids. The enhancement in the heating capacity is necessary in different range of industrial and engineering mechanisms. The availability of increasing heat source is necessary to improve the optimal growth of industrial products. The communication and growing research in nanofluids predict novel thermal applications of such particle in the energy sector, power plants, chemical processes, engineering devices, thermal mechanisms etc. Choi[1] explored the novel thermal aspect of nanofluids in primary research via experimental support. Alharbi[2] studied flow of hybrid nanofluids with heat energy impacts. Uddinet al.[3] scrutinized the radiative slip impact in nanofluid. Hassan et al.[4] focused the thermal availably of hybrid nanofluid by assuming the shear thinning as a base liquid. Tlili et al.[5] premeditated the thermal flow of Oldroyd-B nanofluid with isotropic slip impact. Khan et al.[6] analyzed Darcy-Forchheimer flow in hybrid nanofluid. Haq et al.[7] focused the improve thermal properties of nanofluid with suspension of Casson liquid. Xia et al.[8] addressed the natural convective optimized analysis for Eyring-Powell nanofluid subject to microorganisms. The cross nanofluid flow with entropy generation assessment was utilized by Haq et al.[9]. In another reports, Haq et al.[10] inspected the bioconvection applications for the nanofluid flow with controlled optimized phenomenon. Hussain et al.[11] observed the carbon nanotube thermal outcomes with fluctuation of dynamic viscosity. The melting applications for hybrid nanofluid in addition of variable viscosity has been intended by Hussain et al.[12]. The analysis of Hussain et al.[13] reports the improved heat transfer analysis for Jeffrey material with external heat source. Some more recent research on nanofluids is referred to refs.[14-17].

The Sutterby nanofluid study is another topic of interest to enhance the thermal features of base fluids. Mir et al.[18] addressed the relative improvement of heat transportation phenomenon by endorsing the Sutterby nanofluid in base liquid. The thermal research via interaction of Sutterby nanofluid with progressive thermal change was explored by Nawaz et al.[19]. Bilal et al.[20] observed the role of magnetic force for Sutterby nanofluid flow. The influence of thermal radiation and the inclined magnetic field on the Sutterby fluid by focusing on the Cattaneo-Christov heat flux structure is examined by Sabir et al.[21]. Song et al.[22] performed the Marangoni convection analysis for Sutterby nanofluid with melting and solutal constraints. Abbasi et al.[23] endorsed the Sutterby nanoparticles properties in trapped channel.

The self-induced motile micro-organisms can increase the density of ordinary fluids in specific direction and, as a response a bioconvection phenomenon has evolved. Such nanoparticles motion is independent of microorganism movement and therefore the collective functionality of bioconvection and nanofluids seems to also be vital for microfluidics devices. Waqas et al.[24] reported the bio-convective model for generalized viscoelastic nanofluid by performing the numerical simulations. Aziz et al.[25] tackled a bidirectional bioconvection thermal nanofluid problem subject to accelerating space. Khan et al.[26] discovered the theoretical continuation for couple stress nanofluid with consequences of activation energy and bioconvection phenomenon. Tong et al.[27] explored the suspension of microorganisms subject to slip implementation for nanofluid. Li et al.[28] reported the bio-convection applications for modified second grade fluid. Alwatban et al.[29] expressed the physical onset of bioconvection phenomenon regarding the nanofluid flow.

Although a lot of research have been performed on the nanofluids, however, the thermal applications of Sutterby nanofluid subject to the bioconvection phenomenon with diverse flow features has not been performed yet. One this end, this research presents the thermal flow of Sutterby nanofluid containing the microorganisms due to moving cylinder. The navel aspects of this model are:

The Darcy resistance flow of radiated Sutterby nanofluid with microorganisms due moving cylinder is presented.

The Fourier and Fick’s modified expressions are used for examining the heat and mass transfer phenomenon.

The novel thermal features like nonlinear thermal radiation and activation energy are also incorporated.

The convective boundary conditions are utilized with motivations of enhancing the thermal transport of Sutterby nanofluid.

The shooting technique for the formulated boundary value problem is implemented for obtained numerical simulations.

The obtained simulations may present novel significances in bio-fuels, enzymes, thermal processes, energy systems, heat transfer devices etc.

Flow model

The thermal transport of Sutterby nanofluid with consideration of suspension of microorganisms is taken into consideration. The stretched cylinder is assumed to originate the laminar flow. The velocity of moving cylinder is attributed to be The normal aspect of magnetic force is also utilized[30]. The consequence of activation energy for solutal transport is discussed. A physical schematic of flow model is depicted in Fig. 1. Moreover , and signifies wall surface temperature, concentration and microorganisms respectively. Here ambient temperature, concentration and microorganisms are symbolized by , and correspondingly.

Figure 1

Physical view of flow.

Under the above postulations, the leading governing boundary layer equations of flow are given as follows[18-20]:

Here,wherewith boundary conditions:

In the above Eq. (2) Darcy resistance is defined as[31]:

Introducing the following suitable similarities variables[18-20]:

After introducing the above appropriate transformation (10) in governing PDE’s, we acquireWithwith dimensionless parameter:

Figure 10

Profile of for .

Physical quantities of interest are

the local heat flux, for local mass flux and microorganisms flux, which are addressed as

In the dimensionless forms are

Here the local Reynolds number is symbolized by .

Numerical procedure

The coupled governing odrinary differential Eqs. (11–14) with boundary restrictions (15) is higly nonlinear in nature. The numerical solutions of these system is very diffecullt. There for obtain numerical solutions of model we employ shooting technique via Matlab tool bvp4c. initially, the higher-order ODE’s are converted into first order system, by implementing following procedure

Validation of results

The solution verification and validity has been checked in Table 1 with comparison the numerical with investigation of Fathizadeh et al.[32] and Fang et al.[33]. A fine accuracy of obtained results is noted with these studies.

Table 1

The comparative analysis when

M	Fathizadeh et al.[32]			Fang et al.[33]	Present results
	HPM	MHPM	Exact solution
0	1	1	1	1	1.0000
1.0	− 1.4142	− 1.4142	− 1.4142	− 1.4142	− 1.4145

Results and discussion

This section communicates the physical aspect of Sutterby nanofluid in view of flow parameters. The comparative analysis is performed for flow due to plate and cylinder Figure 2 is drawn to estimate the consequence of Hartmann number on flow velocity . The interaction of magnetic force reports a declining change in velocity due to presence of Lorentz force. Moreover, the declining change in velocity is more progressive for plate as compared to stretched cylinder. Figure 3 characterizes the impact of Sutterby fluid parameter on velocity . The velocity dwindles for increasing change in Sutterby fluid parameter . The outcomes of against Darcy resistance parameter is delineated via Fig. 4. An increasing change in velocity enhanced the for Darcy resistance parameter . Figure 5 is portrayed to understand the impact of buoyancy ratio parameter on It is analyzed that reduces for increasing variation of for both plate and stretched cylinder . Physically, the buoyancy forces play novel contribution to reduce the velocity rate effectively. Figure 6 is plotted to investigate the insight of bioconvection Rayleigh factor on . The velocity field exaggerates for larger bioconvection Rayleigh number . Figure 7 is illustrating the deviation of mixed convection parameter on for both plate and cylinder. The depicted change in show an increasing fluctuation by raise in estimation of mixed convection parameter . Physically, the mixed convection constant explores the ratio between buoyancy to viscous force. The increasing contributions of buoyancy forces results an increment in the velocity. The inspiration of thermal relaxation time parameter against temperature profile is illustrated in Fig. 8. The temperature profile reduces for larger number of thermal relaxation parameter . Figure 9 manifests the effect of Biot number on . It is seen that increased against larger values of Biot number . Physically, the Biot number present the heat transfer coefficient which enhanced the temperature profile. Moreover, the rate of heat transfer is relatively more growing for stretched cylinder as compared to plate. To consequence of sponginess parameter on , Fig. 10 is pictured. It is predicted that upsurges for larger values of sponginess parameter . Fig. 11 is inserted to envision the effect of thermophoresis parameter on thermal field of fluid. The thermophoresis phenomenon is based on the collection of nanoparticles which migrated to the cooler surface because of temperature gradient. This fluctuation in temperature due to thermophoresis phenomenon increase the temperature profile. Figure 12 reports the inspiration of Prandtl number on . The lower temperature changes have been noted due is noted. Physically, the increasing outcomes in Prandtl number declined the thermal diffusivity due to which declined. Figure 13 impacted the change in due to temperature ratio parameter . It is analyzed that shows an enlarging trend for growing values of temperature ratio parameter . The effect of solutal relaxation time parameter against concentration field is deliberated in Fig. 14. The concentration profile reduces for . The features of activation energy parameter on is scrutinized in Fig. 15. It is perceived that enhanced by raising the numbers of activation energy parameter . The activation energy determines the minimum energy supply to start the reaction phenomenon. The presence of activation energy enhanced the concentration change more effectively. The physical aspect of Lewis number on is examined via Fig. 16. The concentration reduces by augmenting the values of Lewis number . Physically, this reduction in concentration is due to low mass diffusivity of nanoparticles associated to the higher values of Lewis number. The physical outcomes of Brownian motion parameter against of is considered in Fig. 17. The diminishes change in against larger is noticed. The characteristics of thermophoresis parameter against is characterized in Fig. 18. The concentration rate is dwindles with larger thermophoresis parameter . Figure 19 is presented to estimate the variation of bioconvection Lewis number on microorganism field for plate and cylinder It is analyzed that microorganism field declines with enlarge numbers of . The physical features of Peclet number on disclosed in Fig. 20. The lower microorganism rate for higher is observed. The lower microorganism profile with higher Peclet number is owing to low motile diffusivity.

Figure 2

Profile of for .

Figure 3

Profile of for .

Figure 4

Profile of for .

Figure 5

Profile of for .

Figure 6

Profile of for .

Figure 7

Profile of for .

Figure 8

Profile of for .

Figure 9

Profile of for .

Figure 11

Profile of for .

Figure 12

Profile of for .

Figure 13

Profile of for .

Figure 14

Profile of for .

Figure 15

Profile of for .

Figure 16

Profile of for .

Figure 17

Profile of for .

Figure 18

Profile of for .

Figure 19

Profile of for .

Figure 20

Profile of for .

The numerical data is achieved for inspecting the heat transfer rate, mass transfer pattern, motile density and wall shear forces in tables (2, 3, 4, 5). Table 2 suggests that the wall shear rate enhance with Hartmann number while it reduces for mixed convection parameter. The numerical observations are relatively growing for stretched cylinder. The numerical outcomes listed in Table 3 explores the Nusselt number variation when different numerical values are assigned to parameters. The decreasing numerical data is achieved for thermophoresis constant while increasing observations are predicted against thermal relaxation constant. Table 4 notify that the Sherwood number numerical variation is larger for Lewis number and concentration relaxation number. From Table 5, the enhanced numerical data is results against the Peclet number.

Table 2

Outcomes of versus , , , , and .

Parameters						-f^″0
M	S^∗	A_1	A_2	α_1	α_2	β=0.0	β=0.3
0.1	0.2	0.1	0.1	0.1	1.0	1.0193	1.3106
0.8						1.3058	1.5981
1.6						1.5705	1.8620
0.5	0.1	0.1	0.1	0.1	1.0	1.1914	1.3598
	1.0					1.1887	1.3469
	2.0					0.1825	1.3376
0.5	0.2	0.2	0.1	0.1	1.0	1.1911	1.4835
		0.8				1.1892	1.4819
		1.6				1.1866	1.4799
0.5	0.2	0.1	0.2	0.1	1.0	1.1922	1.4846
			0.8			1.1969	1.4895
			1.6			1.2030	1.4961
0.5	0.2	0.1	0.1	0.1	1.0	1.1923	1.4787
				0.6		1.1975	1.4867
				1.2		1.2042	1.4939
0.5	0.2	0.1	0.1	0.1	0.1	0.1936	1.4809
					0.6	0.1865	1.4650
					1.2	0.1821	1.4318

Table 3

Outcomes of versus , , , , , , , , and .

Parameters										-θ^′0
M	S^∗	A_2	A_1	Nt	Nb	Pr	Rd	λ_T	Le	β=0.0	β=0.3
0.1	0.2	0.1	0.1	0.3	0.2	2	0.8	0.3	2.0	0.6059	0.5734
0.8										0.5726	0.5399
1.6										0.5423	0.5423
0.5	0.1	0.1	0.1	0.3	0.2	2	0.8	0.3	2.0	0.5724	0.5480
	1.6									0.5836	0.5317
	2.0									0.5980	0.5058
0.5	0.2	0.2	0.1	0.3	0.2	2	0.8	0.3	2.0	0.5858	0.5528
		0.8								0.5849	0.5517
		1.6								0.5836	0.5503
0.5	0.2	0.1	0.2	0.3	0.2	2	0.8	0.3	2.0	0.5859	0.5530
			0.8							0.5960	0.5529
			1.6							0.5861	0.5528
0.5	0.2	0.1	0.1	0.1	0.2	2	0.8	0.3	2.0	0.5954	0.5643
				0.6						0.5714	0.5359
				1.2						0.5411	0.5008
0.5	0.2	0.1	0.1	0.3	0.1	2	0.8	0.3	2.0	0.5952	0.5560
					0.6					0.5711	0.5423
					1.2					0.5352	0.5081
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.8	0.3	2.0	0.4603	0.4295
						2.2				0.6109	0.5782
						3.2				0.7109	0.6802
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.1	0.3	2.0	0.5767	0.5432
							0.6			0.5841	0.5510
							1.2			0.5885	0.5559
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.8	0.1	2.0	0.5811	0.5412
								0.6		0.5865	0.5532
								1.2		0.5889	0.5537
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.8	0.3	1.0	0.5874	0.5545
									2.0	0.5816	0.5487
									3.0	0.5781	0.5453

Table 4

Outcomes of versus , , , , , , , and

Parameters									-ϕ^′0
M	S^∗	A_2	A_1	Nt	Nb	Pr	λ_C	Le	β=0.0	β=0.3
0.1	0.2	0.1	0.1	0.3	0.2	2	0.3	2.0	0.9089	0.8601
0.8									0.8589	0.8098
1.6									0.8135	0.7663
0.5	0.1	0.1	0.1	0.3	0.2	2	0.3	2.0	0.8642	0.8407
	1.6								0.8763	0.8371
	2.0								0.8860	0.8185
0.5	0.2	0.2	0.1	0.3	0.2	2	0.3	2.0	0.8787	0.8293
		0.8							0.8773	0.8276
		1.6							0.8754	0.8254
0.5	0.2	0.1	0.2	0.3	0.2	2	0.3	2.0	0.8789	0.8295
			0.8						0.8790	0.8294
			1.6						0.8791	0.8293
0.5	0.2	0.1	0.1	0.1	0.2	2	0.3	2.0	0.2977	0.2821
				0.6					1.7141	1.6076
				1.2					3.2467	3.0050
0.5	0.2	0.1	0.1	0.3	0.1	2	0.3	2.0	1.8632	1.6590
					0.6				0.8512	0.2767
					1.2				0.8489	0.1383
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.3	2.0	0.6904	0.6443
						2.2			0.9163	0.8673
						3.2			1.0663	1.0204
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.1	2.0	0.8518	0.8421
							0.6		0.8650	0.8322
							1.2		0.8719	0.8272
0.5	0.2	0.1	0.1	0.3	0.2	1.2	0.3	1.0	0.8811	0.5767
								2.0	0.8924	0.5841
								3.0	0.9672	0.5885

Table 5

Outcomes of versus , , , and .

Parameters					-χ^′0
M	A_1	A_2	Pe	Lb	β=0.0	β=0.3
0.1	0.1	0.1	0.1	2	0.9749	0.9031
0.8					0.9062	0.8368
1.6					0.8459	0.7813
0.5	0.2	0.1	0.1	2	0.9334	0.8625
	0.8				0.9335	0.8624
	1.6				0.9337	0.8623
0.5	0.1	0.2	0.1	2	0.9330	0.8622
		0.8			0.9371	0.8601
		1.6			0.9286	0.8572
0.5	0.1	0.1	0.2	2	1.0098	0.9285
			1.2		1.8029	1.6069
			2.2		2.6763	2.3147
0.5	0.1	0.1	0.3	2.4	1.0393	0.9666
				2.8	1.1374	1.0636
				3.2	1.2289	1.1545

Conclusions

The bioconvective thermal determination of Sutterby nanofluid confined via stretched cylinder has been evaluated numerically. The aspect of Darcy resistance for nonlinear radiated flow is also inspected. The numerical outcomes are listed with shooting solver which are further verified to maintain the accuracy. The major results are:

The declining velocity change for Sutterby nanofluid is observed for Sutterby fluid parameter.

The increment in velocity change due to Darcy resistance factor is predicted for both stretched cylinder and plate. However, the reducing change in velocity is comparatively progressive for plate.

The temperature profile for plate and moving cylinder is lower subject to the increasing fluctuation of thermal relaxation constant.

The stronger heat transmission is observed for sponginess parameter and Biot constant.

The increasing concertation change of Sutterby nanofluid is noted for activation energy and thermophoresis factor.

The increasing numerical values of local Nusselt number are predicted for thermal relaxation time constant.