Exploring the magnetohydrodynamic stretched flow of Williamson Maxwell nanofluid through porous matrix over a permeated sheet with bioconvection and activation energy.

The evolution of compact density heat gadgets demands effective thermal transportation. The notion of nanofluid plays active role for this requirements. A comparative account for Maxwell nanofluids and Williamson nanofluid is analyzed. The bioconvection of self motive microorganisms, non Fourier heat flux and activation energy are new aspects of this study. This article elaborates the effects of viscous dissipation, Cattaneo-Christov diffusion for Maxwell and Williamson nanofluid transportation that occurs due to porous stretching sheet. The higher order non-linear partial differential equations are solved by using similarity transformations and a new set of ordinary differential equations is formed. For numerical purpose, Runge-Kutta method with shooting technique is applied. Matlab plateform is used for computational procedure. The graphs for various profiles .i.e. velocity, temperature, concentration and concentration of motile micro-organisms are revealed for specific non-dimensional parameters. It is observed that enhancing the magnetic parameter M, the velocity of fluid decreases but opposite behavior happens for temperature, concentration and motile density profile. Also the motile density profile decrease down for Pe and Lb. The skin friction coefficient is enhanced for both the Williamson and Maxwell fluid.

Introduction

In modelling of fluids flows with shear-dependent viscosity, the power-law model is commonly used. But the consequences of elasticity are impossible to be anticipated. The results of elasticity can be achieved by second grade or third grade fluids. But the viscosity is not shear-dependent in these models. In addition, they are unable to determine the impacts of calming stress. A subdivision of fluids of the rate form, viz., Maxwell model can predict stress relaxation and has thus become more common. As in the Maxwell model, a strictly viscous damper and a purely elastic spring can be described. Many researchers have pondered Maxwell nanofluid flow simulations. Maxwell nanofluids with normal kernel through fractional derivatives studded by Abro et al.[1]. Numerical analysis of Maxwell nanofluid flows past a linearly stretched layer was reviewed by Sharama et al.[2]. Assessment of bioconvection in Maxwell nanofluid configured with nonlinear thermal radiation and activation energy using a Riga surface was conducted by Ramesh et al.[3]. Ahmed et al.[4] reviewed the impact of radiative heat flux in Maxwell nanofluid flow over a chemically reacted spiralling disc. Mathematical analysis of Maxwell nanofluid with hydromagnetic dissipative and radiation was considered by Hussain et al.[5]. Jawad et al.[6] calculated the entropy generation for magnetohydrodynamic (MHD) mixed convection and Maxwell nano-fluid flow over an elongating and penetrable surface in the presence of heat conductivity, velocity slip boundary condition and thermal radiation.

English scientist Williamson presented the Williamson model[7] in 1929 and several other researchers studied it. Williamson fluid has a shear thinning property typical of a non-Newtonian fluid model. Later, perturbation solution of incompressible flows in a rock fracturing with a non-Newtonian Williamson fluid was studied by Dapra et al.[8]. Nadeem et al.[9] studded Wiliamson nanofluid stream on an extended sheet in twenty first century. Even now in 2020, several scholar are rivaling the properties of Williamson nano-fluid such as, analysis of gyrotactic microorganisms with activation energy in the radiated liquid was provoked Haq et al.[10]. The impact of the magnetic field in a curved channel on the peristaltic flow of Williamson fluid was explored by Rashid et al.[11]. Habib et al.[12] carried out a comparison research on micropolar, Williamson, and Maxwell nanofluids flow owing to a stretched surface in the presence of bioconvection and double diffusion along the activation energy. Waqas et al.[13] conducted research for the generalized principle of Fourier and Fick to effect Williamson fluid flow. Analysis of heat transfer for electro-osmotic stream across of the micro-channel for Williamson fluid motion was discussed by Noreen et al.[14].

Hydrodynamic volatility and patterns in suspensions of biased swimming microorganisms are characterized by the term bio-convection. In the early 21st century, biotechnology has developed to incorporate modern and complex disciplines. Bioconvection has various uses in natural systems and biotechnology. Various researchers commit their resources to reveal the characteristics of bio-convection. Bioconvection has various uses in natural systems and biotechnology. Numerical study of bio-convective nanofluid on a magnetohydrodynamic slip flow with stephan blowing was considered by Tuz et al.[15]. The effects of energy on Eyring nanofluid flow containing motile microorganism was studied by Sharief et al.[16]. Sharma et al.[17] adopted non-dimensional measurements to examine the effect of heat and mass flux on the natural convective laminar movement of a viscoelastic immiscible fluid. Sharma et al.[18] studied the impacts of a chemical change and a heat source on magneto-hydrodynamic assorted convective mass and heat transfer flow over a vertical plate. For a Casson nanofluid stream over a stretching/shrinking sheet, magagula et al.[19] tested the influence of an applicable magnetic field, nonlinear thermal radiation and first-order chemical reaction double dispensed bioconvection. Shaw et al.[20] presented a homogeneous model for oxytactic microbial bioconvection in a non-Darcy permeable material. Khan et al.[21] studied the impacts of Arrhenius activation energy, binary chemical reaction and viscous dissipation on the bioconvective micropolar flow of nanofluid over a slim moving needle containing gyrotactic micro-organisms.

A highly promising nanostructure analysis was primarily conducted by Choi et al.[22] to show that the presence of nanoparticles in base fluids can be beneficial in improving the rheological effects of typical base fluids. “Nanofluid is a suspension comprising nanoparticles in a base liquid (water, base fluid mixture, kerosene, bio-fluids and organic liquids) that alters the viscosity, thermal conductivity, density, and mass diffusivity of the base fluid”.In the last few decades, companies have evolved to research fluid mechanics on the nano and micro scales. For improved oil recovery in a heterogeneous two-dimensional anticline geometry, nanofluid flooding investigated by Esfe et al.[23]. Khader et al.[24]utilized the fourth order predictor-corrector finite difference technique to analyse the influence of thermal radiation and a non-uniform heat source on unsteady MHD micropolar fluid flow past a stretching sheet. Verma et al.[25] examined the steady two-dimensional, laminar, viscous and incompressible boundary layer movement of nanoparticles. kumaran et al.[26] employed the Keller box technique to examine the effect of thermal conductivity variation and thermal radiation on chemically reacting and free convective Powell-Eyring nanofluid flow across a cylinder. Sharma et al.[27] explored an irregular magneto-hydrodynamic (MHD) natural convection transmission of mass and thermal over a porous medium sheet under the impact of thermal radiation and thermo-diffusion consequences using the Laplace transform tactic. Jena et al.[28] numerically investigated nanofluid stream transmission through permeable media with a heat source/sink and a chemical process. Nayak et al.[29] discussed slip consequences of the flow of chemically reactive Casson nanofluid flowing over an exponentially stretched electromagnetic sheet. Saranya et al.[30] researched the two-dimensional consistent convective boundary layer fluid motion and heat transmit of Newtonian/non-Newtonian base liquids with magnetic nanoparticles over a flat plate with non-linear thermal radiation and slip implications. Kairi et al.[31] assumed that the bioconvection of Casson nanoparticles over an inclined elongating sheet. Sharma et al.[32] exploited the perturbation approach to analyse the electrically conducting fluid motion over a vertical plane sheet. The diverse applications of nanofluids in various aspects are examined by some other scholars[33-36].

In1942, Hannes Alfvén made the first recorded use of the term “Magnetohydrodynamics”.It originated from the magneto-meaning magnetic field, hydro-meaning water, and dynamics-meaning motion. MHD is utilized in divers fields of engineering, industrial and mostly in the nanotech field. All around the globe, scholars are devoting their time and efforts to reveal the characteristics. In recent time, many scholars worked on the properties of MHD as discussed Rasool et al.[37]. Ahmad et al.[38] applied the homotopy analysis method (HAM) to examine the fluid film flow of an Oldroyd-B fluid past a spinning disc in three-dimensional space under the implications of heat absorption/omission and radiation terminologies. Gul et al.[39] employed the RK-4 method to investigate the MHD (magneto-hydrodynamics) unsteady and immiscible flow of nanofluid flow caused by a stretching rotating disc with the effect of Joule heating and dispersion. Gul et al.[40] investigated the 3-D Darcy-Forchheimer MHD Casson fluid and steady flow between the difference between a disc and a cone in a rotating strategy in the presence of thermal diffusion effect and Brownian movement[41,42] scrutinized the MHD flow attributes in distinct aspects.

The impact of bioconvection of micro-organism for two different flows of Williamson nanofluid and Maxwell nanofluid is rarely studied in the existing literature. It helps to cope with the possible settling of nano-entities and thus improvement of thermal conductivity of the common fluids can be established. The novelty of this work pertains to Cattaneo–Christov heat flux, bioconvection and activation energy. The buoyancy effects, magnetic field strength imparts differentiated influence on flow temperature and concentration characteristics. The role of these useful and practicable physical aspects for two significant non-Newtonian fluids is enumerated and their comparative outcomes help to broaden our understanding and utilization of these flows. By utilization of similarity transforms enable to yield numerical solution of non-linear coupled system with the implementation of Runge–Kutta method and shooting technique. By utilizing the parametric evolution of bioconvection, nanofluid slips and Cattaneo–Christov factor, the findings are evaluated for Maxwell fluid and Williamson fluid. These results are beneficial for thermal management of heat exchangers of the emerging technologies.

Physical model and mathematical formulation

A two-dimensional, steady flow of Maxwell and Williamson nanofluid with activation energy and Cattaneo–Christov diffusion embedded in a porous medium is discussed in the existence of bioconvection of micro-organisms. The fluids flow owing to a stretching sheet, The wall velocity is , D is stretching constant. is the stretching velocity of sheet along x-axis. It is assumed that behavior of the flow is linear due to the impermeable stretching sheet. Here and are the velocities of fluids flow along x-axis and y-axis. In Fig. 1, the temperature of fluids, concentration and density of microorganisms near the sheet are taken as T, C, and N. Now our problem represents the basic equations of continuity, equation of momentum, energy equation, concentration equation and gyrotactic microorganism’s concentration equation in Cartesian coordinates.

Figure 1

Problem description.

The following equation are established in the aspects of flow of fluids. The governing equations for the problems are given below[5,8,16]:

Using the similarity variables:where we have is fluid relaxation time, is Williamson fluid parameter. It is to mentioned that in Eq. (2) when and , it is Williamson fluid flow. Also for and , than it is Maxwell fluid flow. Equation (1) is balanced identically. Equations (2) to (5) are rescaled as below:

The non-dimensional parameters in their respective order:

, M, Kp, , Nr, Rb, Rd, Pr, Le, Nbt, b, Q, Sc, A, , E, Lb, Pe, , and S are Deborah number, non-newtonian Williamson parameter, magnetic parameter, porosity parameter, mixed convection, buoyancy ratio, Rayleigh number, radiation, Prandtl number, Lewis number, diffusivity ratio, thermal relaxation constant, heat source, Schimdt number, dimensionless reaction rate, temperature difference, activation energy, bio-convection Lewis number, bio-convection Peclet number, microorganism concentration difference, suction/injection parameter.

Where: . , , , , , , , , , , , , , , , , , .

The physical quantities are stated as[43]:

(skin friction coefficient), (local Nusselt number), (local Sherwood number) and (local density of microorganism) are given below:where , , and denotes shear stress, surface heat flux, surface mass flux and motile microorganism flux are given by (at ),

On solving these quantities with the help of given similarity transformation, we obtain:where, is the local Reynolds number.

Solution procedure

For numerical results, nonlinear ordinary differential Equations (8) to (11) with the given boundary conditions Eq. (12) are solved by Runge–Kutta method of order four along with shooting technique. To execute this numerical method, the new variable expressed as below[44-46]:along with the boundary conditions[47]:

Results and discussions

In this part, we presented and explained the outcomes as computed from the above mentioned procedure. The validation of current results is established when these are compared with the existing findings of[48-50] in limiting cases.The four sets of results are contained in Table 1 and a close agreement is seen among them. Table 2 indicates that the drag force is supplemented with incremental variation of M, . Nr, Rb in case of Maxwell fluid as well as Williamson fluid and hence the magnitude of (skin friction factor) is enhanced because of the additional resistance to the flow that comes into play with enhancement of magnetic parameters. However, the skin friction is reciprocated against mixed convection parameter , because the rise in mixed convection helps to accelerate the flow. Further, it is seen that the absolute value of is larger for Maxwell nanofluid than that of Williamson nanofluid. Table 3 shows that Nusselt number improves in direct relation to Pr, Nbt and b, but it diminishes against Nc. Table 4 convinced that recedes only when E rises, while is upsurged directly with Sc, A, and n. From Table 5, it is noticed that magnitude of motile density is increased when Lb, Pe and are uplifted. Also, it is observed that Maxwell nanofluid takes larger values than that of Williamson nanofluid for Nusselt number , and for motile density .

Table 1

The comparative outputs for .

Pr	Nadeem et al.[48]	Khan and Pop[49]	Golra and Sidawi[50]	Present results
0.07	0.066	0.066	0.066	0.065
0.2	0.169	0.169	0.169	0.167
0.7	0.454	0.454	0.454	0.435
2.0	0.911	0.911	0.911	0.910

Table 2

Results for .

M	Kp	ω	Nr	Rb	Maxwell fluid β=0.5	Williamson fluid λ=0.1
0.1	0.1	0.1	0.1	0.1	2.5113	1.2342
0.3					2.6257	1.3090
0.5					2.7341	1.3796
0.5	0.2				2.5693	1.2722
	0.4				2.6806	1.3448
	0.5				2.7341	1.3796
	0.5	0.1			2.7341	1.3796
		0.2			2.6512	1.3339
		0.3			2.5712	1.2897
		0.1	0.1		2.7341	1.3796
			0.2		2.7407	1.3834
			0.3		2.7474	1.3873
			0.1	0.1	2.7341	1.3796
				0.2	2.7389	1.3824
				0.3	2.7438	1.3852

Table 3

Results for .

Pr	Nc	Nbt	b	Maxwell fluid β=0.5	Williamson fluid λ=0.1
0.71	0.5	2.0	0.5	0.4045	0.4418
1.0				0.5237	0.5680
1.5				0.7129	0.7639
0.71	0.1			0.4439	0.4843
	0.5			0.4045	0.4418
	1.0			0.3598	0.3934
	0.5	1.0		0.3942	0.4306
		2.0		0.4045	0.4418
		3.0		0.4081	0.4456
		2.0	0.1	0.3769	0.4063
			0.3	0.3905	0.4238
			0.5	0.4045	0.4418

Table 4

Results for .

Sc	A	δ	E	n	Maxwell fluid β=0.5	Williamson fluid λ=0.1
1.0	0.1	0.1	0.3	0.5	0.7333	0.7408
2.0					1.0157	1.0395
3.0					1.2887	1.3218
2.0	0.1				1.0157	1.0395
	0.2				1.1912	1.2084
	0.3				1.3389	1.3526
	0.1	0.1			1.0157	1.0395
		0.2			1.0261	1.0493
		0.3			1.0359	1.0586
		0.1	0.1		1.0546	1.0767
			0.3		1.0157	1.0395
			0.5		0.9818	1.0073
			0.3	0.1	1.0101	1.0342
				0.3	1.0129	1.0369
				0.5	1.0157	1.0395

Table 5

Results for .

Lb	Pe	Ω	Maxwell fluid β=0.5	Williamson fluid λ=0.1
0.4	1.2	0.2	1.6319	1.6748
0.8			1.7782	1.8269
1.2			1.9140	1.9662
1.2	0.4		1.0467	1.0892
	0.8		1.4761	1.5228
	1.2		1.9140	1.9662
	1.2	0.1	1.8102	1.8621
		0.2	1.9140	1.9662
		0.3	2.0178	2.0704

Figure 2a explains the effect of magnetic field M on velocity profile for Maxwell and Williamson fluids. A reduction is produced in the boundary layer thickness by increasing the value of M. It is observed that Williamson fluid is highly affected as compared to the Maxwell fluids. Physically, the magnetic field has ability to enhance the drag force which is said as Lorentz force. This force prevents the fluid from flowing as well as causes an increment in the thickness of boundary layer. From Fig. 2b, it is seen that the enhanced porosity parameter retards the flow velocity. Because has ability to produce resistance for fluid flow. Mixed convection parameter effect on for Williamson and Maxwell nanofluids is discussed in Fig. 3a. The greater values of mixed convection parameter causes an increment in due to the larger buoyancy force. Figure 3b and c both explain the aspects of (Rb) bioconvection Rayleigh number as well as (Nr) buoyancy ratio parameter on . The changing values of both the parameters cause decline in velocity distribution significantly. Physically, both parameters show their relation with buoyancy ratio forces which causes resistance in the motion of the two fluids. Moreover, the velocity for Williamson nanofluids is faster than that of Maxwell nanofluids. It means Maxwell fluid imparts more viscus effects on the flow than the Williamson fluid. The impacts of magnetic parameter M and Cattaneo-Christov parameter b on temperature profile are illustrated in Fig. 4. Physically, M has the ability to intensify the temperature distribution profile due to the involvement of Lorentz force. Moreover by increasing the Cattaneo-Christov parameter b causes a reduction in . Thus Maxwell fluid diffuses more heat as compared to Williamson fluid and hence the temperature curve for Maxwell fluid is elevated than that of Williamson fluid. As it is clear that is a ratio between buoyancy forces to viscous force. Due to which a little change in causes a decrement in the . The impact of Lewis number Le on is depicted in Fig. 5a. Depreciation in temperature profile is due to its direct relation with diffusivity of mass which causes reduction in the temperature distribution. The effect of Pr can be visualized on temperature profile in Fig. 5b. It is clear that by increasing the values of Pr causes a reduction in . Physically, it is due to the reverse relation of Pr with the thermal diffusivity. Fig. 5c explains the effect of Nc on temperature distribution, an increment in Nc causes an accession in temperature distribution. From Fig. 6, it can be visualized that increment in the magnetic number M as well as activation energy E causes an accession in the concentration profile . Increment in M causes an increase in Lorentz force which causes more resistance to the flow of Williamson fluid as compared to Maxwell fluid. A decrease in concentration profile can be observed in Fig. 7 with the exceeding of chemical reaction rate parameter A as well as Lewis number Le. The larger inputs of A means faster chemical reaction to decline the concentration function . Similarly, higher Lewis number Le causes depreciation in because it is related reciprocally to mass diffusivity. While the increment in A generates more resistance in Williamson fluid flow as compared to Maxwell fluids. Moreover decrement in mass diffusivity becomes Williamson fluid less resistive as compared to Maxwell fluid for flow. The motile density profile is increased with the increment in M. It is due to slowing of flow. The curves generated by M results an accession in the motile density profile . It also shows an accession in it by increasing the value of microorganism difference parameter which can be observed in Fig. 8. Figure 9 is drawn to discuss the results of Lewis number Lb as well as Peclet number Pe. The diffusivity of microorganisms is related with Peclet number Pe, which means greater value of Peclet number causes depreciation in diffusivity of microorganisms. The increment in the Peclet number Pe starts a reduction in the motile organisms boundary layer thickness, as well as a decrement in motile density profile . Moreover, it is clear that larger value of Lb causes a lesser microorganism distribution .

Figure 2

Velocity fluctuation with (a) M and (b) .

Figure 3

Velocity fluctuation with (a) , (b) Nr and (c) Rb.

Figure 4

Temperature fluctuation with of (a) M and (b) b.

Figure 5

Temperature fluctuation with (a) Le, (b) Pr and (c) Pr.

Figure 6

Concentration fluctuation with (a) M and (b) E.

Figure 7

Concentration fluctuation with (a) A and (b) Le.

Figure 8

Motile density fluctuation with (a) M and (b) .

Figure 9

Motile density fluctuation with (a) Lb and (b) Pe.

Conclusion

This article explores the effect of magnetic field, bio-convection and activation energy for Williamson and Maxwell nanofluid on a porous stretched sheet with Cattaneo–Christov diffusion. The main findings of this study are given below:

Enhancement in magnetic parameter M retards the flow of Williamson fluid and Maxwell fluids.

Rayleigh number Rb, buoyancy ratio parameter Nr, mixed convection parameter , and Kp parameters highly influence on Williamson fluid than Maxwell fluid.

For temperature profile Cattaneo Christov parameter b, Lewis number Le, and Prandtl number Pr have similar effects. Maxwell fluid is more affected by these parameters than Williamson fluid. On contrary, M and Nc shows opposite behavior for it.

For concentration profile Williamson fluid is more influenced than Maxwell fluid by M as compared to the other various parameters i.e. Activation energy E, A and Le.

For motile density profile various parameters i.e. M, microorganisms difference parameter , Lb and Peclet number Pe express more effect for Maxwell fluids as compared to Williamson fluids.

Skin friction coefficient enhanced for , M, , Nr and Rb while decrease down for and .

Nusselt number reduces when uplifting the parameter Nc, while it increases for Pr, Nbt, and b.

Sherwood number increases for Sc, A, and n while it reduces when uplifting the parameter E.

Motile density number increases for Lb, and Pe.

Future direction

These results can be further modified and improved for flow of hybrid nanofluids with the implementation of finite element or finite difference discretization.