This study theoretically investigates the conjugate effects of radiation flux and magnetohydrodynamic (MHD) on free convection boundary layer flow of a nanofluid over a nonlinear stretching sheet. It is assumed that the magnetic Reynolds number is small enough and the sheet is stretched with a power law velocity under the effects of the magnetic field, the buoyancy parameter, and the solutal buoyancy parameter. The model used for the nanofluid incorporates the effects of Rosseland approximation, Brownian motion, and thermophoresis parameters. By using appropriate similarity transformations, the governing nonlinear partial differential equations are transformed into dimensionless form and numerically solved using an implicit finite difference scheme known as the Keller-box method. It is found that the variations of magnetic field, buoyancy parameter, solutal buoyancy parameter, and the power law velocity parameter have strong influence on the motion.
1. Introduction
Heat and mass transfer (double-diffusion) phenomenon on free convection is driven by two density gradients which have different rates of diffusion and currently is an important fluid dynamics topic. A common example of double diffusive convection appears in oceanography, where heat and salt concentrations exist with different gradients and diffuse at differing rates. Double diffusive convection is also important in understanding the evolution of a number of systems that have multiple causes for density variations. These include convection in the earth’s oceans, in magma chambers, and in the sun where heat and helium diffuse at differing rates [1]. Double diffusive convection flows for Newtonian and non-Newtonian fluids are extensively studied. However, for nanofluids such studies are scarce due to their complicated nature [2, 3]. Moreover, recent developments in the field of fluid dynamics and nanotechnology confirm that nanofluids are industrially more important than other available fluids [4]. There are numerous biomedical applications that involve nanofluids such as magnetic cell separation, drug delivery, hyperthermia, and contrast enhancement in magnetic resonance imaging [5]. Hence, this is the motivation for considering nanofluids in the present work.
On the other hand, stretching sheet problems with double diffusion are important in extrusion process, glass fibber, paper production, hot rolling, wire drawing, electronic chips, crystal growing, plastic manufactures, food processing, and movement of biological fluids [6]. Khan and Pop [7] investigated the laminar flow of a nanofluid on a stretching flat surface by incorporating the effects of Brownian motion, thermophoresis and reported to be the pioneer for this study of stretching sheet in nanofluid. Rana and Bhargava [8] discussed the flow and heat transfer of a nanofluid over a nonlinearly stretching sheet. Furthermore, the MHD flow of nanofluid over a power law stretching sheet plays an important role in various industrial applications including magnetic control of molten iron flow in the steel industry and liquid metal cooling in nuclear reactors [9]. In addition, when thermal radiation is considered, such studies have useful chemical processing applications [10, 11].
The aim of the present work is to numerically investigate the thermal radiation effects for double diffusive MHD free convection flow of a nanofluid over a power law stretching sheet using Keller-box method [12, 13].
2. Problem Formulation
The unsteady two-dimensional boundary layer flow of a nanofluid past a stretching surface is considered. The stretching velocity is assumed to be uw(x)=axm, where a is constant, m(m≥0) is the power law exponent parameter, and x is the coordinate measured along the stretching surface. The flow takes place above the stretching surface at y≥0. Here y is the coordinate axis measured normal to the stretching surface. A uniform stress leading to equal and opposite forces is applied along the x-axis, so that the sheet is stretched keeping the origin fixed. Following the work of Hayat and Qasim [6] and Afify [14], we assume that a uniform transverse magnetic field of strength B is imposed parallel to the y-axis and the induced magnetic field due to the motion of the electrically-conducting fluid is negligible. We also assume that the external electrical field is zero and the electric field due to the polarization of charges is negligible. Moreover, it is further supposed that at the stretching surface, the temperature T and the nanoparticle fraction C take constant values Tw and Cw, whereas the ambient values of temperature T∞ and the nanoparticle fraction C∞are attained as y tends to infinity. The Oberbeck-Boussinesq approximation is employed to the field equations. The governing boundary layer equations that are based on the balance laws of mass, linear momentum, energy, and concentration species for the present problem are given as follows [3, 15]:
(1)∇·V=0,ρf(∂V∂t+(V·∇)V)=-∇p+μ∇2V-σB2(x)V+[(1-C∞)ρf∞βT(T-T∞)+(ρp-ρf∞)βC(C-C∞)]g,(ρc)f(∂T∂t+V·∇T)=k∇2T+(ρc)p[DB∇C·∇T+(DTT∞)∇T·∇T]-∇·qr,∂C∂t+V·∇C=DB∇C+(DTT∞)∇2T,
where V is the velocity vector, σ is the electrical conductivity, gis the acceleration due to gravity, μ is the viscosity, ρf is the density of the base fluid, ρp is the density of the nanoparticle, βT is the coefficient of volumetric thermal expansion, βCis the coefficient of volumetric concentration expansion, DB is the Brownian diffusion coefficient, DT is the thermophoretic diffusion coefficient, kis the thermal conductivity, (ρc)pis the heat capacitance of the nanoparticles, (ρc)f is the heat capacitance of the base fluid, and qr is the radiation flux.
In the following we adopt the Rosseland approximation
(2)qr=-4σ*3k*∂T4∂y,
where σ* is the Stefan-Boltzmann constant and k* is the mean absorption coefficient. Assuming that the temperature difference between the free stream T∞ and the fluid temperature T is small enough, expanding T4 in a Taylor series about T∞, and neglecting higher order terms, it results that
(3)T4=4T∞3T-3T∞4.
The two-dimensional steady-state problem under the assumptions that the external pressure in x-direction having diluted nanoparticles is constant; that is, p=p0 yields to ∂p/∂x=0 and Oberbeck-Boussinesq approximation is valid; (1) are reduced to
(4)∂u∂x+∂v∂y=0,0=μ∂2u∂y2-ρf(u∂u∂x+v∂u∂y)-σB2(x)u+[(1-C∞)ρf∞gβT(T-T∞)m-(ρp-ρf∞)gβC(C-C∞)],u∂T∂x+v∂T∂y=α∇2T+τ[DB∂C∂y∂T∂y+DTT∞(∂T∂y)2]-1(ρc)f∂qr∂y,u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2,
where
(5)α=k(ρc)f,τ=(ρc)p(ρc)f.
Here u and v are the velocity components in the x and y-directions, respectively, α is the thermal diffusivity parameter and τ is the ratio between the effective heat capacity of the nanoparticle material and heat capacity of the fluid.
The associated boundary conditions are
(6)u=uw(x)=axm,v=0,T=Tw,C=Cwaty=0,u→0,v→0,T→T∞,C→C∞asy→∞.
Using the stream function ψ=ψ(x,y), the components of velocity are defined as
(7)u=∂ψ∂y,v=-∂ψ∂x.
By using the following similarity transformations [14],
(8)ψ=2νaxm+1m+1f(η),θ(η)=T-T∞Tw-T∞,ϕ(η)=C-C∞Cw-C∞,η=y(m+1)axm-12ν,
into (4), we get
(9)f′′′+ff′′-2mm+1f′2+2m+1(M2f′+λθ-δϕ)=0,1Pr(1+N)θ′′+fθ′+Nbϕ′θ′+Ntθ′2=0,ϕ′′+Lefϕ′+NtNbθ′′=0,
where
(10)M2=σB02aρf,λ=GrRex3/2,δ=GmRex3/2,Pr=να,N=16σ*T∞33kk*,Le=νDB,ν=μρf,Nb=τDB(Cw-C∞)ν,Nt=τDT(Tw-T∞)νT∞,Rex=uw(x)xν,Gr=(1-C∞)(ρf∞/ρf)gn(Tw-T∞)ν2Rex1/2,Gm=((ρp-ρf∞)/ρf)gn1(Cw-C∞)ν2Rex1/2.
Here Mis the magnetic parameter called Hartmann number, λ is the buoyancy parameter, δ is the solutal buoyancy parameter, Pr is the Prandtl number, N is the radiation parameter, Le is the Lewis number, ν is the kinematic viscosity of the fluid, Nb is the Brownian motion parameter, Nt is the thermophoresis parameter, Rexis the local Reynolds number based on the stretching velocity, Gr is the local thermal Grashof number, Gm is the local concentration Grashof number, and f, θ, and ϕ are the dimensionless stream functions, temperature, and rescaled nanoparticle volume fraction, respectively. Further, in order to eliminate the dependence of M on x, the magnetic field strength has been taken proportional to x, that is, B(x)=B0x(m-1)/2 where B0 is a constant [14]. Also, βT and βC are proportional to x-3, that is, βT=nx-3, βC=n1x-3, where n and n1are constants of proportionality [16].
The corresponding boundary conditions are transformed to
(11)f=0,f′=1,θ=1,ϕ=1atη=0,f′→0,θ→0,ϕ→0asη→∞.
In the present case, we define Nusselt number, Sherwood number, and skin-friction as follows:
(12)Nu=qwxk(Tw-T∞),Sh=qmxDB(Cw-C∞),Cf=τw(1/2)ρU2,
where
(13)qw=-k∂T∂y,qm=-DB∂C∂y,τw=μ∂u∂y,aty=0.
The associated expressions of the dimensionless reduced Nusselt number -θ′(0), reduced Sherwood number -ϕ′(0), and skin-friction coefficientCfx are defined as
(14)Nur=Nu((m+1)/2)Rex,Shr=Sh((m+1)/2)Rex,Cfx=Cf22m+1Rex.
The transformed nonlinear ordinary differential (9) with the boundary conditions (11) are numerically solved by means of the Keller-box method [12, 13].
3. Results and Discussion
Numerical results for some physical parameters of interest are shown in Tables 1 and 2 and in Figures 1–9. In Table 1, the results for reduced Nusselt number -θ′(0) and reduced Sherwood number -ϕ′(0), when the Hartmann number M, the buoyancy parameter λ, the solutal buoyancy parameter δ, and the radiation parameter N are zero and the power law velocity parameter m=1, are compared with those obtained by Khan and Pop [7]. We found that, this comparison shows an excellent agreement for the involved flow parameters and in general their solutions appeared as the special cases of our general solutions. Similarly if we take Nb=0.5,Nt=0.5,Pr=2.0,Le=2.0, and m=2.0, the results of Rana and Bhargava [8] can also be obtained as a special case. Hence we are quite confident that our results are accurate and are more general than the existing studies in the literature.
Comparison of reduced Nusselt number and reduced Sherwood number.
Nb
Nt
Pr
Le
M
λ
δ
N
m
Khan and Pop [7]
Present result
-θ′(0)
-ϕ′(0)
-θ′(0)
-ϕ′(0)
0.1
0.1
10
10
0.0
0.0
0.0
0.0
1.0
0.9524
2.1294
0.9524
2.1294
0.2
0.2
10
10
0.0
0.0
0.0
0.0
1.0
0.3654
2.5152
0.3654
2.5152
0.3
0.3
10
10
0.0
0.0
0.0
0.0
1.0
0.1355
2.6088
0.1355
2.6088
0.4
0.4
10
10
0.0
0.0
0.0
0.0
1.0
0.0495
2.6038
0.0495
2.6038
0.5
0.5
10
10
0.0
0.0
0.0
0.0
1.0
0.0179
2.5731
0.0179
2.5731
Variations of reduced Nusselt number, reduced Sherwood number, and skin-friction coefficient.
Nb
Nt
Pr
Le
M
λ
δ
N
m
-θ′(0)
-ϕ′(0)
Cfx
0.1
0.1
0.71
10
0.1
1.0
1.0
1.0
0.5
0.3399
2.3369
0.6148
0.5
0.1
0.71
10
0.1
1.0
1.0
1.0
0.5
0.3004
2.3976
0.6021
0.1
0.5
0.71
10
0.1
1.0
1.0
1.0
0.5
0.3131
2.1134
0.6473
0.1
0.1
7.00
10
0.1
1.0
1.0
1.0
0.5
0.8824
2.0780
0.8870
0.1
0.1
0.71
25
0.1
1.0
1.0
1.0
0.5
0.3418
3.8356
0.5647
0.1
0.1
0.71
10
2.0
1.0
1.0
1.0
0.5
0.2393
2.1222
1.6688
0.1
0.1
0.71
10
0.1
3.0
1.0
1.0
0.5
0.4042
2.4740
−0.1679
0.1
0.1
0.71
10
0.1
1.0
4.0
1.0
0.5
0.3265
2.2761
1.0759
0.1
0.1
0.71
10
0.1
1.0
1.0
6.0
0.5
0.1817
2.4120
0.4856
0.1
0.1
0.71
10
0.1
1.0
1.0
1.0
5.0
0.2804
2.2294
1.1273
Variations of velocity profiles with η for different values of Pr, M, λ, δ, N, and m.
Variations of temperature profiles with η for different values of Pr, M, λ, δ, N, and m.
Variations of concentration profiles with η for different values of Pr, M, λ, δ, N, and m.
Variations of reduced Nusselt number with λ for different values of Pr, M, λ, δ, N, and m.
Variations of reduced Sherwood number with λ for different values of Pr, M, λ, δ, N, and m.
Variations of skin-friction coefficient with λ for different values of Pr, M, λ, δ, N, and m.
Variations of reduced Nusselt number with δ for different values of Pr, M, λ, δ, N, and m.
Variations of reduced Sherwood number with δ for different values of Pr, M, λ, δ, N, and m.
Variations of skin-friction coefficient with δ for different values of Pr, M, λ, δ, N, and m.
The variations of reduced Nusselt number -θ′(0), reduced Sherwood number -ϕ′(0), and the skin-friction coefficient Cfxfor different values of Nb, Nt, Pr, Le, M, λ, δ, N, and m are shown in Table 2. It is observed that -θ′(0) is a decreasing function of Nb, Nt, M, δ, N, and m, whereas -θ′(0) increases with increasing values of Pr, Le, and λ. However, it is found that -ϕ′(0) decreases for large values of Nt, Pr, M, δ, and m and increases for increasing values of Nb, Le, λ, and N. It is further observed from this table that Cfx is an increasing function of Nt, Pr, M, δ, and m and a decreasing one with respect to Nb, Le, λ, and N.
In order to get physical insight into the problem, convergence of plots for dimensionless velocity, temperature and concentration profiles for different flow parameters are provided in Figures 1–3. Furthermore, the variations of the reduced Nusselt number, the reduced Sherwood number, and skin-friction coefficient as a function of λand δ are shown in Figures 4–9. All these graphs are plotted for fixed values of Nb=Nt=0.1and Le=10.
It is evident from Figure 1 that an increase in Prresults in a decrease in the velocity, whereas velocity increases for the increasing values of λand N. The effect of increasing values of M results in a decrease of momentum boundary layer thickness which reduces the dimensionless stream function and the velocity of nanofluid. Physically, it is possible due to the fact that an increase in Pr makes the fluid to be more viscous, which leads to decrease in the velocity for increasing values of Hartmann number and increase the frictional force which reduce the fluid velocity to slow down the fluid motion. Further, it is found that velocity decreases with increasing values of δ and m.
Figure 2 is plotted to analyze the heat transfer for combined effects of Pr, M, λ, δ, N, and m on the temperature profiles. This figure indicates that an increase in Pr and λ results in a decrease in the temperature profiles, whereas these profiles increase for the increasing values of δ, N, and m. The radiation parameter N is very effective for the thickening of the thermal boundary layer which releases the heat energy from the flow region and causes the cooling of the system. Physically, it is true due to the fact that temperature increases by increasing the Rosseland approximation. Increasing values of Mcause the temperature to rise up. This is due to the fact that nanoparticles dissipate energy in the form of heat which causes the thermal boundary layer thickness to increase in the case of nanofluids and ultimately a localized rise in temperature of the fluid occurs.
Figure 3 is prepared to study the effects of incorporated different flow parameters on concentration profiles. It results that the increasing values of λ and N cause a decrease in the mass transfer. However, the concentration increases for large values of Pr, M, δ, and m. Physically, the variations of Pr, M, λ, δ, and m are important in the boundary layer flow and play a significant role in shortening the concentration boundary layer for the mass fraction. All profiles discussed above descend smoothly in the free stream satisfying boundary conditions. This ensures the accuracy of the obtained numerical results.
The variations of the reduced Nusselt number, reduced Sherwood number, and skin-fiction coefficient as functions of λ for the different values of Pr, M, δ, N, and m are shown in Figures 4–6. From Figure 4, we observed that the reduced Nusselt number increases for increasing values of Pr and decreases for increasing values of M, δ, N, and m. This increase in the heat transfer rate at the surface of the sheet is due to the reason that nanofluids with high viscosities have large values of Prandtl number which results for relatively low thermal conductivities to reduce the conduction phenomenon [7]. Physically, this behaviour is meaningful due to the fact that with increasing viscosity of the base fluid, the thermal boundary layer thickness decreases because of less colloidal forces between the base fluid and suspended nanoparticles. So the thermal conductivities of nanofluids play an important role to increase the heat transfer coefficient between the heat transfer medium and heat transfer surface. Further, Figure 5 shows that the reduced Sherwood number decreases for increasing values of Pr but increases for increasing values of N, M, δ, and m. It is clearly seen from Figure 6 that the dimensionless skin-friction coefficient along λ increases with increasing values of Pr, M, δ, and mbut decreases with increasing values of N.
Finally, Figures 7–9 show the variation in the reduced Nusselt number, the reduced Sherwood number, and skin-fiction coefficient along δ for different values of Pr, M, λ, N, and m. It is found from Figure 7 that the reduced Nusselt number decreases for increasing values of M, N, and m and increases for increasing values of Pr and λ. The reduced Sherwood number increases for large values of λ and N but decreases when the values of Pr, M, and m are increased as shown in Figure 8. Further, Figure 9 shows that skin-friction coefficient increases for the larger values of Pr, M, and mand decreases with regard to λand N. It is clearly seen that the influence of these parameters on the local Sherwood number and skin-fiction coefficient with λ and δ is quite opposite to each other.
4. Conclsions
In the present study we have investigated the conjugate effects of thermal radiation flux on double diffusive MHD free convection boundary layer flow of a nanofluid over a nonlinear stretching sheet. The model used here for nanofluid incorporates the effects of thermophoresis parameter Nb, Brownian motion parameter Nt, Prandtl number Pr, Lewis number Le, Hartmann number M, buoyancy parameter λ, solutal buoyancy parameter δ, radiation parameter N, and nonlinear stretching parameter m. The governing nonlinear partial differential equations are numerically solved using implicit finite difference scheme known as Keller-box method. Numerical solutions depending upon all parameters are presented in tables and figures. The physical interpretations provided by the authors reflect some important characteristics of free convection flow of a nanofluid over a power law stretching sheet and play a significant role in the concentration boundary layer for the mass fraction. The following main results are concluded from this study.
The reduced Nusselt number -θ′(0) is a decreasing function of Nb, Nt, M, δ, N, and mwhereas increases with increasing values of Pr, Le, and λ.
The reduced Sherwood number -ϕ′(0) decreases for large values ofNt, Pr, M, δ, and mand increases for increasing values of Nb, Le, λ, and N.
The skin-friction coefficient Cfx is an increasing function of Nt, Pr, M, δ, and mand a decreasing one with respect to Nb, Le, λ, and N.
The velocity of nanofluid decreases for increasing Pr whereas increases for the increasing values of λ and N.
The temperature of nanofluid decreases when Pr and λ are increased whereas increases for large values of M, δ, and m.
The species concentration of nanofluid decreases for large values of λand N whereas increases for increasing values of Pr, M, δ, and m.
The reduced Nusselt number versus λ increases for Pr and decreases for M, δ, N, and m whereas the reduced Nusselt number against δ decreases for M, N, and m while increases for Pr and λ. However, the reduced Sherwood number and skin-friction coefficient versus λ and δ are found to be opposite in nature.
Acknowledgments
The authors would like to acknowledge the financial support received from MOHE, Research Management Centre-RMC (4F109), UTM, Malaysia, and University Malaysia Pahang (RDU110108).
KhanW. A.AzizA.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 fluxes201150215421602-s2.0-7995973563410.1016/j.ijthermalsci.2011.05.022NieldD. A.KuznetsovA. V.Thermal instability in a porous medium layer saturated by a nanofluid20095225-26579658012-s2.0-7584914159410.1016/j.ijheatmasstransfer.2009.07.023KuznetsovA. V.NieldD. A.The onset of double-diffusive nanofluid convection in a layer of a saturated porous medium20108539419512-s2.0-7814950125410.1007/s11242-010-9600-1ChoiS. U. S.SingerD. A.WangH. P.Enhancing thermal conductivity of fluids with nanoparticles1995231New York, NY, USAAmerican Society of Mechanical Engineers99105ChoiS. U. S.ZhangZ. G.YuW.LockwoodF. E.GrulkeE. A.Anomalous thermal conductivity enhancement in nanotube suspensions20017914225222542-s2.0-003547352910.1063/1.1408272HayatT.QasimM.Radiation and magnetic field effects on the unsteady mixed convection flow of a second grade fluid over a vertical stretching sheet20116678208322-s2.0-7844928947410.1002/fld.2285KhanW. A.PopI.Boundary-layer flow of a nanofluid past a stretching sheet20105311-12247724832-s2.0-7764934161010.1016/j.ijheatmasstransfer.2010.01.032RanaP.BhargavaR.Flow and heat transfer of a nanofluid over a nonlinearly stretching sheet: a numerical study20121712122262-s2.0-7996103006910.1016/j.cnsns.2011.05.009HayatT.QasimM.MesloubS.MHD flow and heat transfer over permeable stretching sheet with slip conditions20116689639752-s2.0-7955166846810.1002/fld.2294MagyariE.PantokratorasA.Note on the effect of thermal radiation in the linearized Rosseland approximation on the heat transfer characteristics of various boundary layer flows20113855545562-s2.0-7995463231410.1016/j.icheatmasstransfer.2011.03.006HamidR. A.ArifinN. M.NazarR. M.PopI.Radiation effects on Marangoni boundary layer flow past a flat plate in nanofluidProceedings of the International MultiConference of Engineers and Computer Scientists 2011 (IMECS '11)March 2011Hong Kong126012632-s2.0-79960584032CebeciT.BradshawP.1977New York, NY, USAHemisphere Publishing CorporationCebeciT.BradshawP.1988New York, NY, USASpringerAfifyA. A.Similarity solution in MHD: effects of thermal diffusion and diffusion thermo on free convective heat and mass transfer over a stretching surface considering suction or injection2009145220222142-s2.0-5604912479510.1016/j.cnsns.2008.07.001KaharR. A.KandasamyR.MuhaiminScaling group transformation for boundary-layer flow of a nanofluid past a porous vertical stretching surface in the presence of chemical reaction with heat radiation2011521521MakindeO. D.OlanrewajuP. O.Buoyancy effects on thermal boundary layer over a vertical plate with a convective surface boundary condition20101324044502