Next Article in Journal
Semantic Image Segmentation with Deep Convolutional Neural Networks and Quick Shift
Previous Article in Journal
Use of Determination of the Importance of Criteria in Business-Friendly Certification of Cities as Sustainable Local Economic Development Planning Tool
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of Entropy Generation with Multi-Slip Effects in MHD Unsteady Flow of Viscous Fluid Past an Exponentially Stretching Surface

1
College of Applied Science, Beijing University of Technology, Beijing 100124, China
2
Department of Mathematics, Govt postgraduate Gordon College, Rawalpindi 46000, Pakistan
3
Department of Mathematics, Comsats University, Islamabad 44000, Pakistan
4
Department of Mathematics, Islamabad Campus, Preston University, Kohat 26000, Pakistan
*
Author to whom correspondence should be addressed.
Submission received: 9 December 2019 / Revised: 23 January 2020 / Accepted: 3 February 2020 / Published: 6 March 2020

Abstract

:
The current study aims to probe the impacts of entropy in a hydromagnetic unsteady slip flow of viscous fluid past an exponentially stretching sheet. Appurtenant similarity variables are employed to transmute the governing partial differential equations into a system of non-linear differential equations, which are analytically solved by utilizing the homotopy analysis method (HAM). Moreover, a shooting technique with fourth–fifth order Runge–Kutta method is deployed to numerically solve the problem. The impact of the physical parameters that influence the flow and heat transmission phenomena are sketched, tabulated and discussed briefly. Additionally, the impact of these parameters on entropy generation is thoroughly discussed by plotting graphs of the local entropy generation number and the Bejan number.

1. Introduction

Initially, the analytical solution for a boundary layer flow problem over a linearly stretching sheet was explored by Crane [1]. Since then, numerous studies have been carried out on this topic [2,3,4,5]. Such investigations have many diverse applications in crystal growth, polymer extrusion, the spinning of fibers, condensation process, metallic sheet cooling, etc. However, the stretching rate of the surface may not always be linear in practical cases. There may be cases where the sheet may stretch with an exponential order. Motivated by this idea, Magyari and Keller [6] investigated the flow and heat transmission of a Newtonian fluid over an exponentially stretching sheet with an exponential temperature distribution. Elbashbeshy [7] discussed the characteristics of the heat transmission of viscous fluid over a permeable exponentially stretching sheet and attained similarity solutions. Al-Odat et al. [8] examined the impacts of magnetic field on the thermal boundary layer of a viscous fluid past an exponentially stretching surface. El-Aziz [9] assumed the impacts of viscous dissipation on a mixed convective micropolar fluid flow that was induced by an exponentially stretching surface. Bidin and Nazar [10] applied the implicit finite difference scheme to obtain to similarity solutions of radiative boundary layer flow of a Newtonian fluid over an exponentially stretching surface. Ishak [11] numerically investigated the steady hydro magnetic flow of a viscous fluid past a stretching surface in the presence of thermal radiation and found that thermal radiation plays a vital role in the augmentation of fluid temperature. Mukhopadhyay and Gorla [12] assumed the velocity and thermal slip conditions at the surface of a exponentially stretching sheet and explored the flow and heat transmission characteristics. They concluded that velocity and thermal slip parameters have decreasing impact on temperature. Sahoo and Poncet [13] explored the influence of partial slip on the flow and heat transmission of third-grade fluid past a continuously stretching surface. Bhattacharyya and Vajravelu [14] examined the stagnation point flow of a viscous fluid past an exponentially contracting sheet and attained similar solutions. Mukhopadhyay et al. [15] studied a Casson fluid flow that was induced by a stretching of sheet with an exponential order. Jyothi [16] examined the effects of thermophoresis on a dissipative Magneto-hydrodynamics viscous fluid flow due to an exponentially stretching surface and observed that thermophoresis had a decreasing effect on a species concentration profile. Reddy [17] considered a hydromagnetic Casson fluid flow past an inclined exponential surface in the presence of thermal radiation and a chemical reaction, and they carried out a comprehensive parametric analysis of the problem. Patil et al. [18] obtained non-similar solutions of a mixed convective flow of a Newtonian fluid past a stretching surface in an exponential order. Jusoh et al. [19] analyzed the MHD rotating flow of a ferrofluid past a permeable, exponentially stretching/shrinking sheet and obtained dual solutions.
Flow and heat transmission procedures are subjected to energy losses that are irreversible. According to Bejan [20], these energy losses can be measured by using the quantity known as the entropy generation. He remarked that by recognizing the elements that produce entropy, energy losses can be reduced. This idea was adopted by engineers and investigators to explicate entropy impacts in various geometrical configurations. However, Tayamol [21] was the first to examine the entropy generation impacts in a viscous flow past a stretching sheet that was immersed in permeable media. Aiboud and Saouli [22] explored the entropy production in a hydromagnetic viscoelastic fluid past a linearly stretching sheet. Butt et al. [23] investigated the entropy generation effects in a mixed convective flow of a second-grade fluid past a stretching surface and found that the viscoelastic parameter has a significant impact on entropy production. Later on, Butt et al. [24] considered the slip effects on a stretchable surface and studied the entropy effects in the magneto hydrodynamic flow of a viscous fluid. They remarked that the presence of a slip at the surface reduces the production of entropy. The effects of slip and heat generation/absorption effects on a nanofluid flow past a linearly stretching sheet were examined by Noghrehabadi et al. [25]. Butt and Ali [26] carried out a study related to entropy effects in a viscous fluid flow past an unsteady stretching surface. Rashidi et al. [27] considered the stagnation point flow of a Newtonian fluid past a porous stretching sheet that was immersed in spongy media and analyzed the entropy impacts. The entropy effects of a three-dimensional Newtonian fluid over an exponential sheet were considered by Afridi and Qasim [28]. A mathematical investigation of an asymmetrical wavy motion of blood under entropy generation effects was presented by Riaz et al. [29]. Zaib, A. et al. explored the aligned magnetic flow comprising of nanoliquid over a radially stretching surface with entropy generation [30]. Recently, Butt et al. [31] theoretically explored the entropy impacts of a Casson nanofluid flow past an unsteady stretching surface. The found that unsteadiness had an enhancing effect on the generation of entropy effects.
Motivated by the stated facts, efforts have been devoted to scrutinize the impacts of entropy generation in the magnetohydrodynamic unsteady slip flow of a viscous fluid passing over an exponentially stretching sheet. To the best of our knowledge, entropy impacts, with the effect logs of an MHD unsteady slip flow have not been presented yet, and the results that are communicated here are new. The homotopy analysis method (HAM) and the shooting technique with the fourth–fifth order Runge–Kutta method were employed to solve the considered problem. The impact of the pertinent parameters on flow and heat and mass transmission features as well as on entropy generation are presented via graphs and tables, and they are explained briefly.

2. Mathematical Modeling

Consider the unsteady, two-dimensional incompressible flow of an electrically conducting viscous fluid due to an exponentially stretching surface as shown in Figure 1. It is assumed that the exponentially stretching surface is placed along the x-axis, and the y-axis is situated normal to the surface. The viscous fluid is considered to be confined in the area y 0 . A variable magnetic field B ( x , t ) = B 0 1 c t e x 2 L is applied normal to the surface where B 0 is the constant. Here, c is the dimensional constant and L represents the characteristics’ length. The effects of the induced magnetic field are neglected as the value of the magnetic Reynolds number is small. The surface is stretched with the exponential velocity U w ( x , t ) = U 0 1 c t e x L , and the surface temperature is kept at T w ( x ) = T + T o 1 c t e x 2 L . The ambient fluid temperature is assumed to be \ T . The effects of viscous and joule dissipations are assumed to be present. Under the application of the boundary layer approximation, the governing equations for continuity, momentum and heat are:
u x + v y = 0 ,
u t + u u x + v u y = ν 2 u y 2 σ B 2 ρ u ,
T t + u T x + v T y = k ρ c p 2 T y 2 + ν c p ( u y ) 2 + σ B 2 ρ c p u 2 ,
The appurtenant boundary constraints are
u =   u w ( x , t ) = U 0 1 c t e x L + α u y , v = 0 , T   = T w ( x , t ) = T + T 0 1 c t e A o x 2 L   at   y = 0 , u 0 , v 0 , T T   as   y . }
Here, ( u , v ) are the velocity components in the ( x , y ) directions, respectively, T is the temperature of the fluid, T 0 is the reference temperature, T is the temperature far away from the surface, ρ represents the fluid density, ν is the kinematic viscosity, k is the thermal conductivity, σ is the electrical conductivity, c p represents the specific heat at constant pressure, U o is the reference velocity, α is the hydrodynamic slip parameter, T o denotes the reference temperature, A 0 is the dimensional constant temperature, and c is the constant with c t < 1 .
The following set of similarity transformation is used to non-dimensionalized Equations (1)–(3):
u = U 0 1 c t e x L f ( ξ ) , v = U 0 ν 2 L ( 1 c t ) e x 2 L [ f ( ξ ) + ξ f ( ξ ) ] , ξ   = U 0 ν 2 ν ( 1 c t ) e x 2 L y ,         θ ( ξ ) = T T T w T . }
By substituting Equation (5) into Equations (1)–(3), the continuity Equation (1) is identically satisfied, and Equations (2) and (3) take the form:
f + f f 2 f 2 A ( ξ f + 2 f ) M f = 0 ,
θ + Pr f θ A 0 Pr f θ A Pr ( 4 θ + ξ θ ) + Pr E c f 2 + M Pr E c f 2 = 0 .
The associated boundary conditions after applying Equation (5) take the form
f ( 0 ) = 0 , f ( 0 ) = 1 + γ f ( 0 ) , θ ( 0 ) = 1 ,   f ( ξ ) 0 , θ ( ξ ) 0    as    ξ . }
Here, M = 2 σ B 0 2 L ρ U 0 represents the magnetic field parameter, A = c L e X U 0 is the unsteadiness parameter, γ = α U 0 e X 2 ν L ( 1 c t ) is the slip parameter, X = x L is the dimensionless parameter, Pr = μ c p k is the Prandtl number, and E c = U w 2 c p ( T w T ) is the Eckert number. The skin friction coefficient C f x and the local Nusselt number N u x are defined as:
C f x = τ w | y = 0 ρ u w 2 ,     N u x = x q w | y = 0 k ( T w T ) ,
where the shear stress τ w and the heat flux q w at the surface are given as:
τ w = μ ( u y ) | y = 0 ,     q w = k ( T y ) | y = 0 .
By substituting Equation (10) into Equation (9) and by using Equation (5), the dimensionless form of the skin friction coefficient C f x and the local Nusselt number N u x are:
Re x 1 / 2 C f x = f ( 0 ) ,     2 X Re x 1 / 2 N u x = θ ( 0 ) ,
where Re x = U w x ν is the local Reynolds number.

3. Entropy Analysis

After the application of the boundary layer approximation, the expression for the local volumetric rate of entropy generation S G for the viscous fluid in the presence of a magnetic field is defined as:
S G = k T ( T y ) 2 Entropy   due   to   Heat   Transfer   + ( u y ) 2   Entropy   due   to   viscous   dissipation + σ B o 2 T u 2 .   Entropy   due   to   magnetic   field
The above expression shows that the key factors that are responsible for the production of entropy effects are heat transfer, viscous dissipation, and the presence of a magnetic field.
The dimensionless form of the expression for the local entropy generation defined in Equation (12) can be obtained by using Equation (5) as:
N G = S G S G o = Re L θ 2 + Re L B r Ω f 2 + Re L B r Ω M   f 2 ,
where N G is the ratio of the local entropy generation rate to the characteristic entropy generation rate and is called the local entropy generation number. Here, S G o = k ( T w T ) 2 T 2 L 2 is the characteristic entropy generation rate, Ω = T w T T w denotes the dimensionless temperature difference, Re L = U L L v is the local Reynolds number, and B r = Pr E c is the Brinkman number. Equation (13) can also be written in the form:
N G = N G H + N G F + N G M ,
where
N G H = Re L θ 2 , N G F = Re L B r Ω f 2 , N G M = Re L B r Ω M   f 2 . }
Here, N G H denotes the entropy generation number due to heat transfer, N G F represents the entropy generation number due to viscous dissipation, and N G M is the entropy generation number due to the presence of the magnetic field.
Another useful quantity is the Bejan number, which is given as:
B e = Entropy   due   to   heat   transfer Local   entropy   generation = N G H N G H + N G F + N G M .
This expression helps to identify whether the heat transfer entropy effects dominate the viscous dissipation and magnetic field or vice versa. It can be seen that the value of the Bejan number varies between 0 and 1. When the value of B e is greater than 0.5, the entropy effects that are caused by heat transfer dominate the fluid friction and magnetic field effects. On the other hand, if the value of the Bejan number is less than 0.5, the entropy effects due to fluid friction and magnetic field dominate the transfer entropy effects. When B e = 0.5 , the contributions of both the entropy effects are equal.

4. Solution to the Problem

4.1. Homotopy Analysis Method

The analytical method known as the homotopy analysis method (HAM) is a powerful technique that is helpful in solving highly nonlinear equations. This method has been employed to obtain the solution of the problems related to fluid dynamics and heat transmission (see [23,24,26,27]). In this article, Equations (6) and (7) with the boundary conditions of Equation (8) were solved by the HAM by considering the following linear operators and linear guesses:
f 0 ( ξ ) = 1 1 + γ ( 1 e ξ ) , θ 0 ( ξ ) = e ξ .   }
f ( f ) = d 3 f d ξ 3 d f d ξ ,   θ ( θ ) = d 2 θ d ξ 2 θ . }
The above operator possesses the following properties:
f ( C 1 + C 2 e ξ + C 3 e ξ ) = 0 , θ ( C 4 e ξ + C 5 e ξ ) = 0 , }
where C j s ,   ( j = 1 , 2 , ...5 ) are arbitrary constants.
The zeroth order deformation problem can be written as:
( 1 q ) f [ F ¯ ( ξ ; q ) f 0 ( ξ ) ] = f q N f [ F ¯ ( ξ ; q ) ] ,
( 1 q ) θ [ Θ ¯ ( ξ ; q ) θ 0 ( ξ ) ] = θ q N θ [ F ¯ ( ξ ; q ) , Θ ¯ ( ξ ; q ) ] ,
F ¯ ( 0 ; q ) = 0 , F ¯ ( 0 ; q ) ξ = 1 + γ 2 F ¯ ξ 2 ( 0 ; q ) , Θ ¯ ( 0 ; q ) = 1 , F ¯ ( ; q ) ξ = 0 , Θ ¯ ( ; q ) = 0 . }
Here, f   and   θ are the convergence control parameters and q [ 0 , 1 ] denotes the embedding parameter that is related to the deformation mappings F ¯ ( ξ ; q ) ,   Θ ¯ ( ξ ; q ) that continuously deform from f 0 ( ξ ) ,   θ 0 ( ξ ) to f ( ξ ) ,   θ ( ξ ) when q varies from 0 to 1.
The non-linear operators are defined as:
N f [ F ¯ ( ξ ; q ) ] = 3 F ¯ ( ξ ; q ) ξ 3 + F ¯ ( ξ ; q ) 2 F ¯ ( ξ ; q ) ξ 2 2 ( F ¯ ( ξ ; q ) ξ ) 2 A ( ξ 2 F ¯ ( ξ ; q ) ξ 2 + 2 F ¯ ( ξ ; q ) ξ ) M F ¯ ( ξ ; q ) ξ , }
N θ [ F ¯ ( ξ ; q ) , Θ ¯ ( ξ ; q ) ] = 2 Θ ¯ ( ξ ; q ) ξ 2 + Pr F ¯ ( ξ ; q ) Θ ¯ ( ξ ; q ) ξ A o Pr F ¯ ( ξ ; q ) ξ Θ ¯ ( ξ ; q ) A Pr ( 4 Θ ¯ ( ξ ; q ) + ξ Θ ¯ ( ξ ; q ) ξ ) + Pr E c ( 2 F ¯ ( ξ ; q ) ξ 2 ) 2 + M Pr E c ( F ¯ ( ξ ; q ) ξ ) 2 , }
The rest of the details of the method can be seen in the literature (see [23,24,26]). The solution of the differential Equations (6) and (7) with the boundary conditions of Equation (8) can be written in the form of infinite series as:
f ( ξ ) = f 0 ( ξ ) + m = 1 f m ( ξ ) , θ ( ξ ) = θ 0 ( ξ ) + m = 1 θ m ( ξ ) . }
From the above procedure, it can be noted that the series solutions that are mentioned in Equation (25) contain the convergence control parameters f   and   θ , which can be adjusted to control the convergence rate of the series solutions. In order to ensure the fast convergence of these series solutions, optimal values of f   and   θ were used. In order to obtain the optimal values of f   and   θ , the averaged root mean square errors were calculated as follows:
E m , f = 1 N + 1 j = 0 k [ N f ( i = 0 m f i ( j Δ ζ ) ) ] 2 , E m , θ = 1 N + 1 j = 0 k [ N f ( i = 0 m f i ( j Δ ζ ) , i = 0 m θ i ( j Δ ζ ) ) ] 2 , }
where f i ( j Δ ζ )   and   θ i ( j Δ ζ ) are the discretization of the continuous functions f ( ξ )   and   θ ( ξ ) into N pieces by selecting Δ ζ = 0.2 and N = 20 . The optimal values of f   and   θ were computed by minimizing the values of E m , f and E m , θ by employing the first derivative law of calculus. Table 1 displays the optimal values of f   and   θ with their corresponding averaged root mean square errors E m , f and E m , θ at different orders of approximation. It was noticed that the values of f   and   θ improved with the simultaneous decrease in E m , f and E m , θ as the order of approximation increased. The convergence rate was further accelerated by employing a Padé approximation on the series solutions that were attained by the homotopy method. Table 2 shows the convergence of the numerical values of f ( 0 )   and   θ ( 0 ) for the [ m / m ] homotopy Padé approximation when f = 0.314784 ,   θ = 0.398461 , M = 1.0 ,   γ = 0.5 ,   A = 1.0 , A 0 = 0.5 , Pr = 1.0 ,   E c = 0.2 . It can be seen that the convergence of the series solutions of Equation (25) was attained up to 5 decimal places at the [ 15 / 15 ] homotopy Padé approximation.

4.2. Numerical Method

In order to counter check the results that were computed by the HAM, the nonlinear differential Equations (6) and (7) were numerically solved by applying shooting technique with fourth–fifth order Runge–Kutta–Fehlberg method. For this, Equations (6) and (7) were transformed in to a set of first order differential equations. For this purpose, let f = y 1 ,   f = y 2 ,   f = y 3 ,   θ = y 4 ,   θ = y 5 . Thus,
y 1 = y 2 , y 2 = y 3 , y 3 = ( y 1 y 3 + 2 y 2 2 + A ( ξ y 3 + 2 y 2 ) + M y 2 ) , y 4 = y 5 , y 5 = Pr y 1 y 5 + A o Pr y 2 y 4 + A Pr ( 4 y 4 + ξ y 5 ) Pr E c y 3 2 M Pr E c y 2 2 . }
and the associated initial conditions are
y 1 ( 0 ) = 0 , y 2 ( 0 ) = 1 + γ r 1 , y 3 ( 0 ) = r 1 , y 4 ( 0 ) = 1 , y 5 ( 0 ) = r 2 , }
where r 1   and   r 2 are the unknown values. In order to solve the above system of Equations (27) and (28) the values of r 1   and   r 2 were to be determined. In the beginning, these values were guessed, and the solution of the system of Equations (27) and (28) was achieved by applying the fourth–fifth-order Runge–Kutta–Fehlberg method. The Secant method is then employed to improve the guessed values of r 1   and   r 2 and the procedure was iteratively repeated. To obtain the solution of the system, the semi-infinite domain [ 0 , [ was truncated to a finite domain, [ 0 , ξ [ where ξ was so chosen that the boundary conditions defined in Equations (8) were asymptotically satisfied by the numerical solution. The value of ξ was set at 20 for all values of parameters, and the calculations were done. The convergence criterion was taken to be 10 5 , and a step size of 0.001 was selected.

5. Results and Discussion

The solutions of non-linear Ordinary differential Equations (6) and (7) with the boundary constraints of Equation (8) were compared with existing literature under limiting conditions. Table 3 was drawn to compare the numerical values of θ ( 0 ) with those mentioned by Magyari and Keller [6], El-Aziz [9], Nazar [10] and Ishak [11] for several values of Pr and M when γ = 0.0 , A = 0.0 , A o = 1.0 , E c = 0.0 . It was witnessed that the numerical values of θ ( 0 ) that were obtained by the HAM and the shooting method in the current study were in excellent agreement with the literature. Table 4 was constructed to note the impacts of the parameters γ , M and A on the magnitude of the dimensionless skin friction coefficient f ( 0 ) . This illustrates that rise in the values of M and A enhanced the skin friction at the surface, whereas the slip parameter had a decreasing impact on f ( 0 ) . Table 5 shows the values of θ ( 0 ) that were obtained by both the HAM and the shooting technique for various values of the different parameters γ , M , A , Pr , E c and temperature exponent A 0 . An increase in the value of θ ( 0 ) was noticed with the rise in the slip parameter γ , the unsteadiness parameter A , the temperature exponent A 0 , and the Prandtl number Pr . On the other hand, increasing the values of the parameters M and E c had a decreasing impact on θ ( 0 ) .
Figure 2a–c explicates the impact of the parameter M , γ and A on the velocity profile f ( ξ ) . As the value of the magnetic field parameter M increased, the Lorentz force became strong, thus resisting the fluid movement and causing a decline in velocity f ( ξ ) . Additionally, for fixed M , the velocity f ( ξ ) achieved a maximal value at the exponentially stretching surface and reduced with rise in distance ξ . In the far away regime, the velocity was asymptotically zero. The impacts of the slip parameter γ on fluid velocity f ( ξ ) are sketched in Figure 2b. A decline in velocity was observed with an increase in the value of the slip parameter γ , which was quite compatible with the physical fact. In Figure 2c, the velocity profile f ( ξ ) is drawn versus ξ for variation in the unsteadiness parameter A . It is quite evident from the graph that increasing the unsteadiness parameter A results in decreasing the velocity f ( ξ ) . Thus, the unsteadiness parameter A had a decreasing impact on the momentum boundary layer thickness.
Figure 3 and Figure 4 are plotted to demonstrate the impacts of various flow parameters on temperature distribution θ ( ξ ) . Figure 3a describes the impact of variation in the magnetic field parameter M on θ ( ξ ) . The resistive force becomes strong with an increase in the value of M , which ultimately results in a rise in the temperature θ ( ξ ) . The temperature distribution θ ( ξ ) is plotted in Figure 3b for various values of the slip parameter γ . It can be seen that the fluid temperature was augmented with increasing values of γ . Figure 3c illustrates the behavior of the fluid temperature by varying the values of the unsteadiness parameter A . A significant reduction in the temperature was seen by increasing the values of A . The impact of the temperature exponent A o on θ ( ξ ) is presented in Figure 4a. A minute boost in the fluid temperature was seen as the value of A o increased. However, this boost was very miniature as compared to the variation in temperature for other flow parameters. Figure 4b presents the effects of Pr on the temperature distribution θ ( ξ ) . A decline in the thermal boundary layer thickness was seen with a rising Pr value. On the other hand, Figure 4c shows that with the rise in values of E c , the viscous dissipation effects were augmented, which caused an increase in the fluid temperature.
As indicated in Equation (15), there are three sources of entropy generation, i.e., entropy production due to heat transfer N G H , entropy due to viscous dissipation N G F , and the entropy effects caused by an applied magnetic field N G M . Figure 5 shows the contribution of each type of entropy source within the boundary layer thickness. It is quite evident from the figure that entropy production was maximum at the surface of the exponentially stretching sheet. Furthermore, the contribution of the heat transfer entropy effects was maximum as compared to the entropy impacts due to viscous dissipation and the magnetic field. These entropy effects were quite significant within the boundary layer region. However, in the far away regime, these effects were insignificant.
Figure 6 and Figure 7 present the impacts of several flow parameters on the local entropy generation number N G . In Figure 6a, the graph of N G is plotted against ξ to illustrate the effects of M . It was noticed from the figure that in the absence of a magnetic field, the entropy effects were not so significant. However, as the value of M started to increase, N G increased near the exponentially stretching surface. Thus, the presence of a magnetic field plays a notable role in the production of entropy effects. Figure 6b elucidates the impacts of the slip parameter γ on N G . As the value of γ increased, the value of N G reduced, which indicates that the slip phenomena at the surface reduced the impacts of entropy. The influence of A on N G is illustrated in Figure 6c. A boost in the value of N G was noticed near the surface with a rise in value of the unsteadiness parameter A increased. However, the situation was reversed within the boundary layer, and the local entropy generation number decreased with an increase in the unsteadiness parameter A . Figure 7a demonstrates that the temperature exponent A o had a decreasing impact on N G near the exponentially stretching surface where the value of N G was augmented with an increase in A o within the boundary layer region. The impact of the group parameter B r / Ω on N G is depicted in Figure 7b. As the group parameter got higher values, the viscous dissipation effects were augmented, which resulted in increase in the value of N G . Figure 7c reveals that the effects of the Reynolds number on N G were increasing.
The impacts of several flow parameters on the Bejan number B e are presented in Figure 8 and Figure 9. Figure 8a presents the effects of M on B e . It can be seen that the entropy effects due to heat transfer were significant near the vicinity of the exponentially stretching sheet. With a rise in the value of M , the heat transfer entropy effects started to decrease a bit. However, within the boundary layer region and the far away regime, the entropy effects due to heat transfer became significantly dominant over the entropy effects due to the viscous dissipation and the magnetic field with an increase in M . Figure 8b shows that in the absence of slip effects (i.e., when γ = 0 ), the impacts of entropy because of the viscous dissipation and the magnetic field were prominent. As the value of the slip parameter γ increased, the value of the Bejan number increased. This indicates that the entropy effects due to heat transfer started to become dominant over the viscous dissipation and magnetic field entropy effects, as illustrated in the figure. The impact of A on B e is displayed in Figure 8c. It can be seen that the value of B e lied within the range 0.5 < B e 1 in the neighborhood of the exponentially stretching surface. This shows that the heat transfer entropy effects were prominent in this region, and these effects became stronger with increase in the value of A . Moreover, it was observed that the entropy effects due to the viscous dissipation and the magnetic field became significantly stronger within the boundary layer and the far away region as the value of A increased. Thus, there was a dominance of the viscous dissipation and magnetic field entropy effects over the heat transfer entropy effects in these regions. Figure 9a depicts the impact of A o on B e . This shows that rise in the value of A o resulted a slight boost in the value of B e . Additionally, it was seen that the heat transfer entropy effects were prominent at the surface of the exponentially stretching sheet, and the viscous dissipation and magnetic field entropy effects were significant in the faraway region. Figure 9b indicates that the entropy effects due to the viscous dissipation and the magnetic field became dominant over the heat transfer entropy effects with an increase in the value of the group parameter B r / Ω . Figure 9c reveals the impact of increasing Re L on B e . The value of B e decreased very slightly near the exponentially stretching sheet with an increase in the Reynolds number Re L . Within the boundary layer, a slight increasing behavior was observed as the value of Re L increased. Again, a decreasing effect was seen on B e in the far away regime with rise in the value of Re L . Moreover, it was witnessed that heat transfer entropy effects were prominent at the surface, whereas the entropy due to the viscous dissipation and the magnetic field was significant in the faraway region.

6. Conclusions

In this study, the impacts of entropy generation in an MHD unsteady slip flow over an exponentially stretching sheet were explored. The problem was solved by employing both the HAM and the shooting method, and the outcomes were briefly discussed. It was observed that the resistive force reduced the fluid velocity with the rise in the magnetic field parameter M and caused the temperature to augment. Consequently, an increase in entropy was witnessed, as it significantly depended upon the fluid friction. Thus, the presence of magnetic field plays a notable role in the production of entropy effects. Moreover, as the value of the magnetic field parameter M increased, the Lorentz force became strong, which resisted the fluid movement and caused the decline in velocity f ( ξ ) . The heat transfer rate from the exponentially stretching surface increased with a rise in the slip parameter γ , the unsteadiness parameter A , the temperature exponent A 0 , and the Prandtl number Pr . On the contrary, reverse impacts were observed in the case of M and E c . It can be deduced that the presence of a slip at the surface of exponentially stretching sheet depreciates the entropy effect, whereas the application of an externally applied magnetic field augments the entropy effects. Moreover, unsteady flow also has an increasing effect on entropy production. The heat transmission entropy effects are dominant at the exponentially stretching sheet, whereas the viscous dissipation and magnetic field entropy impacts dominate the heat transmission entropy impacts in the faraway region. The value of N G rises with an increase in the unsteadiness parameter A near the surface and declines within the boundary layer regime. In the case of the temperature exponent A o , an opposite behavior was observed as compared to the unsteadiness parameter. The heat transmission entropy effects are dominant at the exponentially stretching sheet, whereas the viscous dissipation and magnetic field entropy impacts dominate the heat transmission entropy impacts in the faraway region. Finally, we conclude that the heat transmission entropy impacts are the major cause of entropy production.

Author Contributions

S.H., A.S.B. and Y.-Z.L. wrote the main manuscript text and performed mathematical modeling. S.M.I., S.H. and A.S.B. prepared the figures and tables. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This effort received funding from the National Natural Science Foundation of China, Grant/ Award Number: 11271037.

Conflicts of Interest

The author(s) declare no competing interest.

References

  1. Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. 1970, 21, 645–647. [Google Scholar] [CrossRef]
  2. Gupta, P.S.; Gupta, A.S. Heat and mass transfer on a stretching sheet with suction or blowing. Can. J. Chem. Eng. 1977, 55, 744–746. [Google Scholar] [CrossRef]
  3. Grubka, J.; Bobba, K.M. Heat transfer characteristics of a continuous stretching surface with variable temperature. Trans. ASME J. Heat Trans. 1985, 107, 248–250. [Google Scholar] [CrossRef]
  4. Chen, C.-K.; Char, M.-I. Heat transfer of a continuous, stretching surface with suction or blowing. J. Math. Anal. Appl. 1988, 135, 568–580. [Google Scholar] [CrossRef] [Green Version]
  5. Elbashbeshy, A.E.M. Heat transfer over a stretching surface with variable surface heat flux. J. Phys. D Appl. Phys. 1998, 31, 1951–1954. [Google Scholar] [CrossRef]
  6. Magyari, E.; Keller, B. Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface. J. Phys. D Appl. Phys. 1999, 32, 577–585. [Google Scholar] [CrossRef]
  7. Elbashbeshy, E.M.A. Heat transfer over an exponentially stretching continuous surface with suction. Arch. Mech. 2001, 53, 643–651. [Google Scholar]
  8. Al-Odat, M.Q.; Damseh, R.A.; Al-Azab, T.A. Thermal boundary layer on an exponentially stretching continuous surface in the presence of magnetic field effect. Int. J. Appl. Mech. Eng. 2006, 11, 289–299. [Google Scholar]
  9. El-Aziz, M.A. Viscous dissipation effect on mixed convection flow of a micro polar fluid over an exponentially stretching sheet. Can. J. Phys. 2009, 87, 359–368. [Google Scholar] [CrossRef]
  10. Bidin, B.; Nazar, R. Numerical solution of the boundary layer flow over an exponentially stretching sheet with thermal radiation. Eur. J. Sci. Res. 2009, 4, 710–717. [Google Scholar]
  11. Ishak, A. MHD Boundary layer flow due to an exponentially stretching sheet with radiation effect. Sains Malays. 2011, 40, 391–395. [Google Scholar]
  12. Mukhopadhyay, S.; Gorla, R.S.R. Effects of partial slip on boundary layer flow past a permeable exponential stretching sheet in presence of thermal radiation. Heat Mass Transf. 2012, 48, 1773–1781. [Google Scholar] [CrossRef]
  13. Sahoo, B.; Poncet, S. Flow and heat transfer of a third grade fluid past an exponentially stretching sheet with partial slip boundary condition. Int. J. Heat Mass Transf. 2011, 54, 5010–5019. [Google Scholar] [CrossRef] [Green Version]
  14. Bhattacharyya, K.; Vajravelu, K. Stagnation-point flow and heat transfer over an exponentially shrinking sheet. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 2728–2734. [Google Scholar] [CrossRef]
  15. Mukhopadhyay, S.; Vajravelu, K.; Van Gorder, R.A. Casson Fluid Flow and Heat Transfer at an Exponentially Stretching Permeable Surface. J. Appl. Mech. 2013, 80, 054502. [Google Scholar] [CrossRef]
  16. Jyothi, K.P.; Sreelakshmi, K.; Nagendramma, V. Sarojamma, Study of thermophoresis on the MHD flow due to an exponentially stretching sheet in the presence of viscous dissipation. Proc. Eng. 2015, 127, 340–346. [Google Scholar]
  17. Reddy, P.B.A.; Bareddy, P. Magnetohydrodynamic flow of a Casson fluid over an exponentially inclined permeable stretching surface with thermal radiation and chemical reaction. Ain Shams Eng. J. 2016, 7, 593–602. [Google Scholar] [CrossRef] [Green Version]
  18. Patil, P.; Latha, D.; Roy, S.; Momoniat, E. Non-similar solutions of mixed convection flow from an exponentially stretching surface. Ain Shams Eng. J. 2017, 8, 697–705. [Google Scholar] [CrossRef] [Green Version]
  19. Jusoh, R.; Nazar, R.; Pop, I. Magnetohydrodynamic rotating flow and heat transfer of ferrofluid due to an exponentially permeable stretching/shrinking sheet. J. Magn. Magn. Mater. 2018, 465, 365–374. [Google Scholar] [CrossRef]
  20. Bejan, A. Entropy Generation Minimization; CRC: New York, NY, USA, 1996. [Google Scholar]
  21. Tamayol, A.; Hooman, K.; Bahrami, M. Thermal analysis of flow in a porous medium over a permeable stretching wall. Transp. Porous Med. 2010, 85, 661–676. [Google Scholar] [CrossRef]
  22. Aïboud, S.; Saouli, S. Entropy analysis for viscoelastic magnetohydrodynamic flow over a stretching surface. Int. J. Non-Linear Mech. 2010, 45, 482–489. [Google Scholar] [CrossRef]
  23. Butt, A.S.; Munawar, S.; Ali, A.; Mehmood, A. Entropy analysis of mixed convective magnetohydrodynamic flow of a viscoelastic fluid over a stretching sheet. Z. Naturforsch 2012, 67, 451–459. [Google Scholar] [CrossRef]
  24. Butt, A.S.; Ali, A.; Munawar, S. Slip effects on entropy generation in MHD flow over a stretching surface in the presence of thermal radiation. Int. J. Exergy 2013, 13, 1–20. [Google Scholar] [CrossRef]
  25. Noghrehabadi, A.; Saffarian, M.R.; Pourrajab, R.; Ghalambaz, M. Entropy analysis for nanofluid flow over a stretching sheet in the presence of heat generation/absorption and partial slip. J. Mech. Sci. Technol. 2013, 27, 927–937. [Google Scholar] [CrossRef]
  26. Butt, A.S.; Ali, A. A computational study of entropy generation in magnetohydrodynamic flow and heat transfer over an unsteady stretching permeable sheet. Eur. Phys. J. Plus 2014, 129, 13. [Google Scholar] [CrossRef]
  27. Rashidi, M.; Mohammadi, F.; Abbasbandy, S.; Alhuthali, M.S. Tongji University 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]
  28. Afridi, M.I.; Qasim, M.; Khan, I.; Shafie, S. Entropy generation in magnetohydrodynamic mixed convection flow over an inclined stretching sheet. Entropy 2017, 19, 10. [Google Scholar] [CrossRef] [Green Version]
  29. Butt, A.S.; Ali, A.; Nazimtufail, M.; Mehmood, A. Theoretical Investigation of Entropy Generation Effects in Magnetohydrodynamic Flow of Casson Nanofluid over an Unsteady Permeable Stretching Surface. J. Nanofluids 2019, 8, 103–116. [Google Scholar] [CrossRef]
  30. Riaz, A.; Bhatti, M.M.; Ellahi, R.; Zeeshan, A.; Sait, S.M. Mathematical Analysis on an Asymmetrical Wavy Motion of Blood under the Influence Entropy Generation with Convective Boundary Conditions. Symmetry 2020, 12, 102. [Google Scholar] [CrossRef] [Green Version]
  31. Zaib, A.; Khan, U.; Khan, I.; Seikh, A.H.; Sherif, E.-S.M. Khan Numerical Investigation of Aligned Magnetic Flow Comprising Nanoliquid over a Radial Stretchable Surface with Cattaneo—Christov Heat Flux with Entropy Generation. Symmetry 2019, 11, 1520. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Coordinate system of fluid flow.
Figure 1. Coordinate system of fluid flow.
Symmetry 12 00426 g001
Figure 2. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the velocity profile f ( ξ ) .
Figure 2. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the velocity profile f ( ξ ) .
Symmetry 12 00426 g002aSymmetry 12 00426 g002b
Figure 3. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the temperature profile θ ( ξ ) .
Figure 3. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the temperature profile θ ( ξ ) .
Symmetry 12 00426 g003
Figure 4. Effects of (a) the temperature exponent A o , (b) the Prandtl number Pr , (c) the Eckert number E c on the temperature profile θ ( ξ ) .
Figure 4. Effects of (a) the temperature exponent A o , (b) the Prandtl number Pr , (c) the Eckert number E c on the temperature profile θ ( ξ ) .
Symmetry 12 00426 g004
Figure 5. Variation of relative entropy generation for different entropy generation factors.
Figure 5. Variation of relative entropy generation for different entropy generation factors.
Symmetry 12 00426 g005
Figure 6. Effects of (a) the magnetic field parameter M , (b) the slip paramete γ , and (c) the unsteadiness parameter A on local entropy generation number N G .
Figure 6. Effects of (a) the magnetic field parameter M , (b) the slip paramete γ , and (c) the unsteadiness parameter A on local entropy generation number N G .
Symmetry 12 00426 g006aSymmetry 12 00426 g006b
Figure 7. Effects of (a) the temperature exponent A o , (b) the group parameter B r / Ω , and (c) the Reynolds number Re L on the local entropy generation number N G .
Figure 7. Effects of (a) the temperature exponent A o , (b) the group parameter B r / Ω , and (c) the Reynolds number Re L on the local entropy generation number N G .
Symmetry 12 00426 g007aSymmetry 12 00426 g007b
Figure 8. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the Bejan number B e .
Figure 8. Effects of (a) the magnetic field parameter M , (b) the slip parameter γ , and (c) the unsteadiness parameter A on the Bejan number B e .
Symmetry 12 00426 g008
Figure 9. Effects of (a) the temperature exponent A o , (b) the group parameter B r / Ω , (c) the Reynolds number Re L on the Bejan number.
Figure 9. Effects of (a) the temperature exponent A o , (b) the group parameter B r / Ω , (c) the Reynolds number Re L on the Bejan number.
Symmetry 12 00426 g009
Table 1. Optimal values of f and θ , as well as their corresponding averaged root mean squared residual errors E m , f and E m , θ at different orders of approximation when M = 1.0 ,   γ = 0.5 ,   A = 1.0 , A 0 = 0.5 , Pr = 1.0 ,   E c = 0.2 .
Table 1. Optimal values of f and θ , as well as their corresponding averaged root mean squared residual errors E m , f and E m , θ at different orders of approximation when M = 1.0 ,   γ = 0.5 ,   A = 1.0 , A 0 = 0.5 , Pr = 1.0 ,   E c = 0.2 .
O r d e r   o f   A p p r o x i m a t i o n   m f θ E m , θ E m , θ
3 0.464353 0.396157 1.05905 × 10 2 2.69482 × 10 3
6 0.405309 0.424927 1.81301 × 10 5 1.25275 × 10 4
9 0.382564 0.413877 1.49644 × 10 5 2.93299 × 10 5
12 0.369398 0.368290 8.38391 × 10 7 9.89133 × 10 6
15 0.314784 0.326304 2.48338 × 10 6 3.91650 × 10 6
Table 2. Convergence table for the [ m / m ] homotopy Padé approximation of f ( 0 )   and   θ ( 0 ) when f = 0.314784 ,   θ = 0.398461 , M = 1.0 ,   γ = 0.5 ,   A = 1.0 , A 0 = 0.5 , Pr = 1.0 ,   E c = 0.2 .
Table 2. Convergence table for the [ m / m ] homotopy Padé approximation of f ( 0 )   and   θ ( 0 ) when f = 0.314784 ,   θ = 0.398461 , M = 1.0 ,   γ = 0.5 ,   A = 1.0 , A 0 = 0.5 , Pr = 1.0 ,   E c = 0.2 .
[ m / m ] f ( 0 ) θ ( 0 )
[ 5 / 5 ] 0.95890 1.93308
[ 10 / 10 ] 0.95923 1.92904
[ 15 / 15 ] 0.95921 1.92905
[ 20 / 20 ] 0.95921 1.92905
[ 25 / 25 ] 0.95921 1.92905
[ 30 / 30 ] 0.95921 1.92905
[ 35 / 35 ] 0.95921 1.92905
Table 3. Comparison of present values of θ ( 0 ) with those reported by Magyari and Keller [6], El-Aziz [9], Nazar [10] and Ishak [11] for variation in the Prandtl number Pr and the magnetic field parameter M when γ = 0.0 , A = 0.0 , A o = 1.0 , E c = 0.0 .
Table 3. Comparison of present values of θ ( 0 ) with those reported by Magyari and Keller [6], El-Aziz [9], Nazar [10] and Ishak [11] for variation in the Prandtl number Pr and the magnetic field parameter M when γ = 0.0 , A = 0.0 , A o = 1.0 , E c = 0.0 .
M Pr θ ( 0 ) Magyari and Keller [6] θ ( 0 ) E1-Aziz [9] θ ( 0 ) Bidin and Nazar [10] θ ( 0 ) Anuar Ishak [11]Present
θ ( 0 ) HAM θ ( 0 ) Numerical
0.0 1.0 0.954782 0.954785 0.9548 0.9548 0.95478 0.95478
2.0 1.4714 1.4715 1.47146 1.47146
3.0 1.869075 1.869072 1.8691 1.86907 1.86907
5.0 2.500135 2.500132 2.5001 2.50012 2.50012
10.0 3.660379 3.660372 3.6604 3.66027 3.66027
1.0 1.0 0.8611 0.86109 0.86109
Table 4. Numerical values of f ( 0 ) that were obtained by the homotopy analysis method (HAM) and the shooting method for variation in the values of the slip parameter γ , the magnetic field parameter M , and the unsteadiness parameter A .
Table 4. Numerical values of f ( 0 ) that were obtained by the homotopy analysis method (HAM) and the shooting method for variation in the values of the slip parameter γ , the magnetic field parameter M , and the unsteadiness parameter A .
γ M A f ( 0 ) H A M f ( 0 ) N u m e r i c a l
0.5 1.0 0.2 0.85691 0.85691
0.5 0.89946 0.89946
1.0 0.95920 0.95920
1.5 1.00841 1.00841
0.5 0.0 1.0 0.88048 0.88048
1.0 0.95920 0.95920
2.0 1.02039 1.02039
3.0 1.06998 1.06998
0.0 1.0 1.0 2.04499 2.04499
0.5 0.95920 0.95920
1.0 0.63912 0.63912
2.0 0.38676 0.38676
Table 5. Numerical values of θ ( 0 ) that were obtained by the HAM and the shooting method for variation in the values of the slip parameter γ , the magnetic field parameter M , the unsteadiness parameter A , the temperature exponent A o , the Prandtl number Pr and the Eckert number E c .
Table 5. Numerical values of θ ( 0 ) that were obtained by the HAM and the shooting method for variation in the values of the slip parameter γ , the magnetic field parameter M , the unsteadiness parameter A , the temperature exponent A o , the Prandtl number Pr and the Eckert number E c .
γ M A A 0 Pr E c θ ( 0 )   H A M θ ( 0 ) N u m e r i c a l
0.5 1.0 1.0 0.5 1.0 0.2 1.92905 1.92905
0.5 1.86944 1.86944
0.7 1.82970 1.82970
0.9 1.78996 1.78996
0.5 1.0 1.0 0.5 0.7 0.2 1.61368 1.61368
1.0 1.92905 1.92905
2.0 2.72687 2.72687
3.0 3.33743 3.33743
0.5 1.0 1.0 0.5 1.0 0.2 1.92905 1.92905
1.0 1.97510 1.97510
1.5 2.02057 2.02057
2.0 2.06548 2.06548
0.5 1.0 0.2 0.5 1.0 0.2 0.96247 0.96247
0.5 1.40692 1.40692
1.0 1.92905 1.92905
1.5 2.33731 2.33731
0.5 0.0 1.0 0.5 1.0 0.2 1.95329 1.95329
1.0 1.92905 1.92905
2.0 1.91235 1.91235
3.0 1.90015 1.90015
0.0 1.0 1.0 0.5 1.0 0.2 1.88525 1.88525
0.5 1.92905 1.92905
1.0 1.92378 1.92378
2.0 1.91240 1.91240

Share and Cite

MDPI and ACS Style

Haider, S.; Saeed Butt, A.; Li, Y.-Z.; Imran, S.M.; Ahmad, B.; Tayyaba, A. Study of Entropy Generation with Multi-Slip Effects in MHD Unsteady Flow of Viscous Fluid Past an Exponentially Stretching Surface. Symmetry 2020, 12, 426. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030426

AMA Style

Haider S, Saeed Butt A, Li Y-Z, Imran SM, Ahmad B, Tayyaba A. Study of Entropy Generation with Multi-Slip Effects in MHD Unsteady Flow of Viscous Fluid Past an Exponentially Stretching Surface. Symmetry. 2020; 12(3):426. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030426

Chicago/Turabian Style

Haider, Sajjad, Adnan Saeed Butt, Yun-Zhang Li, Syed Muhammad Imran, Babar Ahmad, and Asia Tayyaba. 2020. "Study of Entropy Generation with Multi-Slip Effects in MHD Unsteady Flow of Viscous Fluid Past an Exponentially Stretching Surface" Symmetry 12, no. 3: 426. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030426

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