Skip to main content

A new numerical approach for MHD laminar boundary layer flow and heat transfer of nanofluids over a moving surface in the presence of thermal radiation

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 skin-friction coefficient and the heat and mass transfer rates are depicted graphically and/or in tabular form.

MSC: 65PXX, 76-XX.

1 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 radiation-convection 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 two-dimensional, steady magnetohydrodynamic flow and heat transfer analysis of a non-Newtonian 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 semi-infinite 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 U w =λU, where λ is the velocity parameter. The flow is assumed to take place at y0, 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 T w and C w , respectively, while the free stream values are taken to be T and C , 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:

u x + v y =0,
(1)
u u x +v u y =ν 2 u y 2 σ B 0 2 ρ f u,
(2)
u T x +v T y =α 2 T y 2 α k q r y +τ [ D B C y T y + ( D T T ) ( T y ) 2 ] ,
(3)
u C x +v C y = D B 2 C y 2 + D B T 2 T y 2 ,
(4)

where u and v are the velocity components along the x-axis and y-axis, respectively, α=k/ ( ρ c ) f is the thermal diffusivity of the fluid, ν is the kinematic viscosity coefficient, k is the thermal conductivity, q r is the heat flux, D B is the Brownian diffusion coefficient, D T is the thermophoresis diffusion coefficient, B 0 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

v=0,u=λU,T= T w ,C= C w at y=0,
(5)
uU,T T ,C C as y.
(6)

The surface moving parameter λ>0 corresponds to the downstream movement of the plate from the origin, while λ<0 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 q r is given by

q r = 4 σ 3 K s T 4 y ,
(7)

where σ and K s are the Stefan-Boltzman constant and the Rosseland mean absorption coefficient, respectively. We assume that the temperature differences within the flow are sufficiently small such that T 4 may be expressed as a linear function of temperature.

T 4 4 T 3 T3 T 4 .
(8)

Using (7) and (8) in the last term of equation (3), we obtain

q r y = 16 σ T 3 3 K s 2 T y 2 .
(9)

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:

ψ = ( 2 U ν x ) 1 2 f ( η ) , θ ( η ) = T T T w T , ϕ ( η ) = C C C w C , η = ( U / 2 ν x ) 1 2 y ,
(10)

where f(η) 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 u= ψ y and v= ψ x . Substituting the similarity variables into equations (2)-(4) gives

f +f f Ha f =0,
(11)
( 3 + 4 R Pr ) θ +f θ +Nb ϕ θ +Nt θ 2 =0,
(12)
ϕ +Lef ϕ + Nt Nb θ =0,
(13)

subject to the boundary conditions

f(0)=0, f (0)=λ,θ(0)=1,ϕ(0)=1,
(14)
f 1,θ0,ϕ0as η,
(15)

where R= 4 α δ T 3 k k is the radiation parameter, Pr= ν α is the Prandtl number, Le= ν D B is the Lewis number, Ha= 2 x B 0 2 U ρ f is the local Hartman number, Nb= ( ρ c ) p D B ( ϕ w ϕ ) ( ρ C ) f ν is the Brownian motion parameter and Nt= ( ρ C ) p ρ T ( T w T ) ( ρ c ) f T ν is the thermophoresis parameter.

The quantities of engineering interest are the skin-friction coefficient C f , the local Nusselt number Nu x , and the local Sherwood number Sh x . These quantities are defined as follows:

C f = τ w ρ u 2 , Nu x = q w k ( T w T ) , Sh x = x q m D B ( C w C ) ,
(16)

where τ w , q w and q m are the shear stress, heat flux and mass flux at the surface, respectively. Upon using the similarity variables into the above expressions, we get

( 2 Re ) 1 2 C f = f (0), ( Re x 2 ) 1 2 Nu x = θ (0), ( Re x 2 ) 1 2 Sh x = ϕ (0),
(17)

where Re x =Ux/ν is the local Reynolds number.

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 pseudo-spectral 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 skin-friction 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, z i (η), i=1,2,,m, where η[a,b] is the independent variable. To solve the resultant iterative scheme, we then use the Chebyshev pseudo-spectral 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 [1,1] on which the spectral method can then be implemented. We use the transformation η=(ba)(τ+1)/2 to map the interval [a,b] to [1,1].

To apply the SRM to the nonlinear ordinary differential equations, we first set f (η)=g(η) and then write the equations as the following set of equations:

f =g,
(18)
g +f g g 2 Hag=0,
(19)
( 3 + 4 R Pr ) θ +f θ +Nb θ ϕ +Nt θ 2 =0,
(20)
ϕ +Lef ϕ + Nt Nb θ =0,
(21)

and the boundary conditions become

f ( 0 ) = 0 , g ( 0 ) = λ , θ ( 0 ) = ϕ ( 0 ) = 1 , g ( ) = 1 , θ ( ) = 0 , ϕ ( ) = 0 .
(22)

In view of the SRM, we obtain the following iterative scheme:

f r + 1 = g r , f r + 1 (0)=0,
(23)
g r + 1 + f r + 1 g r + 1 Ha g r + 1 =0, g r + 1 (0)=λ, g r + 1 ()=1,
(24)
( 3 + 4 R Pr ) θ r + 1 + f r + 1 θ r + 1 +Nb θ r + 1 = Nt r + 1 2 , θ r + 1 (0)=1, θ r + 1 ()=0,
(25)
ϕ r + 1 +Le f r + 1 ϕ r + 1 = Nt Nb θ r + 1 , ϕ r + 1 (0)=1, ϕ r + 1 ()=0.
(26)

The above equations form a system of linear decoupled equations which can be solved iteratively for r=1,2, , starting from initial guesses/approximations ( g 0 (η),θ(η),ϕ(η)).

Applying the Chebyshev pseudo-spectral method to (23)-(26), we obtain

A 1 f r + 1 = B 1 , f r + 1 ( τ N ¯ )=0,
(27)
A 2 g r + 1 = B 2 , g r + 1 ( τ N ¯ )=λ, g r + 1 ( τ 0 )=1,
(28)
A 3 θ r + 1 = B 3 , θ r + 1 ( τ N ¯ )=1, θ r + 1 ( τ 0 )=0,
(29)
A 4 ϕ r + 1 = B 4 , ϕ r + 1 ( τ N ¯ )=1, ϕ r + 1 ( τ 0 )=0,
(30)

where

A 1 = D , B 1 = gr , A 2 = D 2 + diag [ f i + 1 ] D Ha I , B 2 = 0 , A 3 = ( 3 + 4 R Pr ) D 2 + diag [ f i + 1 + Nb ] D , B 3 = Nt θ r + 1 2 , A 4 = D 2 + diag [ Le f i + 1 ] D , B 4 = Nt Nb θ r + 1 ,

where I is the identity matrix of size ( N ¯ +1)×( N ¯ +1), 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

g 0 (η)=λ1+η+(1λ) e η ,
(31)
θ 0 (η)= e η ,
(32)
ϕ 0 (η)= e η ,
(33)

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

Er=Max ( f r + 1 f r ; g r + 1 g r ; θ r + 1 θ r ; ϕ r + 1 ϕ r ) .
(34)

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 N=50 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 Pr=0.71, Nt=Nb=0.3, Le=2, R=2, λ=0.1.

Tables 1 through 5 show a comparison among the SRM results to those previously obtained results as well as the b v p 4 c 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 for θ (0) and ϕ (0) against those of the b v p 4 c as well as those obtained by Olanrewaju et al. [5] for different values of Pr
Table 2 Comparison of SRM solutions for θ (0) and ϕ (0) against those of the b v p 4 c as well as those obtained by Olanrewaju et al. [5] for different values of Le
Table 3 Comparison of SRM solutions for θ (0) and ϕ (0) against those of the b v p 4 c as well as those obtained by Olanrewaju et al. [5] for different values of λ
Table 4 Computations of θ (0) and ϕ (0) at different values of Nb and Nt
Table 5 Comparison of SRM solutions for f (0) against those of the b v p 4 c for 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 1-9. 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
figure 1

Graph of the SRM solutions of the velocity profiles for different values of Ha .

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
figure 2

Graph of the SRM solutions of the velocity profiles for different values of λ .

Figure 3
figure 3

Temperature profiles for different values of Ha .

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
figure 4

Graph of the SRM solutions of the temperature profiles for different values of R .

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
figure 5

Influence of Nt on the temperature profiles.

Figure 6
figure 6

Graph of the SRM solutions of the temperature profiles for different values of Nb .

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
figure 7

Graph of the SRM solutions of the nanoparticle volume fraction profiles for different values of Ha .

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.

Figure 8
figure 8

Influence of Le on the nanoparticle volume fraction profiles.

Figure 9
figure 9

Influence of Nt on the nanoparticle volume fraction profiles.

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 in-built b v p 4 c 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.

Greek letters

Nomenclature

References

  1. Choi SUS: Enhancing thermal conductivity of fluids with nanoparticles. FED 231. Developments and Applications of Non-Newtonian Flows 1995, 99-105. MD, vol. 66

    Google Scholar 

  2. Mnyusiwalla A, Daar AS, Singer PA: Mind the gap: science and ethics in nanotechnology. Nanotechnology 2003, 14: R9-R13. 10.1088/0957-4484/14/3/201

    Article  Google Scholar 

  3. Kakaç S, Pramuanjaroenkij A: Review of convective heat transfer enhancement with nanofluids. Int. J. Heat Mass Transf. 2009, 52: 3187-3196. 10.1016/j.ijheatmasstransfer.2009.02.006

    Article  Google Scholar 

  4. Godson L, Raja B, Lal DM, Wongwises S: Enhancement of heat transfer using nanofluid - An overview. Renew. Sustain. Energy Rev. 2010, 14: 629-641. 10.1016/j.rser.2009.10.004

    Article  Google Scholar 

  5. 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. 2012, 2(1):274-285.

    Google Scholar 

  6. 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. 2011, 5(9):1323-1334.

    Google Scholar 

  7. Khan WA, Aziz A: Natural convection flow of a nanofluid over a vertical plate with uniform surface heat flux. Int. J. Therm. Sci. 2011, 50: 1207-1214. 10.1016/j.ijthermalsci.2011.02.015

    Article  Google Scholar 

  8. Makinde OD, Aziz A: Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition. Int. J. Therm. Sci. 2011, 50: 1326-1332. 10.1016/j.ijthermalsci.2011.02.019

    Article  Google Scholar 

  9. Khan MS, Karim I, Ali LE, Islam A: Unsteady MHD free convection boundary-layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects. Int. Nano Lett. 2012., 2: Article ID 24 10.1186/2228-5326-2-24

    Google Scholar 

  10. 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. 2012, 7: 229-242. 10.1186/1556-276X-7-229

    Article  Google Scholar 

  11. 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. 2011, 61(11):3391-3399. 10.1016/j.camwa.2011.04.053

    Article  MathSciNet  Google Scholar 

  12. 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. 2012, 31(3):125-132.

    Google Scholar 

  13. Khan Y, Wu Q, Faraz N, Yildirim A, Mohyud-Din ST: Heat transfer analysis on the magnetohydrodynamic flow of a non-Newtonian fluid in the presence of thermal radiation: an analytic solution. Z. Naturforsch. A 2012, 67a: 147-152.

    Google Scholar 

  14. Motsa SS, Dlamini P, Khumalo M: A new multistage spectral relaxation method for solving chaotic initial value systems. Nonlinear Dyn. 2013, 72: 265-283. 10.1007/s11071-012-0712-8

    Article  MathSciNet  Google Scholar 

  15. Motsa SS, Dlamini PG, Khumalo M: Solving hyperchaotic systems using the spectral relaxation method. Abstr. Appl. Anal. 2012. 10.1155/2012/203461

    Google Scholar 

  16. Motsa SS, Dlamini P, Khumalo M: A new multistage spectral relaxation method for solving chaotic initial value systems. Nonlinear Dyn. 2013. 10.1007/s11071-012-0712-8

    Google Scholar 

  17. Motsa SS, Makukula ZG: On spectral relaxation method approach for steady von Kármán flow of a Reiner-Rivlin fluid with Joule heating, viscous dissipation and suction/injection. Cent. Eur. J. Phys. 2013. 10.2478/s11534-013-0182-8

    Google Scholar 

  18. Shateyi S, Makinde OD: Hydromagnetic stagnation-point flow towards a radially stretching convectively heated disk. Math. Probl. Eng. 2013. 10.1155/2013/616947

    Google Scholar 

  19. 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. 10.1155/2013/670205

    Google Scholar 

  20. 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. 10.1186/1687-2770-2013-196

    Google Scholar 

  21. Tiwari RK, Das MK: Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50: 2002-2018. 10.1016/j.ijheatmasstransfer.2006.09.034

    Article  Google Scholar 

  22. Hossain MA, Takhar HS: Radiation effect on mixed convection along a vertical plate with uniform surface temperature. Heat Mass Transf. 1996, 31: 243-248. 10.1007/BF02328616

    Article  Google Scholar 

  23. Raptis A: Flow of a micropolar fluid past a continuously moving plate by the presence of radiation. Int. J. Heat Mass Transf. 1998, 41: 2865-2866. 10.1016/S0017-9310(98)00006-4

    Article  Google Scholar 

  24. Canuto C, Hussaini MY, Quarteroni A, Zang TA: Spectral Methods in Fluid Dynamics. Springer, Berlin; 1988.

    Book  Google Scholar 

  25. Trefethen LN: Spectral Methods in MATLAB. SIAM, Philadelphia; 2000.

    Book  Google Scholar 

  26. Motsa SS: A new spectral relaxation method for similarity variable nonlinear boundary layer flow systems. Chem. Eng. Commun. 2013. 10.1080/00986445.2013.766882

    Google Scholar 

Download references

Acknowledgements

The authors wish to acknowledge financial support from the University of Venda, NRF and University of Botswana.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Stanford Shateyi.

Additional information

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.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Shateyi, S., Prakash, J. A new numerical approach for MHD laminar boundary layer flow and heat transfer of nanofluids over a moving surface in the presence of thermal radiation. Bound Value Probl 2014, 2 (2014). https://doi.org/10.1186/1687-2770-2014-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-2

Keywords