Natural bio-convective flow of Maxwell nanofluid over an exponentially stretching surface with slip effect and convective boundary condition.

The under-consideration article mainly focuses an unsteady three-dimensional Maxwell bio-convective nanomaterial liquid flow towards an exponentially expanding surface with the influence of chemical reaction slip condition. The feature of heat transport is achieving in the existenceof convective boundary condition and variable thermal conductivity. With the help of similarity variables, the flow form of equations is turned into a nonlinear form of coupled ODEs. The numerical solutions are calculated by adopting bvp4c function of MATLAB. Impact of distinct characteristics on the temperature, velocity microorganism and concentration field is graphically evaluated. Moreover, physical quantities are observed via graphs and tabulated data in details. It has been seen by the observation that the involvement of unsteadiness parameter restricts the change of laminar to turbulent flow. Further, for increasing velocity slip parameter velocity component in both directions shows lessening behavior. The Nusselt number exhibits diminishing behavior for larger values of Deborah number, and it shows the opposite behavior for larger values of convective parameter.

Introduction

To understand the rheological aspect and mechanism, many models have been established for non-Newtonian fluid in the past. The researchers have given a special attention to the nonlinear differential and rate type models. Maxwell model is the rate type model, and they discussed the characteristics of relaxation time. Due to specific application and special stress relaxation properties the non-Newtonian fluid is a talking point for the researchers. Non-Newtonian fluids are detected at chemical and nuclear industries, foodstuffs, bio engineering, polymeric liquids, and material processing. The Maxwell liquid model was proposed by Maxwell[1] to illustrate the elastic and viscous reaction of air. Zhao et al.[2] explored the Soret and Dufour impacts with Maxwell MHD liquid in porous surface. Zheng et al.[3] addressed Maxwell generalized liquid with oscillatory and constantly accelerating sheet. Shateyi[4] analyzed the MHD flow of Maxwell liquid in the occurrence of chemical reactions and thermophoresis effect on a vertical stretching sheet. Shafiq et al.[5] deliberated the MHD stagnation point flow of non-Newtonian (Walters-B) fluid across a stretching surface with homogenous-heterogenous reactions and Newtonian heating effect. Farooq et al.[6] work out the MHD Maxwell fluid flow with nanomaterials through an exponentially extending surface. Rasool et al.[7] designated the heat and mass transport investigation of Jeffrey MHD nanofluid flow with porous medium over an extending sheet. Khan and Nadeem[8] propose a comparison of linear and exponential stretching sheets of a rotating Maxwell nanomaterials liquid flow with stratification influence. Some latest research associated to non-Newtonian liquid is given in the Refs.[9-11].

The transfer of heat is a natural mechanism which happens with temperature differences within the system. Recently, the heat transfer phenomenon as a wave inspired the researchers from all over the world because heat transfer prevalent biomedical and industrial application. For example, electronic device cooling, power generator, heat conduction in tissues, and nuclear reactor cooling etc. The law of energy conduction to the analysis of transport of heat is suggested by Fourier[12]. Cattaneo[13] modifying the Fourier law to avoid the heat conduction behavior by exerting time relaxation term. Magyari and Keller[14] presented the mass and heat transport analysis of flow over an exponentially extending surface. The heat transfer of unsteady 3D viscous flow of a boundary layer fluid of the series solution passes an impulsively expanding sheet is employed by Xu et al.[15]. Kuznetsov and Nield[16] reviewed the characteristics of heat transport in natural convective flow of a Buongiorn’s model through a vertical plate. Gul et al.[17] presented the characteristic of transport of heat for the 2nd grade time dependent MHD thin film fluid flow analytically using two different ways. Hayat et al.[18] presented the heat transport analysis of stagnation point MHD flow on a vertical sheet. Gkountas et al.[19] analyzed the heat transfer of a viscous nanofluid in the presence of various nanoparticles.

The fluid dynamics by a stretching sheet are valuable in extrusion processes. The sheeting material formed in industrial production processes, and they consist of both polymer and metal sheets. The material region between the die and the collecting mechanism may logically assume that the stretching process alter with distance from the die, while cooling begins to stretch because of the solidification that ultimately happen. The current research concentrates to examine the flows by an exponentially stretching sheet. Such flow is quite widespread in applications such as paper production, crystal growing, continuous casting, glass fiber, metallurgical processes etc. The field of geophysical fluid dynamics that naturally occurring on earth is the main application of such fluid motion. The geophysical fluid dynamics contain a larger scale motion on earth, such as oceanography, meteorology, river flow, cloud’s motion etc. The extensibility of the sheet is a valuable aspect of the flow which can be carried out to boost the machinelike feature of the sheet. Flow on an extending sheet first time analyzed by Crane[20]. Later on, Gupta and Gupta[21] examined the mass and heat transport of liquid flow on an extending surface. Bidin and Nazar[22] incorporated the two-dimensional viscous liquid in the regime of radiation passes an exponentially extending surface. The flow of viscous liquid passes an exponentially extending surface with MHD is premeditated by Ishak[23]. Liu et al.[24] conferred the heat transport of 3D viscous fluid flow pass an exponentially surface. Hayat et al.[25] investigated the transport of heat on stagnant point MHD flow of nanoliquid across the extending surface in the presence of nonlinear radiation. Benos et al.[26] considered the shrinking / stretching surface to elaborate the transfer of heat on flow of MHD in the existence of radiation. Some studies concerning to stretching surfaces is presented in the Refs.[27-30].

Bioconvection is a phenomenon that is used to describe the instability and unstructured pattern formed due to the microorganisms, as a result the lesser density particles are swimming to the uppermost portion of a liquid. These complex microorganisms, such as gyrotactic microorganisms like algae, tend to cluster at the upper section of the fluid layer as they swim upwards, resulting in an unstable top heavy density stratification. Moreover, microorganisms are the microscopic organisms that lived everywhere in the surrounding such as deep sea, rocks, equator, deserts etc. The area of oil recovery and geophysical fact the bio-convection has a notable role. Kuznetsov[31] manifested the oxytactic microorganisms along similarity of finite depth shallow horizontal surface. The micropolar nanofluid with bio-convection recently suggested by Xu and Pop[32]. The mass and heat transport rate of convective flow of Nano liquid on a stretching sheet with microorganism is presented by Shafiq et al.[33]. Nadeem et al.[34] highlighted the 3D bio-convection nanomaterial liquid flow through an exponentially extending surface with micropolar fluid. Rashed and Nabwey[35] scrutinized the mixed bioconvection flow of nanomaterial liquid with convective conditions over a circular cylinder. Amer et al.[36] investigated the dynamical motion of a symmetric rigid body around a principal axis containing the viscous fluid in the existence of gyrostatic moment. Some recent study about gyrotactic microorganisms is found in the Refs.[37-39].

Motivation of the present work is to examine the three-dimensional bio-convective unsteady Maxwell nanomaterial liquid flow with the convective condition past an exponentially extending surface. The mass and heat transport investigation is represented with the influence of variable thermal conductivity and chemical reaction. The main finding of the current problem is to analyze the convective and concentration boundary condition together on the exponential stretching surface of a Maxwell nanofluid, which in not currently investigated in the literature yet. The transferred equations are tackled by applying bvp4c technique. Graphical outcomes of emerging characteristics are sketched and discussed. Physical behaviors of microorganisms, mass, and heat transport rate are analyzed through graphs and tabulated data.

Mathematical formulation

We studied an unsteady, incompressible, and three-dimensional flow of chemically reactive Maxwell bio-convective nanomaterials liquid towards an exponentially extending surface with The convective and slip boundary condition also taken into account to analyze the mass and heat transport. The flow attends the region shown in Fig. 1. Let the extending velocities of an exponentially stretching sheet is in the direction of and respectively. Inside the boundary layer , and denotes the nanoparticle volume concentration, temperature and microorganism density respectively. Furthermore, nanoparticle volume concentration, temperature and microorganism at the wall is defined by and respectively and away from the wall they are and respectively. Using above mentioned assumption the flow model takes the following form,

Figure 1

Physical picture of the paper.

Here is the thermal diffusivity, is the mass diffusivity, microorganism diffusivity, is the chemical reaction constant, is the variable thermal conductivity, is the body forces, is the chemotaxis constant, is the density of fluid, is the cell swimming speed, and is specific heat. The is the extra stress tensor for Maxwell liquid model, which is characterized as, where is the Rivlin-Ericksen tensor , is the viscosity, is the relaxation time, and is the material derivative. The governing equations of mass, momentum, energy, concentration, and microorganism by using the boundary layer approximation and Eqs. (1–6) takes the following form[40],

The related boundary conditions are[41],

The velocity components in and directions are and respectively. The symbols is the concentration slip factor, is the ratio between heat capacity of nanoparticles to the base fluid, and illustrates the heat transport coefficient. Furthermore, [42] is signified the variable thermal conductivity in which indicates the thermal conductivity of the surrounding.

To transform the flow model PDEs into non-dimensionalized form, we introduced the following non dimensional variable[43],

Here , , and all are the constants. Using Eq. (14), the equation of continuity automatically holds, and other Eqs. (8–12) take the following form,

The dimensionless form of the boundary conditions is,

Here prime stand for derivative with respect to . The symbols , , and are represented the unsteadiness parameter, Deborah number, Prandtl number, Schmidt number, stretching ratio characteristic, concentration slip parameter, thermal conductivity parameter and convection characteristic respectively. The symbols , , and denotes Brownian motion parameter, Peclet number, the thermophoresis parameter and bio-convection Schmidt number respectively. These parameters are defined as,

Here we take

Physical quantities

In terms of engineering, the physical quantities are the most important. These physical quantities observed the mass, heat, and microorganism transport rate. These quantities are defined as,

In Eq. (22) , shows the heat, microorganism and mass fluxes, respectively. These are given as,

The physical quantities in the dimensionless from are,

The local Reynolds number is .

Solution procedure

The nonlinear system of Eqs. (15–19) with Eq. (20) are solved numerically with the help of bvp4c MATLAB solution technique. To use a numerical method first we convert the Eqs. (15–20) into the system of first order differential equations. The convergence criteria were assigned as 10−6 [44,45]. Table 1 demonstrated the comparison of a temperature gradient with the previously published result of Nadeem et al.[44].

Table 1

Previous studies comparison of values, when .

	Nadeem et al.[44]	Presents results
Pr	θ^′(0)	θ^′(0)
0.72	0.809400	0.809401
1.0	1.000000	1.000000
3.0	1.923682	1.923683
10.0	3.720673	3.720674
100.0	12.29408	12.29409

The system of equations of first order,

The appropriative conditions are,

Physical description

The central aim of this work is to determine the Maxwell bio-convective nanomaterial liquid flow on an exponentially extending surface subject to the convective condition. Equations (9)–(14) are numerically manipulated by using bvp4c MATLAB technique. Further, graphically conclusions are conducted for different characteristic on the concentration, velocity, temperature and microorganism profile respectively. The parameter values is specified in the range of (0.0–0.5), (0.0–1.0), (0.0–1.0), (1.0–3.0), (0.1–1.5), (0.1–1.0), (1.0–4.0), (2.0–4.0), (0.0–1.0), (0.0–3.0), (0.0–1.5) and (1.0–4.0). In Table 2 the variation of physical quantities like Nusselt number, microorganism number and Sherwood number are observed. It is examined that the heat and microorganism transfer rate depict flourishing behavior by the boosting values of unsteadiness parameter and stretching parameter , but the mass transfer rate shows reverse trend. Further, when growing the amount of the Nusselt number, Deborah number and microorganism number displays lessening behavior, while the Sherwood number shows enlarging trend for . The tabulated data demonstrates that for larger amount of , the microorganism number and Sherwood number are decline, but heat transfer rate rises. Further for greater values of Brownian motion parameter both microorganism number and Sherwood number enhances, while the Nusselt number reduces for . It is clarified from Table 2 that, microorganism number and Sherwood number displays decreasing behavior for distinct values of thermophoresis characteristic . Moreover, it is demonstrated that different amount of Peclet number the tabulated data represents the reduction behavior for microorganism number.

Table 2

Table of , , and for different parameters when

A	λ	β	Pr	Nt	Nb	Pe	-θ^′0	-ϕ^′0	-χ^′0
0.0	0.5	0.3	2.5	0.3	0.3	1.5	0.3722	0.2898	1.600
0.2	–	–	–	–	–	–	0.4124	0.2327	1.8440
0.3	–	–	–	–	–	–	0.4249	0.2155	1.9620
0.3	0.5	0.3	–	0.3	–	–	0.4249	0.2155	1.9620
–	0.7	–	2.5	–	0.3	1.5	0.4279	0.2114	1.9970
–	0.9	–	–	–	–	–	0.4307	0.2076	2.030
0.3	0.5	0.0	–	0.3	–	–	0.4276	0.2119	1.9930
–	–	0.3	2.5	–	–	1.5	0.4249	0.2155	1.9620
0.3	0.5	0.5	–	–	0.3	–	0.4235	0.2174	1.9460
–	–	0.3	2.0	0.3	–	–	0.4119	0.2334	2.004
0.3	–	–	3.0	–	0.3	–	0.4350	0.2019	1.9280
–	–	–	4.0	–	–	1.5	0.4498	0.1820	1.8770
–	0.5	0.3	2.5	0.2	–	–	0.4283	0.1406	1.8620
0.3	–	–	–	0.4	–	–	0.4215	0.2936	2.0740
–	–	–	–	0.6	–	1.5	0.4145	0.4595	2.3340
–	0.5	0.3	–	0.3	0.1	–		0.6466	2.5910
0.3	–	–	2.5	–	0.3	–		0.2155	1.9620
–	–	–	–	–	0.5	–		0.1293	1.8520
–	–	–	–	–	0.3	1.0			1.8700
–	–	–	–	0.3	–	1.5			1.9620
0.2	0.1	0.3	0.3	0.5	0.0	2.0			2.0580

Figure 2 designates that the reduction in velocity field and is occurred by enlarging the values of unsteadiness parameter . Physically, when the amount of increases, the thickness of the boundary layer decreases significantly, and this development restricts the change of laminar to turbulent flow. Hence the flow is stabilized due to the stretching sheet. The Fig. 3 exhibits the diversion in and sketch for various amount of time relaxation characteristic . As seen in the sketched that the and are reduced as the increases. Physically is the ratio of relaxation time to observation time, as enlarging the the relaxation time also enlarges. Hence higher values of inhibit the fluid motion declines the and sketch. It is obvious that when , then viscous fluid is recovered. Further, when , then fluid is non-Newtonian. Figure 4 depict the influence of magnetic parameters on velocity profiles (, ). Form the figures it is observed that both the velocity field diminish with larger the estimation of magnetic parameter because the Lorentz force enhances which enlarge the resistive force in a fluid. As a result, velocity of fluid decreases while the temperature profile enhances (see in Fig. 5). Figure 5 examines the temperature distribution under the impact of the Eckert number. It is seen that the temperature field and their corresponding boundary layer thickness inclines. Figure 6 point out the fluctuation in temperature profile against the various estimation of the Prandtl number. Physically, as amplifying the reduction is occurred in the thermal diffusivity, hence the penetration depth of temperature decays. Further, also control the thermal boundary layers of fluid and relative momentum thickening. It is found in Fig. 6 that the temperature distribution is expanding due to climbing the values of surface convection parameter . Physically, is the proportion of the hot to colder fluid convection resistance. As increases the thermal resistance fall down and hence temperature raises. The impact of unsteadiness characteristic on the concentration and temperature distribution is deliberated in the Fig. 7. The reduction occurs in temperature and concentration plots by boosting the values of . The Fig. 8 examined the outcome of thermophoresis parameter and concentration slip parameter on the temperature sketch and concentration sketch respectively. It is seemed that the temperature sketch enlarged for greater values of , while concentration sketch declines for the higher values of . Physically, thermophoresis force occurs due to the temperature gradient, which leads to fast flow far off the sheet. The thermal boundary layer becomes dense by the large values of . The diversion in concentration sketch against the distinct values of Brownian motion parameter is visualized in Fig. 9. It is represented that the concentration sketch shows decreasing behavior for the greater values of . Figure 9 reported the variation in Schmidt number via concentration sketch. The sketch identified that the concentration sketch diminished by the enhancement of . Physically, the ratio between viscous diffusion to molecular diffusion rate is said to be Schmidt number. The reduction is noted in Schmidt number by rising the mass diffusivity, hence the concentration sketch shrinks for . The Fig. 10 found the impact of bio-convection Schmidt number and Peclet number via microorganism sketch. It is portrayed that the various values of and , the microorganism sketch shows decreasing behavior. Physically, have an inverse relation with microorganism diffusivity, as microorganism diffusivity boosts then the devaluation is occurred in the microorganism density profile. The amplification in unsteadiness parameter shows the diminishing behavior for microorganism sketch, which is illustrated in Fig. 11. The graphical description of Nusselt number, Sherwood number, and microorganism number for several parameters is presented in Fig. 12a–c. The Fig. 12a described that as growing the values of and the Nusselt number shows the increasing behavior in both cases. Figure 12b clearly shows that the Sherwood number decreases for distinct values of and . Figure 12c revealed the impact of microorganism number for distinct values of and . It is portrayed that the microorganism transfer rate enhances for the higher values of and . The heat and mass transfer rate shows the effectiveness of heat and mass convection at the surface.

Figure 2

Variation in and for .

Figure 3

Variation in and for .

Figure 4

Variation in and for .

Figure 5

Variation in for and .

Figure 6

Variation in for and .

Figure 7

Variation in and for .

Figure 8

Variation in for and for

Figure 9

Variation in for and

Figure 10

Variation in for and .

Figure 11

Variation in for and .

Figure 12

Variation of Nusselt number, Sherwood number and microorganism number for distinct parameters.

Concluding remarks

The 3D MHD flow of chemically reactive Maxwell bio-convective nanomaterial liquid embedded by an exponentially extending surface in the presence of viscous dissipation and joule heating effect. The thermal and solutal energy aspect has been addressed with the influence of Brownian motion and thermophoresis effect. The thermal convective and concentration slip boundary conditions are imposed on the boundary of the sheet. The specific observation of distinct characteristics is analyzed and summarized. The main outcomes of paper are:

By the enhancement of the magnetic characteristic and Deborah number the fluid velocity is declines due to occurrence of retardation effect.

By the increment of the unsteadiness parameter, the fluid velocity shows diminishing behavior.

The higher values of reduce the fluid temperature, while opposite trend is noted in the case of stronger .

The lager estimation of the Eckert number and magnetic characteristic boost the fluid temperature.

The effect of on temperature and concentration distribution are quantitatively similar, which is decreasing by enlarging the values of .

The concentration slip parameter shows declining behavior for higher values of . Further slip condition finished for .

The effects of and on microorganism distribution is qualitatively similar, which is decreasing as growing the values of and .

The tabulated results show that the Nusselt number and microorganism number increases, while Sherwood number exhibits declining behavior for higher amount of .

Sherwood and Nusselt number showing opposite behavior for distinct values of .