Nonlinear radiative Maxwell nanofluid flow in a Darcy-Forchheimer permeable media over a stretching cylinder with chemical reaction and bioconvection.

Studies accentuating nanomaterials suspensions and flow traits in the view of their applications are the focus of the present study. Especially, the usage of such materials in biomedical rheological models has achieved great importance. The nanofluids' role is essential in the cooling of small electronic gizmos like microchips and akin devices. Having such exciting and practical applications of nanofluids our goal is to scrutinize the Maxwell MHD nanofluid flow over an extended cylinder with nonlinear thermal radiation amalgamated with chemical reaction in a Darcy-Forchheimer spongy media. The presence of gyrotactic microorganisms is engaged to stabilize the nanoparticles in the fluid. The partial slip condition is considered at the boundary of the stretching cylinder. The Buongiorno nanofluid model is betrothed with impacts of the Brownian motion and thermophoresis. The analysis of entropy generation is also added to the problem. The highly nonlinear system is tackled numerically is addressed by the bvp4c built-in function of the MATLAB procedure. The outcomes of the prominent parameters versus embroiled profiles are portrayed and conversed deeming their physical significance. It is perceived that fluid temperature is augmented for large estimates of the radiation and Darcy parameters. Moreover, it is noticed that the magnetic and wall roughness parameters lower the fluid velocity. To corroborate the presented results, a comparison of the current study with a previously published paper is also executed. An outstanding correlation in this regard is attained.

Introduction

Nanofluid, an arising field of engineering, has caught the eye of numerous researchers who were observing the ways to improve the efficiency of cooling measures in industries. Nanofluids are used to improve rates of heat transfer in an assortment of applications including nuclear reactors, transportation industry, mechanical cooling applications, heat exchangers, micro-electromechanical systems, fiber, and granular insulation, chemical catalytic reactors, packed blood flow in the cardiovascular system engaging the Navier–Stokes equation. Advanced thermal features of the nanofluid are imperative in many fields like pharmaceutical, air-conditioning, micromanufacturing, microelectronics, power generation, thermal therapy for cancer surgery, transportation, chemical, and metallurgical engineering fields, etc. Due to the significant advancement in aerodynamics automotive, there is great importance in breaking down systems by direct heat dissipation. Many investigators have recently added some work to promote solar cells with high digestion of solar radiation. As Choi and Eastman[1] found that the incorporation of nanoparticles to the base liquids significantly enhances their thermal efficiency. The rising demand for highly efficient cooling devices encourages Koo and Kleinstreuer[2] to study the steady laminar nanofluid flow in micro heat sinks. It is noticed that very low nanoparticle concentration in nanofluids results in a higher thermal conductivity that exhibits a remarkable state of nanofluids[3,4]. Bilal et. al[5] scrutinized the numerical study of unsteady Maxwell flow of nanofluid influenced by the magnetic field, melting heat, and the Fourier and Fick laws. This investigation reveals that the liquid temperature is dropped for versus melting heat and unsteadiness parameters. The flow of 3D non-radiative Maxwell nanofluid with thermal and solutal stratification with chemical reaction is analytically studied by Tlili et. al[6]. Here, the noticeable outcome of the model is that the fluid concentration and temperature are declined for solutal and thermal stratifications respectively. Farooq et al.[7] analytically conversed the Maxwell nanofluid flow over an exponentially extended surface. Various researchers revealed the numerous aspects of the Williamson nano liquid[7-22].

The term Darcy–Forchheimer comes from the law of Darcy which interprets the liquid flow along a spongy channel. This law was originated and dependent upon the consequences of analysis on the water flow across the beds of sand. Movements in the spongy medium in which inertial effects are prominent come with the variations of Reynolds numbers. Therefore, this introductory term is adding up to the Darcy equation and is referred to as the Darcy–Forchheimer term. This term represents the non-linear behavior of the flow data versus pressure difference. With wide utilization of grain stockpiling, petroleum technology, frameworks of groundwater and oil assets, this Darcy law is of immense importance in the field of Fluid Mechanics. In places where the porous medium has larger flow rates due to non-uniformity, such as near the wall, Darcy's law is not applicable. The substance with stomata is named as a permeable medium. It includes an application of large numbers such that oil manufacturing, liquid flow in catalytic vessels, and reservoirs, etc. The suggestion of the fluid flow passes a porous surface was first given by Darcy[23] in 1856. However, this idea couldn't be so famous inferable from its restrictions of lower porosity and smaller speed. Afterward, Forchheimer[24] amended the equation of momentum by adding the square velocity condition into the Darcian velocity to convey the undeniable lack. Muskat[25] later call it the "Forchheimer term" which is true of the high Reynolds number. Pal and Mondal[26] addressed the Darcy–Forchheimer model over permeable media past the linearly expanded region and assume that concentration distribution is diminishing function of the electric field parameter. The movement of the hydromagnetic nano liquid past the Darcy–Forchheimer media forum effect on the boundary condition of second order is mathematically evaluated by Ganesh et al.[27]. Alshomrani et al.[28] explained the 3D Darcy–Forchheimer law with carbon nanotubes and homogeneous heterogeneous reactions. The viscous nanofluid with Darcy–Forchheimer effect over a curved area is analyzed by Saif et al.[29]. Seth et al.[30] examined mathematically the movement of carbon nanotubes over a porous Darcy–Forchheimer media in a rotating frame and many therein[31-41].

In numerous processes including dispersion of nutrients in nerves, condensation in mixtures, and thermal insulation, mass transfer plays a vital role. One live example of transfer of mass may be seen in the living matter processes like respiration and sweating. A good number of studies may be quoted where chemical reactions play a vital role in mass transfer procedures.

Recently, Mahmood[42] explored the nanofluid flow with an amalgamation of the CNT’s of both types and the engine oil over a stretched surface with the impact of the activation energy merged with the chemical reaction. The problem is solved numerically and with a surface response statistical technique. It is inferred from this model that the surface drag coefficient is negatively sensitive concerning the magnetic parameter. The numerical solution of the nanofluid flow involving the CNT’s and water over an extended/contracting sheet with quartic autocatalysis chemical reaction and Thompson and Troian slip boundary conditions is discussed by Ramzan et al.[43]. The results exposed that the fluid concentration is enhanced for quartic autocatalysis chemical reaction. Khan et al.[44] numerically tackled the Carreau nanofluid flow over an extended surface in a Homann stagnation region with chemical reaction and modified Fourier law using shooting scheme. The study divulges that the fluid velocity hinders owing to the Hartmann number and the porosity parameter. In the perspective of its clarity, the remarkable work of current researchers, see few studies[45-56].

In thermodynamics, entropy is an essential concept. One of the most powerful methods for investigating the efficiency of thermal systems is entropy generation analysis. The idea of irreversibility is inextricably related to the concept of entropy. All have an instinctive understanding of irreversibility. For example, by playing a video game in both forward and reverse, we can merely explain the irreversibility phenomenon by using backward order. There are numerous forward procedures of daily life that cannot be undone, such as pouring water into a bottle, egg unscrambling, unconstrained expansion of fluids, plastic deformation, gas uprising from the chimney, etc. In this perspective, Yusuf et al.[57] explored the entropy generation in a Maxwell fluid flow over an inclined extended surface in a non-Darcian spongy media with thermal radiation. An analytical solution of the erected mathematical model is attained. The major result of the presented model is that the rate of the entropy generation is boosted for the local inertial coefficient parameter. In a recent study, Adesanya et al.[58] performed the entropy generation appraisal for a couple stress fluid film flow on an inclined heated surface with viscous dissipation impacts. In this analysis, it is comprehended that the fluid temperature and velocity show opposing tendency versus the couple stress parameter. Furthermore, many investigators have worked on entropy generation analysis on a wide number of geometries that may found in[58-65].

The studies deliberated above reveal that abundant literature discussing the nanofluid flow of an extended cylinder is available under the influence of varied impacts. Nevertheless, comparatively less literature can be witnessed that discusses the flow of Maxwell nanofluid over an extended cylinder. But so far one has discussed the Maxwell nanofluid flow past an extended cylinder with thermal radiation, chemical reaction, and gyrotactic microorganisms with partial slip in a Darcy–Forchheimer spongy medium. The problem is solved numerically, and pertinent graphs are plotted versus the involved profiles with logical descriptions. Table 1 illustrates the originality/uniqueness of the stated fluid mathematical model by assessing it with the available researches.

Table 1

A literature analysis for the individuality of the stated model.

Authors	Buongiorno model	Maxwell nanofluid flow over a cylinder	Darcy–Forchheimer impact	Nonlinear thermal radiation	Bioconve-ction	Chemical reaction
Islam et al.[15]	√	√	×	×	×	×
Ahmed et al.[16]	√	√	×	√	×	×
Hayat et al.[66]	√	√	×	×	×	√
Raju et al.[67]	√	√	×	×	×	×
Present	√	√	√	√	√	√

(√) means said effect present, and ( ×) signifies the impact is absent.

Mathematical modeling

We examine an incompressible flow outside a cylinder having radius R and the constant temperature . As the axial direction of a cylinder along the x-axis while radial direction along r-axis. A stretching surface of the cylinder has velocity , where is the characteristic length and shows the reference velocity. The flow situation induced by a magnetic field of intensity is displayed in Fig. 1.

Figure 1

The geometry of the flow.

The resulting boundary layer equations defining the depicted scenario are given as[66]:

With boundary conditions[66]:

The radiative heat flux is given by:

The following transformation is used to obtain the non-dimensional structure of the above-mentioned flow model[66]:

Equation (1) trivially fulfilled and Eqs. (2–6) are as follows:

Particular dimensionless parameters emerging in the above equations are portrayed as:

The dimensional form of drag force coefficients, rate of mass transfer, rate of heat transfer, and Motile microorganisms are appended as below:

With

The Drag force coefficient, mass transfer rate, Motile microorganism, and rate of heat in the dimensionless form are supplemented below:

Rate of entropy generation (EG)

The volumetric equation is represented as:

The characteristics EG is framed as:

The entropy generation is given as the quotient of the and , i.e.,

In dimensionless form:

The parameters used in Eq. (22) are defined as:

Numerical procedure

For the nonlinear arrangement of equations and boundary conditions (9)–(13) the finite difference MATLAB bvp4c procedure is applied which is solid at 4th order and the grid size of 0.01 is viewed as acknowledged 10–6. The numerical plan requires the change of higher-order differential equations into one-order differential equations.

Results with discussion

In this segment, we will examine the effect of distinct parameters on velocity, concentration, temperature, and gyrotactic microorganism fields. The numerous parameters like the magnetic interaction parameter Darcy parameter radiation parameter Schmidt number temperature ratio parameter porosity parameter Deborah number thermophoresis parameter curvature parameter bioconvection Lewis number Prandtl number Brownian motion parameter wall roughness parameter chemical reaction parameter concentration slip parameter thermal slip parameter Peclet number Bioconvection parameter and reaction order are discussed on temperature, velocity and nanoparticles concentration, and gyrotactic microorganism fields. Figure 2 demonstrated the behavior of on the . The strength of the Lorentz force is measured by The increase in enhances the Lorentz force strength and due to the rise in the velocity in axial direction decreases. As a result, the gradient of velocity at the surface is decreased. The wall effect on the velocity profile is described in Fig. 3. The slip decreases the speed close to the disk and this condition enhances by increasing in . Practically, the stretched impact of a cylinder is moderately shifted to the liquid layers which result in a decrease in . The influence of on the velocity field is investigated in Fig. 4. It is examined that by escalating the variations of , the decreasing trend of the velocity field is seen. This is because the higher values of produce resistance in a liquid flow and hence velocity decreases. Figure 5 demonstrated the effect of the on the velocity distribution of . The liquid's velocity diminishes on greater estimations of the . Actually, the movement of the liquid is stalled because of the presence of permeable media, and this results in the falloff of the liquid velocity. The effect of on the velocity field is investigated in Fig. 6. It is examined that by escalating the variations of , the diminishing behavior of the velocity field is seen. The impact of the upon is explained by Fig. 7. A significant increase in is observed. Enhancing signifies the temperature of the wall that causes thicker penetration depth for temperature profile. Likewise, the thermal diffusivity lies in the boundary layer with the relating exchange of heat. The thermal boundary layer corresponds to be larger nearby the region where the hotness is larger while it is lower a long way from cylinder because here temperature is low when compared to others. Subsequently, an intonation point emerges on the region when greater is considered. The temperature field for numerous is shown in Fig. 8. A generous upgrade in the temperature of the liquid is seen when the radius of the cylinder is reduced. The effect of on the thermal profile is described in Fig. 9. It is observed that the existence of melting phenomenon of the liquid temperature increases with rising variations of . Therefore, we can judge that greater variation of enhances the temperature field. Figure 10 defines that the large estimation of the descends the dimensionless liquid's temperature. From the figure, it is noticed that the thermal boundary layer becomes thicker on enhancing the values of the curvature parameter. Figure 11 is outlined to show the plots of for numerous terms of when other variables are fixed. It can be judged that growing values of increase the temperature and its parallel thickness of layer become thicker. Figure 12 elucidates that an increment in decays the nanoparticle concentration distribution . There is an opposite relationship between the and the Brownian diffusion coefficient. Greater the values of Schmidt number lower will be the Brownian diffusion coefficient, which tends to decrease the . Figure 13 portrays the concentration field for different estimations of . Large variations of the tend to smaller the nanoparticle concentration field. The descending behavior of a concentration profile on a is drawn in Fig. 14. Figure 15 demonstrated that for greater values of reaction order the concentration profile becomes higher. Figure 16 depicts the influence of on . Both the concentration and thermal layer thickness are increased by accumulating the variations of the . Greater estimations of the give rise to thermophoresis force which increases the movement of nanoparticles from cold to hot surfaces and also increases in the thermal layer thickness. The descending behavior in concentration distribution against is shown in Fig. 17. An enhancement in increases the Brownian motion due to which there is an escalation in the movement of nanoparticles and hence boundary layer thickness reduces. Figure 18 plotted to draw the curves of for different terms of while other variables are fixed. It is observed that depicts the decreasing behavior for large values of . Figure 19 shows the behavior of on gyrotactic microorganisms’ profile . Here, is an increasing function of that effects to a decrease the diffusivity of microorganisms. Figure 20 indicates the variations in gyrotactic microorganism profile for distinct estimations of the. Large variations of decrease the gyrotactic microorganism field. Figure 21 is drawn to illustrates the impact of the Brinkman numbers on the Entropy generation number. It is found that Entropy escalates for the Brinkman number. In Fig. 22, with rising estimations of , an increase is seen in the entropy generation number. A higher causes more aggravation in the field and expand fluid friction and heat transfer, which eventually increases the rate of entropy in the boundary layer region.

Figure 2

Plot of for .

Figure 3

Plot of for .

Figure 4

Plot of for .

Figure 5

Plot of for .

Figure 6

Plot of for .

Figure 7

Plot of for .

Figure 8

Plot of for .

Figure 9

Plot of for .

Figure 10

Plot of for .

Figure 11

Plot of for .

Figure 12

Plot of for .

Figure 13

Plot of for .

Figure 14

Behaviour of for .

Figure 15

Plot of for .

Figure 16

Plot of for .

Figure 17

Plot of for .

Figure 18

Plot of for .

Figure 19

Plot of for .

Figure 20

Plot of for .

Figure 21

Plot of for .

Figure 22

Plot of for .

Table 2 is generated to substantiate the outlined results in this study by comparing it with Khan and Mustafa[68] and Tamoor et al.[56] in limiting case. A good agreement between the two outcomes is seen. Besides, Tables 3, 4, 5 express the numerical variations of the local Sherwood number , local Nusselt number and density amount of motile microorganism for distinct estimations of , , , , , , , , , , , , and . It is witnessed in Table 3 here that the heat flux rate is escalated for the growing estimates of , but the opposite trend is perceived for the estimations of , , In Table 4, it is witnessed that the mass flux rate is declined for the values of the , and , however, it is enhanced for the increasing estimates of the and The behavior of the varied parameters versus the density amount of motile microorganism is portrayed in Table 5. It is renowned that density amount of motile microorganism is improved for the estimations of , , , , and .

Table 2

Validation of numerical outcomes for with Khan and Mustafa[68] and Tamoor et al.[56] and when  =  = 0.

K	-f^″0	-f^″0	-f^″0
	[68]	[56]	Present
0	1	1	1
0.2	1.0198039	1.01980	1.01981
0.5	1.1180340	1.11803	1.11803
0.8	1.2806248	1.28063	1.28062
1	1.4142136	1.41421	1.41421

Table 3

Computations of for various variations of , , , when  = 7 and  =  = 0.5.

K	M	θ_w	Rd	Re_x^-12Nu_x
0.5	0.2	1.5		0.7654
1				0.6803
1.5				0.5851
2				0.5003
	0.5			0.7745
	0.7			0.7801
	1			0.7883
		2		1.1440
		2.5		1.4871
		3		1.7770
			0.1	1.3736
			0.3	2.2339
			0.5	3.3713

Table 4

Computations of for numerous variations of ,, and when  = 7,  = 0.5 and  =  =  = 0.5.

M	Sc	γ	n	Re_x^-12Sh_x
0.2	5	1	1	1.15817
0.5				1.17480
0.7				1.18683
1				1.20568
	2			0.97565
	3			1.05478
	7			1.22726
		2		1.27036
		3		1.34142
		4		1.39264
			2	0.99641
			3	0.92496
			5	0.87854

Table 5

Computations of for various variations of , , , and .

M	Pr	σ	Pe	Lb	Re_x^-12Nn_x
0.3	1.5	0.2	0.3	0.2	0.935596
0.6					1.059580
0.9					1.193600
	1.0				0.915801
	1.5				0.935596
	2.0				0.940771
		0.2			0.935596
		0.4			1.012580
		0.6			1.089560
			0.4		1.065290
			0.5		1.189010
			0.6		1.307130
				0.3	0.963155
				0.6	1.038210
				0.9	1.103780

Concluding remarks

In the current investigation, we have discussed nonlinear radiative MHD Williamson nano liquid flow via a stretched cylinder in a Darcy–Forchheimer porous media. The flow is assisted by the impacts of the chemical reaction, and gyrotactic microorganisms with partial slip condition at the boundary. The solution to the problem is addressed by the MATLAB scheme of the bvp4c built-in function. The main results of the present investigation are appended below:

and show the opposing nature against the concentration field.

The velocity distribution is lowered for large variations of , and

decreases the gyrotactic microorganism profile.

An increment in , and leads to a lowering concentration profile.

Enhanced variations of and display the decreasing behavior on the temperature profile.

Large values of , and causes an increment in temperature distribution.

Gyrotactic microorganism profile reduces for greater variations of and .

Entropy intensifies for ; however, an opposite tendency is observed for .