Skip to main content

MHD mixed convection slip flow near a stagnation-point on a nonlinearly vertical stretching sheet

Abstract

The problem of magnetohydrodynamic (MHD) mixed convection flow near a stagnation-point region over a nonlinear stretching sheet with velocity slip and prescribed surface heat flux is investigated; this has not been studied before. Using a similarity transformation, the governing equations are transformed into a system of ordinary differential equations, and then are solved by employing a homotopy analysis method. The effects of the nonlinearity parameter, the magnetic field, mixed convection, suction/injection, and the boundary slip on the velocity and temperature profile are analyzed and discussed. The results reveal that the increasing exponent of the power-law stretching velocity increases the heat transfer rate at the surface. It is also found that the velocity slip and magnetic field increase the heat transfer rate when the free stream velocity exceeds the stretching velocity, i.e. \(\varepsilon< 1\), and they suppress the heat transfer rate for \(\varepsilon> 1\).

1 Introduction

The problem of stagnation-point flow and heat transfer on stretching sheet arises in an abundance of practical applications in industry and engineering, such as cooling of electronic devices and nuclear reactors, polymer extrusion, drawing of plastic sheets, etc.; and, moreover, in the magnetohydrodynamic (MHD) flow which has both liquid and magnetic properties and can exhibit particular characteristics in thermal conductivity. Thus the study of MHD stagnation-point flow on stretching sheet has attracted many researchers in recent times, and many problems are discussed as regards different aspects, including the stretching sheet with variable surface temperature [1] or viscous dissipation [2, 3], the effect of slip [4, 5], and the analysis of the unsteady case [6].

Different from the flow induced by a stretching horizontal sheet, the effect of mixed convection due to a buoyancy force could not be neglected for the vertical sheet. There has been increasing interest in studying the problem of MHD with mixed convection boundary layer flow and heat transfer characteristics over a stretching vertical surface [715]. Very recently, Ali et al. [16] studied the MHD mixed convection stagnation-point flow and heat transfer of an incompressible viscous fluid over a vertical stretching sheet, and the MHD boundary layer flow over a vertical stretching/shrinking sheet in a nano-fluid was investigated by Makinde et al. [17] and Das et al. [18].

The above investigations considered the flow on the linearly stretching sheet or vertical sheet, but the real stretching velocity does not always need to be linear or uniform. Some work has been done in this field. The similarity solution of the boundary layer equations for a nonlinearly stretching sheet has been found by Akyildiz et al. [19]. The flow and heat transfer over a nonlinearly stretching sheet has been investigated by Akyildiz and Siginer [20] by using a Legendre spectral method. Recently, Dhanai et al. [21], Ashraf et al. [22], and Mabood et al. [23] analyzed the boundary layer flow and heat transfer on a nonlinearly shrinking/stretching sheet immersed in a nanofluid.

In the present paper, motivated by the above studies, the problem of MHD mixed convection stagnation-point flow on the nonlinearly vertical stretching sheet is discussed in the presence of buoyancy force, suction/injection parameters, and boundary slip. The governing nonlinear coupled partial differential equations are reduced to a set of ordinary differential equations by means of similarity transformations. The reduced equations are solved by the homotopy analysis method (HAM) [24], which has been successfully applied to various interesting complicated fluid problems [2530]. Graphs are plotted to gain physical insight toward the key embedding physical parameters. To the best of our knowledge, the series solutions for this model have not been presented before.

2 Mathematical formulation of the problem

We consider the steady two-dimensional MHD mixed convection flow in the vicinity of a stagnation point over a nonlinear stretching sheet with velocity slip and prescribed surface heat flux. The flow is confined to the region \(y\geq0\), where y is the coordinate measured normal to the stretching surface. A uniform magnetic field of strength is applied in the direction normal to the surface \(y=0\). The flow model along with the coordinate system is shown in Figure 1. It is assumed that the sheet stretching velocity \(u_{w}(x) =cx^{m}\) and the external velocity is prescribed as \(u_{e}(x)= ax^{m}\) where c and a are positive constants. The constant m is the nonlinearity parameter, with \(m=1\) for the linear case and \(m \neq1\) for the nonlinear case. Under the boundary layer approximation and the assumptions that the viscous dissipation and Joule heating are neglected, the basic equations of continuity, momentum, and energy are given by

$$\begin{aligned}& \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0, \end{aligned}$$
(1)
$$\begin{aligned}& u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=u_{e}\frac {d u_{e}}{d x}+\nu \frac{\partial^{2} u}{\partial y^{2}}+\frac{\sigma B^{2}(x)}{\rho} (u_{e}-u)+g\beta(T-T_{\infty}), \end{aligned}$$
(2)
$$\begin{aligned}& u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha\frac{\partial^{2}T}{\partial y^{2}}, \end{aligned}$$
(3)

where u and v are the velocities in the x and y directions, ν is the kinematic viscosity, ρ is the fluid density, σ is the electrical conductivity, \(B(x)\) is the transverse magnetic field, g is the acceleration due to gravity, β is the thermal expansion coefficient, T is the fluid temperature, and α is the thermal diffusivity. The relevant boundary conditions are given by

$$\begin{aligned}& u=u_{w}(x)+\frac{2-\sigma_{v}}{\sigma_{v}}\lambda_{0} \frac{\partial u}{\partial y}, \qquad v=v_{w}(x),\qquad \frac{\partial T}{\partial y}=- \frac{q_{w} (x)}{k}, \quad \mbox{at } y=0, \end{aligned}$$
(4)
$$\begin{aligned}& u\rightarrow u_{e}(x), \qquad T \rightarrow T_{\infty}\quad \mbox{as } y \rightarrow\infty, \end{aligned}$$
(5)

where \(\sigma_{v} \) is the tangential momentum accommodation coefficient, \(\lambda_{0}\) is the mean free path, \(v_{w}(x)\) is the suction (injection) velocity, k is the thermal conductivity, and \(q_{w}(x)\) is the surface heat flux.

Figure 1
figure 1

Physical model and coordinate system.

We introduce the following similarity transformations:

$$ \left \{ \begin{array}{l} \eta=\sqrt{\frac{a}{\nu}}yx^{(m-1)/2}, \\ \psi=\sqrt{a \nu}x^{(m+1)/2} f(\eta), \\ \theta=\sqrt{\frac{a}{\nu}}\frac{k(T-T_{\infty})}{q_{0}x^{2m-1}}. \end{array} \right . $$
(6)

Here ψ is the stream function. Equation (1) is satisfied by introducing ψ such that \(u=\partial\psi/\partial y\) and \(v=-\partial\psi /\partial x\). Employing the similarity variables (6), the velocity components u and v are given by

$$ u=ax^{m}f'(\eta) \quad \mbox{and} \quad v=-\sqrt{a \nu}x^{(m-1)/2} \biggl[\frac{m+1}{2}f(\eta)+ \frac{m-1}{2} \eta f'(\eta) \biggr], $$
(7)

where the prime denotes differentiation with respect to η. To obtain similarity solutions, \(B(x)\), \(v_{w}(x)\), and \(q_{w}(x)\) are taken as

$$ \begin{aligned} &B(x)=B_{0}x^{(m-1)/2},\qquad v_{w}= -\frac{\sqrt{a\nu }(m+1)}{2}x^{(m-1)/2}S, \\ &q_{w}(x)=q_{0}x^{(5m-3)/2}, \end{aligned} $$
(8)

where \(B_{0}\), S, and \(q_{0}\) are constants. It is noted that \(S>0\) corresponds to the injection case and \(S<0\) implies suction. Substituting (6) into (2) and (3), we get the following ordinary differential equations:

$$\begin{aligned}& f'''+\frac{m+1}{2}ff''+m \bigl(1-f^{\prime2}\bigr)+M\bigl(1-f'\bigr)+ \lambda\theta=0, \end{aligned}$$
(9)
$$\begin{aligned}& \theta''+\frac{P_{r}(m+1)}{2}f\theta'-P_{r}(2m-1)f' \theta=0, \end{aligned}$$
(10)

subject to the boundary conditions (4) and (5), which become

$$\begin{aligned}& f(0)=S,\qquad f'(0)=\varepsilon+\delta f''(0), \qquad \theta'(0)=-1, \end{aligned}$$
(11)
$$\begin{aligned}& f'(\infty)=1,\qquad \theta(\infty)=0. \end{aligned}$$
(12)

In the above equations, M is the magnetic parameter, λ is the mixed convection parameter, and \(P_{r}\) is the Prandtl number, and they are given by

$$ M=\frac{\sigma B_{0}^{2}}{a\rho}, \qquad \lambda=\frac{g\beta q_{0}\sqrt{\nu}}{ka^{5/2}}= \frac{\mathit{Gr}_{x}}{\mathit{Re}_{x}^{5/2}} \quad \mbox{and} \quad P_{r}=\frac {\nu}{\alpha} $$
(13)

with \(\mathit{Gr}_{x}=g\beta q_{w}x^{4}/k\nu^{2}\) and \(\mathit{Re}_{x}=u_{e}x/\nu\) being the local Grashof and Reynolds numbers, respectively. It is noticed that λ is a constant with \(\lambda<0\) and \(\lambda>0\) corresponding to the opposing and assisting flows, respectively, while \(\lambda=0 \) is for pure forced convection flow. Further, \(\varepsilon=c/a\) is the velocity ratio parameter, and \(\delta=(2-\sigma_{v})Kn_{x} \mathit{Re}_{x}^{1/2}/\sigma_{v}\) is the velocity slip parameter with the local Knudsen number \(Kn_{x}=\lambda_{0}/\sqrt{\varepsilon}x\).

In addition, the quantities of practical interest in this study are the skin friction coefficient \(C_{f}\) and the local Nusselt number \(\mathit{Nu}_{x}\), which are defined as

$$ C_{f}=\frac{\tau_{w}(x)}{\rho u_{e}^{2}}, \qquad \mathit{Nu}_{x}= \frac{xq_{w}(x)}{k(T_{w} -T_{\infty})}, $$
(14)

where the surface shear stress \(\tau_{w}(x)=\mu(\partial u/\partial y)_{y=0}\) and \(q_{w}(x)\) is the wall heat flux given by (8). Using the similarity variables (6), we obtain

$$ \mathit{Re}_{x}^{1/2}C_{f}=f''(0), \qquad \mathit{Re}_{x}^{-1/2}\mathit{Nu}_{x}=1/ \theta(0). $$
(15)

3 Series solutions of HAM

In the framework of the HAM technique, we select the initial guesses and the linear operator as

$$\begin{aligned}& f_{0}(\eta)=S+\eta+\frac{(\varepsilon-1)}{2\delta+1}\biggl(1-\frac{1}{1+\eta } \biggr),\qquad \theta_{0}(\eta)=\frac{1}{1+\eta}, \end{aligned}$$
(16)
$$\begin{aligned}& \mathcal {L}_{f}=f'''+ \frac{3}{1+\eta}f'', \qquad \mathcal {L}_{\theta}= \theta''+\frac {2}{1+\eta}\theta', \end{aligned}$$
(17)

where \(\mathcal {L}_{f}\) and \(\mathcal {L}_{\theta}\) satisfy

$$\begin{aligned}& \mathcal {L}_{f}\biggl[C_{1}+C_{2}\eta+ \frac{C_{3}}{1+\eta}\biggr]=0, \end{aligned}$$
(18)
$$\begin{aligned}& \mathcal {L}_{\theta}\biggl[C_{4}+\frac{C_{5}}{1+\eta} \biggr]=0, \end{aligned}$$
(19)

in which \(C_{i}\) (\(i=1\mbox{-}5\)) are arbitrary constants.

The zeroth order deformation equations are given by

$$\begin{aligned}& (1-q)\mathcal {L}_{f}\bigl[\bar{f}(\eta,q)-f_{0}(\eta) \bigr]=qh_{f}H_{f}(\eta)N_{f}\bigl[\bar {f}(\eta,q), \bar{\theta}[\eta,q]\bigr], \end{aligned}$$
(20)
$$\begin{aligned}& (1-q)\mathcal {L}_{\theta}\bigl[\bar{\theta}(\eta,q)-\theta_{0}( \eta)\bigr]=qh_{\theta}H_{\theta}(\eta) N_{\theta}\bigl[ \bar{f}(\eta,q),\bar{\theta}[\eta,q]\bigr], \end{aligned}$$
(21)

and they satisfy the following boundary conditions:

$$\begin{aligned}& \bar{f}(\eta,q)=S, \qquad \bar{f}'(\eta,q)=\varepsilon+\delta \bar{f}''(\eta ,q), \qquad \bar{\theta}'( \eta,q)=-1,\quad \mbox{at } \eta=0, \end{aligned}$$
(22)
$$\begin{aligned}& \bar{f}(\eta,q)\rightarrow1, \qquad \bar{\theta}(\eta,q)\rightarrow 0 \quad \mbox{as } \eta\rightarrow\infty. \end{aligned}$$
(23)

Here \(q\in(0,1)\) is an embedding parameter and \(h_{f}\) and \(h_{\theta}\) indicate the non-zero auxiliary parameters, \(H_{f}(\eta)\) and \(H_{\theta}(\eta)\) indicate the non-zero auxiliary functions, the nonlinear operators \(N_{f}\), \(N_{\theta}\) are defined as

$$\begin{aligned}& N_{f}\bigl[\bar{f}(\eta,q),\bar{\theta}(\eta,q)\bigr] = \frac{\partial^{3}\bar{f}(\eta ,q)}{\partial\eta^{3}}+\frac{m+1}{2}\bar{f}(\eta,q)\frac{\partial^{2}\bar {f}(\eta,q)}{\partial\eta^{2}} \\& \hphantom{N_{f}\bigl[\bar{f}(\eta,q),\bar{\theta}(\eta,q)\bigr] =}{} + m\biggl[1-\biggl(\frac{\partial\bar{f}(\eta,q)}{\partial\eta}\biggr)^{2} \biggr]+M\biggl(1-\frac {\partial\bar{f}(\eta,q)}{\partial\eta}\biggr) \\& \hphantom{N_{f}\bigl[\bar{f}(\eta,q),\bar{\theta}(\eta,q)\bigr] =}{} + \lambda\bar{\theta}(\eta,q), \end{aligned}$$
(24)
$$\begin{aligned}& N_{\theta}\bigl[\bar{f}(\eta,q),\bar{\theta}(\eta,q)\bigr] = \frac {\partial^{2}\bar{\theta}(\eta,q)}{\partial\eta^{2}}+ \frac{P_{r}(m+1)}{2}\bar{f}(\eta,q)\frac{\partial\bar{\theta}(\eta ,q)}{\partial\eta} \\& \hphantom{N_{\theta}\bigl[\bar{f}(\eta,q),\bar{\theta}(\eta,q)\bigr] =}{} - P_{r}(2m-1)\bar{\theta}(\eta,q) \frac{\partial\bar{f}(\eta,q)}{\partial \eta}. \end{aligned}$$
(25)

For \(q=0\) and \(q=1\), we have

$$\begin{aligned}& \bar{f}(\eta,0)=f_{0}(\eta), \qquad \bar{f}(\eta,1)=f( \eta), \end{aligned}$$
(26)
$$\begin{aligned}& \bar{\theta}(\eta,0)=\theta_{0}(\eta), \qquad \bar{\theta}(\eta,1)= \theta (\eta). \end{aligned}$$
(27)

By using a Taylor series, it is easy to obtain

$$\begin{aligned}& \bar{f}(\eta,q)=f_{0}(\eta)+\sum_{n=1}^{\infty}f_{n}(\eta)q^{n}, \end{aligned}$$
(28)
$$\begin{aligned}& \bar{\theta}(\eta,q)=\theta_{0}(\eta)+\sum _{n=1}^{\infty}\theta_{n}(\eta)q^{n}, \end{aligned}$$
(29)

in which

$$ f_{n}(\eta)={\biggl.\frac{1}{n!}\frac{\partial^{n}\bar{f}(\eta ,q)}{\partial q^{n}}\biggr|_{q=0}}, \qquad \theta_{n}(\eta)={\biggl.\frac{1}{n!}\frac{\partial^{n}\bar{\theta}(\eta ,q)}{\partial q^{n}}\biggr|_{q=0}}. $$
(30)

The parameters \(h_{f}\) and \(h_{\theta}\) are properly selected such that series solutions converge at \(p = 1\). Substituting \(p = 1\) into (28)-(29) gives

$$\begin{aligned}& f(\eta)=f_{0}(\eta)+\sum_{n=0}^{\infty}f_{n}(\eta), \end{aligned}$$
(31)
$$\begin{aligned}& \theta(\eta)=\theta_{0}(\eta)+\sum_{n=0}^{\infty}\theta_{n}(\eta). \end{aligned}$$
(32)

Then the mth order deformation equations and boundary conditions are

$$\begin{aligned}& \mathcal {L}_{f}\bigl[f_{n}(\eta)-\chi_{n} f_{n-1}(\eta)\bigr]=h_{f}H_{f}(\eta) \Re_{n}^{f}(\eta), \end{aligned}$$
(33)
$$\begin{aligned}& \mathcal {L}_{\theta}\bigl[\theta_{n}(\eta)-\chi_{n} \theta_{n-1}(\eta)\bigr]=h_{\theta}H_{\theta}(\eta) \Re_{n}^{\theta}(\eta), \end{aligned}$$
(34)
$$\begin{aligned}& f_{n}(0)=f_{n}'(0)-\delta f_{n}''(0)=f_{n}'( \infty)=0,\qquad \theta_{n}'(0)=\theta _{n}( \infty)=0, \end{aligned}$$
(35)
$$\begin{aligned}& \Re_{n}^{f}(\eta)=f_{n-1}'''-Mf_{n-1}'+ \lambda\theta_{n-1}+\frac {m+1}{2}\sum_{k=0}^{n-1}f_{n-1-k}f_{k}'' \\& \hphantom{\Re_{n}^{f}(\eta)=}{}-m\sum_{k=0}^{n-1}f_{n-1-k}'f_{k}'+(M+m) (1-\chi_{n}), \end{aligned}$$
(36)
$$\begin{aligned}& \Re_{n}^{\theta}(\eta)=\theta_{n-1}''+ \frac{P_{r}(m+1)}{2}\sum_{k=0}^{n-1}f_{n-1-k} \theta_{k}'-P_{r}(2m-1)\sum _{k=0}^{n-1}f_{n-1-k}' \theta_{k} \end{aligned}$$
(37)

with

$$ \chi_{n}=\left \{ \begin{array}{l@{\quad}l} 0, & n\leq1, \\ 1, & n\geq1. \end{array} \right . $$
(38)

The general solutions of (33)-(35) are

$$\begin{aligned}& f_{n}(\eta)=\tilde{f}_{n}(\eta)+C_{1}+C_{2} \eta+\frac{C_{3}}{1+\eta}, \end{aligned}$$
(39)
$$\begin{aligned}& \theta_{n}(\eta)=\tilde{\theta}_{n}(\eta)+C_{4}+ \frac{C_{5}}{1+\eta}, \end{aligned}$$
(40)

where \(\tilde{f}_{n}(\eta)\) and \(\tilde{\theta}_{n}(\eta)\) are the special solutions of the mth order deformation equations and \(C_{i}\) (\(i=1\mbox{-}5\)) are governed by the boundary conditions (35), which are given by

$$ C_{2}=C_{4}=0, \qquad C_{1}=-C_{3}- \tilde{f}_{n}(0), \qquad C_{3}=\frac{\tilde{f}_{n}'(0)-\delta \tilde{f}_{n}''(0)}{1+2\delta}, \qquad C_{5}=\tilde{\theta}_{n}'(0). $$
(41)

For simplicity, here we take \(\hbar_{f}=\hbar_{\theta}=\hbar\). In addition, according to the rule of the solution expression and the mth order deformation equation, the auxiliary functions \(H_{f}(\eta)\) and \(H_{\theta}(\eta)\) are chosen in the form

$$ H_{f}(\eta)=\frac{1}{(1+\eta)^{4}}, \qquad H_{\theta}( \eta)=\frac{1}{(1+\eta)^{3}}. $$
(42)

By employing the software MATHEMATICA, (33) and (34) can be solved one after the other in the order \(n=1,2,3,\ldots\) .

4 Analysis of the results

4.1 Convergence of the solutions

In order to ensure the convergence of the obtained series solutions, it is necessary to choose the appropriate range for the auxiliary parameter ħ. As pointed out by Liao [24], the interval for the admissible values for ħ corresponds to line segments nearly parallel to the horizontal axis. Figures 2 and 3 are plotted to show the admissible values of ħ for the function \(f''(0)\) and \(\theta(0)\) at 20th order approximation. To ensure the convergence of the series solution by HAM, it is observed from Figures 2 and 3 that the value of ħ should be chosen from \(-0.6\leq\hbar\leq-0.1\).

Figure 2
figure 2

ħ -Curves of \(\pmb{f''(0)}\) and \(\pmb{\theta(0)}\) for the 20th order approximation when \(\pmb{\varepsilon=0.1}\) .

Figure 3
figure 3

ħ -Curves of \(\pmb{f''(0)}\) and \(\pmb{\theta(0)}\) for the 20th order approximation when \(\pmb{\varepsilon=2}\) .

From the computation, it is found that the series solutions (31) and (32) converge in the whole region of η when we take \(\hbar=-0.5\). Table 1 shows the convergence of the solutions for different orders of approximations when \(\mathit{Pr}=0.7\), \(M=1\), \(m=2\), \(S=1\), \(\lambda=-1\), and \(\delta=1\). In order to further validate the present results, we show a comparison with previous work in Table 2. The results indicate that the numerical and analytical solutions are in good agreement. Hence we are confident that our results are accurate.

Table 1 Values of \(\pmb{f''(0)}\) and \(\pmb{1/\theta(0)}\) for different orders of approximations
Table 2 Values of \(\pmb{f''(0)}\) and \(\pmb{1/\theta(0)}\) for different orders of approximations

4.2 Results and discussion

The influence of key parameters on the coefficient of skin friction, the local Nusselt number, the velocity, and the temperature profiles are shown in this section. Figures 4 and 5 present the variation of the skin friction coefficient \(f''(0)\) and the local Nusselt number \(1/ \theta(0)\) with mixed convection parameter λ. The results show that the values of \(f''(0)\) and \(1/ \theta(0)\) increase with λ. It is also observed from Figure 4 that the values of \(f''(0)\) are negative for \(\varepsilon=2 >1\), which means that the sheet exerts a drag force on the fluid. For \(\varepsilon=0.1<1\), the values of \(f''(0)\) become positive, which indicates that the formation of the boundary layer does not depend solely on the stretching sheet. On the other hand, the values of \(1/ \theta(0)\), which represents the heat transfer rate at the surface, increase with the nonlinearity parameter m and they are always positive, i.e. the heat is transferred from the hot sheet to the cold fluid.

Figure 4
figure 4

Variation of the skin friction coefficient \(\pmb{f''(0)}\) with λ for different values of m .

Figure 5
figure 5

Variation of the local Nusselt number \(\pmb{1/ \theta (0)}\) with λ for different values of m .

The effects of mixed convection parameter λ on the velocity and temperature profiles are shown in Figures 6 and 7. The flow has a boundary layer structure, while the velocity ratio parameter \(\varepsilon<1\) and the nonlinearity parameter \(m\geq1\). On the other hand, when \(\varepsilon>1\) and \(m\geq1\), the flow has an inverted boundary layer structure which results from the fact that the stretching velocity \(u_{w}(x)\) of the surface exceeds the velocity \(u_{e}(x)\) of the external stream. From Figure 6, it is observed that the thickness of the velocity boundary layer decreases with λ for \(\varepsilon=0.1<1\), but it increases as λ increases for \(\varepsilon=2>1\). The velocity field \(f'(\eta)\) is always an increasing function of λ for the two cases. Figure 7 shows that an increase in λ corresponds to a decrease in the temperature and the thermal boundary layer thickness. From this figure, we can see that the wall temperature \(\theta(0)\) also decreases with increasing λ, which means that the heat transfer rate \(1/\theta(0)\) increases as λ increases. It is worth pointing out here that the temperature profiles show less deviation for different values of λ.

Figure 6
figure 6

Velocity profile \(\pmb{f'(\eta)}\) for different values of λ .

Figure 7
figure 7

Temperature profile \(\pmb{\theta(\eta)}\) for different values of λ .

Figures 8-11 illustrate the effects of m, S, δ, and M on the velocity profiles for the two cases of \(\varepsilon<1\) and \(\varepsilon>1\). It is observed that the thickness of the velocity boundary layer decreases with increasing values of m, S, δ, and M for both cases, which implies an increasing magnitude of the velocity gradient at the surface. Thus, the skin friction coefficient \(f''(0)\) increases with the increasing values of m, S, δ, and M. From a physical point of view this follows from the fact that with a rise in the strength of magnetic parameter M, the Lorentz force associated with the magnetic field makes the boundary layer thinner. Further, it is seen from Figures 9-11 that the velocity inside the velocity boundary layer increases with S, δ, and M for \(\varepsilon<1\) but decreases with these parameters for \(\varepsilon>1\). In addition, it is indicated that the velocity slip parameter has a significant influence on the velocity compared to the parameters m, S, and M.

Figure 8
figure 8

Velocity profile \(\pmb{f'(\eta)}\) for different values of m .

Figure 9
figure 9

Velocity profile \(\pmb{f'(\eta)}\) for different values of S .

Figure 10
figure 10

Velocity profile \(\pmb{f'(\eta)}\) for different values of δ .

Figure 11
figure 11

Velocity profile \(\pmb{f'(\eta)}\) for different values of M .

Figures 12-15 show the temperature profiles for selected values of the parameters m, S, δ, and M for the two cases of \(\varepsilon<1\) and \(\varepsilon>1\). The temperature is found to decrease to zero monotonically as η increases, which satisfies the far field boundary condition \(\theta (\infty)=0\). The results display that the temperature and thickness of the thermal boundary layer is lower for \(\varepsilon=2\) and higher for \(\varepsilon =0.1\) when the other parameters are constant. Figures 12 and 13 show that the wall temperature \(\theta(0)\) decreases with increasing m, S for the two cases. Thus, the heat transfer rate at the surface increases as m and S increase. Furthermore, it is easy to see that the temperature and thickness of thermal boundary layer also decrease with increasing these two parameters. Different characteristics are observed in Figures 14 and 15. The wall temperature and thickness of thermal boundary layer decrease with increasing δ and M for \(\varepsilon<1\), and the reverse trend is observed for the case of \(\varepsilon>1\). Thus, the heat transfer rate at the surface increases with δ or M for \(\varepsilon<1\), and the opposite behaviors are observed for the effects of δ and M for \(\varepsilon >1\).

Figure 12
figure 12

Temperature profile \(\pmb{\theta(\eta)}\) for different values of m .

Figure 13
figure 13

Temperature profile \(\pmb{\theta(\eta)}\) for different values of S .

Figure 14
figure 14

Temperature profile \(\pmb{\theta(\eta)}\) for different values of δ .

Figure 15
figure 15

Temperature profile \(\pmb{\theta(\eta)}\) for different values of M .

5 Conclusions

In this work, the MHD mixed convection stagnation-point flow toward a nonlinearly stretching vertical sheet is analyzed. Different from the previous works, the current results focus on the effect of nonlinearly vertical stretching for the MHD stagnation-point flow with mixed convection. The analytic solutions for momentum and energy equations have been obtained by the method of HAM. The main conclusions can be summarized as follows.

  • The increase of nonlinearity parameter m leads to an increases of the heat transfer rate at surface \(1/\theta(0)\), and to a decrease of both the thicknesses of the velocity and the thermal boundary layer.

  • The coefficient of the skin friction \(f''(0)\) and the heat transfer rate at the surface \(1/\theta(0)\) increase with increasing mixed convection parameter λ.

  • The heat transfer rate increases as the velocity slip parameter δ and magnetic parameter M increase for \(\varepsilon < 1\), but it decreases with these two parameters for \(\varepsilon> 1\).

  • Inside the velocity boundary layer, the velocity increases with the increasing S, δ, and M for \(\varepsilon<1\), and the opposite trend is displayed for \(\varepsilon>1\).

  • Inside the thermal boundary layer, the temperature always decreases with increasing λ, m, and S.

References

  1. Ishak, A, Jafar, K, Nazar, R, Pop, I: MHD stagnation point flow towards a stretching sheet. Physica A 388(13), 3377-3383 (2009)

    Article  Google Scholar 

  2. Khan, ZH, Khan, WA, Qasim, M, Shah, IA: MHD stagnation point ferrofluid flow and heat transfer toward a stretching sheet. IEEE Trans. Nanotechnol. 13(1), 35-40 (2014)

    Article  Google Scholar 

  3. Shateyi, S, Makinde, OD: Hydromagnetic stagnation-point flow towards a radially stretching convectively heated disk. Math. Probl. Eng. 2013, Article ID 616947 (2013)

    MathSciNet  MATH  Google Scholar 

  4. Aman, F, Ishak, A, Pop, I: Magnetohydrodynamic stagnation-point flow towards a stretching/shrinking sheet with slip effects. Int. Commun. Heat Mass Transf. 47, 68-72 (2013)

    Article  Google Scholar 

  5. Singh, G, Makinde, OD: MHD slip flow of viscous fluid over an isothermal reactive stretching sheet. Ann. Fac. Eng. Hunedoara Int. J. Eng. 11(2), 41-46 (2013)

    Google Scholar 

  6. Shateyi, S, Marewo, GT: Numerical analysis of unsteady MHD flow near a stagnation point of a two-dimensional porous body with heat and mass transfer, thermal transfer, and chemical reaction. Bound. Value Probl. 2014, 218 (2014)

    Article  MATH  Google Scholar 

  7. Ishak, A, Nazar, R, Pop, I: MHD mixed convection boundary layer flow towards a stretching vertical surface with constant wall temperature. Int. J. Heat Mass Transf. 53(23-24), 5330-5334 (2010)

    Article  MATH  Google Scholar 

  8. Kumari, M, Nath, G: Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid. Int. J. Non-Linear Mech. 45(3), 310-319 (2010)

    Article  Google Scholar 

  9. Ali, FM, Nazar, R, Arifin, NM, Pop, I: Effect of Hall current on MHD mixed convection boundary layer flow over a stretched vertical flat plate. Meccanica 46(5), 1103-1112 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  10. Hayat, T, Qasim, M: Radiation and magnetic field effects on the unsteady mixed convection flow of a second grade fluid over a vertical stretching sheet. Int. J. Numer. Methods Fluids 66(7), 820-832 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  11. Turkyilmazoglu, M: The analytical solution of mixed convection heat transfer and fluid flow of a MHD viscoelastic fluid over a permeable stretching surface. Int. J. Mech. Sci. 77, 263-268 (2013)

    Article  Google Scholar 

  12. Chamkha, AJ, El-Kabeir, S: Unsteady heat and mass transfer by MHD mixed convection flow over an impulsively stretched vertical surface with chemical reaction and Soret and Dufour effects. Chem. Eng. Commun. 200(9), 1220-1236 (2013)

    Article  Google Scholar 

  13. Hassanien, IA, El-Hawary, HM, Mahmoud, MAA: Thermal radiation effect on flow and heat transfer of unsteady MHD micropolar fluid over vertical heated nonisothermal stretching surface using group analysis. Appl. Math. Mech. 34(6), 703-720 (2013)

    Article  MathSciNet  Google Scholar 

  14. Shateyi, S: A new numerical approach to MHD flow of a Maxwell fluid past a vertical stretching sheet in the presence of thermophoresis and chemical reaction. Bound. Value Probl. 2013, 196 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  15. Makinde, OD: Heat and mass transfer by MHD mixed convection stagnation point flow toward a vertical plate embedded in a highly porous medium with radiation and internal heat generation. Meccanica 47, 1173-1184 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  16. Ali, FM, Nazar, R, Arifin, NM, Pop, I: Mixed convection stagnation-point flow on vertical stretching sheet with external magnetic field. Appl. Math. Mech. 35(2), 155-166 (2014)

    Article  MathSciNet  Google Scholar 

  17. Makinde, OD, Khan, AH, Khan, ZH: Buoyancy effects on MHD stagnation point flow and heat transfer of a nanofluid past a convectively heated stretching/shrinking sheet. Int. J. Heat Mass Transf. 62, 526-533 (2013)

    Article  Google Scholar 

  18. Das, S, Jana, RN, Makinda, OD: MHD boundary layer slip flow and heat transfer of nanofluid past a vertical stretching sheet with non-uniform heat generation/absorption. Int. J. Nanosci. 13(3), 1450019 (2014)

    Article  Google Scholar 

  19. Akyildiz, FT, Siginer, DA, Vajravelu, K, Cannon, JR, Van Gorder, RA: Similarity solutions of the boundary layer equations for a nonlinearly stretching sheet. Math. Methods Appl. Sci. 33(5), 601-606 (2010)

    MATH  MathSciNet  Google Scholar 

  20. Akyildiz, FT, Siginer, DA: Galerkin-Legendre spectral method for the velocity and thermal boundary layers over a non-linearly stretching sheet. Nonlinear Anal., Real World Appl. 11(2), 735-741 (2010)

    Article  MATH  MathSciNet  Google Scholar 

  21. Dhanai, R, Rana, P, Kumar, L: 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. Powder Technol. 273, 62-70 (2015)

    Article  Google Scholar 

  22. Ashraf, MB, Hayat, T, Alsaedi, A: Three-dimensional flow of Eyring-Powell nanofluid by convectively heated exponentially stretching sheet. Eur. Phys. J. Plus 130(1), 1-16 (2015)

    Article  MathSciNet  Google Scholar 

  23. Mabood, F, Khan, WA, Ismail, AIM: MHD boundary layer flow and heat transfer of nanofluids over a nonlinear stretching sheet: a numerical study. J. Magn. Magn. Mater. 374, 569-576 (2015)

    Article  Google Scholar 

  24. Liao, S: Beyond Perturbation: Introduction to the Homotopy Analysis Method. CRC Press, Boca Raton (2003)

    Book  Google Scholar 

  25. Rashidi, MM, Pour, SAM: Analytic approximate solutions for unsteady boundary-layer flow and heat transfer due to a stretching sheet by homotopy analysis method. Nonlinear Anal., Model. Control 15(1), 83-95 (2010)

    MATH  MathSciNet  Google Scholar 

  26. Nadeem, S, Mehmood, R, Akbar, NS: Non-orthogonal stagnation point flow of a nano non-Newtonian fluid towards a stretching surface with heat transfer. Int. J. Heat Mass Transf. 57(2), 679-689 (2013)

    Article  Google Scholar 

  27. Hayat, T, Iqbal, Z, Mustafa, M, Alsaedi, A: Stagnation-point flow of Jeffrey fluid with melting heat transfer and Soret and Dufour effects. Int. J. Numer. Methods Heat Fluid Flow 24(2), 402-418 (2014)

    Article  MathSciNet  Google Scholar 

  28. Hayat, T, Farooq, M, Alsaedi, A: Melting heat transfer in the stagnation-point flow of Maxwell fluid with double-diffusive convection. Int. J. Numer. Methods Heat Fluid Flow 24(3), 760-774 (2014)

    Article  MathSciNet  Google Scholar 

  29. Hayat, T, Iqbal, Z, Mustafa, M, Alsaedi, A: Unsteady flow and heat transfer of Jeffrey fluid over a stretching sheet. Therm. Sci. 18(4), 1069-1078 (2014)

    Article  Google Scholar 

  30. Mabood, F, Khan, WA: Approximate analytic solutions for influence of heat transfer on MHD stagnation point flow in porous medium. Comput. Fluids 100, 72-78 (2014)

    Article  MathSciNet  Google Scholar 

  31. Yacob, NA, Ishak, A: MHD flow of a micropolar fluid towards a vertical permeable plate with prescribed surface heat flux. Chem. Eng. Res. Des. 89(11), 2291-2297. (2011)

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the referees for their pertinent comments and valuable suggestions. This work is supported by the National Natural Science Foundation of China (Grant No. 51305080).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ming Shen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MS and HC conceived of the study and formulated the problem. FW developed the Mathematica codes and generated the results. All authors participated in the analysis of the results and manuscript coordination. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, M., Wang, F. & Chen, H. MHD mixed convection slip flow near a stagnation-point on a nonlinearly vertical stretching sheet. Bound Value Probl 2015, 78 (2015). https://doi.org/10.1186/s13661-015-0340-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0340-6

Keywords