Next Article in Journal
Single-Qubit Driving Fields and Mathieu Functions
Next Article in Special Issue
Modeling of Subcooled Flow Boiling with Nanoparticles under the Influence of a Magnetic Field
Previous Article in Journal
Symmetry at the Cellular Mesoscale
Previous Article in Special Issue
Magnetohydrodynamics Stagnation-Point Flow of a Nanofluid Past a Stretching/Shrinking Sheet with Induced Magnetic Field: A Revised Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple Slip Effects on Magnetohydrodynamic Axisymmetric Buoyant Nanofluid Flow above a Stretching Sheet with Radiation and Chemical Reaction

Department of Applied Mathematics, School of Science, Northwestern Polytechnical University, Dongxiang Road, Chang’an District, Xi’an 710129, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(9), 1171; https://doi.org/10.3390/sym11091171
Submission received: 25 July 2019 / Revised: 4 September 2019 / Accepted: 9 September 2019 / Published: 16 September 2019
(This article belongs to the Special Issue Future and Prospects in Non-Newtonian and Nanofluids)

Abstract

:
The present article investigates the effect of multiple slips on axisymmetric magnetohydrodynamics (MHD) buoyant nano-fluid flow over a stretching sheet with radiation and chemical effect. The non-linear partial differential equations were transformed to a non-linear control equation using an appropriate similarity transformation. The governing equations were solved through the finite element method. The influence of physical parameters such as multiple slips, magnetic, thermal radiation, Prandtl number, stretching, Brownian motion, thermophoresis, Schmidt number, Lewis number and chemical reaction on the radial velocity, temperature, solutal concentration and nano-fluid volume fraction profile were investigated. We noted that the boundary layers increases in the presence of multiple slip effects whereas, the effect of thermal slip on Nusselt number increases with the increasing values of magnetic and thermal radiation. To verify the convergence of the numerical solution, the computations were made by reducing the mesh size. Finally, our results are parallel to previous scholarly contributions.

1. Introduction

The study of the hydromagnetics flow of non-Newtonian fluids related to nano-particles has attracted many mathematicians, physician, engineers and modelers due to numerous significant applications. These applications spread through a wide range including neural-electronic interfaces and protein engineering, solar systems, nuclear reactors, cooling systems and so forth. As a result, nano-fluids have been examined both hypothetically and experimentally to attain higher degree characteristics of heat transfer in base fluids. Studies have also focused on improved energy efficiency in an assortment of thermal exchange system for various mechanical, physics and engineering applications. These novel kinds of fluids have attracted numerous researchers in view of its improved thermal features. Compared with basic fluids such as water and oil, nano-fluids possess higher thermal conductivity and convection heat transfer properties. Choi and Eastman [1] empirically established that adding a small amount of these particles would significantly improve their effective thermal conductivity. Buongiorno [2] in his non-homogenous equilibrium model reveals that the abnormal increase in the thermal conductivity occurs due to the presence of Brownian motion and thermophoretic diffusion of nanoparticles. A wide range of current and future applications of nano-fluids are prescribed [3,4,5,6,7,8]. Nield and Kuznetsov [9] used the Buongiorno model to study the nano-fluids with respect to free convection flow problem in the permeable medium. Similarly Salleh et al. [10] conducted numerical analysis of the boundary layer heat flow adjacent to a thin needle in the nano-fluid in the presence of heat source and chemical reaction.
Kuznetsov and Nield [11] analyzed boundary layer flow through a plate and later through a porous medium [12]. Likewise, the flow of nano-fluids on a moving plate in uniform free flow is studied by Bachok et al. [13].
The effect of magnetic field on free convection heat transfer of nano-fluid in sinusoidal wall channel is studied by Rashidi et al. [14]. Numerical study on the flow of nano-fluids in rotating systems with permeable films is carried out by Sheikholeslami et al. [15]. Hatami et al. [16] analyzed nano-fluid flow in axisymmetric porous channels with expansion or shrinkage wall.
Magnetohydrodynamics plays a key role in the field of engineering, agriculture and metal processing. The application of mixed convection problem under the influence of the magnetic field in geophysics and astrophysics has attracted the attention of many researchers. Soundalgekar et al. [17] studied the free convection effect of vertical plates with transverse magnetic fields on Stokes problem, while Elbasheshy [18] studied the problems of MHD heat and mass transfer and the role of vertical plates along with the buoyancy effect of heat and species diffusion. In the MHD flow, the inductive current in the fluid creates a force that changes the flow field. Mohanty et al. [19] studied the reaction of the chemical to MHD fluid with heat generation/absorption on the stretching sheet. In another study, M. Ramzan et al. [20] examined the effects of MHD, Joule heating and thermal radiation on an incompressible Jeffrey nano-fluid flow over a linearly stretched surface. MHD and heat transfer flow on a vertical stretching sheet with heat sink or source effect was studied by M. Alfari et al. [21]. Steady two-dimensional laminar incompressible MHD flow over an exponentially shrinking sheet with the effects of slip conditions and viscous dissipation was examined by Ali Lund et al. [8]. Najib et al. [22] examined the stagnation point flow and mass transfer with chemical reaction through stretch/shrink cylinder. The heat and mass transfer flow of MHD micro-polar fluid in the presence of heat sources was investigated by Mishra et al. [23].
Thermal system performance using common base fluids is relatively poor. One way to improve the performance of these systems is to suspend metal nano-particles in the base fluid. Thermal radiation for heat transfer and mass transfer on tensile surfaces plays an important role because of its numerous applications in engineering and physics fields such as satellites, solar energy, geophysics, geophysics, spacecraft and fossil fuels combustion. We know that concentrations and temperature gradients describe mass and energy flux. These effects have significant role in density difference flow. It is of great help to manufacture, pollute and discover drugs of optical fibers and others. Zaimi [24] sought the influence of heat transfer on the nano-fluid flow by means of permeable stretching/shrinking sheet. Ahmed et al. [25] presented the unsteady boundary layer flow and heat transfer of power-law fluid model over a radially stretching sheet. Heat transfer of viscous fluid over a non-linear radially stretching sheet was studied by Azeem et al. [26]. Ilyas Khan et al. [27] introduced a mathematical model of a convection flow of MHD nano-fluid in a channel embedded in a porous medium. The influence of the external magnetic field and thermal radiation on heat transfer intensification of nano-fluid in a porous curved enclosure is simulated by M. Sheikholeslami et al. [5]. Steady 2-D MHD free convective boundary-layer flows of an electrically conducting nano-fluid over a non-linear stretching sheet, taking into account the chemical reaction and heat source/sink, was investigated by D. Makinde et al. [28]. Multiple Slip Effects on MHD Unsteady Flow Heat and Mass transfer impinging on permeable stretching sheet with radiation is studied by Fazle and Stanford Shateyi [29]. Chemical reaction effect of an axisymmetric flow over a radially stretching sheet was considered by B. Nayak [30].
To the author’s best knowledge, no previous studies have thus far examined the multiple slips on hydromagnetics axisymmetric nano-fluid steady flow, heat and mass transfer influenced by radiation, buoyancy and chemical effect in a permeable frame of reference. We provide numerical solutions for some special cases, while the physical interpretation for the various parameters have been explained with the help of graphs. The current study investigates the analysis of mass and heat transfer in the axisymmetric MHD nano-fluid flow over a stretching sheet in the presence of multiple slip conditions with buoyancy effect, chemical reaction and thermal radiation. The governing partial differential equations of the model are converted to ordinary differential equations ODE through proper similarity transformations. The ODE’s are then solved through finite element method. The numerical results of radial velocity distribution, temperature profile, solutal concentration and nano-fluid volume fraction with different physical parametric values are discussed graphically. The present results are compared with those of published papers [31,32,33] and found a good agreement.

2. Mathematical Formulation

We consider the steady axisymmetric MHD flow of an electrically conducting incompressible viscous flow immersed in nano-fluid over a stretching sheet. The geometry of the problem is illustrated in Figure 1. The flow of conductive fluids is caused by the sheet stretching in the U ( r ) = c 0 r direction, where c 0 is a dimensional constant. It is considered that there is no nanoparticle flux across the wall and the surface is stretching along z-direction. Furthermore, it is assumed that T w , ψ w and C w denotes the stretching sheet temperature, solute concentration and the nano-particle fraction at the sheet respectively. At the infinity the T , C and ψ are the ambient temperature, ambient solute concentration and ambient nano-particle volume fraction, respectively. In the current study, B ( r ) = B 0 r summarizes the strength of magnetic field, where B 0 is a uniform magnetic field strength. The magnetic field acts along normal to the sheet in the positive z-direction. When the magnetic field is varying in the radial direction, it also vary in the vertical direction, since the divergence of B should be zero, that is, d i v B = 0 . Under the boundary layer approximation, the equations of mass conservation, linear momentum conservation, energy conservation and nano-particle volume fraction (see References [30,34,35,36]) can be obtained as:
u r + u r + w z = 0 ,
u u r + w u z = ν 2 u z 2 σ B 2 ( r ) u ρ
u w z + w w z = ν 2 w z 2 + g β T ( T T ) + g β C ( C C ) + g β ψ ( ψ ψ )
u T r + w T z = α 2 T z 2 + τ D B ψ z T z + D T T T z 2 1 ρ c p q r z
u C r + w C z = D s 2 C z 2 R ( r ) ( C C )
u ψ r + w ψ z = D B 2 ψ z 2 + D T T 2 T z 2
where u and w are supposed to be the velocity components along r and z respectively, ν , σ and ρ are the kinetic viscosity, electric conductivity and viscosity of fluid, respectively. D B , D T , D s , R ( r ) , are the Brownian diffusion, thermophoretic diffusion, the solutal and chemical reaction, respectively.
The quantity q r is the relative heat flux. Applying the Rosseland approximation [37], q r can be defined as:
q r = 4 σ 3 k T 4 y ,
where σ and k are the Stefan-Boltzmann constant and the mean absorbtion coefficient, respectively. If the variation in temperature across the flow is small, T 4 maybe expended using Tayler series at a free stream temperature T after neglecting higher-order terms as
T 4 = 4 T 4 3 T 4 .
Using Equations (7) and (8), we get
q r y = 16 σ T 4 3 k 2 T 2 y .
To simplify the mass equation, the chemical reaction R ( r ) must be constant. This constraint holds if R ( r ) has the following form [38]:
R ( r ) = R 0 a r n
where R 0 is constant and a > 0 is the stretching constant. The required physical boundary conditions of the above governing Equations are as follows [39]:
u = U ( r ) = c 0 r , w = 0 , T = T w + T s l i p , C = C w + C s l i p , C w = C + c 0 r n , ψ = ψ w + ψ s l i p a t z = 0 u 0 , w 0 , T T , C C , ψ ψ , a s z ,
where T s l i p = D 1 T y is the thermal slip, D 1 is the temperature slip factor, C s l i p = E 1 C y is the solutal slip, E 1 is the solutal concentration slip factor, ψ s l i p = F 1 ψ y is the slip velocity of the nanoparticle, F 1 is the nanoparticle concentration slip factor, C w is the variable surface concentration, c 0 is a constant, n is power-law exponent signifies the change of amount of solute in r –direction.
In order to solve Equations (1)–(6), we introduce the following similarity transformations (see References [30,34]):
u = c 0 r f ( η ) , w = c o ν z f ( η ) , η = c o ν z , θ ( η ) = T T T w T , S ( η ) = C C C w C , ϕ ( η ) = ψ ψ ψ w ψ
In view of the similarity transformation equation, the partial nonlinear differential Equations (1)–(6) transform into the following system of nonlinear ODEs:
f + f f f 2 M f = 0 ,
g + f g g 2 + λ 1 θ + λ 2 S + λ 3 ϕ = 0 ,
( 1 + N r ) 1 P r θ f θ + f θ + N b ϕ θ + N t θ 2 = 0
1 S c S n f S + f S γ S = 0
ϕ L e f ϕ + L e f ϕ + N t N b θ = 0
and the transformed boundary conditions of Equation (10) are:
f ( 0 ) = 0 , f ( 0 ) = 1 + S f f ( 0 ) , θ ( 0 ) = 1 + S t θ ( 0 ) , S ( 0 ) = 1 + S p S ( 0 ) , ϕ ( 0 ) = 1 + S g ϕ ( 0 ) , f ( ) 0 , θ ( ) 0 , S ( ) 0 , ϕ ( ) 0 ,
The primes shows the differentiation with respect to η . The parameters in Equations (12)–(17) are defined as: M = σ B 0 2 ρ c o , P r = κ ρ c p , N b = τ D B ( ψ w ψ ) ν , N t = τ D T ( T w T ) ν T , S c = ν D , L e = σ D B , N r = 16 σ T 3 3 k f K , γ = R ( r ) c o are magnetic, Prandtl number, Brownian, thermophoresis, Schmidt number, Lewis number, thermal radiation and chemical reaction parameter, respectively. λ 1 = g β T T 0 / a ν , λ 2 = g β C C 0 / a ν , λ 3 = g β ψ ψ 0 / a ν are the buoyancy parameters. S f , S t , S p and S g are the hydrodynamic slip, thermal slip, solutal slip and nano-fluid slip, respectively.

3. Finite Element Method Solutions

We applied the finite element method to solve nonlinear differential Equations (12)–(17), with the boundary conditions (12). Finite element method is extremely efficient and was applied to solve different problems in fluid mechanics, heat transfer, chemical processing, solid mechanics and other fields also (see References [40,41,42]). In order to find numerical solutions, we consider
f = h .
substitute Equation (13) in to Equation (8)–(11) the system reduce to:
h + f h h 2 M h + λ 1 θ + λ 2 S + λ 3 ϕ = 0 ,
( 1 + N r ) 1 P r θ h θ + f θ + N b ϕ θ + N t θ 2 = 0 ,
1 S c S n h S + f S γ S = 0 ,
ϕ L e h ϕ + L e f ϕ + N t N b θ = 0 ,
with following reduce boundary condition
f ( 0 ) = 0 , h ( 0 ) = 1 + S f h ( 0 ) , θ ( 0 ) = 1 + S t θ ( 0 ) , S ( 0 ) = 1 + S p S ( 0 ) , ϕ ( 0 ) = 1 + S g ϕ ( 0 ) , h ( ) 0 , θ ( ) 0 , S ( ) 0 , ϕ ( ) 0
For the computation purposes, the η at ∞ is chosen large enough, therefore the numerical solution depicts no perceptible variation for η larger than η m a x . Subject to the boundary condition the η m a x has been fixed at η m a x = 12 .

3.1. Variational Formulations

The variational form associated with Equations (13)–(18) over a linear element Ω e = ( η e , η e + 1 ) is given by
η e η e + 1 s 1 { f h } d η = 0 ,
η e η e + 1 s 2 { h + f h M h h 2 + λ 1 θ + λ 2 S + λ 3 ϕ } d η = 0 ,
η e η e + 1 s 3 { ( 1 + N r ) 1 P r θ h θ + f θ + N b ϕ θ + N t θ 2 } d η = 0 ,
η e η e + 1 s 4 { 1 S c S n h S + f S γ S } d η = 0 ,
η e η e + 1 s 5 { ϕ L e h ϕ + L e f ϕ + N t N b θ } d η = 0 ,
where s 1 , s 2 , s 3 , s 4 and s 5 are arbitrary shape function or trial function.

3.2. Finite Element Formulations

The Equation of finite element model is obtained by replacing finite element approximation of the following form in the Equations (19)–(23).
f = j = 1 2 f j ψ j , h = j = 1 2 h j ψ j , θ = j = 1 2 θ j ψ j , S = j = 1 2 S j ψ j , ϕ = j = 1 2 ϕ j ψ j ,
with s 1 = s 2 = s 3 = s 4 = s 5 = ψ i ( i = 1 , 2 ) , where the shape function ψ i for a line element Ω e = ( η e , η e + 1 ) are given by
ψ 1 = η e + 1 η η e + 1 η e , ψ 2 = η η e η e + 1 η e , η e η η e + 1 .
The finite element method model equations are, therefore given by
F 11 F 12 F 13 F 14 F 15 F 21 F 22 F 23 F 24 F 25 F 31 F 32 F 33 F 34 F 35 F 41 F 42 F 43 F 44 F 42 F 51 F 52 F 53 F 54 F 55 f h θ S ϕ = b 1 b 2 b 3 b 4 b 5
where F m n and b m (m, n = 1, 2, 3, 4, 5) are defined as
F i j 11 = η e η e + 1 ψ i d ψ j d η d η , F i j 12 = η e η e + 1 ψ i ψ j d η , F i j 13 = F i j 14 = 0 , F i j 15 = F i j 21 = 0 , F i j 22 = η e η e + 1 d ψ i d η d ψ j d η d η + η e η e + 1 f ¯ ψ i d ψ j d η d η η e η e + 1 h ¯ ψ i ψ j d η M η e η e + 1 ψ i ψ j d η , F i j 23 = λ 1 η e η e + 1 ψ i ψ j d η , F i j 24 = λ 2 η e η e + 1 ψ i ψ j d η , F i j 25 = λ 3 η e η e + 1 ψ i ψ j d η , F i j 31 = 0 , F i j 32 = 0 ,
F i j 33 = ( 1 + N r ) 1 P r η e η e + 1 d ψ i d η d ψ j d η d η η e η e + 1 h ¯ ψ i ψ j d η + η e η e + 1 f ¯ ψ i d ψ j d η d η + N b η e η e + 1 ϕ ¯ ψ i d ψ j d η d η + N t η e η e + 1 θ ¯ ψ i d ψ j d η d η , F i j 34 = 0 , F i j 35 = 0 , F i j 41 = F i j 42 = 0 , F i j 43 = 0 ,
F i j 44 = 1 S c η e η e + 1 d ψ i d η d ψ j d η d η + η e η e + 1 f ¯ ψ i d ψ j d η d η n η e η e + 1 h ¯ ψ i ψ j d η d η γ η e η e + 1 ψ i ψ j d η d η , F i j 45 = F i j 51 = F i j 52 = F i j 54 = 0 , F i j 53 = N t N b η e η e + 1 d ψ i d η d ψ j d η d η
F i j 55 = η e η e + 1 d ψ i d η d ψ j d η d η + L e η e η e + 1 f ¯ ψ i d ψ j d η d η L e η e η e + 1 h ¯ ψ i ψ j d η d η ,
and
b i 1 = 0 , b i 2 = ψ d h d η η e η e + 1 , b i 3 = ( 1 + N r ) 1 P r ψ d θ d η η e η e + 1 , b i 4 = 1 S c ψ d S d η η e η e + 1 , b i 5 = ψ d ϕ d η η e η e + 1 N t N b ψ d θ d η η e η e + 1 ,
where f ¯ = j = 1 2 f ¯ j ψ j , h ¯ = j = 1 2 h ¯ j ψ j , θ ¯ = j = 1 2 θ ¯ j ψ j and ϕ ¯ = j = 1 2 ϕ ¯ j ψ j are considered to be known. After assembling all the element equations the resulting system of equations obtained are nonlinear, therefore an iterative scheme is used in the solution. The system is linearized by incorporating the function f ¯ , h ¯ , θ ¯ and ψ ¯ known at a lower iteration level. Computations for the functions f , h , θ , S , ϕ are carried out for high level. This process is repeated till the desired accuracy of 0.00005 is obtained.

4. Results and Discussion

Research goal of the current study was to identify the role of mass, heat transfer, chemical reaction and thermal radiation in the axisymmetric MHD boundary layer flow of nano-fluid in the presence of multiple slip effects over a stretching sheet. The finite element method was used to solve the control model equations numerically. In this part, we discuss the results of velocity, temperature, solute concentration and the distribution of the volume fraction of nano-fluid by the proposed numerical scheme. We describe the effect of above profiles on the various dimensionless parametric values such as thermal radiation N r , magnetic M, Prandtl number P r , stretching parameter n, Brownian motion parameter N b , buoyancy parameters λ 1 , λ 2 , λ 3 , thermophoresis N t , Schmidt number S c , Lewis number and chemical reaction γ parameters in the presence of multiple slips paramters ( S f , S t , S p , s g ).
The calculations were carried out for fixed values of the parameters P r = 0.773 , M = 2 , N b = 0.1 , N t = 0.1 , S c = 0.1 , R ( r ) = 0.1 , L e = 3 , N r = 1 n = 1 , S f = S t = S p = S g = 0.5 , λ 1 = λ 2 = λ 3 = 0.4 . Convergence of the finite element method is reflected in Table 1, which explicates no significant variation in the values of f, h, θ , S and ψ when the number of elements passes through 320. In order to prove the validity of the numerical results in Table 2, Table 3 and Table 4, current study archived an excellent correlation and exact solution as in the studies of Crane [31] data, Gireesha et al. [32], Mudassar et al. [33], Mabood and Das [43], Fazle Mabood et al. [29] and Ali [44]. An excellent correlation has been obtained and a grid invariance test has been conducted to maintain 4 decimal point accuracy. Table 5 displays the computed skin friction, heat and mass transfer coefficients which are expressed as | f ( 0 ) | , | θ ( 0 ) | and | S ( 0 ) | , respectively for different values of M, γ , P r , N r , λ 1 , λ 3 and λ 3 . It was noted that the skin friction increases while heat and mass transfer coefficient decreases as M increases. Also mass transfer coefficient increases as γ increases, heat transfer coefficient increases with increase in P r while the heat transfer coefficient decreases as N r accelerates. The Skin friction coefficient retarded whereas, heat and mass transfer coefficients amplified with the enhancement in the buoyancy parameters ( λ 1 , λ 2 , λ 3 ). The velocity profile decreases in both cases of hydrodynamic slip ( S f = 1 ) and no hydrodynamic slip ( S f = 0 ) with the increase in magnetic parameter, the results are depicted in Figure 2. These results are in accordance with the previous contributions [29]. This shows that the Lorentz force slows down the speed of the fluid in a radial direction. However, with existence of the hydrodynamic slip as can be clearly seen in Figure 2 that velocity boundary layer increased. Figure 3 shows that the temperature increases with the increase in the value of the magnetic parameter in the presence and absence of thermal slip conditions. Figure 4 demonstrates that the effect of solutal slip ( S p ) for concentration increased in the boundary layer while increasing in the magnetic parameter M. The velocity distribution increases as the value of buoyancy parameters ( λ 1 , λ 2 , λ 3 ) accelerated in the presence and absence of multiple slip conditions (hydrodynamic slip and no hydrodynamic slip, thermal slip and no thermal slip, solute slip and no solute slip) the results maybe seen in Figure 5, Figure 6 and Figure 7. The variation of temperature profile with the increase in thermal radiation parameter N r thermal slip and no thermal slip conditions are described in Figure 8. The temperature distribution increases with the increase of thermal radiation, which leads to the enormous heat transfer rate through the nano-fluid and the surface. Effect of N b on concentration profile in nano-fluid slip condition ( S g ) is reflected in Figure 9, the concentration profiles elucidate an opposite behavior for growing values of N b that the thickness of nanoparticle concentration boundary layer reduces with the increasing values for N b . It is observed that the increase in the number of Schmidt slows down the concentration distribution of nano-particles, which led to delays in concentration and heat transfer rates. The effects of solutal slip ( S p ) was examined to reduce concentration on the boundary layer with the increase in the stretching parameter n, Schmidt number and chemical reaction parameter, results are graphically depicted in Figure 10, Figure 11 and Figure 12.
Effect of Lewis number ( L e ) in nano-fluid slip condition on nano-fluid concentration profile is shown in Figure 13. Increase in Lewis number makes the boundary layer thinner. Lewis number is directly proportional to the Schmidt number and inversely proportional to the Brownian diffusion coefficient, which will cause the solute boundary layer to become‘thinner. Figure 14 exhibits the nature of the Skin friction coefficient with the magnetic, hydrodynamic slips and buoyancy parameters. We observed in Figure 14 that increase in the magnetic field strength M and λ 1 strongly suppresses the Skin friction coefficient. Skin friction coefficient on the stretching surface is significantly enhanced with the increase in the values of the hydrodynamic slip parameters, thermal buoyancy and solutal buoyancy parameters. Effect of M, N r , P r , λ 2 and thermal slip S t on Nusselt number are shown in Figure 15. We noticed in this figure that the effect of thermal slip ( S t ) on Nusselt number increases with the increasing values of magnetic and thermal radiation parameters, while it decreases with increasing values of buoyancy parameter and Prandtl number. From the figures, above, we see that in the existence of the hydrodynamic, thermal, solute and nano-fluid slip the velocity, temperature, solute concentration and nano-fluid volume fraction boundary layers all increased.

5. Conclusions

The present study investigated the axisymmetric MHD buoyant nano-fluid incompressible viscous electrically conducting flow in the presence of multiple slips, chemical reaction and thermal radiation over a stretching sheet. The current study leads to the following conclusions:
  • In the existence of the hydrodynamic, thermal, solute and nano-fluid slip the velocity, temperature, solute concentration and nano-fluid volume fraction boundary layers increases.
  • The radial velocity profile decreases while temperature and solute concentration increases with the increase in magnetic field and slip conditions.
  • The velocity profile increases with increasing buoyancy parameters ( λ 1 , λ 2 , λ 3 ).
  • After enhancement of magnetic and thermal radiation parameter the effect of thermal slip for temperature profile increased the boundary layer.
  • Solute concentration decreases with an increase in the stretching parameter, solutal slip, Schmidt number and chemical reaction parameter.
  • The nano-fluid volume fraction profile retards with large values of Brownian motion parameter and Lewis number.
  • Nusselt number increases in the presence of thermal slip with the increasing values of magnetic and thermal radiation parameters. Whereas, the effect is reversed for enhancing the values of solutal buoyancy parameter and Prandtl number.

Author Contributions

All authors contributed equally.

Funding

This research was funded by National Natural Science Foundation of China No. [11971386].

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u,wvelocity components
Bmagnetic field strength
B 0 uniform magnetic field strength
N b Brownian parameter
S c Schmidt number
N t thermophoresis parameter
L e Lewis number
q r heat flux
D B Brownian diffusion
D T thermophoretic diffusion
D s solute parameter
R ( r ) chemical reaction
T w constant temperature
C w nanoparticle volume fraction
T ambient temperature
C ambient nanoparticle volume fraction
Mmagnetic parameter
P r Prandtl number
N r thermal radiation
c 0 a constant
npower-law exponent
S f hydrodynamic slip
S t thermal slip
S p solute slip
S g nano-fluid slip
D 1 temperature slip factor
E 1 solutal concentration slip factor
F 1 nanoparticle concentration slip factor
k mean absorbtion coefficient

Greek Symbol

τ ratio of nanoparticle heat capacity
ν kinetic viscosity
σ electric conductivity
σ Stefan-Boltzmann constant
α thermal diffusivity
ρ viscosity of the fluid
ψ w solute concentration
γ chemical reaction
ψ ambient solute concentration
λ 1 , λ 2 , λ 3 Buoyancy parameters

References

  1. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Technical Report; Argonne National Lab.: Springfield, IL, USA, 1995. [Google Scholar]
  2. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  3. Turkyilmazoglu, M. Exact analytical solutions for heat and mass transfer of MHD slip flow in nanofluids. Chem. Eng. Sci. 2012, 84, 182–187. [Google Scholar] [CrossRef]
  4. Kandelousi, M.S. Nanofluid Heat and Mass Transfer in Engineering Problems; BoD-Books on Demand: Norderstedt, Germany, 2017. [Google Scholar]
  5. Sheikholeslami, M.; Shamlooei, M.; Moradi, R. Numerical simulation for heat transfer intensification of nanofluid in a porous curved enclosure considering shape effect of Fe3O4 nanoparticles. Chem. Eng. Process. Intensif. 2018, 124, 71–82. [Google Scholar] [CrossRef]
  6. Sheikholeslami, M.; Shehzad, S. Magnetohydrodynamic nanofluid convection in a porous enclosure considering heat flux boundary condition. Int. J. Heat Mass Transf. 2017, 106, 1261–1269. [Google Scholar] [CrossRef]
  7. Sheikholeslami, M.; Sadoughi, M. Simulation of CuO-water nanofluid heat transfer enhancement in presence of melting surface. Int. J. Heat Mass Transf. 2018, 116, 909–919. [Google Scholar] [CrossRef]
  8. Ali Lund, L.; Omar, Z.; Khan, I.; Raza, J.; Bakouri, M.; Tlili, I. Stability Analysis of Darcy-Forchheimer Flow of Casson Type Nanofluid Over an Exponential Sheet: Investigation of Critical Points. Symmetry 2019, 11, 412. [Google Scholar] [CrossRef]
  9. Kuznetsov, A.; Nield, D. The Cheng–Minkowycz problem for natural convective boundary layer flow in a porous medium saturated by a nanofluid: A revised model. Int. J. Heat Mass Transf. 2013, 65, 682–685. [Google Scholar] [CrossRef]
  10. Salleh, S.N.A.; Bachok, N.; Arifin, N.M.; Ali, F.M. Numerical Analysis of Boundary Layer Flow Adjacent to a Thin Needle in Nanofluid with the Presence of Heat Source and Chemical Reaction. Symmetry 2019, 11, 543. [Google Scholar] [CrossRef]
  11. Kuznetsov, A.; Nield, D. Natural convective boundary-layer flow of a nanofluid past a vertical plate. Int. J. Therm. Sci. 2010, 49, 243–247. [Google Scholar] [CrossRef]
  12. Nield, D.; Kuznetsov, A. The Cheng–Minkowycz problem for natural convective boundary-layer flow in a porous medium saturated by a nanofluid. Int. J. Heat Mass Transf. 2009, 52, 5792–5795. [Google Scholar] [CrossRef]
  13. Bachok, N.; Ishak, A.; Pop, I. Boundary-layer flow of nanofluids over a moving surface in a flowing fluid. Int. J. Therm. Sci. 2010, 49, 1663–1668. [Google Scholar] [CrossRef]
  14. Rashidi, M.M.; Nasiri, M.; Khezerloo, M.; Laraqi, N. Numerical investigation of magnetic field effect on mixed convection heat transfer of nanofluid in a channel with sinusoidal walls. J. Magn. Magn. Mater. 2016, 401, 159–168. [Google Scholar] [CrossRef]
  15. Sheikholeslami, M.; Ganji, D. Numerical investigation for two phase modeling of nanofluid in a rotating system with permeable sheet. J. Mol. Liq. 2014, 194, 13–19. [Google Scholar] [CrossRef]
  16. Hatami, M.; Sheikholeslami, M.; Ganji, D. Nanofluid flow and heat transfer in an asymmetric porous channel with expanding or contracting wall. J. Mol. Liq. 2014, 195, 230–239. [Google Scholar] [CrossRef]
  17. Soundalgekar, V.; Gupta, S.; Aranake, R. Free convection effects on MHD Stokes problem for a vertical plate. Nucl. Eng. Des. 1979, 51, 403–407. [Google Scholar] [CrossRef]
  18. Elbashbeshy, E. Heat and mass transfer along a vertical plate with variable surface tension and concentration in the presence of the magnetic field. Int. J. Eng. Sci. 1997, 35, 515–522. [Google Scholar] [CrossRef]
  19. Jena, S.; Mishra, S.; Dash, G. Chemical reaction effect on MHD Jeffery fluid flow over a stretching sheet through porous media with heat generation/absorption. Int. J. Appl. Comput. Math. 2017, 3, 1225–1238. [Google Scholar] [CrossRef]
  20. Ramzan, M.; Bilal, M.; Chung, J.D.; Mann, A. On MHD radiative Jeffery nanofluid flow with convective heat and mass boundary conditions. Neural Comput. Appl. 2018, 30, 2739–2748. [Google Scholar] [CrossRef]
  21. Alarifi, I.M.; Abokhalil, A.G.; Osman, M.; Lund, L.A.; Ayed, M.B.; Belmabrouk, H.; Tlili, I. MHD Flow and Heat Transfer over Vertical Stretching Sheet with Heat Sink or Source Effect. Symmetry 2019, 11, 297. [Google Scholar] [CrossRef]
  22. Najib, N.; Bachok, N.; Arifin, N.M.; Ishak, A. Stagnation point flow and mass transfer with chemical reaction past a stretching/shrinking cylinder. Sci. Rep. 2014, 4, 4178. [Google Scholar] [CrossRef]
  23. Mishra, S.; Dash, G.; Pattnaik, P. Flow of heat and mass transfer on MHD free convection in a micropolar fluid with heat source. Alex. Eng. J. 2015, 54, 681–689. [Google Scholar] [CrossRef]
  24. Zaimi, K.; Ishak, A.; Pop, I. Boundary layer flow and heat transfer over a nonlinearly permeable stretching/shrinking sheet in a nanofluid. Sci. Rep. 2014, 4, 4404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ahmed, J.; Mahmood, T.; Iqbal, Z.; Shahzad, A.; Ali, R. Axisymmetric flow and heat transfer over an unsteady stretching sheet in power law fluid. J. Mol. Liq. 2016, 221, 386–393. [Google Scholar] [CrossRef]
  26. Shahzad, A.; Ahmed, J.; Khan, M. On heat transfer analysis of axisymmetric flow of viscous fluid over a nonlinear radially stretching sheet. Alex. Eng. J. 2016, 55, 2423–2429. [Google Scholar] [CrossRef] [Green Version]
  27. Khan, I.; Alqahtani, A.M. MHD Nanofluids in a Permeable Channel with Porosity. Symmetry 2019, 11, 378. [Google Scholar] [CrossRef]
  28. Makinde, O.D.; Mabood, F.; Ibrahim, M.S. Chemically reacting on MHD boundary-layer flow of nanofluids over a non-linear stretching sheet with heat source/sink and thermal radiation. Therm. Sci. 2018, 22, 495–506. [Google Scholar] [CrossRef]
  29. Mabood, F.; Shateyi, S. Multiple Slip Effects on MHD Unsteady Flow Heat and Mass Transfer Impinging on Permeable Stretching Sheet with Radiation. Model. Simul. Eng. 2019, 2019, 3052790. [Google Scholar] [CrossRef]
  30. Nayak, B.; Mishra, S.; Krishna, G.G. Chemical reaction effect of an axisymmetric flow over radially stretched sheet. Propuls. Power Res. 2019, 8, 79–84. [Google Scholar] [CrossRef]
  31. Crane, L.J. Flow past a stretching plate. Z. FÜR Angew. Math. Und Phys. Zamp 1970, 21, 645–647. [Google Scholar] [CrossRef]
  32. Gireesha, B.; Ramesh, G.; Bagewadi, C. Heat transfer in MHD flow of a dusty fluid over a stretching sheet with viscous dissipation. J. Appl. Sci. Res. 2012, 3, 2392–2401. [Google Scholar]
  33. Jalil, M.; Asghar, S.; Yasmeen, S. An exact solution of MHD boundary layer flow of dusty fluid over a stretching surface. Math. Probl. Eng. 2017, 2017, 2307469. [Google Scholar] [CrossRef]
  34. Nawaz, M.; Hayat, T. Axisymmetric stagnation-point flow of nanofluid over a stretching surface. Adv. Appl. Math. Mech. 2014, 6, 220–232. [Google Scholar] [CrossRef]
  35. Mustafa, M.; Khan, J.A.; Hayat, T.; Alsaedi, A. Analytical and numerical solutions for axisymmetric flow of nanofluid due to non-linearly stretching sheet. Int. J. Non-Linear Mech. 2015, 71, 22–29. [Google Scholar] [CrossRef]
  36. Mustafa, M.; Hayat, T.; Alsaedi, A. Axisymmetric flow of a nanofluid over a radially stretching sheet with convective boundary conditions. Curr. Nanosci. 2012, 8, 328–334. [Google Scholar] [CrossRef]
  37. Pal, D. Heat and mass transfer in stagnation-point flow towards a stretching surface in the presence of buoyancy force and thermal radiation. Meccanica 2009, 44, 145–158. [Google Scholar] [CrossRef]
  38. Almakki, M.; Dey, S.; Mondal, S.; Sibanda, P. On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects. Entropy 2017, 19, 168. [Google Scholar] [CrossRef]
  39. Mabood, F.; Khan, W.; Ismail, A.M. Multiple slips effects on MHD Casson fluid flow in porous media with radiation and chemical reaction. Can. J. Phys. 2015, 94, 26–34. [Google Scholar] [CrossRef]
  40. Gupta, D.; Kumar, L.; Bég, O.A.; Singh, B. Finite-element simulation of mixed convection flow of micropolar fluid over a shrinking sheet with thermal radiation. Proc. Inst. Mech. Eng. Part J. Process. Mech. Eng. 2014, 228, 61–72. [Google Scholar] [CrossRef]
  41. Gupta, D.; Kumar, L.; Bég, O.A.; Singh, B. Finite element analysis of MHD flow of micropolar fluid over a shrinking sheet with a convective surface boundary condition. J. Eng. Thermophys. 2018, 27, 202–220. [Google Scholar] [CrossRef]
  42. Madhu, M.; Kishan, N. Finite Element Analysis of MHD viscoelastic nanofluid flow over a stretching sheet with radiation. Procedia Eng. 2015, 127, 432–439. [Google Scholar] [CrossRef]
  43. Mabood, F.; Das, K. Melting heat transfer on hydromagnetic flow of a nanofluid over a stretching sheet with radiation and second-order slip. Eur. Phys. J. Plus 2016, 131, 3. [Google Scholar] [CrossRef]
  44. Ali, M.E. Heat transfer characteristics of a continuous stretching surface. Wärme-und Stoffübertragung 1994, 29, 227–234. [Google Scholar] [CrossRef]
Figure 1. Physical configuration and coordinate system.
Figure 1. Physical configuration and coordinate system.
Symmetry 11 01171 g001
Figure 2. Effect of M on velocity distribution in the presence and absence of hydrodynamic slip condition S f .
Figure 2. Effect of M on velocity distribution in the presence and absence of hydrodynamic slip condition S f .
Symmetry 11 01171 g002
Figure 3. Effect of M on temperature distribution in the presence and absence of thermal slip condition S t .
Figure 3. Effect of M on temperature distribution in the presence and absence of thermal slip condition S t .
Symmetry 11 01171 g003
Figure 4. Effect of M on solute concentration in the presence and absence of solute slip condition S p .
Figure 4. Effect of M on solute concentration in the presence and absence of solute slip condition S p .
Symmetry 11 01171 g004
Figure 5. Effect of λ 1 on velocity distribution in the presence and absence of solute slip condition S f .
Figure 5. Effect of λ 1 on velocity distribution in the presence and absence of solute slip condition S f .
Symmetry 11 01171 g005
Figure 6. Effect of λ 2 on velocity distribution in the presence and absence of solute slip condition S f .
Figure 6. Effect of λ 2 on velocity distribution in the presence and absence of solute slip condition S f .
Symmetry 11 01171 g006
Figure 7. Effect of λ 3 on velocity distribution in the presence and absence of solute slip condition S f .
Figure 7. Effect of λ 3 on velocity distribution in the presence and absence of solute slip condition S f .
Symmetry 11 01171 g007
Figure 8. Effect of N r on temperature distribution in the presence and absence of thermal slip condition S t .
Figure 8. Effect of N r on temperature distribution in the presence and absence of thermal slip condition S t .
Symmetry 11 01171 g008
Figure 9. Effect of N b on nano-fluid volume fraction profile in the presence and absence of nano-fluid slip condition S g .
Figure 9. Effect of N b on nano-fluid volume fraction profile in the presence and absence of nano-fluid slip condition S g .
Symmetry 11 01171 g009
Figure 10. Effect of n on solute concentration in the presence and absence of solutal slip condition S p .
Figure 10. Effect of n on solute concentration in the presence and absence of solutal slip condition S p .
Symmetry 11 01171 g010
Figure 11. Effect of S c on solute concentration in the presence and absence of solute slip condition S p .
Figure 11. Effect of S c on solute concentration in the presence and absence of solute slip condition S p .
Symmetry 11 01171 g011
Figure 12. Effect of γ on solute concentration in the presence and absence of solute slip condition S p .
Figure 12. Effect of γ on solute concentration in the presence and absence of solute slip condition S p .
Symmetry 11 01171 g012
Figure 13. Effect of L e on nano-fluid volume fraction profile in the presence and absence of nano-fluid slip condition S g .
Figure 13. Effect of L e on nano-fluid volume fraction profile in the presence and absence of nano-fluid slip condition S g .
Symmetry 11 01171 g013
Figure 14. Effect of M, λ 2 and λ 1 in hydrodynamic slip condition S f on Skin friction coefficient.
Figure 14. Effect of M, λ 2 and λ 1 in hydrodynamic slip condition S f on Skin friction coefficient.
Symmetry 11 01171 g014
Figure 15. Effect of M, λ 2 , P r and N r in thermal slip condition S f on Nusselt number.
Figure 15. Effect of M, λ 2 , P r and N r in thermal slip condition S f on Nusselt number.
Symmetry 11 01171 g015
Table 1. Convergence results of FEM with the variation of number of elements ( P r = 0.773 , M = 2 , N b = 0.1 , N t = 0.1 , S c = 0.1 , R ( r ) = 0.1 , L e = 3 , N r = 1 n = 1 , L 1 = 0.4 , L 2 = 0.4 , L 3 = 0.4 , S f = 0.5 , S t = 0.5 , S p = 0.5 , S g = 0.5 ).
Table 1. Convergence results of FEM with the variation of number of elements ( P r = 0.773 , M = 2 , N b = 0.1 , N t = 0.1 , S c = 0.1 , R ( r ) = 0.1 , L e = 3 , N r = 1 n = 1 , L 1 = 0.4 , L 2 = 0.4 , L 3 = 0.4 , S f = 0.5 , S t = 0.5 , S p = 0.5 , S g = 0.5 ).
Number of Elements f ( 3 ) h ( 3 ) θ ( 3 ) S ( 3 ) ϕ ( 3 )
400.52750.06700.34160.28980.0546
800.53050.06700.34090.28980.0546
1200.53100.06700.34080.28980.0547
1600.53120.06700.34070.28980.0547
2000.53130.06700.34070.28980.0547
2400.53130.06700.34070.28980.0547
2800.53140.06700.34070.28980.0547
3200.53140.06700.34070.28980.0547
Table 2. Comparison of the FEM solution with the exact solution of Crane [31].
Table 2. Comparison of the FEM solution with the exact solution of Crane [31].
η Crane [31]FEM η Crane [31]FEM
01170.00090.0009
10.36790.367980.00030.0003
20.13530.135390.00010.0001
30.04980.0498100.00000.0000
40.01830.0183110.00000.0000
50.00670.0067120.00000.0000
60.00250.002513,14,150.00000.0000
Table 3. Comparison of skin friction coefficient for different values of M.
Table 3. Comparison of skin friction coefficient for different values of M.
MGireesha et al. [32]Mudassar et al. [33]
Exact Solution (a)
FEM (Our Results)
(b)
Error in %
| ( b a a ) | × 100
0.01.0001.0000001.00000800.00080
0.21.0951.0954451.09544580.00007
0.51.2241.2247451.22474460.00003
1.01.4141.4142141.41421320.00006
1.21.4831.4832401.48323930.00005
1.51.5811.5811391.58113840.00004
2.01.7321.7320511.73205040.00003
Table 4. Comparison of f ( 0 ) for various values of M when S f = 0 and P r when M = f w = S f = S t = λ 1 = λ 2 = R = 0 .
Table 4. Comparison of f ( 0 ) for various values of M when S f = 0 and P r when M = f w = S f = S t = λ 1 = λ 2 = R = 0 .
MMabood and Das [43]Fazle [29]FEM (Present) Pr Ali [44]Fazle [29]FEM (Present)
0−1.000008−1.0000084−1.0000082----
11.41421351.414213561.41421353----
52.44948972.449489742.449489630.720.80580.80880.8088
103.31662473.316624793.3166246310.96911.00001.0000
507.14142847.141428437.1414283931.91441.92371.9237
10010.04987510.049875610.0498751103.70063.72073.7207
50022.38302922.383029322.3830283----
100031.63858431.638584031.6385833----
Table 5. Computed values of skin fraction coefficient, heat transfer coefficient and mass transfer coefficient for N b = 0.1 , N t = 0.1 , = 0.1 , S c = 0.1 , σ = 0.1 , L e = 3 , S f = 0.5 , S t = 0.5 , S p = 0.5 , S g = 0.5 ).
Table 5. Computed values of skin fraction coefficient, heat transfer coefficient and mass transfer coefficient for N b = 0.1 , N t = 0.1 , = 0.1 , S c = 0.1 , σ = 0.1 , L e = 3 , S f = 0.5 , S t = 0.5 , S p = 0.5 , S g = 0.5 ).
M γ Pr Nr λ 1 λ 2 λ 3 f ( 0 ) θ ( 0 ) S ( 0 )
1 −0.51871906−0.39914438−0.33469636
1.50.10.77310.40.40.4−0.70007439−0.34330234−0.32039573
2 −0.86037797−0.29152259−0.30916534
0.1 −0.86037797−0.29152259−0.30916534
20.50.77310.40.40.4−0.87456687−0.27774645−0.53903314
1 −0.88233600−0.27260197−0.67751871
2 −0.86809394−0.46620439−0.30740431
20.1310.40.40.4−0.87169999−0.55675321−0.30668459
4 −0.87422111−0.62457605−0.30622757
2 −0.85775278−0.23671739−0.30981895
20.10.77330.40.40.4−0.85618413−0.20495708−0.31022032
4 −0.85513340−0.18408657−0.31049339
0.5 −0.84545297−0.29890014−0.31058477
1 −0.77396002−0.32928379−0.31696545
1.5 −0.70660726−0.35282310−0.32248200
0.5 −0.84570667−0.29845530−0.31051308
0.7 −0.81650214−0.31161198−0.31317081
0.9 −0.78749387−0.32389598−0.31577681
0.4−0.86037797−0.29152259−0.30916534
0.8−0.82636008−0.30157756−0.31127720
1.2−0.79384196−0.31065228−0.31323236

Share and Cite

MDPI and ACS Style

Khan, S.A.; Nie, Y.; Ali, B. Multiple Slip Effects on Magnetohydrodynamic Axisymmetric Buoyant Nanofluid Flow above a Stretching Sheet with Radiation and Chemical Reaction. Symmetry 2019, 11, 1171. https://doi.org/10.3390/sym11091171

AMA Style

Khan SA, Nie Y, Ali B. Multiple Slip Effects on Magnetohydrodynamic Axisymmetric Buoyant Nanofluid Flow above a Stretching Sheet with Radiation and Chemical Reaction. Symmetry. 2019; 11(9):1171. https://doi.org/10.3390/sym11091171

Chicago/Turabian Style

Khan, Shahid Ali, Yufeng Nie, and Bagh Ali. 2019. "Multiple Slip Effects on Magnetohydrodynamic Axisymmetric Buoyant Nanofluid Flow above a Stretching Sheet with Radiation and Chemical Reaction" Symmetry 11, no. 9: 1171. https://doi.org/10.3390/sym11091171

APA Style

Khan, S. A., Nie, Y., & Ali, B. (2019). Multiple Slip Effects on Magnetohydrodynamic Axisymmetric Buoyant Nanofluid Flow above a Stretching Sheet with Radiation and Chemical Reaction. Symmetry, 11(9), 1171. https://doi.org/10.3390/sym11091171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop