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 u and v are the velocity components along the xaxis and yaxis, respectively,
The boundary conditions are taken as
The surface moving parameter
By using the Rosseland diffusion approximation (Hossain et al.[22]) and following Raptis [23] among other researchers, the radiative heat flux
where
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
subject to the boundary conditions
where
The quantities of engineering interest are the skinfriction coefficient
where
where
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,
To apply the SRM to the nonlinear ordinary differential equations, we first set
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
Applying the Chebyshev pseudospectral method to (23)(26), we obtain
where
where I is the identity matrix of size
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
Tables 1 through 5 show a comparison among the SRM results to those previously obtained results as well
as the
Table 1. Comparison of SRM solutions for
Table 2. Comparison of SRM solutions for
Table 3. Comparison of SRM solutions for
Table 4. Computations of
Table 5. Comparison of SRM solutions for
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
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 MD, vol. 66 (1995)

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, Article ID 24 (2012)
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