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.


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  www.nature.com/scientificreports/ 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, biofluids 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 Cu/Ag − H 2 O 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 magnetohydrodynamic (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 twodimensional 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][34][35][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.   Where The physical quantities are stated as 43 : Cf x (skin friction coefficient), Nu x (local Nusselt number), Sh x (local Sherwood number) and Nn x (local density of microorganism) are given below: where τ w , q w , q m and q n denotes shear stress, surface heat flux, surface mass flux and motile microorganism flux are given by (at y = 0), On solving these quantities with the help of given similarity transformation, we obtain: where, (Re x ) = xU w ν 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][45][46] : , ,

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][49][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, K p . Nr, Rb in case of Maxwell fluid (β = 0.5) as well as Williamson fluid ( = 0.1) and hence the magnitude of −f ′′ (0) (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 (17) s 1 = S, s 2 = 1, s 4 = 1, s 6 = 1, s 8 = 1, atη = 0, s 2 → 0, s 4 → 0, s 6 → 0, s 8 → 0 as η → ∞.  Table 3. Results for −θ ′ (0). www.nature.com/scientificreports/ in mixed convection helps to accelerate the flow. Further, it is seen that the absolute value of −f ′′ (0) is larger for Maxwell nanofluid than that of Williamson nanofluid. Table 3 shows that Nusselt number −θ ′ (0) improves in direct relation to Pr, Nbt and b, but it diminishes against Nc. Table 4 convinced that −φ ′ (0) recedes only when E rises, while −φ ′ (0) is upsurged directly with Sc, A, δ and n. From Table 5, it is noticed that magnitude of motile density −χ ′ (0) 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 −θ ′ (0) , −φ ′ (0) and for motile density −χ ′ (0). Figure 2a explains the effect of magnetic field M on velocity profile f ′ (η) 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 K p retards the flow velocity. Because K p has ability to produce resistance for fluid flow. Mixed convection parameter effect on f ′ (η) for Williamson and Maxwell nanofluids is discussed in Fig. 3a. The greater values of mixed convection parameter ω causes an increment in f ′ (η) 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 f ′ (η) . 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  www.nature.com/scientificreports/ 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     www.nature.com/scientificreports/ 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 χ(η).

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:

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.