Skip to main content

Thermodiffusion effects on magneto-nanofluid flow over a stretching sheet

Abstract

We study the effect of thermophoresis on boundary layer magneto-nanofluid flow over a stretching sheet. The model includes the effects of Brownian motion and cross-diffusion effects. The governing partial differential equations are transformed to a system of ordinary differential equations and solved numerically using a spectral linearisation method. The effects of the magnetic influence number, the Prandtl number, Lewis number, the Brownian motion parameter, thermophoresis parameter, the modified Dufour parameter and the Dufour-solutal Lewis number on the fluid properties as well as on the heat, regular and nano mass transfer coefficients are determined and shown graphically.

1 Introduction

Most common fluids such as water, ethylene, glycol, toluene or oil generally have poor heat transfer characteristics owing to their low thermal conductivity. A recent technique to improve the thermal conductivity of these fluids is to suspend nano-sized metallic particles such as aluminum, titanium, gold, copper, iron or their oxides in the fluid to enhance its thermal properties, Choi [1]. The enhancement of thermal conductivity in nanofluids has been studied by, among others, Kakac and Pramuanjaroenkij [2], Choi et al. [3], Masuda et al. [4], Eapen et al. [5] and Fan and Wang [6]. Nield and Kuznetsov [7] analyzed the behaviour of boundary layer flow on the Chen-Minkowycz problem in a porous layer saturated with a nanofluid. Nield and Kuznetsov [8] investigated thermal instability in a porous medium saturated with nanofluid using the Brinkman model. The model incorporated the effects of Brownian motion and thermophoresis of nanoparticles. They found that the critical thermal Rayleigh number can be reduced or increased by a substantial amount depending on whether the nanoparticle distribution is top-heavy or bottom-heavy. Aziz et al. [9] studied steady boundary layer flow past a horizontal flat plate embedded in a porous medium filled with a water-based nanofluid containing gyrotactic microorganisms. Cheng [10] investigated the behaviour of boundary layer flow over a horizontal cylinder of elliptic cross section in a porous medium saturated with a nanofluid. Chamkha et al. [11] investigated the non-similar solutions for natural convective boundary layer flow over a sphere embedded in a porous medium saturated with a nanofluid.

During the last few decades, fluid flow over a stretching surface has received considerable attention because of its engineering applications such as in melt-spinning, hot rolling, wire drawing, glass-fiber production and the manufacture of polymer and rubber sheets, Altan and Gegel [12], Fisher [13], and Tidmore and Klein [14]. Nanofluid flow over a stretching surface has been investigated by many researchers. The first study on a stretching sheet in nanofluids was published by Khan and Pop [15]. Makinde and Aziz [16] performed a numerical study of boundary layer flow over a linear stretching sheet. Both Brownian motion and thermophoresis effects on the transport equations were presented. They reported that stronger Brownian motion and thermophoresis lead to an increase in the rate of heat transfer. However, the opposite was observed in the case of the rate of mass transfer. Recent studies in this area include those of Narayana and Sibanda [17] and Kameswaran et al. [18].

Magnetic nanofluids have numerous uses or potential applications in engineering and medicine. Using magnetic nanofluids has the potential to regulate the flow rate and heat transfer by controlling the thermo-magnetic convection current and the fluid velocity (see Shima et al. [19], Ganguly et al. [20]). The effects of a magnetic field on nanofluid flow over a stretching sheet have been investigated by, among others, Bachok et al. [21] and Hanad and Ferdows [22].

The aim of this study is to analyse Dufour and Soret effects in a magneto-nanofluid flow over a stretching sheet. In addition, we study Brownian motion and thermophoresis effects using a spectral linearisation method to obtain numerical solutions of the momentum, energy, concentration and mass fraction equations. The successive linearisation method (SLM) is an accurate method for solving non-linear coupled equations (see [2325]). Recent studies such as [2628] have suggested that the SLM is accurate and converges rapidly to the numerical results when compared to other semi-analytical methods such as the Adomian decomposition method, the variational iteration method and the homotopy perturbation method.

2 Mathematical formulation

Consider two-dimensional nanofluid flow over a linearly stretching sheet with velocity u w =ax, where a is a real positive number. The coordinate system is assumed to define the x-axis along the surface of the sheet and y is the coordinate normal to the surface of the sheet. The surface temperature T w and nanoparticle concentration ϕ ˆ w are higher than the ambient values T and ϕ ˆ , respectively. The governing equations for the problem can be written in the form

u x + v y =0,
(1)
u u x +v u u +σ B 0 2 u=ν 2 u y 2 ,
(2)
u T x +v T y = α m 2 T y 2 +τ [ D B ϕ ˆ y ˆ T y ˆ + D T T ( T y ) 2 ] D T C 2 C y 2 ,
(3)
u C x +v C y = D S 2 C y 2 + D C T 2 T y 2 ,
(4)
u ϕ ˆ x +v ϕ ˆ y = D B 2 ϕ ˆ y 2 + ( D T T ) 2 T y 2 ,
(5)

with the boundary conditions

v = 0 , u = u w ( = a x ) , T = T w , C = C w and ϕ ˆ = ϕ ˆ w on  y = 0 , u 0 , T T , C C and ϕ ˆ ϕ ˆ when  y ˆ ,
(6)

where u and v are the velocity components along the x and y direction respectively, σ is the electrical conductivity, B 0 is magnetic field flux density, ν kinematic viscosity of the base fluid, α is the thermal diffusivity of the porous medium, D B is the Brownian diffusion coefficient, D T is thermophoresis diffusion coefficient, D C T and D T C are the Soret and Dufour diffusivities, D S is the solutal diffusivity, T is the fluid temperature, C is the solutal concentration, ϕ ˆ is the nanoparticle volume fraction, ( ρ c ) f and ( ρ c ) p are the heat capacity of the fluid and the effective heat capacity of the nanoparticle material respectively, τ is a parameter defined by ( ρ c ) f / ( ρ c ) P . Using the similarity variables

η = y a ν , ψ = ( a ν ) 1 2 f ( η ) , θ ( η ) = T T T w T , S ( η ) = C C C w C , ϕ ( η ) = ϕ ˆ ϕ ˆ ϕ ˆ w ϕ ˆ ,
(7)

equations (1)-(5) reduce to the following non-similar forms where primes denote differentiation with respect to η:

f +f f f 2 M f =0,
(8)
θ +Prf θ +PrNb θ ϕ +PrNt θ 2 +Nd S =0,
(9)
S +Lef S +Ld θ =0,
(10)
ϕ +Lnf ϕ + Nt Nb θ =0,
(11)

subject to the boundary conditions

f = 0 , f = 1 , θ = 1 , S = 1 , ϕ = 1 at  η = 0 , f 0 , θ 0 , S 0 , ϕ 0 as  η .
(12)

The parameters in equations (8)-(11) are the magnetic number M, the Prandtl number Pr, the Lewis number Le, the Brownian motion parameter Nb, the thermophoresis parameter Nt, the nanofluid Lewis number Ln, the modified Dufour parameter Nd and the Dufour-solutal Lewis number Ld. These parameters are defined as

M = σ B 0 2 ρ f a , Pr = ν α , Le = α D S , Nb = τ D B ( ϕ ˆ w ϕ ˆ ) ν , Nt = τ D T ( T w T ) T ν , Ln = ν D B , Nd = D T C ( C w C ) α ( T w T ) , Ld = D C T ( T w T ) D S ( C w C ) .

The parameters of engineering interest in heat and mass transport problems are the local Nusselt number Nu x , the Sherwood number Sh x and the nanofluid Sherwood number Sh x , n . These parameters characterise the wall heat, the regular and nano mass transfer rates, respectively, and are defined by

Nu x = x T w T ( T y ) | y = 0 = Re x 1 2 θ ( 0 ) , Shr x = x C w C ( C y ) | y = 0 = Re x 1 2 S ( 0 ) , Sh x , n = x ϕ ˆ w ϕ ˆ ( ϕ ˆ y ) | y = 0 = Re x 1 2 ϕ ( 0 ) .

Following Khan and Aziz [29], the physical parameters of interest are the reduced Nusselt Nur, the Sherwood number Sh ˆ and the reduced Sherwood Shr defined as

Nur= Nu x / Re x 1 2 ,Shr= Sh x , n / Re x 1 2 andSh= Sh ˆ x / Re x 1 2 .

3 Method of solution

The system of equations (8)-(11) together with the boundary conditions (12) were solved using the successive linearisation method (SLM) (see [25, 26, 30]). The unknown functions f(η), θ(η), S(η) and ϕ(η) are expanded as

f ( η ) = f i ( η ) + m = 0 i 1 F m ( η ) , θ ( η ) = θ i ( η ) + m = 0 i 1 Θ m ( η ) , S ( η ) = S i ( η ) + m = 0 i 1 S ˜ m ( η ) , ϕ ( η ) = ϕ i ( η ) + m = 0 i 1 Φ m ( η ) , }
(13)

where f i , θ i , S i and ϕ i are unknown and F m , Θ m , S ˜ m and Φ m (m1) are successive approximations that are obtained by recursively solving the linear forms of the equation system that results from substituting (13) into equations (8)-(11). In particular, the linearised equations to be solved are

F i + a 1 , i 1 F i + a 2 , i 1 F i + a 3 , i 1 F i = r 1 , i 1 ,
(14)
Θ i + b 1 , i 1 Θ i + b 2 , i 1 F i + b 3 , i 1 S ˜ i + b 4 , i 1 Φ i = r 2 , i 1 ,
(15)
S ˜ i + c 1 , i 1 S ˜ i + c 2 , i 1 F i + c 3 , i 1 Θ i = r 3 , i 1 ,
(16)
Φ i + d 1 , i 1 Φ i + d 2 , i 1 F i + d 3 , i 1 Θ i = r 4 , i 1 ,
(17)

subject to the boundary conditions

F i (0)= F i (0)= F i ()= Θ i (0)= Θ i ()= S ˜ i (0)= S ˜ i ()= Φ i (0)= Φ i ()=0,
(18)

where coefficient parameters a k , i 1 , b k , i 1 , c k , i 1 , d k , i 1 (k=1,,4) and r j , i 1 (j=1,,4) are known constants. The initial guesses F 0 (η), Θ 0 (η), S ˜ 0 (η) and Φ 0 (η) are chosen to satisfy the boundary conditions

F 0 ( η ) = 0 , F 0 ( η ) = 1 , Θ 0 ( η ) = 1 , S ˜ 0 ( η ) = 1 , Φ 0 ( η ) = 1 at  η = 0 , F 0 ( η ) 0 , Θ 0 ( η ) 0 , S ˜ 0 ( η ) 0 , Φ 0 ( η ) 0 as  η }
(19)

and are chosen as

F 0 (η)=1 e η , Θ 0 (η)= e η , S ˜ 0 = e η , Φ 0 (η)= e η .
(20)

Starting from the initial guesses and iterating M 1 times, the functions f(η), θ(η) and ϕ(η) are written as

f ( η ) m = 0 M 1 F m ( η ) , θ ( η ) m = 0 M 1 Θ m ( η ) , S ( η ) m = 0 M 1 S ˜ m ( η ) , Φ ( η ) m = 0 M 1 Φ m ( η ) ,
(21)

where M 1 is the order of the SLM approximation. Equations (14)-(17) are solved using the Chebyshev spectral collocation method. The method is based on the Chebyshev polynomials defined on the interval [1,1]. We first transform the domain of solution [0,) into the domain [1,1] using the domain truncation technique where the problem is solved in the interval [0,L] where L is a scaling parameter used to invoke the boundary condition at infinity. This is achieved by using the mapping

η L = ξ + 1 2 ,1ξ1.
(22)

We discretise the domain [1,1] using the Gauss-Lobatto collocation points given by

ξ=cos π j N ,j=0,1,2,,N,
(23)

where N is the number of collocation points used. The functions F i , Θ i , S ˜ i and Φ i for i1 are approximated at the collocation points as follows:

F i ( ξ j ) k = 0 N F i ( ξ k ) T k ( ξ j ) , Θ i ( ξ j ) k = 0 N Θ i ( ξ k ) T k ( ξ j ) , S ˜ i ( ξ j ) k = 0 N S ˜ i ( ξ k ) T k ( ξ j ) , Φ i ( ξ j ) k = 0 N Φ i ( ξ k ) T k ( ξ j ) , } j=0,1,,N,
(24)

where T k is the k th Chebyshev polynomial given by

T k (ξ)=cos [ k cos 1 ( ξ ) ] .
(25)

The derivatives of the variables evaluated at the collocation points ξ= ξ j are represented as

d r F i d η r = k = 0 N D j k r F i ( ξ k ) , d r Θ i d η r = k = 0 N D j k r Θ i ( ξ k ) , d r Θ i d η r = k = 0 N D j k r Θ i ( ξ k ) , d r Φ i d η r = k = 0 N D j k r Φ i ( ξ k ) , } j=0,1,,N,
(26)

where r is the order of differentiation and D= 2 L D with D being the Chebyshev spectral differentiation matrix (see, for example, [3133]), whose entries are defined as

D 00 = 2 N 2 + 1 6 , D j k = c j c k ( 1 ) j + k ξ j ξ k , j k ; j , k = 0 , 1 , , N , D k k = ξ k 2 ( 1 ξ k 2 ) , k = 1 , 2 , , N 1 , D N N = 2 N 2 + 1 6 . }
(27)

Substituting equations (22)-(26) into equations (14)-(17) leads to the matrix equation

A i 1 X i = R i 1 ,
(28)

where A i 1 is a (4N+4)×(4N+4) square matrix and X i and R i 1 are (4N+4)×1 column vectors defined by

A i 1 = [ A 11 A 12 A 13 A 14 A 21 A 22 A 23 A 24 A 31 A 32 A 33 A 34 A 41 A 42 A 43 A 44 ] , X i = [ F i Θ i S ˜ Φ i ] , R i 1 = [ r 1 , i 1 r 2 , i 1 r 3 , i 1 r 4 , i 1 ] .
(29)

The functions and parameters in equation (29) are

F i = [ f i ( ξ 0 ) , f i ( ξ 1 ) , , f i ( ξ N 1 ) , f i ( ξ N ) ] T , Θ i = [ θ i ( ξ 0 ) , θ i ( ξ 1 ) , , θ i ( ξ N 1 ) , θ i ( ξ N ) ] T , S ˜ i = [ S i ( ξ 0 ) , θ i ( ξ 1 ) , , S i ( ξ N 1 ) , S i ( ξ N ) ] T , Φ i = [ ϕ i ( ξ 0 ) , ϕ i ( ξ 1 ) , , ϕ i ( ξ N 1 ) , ϕ i ( ξ N ) ] T , r 1 , i 1 = [ r 1 , i 1 ( ξ 0 ) , r 1 , i 1 ( ξ 1 ) , , r 1 , i 1 ( ξ N 1 ) , r 1 , i 1 ( ξ N ) ] T , r 2 , i 1 = [ r 2 , i 1 ( ξ 0 ) , r 2 , i 1 ( ξ 1 ) , , r 2 , i 1 ( ξ N 1 ) , r 2 , i 1 ( ξ N ) ] T , r 3 , i 1 = [ r 3 , i 1 ( ξ 0 ) , r 3 , i 1 ( ξ 1 ) , , r 3 , i 1 ( ξ N 1 ) , r 3 , i 1 ( ξ N ) ] T , r 4 , i 1 = [ r 3 , i 1 ( ξ 0 ) , r 3 , i 1 ( ξ 1 ) , , r 3 , i 1 ( ξ N 1 ) , r 3 , i 1 ( ξ N ) ] T , A 11 = D 3 + a 1 , i 1 D 2 + a 2 , i 1 D + a 3 , i 1 I , A 12 = [ 0 ] , A 13 = [ 0 ] , A 14 = [ 0 ] , A 21 = b 2 , i 1 I , A 22 = D 2 + b 1 , i 1 D , A 23 = b 3 , i 1 D 2 , A 24 = b 4 , i 1 D , A 31 = c 2 , i 1 I , A 32 = c 3 , i 1 D 2 , A 33 = D 2 + c 1 , i 1 D , A 34 = [ 0 ] , A 41 = d 2 , i 1 I , A 42 = d 3 , i 1 D 2 , A 44 = [ 0 ] , A 43 = D 2 + d 1 , i 1 D .

In the definitions above, T stands for transpose, a k , i 1 (k=1,,3), b k , i 1 (k=1,,4), c k , i 1 (k=1,,3), d k , i 1 (k=1,,3) and r k , i 1 (k=1,,4) are diagonal matrices of order (N+1)×(N+1), I is an identity matrix of order (N+1)×(N+1) and [ 0 ] is a zero matrix of order (N+1)×(N+1). The solution is obtained as

X i = A i 1 1 R i 1 .
(30)

4 Results and discussion

In this section we present solutions of equations (8)-(11) along with the boundary conditions (12) using the SLM iteration scheme. Tables 1 and 2 give a comparison between the present results and Khan and Pop [15] for the reduced Nusselt and Sherwood numbers respectively. There is a good agreement between the two sets of results with the SLM having converged at the fourth order up to eleven decimal places. The velocity components f(η) and f (η) are plotted in Figures 1(a) and 1(b) for different values of the magnetic field parameter M. As is now well known, the velocity decreases with increases in the magnetic field parameter due to an increase in the Lorentz drag force that opposes the fluid motion.

Figure 1
figure 1

Effect of the magnetic field M on the velocity components (a) f(η) and (b) f (η) .

Table 1 Comparison of results for the reduced Nusselt number θ(0) with M=0 , Pr=10 , Le=10
Table 2 Comparison of results for the reduced Sherwood number ϕ(0) with M=0 , Pr=10 , Le=10

Figures 2(a) and 2(b) show the effect of the thermophoresis parameter on the temperature and mass volume fraction profiles. The thermophoretic force generated by the temperature gradient creates a fast flow away from the stretching surface. In this way more fluid is heated away from the surface, and consequently, as Nt increases, the temperature within the boundary layer increases. The fast flow from the stretching sheet carries with it nanoparticles leading to an increase in the mass volume fraction boundary layer thickness.

Figure 2
figure 2

Effect of the thermophoresis parameter Nt on the temperature θ and nanoparticle ϕ profiles.

Figures 3(a) and 3(b) show the effect of the Lewis number Le, and the Dufour-solutal Lewis number Ld on the species concentration in the boundary layer. The concentration profiles significantly contract as the Lewis number increases. The effect of the random motion of the nanoparticles suspended in the fluid on the temperature and nanoparticle volume fraction is shown in Figures 4(a) and 4(b). As expected, the increased Brownian motion of the nanoparticles carries with it heat and the thickness of the thermal boundary layer increases. The Brownian motion of the nanoparticles increases thermal transport which is an important mechanism for the enhancement of thermal conductivity of nanofluids. However, we note that increasing the Brownian motion parameter leads to a clustering of the nanoparticles near the stretching sheet. An increase in the Brownian motion of the nanoparticles leads to a decrease in the mass volume fraction profiles.

Figure 3
figure 3

Effect of the Lewis number Le and the Dufour-solutal Lewis number Ld on concentration profiles.

Figure 4
figure 4

Effect of the Brownian motion parameter Nb on the temperature and nanoparticle volume fraction profiles.

Figures 5(a) and 5(b) show the temperature profiles for several values of the Prandtl number Pr and mass volume fraction profile for several values of the modified Dufour number Nd. The temperature profiles decrease as the Prandtl number increases since, for high Prandtl numbers, the flow is governed by momentum and viscous diffusion rather than thermal diffusion. On the other hand, the thickness of the mass volume fraction boundary layer increases with an increase in Nd.

Figure 5
figure 5

Effect of Pr and Nd on the temperature θ , profiles respectively.

Figures 6(a) and 6(b) show the effects of the thermophoresis parameter Nt, the Lewis number Le, the magnetic field parameter M, the Prandtl number Pr and the modified Dufour number Nd on the wall heat and mass fraction transfer rates. It can be seen that the thermal boundary layer thickness increases when the thermophoresis parameter Nt increases, thus decreasing the reduced Nusselt number. However, increasing the Lewis number Le leads to a decrease in the reduced Nusselt number. On the other hand, the results show that the reduced Nusselt number increases with increasing Prandtl numbers. Increasing both the magnetic field parameter M and the modified Dufour parameter Nd leads to an increase in the thermal boundary layer thickness, thus reducing the Nusselt number.

Figure 6
figure 6

Effect of Nt , Le , Pr , Nd and M on the heat transfer coefficient Nur .

Figures 7(a) and 7(b) show the effects of the Dufour-solutal Lewis number Ld and the nanofluid Lewis number Ln on the reduced Nusselt number Nur as the Brownian motion parameter Nb increases. We note a decrease in the reduced Nusselt number when Ln increases, and an increase in the reduced Nusselt number when Ld increases.

Figure 7
figure 7

Effect of the Dufour-solutal Lewis number Ld and the nanofluid Lewis number Ln on the reduced Nusselt number Nur .

Figures 8(a) and 8(b) show the graphs of θ (0) and S (0) plotted against the Dufour-solutal Lewis number Ld for different values of the parameters Nt, Nb and Le. We observe that θ (0) increases in the absence of the Brownian motion and the thermophoresis parameter while θ (0) decreases in the presence of Brownian motion and thermophoresis parameters. An increase in S (0) is observed in the presence of both the Brownian motion and the thermophoresis parameter. Figures 9(a) and 9(b) show the effect of increasing Nt and Nb respectively on the reduced Sherwood number ϕ(0).

Figure 8
figure 8

Effect of the Lewis number Le , the thermophoresis parameter Nt and the Brownian motion parameter Nb on (a) the reduced Nusselt number θ (0) and (b) the local Sherwood number S (0) .

Figure 9
figure 9

Effect of the nanofluid Lewis number Ln , the thermophoresis parameter Nt and the Brownian motion parameter Nb on the nanofluid Sherwood number ϕ (0) .

5 Conclusions

A numerical study of the magneto-nanofluid boundary layer flow over a stretching sheet was carried out. We determined the effects of various parameters on the fluid properties as well as on the heat, and the regular and nano mass transfer rates. We have shown that increasing the magnetic field parameter M tends to retard the fluid flow within the boundary layer. The effects of the Prandtl number, the Lewis number, the Brownian motion parameter, the thermophoresis parameter, the nanofluid Lewis number, the modified Dufour parameter and the Dufour-solutal Lewis number on the heat, regular and nano mass transfer coefficients and fluid flow characteristics have been studied. We have shown inter alia that:

– the thermal boundary layer thickness increases with the thermophoresis parameter;

– increasing the Lewis number reduces the heat transfer coefficient;

– the heat transfer coefficient increases in the absence of the Brownian motion and the thermophoresis parameter and decreases in the presence of Brownian motion and thermophoresis parameters.

References

  1. Choi SUS: Enhancing thermal conductivity of fluid with nanoparticles.23. In Developments and Applications of Non-Newtonian Flows. FED, New York; 1995:99–105.

    Google Scholar 

  2. Kakac S, Pramuanjaroenkij A: Review of convective heat transfer enhancement with nanofluid. Int. J. Heat Mass Transf. 2009, 52: 3187–3196. 10.1016/j.ijheatmasstransfer.2009.02.006

    Article  Google Scholar 

  3. Choi SUS, Zhang ZG, Yu W, Lockwood FE, Grulke EA: Anomalously thermal conductivity enhancement in nanotube suspensions. Appl. Phys. Lett. 2001, 79: 2252–2254. 10.1063/1.1408272

    Article  Google Scholar 

  4. Masuda H, Ebata A, Teramae K, Hishinuma N: Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles. Netsu Bussei 1993, 7: 227–233. 10.2963/jjtp.7.227

    Article  Google Scholar 

  5. Eapen J, Rusconi R, Piazza R, Yip S: The classical nature of thermal conduction in nanofluids. J. Heat Transf. 2010., 132: Article ID 102402

    Google Scholar 

  6. Fan J, Wang L: Effective thermal conductivity of nanofluids: the effects of microstructure. J. Phys. D, Appl. Phys. 2010., 43: Article ID 165501

    Google Scholar 

  7. Nield DA, Kuznetsov AV: The Cheng-Minkowycz problem for natural convective boundary-layer flow over a porous medium saturated by a nanofluid. Int. J. Heat Mass Transf. 2010, 52: 5792–5795.

    Article  Google Scholar 

  8. Nield DA, Kuznetsov AV: Thermal instability in a porous medium layer saturated by a nanofluid: Brinkman model. Transp. Porous Media 2010, 81: 409–422. 10.1007/s11242-009-9413-2

    Article  MathSciNet  Google Scholar 

  9. Aziz A, Khan WA, Pop I: Free convection boundary layer flow past a horizontal flat plate embedded in porous medium filled by nanofluid containing gyrotactic microorganisms. Int. J. Therm. Sci. 2012, 56: 48–57.

    Article  Google Scholar 

  10. Cheng C-Y: Free convection boundary layer flow over a horizontal cylinder of elliptic cross section in porus media saturated by a nanofluid. Int. Commun. Heat Mass Transf. 2012, 39: 931–936. 10.1016/j.icheatmasstransfer.2012.05.014

    Article  Google Scholar 

  11. Chamkha A, Gorla RSR, Ghodeswar K: Non-similar solution for natural convective boundary layer flow over a sphere embedded in a porous medium saturated with a nanofluid. Transp. Porous Media 2011, 86: 13–22. 10.1007/s11242-010-9601-0

    Article  MathSciNet  Google Scholar 

  12. Altan T, Oh S, Gegel H: Metal Forming Fundamentals and Applications. Am. Soc. Metals, Metals Park; 1979.

    Google Scholar 

  13. Fisher EG: Extrusion of Plastics. Wiley, New York; 1976.

    Google Scholar 

  14. Tidmore Z, Klein I Polymer Science and Engineering Series. In Engineering Principles of Plasticating Extrusion. Van Norstrand, New York; 1970.

    Google Scholar 

  15. Khan WA, Pop I: Boundary layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53: 2477–2483. 10.1016/j.ijheatmasstransfer.2010.01.032

    Article  Google Scholar 

  16. 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 

  17. Narayana M, Sibanda P: Laminar flow of a nanoliquid film over an unsteady stretching sheet. Int. J. Heat Mass Transf. 2012, 55: 7552–7560. 10.1016/j.ijheatmasstransfer.2012.07.054

    Article  Google Scholar 

  18. Kameswaran PK, Narayana N, Sibanda P, Murthy PVSN: Hydromagnetic nanofluid flow due to a stretching or shrinking sheet with viscous dissipation and chemical reaction effects. Int. J. Heat Mass Transf. 2012, 55: 7587–7595. 10.1016/j.ijheatmasstransfer.2012.07.065

    Article  Google Scholar 

  19. Shima PD, Philip J, Raj B: Magnetically controllable nanofluid with tunable thermal conductivity and viscosity. Appl. Phys. Lett. 2009., 95: Article ID 133112

    Google Scholar 

  20. Ganguly R, Sen S, Puri IK: Heat transfer augmentation using a magnetic fluid under the influence of a line dipole. J. Magn. Magn. Mater. 2004, 271: 63–73. 10.1016/j.jmmm.2003.09.015

    Article  Google Scholar 

  21. 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. 2012, 55: 2102–2109. 10.1016/j.ijheatmasstransfer.2011.12.013

    Article  Google Scholar 

  22. Hamad MAA, Ferdows M: Similarity solutions to viscous flow and heat transfer of nanofluid over nonlinearly stretching sheet. Appl. Math. Mech. 2012, 33: 923–930.

    Article  MathSciNet  Google Scholar 

  23. Makukula ZG, Motsa SS, Sibanda P: On a new solution for the viscoelastic squeezing flow between two parallel plates. J. Adv. Res. Appl. Math. 2010, 2: 31–38.

    Article  MathSciNet  Google Scholar 

  24. Awad FG, Sibanda P, Motsa SS, Makinde OD: Convection from an inverted cone in a porous medium with cross-diffusion effects. Comput. Math. Appl. 2011, 61: 1431–1441. 10.1016/j.camwa.2011.01.015

    Article  MathSciNet  Google Scholar 

  25. Makukula ZG, Sibanda P, Motsa SS: A novel numerical technique for two-dimensional laminar flow between two moving porous walls. Math. Probl. Eng. 2010., 2010: Article ID 528956. doi:10.1155/2010/528956

    Google Scholar 

  26. Makukula ZG, Motsa SS, Sibanda P: A novel numerical technique for two-dimensional laminar flow between two moving porous walls. Math. Probl. Eng. 2010., 2010: Article ID 528956. doi:10.1155/2010/528956

    Google Scholar 

  27. Awad FG, Sibanda P, Narayana M, Motsa SS: Convection from a semi-finite plate in a fluid saturated porous medium with cross-diffusion and radiative heat transfer. Int. J. Phys. Sci. 2011, 6: 4910–4923.

    Google Scholar 

  28. Motsa SS, Sibanda P, Shateyi S: On a new quasi-linearization method for systems of nonlinear boundary value problems. Math. Methods Appl. Sci. 2011, 34: 1406–1413. 10.1002/mma.1449

    Article  MathSciNet  Google Scholar 

  29. Khan WA, Aziz A: Double-diffusive natural convective boundary layer flow in a porous medium saturated with a nanofluid over a vertical plate: prescribed surface heat, solute and nanoparticle fluxes. Int. J. Therm. Sci. 2011, 50: 2154–2160. 10.1016/j.ijthermalsci.2011.05.022

    Article  Google Scholar 

  30. Shateyi S, Motsa SS: Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with hall effect. Bound. Value Probl. 2010., 2010: Article ID 257568. doi:10.1155/2010/257568

    Google Scholar 

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

    Book  Google Scholar 

  32. Don WS, Solomonoff A: Accuracy and speed in computing the Chebyshev collocation derivative. SIAM J. Sci. Comput. 1995, 16: 1253–1268. 10.1137/0916073

    Article  MathSciNet  Google Scholar 

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

    Book  Google Scholar 

Download references

Acknowledgements

The authors wish to thank the University of KwaZulu-Natal for financial support.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Precious Sibanda.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The work including proofreading was done by all the authors.

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

Awad, F.G., Sibanda, P. & Khidir, A.A. Thermodiffusion effects on magneto-nanofluid flow over a stretching sheet. Bound Value Probl 2013, 136 (2013). https://doi.org/10.1186/1687-2770-2013-136

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-136

Keywords