We study the effect of thermophoresis on boundary layer magneto-nanofluid flow over a stretching sheet. The model includes the effects of Brownian motion and cross-diffusion effects. The governing partial differential equations are transformed to a system of ordinary differential equations and solved numerically using a spectral linearisation method. The effects of the magnetic influence number, the Prandtl number, Lewis number, the Brownian motion parameter, thermophoresis parameter, the modified Dufour parameter and the Dufour-solutal Lewis number on the fluid properties as well as on the heat, regular and nano mass transfer coefficients are determined and shown graphically.
Most common fluids such as water, ethylene, glycol, toluene or oil generally have poor heat transfer characteristics owing to their low thermal conductivity. A recent technique to improve the thermal conductivity of these fluids is to suspend nano-sized metallic particles such as aluminum, titanium, gold, copper, iron or their oxides in the fluid to enhance its thermal properties, Choi . The enhancement of thermal conductivity in nanofluids has been studied by, among others, Kakac and Pramuanjaroenkij , Choi et al., Masuda et al., Eapen et al. and Fan and Wang . Nield and Kuznetsov  analyzed the behaviour of boundary layer flow on the Chen-Minkowycz problem in a porous layer saturated with a nanofluid. Nield and Kuznetsov  investigated thermal instability in a porous medium saturated with nanofluid using the Brinkman model. The model incorporated the effects of Brownian motion and thermophoresis of nanoparticles. They found that the critical thermal Rayleigh number can be reduced or increased by a substantial amount depending on whether the nanoparticle distribution is top-heavy or bottom-heavy. Aziz et al. studied steady boundary layer flow past a horizontal flat plate embedded in a porous medium filled with a water-based nanofluid containing gyrotactic microorganisms. Cheng  investigated the behaviour of boundary layer flow over a horizontal cylinder of elliptic cross section in a porous medium saturated with a nanofluid. Chamkha et al. investigated the non-similar solutions for natural convective boundary layer flow over a sphere embedded in a porous medium saturated with a nanofluid.
During the last few decades, fluid flow over a stretching surface has received considerable attention because of its engineering applications such as in melt-spinning, hot rolling, wire drawing, glass-fiber production and the manufacture of polymer and rubber sheets, Altan and Gegel , Fisher , and Tidmore and Klein . Nanofluid flow over a stretching surface has been investigated by many researchers. The first study on a stretching sheet in nanofluids was published by Khan and Pop . Makinde and Aziz  performed a numerical study of boundary layer flow over a linear stretching sheet. Both Brownian motion and thermophoresis effects on the transport equations were presented. They reported that stronger Brownian motion and thermophoresis lead to an increase in the rate of heat transfer. However, the opposite was observed in the case of the rate of mass transfer. Recent studies in this area include those of Narayana and Sibanda  and Kameswaran et al..
Magnetic nanofluids have numerous uses or potential applications in engineering and medicine. Using magnetic nanofluids has the potential to regulate the flow rate and heat transfer by controlling the thermo-magnetic convection current and the fluid velocity (see Shima et al., Ganguly et al.). The effects of a magnetic field on nanofluid flow over a stretching sheet have been investigated by, among others, Bachok et al. and Hanad and Ferdows .
The aim of this study is to analyse Dufour and Soret effects in a magneto-nanofluid flow over a stretching sheet. In addition, we study Brownian motion and thermophoresis effects using a spectral linearisation method to obtain numerical solutions of the momentum, energy, concentration and mass fraction equations. The successive linearisation method (SLM) is an accurate method for solving non-linear coupled equations (see [23-25]). Recent studies such as [26-28] have suggested that the SLM is accurate and converges rapidly to the numerical results when compared to other semi-analytical methods such as the Adomian decomposition method, the variational iteration method and the homotopy perturbation method.
2 Mathematical formulation
Consider two-dimensional nanofluid flow over a linearly stretching sheet with velocity , where a is a real positive number. The coordinate system is assumed to define the x-axis along the surface of the sheet and y is the coordinate normal to the surface of the sheet. The surface temperature and nanoparticle concentration are higher than the ambient values and , respectively. The governing equations for the problem can be written in the form
with the boundary conditions
where u and v are the velocity components along the x and y direction respectively, σ is the electrical conductivity, is magnetic field flux density, ν kinematic viscosity of the base fluid, α is the thermal diffusivity of the porous medium, is the Brownian diffusion coefficient, is thermophoresis diffusion coefficient, and are the Soret and Dufour diffusivities, is the solutal diffusivity, T is the fluid temperature, C is the solutal concentration, is the nanoparticle volume fraction, and are the heat capacity of the fluid and the effective heat capacity of the nanoparticle material respectively, τ is a parameter defined by . Using the similarity variables
equations (1)-(5) reduce to the following non-similar forms where primes denote differentiation with respect to η:
subject to the boundary conditions
The parameters in equations (8)-(11) are the magnetic number M, the Prandtl number Pr, the Lewis number Le, the Brownian motion parameter Nb, the thermophoresis parameter Nt, the nanofluid Lewis number Ln, the modified Dufour parameter Nd and the Dufour-solutal Lewis number Ld. These parameters are defined as
The parameters of engineering interest in heat and mass transport problems are the local Nusselt number , the Sherwood number and the nanofluid Sherwood number . These parameters characterise the wall heat, the regular and nano mass transfer rates, respectively, and are defined by
Following Khan and Aziz , the physical parameters of interest are the reduced Nusselt Nur, the Sherwood number and the reduced Sherwood Shr defined as
3 Method of solution
where , , and are unknown and , , and ( ) are successive approximations that are obtained by recursively solving the linear forms of the equation system that results from substituting (13) into equations (8)-(11). In particular, the linearised equations to be solved are
subject to the boundary conditions
where coefficient parameters , , , ( ) and ( ) are known constants. The initial guesses , , and are chosen to satisfy the boundary conditions
and are chosen as
Starting from the initial guesses and iterating times, the functions , and are written as
where is the order of the SLM approximation. Equations (14)-(17) are solved using the Chebyshev spectral collocation method. The method is based on the Chebyshev polynomials defined on the interval . We first transform the domain of solution into the domain using the domain truncation technique where the problem is solved in the interval where L is a scaling parameter used to invoke the boundary condition at infinity. This is achieved by using the mapping
We discretise the domain using the Gauss-Lobatto collocation points given by
where N is the number of collocation points used. The functions , , and for are approximated at the collocation points as follows:
where is the kth Chebyshev polynomial given by
The derivatives of the variables evaluated at the collocation points are represented as
Substituting equations (22)-(26) into equations (14)-(17) leads to the matrix equation
where is a square matrix and and are column vectors defined by
The functions and parameters in equation (29) are
In the definitions above, T stands for transpose, ( ), ( ), ( ), ( ) and ( ) are diagonal matrices of order , I is an identity matrix of order and is a zero matrix of order . The solution is obtained as
4 Results and discussion
In this section we present solutions of equations (8)-(11) along with the boundary conditions (12) using the SLM iteration scheme. Tables 1 and 2 give a comparison between the present results and Khan and Pop  for the reduced Nusselt and Sherwood numbers respectively. There is a good agreement between the two sets of results with the SLM having converged at the fourth order up to eleven decimal places. The velocity components and are plotted in Figures 1(a) and 1(b) for different values of the magnetic field parameter M. As is now well known, the velocity decreases with increases in the magnetic field parameter due to an increase in the Lorentz drag force that opposes the fluid motion.
Figure 1 . Effect of the magnetic fieldMon the velocity components (a) and (b) .
Table 1 . Comparison of results for the reduced Nusselt number with , ,
Table 2 . Comparison of results for the reduced Sherwood number with , ,
Figures 2(a) and 2(b) show the effect of the thermophoresis parameter on the temperature and mass volume fraction profiles. The thermophoretic force generated by the temperature gradient creates a fast flow away from the stretching surface. In this way more fluid is heated away from the surface, and consequently, as Nt increases, the temperature within the boundary layer increases. The fast flow from the stretching sheet carries with it nanoparticles leading to an increase in the mass volume fraction boundary layer thickness.
Figure 2 . Effect of the thermophoresis parameterNton the temperatureθand nanoparticleϕprofiles.
Figures 3(a) and 3(b) show the effect of the Lewis number Le, and the Dufour-solutal Lewis number Ld on the species concentration in the boundary layer. The concentration profiles significantly contract as the Lewis number increases. The effect of the random motion of the nanoparticles suspended in the fluid on the temperature and nanoparticle volume fraction is shown in Figures 4(a) and 4(b). As expected, the increased Brownian motion of the nanoparticles carries with it heat and the thickness of the thermal boundary layer increases. The Brownian motion of the nanoparticles increases thermal transport which is an important mechanism for the enhancement of thermal conductivity of nanofluids. However, we note that increasing the Brownian motion parameter leads to a clustering of the nanoparticles near the stretching sheet. An increase in the Brownian motion of the nanoparticles leads to a decrease in the mass volume fraction profiles.
Figure 3 . Effect of the Lewis numberLeand the Dufour-solutal Lewis numberLdon concentration profiles.
Figure 4 . Effect of the Brownian motion parameterNbon the temperature and nanoparticle volume fraction profiles.
Figures 5(a) and 5(b) show the temperature profiles for several values of the Prandtl number Pr and mass volume fraction profile for several values of the modified Dufour number Nd. The temperature profiles decrease as the Prandtl number increases since, for high Prandtl numbers, the flow is governed by momentum and viscous diffusion rather than thermal diffusion. On the other hand, the thickness of the mass volume fraction boundary layer increases with an increase in Nd.
Figure 5 . Effect ofPrandNdon the temperatureθ, profiles respectively.
Figures 6(a) and 6(b) show the effects of the thermophoresis parameter Nt, the Lewis number Le, the magnetic field parameter M, the Prandtl number Pr and the modified Dufour number Nd on the wall heat and mass fraction transfer rates. It can be seen that the thermal boundary layer thickness increases when the thermophoresis parameter Nt increases, thus decreasing the reduced Nusselt number. However, increasing the Lewis number Le leads to a decrease in the reduced Nusselt number. On the other hand, the results show that the reduced Nusselt number increases with increasing Prandtl numbers. Increasing both the magnetic field parameter M and the modified Dufour parameter Nd leads to an increase in the thermal boundary layer thickness, thus reducing the Nusselt number.
Figure 6 . Effect ofNt,Le,Pr,NdandMon the heat transfer coefficientNur.
Figures 7(a) and 7(b) show the effects of the Dufour-solutal Lewis number Ld and the nanofluid Lewis number Ln on the reduced Nusselt number Nur as the Brownian motion parameter Nb increases. We note a decrease in the reduced Nusselt number when Ln increases, and an increase in the reduced Nusselt number when Ld increases.
Figure 7 . Effect of the Dufour-solutal Lewis numberLdand the nanofluid Lewis numberLnon the reduced Nusselt numberNur.
Figures 8(a) and 8(b) show the graphs of and plotted against the Dufour-solutal Lewis number Ld for different values of the parameters Nt, Nb and Le. We observe that increases in the absence of the Brownian motion and the thermophoresis parameter while decreases in the presence of Brownian motion and thermophoresis parameters. An increase in is observed in the presence of both the Brownian motion and the thermophoresis parameter. Figures 9(a) and 9(b) show the effect of increasing Nt and Nb respectively on the reduced Sherwood number .
Figure 8 . Effect of the Lewis numberLe, the thermophoresis parameterNtand the Brownian motion parameterNbon (a) the reduced Nusselt number and (b) the local Sherwood number .
Figure 9 . Effect of the nanofluid Lewis numberLn, the thermophoresis parameterNtand the Brownian motion parameterNbon the nanofluid Sherwood number .
A numerical study of the magneto-nanofluid boundary layer flow over a stretching sheet was carried out. We determined the effects of various parameters on the fluid properties as well as on the heat, and the regular and nano mass transfer rates. We have shown that increasing the magnetic field parameter M tends to retard the fluid flow within the boundary layer. The effects of the Prandtl number, the Lewis number, the Brownian motion parameter, the thermophoresis parameter, the nanofluid Lewis number, the modified Dufour parameter and the Dufour-solutal Lewis number on the heat, regular and nano mass transfer coefficients and fluid flow characteristics have been studied. We have shown inter alia that:
– the thermal boundary layer thickness increases with the thermophoresis parameter;
– increasing the Lewis number reduces the heat transfer coefficient;
– the heat transfer coefficient increases in the absence of the Brownian motion and the thermophoresis parameter and decreases in the presence of Brownian motion and thermophoresis parameters.
The authors declare that they have no competing interests.
The work including proofreading was done by all the authors.
The authors wish to thank the University of KwaZulu-Natal for financial support.
Kakac, S, Pramuanjaroenkij, A: Review of convective heat transfer enhancement with nanofluid . Int. J. Heat Mass Transf.. 52, 3187–3196 (2009). Publisher Full Text
Choi, SUS, Zhang, ZG, Yu, W, Lockwood, FE, Grulke, EA: Anomalously thermal conductivity enhancement in nanotube suspensions . Appl. Phys. Lett.. 79, 2252–2254 (2001). Publisher Full Text
Masuda, H, Ebata, A, Teramae, K, Hishinuma, N: Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles . Netsu Bussei. 7, 227–233 (1993). Publisher Full Text
Nield, DA, Kuznetsov, AV: Thermal instability in a porous medium layer saturated by a nanofluid: Brinkman model . Transp. Porous Media. 81, 409–422 (2010). Publisher Full Text
Aziz, A, Khan, WA, Pop, I: Free convection boundary layer flow past a horizontal flat plate embedded in porous medium filled by nanofluid containing gyrotactic microorganisms . Int. J. Therm. Sci.. 56, 48–57 (2012)
Cheng, C-Y: Free convection boundary layer flow over a horizontal cylinder of elliptic cross section in porus media saturated by a nanofluid . Int. Commun. Heat Mass Transf.. 39, 931–936 (2012). Publisher Full Text
Chamkha, A, Gorla, RSR, Ghodeswar, K: Non-similar solution for natural convective boundary layer flow over a sphere embedded in a porous medium saturated with a nanofluid . Transp. Porous Media. 86, 13–22 (2011). Publisher Full Text
Khan, WA, Pop, I: Boundary layer flow of a nanofluid past a stretching sheet . Int. J. Heat Mass Transf.. 53, 2477–2483 (2010). Publisher Full Text
Makinde, OD, Aziz, A: Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition . Int. J. Therm. Sci.. 50, 1326–1332 (2011). Publisher Full Text
Narayana, M, Sibanda, P: Laminar flow of a nanoliquid film over an unsteady stretching sheet . Int. J. Heat Mass Transf.. 55, 7552–7560 (2012). Publisher Full Text
Kameswaran, PK, Narayana, N, Sibanda, P, Murthy, PVSN: Hydromagnetic nanofluid flow due to a stretching or shrinking sheet with viscous dissipation and chemical reaction effects . Int. J. Heat Mass Transf.. 55, 7587–7595 (2012). Publisher Full Text
Bachok, N, Ishak, A, Pop, I: Unsteady boundary-layer flow and heat transfer of a nanofluid over a permeable stretching/shrinking sheet . Int. J. Heat Mass Transf.. 55, 2102–2109 (2012). Publisher Full Text
Awad, FG, Sibanda, P, Motsa, SS, Makinde, OD: Convection from an inverted cone in a porous medium with cross-diffusion effects . Comput. Math. Appl.. 61, 1431–1441 (2011). Publisher Full Text
Makukula, ZG, Sibanda, P, Motsa, SS: A novel numerical technique for two-dimensional laminar flow between two moving porous walls . Math. Probl. Eng.. 2010, (2010) Article ID 528956. doi:10.1155/2010/528956
Makukula, ZG, Motsa, SS, Sibanda, P: A novel numerical technique for two-dimensional laminar flow between two moving porous walls . Math. Probl. Eng.. 2010, (2010) Article ID 528956. doi:10.1155/2010/528956
Awad, FG, Sibanda, P, Narayana, M, Motsa, SS: Convection from a semi-finite plate in a fluid saturated porous medium with cross-diffusion and radiative heat transfer . Int. J. Phys. Sci.. 6, 4910–4923 (2011)
Motsa, SS, Sibanda, P, Shateyi, S: On a new quasi-linearization method for systems of nonlinear boundary value problems . Math. Methods Appl. Sci.. 34, 1406–1413 (2011). Publisher Full Text
Khan, WA, Aziz, A: Double-diffusive natural convective boundary layer flow in a porous medium saturated with a nanofluid over a vertical plate: prescribed surface heat, solute and nanoparticle fluxes . Int. J. Therm. Sci.. 50, 2154–2160 (2011). Publisher Full Text
Shateyi, S, Motsa, SS: Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with hall effect . Bound. Value Probl.. 2010, (2010) Article ID 257568. doi:10.1155/2010/257568
Don, WS, Solomonoff, A: Accuracy and speed in computing the Chebyshev collocation derivative . SIAM J. Sci. Comput.. 16, 1253–1268 (1995). Publisher Full Text