Investigation of binary chemical reaction in magnetohydrodynamic nanofluid flow with double stratification

Abstract

This article addresses MHD nanofluid flow induced by stretched surface. Heat transport features are elaborated by implementing double diffusive stratification. Chemically reactive species is implemented in order to explore the properties of nanofluid through Brownian motion and thermophoresis. Activation energy concept is utilized for nano liquid. Further zero mass flux is assumed at the sheet’s surface for better and high accuracy of the out-turn. Trasnformations are used to reconstruct the partial differential equations into ordinary differential equations. Homotopy analysis method is utilized to obtain the solution. Physical features like flow, heat and mass are elaborated through graphs. Thermal stratified parameter reduces the temperature as well as concentration profile. Also decay in concentration field is noticed for larger reaction rate parameter. Both temperature and concentration grows for Thermophoresis parameter. To check the heat transfer rate, graphical exposition of Nusselt number are also discussed and interpret. It is noticed that amount of heat transfer decreases with the increment in Hartmann number. Numerical results shows that drag force increased for enlarged Hartmann number.

Introduction

Having amazing heat transfer characteristics in contrast with usual heat transit-fluids, nanofluids with astonishing heat transit characteristics is the most discoursed topic of present time. Nanofluids contain nanoparticles with size under 100 nm. Over traditional heat transit fluids, nanofluids offer special influence. Sensational development in the thermal effects of host fluids is produced when a very slight quantity of nano-particles suspended thoroughly and dissipates constantly in the base fluids. To generate steady and highly conducted nano-fluids one step and two step methods are applied but creating nanoparticles with both methods experience cluster of nanoparticles. This is the pivotal issue in industrial science including nano-powder. Nanofluid was introduces by Choi and Eastman. ¹ After that different researchers analysed many properties of nanofluids. Evans et al. ² discussed the Brownian motion effect on nanofluid thermal conductivity. He investigated that thermal conductivity can be increased. Yacob et al. ³ presented nanofluid flow along with convective boundary condition passing through a stretching surface. Nanofluid flow accompanied by a nonlinear stretchable surface was elaborated by Nadeem and Lee. ⁴ Influence of viscous dissipatation of Copper-water nanofluid along with side by side plates was discussed by Sheikholeslami and Ganji. ⁵ Hatami and Ganji. ⁶ examined the MHD nanofluid flow with suction between parallel disks. Hayat et al.^7,8 discussed the thermal radiation effect of Powell-Eyring nanofluid and squeezing flow of carbon nanotubes. Ramana Reddy et al. ⁹ studied the inclined MHD unsteady flow of a nanofluid with Hall current and soret effects. They observed that enhanced soret number strengthens the momentum boundary layer thickness. Ayub et al. ¹⁰ inspected the nanofluid flow with MHD, slip effects and Riga plate. Impact of Brownian motion on MHD nanofluid over variable stretching sheet was discussed by Jayachandra Babu and Sandeep. ¹¹ It is found that velocity slip parameter diminishes the velocity field. Ramana Reddy et al. ¹² considered the thermophoresis and slip effects on MHD nanofluid. They used R.K. Fehlberg method for solving this problem and found that unsteadiness parameter decays the concentration profile. Anantha Kumar et al. ¹³ investigates about slip effects on bioconvective flow of nanofluids passing through a stretching sheet. Surface thickness parameter minimizes the temperature as well as concentration profile. Shah et al. ¹⁴ considered the Cattaneo-Christov heat flux model for second grade nanofluid having carbon nanotubes. He investigated that how entropy plays a role in our daily life. Impact of nonlinear radiation on MHD Casson fluid flow was studied by Ramudu et al. ¹⁵ They used numerical technique for solving differential equations. It is concluded that brownian motion and thermophoresis parameters strengthens the temperature profile. Consequences of Joule heating on copper and silver nanofluid flow with porous medium was discussed by Shah et al. ¹⁶ HAM and shooting methods are compared with the aid of numerical tables.

The process of stratification is illustrated as the formation of assorted layers possessing different densities. Stratification process occurs because of temperature changes, diversity of fluids and differences in concentration. Stratification is the interesting process in the mechanics of convective transportation because it has wide applications in the domains of industrial, natural and engineering processes. These applications include manufacturing processes, in atmosphere involving heterogeneous mixture, industrial food and salinity and thermal stratification mechanisms in oceans, rivers, reservoirs and reservoirs of ground water. Moreover, thermal stratification impedes penetration of oxygen between upper and lower layers of water and thus water be nominated as anoxic by the involvement of biological process and hence, this is illustrated as the disadvantage of it. RamReddy et al. ¹⁷ described the impact of non-Darcy porous mediumon nanofluid flow with thermal stratification. Hayat et al. ¹⁸ scrutinized the flow of an Oldroyd-B fluid with thermal stratification and stagnation point. Sheremet et al. ¹⁹ looked into the consequences of nanofluid filled in a square porous cavity under thermal stratification. Abbasi et al. ²⁰ reported the upshot of double stratification and radiation effects on Jeffrey nanofluid. Muhammad et al. ²¹ elaborated the features of thermal stratification in ferromagnetic fluid with stagnation point. Rehman et al. ²² flashes the characteristics of flat and cylindrical surfaces on tangent hyperbolic fluid with thermal radiations. Impact of thermal stratification and Joule heating on MHD nanofluid flow was considered by Daniel et al. ²³ Rehman et al. ²⁴ investigated the double stratification phenomenon on tangent hyperbolic fluid. Hayat et al. ²⁵ studied the characteristics of heat absorption in Oldroyd-B fluid with convective boundary conditions.

The process of mass transfer along with chemical reaction energy has been observed by many researchers due to its consequences in chemical engineering, geothermal reserviors, cooling of nuclear reactor and in oil recovery. Activation energy is the most important factor in chemical reaction. It is the fewer amount of energy for atoms or molecules due to which they experience a chemical reactionca. Bestman ²⁶ initially investigated the mass transfer of boundary layer flow in a porous medium with suction. Mass transfer flow of MHD thermally radiative fluid with activation energy was discussed by Maleque. ²⁷ Shafique et al. ²⁸ declared the behaviour of activation energy on Maxwell fluid. Influence of thermal reaction on Casson fluid with stagnation point was studied by Abbas et al. ²⁹ Zulfiqar et al. ³⁰ analysed the soret, dufour effects on MHD fluid flow with uniform suction/injection. Mustafa et al. ³¹ explored the features of chemical reaction on the MHD nanofluid flow passing through a vertical surface. Outcomes of chemical reaction on Cassonfluid flow over a rotating cone with Halll effect was inspected by Deebani et al. ³² He concluded that enhanced Brinkman number decays the Bejan number. Dawar et al. ³³ considered the Williamson nanofluid with activation energy. It is noted that Williamson parameter decays the drag force.

Magnetohydrodynamic nanofluid flows are very useful in engineering as well as in biomedicines. They have wide applications in nuclear power plant, nanofluid as a coolant, in fuels like Copper-oxide brake nanofluid (CBN), AOBN etc, magnetic nanoparticles for cancer therapy and many more. Aforeknown literature surveys designates that researchers highlighted the characteristics of Magnetohydrodynamic nanofluid flows. Although peculiarity of thermally stratified flows because of stretching sheet with chemical reaction have not been explored until now. Therefore our main objective is to investigate the behaviour of chemical reaction on MHD nanofluid flow over a stretching sheet. Thermal stratification is also contemplated. We considered variable temperature both at the surface and away from the wall. The emerging dominating equations are figure out by homotopic technique.^31,34–42 The consequences of demanding parameters are inspected through graphs. Additionally, drag force and heat transfer rate has been elaborated graphically.

Formulation

Two dimensional and steady state flow of MHD nanofluid passing through a stretching sheet is taken with thermal stratification. Energy and concentration equations are used to discuss the temperature and concentration profiles. Velocity is contemplated at the surface of the plate. Variable temperature is considered at wall that is, $T_w=T_0+\mathrm{bx}$ Similarly, surrounding temperature that is, $T_{\infty }=T_0+\mathrm{dx}$ is also taken variable. Concentration is assumed to be variable at surface and ambient fluid. After utilizing boundary layer approximations the ruling equations takes the succeeding forms: ³¹ see Figure 1 below.

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0,$$

(1)

$$u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}={\nu }_f\frac{{\partial }^{2}u}{\partial y^2}-\frac{{\sigma }^{*}}{{\rho }_f}{B}_0^{2}u,$$

(2)

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{\partial }^{2}T}{\partial y^2}+\tau [D_B\frac{\partial T}{\partial y}\frac{\partial C}{\partial y}+\frac{D_T}{T_{\infty }}{\left(\frac{\partial T}{\partial y}\right)}^2],$$

(3)

$$\begin{array}{c}\multicolumn{1}{c}{u\frac{\partial C}{\partial x}+v\frac{\partial C}{\partial y}=D_B\frac{{\partial }^{2}C}{\partial y^2}+\frac{D_T}{T_{\infty }}\frac{{\partial }^{2}T}{\partial y^2}}\\ -k_r^2(C-C_{\infty }){\left(\frac{T}{T_{\infty }}\right)}^n\exp (-\frac{E_a}{\mathrm{KT}}),\hfill \end{array}$$

(4)

Figure 1.

Physical model of the problem.

exposed boundary conditions are

$$\begin{array}{c}u=U_w\left(x\right)=\mathrm{ax},v=0,T=T_w=T_0+\mathrm{bx},\hfill \\ D_B\frac{\partial C}{\partial y}+\frac{D_T}{T_{\infty }}\frac{\partial T}{\partial y}=0\mathrm{at}y=0,\hfill \end{array}$$

(5)

$$\begin{array}{c}u\to \infty ,T\to T_{\infty }=T_0+\mathrm{dx},\hfill \\ C\to C_{\infty }=C_0+\mathrm{ex},\mathrm{at}y\to \infty .\hfill \end{array}$$

(6)

Here velocity components are denoted by $u$ and $v$ in $x-$ and $y-$ directions respectively, $U_w(x)$ $(=\mathrm{ax})$ is the stretching velocity, Magnetic field is denoted by $B_0$ , $T$ corresponds to fluid temperature, Wall temperature is $T_w(x)(=T_0+\mathrm{bx})$ , $T_{\infty }(x)(=T_0+\mathrm{dx})$ denotes the ambient fluid temperature, Brownian motion parameter is set as $D_B$ , $D_T$ figure out thermophoresis parameter, $C$ represents fluid concentration, $k_r$ exhibits reaction rate, $E_a$ corresponds to activation energy, $b,d$ and $e$ are the dimensional constants, Boltzmann constant is denoted by $K(=8.61\times {10}^{-5}\mathrm{eV}/K)$ , fitted rate constant is represented by $n$ and normally its range is $-1<n<1.$

Using

$$\begin{array}{c}u=\mathrm{axf}\prime (\eta ),v=-\sqrt{a\nu }f(\eta ),\eta =y\sqrt{\frac{c}{\nu }},\hfill \\ \theta \left(\eta \right)=\frac{T-T_{\infty }}{T_w-T_0},\mathrm{\Phi }\left(\eta \right)=\frac{C-C_{\infty }}{C_{\infty }},\hfill \end{array}$$

(7)

law of conservation of mass (1) is fulfilled consonantly. Therefore equations (2)–(6) becomes:

$$f″\prime +\mathrm{ff}″-(f\prime {)}^2-\mathrm{Mf}\prime =0,$$

(8)

$$\frac{1}{\mathrm{Pr}}\theta ″+f\theta \prime -f\prime \theta -S_{1}f\prime +N_b\theta \prime \mathrm{\Phi }\prime +N_t(\theta \prime {)}^2=0,$$

(9)

$$\begin{array}{c}\multicolumn{1}{c}{\mathrm{\Phi }″+\mathrm{Scf}\mathrm{\Phi }\prime +\left(\frac{N_t}{N_b}\right)\theta ″}\\ -\mathrm{Sc}\sigma {(1+\delta \theta )}^n\mathrm{\Phi }\exp \left(\frac{-E}{1+\delta \theta }\right)=0,\hfill \end{array}$$

(10)

$$f\prime (0)=1,f(0)=0,f\prime (\infty )=0,$$

(11)

$$\theta (0)=1-S_1,\theta (\infty )=0,$$

(12)

$$N_b\mathrm{\Phi }\prime \left(0\right)+N_t\theta \prime \left(0\right)=0,\mathrm{\Phi }\left(\infty \right)=0,$$

(13)

In the above equations Hartmann number is represented by $M$ , Stratification parameter is denoted by $S_1,\mathrm{Pr}$ denotes Prandtl number, $\mathrm{Sc}$ expresses the Schmidth number, $N_b$ flashes the Brownian motion parameter, $N_t$ stands for Thermophoresis parameter, $\sigma$ appears as dimensionless reaction rate, Non dimensional activation energy corresponds to $E,$ $\delta$ exhibits the temperature difference parameter. Mathematical form of these parameters are as follows:

$$\begin{array}{cc}\multicolumn{1}{c}{M=}& \frac{{\sigma }^{*}{B}_0^2}{{\rho }_{f}c},\mathrm{Pr}=\frac{v}{\alpha },S_1=\frac{d}{b},\mathrm{Sc}=\frac{v_f}{D_B},N_b=\frac{\tau D_B{C}_{\infty }}{v},\hfill \\ \multicolumn{1}{c}{N_t=}& \frac{\tau D_T(T_w-T_{\infty })}{v{T}_{\infty }},\sigma =\frac{k_r^2}{c},E=\frac{E_a}{T_{\infty }},\delta =\frac{(T_w-T_{\infty })}{T_{\infty }}.\hfill \end{array}$$

(14)

Mathematical expression of skin friction coefficient is disclosed as follows:

$$C_f=\frac{{\tau }_w}{\rho U_w^2},$$

(15)

Undimensioned mode is given as

$$C_f{\mathrm{Re}}_x^{-1/2}=f″(0),$$

(16)

Here, $R{e}_x=U_w(x)x/v$ denotes local Reynolds number.

Nusselt number for the present analysis is given as

$$\mathrm{Nu}{\mathrm{Re}}_x^{-1/2}=-\frac{\theta \prime \left(0\right)}{1-S_1}.$$

(17)

Homotopic solutions

Liao ³⁴ was the first who presented Homotopy analysis method in 1992. We get series solution of highly nonlinear problems with this method. This gave us great opportunity to choose initial guesses and linear operators. They are declared as given below:

$$f_0\left(\eta \right)=\left(1-\exp \left(-\eta \right)\right),$$

(18)

$${\theta }_0(\eta )=-\exp (-\eta )(-1+s_1),$$

(19)

$${\mathrm{\Phi }}_0(\eta )=\frac{\exp (-\eta )N_t(-1+s_1)}{N_b},$$

(20)

$$\begin{array}{c}{\mathbf{L}}_f\left(f\right)=\frac{d^{3}f}{d{\eta }^3}-\frac{\mathrm{df}}{d\eta },{\mathbf{L}}_{\theta }\left(\theta \right)=\frac{d^2\theta }{d{\eta }^2}-\theta ,\hfill \\ {\mathbf{L}}_{\mathrm{\Phi }}\left(\mathrm{\Phi }\right)=\frac{d^2\mathrm{\Phi }}{d{\eta }^2}-\mathrm{\Phi }=0,\hfill \end{array}$$

(21)

With

$${\mathbf{L}}_f[C_1+C_2\exp (\eta )+C_3\exp (-\eta )]=0,$$

(22)

$${\mathbf{L}}_{\theta }[C_4\exp (\eta )+C_5\exp (-\eta )]=0,$$

(23)

$$\mathbf{L}\mathrm{\Phi }[C_6\exp (\eta )+C_7\exp (-\eta )]=0,$$

(24)

here $C_i$ $(i=1,...,7)$ denotes the capricious constants.

The final solutions $(f_m,{\theta }_m,{\mathrm{\Phi }}_m)$ with the fusion of special solutions $(f_m^{*},{\theta }_m^{*},{\mathrm{\Phi }}_m^{?}\left(\eta \right))$ are expressed by

$$f_m\left(\eta \right)=f_m^{?}\left(\eta \right)+C_1+C_2{e}^{\eta }+C_3{e}^{-\eta },$$

(25)

$${\theta }_m\left(\eta \right)={\theta }_m^{?}\left(\eta \right)+C_4{e}^{\eta }+C_5{e}^{-\eta },$$

(26)

$${\mathrm{\Phi }}_m\left(\eta \right)={\mathrm{\Phi }}_m^{?}\left(\eta \right)+C_6{e}^{\eta }+C_7{e}^{-\eta },$$

(27)

Convergence analysis

Convergence of series solution can be secured by HAM. It depends over an auxiliary parameters ${\hslash }_f,$ ${\hslash }_{\theta }$ and ${\hslash }_{\mathrm{\Phi }}$ . Figures 2 and 3 declares the h-curves for velocity, temperature and mass equations. The permissible ranges of ${\hslash }_f$ , ${\hslash }_{\theta }$ and ${\hslash }_{\mathrm{\Phi }}$ are $-1.6\le {\hslash }_f\le -1.0,$ $-0.5\le {\hslash }_{\theta }\le -0.2$ and $-0.6\le {\hslash }_{\mathrm{\Phi }}\le -0.1.$

Figure 2.

$h$ -curve for $f\prime (\eta )$ .

Figure 3.

$h$ -curve for $\theta (\eta )$ and $\Phi (\eta )$ .

Discussion

This section contains graphical results of various opposite parameters for velocity, temperature and concentration distributions. Influence of Hartmann number $M$ on the velocity field is sketched in Figure 4. Decrement in velocity field is noticed for enlarged values of Hartmann number $M$ . Also momentum boundary layer thickness decays. It is because with enhanced Hatrmann number $M$ produces a large amount of Lorentz force which is subjected to reduction of velocity. The outcomes of Hartmann number on temperature field is displayed in Figure 5. Temperature profile grows clearly when Hartmann number $M$ grows. Thickness of thermal boundary layer also expands. The reason behind this is the production of Lorentz force. Maximum resistance produces because of Lorentz force due to which more heat will be generated and ultimately temperature profile raises. Analysis of Brownian motion parameter $N_b$ on the temperature field is illustrated in Figure 6. Both the temperature field and thermal boundary layer thickness grows with an increment of Brownian motion parameter $N_b.$ With the increase in $N_b$ collisions between fluid particles will increase. Due to this more heat will be generated. Figure 7 is plotted for Brownian motion parameter $N_b$ versus concentration field. It is remarked that concentration field shows decaying behaviour for Brownian motion parameter $N_b.$ With the increase in Brownian motion collisions between fluid particles will be maximum. Therefore small quantity of mass is relocated and thus downfall in concentration distribution is noticed. Figure 8 is plotted for noticing the consequences of thermophoresis parameter $N_t$ on temperature distribution. Intensified temperature profile is observed when Thermophoresis parameter $N_t$ increases. Increase in thermal boundary layer thickness is noticed too. Actually the particles near the heated plate have maximum temperature rather than the particles away from the plate. Therefore heated particles shifts heat to cold particles and hence temperature enhanced. Figure 9 is portrayed to display the response of Thermophoresis parameter $N_t$ on concentration profile. Both concentration field and solutal boundary layer thickness grows when Thermophoresis parameter $N_t$ expands. Physically, thermophoresis force grows for enhanced thermophoresis parameter. So more heat transfer will occur and leads to higher diffusive effects. Figure 10 describes the significance of thermal stratification parameter $S_1$ on the temperature field. Decrease in thermal stratification parameter is noticed when temperature field increases. Lower region has maximum density for enlarged thermal stratification parameter $S_1.$ As a outcome, heated wall produces resistance in the flow towards the surrounding wall. Consequences of thermal stratification parameter $S_1$ on concentration field is exposed in Figure 11. Concentration field is detected to increase with the enhanced thermal stratification parameter $S_1$ . Both thermal and solutal boundary layer thickness decays when thermal stratification parameter $S_1$ increases. With the increase in $S_1$ density difference increases between the layers of fluid particles. This density difference creates a barrier for mass transfer. Hence decrement occurs in concentration field. Figure 12 reflects the change in Prandtl number $\mathrm{Pr}$ versus temperature field. Here we examine that temperature field decreases when Prandtl number $\mathrm{Pr}$ grows. Thermal diffusivity decreases as Prandtl number enhances and temperature field reduces. Figure 13 flashes the effect of dimensionless reaction rate parameter $\sigma$ on concentration field. It is esteemed that concentration field decreases for growing dimensionless reaction rate parameter $\sigma$ . It is revealed that the factor $\sigma (1+\delta \theta {)}^n\exp (\frac{-E}{1+\delta \theta })$ enhances for dominant $\sigma .$ Therefore, concentration gradient enhances at the wall. Hence concentration profile diminishes. Figure 14 reflects the deviation of Schmidth number $\mathrm{Sc}$ on concentration field. It is noticed that concentration profile raises for larger values of Schmidth number. Momentum diffusivity enhanced for increment in Schmidth number $\mathrm{Sc}.$ Therefore concentration distribution raises. Impact of wall temperature difference parameter $\delta$ on the concentration field is depicted in Figure 15. Concentration field shows increasing trend when temperature difference parameter $\delta$ gradually increases. Dominant values of $\delta$ result in higher energetic fluid particles which consequently grows the mass transfer. Thus concentration field enhances. Figures 16 and 17 are plotted for numerous values of activation energy $E$ on temperature and concentration distributions. It is noticeable that temperature distribution enhances and concentration profile decays for dominant activation energy $E$ . Enhanced activation energy results in destructive chemical reaction due to which temperature raises and concentration decreases. Figures 18 and 19 are plotted to check the heat transfer rate corresponding to Hartmann number $M,$ thermophoresis parameter $N_t,$ Prandtl number $\mathrm{Pr}$ and stratification parameter $S_1$ . It is perceived that heat transfer rate reduces for dominant Hartmann number $M$ and it grows for higher thermophoresis parameter $N_t$ . Physically, Lorentz force produces maximum resistance. Due to this temperature field increases and hence less heat will be transferred to the environment. Therefore heat transfer rate diminishes. Similarly heat transfer rate increases for growing Prandtl number $\mathrm{Pr}$ and stratification parameter $S_1.$ This is because enhanced Prandtl number and stratification parameter weakens the temperature field. So maximum heat will be transferred from fluid to environment. Therefore, heat transfer rate increases.

Figure 4.

Description of $M$ on $f\prime (\eta ).$

Figure 5.

Description of $M$ on $\theta (\eta ).$

Figure 6.

Description of $N_b$ on $\theta (\eta ).$

Figure 7.

Description of $N_b$ on $\mathrm{\Phi }(\eta ).$

Figure 8.

Description of $N_t$ on $\theta (\eta ).$

Figure 9.

Description of $N_t$ on $\mathrm{\Phi }(\eta ).$

Figure 10.

Description of $S_1$ on $\theta (\eta ).$

Figure 11.

Description of $S_1$ on $\mathrm{\Phi }(\eta ).$

Figure 12.

Description of $\mathrm{Pr}$ on $\theta (\eta ).$

Figure 13.

Description of $\sigma$ on $\mathrm{\Phi }(\eta ).$

Figure 14.

Description of $\mathrm{Sc}$ on $\mathrm{\Phi }(\eta ).$

Figure 15.

Description of $\delta$ on $\mathrm{\Phi }(\eta ).$

Figure 16.

Description of $E$ on $\theta (\eta ).$

Figure 17.

Description of $E$ on $\mathrm{\Phi }(\eta ).$

Figure 18.

Impact of $M$ and $N_t$ on $\mathrm{Nu}{\mathrm{Re}}^{\frac{-1}{2}}.$

Figure 19.

Impact of $\mathrm{Pr}$ and $S_1$ on $\mathrm{Nu}{\mathrm{Re}}^{\frac{-1}{2}}.$

Table 1 is constructed for the inspection of skin friction $(f″(0))$ with published work of Hayat et al. ⁴³ in limiting case. Both results are in good manner.

Table 1.

Comparison of drag force for distinct values of Hartman number $M$ when $\mathrm{Pr}=1,$ $N_b$ = $N_t$ = $\mathrm{Sc}$ = $S_1$ = $E$ = $n$ = $\delta$ = $\sigma$ = $0.$

$M$	Hayat et al. 43	Present
0	1.00000	1.00000
1	−1.41421	−1.41421
5	−2.44948	−2.44949
10	−3.31662	−3.31662

Concluding remarks

In this article we analyse the characteristics of chemical reaction of MHD nanofluid with thermal stratification. The main outcomes are summarized as:

• Increasing Hartmann number $M$ diminishes the velocity field.

• Dominant Brownian motion $N_b$ parameter results in increment of temperature profile.

• Higher Prandtl number $\mathrm{Pr}$ decays temperature profile.

• Thermal stratified parameter $S_1$ reduces the temperature as well as concentration profiles.

• Concentration field decays for growing reaction rate parameter $\sigma$ .

• Temperature difference parameter $\delta$ strengthens the concentration field.

• Enlarged Hartmann number $M$ reduces the heat transfer rate.

It is hoped that the present study subsidized as a motivation for representing supplementary Magnetohydrodynamic nanofluid flows mainly in nuclear reactors, accelerators, flow metres, MHD generators, designing of cooling systems and biomedicines. This article may be used in nano drug delivery, cancer therapy, cell separation, automotive industry and solar energy harvesting.