Unsteady MHD Mixed Convection Flow with Slip of a Nanofluid in the Stagnation Region of an Impulsively Rotating Sphere with Effects of Thermal Radiation and Convective Boundary Conditions ()
1. Introduction
Nanofluids are new types of fluids consist of a base fluid and suspended nanoparticles (1 - 100 nm). These fluids have a higher thermal conductivity and a single-phase heat transfer coefficient compared with their base fluids [1] . Buongiorno [2] made a test for seven slip mechanisms to determine which of them can produce a relative velocity between the base fluid and nanoparticles. These effects are gravity, fluid drainage, Magnus effect, diffusiophoresis, inertia Brownian diffusion and thermophoresis. His results indicated that only Brownian diffusion and thermophoresis are important in nanofluids. Taken into account this observation, he suggests a model for nanofluids consist of a two-component four-equations nonhomogeneous model for momentum, mass and heat transport. Also, he attributed this abnormal improvement in heat transfer coefficient to the significantly variations of nanofluid properties within the boundary layer due to the effects of thermophoresis and temperature gradient. This mentioned effects can produce a reduction in the fluid viscosity and consequently the heat transfer increases. Kuznetsov and Nield [3] applied the Buongiorno’s model [2] to discuss the Pohlhausen, Kuiken and Bejan [4] [5] [6] [7] problems. In this investigation, they assumed that the boundary conditions for the nanoparticle volume fraction are similar to the boundary conditions for the temperature. Kuznetsov and Nield [8] revised their problem [3] to be the nanoparticle volume fraction on the boundary which is passively rather than actively controlled. Makinde and Aziz [9] examined the influence of convective boundary conditions on the convective boundary layer flow of a nanofluid past a stretching sheet. They found that the local concentration of nanoparticles is enhanced as the Biot number increases but decreases as the Lewis number increases. Vajraveu et al. [10] analyzed the convective heat transfer and fluid flow over a stretching surface using nanofluids. It is found that the presence of nanoparticles leads to an improvement in the skin friction coefficient. More recently, Ahmed et al. [11] discussed the entropy generation analysis due to convective flow of nanofluids over a stretching cylinder. They observed that the increase in suction parameter causes an increase in the entropy generation near the surface. Aly and Ahmed [12] performed an investigation for double diffusive boundary layer flow of a nanofluid over a vertical cylinder. It is found that the local skin friction and the local Sherwood number are enhanced as the buoyancy ratio increases but the local Nusselt number takes the inverse behavior. The convective boundary layer flow of nanofluids over a flat plate under the influence of Newtonian heating and viscous dissipation was presented by Makinde [13] . His results indicated that the increase in nanoparticle volume fraction and Biot number enhances the heat transfer rate. Ahmed and Mahdy [14] studied the magnetohydrodynamic bio-convective flow of nanofluids over a vertical wavy surface. They observed that the increase in either Grashof number or magnetic field parameter results in a reduction in the local Nusselt and Sherwood numbers.
On the other hand, the study of convective boundary layer flow over rotating bodies of revolution has several engineering applications such as fiber coating, design of rotating machinery, re-entry missile, projectile motion etc. [15] . Because of this importance, many researchers were interest with this subject [16] - [22] . Chamkha and Ahmed [16] reported on MHD double diffusive mixed convection flow in a forward stagnation region of a rotating sphere under the effect of uniform chemical reaction and heat generation. They observed that regardless the cases of CWT and CHMF, the velocity component in the x-direction is enhanced as the acceleration parameter increases. Chamkha and Ahmed [17] studied the unsteady MHD double diffusive boundary layer flow near the stagnation point over a 3D porous body. They observed that the increase in either transpiration parameter or ratio of velocity gradient parameter improves the local Nusselt number. In Chamkha and Ahmed [18] , the mixed convection boundary layer flow in the forward stagnation region of a rotating sphere was investigated. The results indicated that the increase in chemical reaction parameter causes an enhancement in the local shear stress and Sherwood number. Faltas and Saad [19] examined the slip effect on Stokes flow between eccentric rotating sphere. They concluded that the effect of slip on the couple is valuable and should not be ignored. Chamkha and Ahmed [20] used the similarity solutions to discuss the unsteady MHD boundary layer flow in a stagnation region of a 3D porous body in the presence of chemical reaction and heat generation/absorption effects. They noted that the increase in either the heat generation/absorption parameter or chemical reaction decreases the local Nusselt number. In [21] , the similarity solution was applied on unsteady boundary layer flow over a rotating sphere with a stagnation point. They concluded that the heat transfer obtained by CHF (constant heat flux) case is higher than the case of CWT (constant wall temperature). Takhar and Nath [22] applied a self-similar solution on the unsteady MHD flow over a rotating sphere with a stagnation point region. They showed that the surface shear stress in the x-direction vanishes at a certain value of the acceleration parameter, while the shear stress in the y-direction remains finite.
Ahmed and Mahdy [23] discussed the unsteady magnetohydrodynamic, heat and mass transfer over a rotating sphere with a stagnation point region in the presence of radiation effect. They found that the velocity profiles in x-direction increases as the buoyancy parameter increases. However, they neglected the case of nanofluids and the effects of slip and convective boundary conditions. So, this paper aims to generalize Ahmed and Mahdy [23] using Buongiorno’s model [2] . Effects of the slip and convective boundary conditions are taken into account in this study. Non-similar solutions and MATLAB function bvp4c are used to solve the governing equations. Comprehensive discussion with the physical meaning of the findings is performed. To best knowledge of the authors, this problem did not consider before.
2. Problem Formulation
Figure 1 displays the physical model for the unsteady magnetohydrodynamic
Figure 1. Physical model of the problem.
flow of an incompressible, laminar nanofluid over a rotating sphere with a constant angular velocity Ω in the stagnation region. It is assumed that there is a uniform magnetic field with strength B in the flow domain. A two-phase model (Buongiorno’s nanofluid model) is considered to simulate the nanofluid. The ambient velocity, temperature and nanoparticles volume fraction are denoted by
and
, respectively, while the noral fluxes of the nanoparticles vanish at the surface of the sphere. Also, a uniform slip and convective boundary conditions are imposed to the wall of the sphere. The thermal radiation is taken into account, but the Joule heating and viscous dissipation effects are neglected. All the fluid properties are constant except the nanofluid density which is approximated using a well-known Boussinesq approximation. Under all these assumptions, the continuity, x-component momentum, y-component momentum, energy and nanoparticles volume fractions equations are expressed as:
(1)
(2)
(3)
(4)
(5)
In Equations (1)-(4), u, v and w are the velocity components in x, y and z-direction, respectively, ϕ is the nanoparticle volume fraction, ρ is the density, ν is the kinematic viscosity, T is the temperature, βT is the thermal expansion, g is the gravity acceleration, σ is the electrical conductivity, B0 is the magnetic strength, k is the thermal conductivity, DB is Brownian diffusion coefficient, DT is thermophoretic diffusion, C is the specific heat, q is the radiative heat flux and the subscripts f, p and ∞ refer to base fluid, nanoparticle and free stream condition, respectively.
The radiative flux q can be expressed using the Rossel and approximation [24] as:
(6)
where,
is the Stefan-Boltzmann constant and
is the absorption coefficient. Expanding
using Taylor series about
and neglecting the higher terms, we obtain:
(7)
Substituting Equations (6) and (7) in Equation (4), we get:
(8)
The corresponding initial and boundary conditions are:
(9a)
(9b)
(9c)
where, is the dimensional slip parameter, μ is the dynamic viscosity, Ω is the angular velocity of the sphere and hf is the convective heat transfer.
Introducing the following dimensionless variables:
(10)
Substituting Equation (10) into Equations (2), (3), (5) and (8), the following non-similar equations are obtained:
(11)
(12)
(13)
(14)
where,
is the Grashof number,
is the buoyancy parameter,
is the Reynolds number,
is the magnetic field parameter ,
is the buoyancy ratio,
is the rotation parameter,
is the Brownian motion parameter,
is the thermophoreses parameter,
is the Prandtl number,
is the Lewis number.
The boundary conditions are converted to:
(15)
where,
is the slip parameter and
is the Biot number.
3. Numerical Solutions
In this section the MATLAB function bvp4c is applied to solve Equations (11)-(14) with the boundary conditions (15). The first step in this technique is the converting the system (11)-(14) to an ordinary differential equations system, so all the derivatives with respect to ξ in this system are dropped down, thus:
(16)
(17)
(18)
(19)
Subject to the boundary conditions:
(20)
In the second step, all the neglected terms in previous system are restored and the following new variables are introduced:
(21)
Thus, the governing equations are:
(22)
(23)
(24)
(25)
Subject to the boundary conditions:
(26)
There are four additional equations are needed with appropriate boundary conditions to evaluate the four new variables
and
. These can be obtained by differentiating Equations (22)-(26) with respect to ξ and neglecting
the terms
and
, it can be obtained:
(27)
(28)
(29)
(30)
Subjected to the boundary conditions:
(31)
Equations (16)-(19), with the boundary conditions (20) and Equations (27)-(30) with the boundary conditions (31) are solved numerically using MATLAB software. The value of
is set to be equal 10 and the step size is
. This method is found to be suitable and gives results those are very close to Ahmed and Mahdy [23] . As it can be noted from Figure 2, the present method compares very well with Ahmed and Mahdy [23] .
4. Results and Discussion
Figure 3 and Figure 4 display effects of the slip parameter δ on the dimensionless velocity in the x-direction and temperature. The results indicated that the velocity value increases as the slip parameter increases but the temperature distributions decreases hence the thickness of the thermal boundary layer increases. Physically, the increase in slip parameters causes an enhancement in the velocity gradients on the sphere surface which in turn increases the fluid activity, rate of heat transfer and decreases the temperature as well.
Figure 5 shows that the nanoparticle volume fraction profiles increase with increasing δ near the rotating sphere surface up to a certain value of η but beyond this point, the opposite trend is observed. In Figures 6-8, the influence of thermophoresis parameter Nt on the velocity in the x-direction, temperature
Figure 2. Comparison of the present results and those obtained by Ahmed and Mahdy [23] .
Figure 3. Profiles of the velocity in the x-direction for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 4. Temperature distributions for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 5. Profiles of the nanoparticle volume fraction for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 6. Profiles of the velocity in the x-direction for different values of Nt at Bi = 10, M = 5, R = 5, Le = 10, Pr = 10, Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 7. Temperature distributions for different values of Nt at Bi = 10, M = 5, R = 5, Le = 10, Pr = 10, Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 8. Profiles of the nanoparticle volume fraction for different values of Nt at Bi = 10, M = 5, R = 5, Le = 10, Pr = 10, Nb = Nr = 0.5,δ = 0.1, γ = λ = 5, ξ = 0.5.
and nanoparticle volume fraction profiles are illustrated. It is found that, the increase in Nt leads to an increase in the velocity profiles and temperature distributions but the nanoparticle volume fraction profiles decrease with the increase in δ near the sphere surface, whereas at
, the inverse behavior is observed. Also, the thermal boundary layer increases, slightly, along with the increase in Nt. In Figures 9-11, the effects of Biot number Bi on the velocity profiles in the x-direction, temperature distribution and the nanoparticle volume fraction profiles are presented. It can be noted that the increase in Bi enhances the velocity profiles and temperature distributions, while on contrary the profiles of nanoparticle volume fraction are reduced. Additionally, this mentioned behavior of ϕ is noted near the surface but far away from the sphere surface, the increase in Bi improves the distributions of ϕ. Here, it should be mentioned that when
, the constant wall temperature case (CWT) is obtained. From the
Figure 9. Profiles the velocity in x-direction for different values of Bi at M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 10. Temperature distributions for different values of Bi at M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 11. Profiles of the nanoparticle volume fraction for different values of Bi at M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
physical view, the increase in Bi enhances the temperature differences at the wall which in turn augments the mixed convection and consequently increases the velocity temperature distributions.
Figures 12-14 display effect of the radiation parameter R on the velocity in the x-direction, temperature and the nanoparticle volume fraction profiles. It is found that the increase in R leads to a reduction in the velocity and temperature profiles. Further, the nanoparticle volume fraction profiles decrease with increasing R near the surface up to a certain value of η but beyond this point, the opposite tendency is observed.
The effect of variations of Lewis Le on the velocity profiles in the x-direction for different values of slip parameter δ is examined with the help of Figure 15. It is found that the increase in Le resists the fluid flow in this direction regardless of the value of δ. This can be attributed to the velocity boundary layer thickness which increases as Le increases and therefore the velocity profiles decrease.
Figure 12. Profiles of the velocity in x-directions for different values of R at Bi = 10, M = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 13. Temperature distributions for different values of R at Bi = 10, M = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 14. Profiles of the nanoparticle volume fraction for different values of R at Bi = 10, M = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 15. Profiles of the velocity in the x-direction for different values of Le and δ at at Bi = 0.5, M = 5, R = 5, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 16 shows that as the Prandtl number Pr increases, the temperature distribution decreases. Here it should be mentioned that if
denoted the velocity boundary layer thickness and
represented the thermal boundary layer thickness, so if
, then
, if
then
and if
, then
. Consequently, the thermal boundary layer is inversely proportional with the values of Pr, so the temperature distributions are reduced as Pr increases.
Also, the same effect is observed for the increase in δ, moreover the thickness of the thermal boundary layer decreases. As can be observed in Figure 17 increasing value of the Pr the nanoparticle volume fraction profiles increase near the stretching sheet wall up to a certain value of η but beyond this point, the opposite trend is observed.
Figure 16. Temperature distributions for different values of Pr and δ at Bi = 0.5, M = 5, R = 5, Le = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 17. Profiles of the nanoparticle volume fraction for different values of Pr at Bi = 0.5, M = 5, R = 5, Le = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
In Figure 18 and Figure 19 respectively shows the effects of the rotation parameter λ and the buoyancy parameter γ on the velocity profiles in the x-direction. It clear increase in λ, γ leads to an increase in the velocity profiles .The physical explanation for this behavior is attributed to the fluid momentum which increases by increase rotation parameter.
The effects of the rotation parameter λ and the slip parameter δ on the local skin friction coefficient
are presented in Figure 20 and Figure 21. It is observed from these figures that as the rotation parameter λ increases of
increase Also, the same effect is observed for the increase in δ , but at
, (
) the
has different behavior.
Figure 18. Profiles of the velocity in x-directions for different values of λ at Bi = 0.5, M = R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 19. Profiles of the velocity in x-directions for different values of γ at Bi = 0.5, M = R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, λ = 5, ξ = 0.5.
Figure 20. Local skin friction coefficient for different values of λ and δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = 5, ξ = 0.5.
Figure 21. Local skin friction coefficient for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
The effects of slip parameter δ on the local Nusselt and Sherwood numbers are shown in Figures 22 and Figure 23 respectively. It is seen that the local Nusselt number increases with increasing δ, while local Sherwood number decrease. These behaviors is due to the fact that the increase in slip parameter improves the mixed convection case and thus the rate of heat transfer is enhanced
In Figure 24 and Figure 25, the behavior of local Nusselt number and Sherwood number for various values of Biot number Bi is noted. The results indicate that the increase in Bi number cause an enhancement in the temperature gradients at the surface and consequently the rate of heat transfer is augmented, while the opposite behavior
. Also thermal boundary layer thickness increases.
Figure 26 shows the influences of Brownian motion parameter Nb on the local Sherwood number
. It is observed that as the Brownian motion parameter increases of of the local Sherwood number increase. Physically, the increase in Brownian motion parameter increases the net movement of atoms or molecules from the region near the surface where the high concentration of nanoparticle to the region of low concentration of nanoparticle (free stream region) which results in the increase in Sherwood number.
Figure 27 and Figure 28 exhibit the effect of thermophoresis parameter Nt on the local Nusselt and Sherwood numbers, respectively. It is clear that the increase in Nt leads to a decrease in the local Nusselt and Sherwood numbers.
5. Conclusions
This paper discussed the unsteady MHD mixed convection flow along rotating spheres with a stagnation point in the presence of thermal radiation, slip and convective boundary conditions. The problem was formulated using Navier stokes equations with boundary layer approximations and then solved numerically using MATLAB software. The obtained results were presented in terms of
Figure 22. Local Nusselt number for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 23. Local Sherwood number for different values of δ at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, γ = λ = 5, ξ = 0.5.
Figure 24. Local Nusselt number for different values of Bi at M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 25. Local Sherwood number for different values of Bi at M = 5, R = 5, Le = 10, Pr = 10, Nt = Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 26. Local Sherwood number for different values of Nb at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nt = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 27. Local Nusselt number for different values of Nt at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
Figure 28. Local Sherwood number for different values of Nt at Bi = 0.5, M = 5, R = 5, Le = 10, Pr = 10, Nb = Nr = 0.5, δ = 0.1, γ = λ = 5, ξ = 0.5.
profiles of velocity in the x-direction, temperature and nanoparticle volume fraction as well as skin friction coefficient, local Nusselt number and local Sherwood number. The key findings of the present investigation are summarized as follows:
1) The increase in slip parameter enhances the velocity profiles, skin friction coefficient and local Nusselt number, while the temperature distributions and local Sherwood number take the inverse behaviors.
2) An improvement in profiles of velocity and temperature can be obtained by increasing the thermophores parameter, while the profiles of nanoparticle volume fraction near the surface, local Nusselt and Sherwood numbers are reduced.
3) The increase in Biot number causes an increase in the hot nanofluid at the surface and consequently the mixed convection is improved.
4) The increase in radiation parameter resists the fluid flow and decreases the fluid temperature.
5) Regardless the values of slip parameter, the increase in Lewis number reduces the fluid velocity and the fluid temperature decreases, gradually, as the Prandtl number increases.
Acknowledgements
The authors extend their appreciation to the editors and reviewers for their comments to improve the quality of the paper.
Nomenclature
a Velocity gradient the edge
B Magnetic field
Bi Biot number
cp Specific heat
DB Brownian diffusion coefficient
DT Thermophoretic diffusion coefficient
FD Imensionless stream function
g Gravity acceleration
GrR Grashof number
hf Convective heat transfer coefficient
k Thermal conductivity of base fluid
k* Absorption coefficient
Le Lewis number
M Magnetic field parameter
Nr Buoyancy ratio
Nb Brownian motion parameter
Nt Thermophoresis parameter
Nu Nusselt number
r Radius of the sphere
R Radiation parameter
ReR Reynolds number
S Velocity component in-direction y
Sh Sherwood numbers
t Time
T Local fluid temperature
Tf Temperature of the hot fluid
(u, v, w) Velocity components in the x, y and z directions
U Ambient velocity
(x, y, z) Cartesian coordinates
Greek Symbols
βT Coefficient of thermal expansion
μ Dynamic viscosity
ν Kinematic viscosity
λ Rotation parameter
σ Electrical conductivity
σ1 Stefan-Boltzman constant
ρf Density of base fluid
ρp Nanoparticle mass density
θ Dimensionless temperature
ϕ Dimensionless nanoparticle volume fraction
(ρc)f Heat capacity of the base fluid
(ρc)p Heat capacity of the nanoparticle material
γ Buoyancy parameter
δ Slip parameter
Ω Angular velocity
(ξ, η) Non-similarity variables
Subscripts
w Conditions at the surface
∞ Conditions in the free stream