Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 56, 1, pp. 81-94, Warsaw 2018 DOI: 10.15632/jtam-pl.56.1.81 A LOCALLY MODIFIED SINGLE-PHASE MODEL FOR ANALYZING MAGNETOHYDRODYNAMIC BOUNDARY LAYER FLOW AND HEAT TRANSFER OF NANOFLUIDS OVER A NONLINEARLY STRETCHING SHEET WITH CHEMICAL REACTION Pooria Akbarzadeh School of Mechanical Engineering, Shahrood University of Technology, Shahrood, Iran e-mail: akbarzad@ut.ac.ir; p.akbarzadeh@shahroodut.ac.ir Theproblemof boundary layer flowandheat transfer of nanofluids overnonlinear stretching of a flat sheet in the presence of amagnetic field and chemical reaction is investigated nume- rically. In this paper, anew locallymodified single-phasemodel for the analysis is introduced. In thismodel, the effective viscosity, density and thermal conductivity of the solid-liquidmi- xtures (nanofluids) which are commonly utilized in the homogenous single-phasemodel, are locally combinedwith the prevalent single-phasemodel. Similarity transformation is used to convert the governing equations into three coupled nonlinear ordinarydifferential equations. These equations depend onfive local functions of the nanoparticle volume fractionviz., local viscosity ratio,magnetic, Prandtl,Brownianmotion and thermophoresis functions.The equ- ations are solved using Newton’s method and a block tridiagonalmatrix solver. The results are compared to the prevalent single-phase model. In addition, the effect of important go- verning parameters on the velocity, temperature, volume fraction distribution and the heat andmass transfer rates are examined. Keywords: nanofluid,magnetohydrodynamic, nonlinear stretching sheet, similarity transfor- mation, locally modified single-phasemodel 1. Introduction The enhancement of heat transfer in an industrial process such as cooling of electronic devices, guide or thrust bearings, high-speed sliding surfaces, strips or filaments, and etc., may create energy savings, reduce the process time, increase efficiency and lengthen the working life of equipment. Convective heat transfer can be improved by changing flow geometry, boundary conditions, or by enhancing thermal conductivity of the fluid (Wang and Mujumdar, 2007). One way to increase the thermal conductivity of base fluids is the use of nanoparticles (such as Cu, Ag, Al, and etc. as metallic solid, CuO, Al2O3, TiO2 and etc., as metallic oxide particles, multi walled nanotubes (MWNTs) and single wall nanotubes (SWNTs) (asmetallic nanotubes) suspended in the fluids. Indeed, the Brownian motion of the nanoparticles in these suspensions is one of the potential contributors to this enhancement (Choi, 1995; Choi et al., 2001; Xie et al., 2003; Li and Peterson, 2007; Aminossadati and Ghasemi, 2009; Hamad and Ferddows, 2012). Apparently, Choi (1995) was the first researcher who introduced the term nanofluids to refer to the fluid with suspended nanoparticles, and up from this time, many researchers focused on numerical or experimental studies in the field of thermophysical and also heat and fluid flow properties of nanofluids such by Choi (1995), Choi et al. (2001), Xie et al. (2003), Li andPeterson (2007), Aminossadati andGhasemi (2009), Hamad andFerddows (2012),Wen and Ding Y. (2004), Heris et al. (2006), Abu-Nada andChamkhac (2010), Sheikhzadeh et al. (2011) and to name but a few. 82 P. Akbarzadeh As aforesaid, cooling of continuous strips or filaments by drawing them through amotionless fluid is one of the heat transfer processes in the industry. In these respects, the feature of the final product (e.g. plastic and polymer sheets manufactured by extrusion, glass-fiber and paper production,metal spinning, etc.) dependson the cooling rate and theprocess of stretching (Rana and Bhargava, 2012). Therefore, simulation of boundary layer behavior of fluid flow over the stretching surface (strips or filaments) can be useful for predicting heat transfer characteristics of such a process. The basic and elementary works on Boundary Layer Flow (BLF) over a continuousmoving surface (with constant velocity) were reported by Sakiadis (1961a,b,c). After that, Crane (1970) published a report about the exact temperature distribution of the steady BLF of a viscous fluid caused by stretching a flat isothermal sheet (with a linearly variable velocity). The temperature and concentration distributions of BLF over an isothermal moving plate with blowing or suction were found by Gupta and Gupta (1977). The same problem was studied byChen andChar (1988) in which the sheet was subjected to a prescribed temperature and heat flux. Theboundary layer and thermal characteristics of fluidflowover a nonlinear stretching sheet were examined for the first time by Vajravelu (2001). Then, Vajravelu and Cannon (2006) in- vestigated the existence of a solution for the nonlinear problem by using the Schauder theory. Afterwards, Cortell (2007) analyzed flow and heat transfer over a nonlinear stretching sheet for the prescribed and constant surface temperature as boundary conditions. Prasad et al. (2010) presented a numerical solution forMagneto-Hydrodynamics (MHD) flow of an electrically con- ductingviscousfluidover a stretching sheetwithvariablefluidproperties.Theyassumed that the stretching velocity and the transverse magnetic field varied as a power function of the distance from the origin. Postelnicu andPop (2011) explored steady 2D laminarBLFof a non-Newtonian power-law fluid past a permeable non-linear stretching wedge. And recently, Vajravelu et al. (2014) published a paper on the subject of MHD flow and heat transfer of a non-Newtonian power-law fluid over an unsteady stretching isothermal sheet. Since using nanofluids may improve heat transfer behavior of an engineering process, many researchers examined the effect of those fluids on the boundary layer and thermal characteristics of flow caused by stretching a surface. In this respect, Bachok et al. (2010) studied the steady BLF of nanofluids over a continuous moving surface (with constant velocity). Khan and Pop (2010) obtained the temperaturedistributionof steadyBLFofnanofluidsover a linear isothermal stretching flat sheet. Gorder et al. (2010) presented the similarity solution for the nano-BLFs over a linearly stretching sheet, in which on the surface sheet the velocity slip was assumed to be proportional to the local shear stress. Hassani et al. (2011) obtained an analytical solution for BLF of a nanofluid past a linearly stretching sheet using the homotopy analysis method (HAM). Rana and Bhargava (2012) studied numerically flow and heat transfer of a nanofluid over anonlinearly stretching sheet.Theyused thefinite difference, finite element, andvariational method in their computations. Nadeem et al. (2014) analyzed flow of a three-dimensional water- -based nanofluid over an exponentially stretching sheet. And recently, Das (2015) investigated the problem of BLF of a nanofluid over a non-linear permeable stretching sheet at a predestined surface temperature in the presence of partial slip. The study of Magneto-Hydrodynamics (MHD) flow of an electrically conducting fluid due to a stretching sheet is important in modern engineering processes. Metallurgy, metalworking, metal fusion in an electrical furnace, and etc., are some examples of such processes (Ibrahim and Shanker, 2014). The MHD flow analysis over a stretching sheet with various aspects such as types of fluids, magnetic effect, stagnation geometry and the temperature effect can be seen in works of many researchers such as Ishak et al. (2008), Fadzilah et al. (2011), Mahapatra et al. (2009), Prasad et al. (2010), and to name but a few. Ibrahim et al. (2013) studied the MHD stagnation-point flow and heat transfer due to a nanofluid towards a stretching sheet using the Runge-Kutta fourth order numerical method. They analyzed the effect of the velocity A locally modified single-phase model for analyzing magnetohydrodynamic... 83 ratio parameter on both the local Nusselt number and local Sherwood number. In the same year, Ibrahim and Shankar (2013) investigated MHD boundary layer flow and heat transfer of a nanofluid past a permeable stretching sheet with velocity, thermal and solutal slip boundary conditions. Recently, Ibrahim and Shankar (2014) studied MHD boundary layer flow and heat transfer of a nanofluid over a non-isothermal stretching sheet. Their analysis was done for two different cases, namely a prescribed surface temperature and prescribed heat flux. The above mentioned literature review shows that in order to investigate the heat transfer characteristics of nanoparticles (or generally small solid particles) suspended in a fluid, twomain approaches have been adopted by researchers. The first approach is the two-phase model which considers both the fluid phase and the solid particles behavior in the heat transfer process. The second one is the single-phase model in which both phases are in thermal and hydrodynamic equilibrium state (this approach is simpler and more computationally efficient). Generally in nanofluids, there are several factors that affect heat transfer enhancement. Some of the mo- re important factors are Brownian motion (including diffusion, sedimentation, and dispersion), gravity, layering at the solid/liquid interface, particle clustering, friction between the fluid and the solid particles, etc. Therefore, in the absence of any experimental data and suitable theore- tical studies, the existing macroscopic two-phase model has not enough precision for analyzing nanofluids. Consequently, the modified single-phase, considering some of the above factors, is more convenient than the two-phase model if the main interest of analysis is the heat transfer process (Khanafer et al., 2003). Therefore, in order to improve the results of the single-phase model for analyzing the nanofluids flow, some modifications are needed. In this paper, a new modified single-phase model for analyzing flow and heat transfer in nanofluids is introduced for the first time. In this model, all effective properties of nanofluids such as density, viscosity and thermal conductivity, which are normally used for the effective single-phase model (as constant values), are incorporated locally with the governing equations (as non-constant values). This approach is used for examining the BLF behavior and thermal characteristics of nanofluids flow over a nonlinear stretching sheet in the presence of a magnetic field and chemical reaction. The results for Cu and Al2O3 nanoparticles are compared to the prevalent single-phase model. This comparison depicts that the prevalent single-phase model has a considerable deviation for predicting the behavior of nanofluids flow, especially in dimensionless temperature and nano- particle volume fraction. In addition, the effects of important governing parameters such as the transverse magnetic field, chemical reaction strength, thermophoresis parameter, nanoparticle volume fraction near the surface, etc., on the velocity, temperature, volume fraction distribution and dimensionless heat andmass transfer rates are examined. 2. Mathematical governing equations The steady-state two-dimensional BLF of a nanofluid past a stretching sheet is considered with the nonlinear velocity uw = aux n, where au is a positive constant, n is the nonlinear stretching parameter and x is the coordinate measured along the stretching surface, as shown in Fig. 1. Because of the impermeability characteristic of the sheet, the vertical velocity of the fluid on the surface is vw = 0. In this problem, the stretching surface is a non-isothermal face with the relation of Tw(x) = T∞ +aTx r, where aT is a positive constant, r is the surface temperature parameter in the prescribed surface temperature boundary condition andT∞ is the temperature of the fluid far away from the stretching sheet. In addition, the fluid is under a transverse magnetic field with strength B(x) which is applied in the vertical direction, given by the special form of B(x)= B0x n−1, where B0 is a positive constant. Also, it is assumed that the first-order homogeneous chemical reaction with the non-linear rate of K = K0x n−1 occurs in the fluid (K0 is the constant chemical reaction parameter). According to (Ibrahim and Shanker, 2014; 84 P. Akbarzadeh Fig. 1. Physical model of the stretching sheet and the coordinate system Buongiorno, 2006; Yazdi et al., 2011), the following four governing equations of BLF include the continuity, momentum, energy, and nanoparticles concentration, respectively ∂u ∂x + ∂v ∂y =0 u ∂u ∂x +v ∂u ∂y = νnf ∂2u ∂y2 − σB ρnf u u ∂T ∂x +v ∂T ∂y = αnf ∂2T ∂y2 + (ρc)p (ρc)nf [ DB ∂ϕ ∂y ∂T ∂y + DT T∞ (∂T ∂y )2] u ∂ϕ ∂x +v ∂ϕ ∂y = DB ∂2ϕ ∂y2 + DT T∞ ∂2T ∂y2 −Kϕ (2.1) here, T is temperature of the fluid, u and v are velocity components along the x and y-axis respectively, ϕ is the nanoparticle concentration (or volume fraction), DT is the thermophoretic diffusion coefficient, DB is the Brownian diffusion coefficient, c is specific heat capacity, ρ is density, α = k/(ρc) is thermal diffusivity, ν is kinematic viscosity, k is thermal conductivity, and σ is electrical conductivity of the fluid. The subscripts p and nf refer to the nanoparticles and nanofluid, respectively. In this study, a new local Modified Single-Phase Model (here it is named asMSPM) for analyzing nanofluids flow and heat transfer is introduced. In thismanner, the parameters νnf, knf, ρnf, and cnf of the above governing equations may be introduced by the following four relations (Oztop and Abu-Nada, 2008) (ρc)nf (ρc)f =1−ϕ+ϕ(ρc)p (ρc)f ρnf ρf =1−ϕ+ϕρp ρf νnf νf = 1 (1−ϕ)2.5(1−ϕ+ϕρp/ρf) knf kf = kp+2kf −2ϕ(kf −kp) kp+2kf +ϕ(kf −kp) (2.2) where, the subscript f refers to the base fluid. The associated boundary conditions for the problem are u(x,0)= uw(x)= aux n v(x,0)= vw(x)= 0 T(x,0)= Tw(x)= aTx r+T∞ ϕ(x,0)= ϕw u(x,∞) = 0 v(x,∞) = 0 T(x,∞)= T∞ ϕ(x,∞) = ϕ∞ (2.3) whereϕw andϕ∞ are thenanoparticle volume fraction neighborhood to the surface and far away from the stretching sheet, respectively. By considering the following similarity transformations (Ibrahim and Shankar, 2013) η = yβx n−1 2 u = aux nf ′(η) v =−γx n−1 2 ( f + n−1 n+1 ηf ′ ) θ(η)= T −T∞ Tw−T∞ ψ(η) = ϕ−ϕ∞ ϕw −ϕ∞ (2.4) A locally modified single-phase model for analyzing magnetohydrodynamic... 85 where, β = √ au(n+1)/(2νf) and γ = √ auνf(n+1)/2, BLF governing equations (2.1) are transformed into three non-linear ordinary differential equations as follows f ′′′+F(ϕ) ( ff ′′− 2n n+1 f ′2 ) −G(ϕ) 2 n+1 f ′ =0 1 Pr(ϕ) θ′′+fθ′− 2r n+1 f ′θ+Nb(ϕ)ψ′θ′+Nt(ϕ)θ ′2 =0 ψ′′+Lefψ′+ Nt Nb θ′′− 2R n+1 Leψ =0 (2.5) where F(ϕ) is the local viscosity ratio function, G(ϕ), Pr(ϕ), Nb(ϕ) and Nt(ϕ) are the local magnetic, local Prandtl, local Brownian and local thermophoresis functions are defined as F(ϕ) = (1−ϕ)2.5 ( 1−ϕ+ϕρp ρf ) G(ϕ) = M(1−ϕ)2.5 Nb(ϕ) = Nb 1−ϕ+vp(ρc)p/(ρc)f Nt(ϕ)= Nt 1−ϕ+ϕ(ρc)p/(ρc)f Pr(ϕ)= νf αnf =Prf kp/kf +2+2ϕ(1−kp/kf) kp/kf +2−ϕ(1−kp/kf) ( 1−ϕ+ϕ(ρc)p (ρc)f ) (2.6) In the set of equations (2.6), M = σB0/(auρf) is the magnetic parameter, Prf = νf/αf is the fluid Prandtl number, Le = νf/DB is the Lewis number, Nb = [(ρc)pDB(ϕw −ϕ∞)]/[νf(ρc)f] is the Brownian motion number, Nt = [(ρc)pDT(Tw −T∞)]/[T∞νf(ρc)f ] is the thermophoresis parameter and, finally, R = K0/au. It is clear that Nb is related to Le as follows Nb= (ρc)p(ϕw −ϕ∞) Le(ρc)f (2.7) It should be noted that the enhancement of nanofluids thermal conductivity is because of four major mechanisms: (1) Brownian motion of nanoparticles, (2) nanolayer, (3) clustering, and (4) nature of heat transport in the nanoparticles. Also, some important parameters which affect the thermal conductivity of nanofluids are particle volume fraction, temperature, particles size and the size and property of the nanolayer. Therefore, a considerable number of studies can be found on the modelling of thermal conductivity of nanofluids, such as the Maxwell model, Hamilton-Crossermodel,Bruggemanmodel,Waspmodel,Maxwell-Garnett (MG)model,Gupta model, and to name but a few. However, there is a lack of reliable and comprehensive model which includes all mechanisms and influenced parameters for thermal conductivity of nanofluids (Esfe et al., 2014; Kumar et al., 2015). As amatter of fact, the present study does not focus on a comprehensive thermal conductivity model of nanofluids (so the Maxwell model (Eqs. (2.5) is selected as the model). The novelty of this study can be found in Eqs. (2.5) and (2.6). In Eqs. (2.5), F(ϕ), Pr(ϕ), Nb(ϕ), and Nt(ϕ) are not constant during the calculation and they are updated locally based on the third relation of Eq. (2.5), which is more physical. It means that three relations of Eqs. (2.5) should be solved together. By consideringEqs. (2.3), non-linear ordinary differential equations (2.5) are solved subject to the following boundary conditions f ′(0)= 1.0 θ(0)= 1.0 ψ(0)= 1.0 f(0)= 0 f ′(∞)= 0 θ(∞)=0 ψ(∞)= 0 (2.8) By introducing the local Reynolds number via Rex = uw(x)x/νf = aux n+1/νf, the surface heat flux through qw = knf(∂T/∂y)y=0, the surface mass flux via qm = DB(∂ϕ/∂y)y=0, the local nanofluid Nusselt number through Nunf = xqw/[knf(Tw −T∞)] and the local nanofluid 86 P. Akbarzadeh Sherwood number by means of Shnf = xqm/[DB(ϕw − ϕ∞)], the following relation can be established Nunf√ Rex = √ n+1 2 |θ′(0)| Shnf√ Rex = √ n+1 2 |ψ′(0)| (2.9) where Nunf/ √ Rex and Shnf/ √ Rex are known as dimensionless heat and mass transfer rates, respectively. Without loss of generality, one can assume that the nanoparticle volume fraction far away from the stretching sheet is zero or ϕ∞ =0. In this paper, the fluid is awater based na- nofluid containing different types of prevalent nanoparticles: copper (Cu) and alumina (Al2O3). Based on (Oztop and Abu-Nada, 2008), the thermophysical properties of the fluid (water) and the mentioned nanoparticles are given in Table 1. Table 1.Thermophysical properties of the fluid (water) and three nanoparticles Properties Fluid phase (water) Cu Al2O3 c [J/kgK] 4179 385 765 ρ [kg/m3] 997.1 8933 3970 k [W/mK] 0.613 400 40 Prf 6.83 – – (ρc)p/(ρc)f – 0.825 0.729 kp/kf – 652.53 65.25 3. Algorithm of numerical solution The algorithm of numerical calculation (computer programming procedure) of non-linear dif- ferential equations (2.5)-(2.7) with boundary conditions (2.8) is given via items (i) to (viii). It should be noted that for step (ii), by assuming ϕ∞ = 0, the initial guess can be ψ(η) = ϕ(η)/ϕw = [1− exp(η − ηmax)]/[1− exp(−ηmax)] as an example, where ηmax is the maximum value of η in the numerical calculation. (i) Specify input data of the base fluid and nanoparticles: ρp/ρf, (ρc)p/(ρc)f , kp/kf, Prf, Nb, Nt, Le, ϕw (ii) Set an initial guess for ϕ(η) or its dimensionless parameter, i.e. ψ(η) (iii) Calculate F(ϕ), G(ϕ), Nb(ϕ), Nt(ϕ), and Pr(ϕ) based on step (i) and profile of ϕ (iv) Solve the 1st relation of equation (2.5) to obtain f(η) by using Newton’s method (v) Solve the 2nd relation of equation (2.5) to obtain θ(η) by using Newton’s method (vi) Solve the 3rd relation of equation (2.5) to update ψ(η) or ϕ(η) by usingNewton’s method (vii) Check the error between the updated ψ(η) and the guessed one (viii) If the maximum error¬ 10−5 then the procedure is finished; else go to step (iii) 4. Validation of the numerical solution With reference tomany studies such asKhan andPop (2010), Das (2015), Ibrahimand Shankar (2013), Ibrahim and Shankar (2014) and etc., it is evident that theMSPM is reduced to the old Prevalent Single-Phase Model (PSPM) by considering the relations: F(ϕ) = 1.0, G(ϕ) = M, Nb(ϕ)=Nb, Nt(ϕ)=Nt, and Pr(ϕ)=Pr. In order to check the validity of the present compu- tational programming code, governing equations (2.5)-(2.7) subject to boundaryconditions (2.8) are solved numerically for some values of the governing parameters of the PSPM. The results A locally modified single-phase model for analyzing magnetohydrodynamic... 87 for the reduced Nusselt number |θ′(0)| and the reduced Sherwood number |ψ′(0)| are compared with those obtained by Khan and Pop (2010) in Table 2. In this simulation, the default values of the parameters are considered as Pr = 10, Le = 10, n = 1, r = 0, M = 0, R = 0 and ηmax = 20. Also, Table 3 shows the comparison of the magnitude of the velocity gradient at the wall, |f ′′(0)|, between the present code results and that obtained previously by Fang et al. (2009) for the case of n =1, ηmax =20, and M =0.25-4.0. It can be seen from Tables 2 and 3 that the present results are in very good agreement with those reported by other researchers. Therefore, it is clear that the results obtained in this study are accurate. Table 2. Comparison of |θ′(0)| and |ψ′(0)| for Pr = 10, Le= 10, n =1, r =0, R =0, M =0, ηmax =20 Parameter Present result Khan and Pop (2010) |θ′(0)| |ψ′(0)| |θ′(0)| |ψ′(0)| Nb=0.1, Nt=0.1 0.95238 2.12939 0.9524 2.1294 Nb=0.1, Nt=0.3 0.52007 2.52863 0.5201 2.5286 Nb=0.1, Nt=0.5 0.32105 3.03514 0.3211 3.0351 Nb=0.3, Nt=0.1 0.25215 2.41002 0.2522 2.4100 Nb=0.3, Nt=0.3 0.13551 2.60882 0.1355 2.6088 Nb=0.3, Nt=0.5 0.08329 2.75187 0.0833 2.7519 Nb=0.5, Nt=0.1 0.05425 2.38357 0.0.543 2.3836 Nb=0.5, Nt=0.3 0.02913 2.49837 0.0291 2.4984 Nb=0.5, Nt=0.5 0.01792 2.57310 0.0179 2.5731 Table 3.Comparison of results for |f ′′(0)| when n =1, and ηmax =20 Parameter Present result Fang et al. (2009) |f ′′(0)| |f ′′(0)| M =0.25 1.11803 1.1180 M =4.0 2.23606 2.2361 5. Results and discussion In this Section, the numerical results for profiles of dimensionless velocity f ′(η), temperatu- re θ(η), nanoparticle concentration ψ(η), and etc., are presented for different values of the governing parameters. The obtained results are displayed through graphs in Figs. 2-7. For all simulations and their corresponding figures, the used governing parameters are given inTable 4. It should be noted that the value of the Brownian motion number Nb for the MSPM is cal- culated by equation (2.7), while for PSPM this is considered as the average of Nb for Cu and Al2O3 nanoparticles. In order to compare the effect of usingMSPM and PSPM on the results, the linear and non-linear stretching sheet problems for two different nanoparticles (i.e. Cu and Al2O3) are considered.These comparisons for sixdimensionless profiles (i.e. f,f ′, θ,θ′,ψ andψ′) are presented in Fig. 2 (linear stretching sheet problem) and Fig. 3 (non-linear stretching sheet problem). The illustrated results in both Figs. 2 and 3 confirm that the temperature and na- noparticle volume fraction profiles (i.e. θ and ψ) converge quicker than the horizontal velocity profiles (i.e. f ′) for both PSPMandMSPM.Nevertheless, it is evident from the figures that the PSPM has a remarkable deviation for predicting the behavior of nanofluids flow, especially in dimensionless temperature θ and nanoparticle concentration ψ. Also, the comparison between the results presented in Figs. 2 and 3 displays that this deviation is intensifiedwhen the sheet is 88 P. Akbarzadeh stretched nonlinearly. It should be noted that discussions about the PSPMand its effect on the results of the stretching sheet problemwere recently reported inmany studies such asKhan and Pop (2010), Hassani et al. (2011), Ibrahim and Shankar (2013), Ibrahim and Shanker (2014), and to name but a few. Therefore, for the next simulations, theMSPM is only considered. Also, based on the results shown in Figs. 2 and 3, one can observe that the boundary layer beha- vior and thermal characteristics of both nanoparticle flows are approximately the same. Hence, the further results and discussions are focused only on the Cu-water nanofluid flow (using the MSPM). Fig. 2. Effect of theMSPM and PSPM on the dimensionless profiles for n =1 and r =0 Table 4.The governing parameters for nanofluid flow simulations Nb Fig. n r Nt Le ϕw M R MSPM MSPM PSPM (Cu) (Al2O3) 2 1 0 0.1 5 0.5 0.25 0.1 0.0825 0.0729 0.0777 3 5 5 0.1 5 0.5 0.25 0.1 0.0825 0.0729 0.0777 4 5 5 0.1 5 0.5 0.0-4.0 0.1 0.0825 – – 5 5 5 0.1 5 0.5 2.0 0.0-10 0.0825 – – 6 5 5 0.1 1 0.1-0.5 2.0 5 0.0825-0.4125 – – 7a 5 5 0.1 5 0.5 0.0-4.0 0.1 0.0825 – – 7b 5 5 0.1 5 0.5 2.0 0.0-10 0.0825 – – 7c 5 5 0.1 1 0.1-0.5 2.0 5 0.0825-0.4125 – – A locally modified single-phase model for analyzing magnetohydrodynamic... 89 Fig. 3. Effect of theMSPM and PSPM on the dimensionless profiles for n =5 and r =5 Fig. 4. Effect of M on the dimensionless profiles for n =5 and r =5 The first analysis is related to the effects of the transverse magnetic field M on the flow and thermal characteristics. By applying a transversemagnetic field, a Lorentz force is created, which results in a retarding force on the velocity of the flow (Ibrahim and Shanker, 2014). Therefore, as M increases (and, consequently, increasing the retarding force), the velocity of the fluid decreases. This fact can be observed in Fig. 4a. In addition, Fig. 4b illustrates the impact of the transverse magnetic field on the temperature profile. The results show that as the magnetic parameter M increases, the temperature profile and the thermal boundary layer thickness increase. Also, fromFig. 4c it is seen that the behavior of dimensionless concentration is the same as the temperature profile when the values of the magnetic parameter M increase. 90 P. Akbarzadeh Chemical reaction is an important process that should be considered in micro-mixing of biological systems such as cell-activation and protein-folding, particularly when the mixing of reactants for initiation is necessary. Lorentz forces produced by a magnetic field are able to move nano-liquids in a mixing process. Hence, one of the active methods of micro-mixing of biological samples is using a magnetic field in the presence of chemical reaction (Yazdi et al., 2011). Therefore, the second examination belongs to the effects of the chemical reaction strength R on the nanofluids BLF characteristics. Figure 5 demonstrates the influence of the chemical reaction strength on the velocity, temperature and concentration graph. FromFig. 5a, it is evident that the values of R has not any significant effect on the dimensionless veloci- ty profile. But, by increasing the values of the chemical reaction parameter, the temperatu- re gradient (Fig. 5b) and concentration gradient (Fig. 5c) on the wall intensify. In fact, for the case of the higher nanoparticle concentration gradient near the wall, the nanoparticle con- centration decreases rapidly as η increases. Therefore, the conductive heat transfer process is weakened. Fig. 5. Effect of R on the dimensionless profiles for n =5 and r =5 As theMSPMdepends on the value of the nanoparticle volume fraction near the surface ϕw, it is of great worth to investigate the effects of ϕw on dynamic and thermal characteristics of the nanofluidsflow.As it is predicted,when thevalue of themagnetic parameterM does not change, variation of ϕw has no significant influence on the velocity graph (see Fig. 6a). However, as it is noticed from Fig. 6b, increasing the nanoparticle volume fraction near the surface increases the thermal boundary layer thickness. And finally for Fig. 6c, since ψ = ϕ/ϕw, comparing the concentration graphs for various ϕw has no any important consequence. Fig. 6. Effect of ϕ w on the dimensionless profiles for n=5 and r =5 A locally modified single-phase model for analyzing magnetohydrodynamic... 91 The last results are about the effect of M, R, and ϕw on the dimensionless heatNunf/ √ Rex andmass transfer rates Shnf/ √ Rex on thewall.This information ispresented inFig. 7.Figure 7a illustrates that the heat andmass transfer rates on the stretching wall decrease with an increase in the transverse magnetic field M. In fact, the magnetic field generates more heat in the boundary layer region and, hence, this reduces the wall heat transfer rate. Figure 7b depicts that the wall mass transfer rate and the wall heat transfer rate increase with an increase in the chemical reaction parameter R. In addition, a decrease in the dimensionless heat transfer and an increase in the dimensionless mass transfer on the wall are observed with an increase in ϕw. These are shown in Fig. 7c. Fig. 7. Effect of M, R, and ϕ w on the dimensionless heat andmass transfer rates for n =5 and r =5 6. Conclusion The problem of boundary layer flow and heat transfer of nanofluids over a nonlinear stretching sheet in the presence of a magnetic field and chemical reaction is examined numerically. In this study, a modified single-phase model for analyzing nanofluids flow and heat transfer is initiated. In thismodifiedmodel, the effective density andviscosity ofnanofluidsand the effective thermal conductivity of the solid-liquidmixturewhichareprevalently used in the effective single- phase model (as constant values) are incorporated locally with the governing equations (as no- constant values). A similarity solution is proposed which depends on the local Prandtl number, local Brownian motion number, local Lewis number and the local thermophoresis number. The results for Cu andAl2O3 nanoparticles are compared to the prevalent single-phase model. This comparison shows that theprevalent single-phasemodel has anoticeable deviation for predicting the behavior of the nanofluids flow, especially in dimensionless temperature and nanoparticle volume fraction. In addition, the results exhibit that the heat and mass transfer rates on the stretching surface decrease with increase in the transverse magnetic field. Also, the wall mass transfer rate and the wall heat transfer rate increase with an increase in the chemical reaction parameter. Acknowledgments The author would like to acknowledge the Shahrood University of Technology which supported this project. References 1. Abu-Nada E., Chamkhac A.J., 2010, Mixed convection flow in a lid-driven inclined square enclosure filled with a nanofluid,European Journal of Mechanics B/Fluids, 29, 472-482 92 P. Akbarzadeh 2. Aminossadati S.M., Ghasemi B., 2009, Natural convection cooling of a localised heat source at the bottom of a nanofluid-filled enclosure,European Journal of Mechanics B/Fluids, 28, 630-640 3. Bachok N., Ishak A., Pop, I., 2010, Boundary-layer flow of nanofluids over a moving surface in a flowing fluid, International Journal of Thermal Sciences, 49, 1663-1668 4. Buongiorno J., 2006, Convective transport in nanofluids,ASME Journal of Heat Transfer, 128, 240-250 5. Chen C.K., Char M. I., 1988, Heat transfer of a continuous stretching surface with suction or blowing, Journal Mathematical Analysis and Applications, 135, 568-580 6. Choi S.U.S., 1995, Enhancing thermal conductivity of fluids with nanoparticles, ASME Fluids Engineering Division, 231, 99-105 7. Choi S.U.S., Zhang Z.G., YuW., LockwoodF.E., GrulkeE.A., 2001,Anomalous thermal conductivity enhancement in nanotube suspensions,Applied Physics Letters, 79, 14, 2252-2254 8. Cortell R., 2007, Viscous flow and heat transfer over a nonlinearly stretching sheet, Applied Mathematics and Computation, 184, 864-873 9. CraneL.J., 1970,Flowpasta stretchingplate,KurzeMitteilungen-BriefReports-Communications Breves, 21, 645-647 10. DasK., 2015,Nanofluidflowovera non-linearpermeable stretching sheetwithpartial slip,Journal of the Egyptian Mathematical Society, 23, 2, 451-456 11. Esfe M.H., Saedodin S., Mahian O., Wongwises S., 2014, Thermal conductivity of Al2O3/water nanofluids – Measurement, correlation, sensitivity analysis, and comparisons with literature reports, Journal of Thermal Analysis and Calorimetry, 117, 2, 675-681 12. FadzilahM.,NazarR.,NorihanM.,Pop I., 2011,MHDboundary-layerflowandheat transfer over a stretching sheetwith inducedmagnetic field,Journal of Heat andMassTransfer,47, 155-162 13. Fang T., Zhang J., Yao S., 2009, Slip MHD viscous flow over a stretching sheet – An exact solution,Communications in Nonlinear Science and Numerical Simulation, 14, 3731-3737 14. GorderR.A.V.,SweetE.,VajraveluK., 2010,Nanoboundary layersover stretching surfaces, Communications in Nonlinear Science and Numerical Simulation, 15, 1494-1500 15. Gupta P.S., Gupta A.S., 1977, Heat and mass transfer on a stretching sheet with suction or blowing,The Canadian Journal of Chemical Engineering, 55, 744-746 16. Hamad M.A.A., Ferddows M., 2012, Similarity solutions to viscous flow and heat transfer of nanofluid over nonlinearly stretching sheet,AppliedMathematics andMechanics (English Edition), 33, 7, 923-930 17. Hassani M., Mohammad Tabar Nemati H., Domairry G., Noori F., 2011, An analytical solution for boundary layer flow of a nanofluid past a stretching sheet, International Journal of Thermal Sciences, 50, 2256-2263 18. Heris S.Z.,EtemadS.G.,EsfahanyM.N., 2006,Experimental investigationof oxidenanofluids laminar flow convective heat transfer, International Communications in Heat and Mass Transfer, 33, 529-535 19. Ibrahim W., ShankarB., 2013,MHDboundary layer flow and heat transfer of a nanofluid past a permeable stretching sheetwith velocity, thermal and solutal slip boundary conditions,Computer and Fluids, 75, 1-10 20. Ibrahim W., Shankar B., Nandeppanavar M.M., 2013,MHD stagnation point flow and heat transfer due to nanofluid towards a stretching sheet, International Journal of Heat and Mass Transfer, 56, 1-9 21. IbrahimW., ShankerB., 2014,Magnetohydrodynamic boundary layer flow andheat transfer of a nanofluid over non-isothermal stretching sheet,ASME Journal of Heat Transfer, 136, 051701-9 A locally modified single-phase model for analyzing magnetohydrodynamic... 93 22. IshakA.,NazarR.,Pop, I., 2008,Hydromagneticflowandheat transferadjacent toa stretching vertical sheet, Journal of Heat and Mass Transfer, 44, 921-927 23. KhanW.A.,Pop I., 2010,Boundary-layerflowofananofluidpasta stretching sheet, International Journal of Heat and Mass Transfer, 53, 2477-2483 24. Khanafer K., Vafai K., LightstoneM., 2003, Buoyancy-drivenheat transfer enhancement in a two-dimensional enclosure utilizing nanofluids, International Journal of Heat andMass Transfer, 46, 3639-3653 25. Kumar P.M., Kumar J., Tamilarasan R., Sendhilnathan S., Suresh S., 2015, Review on nanofluids theoretical thermal conductivity models,Engineering Journal, 19, 1, 67-83 26. LiC.H.,PetersonG.P., 2007,Mixing effect on the enhancement of the effective thermal conduc- tivity of nanoparticle suspensions (nanofluids), International Journal of Heat and Mass Transfer, 50, 4668-4677 27. MahapatraT.R.,NandyS.K.,GuptaA.S., 2009,Magnetohydrodynamicstagnationpointflow of a power-law fluid towards a stretching sheet, International Journal of Non-Linear Mechanics, 44, 124-129 28. Nadeem S., Haq R.U., Khan Z.H., 2014, Heat transfer analysis of water-based nanofluid over an exponentially stretching sheet,Alexandria Engineering Journal, 53, 219-224 29. Oztop H.F., Abu-Nada E., 2008, Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids, International Journal of Heat and Fluid Flow, 29, 1326-1336 30. Postelnicu A., Pop I., 2011, Falkner-Skan boundary layer flow of a power-law fluid past a stretching wedge,Applied Mathematics and Computation, 217, 4359-4368 31. PrasadK.V.,VajraveluK.,DattiP.S., 2010,Mixedconvectionheat transfer overanon-linear stretching surface with variable fluid properties, International Journal of Non-Linear Mechanics, 45, 320-330 32. RanaP., BhargavaR., 2012,Flow andheat transfer of a nanofluid over a nonlinearly stretching sheet: A numerical study, Communications in Nonlinear Science and Numerical Simulation, 17, 212-226 33. Sakiadis B.C., 1961a, Boundary layer behavior on continuous moving solid surfaces. I. Bounda- ry layer equations for two-dimensional and axis-symmetric flow, American Institute of Chemical Engineer Journal, 7, 1, 26-28 34. SakiadisB.C., 1961b,Boundary layer behavior on continuousmoving solid surfaces. II. Boundary layer on a continuous flat surface,American Institute of Chemical Engineer Journal, 7, 2, 221-225 35. SakiadisB.C., 1961c,Boundary layerbehavior on continuousmoving solid surfaces. III.Boundary layer on a continuous cylindrical surface, American Institute of Chemical Engineer Journal, 7, 3, 467-472 36. Sheikhzadeh G.A., Arefmanesh A., Kheirkhah M.H., Abdollahi R., 2011, Natural co- nvection of Cu-water nanofluid in a cavity with partially active side walls, European Journal of Mechanics B/Fluids, 30, 2, 166-176 37. Vajravelu K., 2001, Viscous flow over a nonlinearly stretching sheet, Applied Mathematics and Computation, 124, 281-288 38. Vajravelu K., Cannon J.R., 2006, Fluid flow over a nonlinearly stretching sheet, Applied Ma- thematics and Computation, 181, 609-618 39. Vajravelu K., Prasad K.V., Datti P.S., Raju B.T., 2014, MHD flow and heat transfer of an Ostwald-de Waele fluid over an unsteady stretching surface, Ain Shams Engineering Journal, 5, 157-167 40. Wang X.Q., Mujumdar A.S., 2007, Heat transfer characteristics of nanofluids: a review, Inter- national Journal of Thermal Sciences, 46, 1-19 94 P. Akbarzadeh 41. Wen D., Ding Y., 2004, Experimental investigation into convective heat transfer of nanofluids at the entrance region under laminar flow conditions, International Journal of Heat and Mass Transfer, 47, 5181-5188 42. Xie H.Q., Lee H., Youn W., Choi M., 2003, Nanofluids containing multiwalled carbon nano- tubes and their enhanced thermal conductivities, Journal of Applied Physics, 94, 8, 4967-4971 43. YazdiM.H.,Abdullah S.,Hashim I., SopianK., 2011, SlipMHD liquid flowandheat transfer over non-linear permeable stretching surfacewith chemical reaction, International Journal of Heat and Mass Transfer, 54, 3214-3225 Manuscript received May 15, 2017; accepted for print July 25, 2017