Effects of induced magnetic field and homogeneous–heterogeneous reactions on stagnation flow of a Casson fluid

==Abstract==

In this study, we analyzed the induced magnetic field effect on the stagnation-point flow of a non-Newtonian fluid over a stretching sheet with homogeneous–heterogeneous reactions and non-uniform heat source or sink. The transformed ordinary differential equations are solved numerically using Runge–Kutta and Newtons method. For physical relevance we analyzed the behavior of homogeneous and heterogeneous profiles individually in the presence of induced magnetic field. The effects of different non-dimensional governing parameters on velocity, induced magnetic field, temperature and concentration profiles, along with the skin friction coefficient and local Nusselt number, are discussed and presented through graphs. The results of the present study are validated by comparing with the existed literature. Results indicate that induced magnetic field parameter and stretching ratio parameter have the tendency to enhance the heat transfer rate.

==Keywords==

Induced magnetic field ; Homogeneous–heterogeneous reactions ; Non-uniform heat source/sink ; Casson fluid ; Stagnation flow

==Nomenclature==

<math display="inline">u,\mbox{ }v</math>- Velocity components in <math display="inline">x\mbox{ }\mbox{and}\mbox{ }y</math>  directions                          

<math display="inline">I_1\mbox{},\mbox{ }I_2</math>- Magnetic components in <math display="inline">x\mbox{ }\mbox{and}\mbox{ }y</math>  directions                          

<math display="inline">I_0</math>- Magnetic component in free stream flow

<math display="inline">x</math>- Distance along the surface

<math display="inline">y</math>- Distance normal to the surface

<math display="inline">c_p</math>- Specific heat capacity at constant pressure

<math display="inline">f</math>- Dimensionless velocity

<math display="inline">u_e</math>- Free stream velocity

<math display="inline">T</math>- Temperature of the fluid

<math display="inline">q^{{'´´}}</math>- Dimensionless non-uniform heat source/sink

<math display="inline">A\ast </math> ,<math display="inline">B\ast </math>- Non-uniform heat source/sink

<math display="inline">C_f</math>- Skin friction coefficient

<math display="inline">Nu</math>- Local Nusselt number

<math display="inline">{Re}_x</math>- Local Reynolds number

<math display="inline">Pr</math>- Prandtl number

<math display="inline">a,\mbox{ }d</math>- Positive constants

<math display="inline">k</math>- Thermal conductivity

<math display="inline">K</math>- Homogeneous reaction

<math display="inline">K_s</math>- Heterogeneous reaction

===Greek symbols===

<math display="inline">a,\mbox{ }b</math>- Rate constants

<math display="inline">\eta </math>- Similarity variable

<math display="inline">\sigma </math>- Electrical conductivity

<math display="inline">\theta </math>- Dimensionless temperature

<math display="inline">\rho </math>- Density

<math display="inline">\beta </math>- Casson parameter

<math display="inline">\lambda </math>- Reciprocal magnetic Prandtl number

<math display="inline">\Lambda </math>- Magnetic parameter

<math display="inline">{\mu }_e</math>- Magnetic permeability

<math display="inline">\mu </math>- Dynamic viscosity

<math display="inline">\upsilon </math>- Kinematic viscosity

<math display="inline">D_A\mbox{},\mbox{ }D_B</math>- Diffusion coefficients

<math display="inline">\alpha </math>- Thermal diffusivity

<math display="inline">\delta </math>- <math display="inline">a/c</math>  stretching ratio parameter                          

===Subscripts===

<math display="inline">f</math>- Fluid

<math display="inline">w</math>- Condition at the wall

<math display="inline">\infty </math>- Condition at the free stream

==1. Introduction==

The study of convection flows over stretching sheet has extensive applications in the field of manufacturing production of plastic, polythene, paper, polymer extrusion, cooling of elastic sheets, and fiber technology, science and engineering technology. The applications of stretching surface were broadly discussed by Chaim [[#bib0010|[1]]] . The flow of a non-Newtonian fluid over a stretching sheet has potential applications in the biosciences, blood flows, jelly and engineering. Casson is a shear thinning liquid. It can exhibit yield stress. If applied yield stress is greater than the shear stress, then it acts as solid, whereas if yield stress is lesser than the applied shear stress then it starts to move. The application of Casson fluid was discussed by Hayat et al. [[#bib0015|[2]]] . The homogeneous–heterogeneous reactions on stagnation-point flow past a stretching surface was discussed by Bachok et al. [[#bib0020|[3]]] . Peristaltic flow of a Carreau fluid flow with convective boundary conditions was analyzed by Hayat et al. [[#bib0025|[4]]] , and concluded that brinkman number and hall parameters help to enhance the thermal boundary layer thickness. Khan and Pop [[#bib0030|[5]]]  investigated the viscoelastic fluid flow through a stretching surface with homogeneous–heterogeneous effects. Hayat et al. [[#bib0035|[6]]]  illustrated the three-dimensional MHD bidirectional flow of a nanofluid past a permeable stretching sheet in the presence of slip condition and homogeneous–heterogeneous reactions.      

Heat transfer analysis of magneto hydrodynamic stagnation-point flow toward a stretching surface in the presence of induced magnetic field was done by Ali et al. [[#bib0040|[7]]] . Abbas et al. [[#bib0045|[8]]]  considered the slip conditions and homogeneous–heterogeneous effects on MHD stagnation-point flow of viscous fluid past a stretching or shrinking surface. Sandeep and Sulochana [[#bib0050|[9]]]  discussed the mixed convection flow of an unsteady micropolar fluid flow over a stretching or shrinking sheet in the presence of magnetic and non-uniform heat source/sink. Raju et al. [[#bib0055|[10]]]  analyzed the non-uniform heat source or sink effects on ferrofluid flow through a flat plate with thermal radiation and aligned magnetic field effects. The stagnation-point flow toward a vertical stretching surface was explained by Ishak et al. [[#bib0060|[11]]] . Mallikarjuna et al. [[#bib0065|[12]]]  investigated the variable porosity regime on convection flow over a vertical cone with chemical reaction and magnetic field effect. An unsteady MHD nanofluid flow through a stretching surface in the presence of non-uniform heat source/sink and magnetic field effect was studied by Sandeep et al. [[#bib0070|[13]]] . Pal and Mandal [[#bib0075|[14]]]  analyzed the induced magnetic field effect on stagnation point flow of a nanofluid past a non-isothermal stretching surface.      

Heat transfer analysis of a stagnation point viscoelastic fluid flow toward a stretching sheet was investigated by Mahapatra and Gupta [[#bib0080|[15]]] . Gorla [[#bib0085|[16]]]  depicted the stagnation point flow of non-Newtonian fluid flow in the presence of transverse magnetic field. Heat transfer characteristics of MHD Casson fluid flow past a permeable an exponentially stretching sheet was considered by Raju et al. [[#bib0090|[17]]] . Hayat et al. [[#bib0095|[18]]]  analyzed the two-dimensional flow of Eyring–Powell fluid flow toward a stretching sheet with homogeneous–heterogeneous reactions and concluded that the Eyring–Powell fluid material parameter enhances the boundary layer thickness. Heat transfer characteristics of Casson fluid flow past a permeable exponentially stretching sheet in the presence of thermal radiation was discussed by Pramanik [[#bib0100|[19]]]  and concluded that an increasing value of Casson parameter suppresses the velocity profiles. Animasaun [[#bib0105|[20]]]  discussed the MHD Casson fluid flow in the presence of viscous dissipation. Stagnation point flow of optimized exact solution for oblique flow of Casson-nanofluid with convective boundary conditions was discussed by Nadeem et al. [[#bib0110|[21]]] . Makinde and Aziz [[#bib0115|[22]]]  and [[#bib0120|[23]]]  examined the boundary layer analysis of a nanofluid flow through a stretching surface in the presence of convective boundary conditions and buoyancy effects. MHD heat transfer characteristics of an electrically conducting nanofluid flow through nonlinearly stretching surface in the presence of viscous dissipation was examined numerically by Mabood et al. [[#bib0125|[24]]] . Nadeem and Saleem [[#bib0130|[25]]]  proposed an analytical solution for an unsteady mixed convection MHD flow on a rotating cone. Dissipation and magnetic field effects on the boundary layer flow of a power-law nanofluid induced by a permeable stretching/shrinking sheet were investigated numerically by Dhanai et al. [[#bib0135|[26]]] . Chaudhary and Merkin [[#bib0140|[27]]]  proposed a simple isothermal model for homogeneous–heterogeneous reactions for boundary layer flow. Very recently, Hayat et al. [[#bib0145|[28]]]  and Mahanta and Shaw [[#bib0150|[29]]]  investigated the heat transfer characteristics of Casson fluid through different channels. Hayat et al. [[#bib0155|[30]]]  studied the homogeneous and heterogeneous reactions on the flow over a nanotube with homogeneous heating. Further, Farooq et al. [[#bib0160|[31]]]  discussed the homogeneous–heterogeneous characteristics on the flow of Jeffrey fluid.      

In this study, we investigated the induced magnetic field and non-uniform heat source or sink effects on the stagnation-point flow of a non-Newtonian fluid toward a stretching sheet in the presence of homogeneous–heterogeneous reactions. The emerging set of governing nonlinear partial differential equations is transformed into a set of ordinary differential equations using similarity transformation, which are then solved numerically using Runge–Kutta and Newtons method. The effects of different non-dimensional governing parameters on velocity, induced magnetic field, temperature and concentration profiles, along with the friction factor and local Nusselt number, are discussed and presented through graphs and tables.

==2. Flow analysis==

Consider a steady, incompressible, electrically conducting stagnation point flow of a Casson fluid toward a stretching sheet in the presence of induced magnetic field, non-uniform heat source/sink and homogeneous–heterogeneous reactions. The stretching sheet is considered along the <math display="inline">x-</math> axis and <math display="inline">y-</math> axis is normal to it as displayed in [[#f0010|Fig. 1]] .

<span id='f0010'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0001.jpg|center|395px|Physical model and coordinate system.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 1.

Physical model and coordinate system.

</span>
|}

The governing boundary layer equations for the flow and induced magnetic field profiles are given by [[#bib0040|[7]]] .

<span id='e0010'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=</math><math>0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 1)
|}

<span id='e0015'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\frac{\partial I_1}{\partial x}+\frac{\partial I_2}{\partial y}=</math><math>0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 2)
|}

<span id='e0020'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\rho \left(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}\right)=</math><math>u_e\mbox{ }\left(x\right)\frac{\partial u_e\mbox{ }\left(x\right)}{\partial x}+</math><math>\nu \left(1+\frac{1}{\beta }\right)\frac{{\partial }^2u}{\partial y^2}</math>
|-
|<math>\mbox{ }+\mbox{ }\frac{\mu }{4\pi }\left(I_1\frac{\partial I_1}{\partial x}+\right. </math><math>\left. I_2\frac{\partial I_1}{\partial y}-I_e\frac{\partial I_e}{\partial x}\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 3)
|}

<span id='e0025'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u\frac{\partial I_1}{\partial x}+v\frac{\partial I_1}{\partial y}=</math><math>I_1\frac{\partial u}{\partial x}+I_2\frac{\partial u}{\partial y}+</math><math>{\mu }_e\frac{{\partial }^2I_1}{\partial y^2},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 4)
|}

Subjected to the boundary conditions

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
u=u_w\mbox{ }\left(x\right)=cx,\mbox{ }v=0,\mbox{ }\frac{\partial I_1}{\partial y}=I_2=0,\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{at}\mbox{ }y=0,\\
u=u_e\mbox{ }\left(x\right)=ax,\mbox{ }\frac{\partial u}{\partial y}\rightarrow 0,\mbox{ }\frac{\partial v}{\partial y}\rightarrow 0,\mbox{ }I_1=I_e\mbox{ }\left(x\right)=I_0x,\mbox{ }\mbox{as}\mbox{ }y\rightarrow \infty ,
\end{array}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 5)
|}

where <math display="inline">{\mu }_e</math>  is the magnetic diffusivity of the fluid, which is given by <math display="inline">{\mu }_e=1/4\pi \sigma </math> , <math display="inline">u</math>  and <math display="inline">v</math>  are the velocity components along the <math display="inline">x</math>  and <math display="inline">y</math> directions, respectively, <math display="inline">\rho </math>  is the density, and <math display="inline">\nu </math>  is the viscosity.      

To convert the governing equations of the flow into a set of nonlinear ordinary differential equations, we now introduce the following similarity transformation:

<span id='e0035'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
u=cxf^{{'}}\left(\eta \right),\mbox{ }v=-{\upsilon }_f{}^{1/2}c^{1/2}f\left(\eta \right),\mbox{ }\eta ={\upsilon }_f{}^{-1/2}c^{1/2}y,\\
I_1=I_0xg^{{'}}\left(\eta \right),\mbox{ }I_2=-I_0{\upsilon }_f^{1/2}c^{-1/2}g\left(\eta \right),
\end{array}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 6)
|}

Equation [[#e0035|(6)]]  identically satisfies the continuity equations [[#e0010|(1)]]  and [[#e0015|(2)]] , and equations [[#e0020|(3)]]  and [[#e0025|(4)]]  become

<span id='e0040'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left(1+\frac{1}{\beta }\right)f^{{'''}}-\left({f^{{'}}}^2-\right. </math><math>\left. ff^{{''}}\right)+\Lambda \left({g^{{'}}}^2-\right. </math><math>\left. gg^{{''}}-1\right)+{\left(\frac{a}{c}\right)}^2=</math><math>0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 7)
|}

<span id='e0045'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\lambda g^{{'''}}+fg^{{''}}-f^{{''}}g=0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 8)
|}

Subjected to the transformed boundary conditions

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
f=g=0,\mbox{ }f^{{'}}=1,\mbox{ }g^{{''}}=0,\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{at}\mbox{ }\eta =0,\\
f^{{'}}=a/c,\mbox{ }g^{{'}}=1,\mbox{ }g^{{''}}=f^{{''}}=0,\mbox{ }\mbox{as}\mbox{ }\eta \rightarrow \infty ,
\end{array}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 9)
|}

where <math display="inline">\beta </math>  is the Casson parameter, <math display="inline">\Lambda </math>  is the magnetic parameter, and <math display="inline">\lambda </math>  is the reciprocal magnetic Prandtl number, which are given by <math display="inline">\Lambda =\frac{\mu H_0^2}{4\pi \rho c^2},\mbox{ }\lambda =</math><math>\frac{1}{4\pi \sigma \upsilon }</math> .      

==3. Heat transfer analysis==

The governing equation for energy in the presence of no-uniform heat source/sink is given by

<span id='e0055'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}\right)=</math><math>\alpha \frac{{\partial }^2T}{\partial y^2}+\frac{1}{\rho c_p}q^{{'''}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 10)
|}

With the boundary conditions

<span id='e0060'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>T=T_w\mbox{},\mbox{ }\mbox{at}\mbox{ }y=0,\mbox{ }T=</math><math>T_{\infty }\mbox{},\mbox{ }\mbox{as}\mbox{ }y\rightarrow \infty ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 11)
|}

The space and temperature-dependent heat generation/absorption (non-uniform heat source/sink) <math display="inline">q^{{'´´}}</math> is defined by Sandeep and Sulochana [[#bib0050|[9]]] :

<span id='e0065'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>q^{{'''}}=\left(\frac{ku_w\mbox{ }\left(x\right)}{x\upsilon }\right)\left(A\ast \left(T_w-\right. \right. </math><math>\left. \left. T_{\infty }\right)f^{{'}}+B\ast \left(T-\right. \right. </math><math>\left. \left. T_{\infty }\right)\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 12)
|}

where ''T  ''  is the fluid temperature, <math display="inline">\alpha </math>  is the thermal diffusivity coefficient, <math display="inline">\upsilon </math>  is the kinematic viscosity of the base fluid, and <math display="inline">A\ast </math> and <math display="inline">B\ast </math> are parameters of the space and temperature-dependent internal heat generation/absorption. The positive and negative values of <math display="inline">A\ast </math>  and <math display="inline">B\ast </math> represent heat generation and absorption parameters respectively. Now introducing the similarity transformation

<span id='e0070'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\eta ={\upsilon }_f{}^{-1/2}c^{1/2}y,\mbox{ }\theta \left(\eta \right)=</math><math>\left(T-T_{\infty }\right)/\left(T_w-T_{\infty }\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 13)
|}

By using equation [[#e0070|(13)]] , equations [[#e0055|(10)]] , [[#e0060|(11)]]  and [[#e0065|(12)]]  transformed to

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\theta }^{{''}}+\left(A\ast f^{{'}}+B\ast \theta \right)+</math><math>Prf{\theta }^{{'}}=0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 14)
|}

With the transformed boundary conditions

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\theta =1,\mbox{ }\mbox{at}\mbox{ }\eta =0,\mbox{ }\theta =</math><math>0,\mbox{ }\mbox{as}\mbox{ }\eta \rightarrow \infty ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 15)
|}

where <math display="inline">Pr=\frac{v}{\alpha }</math> is the Prandtl number.      

==4. Mass transfer analysis==

We considered a simple homogeneous–heterogeneous reaction model that exists as proposed by Chaudhary and Merkin [[#bib0140|[27]]]  in the following form:

<span id='e0085'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left(u\frac{\partial a}{\partial x}+v\frac{\partial a}{\partial y}\right)=</math><math>D_A\frac{{\partial }^2a}{\partial y^2}-k_cab^2\mbox{},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 16)
|}

<span id='e0090'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left(u\frac{\partial b}{\partial x}+v\frac{\partial b}{\partial y}\right)=</math><math>D_B\frac{{\partial }^2b}{\partial y^2}+k_cab^2\mbox{},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 17)
|}

Subject to the boundary conditions

<span id='e0095'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
D_A\frac{\partial a}{\partial y}=k_sa\left(y\right),\mbox{ }D_B\frac{\partial b}{\partial y}=-k_sa\left(y\right)\mbox{ }\mbox{at}\mbox{ }y=0,\\
a\left(y\right)=a_0\mbox{},\mbox{ }b\left(y\right)=0\mbox{ }\mbox{as}\mbox{ }y\rightarrow \infty ,
\end{array}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 18)
|}

where<math display="inline">D_A\mbox{ }\mbox{and}\mbox{ }D_B</math>  are the diffusion coefficients. We now introduced the similarity transformations as

<span id='e0100'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>G\left(\eta \right)=a/a_0\mbox{},\mbox{ }H\left(\eta \right)=</math><math>b/a_0\mbox{},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 19)
|}

Using equation [[#e0100|(19)]] , equations [[#e0085|(16)]] , [[#e0090|(17)]]  and [[#e0095|(18)]]  transformed to

<span id='e0105'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>G^{{''}}+Sc\left(fG^{{'}}-KGH^2\right)=0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 20)
|}

<span id='e0110'></span>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\delta H^{{''}}+Sc\left(fH^{{'}}+KGH^2\right)=0,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 21)
|}

With the transformed boundary conditions

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
G^{{'}}=K_sG,\mbox{ }\delta H^{{'}}=-K_sG\mbox{ }\mbox{at}\mbox{ }\eta =0,\\
G=1,\mbox{ }H=0\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{as}\mbox{ }\eta \rightarrow \infty ,
\end{array}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 22)
|}

<math display="inline">Sc</math>  is the Schmidt number, <math display="inline">K</math>  is the measure of strength of homogeneous reaction, <math display="inline">K_s</math>  is the strength of heterogeneous reaction, <math display="inline">Re=c/\nu </math>  is the Reynolds number and <math display="inline">\delta </math>  is the ratio of diffusion coefficient, which are represented below.

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\delta =\frac{D_B}{D_A},\mbox{ }Sc=\frac{\nu }{D_A},\mbox{ }K=</math><math>\frac{k_ca_0{}^2}{c},\mbox{ }K_s=\frac{k_s}{D_A\sqrt{Re}}\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 23)
|}

In most of the applications, we expect the diffusion coefficient of chemical species A and B to be of comparable size. This leads us to make an assumption that the diffusion coefficients D<sub>A</sub>  and D<sub>B</sub>  are equal. That is the case to <math display="inline">\delta =1</math> .      

For engineering interest, the shear stress coefficient <math display="inline">C_f</math>  (friction factor) and the local Nusselt number <math display="inline">Nu_x</math>  are given by

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>C_f=\frac{{\tau }_w}{{\rho }_fU_w^2},\mbox{ }\mbox{ }Nu_x=</math><math>\frac{xq_w}{k_f\mbox{ }\left(T_w-T_{\infty }\right)},\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 24)
|}

Here, <math display="inline">{\tau }_w</math>  is the wall shear stress and <math display="inline">q_w</math> is the wall heat flux, which are given by

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\tau }_w=\mu {\left(\frac{\partial u}{\partial y}\right)}_{y=0}\mbox{},\mbox{ }q_w=</math><math>-k{\left(\frac{\partial T}{\partial y}\right)}_{y=0}\mbox{},\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 25)
|}

Using (6) and (13) we get

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{Re}_x^{1/2}C_f=\left(1+\frac{1}{\beta }\right)f^{{''}}\left(0\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 26)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{Re}_x^{-1/2}Nu_x=-{\theta }^{{'}}\left(0\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 27)
|}

where <math display="inline">{Re}_x=u_wx/v_f</math>  is the local Reynolds number.      

==5. Results and discussion==

The set of nonlinear ordinary differential equations [[#e0040|(7)]] , [[#e0045|(8)]] , [[#e0055|(10)]] , [[#e0105|(20)]]  and [[#e0110|(21)]]  subjected to the boundary conditions (9), (11) and (22) are solved numerically using Runge–Kutta and Newtons method (Mallikarjuna et al. [[#bib0065|[12]]] ). The influence of non-dimensional governing parameters on velocity, induced magnetic field, temperature and concentration profiles, along with the friction factor coefficient and local Nusselt number, are discussed and presented through graphs and tables. For numerical calculations, we considered <math display="inline">K=K_S=1,\mbox{ }a/c=0.5,\mbox{ }\delta =</math><math>3,\mbox{ }\lambda =0.2,\mbox{ }Pr=1</math> , <math display="inline">\Lambda =A\ast =B\ast =0.1,\mbox{ }Sc=</math><math>0.6</math> . In the entire study we kept these values as common except the variations in the corresponding graphs and tables. In this paper green color indicates the Newtonian fluid case and red color indicates the non-Newtonian fluid case.      

[[#f0015|Fig. 2]] , [[#f0020|Fig. 3]] , [[#f0025|Fig. 4]]  and [[#f0030|Fig. 5]]  depict the influence of stretching ratio parameter on velocity, induced magnetic field and concentration profiles for both Newtonian and non-Newtonian fluid cases. It is evident that the velocity and homogeneous concentration profiles are encouraged with increasing values of stretching ratio parameter, and reverse profiles are observed in induced magnetic field and heterogeneous concentration profiles. Basically, an increase in the stretching ratio parameter initially develops the more pressure on the flow; due to this reason, we have seen an enhancement in the velocity profiles.

<span id='f0015'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0002.jpg|center|546px|Velocity profiles for different values of stretching ratio parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 2.

Velocity profiles for different values of stretching ratio parameter.

</span>
|}

<span id='f0020'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0003.jpg|center|546px|Induced magnetic field profiles for different values of stretching ratio ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 3.

Induced magnetic field profiles for different values of stretching ratio parameter.

</span>
|}

<span id='f0025'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0004.jpg|center|546px|Concentration profiles for different values of stretching ratio parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 4.

Concentration profiles for different values of stretching ratio parameter.

</span>
|}

<span id='f0030'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0005.jpg|center|546px|Concentration profiles for different values of stretching ratio parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 5.

Concentration profiles for different values of stretching ratio parameter.

</span>
|}

The effects of induced magnetic field parameter on velocity, induced magnetic field and concentration profiles for both Newtonian and non-Newtonian cases are plotted in [[#f0035|Fig. 6]] , [[#f0040|Fig. 7]] , [[#f0045|Fig. 8]]  and [[#f0050|Fig. 9]] . It is noticed that the momentum, thermal and concentration boundary layer thicknesses are enhanced with an increase in the induced magnetic field parameter. Usually an increase in the induced magnetic field develops the electric current. This electric force can help to enhance the momentum and thermal boundary layer thickness. This leads to an increase in the momentum and thermal boundary layer thickness.

<span id='f0035'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0006.jpg|center|546px|Velocity profiles for different values of induced magnetic field parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 6.

Velocity profiles for different values of induced magnetic field parameter.

</span>
|}

<span id='f0040'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0007.jpg|center|546px|Induced velocity profiles for different values of induced magnetic field ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 7.

Induced velocity profiles for different values of induced magnetic field parameter.

</span>
|}

<span id='f0045'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0008.jpg|center|546px|Concentration profiles for different values of induced magnetic field parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 8.

Concentration profiles for different values of induced magnetic field parameter.

</span>
|}

<span id='f0050'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0009.jpg|center|546px|Concentration profiles for different values of induced magnetic field parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 9.

Concentration profiles for different values of induced magnetic field parameter.

</span>
|}

The influence of homogeneous reaction parameter on concentration profiles for both Newtonian and non-Newtonian cases is depicted in [[#f0055|Fig. 10]]  and [[#f0060|Fig. 11]] . It displays that an increasing value of <math display="inline">K</math> enhances the <math display="inline">H</math> concentration profiles and depreciates the<math display="inline">G</math>  concentration profiles. It may be due to the fact that the reaction rates dominate diffusion coefficients. Quite similar results are displayed in [[#f0065|Fig. 12]]  and [[#f0070|Fig. 13]] , with an increase in the heterogeneous parameter. This agrees with the general physical behavior of homogeneous and heterogeneous reactions.

<span id='f0055'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0010.jpg|center|546px|Concentration profiles for different values of strength of homogeneous reaction ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 10.

Concentration profiles for different values of strength of homogeneous reaction parameter.

</span>
|}

<span id='f0060'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0011.jpg|center|546px|Concentration profiles for different values of strength of homogeneous reaction ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 11.

Concentration profiles for different values of strength of homogeneous reaction parameter.

</span>
|}

<span id='f0065'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0012.jpg|center|546px|Concentration profiles for different values of heterogeneous reaction parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 12.

Concentration profiles for different values of heterogeneous reaction parameter.

</span>
|}

<span id='f0070'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0013.jpg|center|546px|Concentration profiles for different values of strength of heterogeneous ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 13.

Concentration profiles for different values of strength of heterogeneous reaction parameter.

</span>
|}

The effect of the ratio of the diffusion coefficient on the concentration profiles is depicted in [[#f0075|Fig. 14]] . Mixed profiles in the concentration are observed; initially, less diffused particles are developed in the concentration. Due to this reason, initially we have noticed depreciation in the concentration profiles. [[#f0080|Fig. 15]]  and [[#f0085|Fig. 16]]  illustrate the effect of non-uniform heat source or sink parameter on temperature profiles. Generally, higher values of non-uniform heat source or sink parameter enhance temperature profiles. Because it generates an internal energy in the flow, we have seen an improvement in temperature profiles.

<span id='f0075'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0014.jpg|center|546px|Concentration profiles for different values of ratio of diffusion parameter.]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 14.

Concentration profiles for different values of ratio of diffusion parameter.

</span>
|}

<span id='f0080'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0015.jpg|center|546px|Temperature profiles for different values of non-uniform heat source/sink ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 15.

Temperature profiles for different values of non-uniform heat source/sink parameter.

</span>
|}

<span id='f0085'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_920191755-1-s2.0-S221509861530121X-jestch199-fig-0016.jpg|center|546px|Temperature profiles for different values of non-uniform heat source/sink ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 16.

Temperature profiles for different values of non-uniform heat source/sink parameter.

</span>
|}

[[#t0010|Table 1]]  and [[#t0015|Table 2]]  depict the comparison of the present results with previously published work under some special limited cases. We found an excellent agreement of the present results with existing results. This proves the validity of the present results along with the accuracy of the numerical technique we used in this study. [[#t0020|Table 3]]  and [[#t0025|Table 4]]  presented the effect of various emerging thermo physical parameters for skin friction coefficients and local Nusselt number for both Newtonian and non-Newtonian fluid cases. It is evident that the friction factor coefficients and heat transfer rate are enhanced with an increase in the induced magnetic field parameter. The quite similar types of results are displayed in stretching ratio parameter. It is also observed that an increase in the non-uniform heat source or sink parameter declines the rate of heat transfer.

<span id='t0010'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Table 1.

Comparison of the skin friction coefficient when <math display="inline">K_s=K=\Lambda =K=A\ast =B\ast =</math><math>Sc=0\mbox{ }and\mbox{ }\beta \rightarrow \infty </math> .                  

|-

!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>a/c</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


! Ishak et al. [[#bib0060<nowiki>|</nowiki>[11]]]
! Mahapatra and Gupta [[#bib0080<nowiki>|</nowiki>[15]]]
! Ali et al. [[#bib0040<nowiki>|</nowiki>[7]]]
! Present
|-

| 0.1
| −0.9694
| −0.9694
| −0.9694
| −0.96941
|-

| 0.2
| −0.9181
| −0.9181
| −0.9181
| −0.91810
|-

| 0.5
| −0.6673
| −0.6673
| −0.6673
| −0.66731
|-

| 2
| 2.0175
| 2.0175
| 2.0175
| 2.01750
|-

| 3
| 4.7294
| 4.7293
| 4.7293
| 4.72931
|}

<span id='t0015'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Table 2.

Comparison of local Nusselt number when <math display="inline">K=K_s=Sc=\Lambda =K=A\ast =</math><math>B\ast =0,\mbox{ }\beta \rightarrow \infty </math> .                  

|-

! rowspan="3" |

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>a/c</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


! colspan="8" |

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>Pr</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


|-

! colspan="4" | Mahapatra and Gupta [[#bib0080<nowiki>|</nowiki>[15]]]
! colspan="4" | Present results
|-

! 0.05
! 0.5
! 1
! 1.5
! 0.05
! 0.5
! 1
! 1.5
|-

| 0.1
| −0.081
| −0.383
| −0.603
| −0.777
| −0.0811
| −0.3833
| −0.6031
| −0.7772
|-

| 0.2
| −0.099
| −0.408
| −0.625
| −0.797
| −0.0991
| −0.4082
| −0.6252
| −0.7971
|-

| 0.5
| −0.136
| −0.473
| −0.692
| −0.863
| −0.1360
| −0.4731
| −0.6921
| −0.8632
|-

| 1
| −0.178
| −0.563
| −0.796
| −0.974
| −0.1782
| −0.5632
| −0.7962
| −0.9741
|-

| 2
| −0.241
| −0.709
| −0.974
| −1.171
| −0.2413
| −0.7094
| −0.9743
| −1.1713
|-

| 3
| −0.289
| −0.829
| −1.124
| −1.341
| −0.2891
| −0.8290
| −1.1241
| −1.3411
|}

<span id='t0020'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Table 3.

The physical parameter values for Newtonian fluid case.

|-

!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>a/c</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\Lambda </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>K</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>K_s</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>A\ast </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>B\ast </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\delta </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f^{{''}}\left(0\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\theta }^{{'}}\left(0\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


|-

| 0.5
| 
| 
| 
| 
| 
| 
| −0.510405
| 1.731575
|-

| 0.6
| 
| 
| 
| 
| 
| 
| −0.458589
| 1.739244
|-

| 0.7
| 
| 
| 
| 
| 
| 
| −0.371571
| 1.753433
|-

| 
| 0.1
| 
| 
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 0.15
| 
| 
| 
| 
| 
| −0.435305
| 1.739196
|-

| 
| 0.2
| 
| 
| 
| 
| 
| −0.321719
| 1.760058
|-

| 
| 
| 1
| 
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 1.4
| 
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 1.8
| 
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 
| 0.4
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 
| 0.6
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 
| 0.8
| 
| 
| 
| −0.516500
| 1.723825
|-

| 
| 
| 
| 
| 0.5
| 
| 
| −0.520685
| 1.268812
|-

| 
| 
| 
| 
| 1
| 
| 
| −0.520685
| 0.707954
|-

| 
| 
| 
| 
| 1.5
| 
| 
| −0.520685
| 0.147097
|-

| 
| 
| 
| 
| 
| 1
| 
| −0.475514
| 1.521289
|-

| 
| 
| 
| 
| 
| 2
| 
| −0.475514
| 1.216227
|-

| 
| 
| 
| 
| 
| 3
| 
| −0.475514
| 0.860203
|-

| 
| 
| 
| 
| 
| 
| 1.5
| −0.516500
| 1.723825
|-

| 
| 
| 
| 
| 
| 
| 2
| −0.516500
| 1.723825
|-

| 
| 
| 
| 
| 
| 
| 2.5
| −0.516500
| 1.723825
|}

<span id='t0025'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Table 4.

The physical parameter values in non-Newtonian fluid case.

|-

!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>a/c</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\Lambda </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>K</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>K_s</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>A\ast </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>B\ast </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\delta </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f^{{''}}\left(0\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


!

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>-{\theta }^{{'}}\left(0\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | 
|}


|-

| 0.5
| 
| 
| 
| 
| 
| 
| −0.305742
| 1.766651
|-

| 0.6
| 
| 
| 
| 
| 
| 
| −0.277062
| 1.769630
|-

| 0.7
| 
| 
| 
| 
| 
| 
| −0.225497
| 1.777041
|-

| 
| 0.1
| 
| 
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 0.15
| 
| 
| 
| 
| 
| −0.253464
| 1.769873
|-

| 
| 0.2
| 
| 
| 
| 
| 
| −0.183574
| 1.782195
|-

| 
| 
| 1
| 
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 
| 1.4
| 
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 
| 1.8
| 
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 
| 
| 0.4
| 
| 
| 
| −0.305683
| 1.760186
|-

| 
| 
| 
| 0.6
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 
| 
| 0.8
| 
| 
| 
| −0.305683
| 1.760185
|-

| 
| 
| 
| 
| 0.5
| 
| 
| −0.305679
| 1.297576
|-

| 
| 
| 
| 
| 1
| 
| 
| −0.305679
| 0.726128
|-

| 
| 
| 
| 
| 1.5
| 
| 
| −0.305679
| 0.154680
|-

| 
| 
| 
| 
| 
| 1
| 
| −0.328580
| 1.551278
|-

| 
| 
| 
| 
| 
| 2
| 
| −0.328580
| 1.255117
|-

| 
| 
| 
| 
| 
| 3
| 
| −0.328580
| 0.913042
|-

| 
| 
| 
| 
| 
| 
| 1.5
| −0.305683
| 1.760185
|-

| 
| 
| 
| 
| 
| 
| 2
| −0.305683
| 1.760185
|-

| 
| 
| 
| 
| 
| 
| 2.5
| −0.305683
| 1.760185
|}

==6. Conclusions==

This study deals with the stagnation-point flow of a Casson fluid past a stretching sheet in the presence of induced magnetic field and homogeneous–heterogeneous reactions. The non-uniform heat source/sink is also taken into account. The arising sets of governing nonlinear coupled partial differential equations are transformed into a set of nonlinear coupled ordinary differential equations using self-similarity transformations, which are then solved numerically. The conclusions are as follows:
* Heat transfer performance of Newtonian fluids was comparatively high when compared with the transfer performance of non-Newtonian fluids.
* The non-uniform heat source/sink parameter enhances the temperature profiles for both Newtonian and non-Newtonian fluid cases.
* Stretching ratio parameter improves the local Nusselt number and friction factor coefficients.
* An increase in the induced magnetic field parameter enhances the heat transfer rate.
* A homogeneous–heterogeneous reaction does not show significant variations in the heat transfer rate.

==References==

<ol style='list-style-type: none;margin-left: 0px;'><li><span id='bib0010'></span>
[[#bib0010|[1]]] T.C. Chiam; Micropolar fluid flow over a stretching sheet; Zeitschrift für Angewandte Mathematik und Mechanik, 62 (1982), pp. 565–568</li>
<li><span id='bib0015'></span>
[[#bib0015|[2]]] T. Hayat, M. Awais, M. Sajid; Mass transfer effects on the unsteady flow of UCM fluid over a stretching sheet; Int. J. Mod. Phys. B, 25 (2011), pp. 2863–2878</li>
<li><span id='bib0020'></span>
[[#bib0020|[3]]] N. Bachok, A. Ishak, I. Pop; On the stagnation-point flow towards a stretching sheet with homogeneous-heterogeneous reactions effects; Commun. Nonlinear Sci. Numer. Simul, 16 (2011), pp. 4296–4302</li>
<li><span id='bib0025'></span>
[[#bib0025|[4]]] T. Hayat, A. Tanveer, H. Yasmin, A. Alsaedi; Homogeneous-heterogeneous reactions in peristaltic flow with convective conditions; PLoS ONE, 9 (12) (2014), p. e113851</li>
<li><span id='bib0030'></span>
[[#bib0030|[5]]] W.A. Khan, I.M. Pop; Effects of homogeneous-heterogeneous reactions on the viscoelastic fluid towards a stretching sheet; J. Heat Transfer, 134 (2012), pp. 1–5</li>
<li><span id='bib0035'></span>
[[#bib0035|[6]]] T. Hayat, M. Imtiaz, A. Alsaedi; Impact of magneto hydrodynamics in bidirectional flow of nanofluid subject to second order slip velocity and heterogeneous-homogeneous reactions; J. Magn. Mater, 395 (2015), pp. 294–302</li>
<li><span id='bib0040'></span>
[[#bib0040|[7]]] F.M. Ali, R. Nazar, N.M. Arfin, I. Pop; MHD stagnation-point flow and heat transfer towards stretching sheet with induced magnetic field; Appl. Math. Mech. Engl. Ed, 32 (4) (2011), pp. 409–418</li>
<li><span id='bib0045'></span>
[[#bib0045|[8]]] Z. Abbas, M. Shekh, I. Pop; Stagnation-point flow of hydromagnetic viscous fluid over stretching/shrinking sheet with generalized slip condition in the presence of homogeneous-heterogeneous reactions; J. Taiwan Inst. Chem. Eng, 55 (2015), pp. 69–75</li>
<li><span id='bib0050'></span>
[[#bib0050|[9]]] N. Sandeep, C. Sulochana; Dual solutions for unsteady mixed flow of MHD micropolar fluid over a stretching/shrinking sheet with non-uniform heat source/sink; Eng. Sci. Technol. Int. J., 8 (2015), pp. 1–8</li>
<li><span id='bib0055'></span>
[[#bib0055|[10]]] C.S.K. Raju, N. Sandeep, C. Sulochana, V. Sugunamma; Effects of aligned magnetic field and radiation on the flow of ferrofluids over a flat plate with non-uniform heat source/sink; Int. J. Sci. Eng, 8 (2) (2015), pp. 151–158</li>
<li><span id='bib0060'></span>
[[#bib0060|[11]]] A. Ishak, R. Nazar, I. Pop; Mixed convection boundary layers in the stagnation-point flow towards a vertical stretching sheet; Meccanica, 41 (2006), pp. 509–518</li>
<li><span id='bib0065'></span>
[[#bib0065|[12]]] B. Mallikarjuna, A.M. Rashad, A.J. Chamka, S. Hariprasad Raju; Chemical reaction effects on MHD convective heat and mass transfer flow past a rotating vertical cone embedded in a variable porosity regime; Afr. Math (2015) [http://dx.doi.org/10.1007/s13370-015-0372-1 http://dx.doi.org/10.1007/s13370-015-0372-1]</li>
<li><span id='bib0070'></span>
[[#bib0070|[13]]] N. Sandeep, C. Sulochana, C.S.K. Raju, M. Jayachandra Babu, V. Sugunamma; Unsteady boundary layer flow of thermophoretic MHD nanofluid past a stretching sheet with space and time dependent internal heat source/sink; Appl. Appl. Math, 10 (1) (2015), pp. 312–327</li>
<li><span id='bib0075'></span>
[[#bib0075|[14]]] D. Pal, G. Mandal; MHD convective stagnation-point flow of nanofluids over a non-isothermal stretching sheet with induced magnetic field; Meccanica, 50 (8) (2015), pp. 2023–2035</li>
<li><span id='bib0080'></span>
[[#bib0080|[15]]] T.R. Mahapatra, A.S. Gupta; Stagnation-point flow of a viscoelastic fluid towards a stretching surface; Int. J. Non-Linear Mech, 39 (2004), pp. 811–820</li>
<li><span id='bib0085'></span>
[[#bib0085|[16]]] R.S.R. Gorla; Non-Newtonian fluid at a stagnation point in presence of a transverse magnetic field; Mech. Res. Commun, 3 (1976), pp. 1–6</li>
<li><span id='bib0090'></span>
[[#bib0090|[17]]] C.S.K. Raju, N. Sandeep, V. Sugunamma, M. Jayachandrababu, J.V. Ramanareddy; Heat and mass transfer in magneto hydrodynamic Casson fluid over an exponentially permeable stretching surface; Eng. Sci. Technol. Int. J. (2015) [http://dx.doi.org/10.1016/j.jestch.2015.05.010 http://dx.doi.org/10.1016/j.jestch.2015.05.010]</li>
<li><span id='bib0095'></span>
[[#bib0095|[18]]] T. Hayat, M. Imtiaz, A. Alsaedi; Effects of homogeneous-heterogeneous reactions in flow of Powell-Eyring fluid; J. Cent. South Univ, 22 (2015), pp. 3211–3216</li>
<li><span id='bib0100'></span>
[[#bib0100|[19]]] S. Pramanik; Casson fluid flow and Heat transfer past an exponentially porous stretching sheet in presence of thermal radiation; Ain Shams Eng. J., 5 (2014), pp. 205–212</li>
<li><span id='bib0105'></span>
[[#bib0105|[20]]] I.L. Animasaun; Effects of thermophoresis, variable viscosity and thermal conductivity on free convective heat and mass transfer of non-Darcian MHD dissipative Casson fluid flow with suction and nth order of chemical reaction; J. Niger. Math. Soc, 34 (2015), pp. 11–31</li>
<li><span id='bib0110'></span>
[[#bib0110|[21]]] S. Nadeem, R. Mehmood, N.S. Akbar; Optimized analytical solution for oblique flow of a Casson-nanofluid with convective boundary conditions; Int. J. Therm. Sci, 78 (2014), pp. 90–100</li>
<li><span id='bib0115'></span>
[[#bib0115|[22]]] O.D. Makinde, A. Aziz; Boundary layer flow of a nanofluid past a stretching sheet with convective boundary condition; Int. J. Therm. Sci, 50 (2011), pp. 1326–1332</li>
<li><span id='bib0120'></span>
[[#bib0120|[23]]] O.D. Makinde, W.A. Khan, Z.H. Khan; Buoyancy effects on MHD stagnation point flow and heat transfer of a nanofluid past a convectively heated stretching/shrinking sheet; Int. J. Heat Mass Transfer, 62 (2013), pp. 526–533</li>
<li><span id='bib0125'></span>
[[#bib0125|[24]]] F. Mabood, W.A. Khan, A.I.M. Ismail; MHD boundary layer flow and heat transfer of nanofluids over a non-linear stretching sheet: a numerical study; J. Magn. Magn. Mater, 374 (2015), pp. 569–576</li>
<li><span id='bib0130'></span>
[[#bib0130|[25]]] S. Nadeem, S. Saleem; Analytical treatment of unsteady mixed convection MHD flow on a rotating cone in a rotating frame; J. Taiwan Inst. Chem. Eng, 44 (2013), pp. 596–604</li>
<li><span id='bib0135'></span>
[[#bib0135|[26]]] R. Dhanai, P. Rana, L. Kumar; Multiple solutions of MHD boundary layer flow and heat transfer behavior of nanofluids induced by a power-law stretching/shrinking permeable sheet with viscous dissipation; Power Technol, 273 (2015), pp. 62–70</li>
<li><span id='bib0140'></span>
[[#bib0140|[27]]] M.A. Chaudhary, J.H. Merkin; A simple isothermal model for homogeneous-heterogeneous reactions in boundary layer flow: I. Equal diffusivities; Fluid Dyn. Res, 16 (1995), pp. 311–333</li>
<li><span id='bib0145'></span>
[[#bib0145|[28]]] T. Hayat, M. Farooq, A. Alsaedi; Thermally stratified stagnation point flow of Casson fluid with slip conditions; Int. J. Numer. Methods Heat Fluid Flow, 25 (4) (2015), pp. 724–748</li>
<li><span id='bib0150'></span>
[[#bib0150|[29]]] G. Mahanta, S. Shaw; 3D Casson fluid flow past a porous linearly stretching sheet with convective boundary condition; Alex. Eng. J., 54 (3) (2015), pp. 653–659</li>
<li><span id='bib0155'></span>
[[#bib0155|[30]]] T. Hayat, M. Farooq, A. Alsaedi; Homogeneous-heterogeneous reactions in the stagnation point flow of carbon nano tubes with Newtonian heating; AIP Adv, 5 (2015), p. 027130 [http://dx.doi.org/10.1063/1.4908602 http://dx.doi.org/10.1063/1.4908602]</li>
<li><span id='bib0160'></span>
[[#bib0160|[31]]] M. Farooq, A. Alsaedi, T. Hayat; Note on characteristics of homogeneous-heterogeneous reaction in flow of Jeffrey fluid; Appl. Math. Mech. Engl. Ed, 36 (10) (2015), pp. 1319–1328</li>
</ol>