- Research
- Open access
- Published:
On unsteady MHD mixed convection in a nanofluid due to a stretching/shrinking surface with suction/injection using the spectral relaxation method
Boundary Value Problems volume 2015, Article number: 24 (2015)
Abstract
In this study we investigate heat and mass transfer in magnetohydrodynamic mixed convection flow of a nanofluid over an unsteady stretching/shrinking sheet. The flow is subject to a heat source, viscous dissipation and Soret and Dufour effects are assumed to be significant. We have further assumed that the nanoparticle volume fraction at the wall may be actively controlled. The physical problem is modeled using systems of nonlinear differential equations which we have solved numerically using the recent spectral relaxation method. In addition to the discussion on physical heat and mass transfer processes, we also show that the spectral relaxation technique is an accurate technique for solving nonlinear boundary value problems.
1 Introduction
Nanofluids are suspensions of metallic, non-metallic or polymeric nano-sized powders in a base liquid which are used to increase the heat transfer rate in various applications. In recent years, the concept of nanofluid has been proposed as a route for increasing the performance of heat transfer liquids. Due to the increasing importance of nanofluids, there is a large amount of literature on convective heat transport in nanofluids and problems linked to a stretching surface. An excellent collection of articles on this topic can be found in [1–4]. The majority of the previous studies have been restricted to boundary layer flow and heat transfer in nanofluids. Following the early work by Crane [5], Khan and Pop [6] were among the first researchers to study nanofluid flow due to a stretching sheet. Other researchers studied various aspects of flow and heat transfer in a fluid of infinite extent; see, for instance, Chen [7] and Abo-Eldahab and Abd El-Aziz [8]. A mathematical analysis of momentum and heat transfer characteristics of the boundary layer flow of an incompressible and electrically conducting viscoelastic fluid over a linear stretching sheet was carried out by Abd El-Aziz [9]. In addition, radiation effects on viscous flow of a nanofluid and heat transfer over a nonlinearly stretching sheet were studied by Hady et al. [10]. Theoretical studies include, for example, modeling unsteady boundary layer flow of a nanofluid over a permeable stretching/shrinking sheet by Bachok et al. [11]. Rohni et al. [12] developed a numerical solution for the unsteady flow over a continuously shrinking surface with wall mass suction using the nanofluid model proposed by Buongiorno [13].
The effect of an applied magnetic field on nanofluids has substantial applications in chemistry, physics and engineering. These include cooling of continuous filaments, in the process of drawing, annealing and thinning of copper wire. Drawing such strips through an electrically conducting fluid subject to a magnetic field can control the rate of cooling and stretching, thereby furthering the desired characteristics of the final product. Such an application of a linearly stretching sheet of incompressible viscous flow of MHD was discussed by Pavlov [14]. In other work, Jafar et al. [15] studied the effects of magnetohydrodynamic (MHD) flow and heat transfer due to a stretching/shrinking sheet with an external magnetic field, viscous dissipation and Joule effects.
A model for magnetohydrodynamic flow over a uniformly stretched vertical permeable surface subject to a chemical reaction was suggested by Chamkha [16]. An analysis of the effects of a chemical reaction on heat and mass transfer on a magnetohydrodynamic boundary layer flow over a wedge with ohmic heating and viscous dissipation in a porous medium was done by Kandasamy and Palanimani [17]. Rashidi and Erfani [18] studied the steady MHD convective and slip flow due to a rotating disk with viscous dissipation and ohmic heating. Rashidi et al. [19] found approximate analytic solutions for an MHD boundary-layer viscoelastic fluid flow over a continuously moving stretching surface using the homotopy analysis method. Rashidi and Keimanesh [20] used the differential transform method and Padé approximants to solve the equations that model MHD flow in a laminar liquid film from a horizontal stretching surface. The effect of a transverse magnetic field on the flow and heat transfer over a stretching surface were examined by Anjali-Devi and Thiyagarajan [21]. The influence of a chemical reaction on heat and mass transfer due to natural convection from vertical surfaces in porous media subject to Soret and Dufour effects was also studied by Postelnicu [22].
Despite all this previous work, there is still a lot that is unknown about the flow and heat and mass transfer properties of different nanofluids. For instance, the composition and make of the nanoparticles may have an impact on the performance of the nanofluid as a heat transfer medium. In this paper we investigate unsteady MHD mixed convection boundary layer with suction/injection subject to a number of source terms including Dufour and Soret effects, heat generation, an applied magnetic field and viscous dissipation. Various numerical and or semi-numerical methods can and have been used to solve the equations that model this type of boundary layer flow. These equations are non-similar and coupled. In this paper we use the spectral relaxation method (SRM) that was recently proposed by Motsa [23]. This spectral relaxation method promises fast convergence with good accuracy, has been successfully used in a limited number of boundary layer flow and heat transfer problems (see [24, 25]). In this paper we discuss the fluid flow and heat transfer as well as highlight the strengths of the solution method.
2 Governing equations
Consider the two-dimensional unsteady laminar MHD mixed convective flow of a nanofluid due to a stretching sheet situated at \(y = 0\) with stretching velocity \(u=ax\), where a is a constant. The temperature and nanoparticle volume fraction at the stretching surface are \(T_{w}\) and \(C_{w}\), respectively, and those of the ambient nanofluid are \(T_{\infty}\) and \(C_{\infty}\), respectively. The x and y directions are in the plane of and perpendicular to the sheet, respectively. The continuity, momentum, energy and concentration equations of unsteady, incompressible nanofluid boundary layer flow are as follows (see Yang [26]):
where t, u and v are the time, the fluid velocity and the normal velocity components in the x and y orientations, respectively; \(\nu_{nf}\), p, \(\rho_{nf}\), σ, \(B_{0}\), \(\mu_{nf}\), g are the nanofluid kinematic viscosity, the pressure, nanofluid density, electrical conductivity, the uniform magnetic field in the y direction, the effective dynamic viscosity of the nanofluid and gravitational acceleration, respectively; \(\beta_{T}\), \(\beta_{C}\), T, C, \(\alpha _{nf}\), \((\rho cp)_{nf}\), Q are the volumetric thermal expansion coefficient, the solutal expansion coefficient, the temperature of the fluid in the boundary layer, fluid solutal concentration, the thermal diffusivity of the nanofluid, the nanofluid heat capacitance and the volumetric rate of heat generation, respectively; \(\rho_{f}\), \(D_{m}\), \(K_{T}\), \(C_{s}\), \((cp)_{nf}\), \(T_{m}\), R are the density of the base fluid, the mass diffusivity of concentration, thermal diffusion ratio, concentration susceptibility, specific heat of the fluid at constant pressure, mean fluid temperature and the chemical reaction parameter, respectively.
The boundary conditions are as follows:
and the initial conditions are
where \(a_{\infty}\) (>0) is the stagnation flow rate parameter, \(a< 0\) for a shrinking surface and \(a > 0\) for a stretching surface. Here \(v_{w}\) is prescribed suction velocity (\(v_{w} < 0\)) or blowing velocity (\(v_{w} >0\)).
In the free stream the momentum equation (2.2) becomes
Substituting (2.7) in (2.2) the momentum equation is written as
The effective dynamic viscosity of the nanofluid was given by Brinkman [27] as
where Ï• is the solid volume fraction of nanoparticles, \(\mu_{f}\) is the dynamic viscosity of the base fluid. In equations (2.1)-(2.4),
where \(k_{nf}\) is the thermal conductivity of the nanofluid, \(k_{f}\) and \(k_{s}\) are the thermal conductivities of the fluid and of solid fractions, respectively, and \(\rho_{s}\) is the density of solid fractions, \((\rho c_{p})_{f}\) and \((\rho c_{p})_{s}\) are the heat capacity of the base fluid and the effective heat capacity of a nanoparticle, respectively, \(k_{nf}\) is the thermal conductivity of the nanofluid.
The continuity equation (2.1) is satisfied by introducing a stream function \(\psi(x,y)\) such that
We introduce the following non-dimensional variables (see Liao [28]):
where \(f(\xi,\eta)\) is a dimensionless stream function, \(\theta(\xi,\eta)\) is the dimensionless temperature and \(\phi(\xi,\eta)\) is the dimensionless solute concentration. By using (2.11) and (2.12), the governing equations (2.3), (2.4) and (2.8) along with the boundary conditions (2.5) are reduced to the following two-point boundary value problem:
The boundary conditions are as follows:
where primes denote differentiation with respect to η, \(\alpha_{f}=k_{f}/(\rho c_{p})_{f}\) and \(\nu_{f}=\mu_{f}/\rho_{f}\) are the thermal diffusivity and kinetic viscosity of the base fluid, respectively. Other non-dimensional parameters appearing in equations (2.13) to (2.15) are Ha, \(\mathit{Gr}_{t}\), \(\mathit{Gr}_{c}\), \(\mathit{Gr}_{c}\), Pr, δ, \(D_{f}\), Sc, γ and Sr, and they denote the Hartmann number, the local temperature Grashof number, the local concentration Grashof number, the Prandtl number, the dimensionless heat generation parameter, the Dufour number, the Schmidt number, the scaled chemical reaction parameter and the Soret number, respectively. These parameters are defined mathematically as
The boundary conditions are as follows:
The nanoparticle volume fractions \(\phi_{1}\) and \(\phi_{2}\) are defined as
In equations (2.18), \(f_{w} =-v_{w}/\sqrt{a_{\infty}\nu_{f} \xi}\) represents suction (\(f_{w}>0\)) or injection (\(f_{w}<0\)) and λ (\(=a/a_{\infty}\)) is the stretching/shrinking parameter.
3 Skin friction, heat and mass transfer coefficients
The skin friction coefficient \(C_{f}\), the local Nusselt number \(\mathit{Nu}_{x}\) and the local Sherwood number \(\mathit{Sh}_{x}\) characterize the surface drag, wall heat and mass transfer rates, respectively.
The shearing stress at the surface of the wall \(\tau_{w}\) is defined as
where \(\mu_{nf}\) is the coefficient of viscosity.
The skin friction coefficient is obtained as
and using equation (3.1) in (3.2) we obtain
The heat transfer rate at the surface flux at the wall is defined as
where \(k_{nf}\) is the thermal conductivity of the nanofluid. The local Nusselt number is defined as
Using equation (3.4) in equation (3.5), the dimensionless wall heat transfer rate is obtained as
The mass flux at the wall surface is defined as
and the local Sherwood number (mass transfer coefficient) is obtained as
The dimensionless wall mass transfer rate is obtained as
where \(\mathit{Re}_{x}\) represents the local Reynolds number and is defined as
4 Cases of special interest
In this section we highlight two particular cases where equations (2.12) to (2.14) reduce to ordinary differential equations.
4.1 Initial steady flow
For steady flow and a regular fluid, if we assume that \(\xi\to0\), where \(0< \xi\le1\), then \(t \approx0\). Thus \(f(\eta,\xi)\approx f(\eta)\), \(\theta(\eta,\xi) \approx\theta(\eta)\) and \(\varPhi (\eta,\xi)\approx \varPhi (\eta)\). In this case equations (2.12) to (2.14) reduce to
subject to the appropriately modified boundary conditions (2.18). The exact solutions of these equations cannot be easily obtained. The numerical solutions were obtained using the spectral relaxation method (SRM).
4.2 Final steady state flow
In this case, we have \(\xi= 1\) when \(t \rightarrow\infty\), corresponding to \(f(\eta,1)= f(\eta)\), \(\theta(\eta,1) = \theta(\eta)\) and \(\varPhi (\eta,1)= \varPhi (\eta)\). Equations (2.12) to (2.14) reduce to the following forms:
subject to the boundary conditions (2.18). Equations (4.1) to (4.6) were solved using the SRM, Motsa [23].
The spectral relaxation method (SRM) is an iterative procedure that employs the Gauss-Seidel type of relaxation approach to linearize and decouple the system of differential equations. Further details of the rules of the SRM can be found in [24, 25]. The linear terms in each equation are evaluated at the current iteration level (denoted by \(r+1\)) and the non-linear terms are assumed to be known from the previous iteration level (denoted by r). The linearized form of (2.13)-(2.15) is
where
Equations (4.7)-(4.9) are now linear and decoupled. The equations can be solved sequentially to obtain approximate solutions for \(f(\eta,\xi)\), \(\theta(\eta,\xi)\) and \(\phi(\eta,\xi)\). In this study, the Chebyshev spectral collocation method was used to discretize in η and finite differences used to discretize in ξ directions. Starting from initial guesses for f, θ and ϕ, equations (4.7)-(4.9) were solved iteratively until the approximate solutions converged within a certain prescribed tolerance level. The accuracy of the results was validated against results from the literature for some special cases of the governing equations.
5 Results and discussion
The system of partial differential equations (2.13) to (2.15) subject to boundary conditions (2.18) were solved numerically using the spectral relaxation method (SRM) for Cu-water and Ag-water nanofluids. The thermophysical properties of the nanofluids used in the numerical simulations are given in Table 1.
To determine the accuracy of our numerical results, the skin friction coefficient is compared with the published results of Jafar et al. [15], Wang [31] and Suali et al. [32] in Tables 2-4. Here we have varied the stretching parameter while keeping other physical parameters fixed. Table 2 gives a comparison of the SRM results with those obtained by Jafar et al. [15] and Wang [31] when \(\mathit{Ha} = \mathit{Gr}_{t} = \mathit{Gr}_{c} = \delta= D_{f} = \mathit{Sc} = \mathit{Sr} = \gamma= \phi= 0\), \(\mathit{Pr} = 1\) and \(\xi=1\) for different values of the stretching/shrinking parameter. It is observed that for increasing λ, the present results are in good agreement with results in the literature.
Table 4 gives the skin friction coefficient for selected stretching λ parameter values. Here we note that as the stretching rate decreases, the skin friction coefficient increases. These results are in good agreement with those obtained by Suali et al. [32].
The effects of the nanoparticle volume fraction on the fluid velocity, temperature, concentration profiles as well as skin friction, local Nusselt and Sherwood numbers are given in Figures 1-4. It is evident that the solute concentration, skin friction and the local Nusselt number decrease with increasing nanoparticle volume fraction while the velocity, temperature, and the local Sherwood number increase. This is because with an increase in nanoparticles volume fraction, the thermal conductivity of the nanofluid increases, which reduces the thermal boundary layer thickness and the temperature gradient at the wall.
The axial velocity in the case of an Ag-water nanofluid is comparatively higher than that in the case of a Cu-water nanofluid. The temperature distribution in an Ag-water nanofluid is higher than that in a Cu-water nanofluid and this is explained by the observation that the thermal conductivity of silver is higher than that of copper. The concentration boundary layer thickness is higher for the case of a Cu-water than that for the case of an Ag-water nanofluid.
Figure 3 shows that the skin friction coefficient decreases monotonically with increasing ξ. The result is true for both types of fluids. The maximum value of the skin friction in the case of a Cu-water nanofluid is achieved at a smaller value of ξ in comparison with an Ag-water nanofluid. Furthermore, in this paper it is found that the Ag-water nanofluid shows less drag as compared to the Cu-water nanofluid. The dimensionless wall heat transfer rate and the dimensionless wall mass transfer rate are shown as functions of ξ in Figure 4(a) and (b), respectively. We observe that the wall heat transfer rate decreases while the opposite is true in case of the wall mass transfer rate. The Cu-water nanofluid exhibits higher wall heat transfer rate as compared to the Ag-water nanofluid, while the Cu-water nanofluid exhibits less than the Ag-water nanofluid. The presence of nanoparticle tends to increase the wall heat transfer rate and to reduce the wall mass transfer rates with increasing the values of ξ.
Figures 5-7 show the influence of the Hartmann number on the velocity, temperature, skin friction, the local Nusselt number and the local Sherwood number. The effect of the Hartmann number Ha is to increase the nanofluid velocity and the wall heat transfer rate,whereas it reduced the skin friction coefficient and the wall mass transfer rate. A similar observation was made by Jafar et al. [15]. The momentum boundary layer thickness increases with increase in the Hartmann number.
Figure 6 shows the skin friction coefficient as a function of ξ. It is clear that for Ag-water and Cu-water nanofluids, the skin friction reduces when ξ increases. We note that the Cu-water nanofluid exhibits higher drag to the flow as compared to the Ag-water nanofluid. Figure 7 shows the wall heat and mass transfer rates for a different Hartmann number Ha, it is clear that the value of wall heat transfer rate increases as ξ increases, in the case of an Ag-water nanofluid it is less than in the case of a Cu-water nanofluid. Further, the wall mass transfer rate increases up to the value of ξ before reducing.
Figures 8-11 show the velocity, temperature, concentration of nanofluid with skin friction, the wall heat and mass transfer rates for various values of the suction/injection parameter. We observe that the velocity boundary layer thickness decreases with increasing values of the suction parameter. This is because due to suction, the fluid is removed from the system which reduces the momentum boundary layer thickness. Similarly, the boundary layer thickness increases with increase of the injection parameter as injection allows the fluid to enter the system. The thermal boundary layer thickness decreases due to injection, while it increases with suction. The effect of the suction/injection parameter is to increase the concentration profile at the surface. Beyond this critical value, the concentration profile decreases with increasing suction/injection. The solute concentration boundary layer thickness is larger for the case of a Cu-water nanofluid than that for the case of an Ag-water nanofluid (see Figure 9). The skin friction coefficient decreases with increasing the values of ξ. It is obvious that the skin friction for the case of an Ag-water nanofluid is relatively less than that for the case of a Cu-water nanofluid (see Figure 10).
The axial distributions of the wall heat and mass transfer rates are shown in Figure 11(a) and (b), respectively. The wall heat transfer rate increased with ξ, and we observe that the heat transfer rate is higher for a Cu-water nanofluid than for an Ag-water nanofluid. It is interesting to note that with suction (\(f_{w}=-1\)), the heat transfer rate is less for an Ag-water nanofluid than for a Cu-water nanofluid up to a certain value of ξ. Beyond this point, the heat transfer rate is higher for an Ag-water nanofluid as compared to a Cu-water nanofluid, while the wall mass transfer rate increases monotonically with ξ to a maximum values before reducing. It is shown that the mass transfer rate is higher for a Cu-water nanofluid than for an Ag-water nanofluid. The opposite behavior is observed in the case of suction when \(f_{w} = -1\). The mass transfer rate for an Ag-water nanofluid is higher than that for a Cu-water nanofluid up to a certain value of ξ, and beyond this critical value, the mass transfer in an Ag-water nanofluid is less than that in a Cu-water nanofluid, Figure 11(b).
The influence of stretching/shrinking on velocity, temperature, solutal concentration profiles, the skin friction coefficient, wall heat and mass transfer rates are shown in Figures 12-15. Figure 12 shows that the momentum boundary layer thickness increases with the stretching/shrinking rate. This may be attributed to the fact than an increase in the stretching parameter enhances the velocity of the nanofluid which in turn enhances the momentum boundary layer thickness.
Figure 13 shows that the thermal and concentration boundary layer thicknesses decrease as the stretching rate increases. For the shrinking case, when \(\lambda= -2\), the momentum boundary layer for an Ag-water nanofluid is greater than that for a Cu-water nanofluid, while the opposite is observed for the stretching case when \(\lambda = 2\). The Ag-water nanofluid thermal boundary is higher than that of a Cu-water nanofluid (see Figure 13(a)). The solutal concentration increases up to a critical η, and beyond this critical value the concentration profile decreases (see Figure 13(b)). We observe that solute concentration profiles are larger for the case of a Cu-water than those for the case of an Ag-water nanofluid for the shrinking sheet with \(\lambda= -2\), while the opposite is true for the stretching sheet with \(\lambda= 2\).
Figure 14 shows the effect of stretching/shrinking on the shear stress, while Figure 15 shows the effect of the stretching rate on the wall heat and mass transfer rates. From Figure 14 we note that the shear stress increases with the stretching/shrinking parameter. The shear stress decreases with ξ. Figure 15(a) shows that the heat transfer rate increases with increasing λ. The mass transfer at the wall decreases with the increase in λ. The heat transfer rate is larger for the case of an Ag-water nanofluid compared to that of a Cu-water nanofluid, while the opposite is true for the mass transfer rate (see Figure 15).
6 Conclusions
We have investigated heat and mass transfer in unsteady MHD mixed convection in a nanofluid due to a stretching/shrinking sheet with heat generation and viscous dissipation. Other parameters of interest in this study included the Soret and Dufour effects. In this paper we considered Cu-water and Ag-water nanofluids and assumed that the nanoparticle volume fraction can be actively controlled at the boundary surface. We have solved the model equations using the spectral relaxation method, and to benchmark our solutions, we compared our results with some limiting cases from the literature. These results were found to be in good agreement.
The numerical simulations show, inter alia, that the skin friction factor increases with both an increase in the nanoparticle volume fraction and the stretching rate and that an increase in the nanoparticle volume fraction leads to a reduction in the wall mass transfer rate.
References
Kuznetsov, AV, Nield, DA: Natural convective boundary layer flow of a nanofluid past a vertical plate. Int. J. Therm. Sci. 49, 243-247 (2010)
Ishak, A, Nazar, R, Pop, I: Hydromagnetic flow and heat transfer adjacent to a stretching vertical sheet. Heat Mass Transf. 44, 921-927 (2008)
Mahapatra, TR, Mondal, S, Pal, D: Heat transfer due to magnetohydrodynamic stagnation-point flow of a power-law fluid towards a stretching surface in the presence of thermal radiation and suction/injection. ISRN Thermodyn. 9, 1-9 (2012)
Das, SK, Choi, SUS, Yu, W, Pradeep, T: Nanofluids: Science and Technology. Wiley, New York (2007)
Crane, LJ: Flow past a stretching plate. Z. Angew. Math. Phys. 21, 645-647 (1970)
Khan, WA, Pop, I: Boundary layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 53, 2477-2483 (2010)
Chen, C-H: Laminar mixed convection adjacent to vertical, continuously stretching sheets. Heat Mass Transf. 33, 471-476 (1988)
Abo-Eldahab, EM, Abd El-Aziz, M: Blowing/suction effect on hydromagnetic heat transfer by mixed convection from an inclined continuously stretching surface with internal heat generation/absorption. Int. J. Therm. Sci. 43, 709-719 (2004)
Abd El-Aziz, M: Thermal-diffusion and diffusion-thermo effects on combined heat and mass transfer by hydromagnetic three-dimensional free convection over a permeable stretching surface with radiation. Phys. Lett. 372(3), 263-272 (2007)
Hady, FM, Ibrahim, FS, Abdel-Gaied, SM, Eid, MR: Radiation effect on viscous flow of a nanofluid and heat transfer over a nonlinearly stretching sheet. Nanoscale Res. Lett. 7, 229 (2012)
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)
Rohni, AM, Ahmad, S, Ismail, AIM, Pop, I: Flow and heat transfer over an unsteady shrinking sheet with suction in a nanofluid using Buongiorno’s model. Int. Commun. Heat Mass Transf. 43, 75-80 (2013)
Buongiorno, J: Convective transport in nanofluids. J. Heat Transf. 128, 240-250 (2006)
Pavlov, KB: Magnetohydromagnetic flow of an incompressible viscous fluid caused by deformation of a surface. Magn. Gidrodin. 4, 146-147 (1974)
Jafar, K, Nazar, R, Ishak, A, Pop, I: MHD flow and heat transfer over stretching/shrinking sheets with external magnetic field, viscous dissipation and Joule effects. Can. J. Chem. Eng. 99, 1-11 (2011)
Chamka, AJ: MHD flow of a uniformly stretched vertical permeable surface in the presence of heat generation/absorption and a chemical reaction. Int. Commun. Heat Mass Transf. 30, 413-422 (2003)
Kandasamy, R, Palanimani, PG: Effects of chemical reactions, heat, and mass transfer on nonlinear magnetohydrodynamic boundary layer flow over a wedge with a porous medium in the presence of ohmic heating and viscous dissipation. J. Porous Media 10, 489-502 (2007)
Rashidi, MM, Erfani, E: Analytical method for solving steady MHD convective and slip flow due to a rotating disk with viscous dissipation and Ohmic heating. Eng. Comput. 29, 562-579 (2012)
Rashidi, MM, Momoniat, E, Rostami, B: Analytic approximate solutions for MHD boundary-layer viscoelastic fluid flow over continuously moving stretching surface by homotopy analysis method with two auxiliary parameters. J. Appl. Math. 2012, Article ID 780415 (2012)
Rashidi, MM, Keimanesh, M: Using differential transform method and Padé approximant for solving MHD flow in a laminar liquid film from a horizontal stretching surface. Math. Probl. Eng. 2010, Article ID 491319 (2010)
Anjali-Devi, SP, Thiyagarajan, M: Steady nonlinear hydromagnetic flow and heat transfer over a stretching surface of variable temperature. Heat Mass Transf. 42, 671-677 (2006)
Postelnicu, A: Influence of chemical reaction on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects. Heat Mass Transf. 43, 595-602 (2007)
Motsa, SS: A new spectral relaxation method for similarity variable nonlinear boundary layer flow systems. Chem. Eng. Commun. 201, 241-256 (2014)
Motsa, SS, Dlamini, PG, Khumalo, M: Spectral relaxation method and spectral quasilinearization method for solving unsteady boundary layer flow problems. Adv. Math. Phys. 2014, Article ID 341964 (2014). doi:10.1155/2014/341964
Motsa, SS, Makukula, ZG: On spectral relaxation method approach for steady von Karman flow of a Reiner-Rivlin fluid with Joule heating and viscous dissipation. Cent. Eur. J. Phys. 11, 363-374 (2013)
Yang, KT: Unsteady laminar boundary layers in an incompressible stagnation flow. J. Appl. Mech. 25, 421-427 (1958)
Brinkman, HC: The viscosity of concentrated suspensions and solution. J. Chem. Phys. 20, 571-581 (1952)
Liao, SJ: An analytic solution of unsteady boundary layer flows caused by an impulsively stretching plate. Commun. Nonlinear Sci. Numer. Simul. 11, 326-329 (2006)
Sheikholeslami, M, Bandpy, MG, Ganji, DD, Soleimani, S, Seyyedi, SM: Natural convection of nanofluids in an enclosure between a circular and a sinusoidal cylinder in the presence of magnetic field. Int. Commun. Heat Mass Transf. 39, 1435-1443 (2012)
Oztop, HF, Abu-Nada, E: Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 29, 1326-1336 (2008)
Wang, CY: Stagnation flow towards a shrinking sheet. Int. J. Non-Linear Mech. 43, 377-382 (2008)
Suali, M, Long, NMAN, Ishak, A: Unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with prescribed surface heat flux. Appl. Math. Comput. Intell. 1, 1-11 (2012)
Acknowledgements
The authors wish to thank the University of KwaZulu-Natal for financial support.
Author information
Authors and Affiliations
Corresponding author
Additional information
Abbreviations
a: a positive constant; \(B_{0}\): magnetic field; \(C_{\infty}\): ambient concentration; \(C_{fx}\): skin friction coefficient; C: solutal concentration; \(C_{w}\): surface concentration; \(C_{s}\): concentration susceptibility; \(D_{m}\): concentration mass diffusivity; \(D_{f}\): Dufour number; f: dimensionless stream function; g: acceleration due to gravity; \(\mathit{Gr}_{c}\): solutal Grashof number; \(\mathit{Gr}_{t}\): thermal Grashof number; Ha: Hartmann number; \(k_{f}\): fluid thermal conductivity; \(k_{s}\): solid volume fraction; \(k_{nf}\): nanofluid thermal conductivity; \(K_{T}\): thermal diffusion ratio; \(\mathit{Nu}_{x}\): local Nusselt number; p: fluid pressure; Pr: Prandtl number; \(q_{w}\): wall heat flux; \(q_{m}\): wall mass flux; Q: volumetric rate of heat generation; R: chemical reaction parameter; \(\mathit{Re}_{x}\): local Reynolds number; Sc: Schmidt number; \(\mathit{Sh}_{x}\): local Sherwood number; Sr: Soret number; t: time; T: fluid temperature; \(T_{m}\): mean fluid temperature; \(T_{w}\): wall temperature; \(T_{\infty}\): ambient temperature; u, v: fluid velocity components; \(v_{w}\): suction velocity. Greek symbols: \(\alpha_{nf}\): nanofluid thermal diffusivity; \(\beta_{C}\): volumetric solutal expansion coefficient; \(\beta_{T}\): volumetric thermal expansion coefficient; δ: heat generation parameter; \(\mu_{f}\): base fluid dynamic viscosity; \(\mu_{nf}\): nanofluid effective dynamic viscosity; \(\nu_{nf}\): nanofluid kinematic viscosity; \(\rho_{f}\): density of the base fluid; \(\rho_{nf}\): nanofluid density; \(\rho_{s}\): density of the solid fractions; σ: electrical conductivity; \(\tau_{w}\): wall shear stress; ϕ: nanoparticle solid volume fraction; ψ: stream function.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All the authors participated in the design of this study and helped to draft and proofread the manuscript. NAHH produced the initial draft.
Rights and permissions
Open Access This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.
About this article
Cite this article
Haroun, N.A., Sibanda, P., Mondal, S. et al. On unsteady MHD mixed convection in a nanofluid due to a stretching/shrinking surface with suction/injection using the spectral relaxation method. Bound Value Probl 2015, 24 (2015). https://doi.org/10.1186/s13661-015-0289-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-015-0289-5