Next Article in Journal
Some Remarks on Classical and Classical-Quantum Sphere Packing Bounds: Rényi vs. Kullback–Leibler
Next Article in Special Issue
Investigation of Oriented Magnetic Field Effects on Entropy Generation in an Inclined Channel Filled with Ferrofluids
Previous Article in Journal
Fourier’s Law in a Generalized Piston Model
Previous Article in Special Issue
Natural Convection and Entropy Generation in a Square Cavity with Variable Temperature Side Walls Filled with a Nanofluid: Buongiorno’s Mathematical Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects on a Non-Linear Stretching Sheet

1
School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X01, Scottsville, 3209 Pietermaritzburg, South Africa
2
DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), Private Bag 3, Wits, 2050 Johannesburg, South Africa
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(7), 168; https://doi.org/10.3390/e19070168
Submission received: 22 February 2017 / Revised: 8 April 2017 / Accepted: 12 April 2017 / Published: 12 July 2017
(This article belongs to the Special Issue Entropy Generation in Nanofluid Flows)

Abstract

:
The entropy generation in unsteady three-dimensional axisymmetric magnetohydrodynamics (MHD) nanofluid flow over a non-linearly stretching sheet is investigated. The flow is subject to thermal radiation and a chemical reaction. The conservation equations are solved using the spectral quasi-linearization method. The novelty of the work is in the study of entropy generation in three-dimensional axisymmetric MHD nanofluid and the choice of the spectral quasi-linearization method as the solution method. The effects of Brownian motion and thermophoresis are also taken into account. The nanofluid particle volume fraction on the boundary is passively controlled. The results show that as the Hartmann number increases, both the Nusselt number and the Sherwood number decrease, whereas the skin friction increases. It is further shown that an increase in the thermal radiation parameter corresponds to a decrease in the Nusselt number. Moreover, entropy generation increases with respect to some physical parameters.

1. Introduction

The study of unsteady nanofluid flow, heat and mass transfers along a nonlinear stretching surface has received considerable attention during the last few years because of several applications in engineering processes, such as in materials manufacturing through extrusion, glass-fiber and paper production. Similarly, unsteady mixed convection in boundary layer flows have received attention with a large number of studies focusing on heat and mass transfer characteristics in nanofluids (Dessie et al. [1]). Nanofluids have increased thermal conductivity and convective heat transfer performance as compared to base fluids such as water and oils. The notion of a nanofluid was introduced by Choi [2] when he proposed the suspension of nanoparticles in a base fluid like water, oil or an ethylene-glycol mixture. These common base fluids have lower thermal conductivity, which is increased when nanoparticles are added. The increase in the thermal conductivity was explained by Buongiorno [3] in terms of the effect of particle Brownian motion and thermophoresis. Nanofluids have several applications due to the unique chemical and physical properties of the constituent nanoparticles. For instance, nanofluids have been used in applications that require high-performance cooling systems such as hot rolling, glass fiber production, rubber and the manufacture of metallic sheets [4].
Heat transfers due to free and mixed convection have several applications, for instance in electronic cooling, heat exchangers, etc. The study of axisymmetric magnetohydrodynamics (MHD) flow and heat transfer of power law fluid along an unsteady radially-stretching sheet was carried out by Ahmed et al. [5]. The conservation equations are solved analytically and numerically. Mohammadiun et al. [6] derived an exact solution of axisymmetric stagnation-point flow and heat transfer along a stationary infinite circular cylinder of a steady viscous compressible fluid with constant heat flux. The general self-similar solution was obtained with constant wall heat flux. The solutions of the system were obtained for different Reynolds numbers, compressibility factors and Prandtl numbers. Xiao et al. [7] proposed a mathematical model for heat convection in the presence of Brownian motion of nanoparticles where the physical properties of the fluid are predicted by the mean of fractal geometry. Thermal conductivity was found to have a positive and negative relation with concentration and the size of nanoparticles, respectively. The model output was found to be in good agreement with experimental data. Cai et al. [8] presented a review of research conducted on thermal conductivity and convective heat transfer in nanofluids by using fractal models and fractal-based techniques. Cai et al. [8] presented a model that considered the fractal distribution of nanoparticle sizes and heat convection between the fluid and the nanoparticles. Additionally, the heat transfer in nanofluids was modeled using three fractal-based models. The models were used to derive formulas for predicting the heat flux for boiling heat transfer. The predictions were found to be in a good agreement with the experimental data.
Shankar et al. [9] studied MHD flow, heat and mass transfer in a nanofluid flow along a stretching boundary with a non-uniform heat source/sink. The system of equations was solved using the Keller-box method. MHD effects in heat transfer have applications in science and technology. In MHD flow, the induced currents in the fluid generate forces, which in turn modify the flow field. Shahzad et al. [10] presented the exact solution for the MHD flow and heat transfer for a viscous incompressible fluid and a nonlinear sheet stretching radially in a porous medium.
Thermal radiation is important in the flow of a fluid, and consequently, the effects of thermal radiation on heat and mass transfer have been extensively studied. Ahmad et al. [11] considered the effect of thermal radiation on the steady MHD axisymmetric stagnation point flow of a viscous incompressible micropolar fluid along a shrinking sheet. The system of equations was solved numerically using an algorithm based on finite difference approximations. They found that the thickness of the thermal boundary layer becomes thinner as the thermal radiation parameter increased. The effect of melting heat transfer and second order slip with thermal radiation for a stagnation point flow was examined by Mabood et al. [12]. Singh et al. [13] analyzed unsteady MHD flow of a viscous incompressible fluid over a permeable stretching sheet and took into account the effect of thermal radiation.
In practical applications, mass transfer occurs due to molecular diffusion of species in homogeneous and heterogeneous chemical reactions (Sarada et al. [14]). The properties of the fluid can be affected by the diffusion of the species, which can either be generated or absorbed by the fluid. For this reason, the study of a chemical reaction in a fluid flow has attracted many researchers. Sarada et al. [14] analyzed the influence of a chemical reaction on unsteady MHD flow of a viscous incompressible fluid passing through an infinite vertical porous plate with varying suction.
The study of the unsteady flow of a viscous incompressible fluid along a linear stretching with a chemical reaction was investigated by Hunegnaw et al. [15]. They used a shooting technique and a fourth-order Runge–Kutta integration scheme combined with the Newton–Raphson method to solve the conservation equations. Barik [16] presented a study of the effects of a chemical reaction on unsteady rotating MHD flow in a porous medium with a heat source.
Entropy generation plays a vital role in the study of heat transfer processes. Entropy is a measure of the randomness or molecular disorder of a system. In accordance with the second law of thermodynamics, the entropy of a system always increases during an irreversible process and remains constant during a reversible process, that is entropy generation ( E g e n ) is always positive for an irreversible process and zero for reversible process. The performance of any engineering system is degraded by irreversibility, and entropy generation is a measure of the magnitude of the irreversibility of the process.
Entropy generation is disregarded in most of the studies reviewed earlier. In this work, the mechanisms for generating entropy are connected to heat transfer, fluid friction irreversibility, magnetic field and mass transfer. The pioneering work in the analysis of entropy generation was done by Bejan [17]. Subsequently, entropy generation in MHD Casson nanofluid flow in the proximity of a stagnation point was investigated by Qing et al. [18]. The findings suggested a positive correlation between entropy generation and an increase in the Brinkman number, Reynolds number, Hartmann number and porosity. The study of entropy generation with MHD flow peristaltic blood of nanofluid across a porous medium was considered by Rashidi et al. [19]. The system was solved numerically using the homotopy perturbation method (HPM). The entropy generation of MHD an Eyring–Powell nanofluid flow towards a permeable stretching surface in the presence of nonlinear thermal radiation was investigated by Bhatti et al. [20]. The successive linearization method (SLM) and Chebyshev spectral collocation method are used. Rashidi et al. [21] studied entropy generation in nanofluid flow along a permeable stretching surface near the stagnation point with heat generation/absorption and a convective boundary condition. Bhatti et al. [22] used the successive linearization method (SLM) to study the entropy generation on non-Newtonian Carreau nanofluid over a shrinking sheet. The thermal radiation and magnetohydrodynamics (MHD) are taken into account. The idea of entropy generation in nanofluids is a growing area of research, and recent studies have examined different source terms and flow geometries [23].
The spectral quasi-linearization method (SQLM) has not been previously used to solve equations for three-dimensional axisymmetric MHD nanofluid flow. Moreover, to the best of the authors’ knowledge, entropy generation in this type of fluid flow has not been previously studied. The aim of the study is to analyze thermo-diffusion effects in three-dimensional axisymmetric MHD nanofluid flow, heat and mass transfers over a nonlinearly circular stretching sheet with entropy generation. The flow is subjected to thermal radiation and a chemical reaction. The conservation equations are solved numerically using the spectral quasi-linearization method [24]. The SQLM combines fast convergence with accuracy. The method has been used in recent boundary layer flow and heat transfer studies, such as [24,25]. In addition, the nanofluid boundary condition suggested by Nield and Kuznetsov [26] is adopted where it is assumed that the nanoparticle mass flux at the wall vanishes. The results are found to be in good agreement with previously published work, such as Mustafa et al. [27].

2. Problem Formulation

Consider the unsteady three-dimensional MHD flow of an incompressible viscous flow of nanofluid. A cylindrical coordinate system ( r , θ , z ) is used. The velocity of the stretching sheet is assumed to be nonlinear along the radial direction. It is assumed that there is no nanoparticle flux across the wall, and the surface is stretching along the z-direction. In the ambient fluid, the temperature, solute concentration and nanoparticle concentration are denoted by T , C and ψ , respectively (see Figure 1). Here, C w , T w and ψ w represent the solute concentration, constant temperature and nanoparticle concentration respectively on the wall. The variable magnetic field intensity is denoted by B ( r , t ) where t represents time. The magnetic field acts in the positive z-direction normal to the sheet. In this study, B ( r , t ) generalizes the magnetic field term provided previously in [9,10] to:
B ( r , t ) = B 0 r ( n 1 ) / 2 1 λ t ,
where B 0 is the uniform magnetic field strength, n > 0 is the power-law index or stretching sheet parameter and λ represents the unsteadiness parameter.
The equations for the conservation of momentum, energy, mass and nanoparticle volume fraction, considering the boundary layer assumptions presented above, can be obtained [27,28] as:
u r + u r + w z = 0 ,
u t + u u r + w u z = ν f 2 u z 2 σ B 2 ( r , t ) u ρ f ,
T t + u T r + w T z = α f 2 T z 2 + τ D B ψ z T z + D T T T z 2 + D T C 2 C z 1 ρ f c p q r z ,
C t + u C r + w C z = D s 2 C z 2 + D C T 2 T z R ( r , t ) ( C C ) ,
ψ t + u ψ r + w ψ z = D B 2 ψ z 2 + D T T 2 T z 2 ,
where u and w represent the fluid velocity components in r and z directions, respectively; ν f , σ and ρ f are the kinematic viscosity, electrical conductivity and the density of the fluid, respectively; α f = k f / ( ρ c p ) f is the thermal diffusivity, c p is the heat capacity, D B is the Brownian diffusion coefficient, D T is the thermophoretic diffusion coefficient; τ = ( ρ c p ) s / ( ρ c p ) f the ratio of effective heat capacity of the nanoparticle material to heat capacity of the fluid; D C T and D T C are the Soret and Dufour diffusivities, D s is the solute diffusivity and R ( r , t ) is the chemical reaction.
The quantity q r is the relative heat flux. Applying the Rosseland approximation, [29], q r can be expressed as:
q r = 4 σ 3 K T 4 z ,
In Equation (7), σ is the Stefan–Boltzmann constant while K is the Rosseland mean spectral absorption coefficient. If the variation in temperature across the flow is small, T 4 may be expanded using the Taylor series about T , which yields:
T 4 4 T 3 T 3 T 4 .
Substituting Equation (8) in Equation (7), the relative heat flux becomes:
q r = 16 σ T 3 3 K T z .
To simplify the mass equation, the chemical reaction R ( r , t ) must be a constant. This condition holds if R ( r , t ) has the following form:
R ( r , t ) = R 0 a r n 1 1 λ t ,
where R 0 is constant and a > 0 is the stretching constant.
The boundary conditions considered (see [27]) are:
u = u w ( r ) = a r n 1 λ t , k f T z = h f ( T w T ) , D B ψ z + D T T T z = 0 , C = C w a t z = 0 , u 0 , T T , C C and ψ ψ as z ,
where k f = k 0 1 λ t is the thermal conductivity of the base fluid where k 0 is a constant. A similarity solution of the energy equation can be obtained if the Biot number,
B i = h f k 0 ν f a r ( n 1 ) ,
is a constant. This condition is satisfied if the heat transfer coefficient, h f , is proportional to r ( n 1 ) / 2 . Thus, the heat transfer coefficient is expressed as h f = b 0 r ( n 1 ) / 2 , where b 0 is a constant. The Biot number can be obtained as B i = h 0 ( ν f / a ) / k 0 .
Equations (2)–(6) are converted into coupled ordinary differential equations by using the following similarity variables (see [27]):
η = a ν f ( 1 λ t ) r ( n 1 ) / 2 z , θ ( η ) = T T T w T , S ( η ) = C C C w C , ϕ ( η ) = ψ ψ ψ .
The system of Equations (2)–(6) is transformed to:
f + ( n + 3 ) 2 f f n f 2 A f + 1 2 η f H a f = 0 ,
1 + N r 1 P r θ A 2 η θ + n + 3 2 f θ + N b ϕ θ + N t θ 2 + N d S = 0 ,
1 S c S A 2 η S + ( n + 3 ) 2 f S R 0 S + L d θ = 0 ,
ϕ A 2 S c η ϕ + n + 3 2 S c f ϕ + N t N b θ = 0 ,
subject to the boundary conditions:
f ( 0 ) = 0 , f ( 0 ) = 1 and f ( ) 0 , θ ( 0 ) = B i ( 1 θ ( 0 ) ) , θ ( ) 0 , S ( 0 ) = 1 , S ( ) 0 , N b ϕ ( 0 ) + N t θ ( 0 ) = 0 and ϕ ( ) 0 .
The parameters in Equations (14)–(17) are given by:
A = λ a r n 1 , H a = σ B 0 2 a ρ f , P r = ν f α , N r = 16 σ T 3 3 k f K , N b = τ D B ϕ ν f , N t = τ ( T w T ) D T ν f T , N d = D T C ( C w C ) ν f ( T w T ) , S c = ν f D B , R 0 = ( 1 λ t ) R ( x , t ) a r n 1 and L d = D C T ( T w T ) ν f ( C w C ) ,
where A is the unsteadiness parameter, H a is the Hartmann number, P r is the Prandtl number, N r is the thermal radiation parameter, N b is the Brownian motion parameter, N t is the thermophoresis parameter, N d is the modified Dufour parameter, S c is the Schmidt number, R 0 is a constant and L d is the modified Soret parameter.
The skin friction coefficient C f , the Nusselt number N u r and the Sherwood number S h r can be expressed as:
(i)
Skin-fraction coefficient:
C f R e 1 / 2 = f ( 0 ) ,
where R e = u w ( r ) r ν f is the Reynolds number.
(ii)
Nusselt number:
N u = r h w k f ( T w T ) ,
where h w represents surface heat flux, which can be obtained by:
h w = k f + 16 σ T 3 3 K T z z = 0 .
Equation (21) becomes:
1 + N r θ ( 0 ) = R e 1 / 2 N u r
(iii)
The Sherwood number for solute concentration equation is:
S h r = r h m D s ( C w C ) ,
where:
h m = D s C z z = 0 ,
where h m is the surface mass flux. Equation (24) can be written as:
S ( 0 ) = R e r ( 1 / 2 ) S h r
The mass flux is zero due to the nanoparticle boundary condition used, and for this reason, it is not possible to define the Sherwood number for nanoparticle concentration at the wall.

3. Entropy Generation Analysis

Entropy generation suggests a wastage of energy; thus, minimization of entropy production is often a major objective. Entropy generation analysis can be used as an effective tool for the identification of causes of inefficiency in any system and offers scope for the improvement in the design of any device or process. The limitation of global energy resources provides an impetus for re-examining energy production systems and consumption patterns (Arikoglu et al. [30]). From a theoretical perspective, the second law of thermodynamics is utilized to study energy producing, converting and consuming systems. The volumetric rate of local entropy generation for a nanofluid along with thermal radiation and a magnetic field can be expressed as (see [18,22,23]:
E g e n = k T ˜ 2 T ˜ z 2 + 16 σ T ˜ 3 3 K T ˜ z 2 1 st   part + μ T ˜ u ˜ z 2 2 nd   part + σ B 0 2 T ˜ u ˜ 2 3 rd   part + R D B ψ ψ z 2 + R D B T ˜ T z ψ z 4 th   part
In Equation (27), the entropy generation is stated in four parts. The first part is entropy generation due to heat transfer irreversibility; the second part is the entropy generation due to the viscous dissipation irreversibility; the third part is the entropy generation due to the applied magnetic field irreversibility; and the fourth part is due to the diffusive irreversibility.
We define the entropy generation number as the ratio of the local volumetric entropy generation rate E g e n to a characteristic rate of entropy generation E 0 , that is,
N G = E g e n E 0 ,
where:
E 0 = k f ( T w T ) 2 T 2 r 2 .
Using Equations (27)–(29), the non-dimensional entropy generation can be represented as:
N G = R e 1 + N r θ 2 ( η ) + B r R e Ω f 2 + B r H a 2 Ω f 2 + R e Σ Ω 2 ϕ 2 ( η ) + R e Σ Ω θ ( η ) ϕ ( η ) ,
where R e , N r , B r , Ω , H a , Σ are the Reynolds number, the thermal radiation parameter, Brinkman number, dimensionless temperature difference, Hartmann number and the diffusive parameter, respectively. These parameters are expressed by the following relations:
B r = μ u w 2 ( r ) k f Δ T , Ω = Δ T T = T w T T , Σ = R D B ψ k f .

4. Method of Solution

The quasi-linearization method (QLM) is a generalization of the Newton–Raphson method (Bellman and Kalaba [31]). The derivation of the QLM is based on the linearization of the nonlinear components of the governing equations using the Taylor series assuming that the difference between the value of the unknown function is negligible between the current iteration, r + 1 , and the previous iteration, r.
Applying the quasi-linearization scheme to Equations (14)–(17) with the boundary conditions, Equation (18), yields the following iterative schemes:
a 0 , r f r + 1 + a 1 , r f r + 1 + a 2 , r f r + 1 + a 3 , r f r + 1 = R f ,
b 0 , r θ r + 1 + b 1 , r θ r + 1 + b 2 , r θ r + 1 + b 3 , r f r + 1 + b 4 , r S r + 1 + b 5 , r ϕ r + 1 = R θ ,
c 0 , r S r + 1 + c 1 , r S r + 1 + c 2 , r S r + 1 + c 3 , r f r + 1 + c 4 , r θ r + 1 = R S ,
d 0 , r ϕ r + 1 + d 1 , r ϕ r + 1 + d 2 , r f r + 1 + d 3 , r θ r + 1 = R ϕ ,
subject to
f r + 1 ( 0 ) = 0 , f r + 1 ( 0 ) = 1 , f r + 1 ( ) 0 , θ r + 1 ( 0 ) = B i ( 1 θ r + 1 ( 0 ) ) , θ r + 1 ( ) 0 , S r + 1 ( 0 ) = 1 , S r + 1 ( ) 0 , N b ϕ r + 1 ( 0 ) + N t θ r + 1 ( 0 ) = 0 , and ϕ r + 1 ( ) 0 ,
where the coefficients in Equations (32)–(35) are obtained as:
a 0 , r = 1 , a 1 , r = n + 3 2 f r A 2 η ,
a 2 , r = 2 n f r A H a , a 3 , r = n + 3 2 f r , b 0 , r = 1 + N r 1 P r , b 1 , r = n + 3 2 f r A 2 η + N b ϕ r + 2 N t θ r , b 2 , r = 0 , b 3 , r = n + 3 2 θ r ,
b 4 , r = N d , b 5 , r = N b θ r , c 0 , r = 1 S c , c 1 , r = n + 3 2 f r A 2 η , c 2 , r = R 0 ,
c 3 , r = n + 3 2 S r , c 4 , r = L d , d 0 , r = 1 , d 1 , r = n + 3 2 S c f r A 2 S c η ,
d 2 , r = n + 3 2 S c ϕ r , d 3 , r = N t N b .
A Chebyshev pseudo-spectral method [24] was used to solve Equations (32)–(35). The Chebyshev interpolating polynomial defined by Equation (41) is used with Gauss–Lobatto points [32,33] to define the unknown functions where:
x i = cos π i N , i = 0 , 1 , . . . , N ; 1 x i 1 .
The variable N in Equation (41) is the number of collocation points. A truncated domain [ 0 , L ] is used to approximate the semi-infinite domain to facilitate computations. The parameter L represents the boundary condition at infinity. In order to model the behavior of the flow at infinity, the parameter L should be a large number. The domain [ 0 , L ] is transformed into [ 1 , 1 ] using the linear transformation η = L ( x + 1 ) 2 .
The spectral collocation method is used to construct a differentiation matrix to approximate the derivative of unknown variables at the collocation points as a matrix vector product:
d F r ( 1 ) d η ( η j ) = k = 0 N D j k f ( η k ) = D F m , j = 0 , 1 , 2 , . . . , N ,
where D = 2 D / L and F = f ( η 0 ) , f ( η 1 ) , f ( η 1 ) , . . . , f ( η N ) T represent the vector function at the collocation points. The high order derivatives are given as powers of D, such as:
F r ( p ) = D p F r ,
where p is the order of the derivative. Spectral collocation is applied at r using the differentiation matrix D in order to approximate derivatives of unknown functions in Equations(32)–(35) with Equation (36), which yields:
A 1 , 1 f + A 1 , 2 θ + A 1 , 3 S + A 1 , 4 ϕ = R f ,
A 1 , 2 f + A 2 , 2 θ + A 2 , 3 S + A 2 , 4 ϕ = R θ ,
A 1 , 3 f + A 2 , 3 θ + A 3 , 3 S + A 3 , 4 ϕ = R S ,
A 1 , 4 f + A 2 , 4 θ + A 3 , 4 S + A 4 , 4 ϕ = R ϕ .
Here:
A 1 , 1 = a 0 , r D 3 + diag a 1 , r D 2 + diag a 2 , r D + diag a 3 , r I , A 1 , 2 = 0 , A 1 , 3 = 0 , A 1 , 4 = 0 ,
A 2 , 1 = diag b 3 , 1 I , A 2 , 2 = diag b 0 , r D 2 + diag b 1 , r D + diag b 2 , r I , A 2 , 3 = diag b 4 , r D 2 , A 2 , 4 = diag b 5 , r D ,
A 3 , 1 = diag c 3 , r I , A 3 , 2 = diag c 4 , r D 2 , A 3 , 3 = c 0 , r D 2 + diag c 1 , r D + c 2 , r I and A 3 , 4 = 0 ,
A 4 , 1 = diag d 2 , r I , A 4 , 2 = diag d 3 , r D 2 ,
A 4 , 3 = 0 , A 4 , 4 = diag d 0 , r D + diag d 1 , r D ,
R f = n + 3 2 f r f r n f r 2 , R ϕ = n + 3 2 S c f r ϕ r ,
R s = n + 3 2 f r S r a n d R ϕ = n + 3 2 S c f r ϕ r ,
where diag ( ) represents diagonal matrices of order ( N + 1 ) × ( N + 1 ) , I is an ( N + 1 ) × ( N + 1 ) identity matrix and f , θ , S and ϕ are the values of functions f , θ , S and ϕ , respectively. Equations (44)–(47) are solved as a matrix system using the SQLM scheme where the iteration is started with initial approximate solutions obtained as:
f 0 ( η ) = 1 exp ( η ) , θ 0 ( η ) = B i 1 + B i exp ( η ) , S 0 ( η ) = exp ( η ) and ϕ 0 ( η ) = N t N b B i 1 + B i exp ( η ) .
The above equations can be expressed in matrix form as follows:
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 F r + 1 Θ r + 1 S r + 1 Φ r + 1 = R f R θ R s R ϕ

5. Results and Discussion

The study investigated entropy generation in an unsteady three-dimensional MHD nanofluid along a nonlinear stretching sheet and considered the influence of thermal radiation and a chemical reaction. The conservation equations are solved numerically using the spectral quasi-linearization method (SQLM). The SQLM has been used in a limited number of studies to solve boundary layer flow, heat and mass transfer problems [34]. The physical parameters in this paper are mostly chosen from the literature, for example in the papers [21,26,27,35].
A comparison with previously published results is shown in Table 1 when A = 0 , H a = 0 , Nr = 0 , N d = 0 and L d = 0 (i.e., in the absence of unsteadiness parameter, Hartmann number, the thermal radiation number, Dufour parameter and Soret parameter, respectively). These results are comparable to those of Mustafa et al. [27], which validates this numerical method.
To demonstrate further the convergence of the SQLM, Figure 2, Figure 3, Figure 4 and Figure 5 show the effect of the number of collocation points on the accuracy of the solutions obtained for the velocity, temperature, solute concentration and nanoparticle concentration profiles, respectively.
It was found that an increase in the number of collocation points results in a reduction of residual error. However, after a certain point (the point at which an optimal residual error is obtained), an increase in the number of collocation points does not have a significant effect on residual error.
In Figure 2, it is observed that the optimal residual error for the velocity profile is around 10 9 , which is achieved when the number of the collocation points is within the range 30 and 40; beyond that range, the accuracy starts declining. From Figure 3 and Figure 5, it is noted that the optimal residual errors for the temperature and nanoparticle concentration profiles, respectively, are found around 10 12 when the number of collocation points is within 100 and 120. Figure 4 shows that the optimal residual of the solute concentration profile occurs approximately at 10 11 when the number of collocation points is between 80 and 100. It is interesting to note that there is no magic number of collocation points that induces the optimal residual errors for all of the physical properties. Additionally, there is variability in the behavior of the decline of the accuracy after the range at which the optimal residuals are obtained; for instance, in Figure 2, the accuracy of the solution of the velocity profile declines much faster as compared to the accuracy of the solution of temperature, solute concentration and nanoparticle concentration profiles respectively shown in Figure 3, Figure 4 and Figure 5. Overall, these results demonstrate the accuracy and convergence of the SQLM technique.
Table 2 shows the computed skin friction, heat transfer and the mass transfer coefficients, which are represented by f ( 0 ) , θ ( 0 ) and S ( 0 ) , respectively, for various values of n , A , H a , N r and P r . It is observed that the skin friction decreases as n increases, whereas the heat transfer coefficient and mass transfer increase with increasing n with other parameters fixed. It is also noted that the skin friction reduces as A increases while the heat and mass transfer rates decrease. The skin friction, heat and mass transfer rates vary inversely with H a and appear to be independent of changes in N r , while the heat transfer rate decreases with the increase in N r . The same result also holds for the skin friction and the mass transfer with respect to P r . The heat transfer increases when P r increases.
In Table 3 it is noted that for small values of the Biot number, small changes in the Biot number correspond to a large changes in the maximum temperature. In contrast, with high values of the Biot number, a large change in the Biot number does not appear to affect the maximum temperature significantly.
The parameter n plays an important role in fluid flow. A fluid is called pseudoplastic fluid when n < 1 and dilatant fluid when n > 1 . It is a Newtonian fluid when n = 1 . A fluid is shear-thinning or shear-thickening according as n < 1 or n > 1 . The effect of n on the fluid velocity, temperature and solute concentration profiles, respectively, is shown in Figure 6, Figure 7 and Figure 8. Figure 6 illustrates the velocity profile ( f ) with the variation in the stretching sheet parameter keeping other parameters fixed. The magnitude of the radial component of the velocity decreases with increasing n. In contrast, n has a significant impact on η and the velocity components u and w, respectively. Figure 7 shows that an increases in n leads to a decrease in the thickness of the thermal boundary layer. As a result, the temperature decreases with increasing n and an increase in the heat transfer rate from the sheet. Increasing n increases the deformation rate from the wall to the fluid. Figure 8 shows the effect of n on the solute concentration when the other parameters take fixed values. It is noticed that increases in the non-linear stretching lead to a decrease in the solutal concentration.
Figure 9 and Figure 10 describe the effects of unsteadiness parameter A on the temperature and the solute concentration profiles, respectively. Two cases have been studied, namely A < 0 and A > 0 . From Figure 9, it is observed that the temperature increases with increasing values of A > 0 . A similar effect is observed for A < 0 . That the temperature decreases, as moving far away from the stretching sheet in the dynamics region due to the increase of the unsteadiness parameter. Physically, the decrease of temperature can be attributed to heat lostfrom the sheet as a result of the increase in the unsteadiness. Similar properties of solute concentration profiles are observed in Figure 10.
Figure 11 and Figure 12 demonstrate the behavior of several values of the Hartmann number on the velocity and temperature profiles when the other physical parameters are fixed. Figure 11 shows that the dimensionless radial component of the velocity decreases as the Hartmann number increases. This is an indication of a Lorentz force slowing the motion of the fluid in the radial direction. Consequently, the boundary layer thickness increases. The physical implication is that the momentum of the fluid flow may be controlled through an applied magnetic field. Figure 12 shows that the temperature increases with an increase in the Hartmann number due to the effect of the transverse magnetic field in the fluid. It is also worth noting that the application of the magnetic field affects the thermal boundary layer thickness positively in the sense that the thickness of the thermal boundary increases in the presence of a magnetic field.
Figure 13 shows how the temperature profiles vary with the thermal radiation parameter. The temperature increases with increased thermal radiation, which is in line with physical observations regarding the impact of increasing thermal radiation.
Figure 14 presents the influence of the Prandtl number on the temperature profiles. As the Prandtl number increases, the temperature decreases. Smaller Prandtl numbers suggest a fluid with higher thermal conductivities; hence, heat diffuses more rapidly from the heated surface than in the case of fluids with higher Prandtl numbers.
Figure 15 illustrates the behavior of the temperature profiles with different values of the thermophoresis parameter, N t , when the other parameter are fixed. It is observed that the temperature increases as the thermophoresis parameter increases because the thermal boundary layer increases with increases in the thermophoresis parameter.
Figure 16 demonstrates the influence of the Schmidt number on the solute concentration profile with fixing the other parameters. It is noticed that the solute concentration decreases as the Schmidt number increases.
Figure 17 shows the relation between the Biot number B i and the temperature profiles. The Biot number can be perceived as the ratio of internal (conductive) resistance of solid to external (convective resistance). Here, it is found that temperature increases with the increase in the Biot number. An increase in the Biot number causes a stronger convection, and this leads to higher surface temperatures. Further, for sufficiently large values of the Biot number, the temperature approaches its maximum value.
Figure 18 shows that high values of the stretching sheet parameter n lead to a reduction in the thickness of the nanoparticle concentration boundary layer thickness up to a certain value of η . Beyond this critical distance, the opposite trend is observed. In Figure 19, we observe a similar pattern of nanoparticle concentration profiles as in Figure 9 and Figure 10. In Figure 20, the nanoparticle concentration profiles increase with the thermophoresis parameter. Similar results are shown in Figure 21 for the Brownian motion parameter. In Figure 22, it is observed that increasing the Schmidt number reduces the nanoparticle concentration profile. The Schmidt number S c has an opposite effect on the Brownian diffusion coefficient. When D B is small, the penetration depth of the nanoparticle concentration profiles becomes shorter. Thus, as the Schmidt number increases, the penetration depth becomes shorter. Moreover, the increase in S c results in a decrease in the concentration rate, as well as the heat transfer rate.
It is interesting to note that Figure 23, Figure 24, Figure 25, Figure 26 and Figure 27 suggest different ways of controlling the entropy generation. Figure 23 shows the effect of the Reynolds number on the entropy generation N G . It is observed that entropy generation increases with the increase in the Reynolds number. This observation could be explained by the increase in random eddies, vortices and flow fluctuation in the fluid, which are associated with increasing R e . As a consequence, the heat transfer increases, which, in turn, enhances the entropy generation. The entropy generation number increases in the proximity of the sheet due to a decrease in the fluid friction. Reducing the Reynolds number induces a reduction in the entropy generation number.
Figure 24 shows entropy generation with the Brinkman number. The Brinkman number is a measure of the importance of heat produced by viscous heating relative to heat transported by molecular conduction. An increase in the Brinkman number causes an increase in the entropy generation in the vicinity of the sheet. Heat produced by viscous dissipation dominates over heat transported by molecular conduction in the proximity of the sheet. In the vicinity of the sheet, significant heat generation occurs within the layers of the moving fluid particles, which in turn increases the entropy generation number by enhancing the degree of disorder of the system. This gradually fades with the distance from the sheet.
The influence of the Hartmann number on the entropy generation number is given in Figure 25. In close proximity to the sheet, an increase in the Hartmann number corresponds to a noticeable increase in the entropy generation, whereas far from the sheet, the increase in the Hartmann number has rather negligible impact on entropy generation. This change in entropy generation number in response to the Hartmann number is related to the increase in the resistance of the fluid motion due to an increase in the Hartmann number. This variation in entropy generation with the Hartmann number can be attributed to higher heat transfer rates, which result in an increase of the entropy generation number. However, far from the sheet, the effect of the Hartmann number is insignificant, causing a corresponding insignificant entropy generation.
Figure 26 displays the entropy generation with the Brinkman parameter ( B r Ω 1 ). An increase in B r Ω 1 results in a remarkable increase in the entropy generation number in the area surrounding the sheet. The increase in the entropy generation number around the sheet due to the increase in B r Ω 1 can be related to a decrease in the fluid friction caused by increasing B r Ω 1 . Moreover, Equation (31) suggests that the increase in B r Ω 1 results in an increase in the sheet velocity, which affects the fluid in the surroundings of the sheet and that explains, in part, the noticeable effect of B r Ω 1 on the entropy generation proximity of the sheet and the negligible effect in the further region. Figure 27 demonstrates the relation between the Biot number and the entropy generation number. The entropy generation number increases prominently in the proximity of the sheet as the Biot number increases. On the contrary, the Biot number has an insignificant effect on the entropy generation number far away from the sheet.
Figure 28 demonstrates how the skin friction f ( 0 ) varies according to different values of the magnetic field parameter H a and the stretching sheet parameter n. It is observed that f ( 0 ) increases with increasing values of H a and n.
The variation in the Nusselt number θ ( 0 ) due to the thermophoresis number N t and the stretching sheet parameter n is graphed in Figure 29. It is noticed that with increasing values of N t and n, θ ( 0 ) decreases.
Figure 30 displays the change in the Sherwood number S ( 0 ) caused by the effect of the Schmidt number S c and the stretching sheet parameter n. It is observed that S ( 0 ) increases as S c and n increase.
The changing of the Nusselt number θ ( 0 ) due to the effect of the Prandtl number P r , and the thermal radiation parameter N r is exhibited in Figure 31. The figure depicts that the θ ( 0 ) is an increasing function of P r and N r .

6. Conclusions

A mathematical formulation and analysis of the entropy generation rate in an unsteady three-dimensional axisymmetric MHD nanofluid flow over a non-linear stretching sheet subject to thermal radiation and a chemical reaction has been presented. The influence of different physical parameters on the entropy generation number has been demonstrated and discussed in detail. From the discussion, the following outcomes may be inferred:
  • The heat transfer rate increases with increasing sheet stretching.
  • An increase in the Reynolds number and the Brinkman number corresponds to a significant increase in the entropy generation number. Therefore, it can be ascertained that the entropy generation number is highly affected by viscous dissipation when the nanofluid flow has a large Reynolds number.
  • An increase in the Biot and Hartmann numbers corresponds to a significant increase in the entropy generation number in the vicinity of the sheet surface. The significance of the Biot and Hartmann numbers gradually fades with distance from the sheet.
  • The entropy generation rate can be minimized by controlling the physical parameters.
  • The number of collocation points has a significant influence on the accuracy of the solutions.

Acknowledgments

The authors would like to thank the University of KwaZulu-Natal, DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS) and Claude Leon Foundation, South Africa, for the financial support. Furthermore, the authors are thankful to the reviewers for their valuable comments and suggestions, which have improved the quality of this paper.

Author Contributions

Sabyasachi Mondal proposed the main idea and design framework. Mohammed Almakki carried out the numerical simulations and the first draft of this paper. Sharadia Dey and Precious Sibanda revised the draft, done some write-up and checked the calculations. All authors contributed equally in the analysis of the results.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
r , θ , z cylindrical polar coordinate axes
B 0 magnetic field strength (kg·s 2 A 1 )
uvelocity in radial direction (m·s 1 )
n , λ constants
wvelocity in axial direction (m·s 1 )
ν f kinematic viscosity (m 2 s 1 )
Ttemperature variable (K)
σ electrical conductivity ( Ω 1 m 1 )
T w temperature of fluid at sheet (K)
ρ f density of fluid (kg·m 3 )
T ambient temperature of fluid (K)
α f thermal diffusion of fluid (m 2 s 1 )
Csolute concentration of fluid (kg·m 3 )
C W solute concentration at wall (kg·m 3 )
C solute concentration far away from the disk (kg·m 3 )
ψ nanoparticle concentration
ψ W nanoparticle concentration at wall
ψ nanoparticle concentration far away from the disk
τ ratio of nanoparticle heat capacity
D B Brownian mention coefficient (kg·m 1 s 1 )
D T Thermophoretic diffusion coefficient (kg·m 1 s 1 K 1 )
D T C Dufour diffusion coefficient
c p specific heat (m 2 s 2 K 1 )
k f thermal conductivity (W·m 1 K 1 )
q r radiative heat flux (kg·m 2 )
D s solute diffusion coefficient
D C T Soret diffusion coefficient
Rchemical reaction parameter
( ρ c p ) s heat capacity of the nanoparticle
( ρ c p ) f heat capacity of fluid
σ Stefan–Boltzmann constant
K mean absorption coefficient
R 0 , a constants
k 0 , b 0 constants
h f heat transfer coefficient (W·m 2 K 1 )
B i Biot number
η dimensionless variable
θ dimensionless temperature
Sdimensionless solute concentration
ϕ dimensionless nanoparticle concentration
fdimensionless velocity
Aunsteadiness parameter
H a Hartmann number
N r thermal radiation parameter
P r Prandtl number
N b Brownian mention parameter
N t Thermophoresis parameter
N d Dufour parameter
S c Schmidt number
L d Soret parameter
C f skin friction coefficient
N u r Nusselt number
S h r Sherwood number
u w ( r ) velocity of the stretching sheet
R e Reynolds number
h w heat flux (W·m 2 )
h m mass flux (kg·m 2 s 1 )
E g e n volumetric entropy generation per unit length (W·m 3 K 1 )
E 0 dimensionless entropy generation rate
B r Brinkman number
Ω dimensionless parameter
H a Hartmann number
Σ diffusive constant parameter
Δ T difference between ( T W T )(K)

References

  1. Dessie, H.; Kishan, N. Unsteady MHD flow of heat and mass transfer of nanofluids over stretching sheet with a non-uniform heat/source/sink considering viscous dissipation and chemical reaction. Int. J. Eng. Res. Afr. 2015, 14. [Google Scholar] [CrossRef]
  2. Chol, S. Enhancing thermal conductivity of fluids with nanoparticles. ASME Publ. FED 1995, 231, 99–106. [Google Scholar]
  3. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  4. Ahmed, N.; Goswami, J.; Barua, D. Effects of chemical reaction and radiation on an unsteady MHD flow past an accelerated infinite vertical plate with variable temperature and mass transfer. Indian J. Pure Appl. Math. 2013, 44, 443–466. [Google Scholar] [CrossRef]
  5. 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]
  6. Mohammadiun, H.; Rahimi, A.; Kianifar, A. Axisymmetric stagnation-point flow and heat transfer of a viscous, compressible fluid on a cylinder with constant heat flux. Sci. Iran. 2013, 20, 185–194. [Google Scholar] [CrossRef]
  7. Xiao, B.; Yang, Y.; Chen, L. Developing a novel form of thermal conductivity of nanofluids with Brownian motion effect by means of fractal geometry. Powder Technol. 2013, 239, 409–414. [Google Scholar] [CrossRef]
  8. Cai, J.; Hu, X.; Xiao, B.; Zhou, Y.; Wei, W. Recent developments on fractal-based approaches to nanofluids and nanoparticle aggregation. Int. J. Heat Mass Transf. 2017, 105, 623–637. [Google Scholar] [CrossRef]
  9. Shankar, B.; Yirga, Y. Unsteady heat and mass transfer in MHD flow of nanofluids over stretching sheet with a non-uniform heat source/sink. Int. J. Math. Comput. Sci. Eng. 2013, 7, 1267–1275. [Google Scholar]
  10. Shahzad, A.; Ali, R.; Khan, M. On the exact solution for axisymmetric flow and heat transfer over a nonlinear radially stretching sheet. Chin. Phys. Lett. 2012, 29, 084705. [Google Scholar] [CrossRef]
  11. Ahmad, S.; Ashraf, M.; Syed, K. Effects of thermal radiation on MHD axisymmetric stagnation point flow and heat transfer of a micropolar fluid over a shrinking sheet. World Appl. Sci. J. 2011, 15, 835–848. [Google Scholar]
  12. Mabood, F.; Shafiq, A.; Hayat, T.; Abelman, S. Radiation effects on stagnation point flow with melting heat transfer and second order slip. Results Phys. 2017, 7, 31–42. [Google Scholar] [CrossRef]
  13. Prasad, K.; Vaidya, H.; Vajravelu, K.; Datti, P.; Umesh, V. Axisymmetric mixed convective MHD flow over a slender cylinder in the presence of chemically reaction. Int. J. Appl. Mech. Eng. 2016, 21, 121–141. [Google Scholar] [CrossRef]
  14. Sarada, K.; Shanker, B. The effect of chemical reaction on an unsteady MHD free convection flow past an infinite vertical porous plate with variable suction. Int. J. Eng. Mod. Res. 2013, 3, 725–735. [Google Scholar]
  15. Hunegnaw, A.; Kishan, N. Unsteady MHD heat and mass transfer flow over stretching sheet in porous medium with variable properties considering viscous dissipation and chemical reaction. Am. Chem. Sci. J. 2014, 4, 901–917. [Google Scholar] [CrossRef]
  16. Barik, R. Heat and Mass Transfer Effects on Unsteady MHD Flow through an Accelerated Isothermal Vertical Plate Embedded in Porous Medium in the Presence of Heat Source and Chemical Reaction. Eur. J. Adv. Eng. Technol. 2016, 3, 56–61. [Google Scholar]
  17. Bejan, A. Entropy Generation Minimization: The Method of Thermodynamic Optimization of Finite-Size Systems and Finite-Time Processes; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  18. Qing, J.; Bhatti, M.M.; Abbas, M.A.; Rashidi, M.M.; Ali, M.E.S. Entropy generation on MHD Casson nanofluid flow over a porous stretching/shrinking surface. Entropy 2016, 18, 123. [Google Scholar] [CrossRef]
  19. Rashidi, M.M.; Bhatti, M.M.; Abbas, M.A.; Ali, M.E.S. Entropy generation on MHD blood flow of nanofluid due to peristaltic waves. Entropy 2016, 18, 117. [Google Scholar] [CrossRef]
  20. Bhatti, M.M.; Abbas, T.; Rashidi, M.M.; Ali, M.E.S.; Yang, Z. Entropy generation on MHD eyring—Powell nanofluid through a permeable stretching surface. Entropy 2016, 18, 224. [Google Scholar] [CrossRef]
  21. Rashidi, M.; Mohammadi, F.; Abbasbandy, S.; Alhuthali, M. Entropy generation analysis for stagnation point flow in a porous medium over a permeable stretching surface. J. Appl. Fluid Mech. 2015, 8, 753–765. [Google Scholar] [CrossRef]
  22. Bhatti, M.M.; Abbas, T.; Rashidi, M.M.; Ali, M.E.S. Numerical simulation of entropy generation with thermal radiation on MHD Carreau nanofluid towards a shrinking sheet. Entropy 2016, 18, 200. [Google Scholar] [CrossRef]
  23. Freidoonimehr, N.; Rahimi, A.B. Comment on “Effects of thermophoresis and Brownian motion on nanofluid heat transfer and entropy generation” by M. Mahmoodi, Sh. Kandelousi, Journal of Molecular Liquids, 211 (2015) 15–24. J. Mol. Liq. 2016, 216, 99–102. [Google Scholar] [CrossRef]
  24. Agbaje, T.; Mondal, S.; Makukula, Z.; Motsa, S.; Sibanda, P. A new numerical approach to MHD stagnation point flow and heat transfer towards a stretching sheet. Ain Shams Eng. J. 2016, in press. [Google Scholar] [CrossRef]
  25. Kaladhar, K.; Motsa, S.; Srinivasacharya, D. Mixed Convection Flow of Couple Stress Fluid in a Vertical Channel with Radiation and Soret Effects. J. Appl. Fluid Mech. 2016, 9, 43–50. [Google Scholar]
  26. Nield, D.; Kuznetsov, A. The onset of convection in a horizontal nanofluid layer of finite depth: A revised model. Int. J. Heat Mass Transf. 2014, 77, 915–918. [Google Scholar] [CrossRef]
  27. 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]
  28. Sajid, M.; Hayat, T.; Asghar, S.; Vajravelu, K. Analytic solution for axisymmetric flow over a nonlinearly stretching sheet. Arch. Appl. Mech. 2008, 78, 127–134. [Google Scholar] [CrossRef]
  29. 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]
  30. Arikoglu, A.; Ozkol, I.; Komurgoz, G. Effect of slip on entropy generation in a single rotating disk in MHD flow. Appl. Energy 2008, 85, 1225–1236. [Google Scholar] [CrossRef]
  31. Bellman, R.E.; Kalaba, R.E. Quasilinearization and Nonlinear Boundary-Value Problems; RAND Corporation: New York, NY, USA, 1965. [Google Scholar]
  32. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.M.; Thomas, A., Jr. Spectral Methods in Fluid Dynamics; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  33. Trefethen, L.N. Spectral Methods in MATLAB; SIAM: Oxford, UK, 2000. [Google Scholar]
  34. Motsa, S.; 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. [Google Scholar] [CrossRef]
  35. Rashidi, M.; Abelman, S.; Mehr, N.F. Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. Int. J. Heat Mass Transf. 2013, 62, 515–525. [Google Scholar] [CrossRef]
Figure 1. Physical configuration and coordinate system.
Figure 1. Physical configuration and coordinate system.
Entropy 19 00168 g001
Figure 2. Effects of the number of collocation points number on the residual error of the velocity profile | | R e s ( f ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 2. Effects of the number of collocation points number on the residual error of the velocity profile | | R e s ( f ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g002
Figure 3. Effects of the number of collocation points number on the residual error of the temperature profile | | R e s ( θ ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 3. Effects of the number of collocation points number on the residual error of the temperature profile | | R e s ( θ ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g003
Figure 4. Effects of the number of collocation points number on the residual error of the solute concentration profile | | R e s ( S ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 4. Effects of the number of collocation points number on the residual error of the solute concentration profile | | R e s ( S ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g004
Figure 5. Effects of the number of collocation points number on the residual error of the nanoparticle concentration profile | | R e s ( ϕ ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 5. Effects of the number of collocation points number on the residual error of the nanoparticle concentration profile | | R e s ( ϕ ) | | when n = 3 , A = 0 . 3 , P r = 7 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g005
Figure 6. Effects of the stretching parameter n on f ( η ) when A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 6. Effects of the stretching parameter n on f ( η ) when A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g006
Figure 7. Effects of the stretching parameter n on θ ( η ) where A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 10 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 7. Effects of the stretching parameter n on θ ( η ) where A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 10 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g007
Figure 8. Effects of the stretching parameter n on S ( η ) for A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 8. Effects of the stretching parameter n on S ( η ) for A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g008
Figure 9. Effects of the unsteadiness parameter (A) on θ ( η ) for n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 9. Effects of the unsteadiness parameter (A) on θ ( η ) for n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g009
Figure 10. Effects of the unsteadiness parameter (A) on S ( η ) when n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 10. Effects of the unsteadiness parameter (A) on S ( η ) when n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g010
Figure 11. Effects of the Hartmann number ( H a ) on f ( η ) for n = 3 , A = 0 . 3 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 11. Effects of the Hartmann number ( H a ) on f ( η ) for n = 3 , A = 0 . 3 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g011
Figure 12. Effects of the Hartmann number ( H a ) on θ ( η ) when n = 3 , A = 0 . 3 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 12. Effects of the Hartmann number ( H a ) on θ ( η ) when n = 3 , A = 0 . 3 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g012
Figure 13. Effects of the thermal radiation parameter ( N r ) on θ ( η ) where n = 3 , A = 0 . 3 , H a = 0 . 5 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 13. Effects of the thermal radiation parameter ( N r ) on θ ( η ) where n = 3 , A = 0 . 3 , H a = 0 . 5 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g013
Figure 14. Effects of the Prandtl number ( P r ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , β i = 0 . 2 .
Figure 14. Effects of the Prandtl number ( P r ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , β i = 0 . 2 .
Entropy 19 00168 g014
Figure 15. Effects of the thermophoresis parameter ( N t ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 15. Effects of the thermophoresis parameter ( N t ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g015
Figure 16. Effects of the Schmidt number ( S c ) on S ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 16. Effects of the Schmidt number ( S c ) on S ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g016
Figure 17. Effects of the Biot number ( B i ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 .
Figure 17. Effects of the Biot number ( B i ) on θ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 .
Entropy 19 00168 g017
Figure 18. Effects of the stretching parameter (n) on ϕ ( η ) when A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 18. Effects of the stretching parameter (n) on ϕ ( η ) when A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g018
Figure 19. Effects of the unsteadiness parameter (A) on ϕ ( η ) for n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 19. Effects of the unsteadiness parameter (A) on ϕ ( η ) for n = 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g019
Figure 20. Effects of the thermophoresis parameter ( N t ) on ϕ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 20. Effects of the thermophoresis parameter ( N t ) on ϕ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g020
Figure 21. Effects of the Brownian motion parameter ( N b ) on ϕ ( η ) when n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Figure 21. Effects of the Brownian motion parameter ( N b ) on ϕ ( η ) when n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N t = 0 . 5 , N d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , L d = 0 . 02 , B i = 0 . 2 .
Entropy 19 00168 g021
Figure 22. Effects of the Schmidt number ( S c ) on ϕ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 22. Effects of the Schmidt number ( S c ) on ϕ ( η ) for n = 3 , A = 0 . 3 , H a = 0 . 5 , N r = 0 . 2 , P r = 7 , N b = 0 . 5 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g022
Figure 23. Effects of the Reynolds number R e on the entropy generation number N G when N r = 0 . 2 , B r = 1 , Ω = 1 , H a = 0 . 5 and Σ = 0 . 5 .
Figure 23. Effects of the Reynolds number R e on the entropy generation number N G when N r = 0 . 2 , B r = 1 , Ω = 1 , H a = 0 . 5 and Σ = 0 . 5 .
Entropy 19 00168 g023
Figure 24. Effects of Brinkman number B r on the entropy generation number N G when R e = 2 , N r = 0 . 2 , Ω = 1 , H a = 0 . 5 and Σ = 0 . 5 .
Figure 24. Effects of Brinkman number B r on the entropy generation number N G when R e = 2 , N r = 0 . 2 , Ω = 1 , H a = 0 . 5 and Σ = 0 . 5 .
Entropy 19 00168 g024
Figure 25. Effects of Hartmann number H a on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 , Ω = 1 and Σ = 0 . 5 .
Figure 25. Effects of Hartmann number H a on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 , Ω = 1 and Σ = 0 . 5 .
Entropy 19 00168 g025
Figure 26. Effects of the Brinkman group parameter B r Ω 1 on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 and Σ = 0 . 5 .
Figure 26. Effects of the Brinkman group parameter B r Ω 1 on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 and Σ = 0 . 5 .
Entropy 19 00168 g026
Figure 27. Effects of Biot number B i on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 , H a = 0 . 5 , Ω = 1 and Σ = 0 . 5 .
Figure 27. Effects of Biot number B i on the entropy generation number N G when R e = 2 , N r = 0 . 2 , B r = 1 , H a = 0 . 5 , Ω = 1 and Σ = 0 . 5 .
Entropy 19 00168 g027
Figure 28. Effect of the magnetic field H a and the stretching sheet parameter n on the skin friction f ( 0 ) when A = 0 . 3 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 28. Effect of the magnetic field H a and the stretching sheet parameter n on the skin friction f ( 0 ) when A = 0 . 3 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g028
Figure 29. Effect of the thermophoresis number ( N t ) and the stretching sheet parameter (n) on the Nusselt number ( θ ( 0 ) ) when A = 0 . 3 , H a = 0 . 5 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 29. Effect of the thermophoresis number ( N t ) and the stretching sheet parameter (n) on the Nusselt number ( θ ( 0 ) ) when A = 0 . 3 , H a = 0 . 5 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g029
Figure 30. Effect of the Schmidt number S c and the stretching sheet parameter n on the Sherwood number S ( 0 ) when A = 0 . 3 , H a = 0 . 5 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 30. Effect of the Schmidt number S c and the stretching sheet parameter n on the Sherwood number S ( 0 ) when A = 0 . 3 , H a = 0 . 5 , P r = 7 , N r = 0 . 2 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g030
Figure 31. Effect of the Prandtl number R e and the thermal radiation parameter N r on the Nusselt number θ ( 0 ) when n = 3 , A = 0 . 3 , H a = 0 . 5 , S c = 7 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Figure 31. Effect of the Prandtl number R e and the thermal radiation parameter N r on the Nusselt number θ ( 0 ) when n = 3 , A = 0 . 3 , H a = 0 . 5 , S c = 7 , N b = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , N t = 0 . 5 , R 0 = 0 . 3 , B i = 0 . 2 .
Entropy 19 00168 g031
Table 1. Current Nusselt number θ ( 0 ) compared with Mustafa et al. [27].
Table 1. Current Nusselt number θ ( 0 ) compared with Mustafa et al. [27].
n N t S c P r Mustafa et al. [27]Present Results
0.50.1205 1 . 9112911 1 . 91068095
0.5 1 . 2170065 1 . 21659065
0.7 0 . 9815765 0 . 98122822
1.00.555 1 . 6914582 1 . 69104675
10 1 . 4740787 1 . 47375172
20 1 . 2861370 1 . 28590965
2.50.5200.7 0 . 6619164 0 . 66986678
5 1 . 4784288 1 . 47847763
7 1 . 5758736 1 . 57604858
Table 2. Skin friction coefficient, heat transfer coefficient and mass transfer coefficient for N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
Table 2. Skin friction coefficient, heat transfer coefficient and mass transfer coefficient for N b = 0 . 5 , N t = 0 . 5 , N d = 0 . 02 , L d = 0 . 02 , S c = 7 , R 0 = 0 . 3 , B i = 0 . 2 .
nA H a N r P r f ( 0 ) θ ( 0 ) S ( 0 )
1 1 . 43922866 0 . 16417529 2 . 89107257
20.30.50.27 1 . 70034047 0 . 16661752 3 . 16939593
4 2 . 12897896 0 . 16969131 3 . 66539915
−0.5 1 . 77701012 0 . 17085800 3 . 72686482
300.50.27 1 . 87092985 0 . 16944779 3 . 54329048
0.5 1 . 96353732 0 . 16749470 3 . 34434681
1.5 2 . 17151786 0 . 16755905 3 . 37606776
30.32.50.27 2 . 39114443 0 . 16678981 3 . 33136800
5 2 . 86708560 0 . 16495200 3 . 23526081
1 1 . 92657387 0 . 16615306 3 . 42533294
30.30.51.57 1 . 92657387 0 . 16464481 3 . 42514164
2 1 . 92657387 0 . 16310964 3 . 42504767
4 1 . 92657387 0 . 16585371 3 . 42553409
30.30.50.25 1 . 92657387 0 . 16708464 3 . 42554096
9 1 . 92657387 0 . 16894896 3 . 42668462
Table 3. Effects of the Biot number on the maximum temperature.
Table 3. Effects of the Biot number on the maximum temperature.
Biot NumberMaximum TemperatureChange in Maximum Temperature
0.10.1264
0.90.3401169.07
20.521853.43
50.737741.38
100.85315.63
500.967813.46

Share and Cite

MDPI and ACS Style

Almakki, M.; Dey, S.; Mondal, S.; Sibanda, P. On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects on a Non-Linear Stretching Sheet. Entropy 2017, 19, 168. https://doi.org/10.3390/e19070168

AMA Style

Almakki M, Dey S, Mondal S, Sibanda P. On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects on a Non-Linear Stretching Sheet. Entropy. 2017; 19(7):168. https://doi.org/10.3390/e19070168

Chicago/Turabian Style

Almakki, Mohammed, Sharadia Dey, Sabyasachi Mondal, and Precious Sibanda. 2017. "On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects on a Non-Linear Stretching Sheet" Entropy 19, no. 7: 168. https://doi.org/10.3390/e19070168

APA Style

Almakki, M., Dey, S., Mondal, S., & Sibanda, P. (2017). On Unsteady Three-Dimensional Axisymmetric MHD Nanofluid Flow with Entropy Generation and Thermo-Diffusion Effects on a Non-Linear Stretching Sheet. Entropy, 19(7), 168. https://doi.org/10.3390/e19070168

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