Abstract
Numerical analysis has been carried out on the problem of magnetohydrodynamic boundary layer flow of a nanofluid over a moving surface in the presence of thermal radiation. The governing partial differential equations were transformed into a system of ordinary differential equations using suitable similarity transformations. The resultant ordinary equations were then solved using the spectral relaxation method. Effects of the physical parameters on the velocity, temperature and concentration profiles as well as the local skinfriction coefficient and the heat and mass transfer rates are depicted graphically and/or in tabular form.
MSC: 65PXX, 76XX.
Keywords:
MHD boundary layer flow; nanofluids; thermal radiation; moving surface1 Introduction
Many engineering and industrial processes involve heat transfer by means of a flowing fluid in either laminar or turbulent regimes. A decrease in thermal resistance of heat transfer in the fluids would significantly benefit many of these applications/processes. Nanofluids have the potential to reduce thermal resistances, and industrial groups such as electronics, medical, food and manufacturing would benefit from such improved heat transfer. It is well known that conventional heat transfer fluids, such as oil, water and ethylene glycol mixture, are poor heat transfer fluids. Choi [1] introduced the technique of nanofluids by using a mixture of nanoparticles and the base fluids. The presence of the nanoparticles in the nanofluid increases the thermal conductivity and therefore substantially enhances the heat transfer characteristics of the nanofluid. Nanotechnology has been an ongoing hot topic of discussion in public health as researchers claim that nanoparticles could present possible dangers in health and environment, Mnyusiwella et al.[2]. There are several numerical studies on the modelling of natural convection heat transfer in nanofluids (Kaka and Pramuanjaroekij [3], Godson et al.[4], Olanrewaju et al.[5]). Gbadeyan et al.[6] numerically studied boundary layer flow induced in a nanofluid due to a linearly stretching sheet in the presence of thermal radiation and induced magnetic field. Khan and Aziz [7] studied natural convection flow of a nanofluid over a vertical plate with a uniform surface heat flux. Makinde and Aziz [8] investigated boundary layer flow of a nanofluid past a stretching sheet with a convective boundary conditions. Khan et al.[9] studied the unsteady free convection boundary layer flow of a nanofluid along a stretching sheet with thermal radiation in the presence of a magnetic field.
Thermal radiationconvection interaction problems are found in the cooling of high temperature components design where heat transfer from surfaces occurs by parallel radiation and convection, the interaction of incident solar radiation with the earth’s surface to produce complex free convection patterns. Hadey et al.[10] studied the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation. Khan et al.[11] analyzed the effects of variable viscosity and thermal conductivity on the flow and heat transfer in a laminar liquid film on a horizontal shrinking/stretching sheet. Khan et al.[12] investigated, by employing the homotopy perturbation transform method (HPTM) and the Padé approximation, the problem of magnetohydrodynamic (MHD) boundary layer flow over a nonlinear porous stretching sheet. Khan et al.[13] considered a twodimensional, steady magnetohydrodynamic flow and heat transfer analysis of a nonNewtonian fluid in a channel with a constant wall temperature in the presence of thermal radiation.
As the governing equations modelling MHD flow and heat transfer of nanofluids are highly nonlinear, exact solutions are impossible to obtain. Over the years, numerical methods have been developed, improved, and hybred as a way of getting more accurate solutions. The current study seeks to employ a recently developed numerical technique known as spectral relaxation method [14] to solve the problem of magnetohydrodynamic boundary layer flow of nanofluids over a moving surface in the presence of thermal radiation. The current method has been successfully employed in [1520], among others. We apply this method to the problem of MHD flow of a nanofluid past a stretching sheet in the presence of thermal radiation. The governing boundary layer equations are transformed to a system of nonlinear ordinary differential by using suitable local similarity variables.
2 Mathematical formulation
We consider the steady MHD boundary layer flow of a nanofluid past a moving semiinfinite flat plate in a uniform free stream in the presence of thermal radiation. We assume that the velocity of the uniform stream is U and that of the plate is , where λ is the velocity parameter. The flow is assumed to take place at , with y being the coordinate measured normal to the moving surface. A uniform magnetic field is applied in the y direction. At moving surface, the temperature and the nanoparticles take constant values and , respectively, while the free stream values are taken to be and , respectively. Following the nanofluid model proposed by Tiwari and Das [21] along with Boussinesq and boundary layer approximations, the governing equations for the present problem are:
where u and v are the velocity components along the xaxis and yaxis, respectively, is the thermal diffusivity of the fluid, ν is the kinematic viscosity coefficient, k is the thermal conductivity, is the heat flux, is the Brownian diffusion coefficient, is the thermophoresis diffusion coefficient, is the uniform magnetic field strength of the base fluid, σ is the electrical conductivity of the base fluid, τ is the ratio of the nanoparticle heat capacity and the base fluid heat capacity.
The boundary conditions are taken as
The surface moving parameter corresponds to the downstream movement of the plate from the origin, while corresponds to the upstream movement of the plate.
By using the Rosseland diffusion approximation (Hossain et al.[22]) and following Raptis [23] among other researchers, the radiative heat flux is given by
where and are the StefanBoltzman constant and the Rosseland mean absorption coefficient, respectively. We assume that the temperature differences within the flow are sufficiently small such that may be expressed as a linear function of temperature.
Using (7) and (8) in the last term of equation (3), we obtain
2.1 Similarity transformations
In order to reduce the governing equations into a system of ordinary differential equations, we introduce the following local similarity variables:
where is the dimensionless stream function, η is the similarity variable, is the dimensionless temperature, is the dimensionless nanoparticles concentration. It is worth mentioning that the continuity equation (1) is identically satisfied from our choice of the stream function with and . Substituting the similarity variables into equations (2)(4) gives
subject to the boundary conditions
where is the radiation parameter, is the Prandtl number, is the Lewis number, is the local Hartman number, is the Brownian motion parameter and is the thermophoresis parameter.
The quantities of engineering interest are the skinfriction coefficient , the local Nusselt number , and the local Sherwood number . These quantities are defined as follows:
where , and are the shear stress, heat flux and mass flux at the surface, respectively. Upon using the similarity variables into the above expressions, we get
3 Method of solution
To solve the set of ordinary differential equations (11)(13) together with the boundary conditions (14) and (15), we employ the Chebyshev pseudospectral method known as spectral relaxation method (SRM). This method transforms sets of nonlinear ordinary differential into sets of linear ordinary differential equations. The entire computational procedure is implemented using a program written in MATLAB computer language. The fluid velocity, temperature, the local skinfriction coefficient and the local Nusselt and Sherwood numbers are determined from these numerical computations. The SRM algorithm starts with the assumption of having a system of m nonlinear ordinary differential equations in m unknowns functions, , , where is the independent variable. To solve the resultant iterative scheme, we then use the Chebyshev pseudospectral method. The details of the spectral methods can be found in Canuto et al.[24], Trefethen [25]. Before applying the spectral method, the domain on which the governing equation is defined is transformed to the interval on which the spectral method can then be implemented. We use the transformation to map the interval to .
To apply the SRM to the nonlinear ordinary differential equations, we first set and then write the equations as the following set of equations:
and the boundary conditions become
In view of the SRM, we obtain the following iterative scheme:
The above equations form a system of linear decoupled equations which can be solved iteratively for , starting from initial guesses/approximations .
Applying the Chebyshev pseudospectral method to (23)(26), we obtain
where
where I is the identity matrix of size , f, g, ϕ and θ are the values of f, g, ϕ and θ, respectively, when evaluated at the grid points. Equations (27)(30) constitute the SRM scheme. The initial approximation required to start the iterative process is
which are randomly chosen functions that satisfy the boundary conditions. The iteration is repeated until convergence is achieved. The convergence of the SRM scheme is defined in terms of the infinity norm as
Accuracy of the scheme is established by increasing the number of collocation points N until the solutions are consistent and further increases do not change the value of the solutions.
4 Results and discussion
The system of ordinary differential equations (12)(13) subject to the boundary conditions (14)(15) are numerically solved by using spectral relaxation method (SRM). This is a recently developed method, and details of the method are found in [26]. The SRM results presented in this work were obtained using collocation points, and also the convergence was achieved after as few as five iterations. We also take 15 to be the infinity value . We use these default values for the parameters , , , , .
Tables 1 through 5 show a comparison among the SRM results to those previously obtained results as well as the results. As can be clearly observed from these tables, there are excellent agreements in all the tables giving us confidence in the findings of this study. Table 1 depicts the effects of increasing the values of the Prandtl number on heat and mass transfers. Physically, the Prandtl number is a dimensionless number which is the ratio of momentum diffusivity (kinematic viscosity) to thermal diffusivity. Increasing values of the Prandtl number means that momentum diffusivity dominates thermal diffusivity. Thus, the rate of heat transfer at the surface increases with increasing values of the Prandtl number while the rate of mass transfer is reduced as the Prandtl number increases. Table 2 shows the effect of increasing the Lewis number Le on the heat and mass transfer coefficients. The Lewis number is defined as the ratio of thermal diffusivity to mass diffusivity or the ratio of the Schmidt number to the Prandtl number. As can be clearly seen in Table 2, increasing the Lewis number means that the fluid becomes more viscous therefore causing increases in the rate of mass transfer at the surface. As expected, increasing the Lewis number reduces the rate of heat transfer on the surface as this corresponds to the reduced Prandtl number. In Table 3, we display the effect of increasing the values of the moving surface parameter λ on the rates of heat and mass transfer. These rates both increase with increasing values of the surface moving parameter. Table 4 depicts the effects of Brownian motion parameter Nb and the thermophoresis parameter on the heat and mass transfer coefficients. The rate of heat transfer on the surface as expected is reduced as the Brownian motion parameter increases, while the rate of mass transfer slightly increases as this parameter increases. Also in Table 4, we depict the effect of the thermophoresis parameter Nt on the Nusselt and Sherwood numbers. By definition thermophoresis is the migration of a colloidal particle in a solution in response to a macroscopic temperature gradient. Both the Nusselt number and the Sherwood number are reduced as the values of this parameter increases.
Table 1. Comparison of SRM solutions forandagainst those of theas well as those obtained by Olanrewaju et al.[5]for different values ofPr
Table 2. Comparison of SRM solutions forandagainst those of theas well as those obtained by Olanrewaju et al.[5]for different values ofLe
Table 3. Comparison of SRM solutions forandagainst those of theas well as those obtained by Olanrewaju et al.[5]for different values ofλ
Table 4. Computations ofandat different values ofNbandNt
Table 5. Comparison of SRM solutions foragainst those of thefor different values ofλ
In Table 5 we show the influence of the velocity parameter on the skin friction. Increasing positive values of λ corresponds to accelerated movement of the plate. This in turn results in reduction of skin friction. The effects of physical parameters of importance in this study on the velocity, temperature and nanoparticles distributions are depicted in Figures 19. Figure 1 displays the effect of the Hartman number Ha on the dimensionless velocity. The presence of a magnetic field normal to the flow in an electrically conducting fluid gives rise to a Lorentz force, which acts against the flow. Thus the velocity distributions of the nanofluid are greatly reduced with increasing values of the Hartman number.
Figure 1. Graph of the SRM solutions of the velocity profiles for different values ofHa.
Figure 2 presents the effect of λ on the velocity profiles. Increasing the down stream movement of the plate from the origin, as expected, enhances the fluid velocity. The influence of the magnetic field parameter on the temperature distribution is shown in Figure 3. From this figure, we observe that the temperature profiles increase with increasing values of the magnetic field parameter. This implies that the applied magnetic field tends to heat the fluid, thus reducing the heat transfer from the wall.
Figure 2. Graph of the SRM solutions of the velocity profiles for different values ofλ.
Figure 3. Temperature profiles for different values ofHa.
Figure 4 shows the effect of increasing the radiation parameter R on the temperature distribution. We observe in this figure that the thickness of the thermal boundary layer increases as the values of the radiation parameter are increased. Thus thermal radiation enhances thermal diffusion.
Figure 4. Graph of the SRM solutions of the temperature profiles for different values ofR.
In Figure 5, we have the influence of the thermophoresis parameter on the temperature distribution. Increases in this parameter physically imply high temperature gradients. This in turn results in high temperature distributions within the nanofluid flow. Increasing the values of the Brownian motion parameter Nb results in thickening of the thermal boundary layer. Thus enhancing the temperature of the nanofluid as can be easily seen in Figure 6.
Figure 5. Influence ofNton the temperature profiles.
Figure 6. Graph of the SRM solutions of the temperature profiles for different values ofNb.
The effect of the magnetic field strength Ha on the nanoparticle volume fraction profiles is shown in Figure 7. Application of a normal magnetic field results in slowing down of the nanofluid flow thereby resulting in volumetric increase of the nanoparticles within the flow as less particles are now transported per given time.
Figure 7. Graph of the SRM solutions of the nanoparticle volume fraction profiles for different values ofHa.
In Figure 8 we depict the effect of the Lewis number on the nanoparticles. Large values of the Lewis number implies increased values of the Schmidt number which in turn results in thinning of the solutal boundary layer as can be clearly seen in Figure 8. Increasing the thermophoresis as expected results in increase of the nanoparticle volume fraction profiles. This is displayed in Figure 9.
5 Conclusion
Numerical analysis has been carried out on MHD boundary layer flow of nanofluids over a moving surface in the presence of thermal radiation. The governing partial differential equations were transformed into a system of ordinary differential equations using suitable similarity transformations. The resultant equations were then solved using the spectral relaxation method. The accuracy of the SRM is validated against the MATLAB inbuilt routine for solving boundary value problems as well as previously obtained results. An excellent agreement was observed between our results and those obtained using other methods giving confidence to our present results. We observed that the local temperature rises as the Brownian motion, thermophoresis and radiation effects intensify. Lastly, the Nusselt number decreases while the Sherwood number increases as the Brownian motion and thermophoresis effects are increased.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
SS formulated, generated and discussed the results, and JP discussed the results and proofread the manuscript. All authors read and approved the final manuscript.
Greek letters
Nomenclature
Acknowledgements
The authors wish to acknowledge financial support from the University of Venda, NRF and University of Botswana.
References

Choi, SUS: Enhancing thermal conductivity of fluids with nanoparticles. Developments and Applications of NonNewtonian Flows. 99–105 (1995) MD, vol. 66

Mnyusiwalla, A, Daar, AS, Singer, PA: Mind the gap: science and ethics in nanotechnology. Nanotechnology. 14, R9–R13 (2003). Publisher Full Text

Kakaç, S, Pramuanjaroenkij, A: Review of convective heat transfer enhancement with nanofluids. Int. J. Heat Mass Transf.. 52, 3187–3196 (2009). Publisher Full Text

Godson, L, Raja, B, Lal, DM, Wongwises, S: Enhancement of heat transfer using nanofluid  An overview. Renew. Sustain. Energy Rev.. 14, 629–641 (2010). Publisher Full Text

Olanrewaju, PO, Olanrewaju, MA, Adesanya, AO: Boundary layer flow of nanofluids over a moving surface in a flowing fluid in the presence of radiation. Int. J. Appl. Sci. Technol.. 2(1), 274–285 (2012)

Gbadeyan, JA, Olanrewaju, MA, Olanrewaju, PO: Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition in the presence of magnetic field and thermal radiation. Aust. J. Basic Appl. Sci.. 5(9), 1323–1334 (2011)

Khan, WA, Aziz, A: Natural convection flow of a nanofluid over a vertical plate with uniform surface heat flux. Int. J. Therm. Sci.. 50, 1207–1214 (2011). 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

Khan, MS, Karim, I, Ali, LE, Islam, A: Unsteady MHD free convection boundarylayer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects. Int. Nano Lett.. 2, (2012) Article ID 24
Article ID 24
BioMed Central Full Text 
Hady, FM, Ibrahim, FS, AbdelGaied, SM, Eid, MR: Radiation effect on viscous flow of a nanofluid and heat transfer over a nonlinearly stretching sheet. Nanoscale Res. Lett.. 7, 229–242 (2012). PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Khan, Y, Wua, Q, Faraz, N, Yildirim, A: The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet. Comput. Math. Appl.. 61(11), 3391–3399 (2011). Publisher Full Text

Khan, Y, Abdou, MA, Faraz, N, Yildirim, A, Wu, Q: Numerical solution of MHD flow over a nonlinear porous stretching sheet. Iran. J. Chem. Chem. Eng.. 31(3), 125–132 (2012)

Khan, Y, Wu, Q, Faraz, N, Yildirim, A, MohyudDin, ST: Heat transfer analysis on the magnetohydrodynamic flow of a nonNewtonian fluid in the presence of thermal radiation: an analytic solution. Z. Naturforsch. A. 67a, 147–152 (2012)

Motsa, SS, Dlamini, P, Khumalo, M: A new multistage spectral relaxation method for solving chaotic initial value systems. Nonlinear Dyn.. 72, 265–283 (2013). Publisher Full Text

Motsa, SS, Dlamini, PG, Khumalo, M: Solving hyperchaotic systems using the spectral relaxation method. Abstr. Appl. Anal. (2012). Publisher Full Text

Motsa, SS, Dlamini, P, Khumalo, M: A new multistage spectral relaxation method for solving chaotic initial value systems. Nonlinear Dyn. (2013). Publisher Full Text

Motsa, SS, Makukula, ZG: On spectral relaxation method approach for steady von Kármán flow of a ReinerRivlin fluid with Joule heating, viscous dissipation and suction/injection. Cent. Eur. J. Phys. (2013). Publisher Full Text

Shateyi, S, Makinde, OD: Hydromagnetic stagnationpoint flow towards a radially stretching convectively heated disk. Math. Probl. Eng. (2013). Publisher Full Text

Shateyi, S, Marewo, GT: A new numerical approach of MHD flow with heat and mass transfer for the UCM fluid over a stretching surface in the presence of thermal radiation. Math. Probl. Eng. (2013). Publisher Full Text

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

Tiwari, RK, Das, MK: Heat transfer augmentation in a twosided liddriven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf.. 50, 2002–2018 (2007). Publisher Full Text

Hossain, MA, Takhar, HS: Radiation effect on mixed convection along a vertical plate with uniform surface temperature. Heat Mass Transf.. 31, 243–248 (1996). Publisher Full Text

Raptis, A: Flow of a micropolar fluid past a continuously moving plate by the presence of radiation. Int. J. Heat Mass Transf.. 41, 2865–2866 (1998). Publisher Full Text

Canuto, C, Hussaini, MY, Quarteroni, A, Zang, TA: Spectral Methods in Fluid Dynamics, Springer, Berlin (1988)

Trefethen, LN: Spectral Methods in MATLAB, SIAM, Philadelphia (2000)

Motsa, SS: A new spectral relaxation method for similarity variable nonlinear boundary layer flow systems. Chem. Eng. Commun. (2013). Publisher Full Text