Next Article in Journal
Bearing Remaining Useful Life Prediction Based on Naive Bayes and Weibull Distributions
Next Article in Special Issue
A Simple Thermodynamic Model of the Internal Convective Zone of the Earth
Previous Article in Journal
Morphological Computation: Nothing but Physical Computation
Previous Article in Special Issue
Entropy Analysis of 3D Non-Newtonian MHD Nanofluid Flow with Nonlinear Thermal Radiation Past over Exponential Stretched Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Irreversibility Analysis of Dissipative Fluid Flow Over A Curved Surface Stimulated by Variable Thermal Conductivity and Uniform Magnetic Field: Utilization of Generalized Differential Quadrature Method

by
Muhammad Idrees Afridi
1,
Abderrahim Wakif
2,
Muhammad Qasim
1 and
Abid Hussanan
3,4,*
1
Department of Mathematics, COMSATS University Islamabad (CUI), Park Road, Tarlai Kalan, Islamabad 455000, Pakistan
2
Laboratory of Mechanics, Faculty of Sciences Aïn Chock, Hassan II University, B.P. 5366 Mâarif, Casablanca 20000, Morocco
3
Division of Computational Mathematics and Engineering, Institute for Computational Science, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
4
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(12), 943; https://doi.org/10.3390/e20120943
Submission received: 10 October 2018 / Revised: 18 November 2018 / Accepted: 3 December 2018 / Published: 7 December 2018
(This article belongs to the Special Issue Entropy Generation and Heat Transfer)

Abstract

:
The effects of variable thermal conductivity on heat transfer and entropy generation in a flow over a curved surface are investigated in the present study. In addition, the effects of energy dissipation and Ohmic heating are also incorporated in the modelling of the energy equation. Appropriate transformations are used to develop the self-similar equations from the governing equations of momentum and energy. The resulting self-similar equations are then solved by the Generalized Differential Quadrature Method (GDQM). For the validation and precision of the developed numerical solution, the resulting equations are also solved numerically using the Runge-Kutta-Fehlberg method (RKFM). An excellent agreement is found between the numerical results of the two methods. To examine the impacts of emerging physical parameters on velocity, temperature distribution and entropy generation, the numerical results are plotted against the various values of physical flow parameters and discussed physically in detail.

1. Introduction

The first law of thermodynamics gives a quantitative estimate of heat and work interactions between some system and surroundings if the system undergoes a thermodynamic process or a cycle. However, it does not say whether the process or the cycle in a particular direction would occur or not. Further, the first law does not indicate whether conversion of energy from one form to another are performed perfectly or whether some forms are completely converted to others. The first law establishes the equivalence of heat and work and shows there is a fixed rate of exchange between heat and work. It does not talk about the conditions under which the transformations of energy are possible. Second law of thermodynamics puts a condition on the conversion of energy from one form to another. The second law states that it is not possible to convert heat energy completely into work. That part of heat energy which cannot be converted into work is known as unavailable energy and must be rejected as a low-grade energy. This means that, the availability of energy in a thermal system decreases. The phenomenon of increasing unavailable energy in a thermal process is called entropy generation. The thermal efficiency reduces with increasing entropy generation and therefore needs to examine the factors that reduces the entropy generation. Bejan [1] introduced the innovative idea of reducing the unavailable energy (entropy generation) in a convective heat transfer problem. After the pioneering work of Bejan [1], analysis of entropy generation in fluid flow problems are reported by many researchers. Recently, Afridi et al. [2] studied the entropy generation of carbon nanotubes CNTs nanofluids in a flow over a thin needle by incorporating the effects of nonlinear thermal radiation and viscous dissipation. Makinde [3] reported the effects of variable viscosity on inherent irreversibility in a flow over a flat plate with Newtonian heating and uniform magnetic field. The combined effects of linear thermal radiation and viscous dissipation on entropy generation in a Blasius flow are reported by Butt et al. [4]. Recently, Afridi and Qasim [5] examined the influences of frictional heating on entropy production rate in a three-dimensional flow. Entropy generation in a non-linear convection flow over a vertical plate with convective boundary condition and porous medium is studied by Vasu et al. [6]. Afridi et al. [7] studied the exact effects of viscous dissipation, Joule heating and heat transfer on entropy generation in a flow of Newtonian fluid over an elastic stretching boundary. The heat transfer and irreversibility analysis of nanofluid flow containing copper nanoparticles with water as a base fluid is reported by Butt et al. [8]. Makinde and Eegunjobi [9] studied the heat transfer effects in a couple stress fluid flow inside a vertical channel filled with porous medium with entropy analysis. Rashidi et al. [10] investigated the entropy generation in a nanofluid flow over a rotating disk under the influence of uniform magnetic field.
Boundary layer flows have significant number of applications in industrial and engineering processes such as extraction of polymer sheet, glass-fiber production, spinning of fibers, wire drawing, extruded plastic sheets, paper production, hot rolling and materials handling conveyors. Boundary layer flow is the corner stone of modern fluid dynamics due to vital application in manufacturing processes [11]. After the seminal work of Sakiadis [12], many researchers examined the boundary layer flow either by using analytic techniques or by utilizing various numerical methods. Crane [13] reported the exact solution of boundary layer flow driven by a stretching surface. Gupta and Gupta [14] reported the effects of transpiration on heat and mass transfer in a flow over a stretching boundary. Flow driven by a stretching cylinder is reported by Wang [15]. Wang [16] also reported the rotating fluid flow over a stretching boundary. The thin film flow over a stretching surface with variable fluid properties is studied by Dandapat et al. [17]. Vajravelu and Rollins [18] reported the flow of second grade fluid under the effect of magnetic field. Pal and Mondal [19] examined the effects of prescribed surface temperature (PST) and prescribed heat flux (PHF) on the mixed convection fluid with variable transport properties. The effects of energy dissipation with prescribed surface temperature on a power-law fluid flow driven by a permeable boundary are examined by Yazdi et al. [20]. Hsiao [21] analyzed the influences of viscous dissipation on viscoelastic fluid. The flow over a slender stretching sheet with Hall effects and variable fluid properties is studied by Vajravelu et al. [22]. Mixed convection flow of Casson fluid over a vertical flat surface is investigated by Vajravelu et al. [23]. Some of other recent studies on a boundary layer flow over a flat surface and curved surface are referenced in [24,25,26,27,28,29]. The effects of magnetic field on fluid flow are investigated in [30,31,32].
The objective of the present article is to study the heat transfer effects in a dissipative fluid flow over an elastic curved surface with variable thermal conductivity. One of the most important effect in the boundary layer flow knowns as viscous dissipation is also considered by adding the dissipation function the energy equation. Besides this, the entropy generation analysis is also performed. Numerical solutions of the transformed set of highly nonlinear differential equations are obtained and utilized to compute the entropy generation number. Generalized Differential Quadrature Method (GDQM) is used to get the numerical solutions. The numerical results are presented graphically and discussed.

2. Problem Formulation

Let us consider a two-dimensional incompressible steady flow of a Newtonian electrically conducting fluid over a linearly curved stretching surface as schematically described in Figure 1. Moreover, this curved boundary (i.e., r = R) is taken to be coiled in a circle of radius R. In this investigation, we consider that the induced boundary layer flow is significantly affected by a radial magnetic field with constant strength B0, in which the fluid is flowing under the combined effects of viscous dissipation and Joule heating with the presence of a temperature dependent behavior for the thermal conductivity in the form k* = k ω(T). In order to define the flow geometry, we choose to employ the curvilinear coordinates (r,s) as the best way for modelling the present problem, where r is the radial coordinate measured from the center O of the curved surface, whereas s is the coordinate of the arc length along the flow direction. Furthermore, the stretching velocity and temperature of the curved sheet are taken to be in the form uw(s) = uos and Tw(s) = Tb + Tos2, respectively, where Tb represents the temperature of the bulk fluid, so that uo and T0 are two dimensional constants.
Based on the Prandtl boundary layer approximations and the afore-mentioned assumptions, the basic equations of continuity, momentum and thermal energy are written as follows:
  r ( r ¯ u r ) + R u s s = 0   ,
1 r ¯ u s 2 = 1 ρ p r   ,
ρ ( u r u s r + R r ¯ u s u s s + 1 r ¯ u s u r ) = R r ¯ p s + μ ( 2 u s r 2 + 1 r ¯ u s r 1 r ¯ 2 u s ) σ B o 2 u s   ,
ρ c p ( u r T r + R r ¯ u s T s ) =   1   r * r ( r ¯ k * T r ) + μ ( u s r 1 r ¯ u s ) 2 + σ B o 2 u s 2 .
Here, the symbols ρ, μ, σ and k* represents the thermo-physical properties of the electrically conducting fluid. These quantities denote the density, dynamic viscosity, electrical conductivity and thermal conductivity of the studied fluid, where r ¯ is the reduced radial variable, where r ¯ = r + R .
Additionally, it is worth noting here that the characteristic function ω(T) mentioned above is taken in the form:
ω ( T ) = 1 + τ ( T T b T w T b ) ,
where τ is an adjusted parameter showing the importance of the temperature dependence in the thermal conductivity k*.
For the present two-dimensional steady flow model, the associated physical boundary conditions are written in curvilinear coordinates as follows:
u s = u w ,   T = T w   at   r = 0 ,
u s 0 ,   u s r 0 ,   T T b   as   r .
Here, u s and u r are the velocity components in s and r directions, respectively, Bo denotes the magnetic field strength, p and T show the pressure and temperature of the fluid, respectively, Tw represents the temperature of the curved sheet, whereas Tb indicates the temperature of the fluid in the stress free region.
Now, by introducing the following similarity transformations:
χ = ( u o ρ μ   ) 0.5 r ,   g ( χ ) = u s ( r , s ) u w ,   g ( χ ) = ( ρ u o μ ) 0.5 ( r ¯ R ) u r ( r , s ) ,   θ = T T b T w T b ,   P ( χ ) = 1 ρ u o 2 s 2 p ,
Equations (1)–(4) reduce to
P χ = g 2   h   ,
2 κ h P = g + g h ( 1 h 2 + M ) g +   κ h g g + κ h 2 g g   κ h g 2 ,
1 Pr ( 1 + τ θ ) ( θ + θ h ) + τ Pr θ 2 +   κ h ( g θ 2 g θ ) + E c ( g g h ) 2 + E c M g 2 = 0   .
Here, κ denotes the curvature parameter, Ec and Pr represent the Eckert and Prandtl numbers, respectively, whereas M represents the magnetic parameter, where:
κ = R ( u o ρ μ   ) 0.5 ,   E c = u w 2 c p ( T w T b ) ,   Pr = c p μ k ,   M = B o 2 σ u o ρ .
In this paper, we use the prime notation (e.g., g’ or θ’) in Equations (8)–(11) as subscript to denote the derivative with respect to χ . Furthermore, the function h shown in Equations (9)–(11) is a linear function of the similarity variable χ , which is given by:
h ( χ ) = χ + κ   .
Also, after combining Equations (9) and (10), we get:
g + 2 h g g ¯ 1 g + g ¯ 2 g +   κ h ( g g + g h g 1 h g 2 g g g h 2 g ) = 0 ,
where:
g ¯ 1 = 1   h 2 + M ,
g ¯ 2 = 1   h 3 M h .
Upon making use of the transformations (8), the boundary conditions (6) and (7) become:
g ( 0 ) = 0 ,   g ( 0 ) = 1 ,   θ ( 0 ) = 1 ,
g ( χ ) 0 ,   g ( χ ) 0 ,   θ ( χ ) 0 .
Furthermore, the important physical quantities of practical interest arising from this investigation are the local skin friction coefficient Cfs and the local Nusselt number Nus, which are expressed by:
Re s 0.5 C f s = g ( 0 ) 1 κ g ( 0 ) ,
Re s 0.5 N u s = θ ( 0 ) τ θ 2 ( 0 ) .
Here, Res represents the local Reynolds number, where Re s = u o s 2 / ν .
For more helpful simplifications, we can put:
{ χ = χ η , h ( χ ) = h ( χ η ) = H ( η ) , g ¯ 1 ( χ ) = g ¯ 1 ( χ η ) = G ¯ 1 ( η ) , g ¯ 2 ( χ ) = g ¯ 1 ( χ η ) = G ¯ 2 ( η ) , g ( χ ) = g ( χ η ) = G ( η ) , θ ( χ ) = θ ( χ η ) = Θ ( η ) .
Keeping in mind the above considerations, the derivatives of g(χ) and θ(χ) can be expressed as function of the derivatives of G(η) and Θ(η), respectively, as follows:
{ g ( n ) ( χ ) = G ( n ) ( η ) χ n , θ ( n ) ( χ ) = Θ ( n ) ( η ) χ n ,
where n denotes the integer-order derivative with respect to the spatial variables χ or η.
Accordingly, Equations (11) and (14) with the boundary conditions (17) and (18) can be written in the following general form:
L G ( G ) + N G ( G , Θ )   = 0 ,  
L Θ ( Θ ) + N Θ ( G , Θ )   = 0 ,  
G ( η ) = 0 ,   G ( η ) = χ ,   Θ ( η ) = 1 ,   at   η = 0 ,
G ( η ) 0 ,   G ( η ) 0 ,   Θ ( η ) 0 ,   as   η 1 ,
in which:
L G ( G ) =   G + 2 χ H G χ 2 G ¯ 1 G + χ 3 G ¯ 2 G ,
L Θ ( Θ ) = χ 2 Pr Θ + χ 3 Pr H Θ ,
N G ( G , Θ ) = χ κ H ( G G + χ G H G χ H G 2 G G χ 2 G H 2 G ) ,
N Θ ( G , Θ ) = { χ 2 τ Pr Θ 2 + χ 3 τ Θ Pr H Θ + χ 2 τ Pr Θ Θ + χ 3 κ G H Θ 2 χ 3 κ Θ H G + E c G 2 2 χ E c H G G + χ 2 E c G ¯ 1 G 2 } .
Here, χ represents the optimum value of the boundary layer thickness, which ensures our numerical findings are approached asymptotically to their exact values.
After substituting Equation (21) into Equations (19) and (20), the physical quantities (Res)0.5 Cfs and (Res)−0.5Nus become:
Re s 0.5 C f s = 1   χ 2 G ( 0 ) 1 χ κ G ( 0 ) ,
Re s 0.5 N u s = 1 χ Θ ( 0 ) τ χ 2 Θ 2 ( 0 ) .

3. Analysis of Entropy Production

As is well known, the local volumetric rate of entropy generation ξ g of a fluidic system in the presence of viscous dissipation and Ohmic heating is given by:
ξ g = ξ t + ξ f + ξ m ,
where ξ t represents the entropy generation due to heat transfer across a finite temperature difference, ξ f shows the local entropy generation due to viscous dissipation and ξ m characterizes the local entropy generation due to the presence of Lorentz force, where:
ξ t = k ω ( T ) T 2 ( T r ) 2 ,
ξ f = μ T ( u s r u s r * ) 2 ,
ξ m = σ B o 2 T u s 2 .
By considering the following characteristic entropy generation:
ξ c = k u o ν ,
the entropy generation in dimensionless form can be written as follows:
N s = ξ g ξ c = ( 1 + τ θ ) ( θ + λ ) 2 θ 2 T h e r m a l   c o n t r i b u t i o n + E c Pr ( θ + λ ) ( g g h ( ξ ) ) 2 F r i c t i o n a l   c o n t r i b u t i o n + M E c Pr ( θ + λ ) g 2 M a g n e t i c   c o n t r i b u t i o n ,
where λ = T b / ( T w T b ) denotes the temperature difference parameter.
By virtue of Equation (21), the entropy generation Ns reduces to:
N s = ( 1 + τ Θ ) χ 2 ( Θ + λ ) 2 Θ 2 + E c Pr χ 4 ( Θ + λ ) ( G χ G H   ) 2 + M E c Pr χ 2 ( Θ + λ ) G 2 .

4. Solution Methodology

The boundary layer flow model induced over the curved stretching surface r = R can be regarded as a complicated two-point boundary value problem. As mentioned previously in Equations (23) and (24), the present physical problem is governed by a set of ordinary differential equations ODEs, which are highly nonlinear. From the mathematical point of view, the flow and heat transfer characteristics of the studied fluid are extremely difficult to be found analytically as closed form solutions. Therefore, for solving this kind of physical problems, it is more recommended here to adopt a powerful numerical method in terms of accuracy and efficiency to predict approximate numerical solutions for Equations (23) and (24) along with the boundary conditions (25) and (26). Hence, in order to achieve this objective and ensure enough accuracy, the resulting ODEs are handled numerically by discretizing the present boundary layer equations using generalized differential quadrature method (GDQM) with the following non-uniform grid points:
η i = 1 2 1 2 c o s ( π i π N 1 )   .
Here, 1 i N and η1η1ηN, where η1 = 0 and ηN = 1.
Accordingly, the discretized form of the derivatives of the functions G(η) and Θ(η) with respect to the variable η at a collocation point η i are defined as follows:
{ G ( n ) ( η i ) = j = 1 N d i j ( n ) G ( η j ) = j = 1 N d i j ( n ) G j   , Θ ( n ) ( η i ) =   j = 1 N d i j ( n ) Θ ( η j ) = j = 1 N d i j ( n ) Θ j ,
Here, d i j ( n ) are the weighting coefficients for the nth-order derivative and N is the total number of collocation points, where i and j are integers varying from 1 to N.
In order to provide enough information about the proposed numerical method, the interested readers can refer to [33,34] and the reference therein. According to the pioneer work of Shu [33], the weighting coefficients d i j ( 1 ) for the first-order derivative discretization can be expressed as follows:
{ d i j ( 1 ) = k = 1 , k i N ( η i η k ) ( η i η j ) k = 1 , k j N ( η j η k ) , for   i j   ,   d i j ( 1 ) = j = 1 , j i N d i j ( 1 ) , for   i = j   ,  
where 1 ≤ i, j ≤ N.
Similarly, the weighting coefficients d i j ( n ) for the higher-order derivatives can be found numerically using the following recurrence relations:
{ d i j ( n ) = n ( d i i ( n 1 ) d i j ( 1 ) d i j ( n 1 ) η i η j ) , for   i j   , d i j ( n ) = j = 1 , j i N d i j ( n ) , for   i = j   ,  
where 1 ≤ I, j ≤ N and n ≥ 2.
Therefore, after discretization of Equations (23)–(26), the functions H ( η ) , G ¯ 1 ( η ) , G ¯ 2 ( η ) , G ( η ) and Θ ( η ) are approximated in each collocation point η i by H i , G ¯ 1 i , G ¯ 2 i , G i and Θ i , respectively. Consequently, the discretized form of our problem is given as follows:
( S ) { G 1 = 0 ,   j = 1 N d 1 j ( 1 ) G j χ = 0 , L G i ( G i ) + N G i ( G i , Θ i )   = 0   ,   for   3 i N 2 , j = 1 N d N j ( 1 ) G j = 0 , j = 1 N d N j ( 2 ) G j = 0 , Θ 1 1 = 0 . L Θ i ( Θ i ) + N Θ i ( G i , Θ i )   = 0   ,   for   2 i N 1 Θ N = 0 ,  
in which:
L G i ( G i ) = { j = 1 N d i j ( 4 ) G j + 2 χ H i ( j = 1 N d i j ( 3 ) G j ) χ 2 G ¯ 1 i ( j = 1 N d i j ( 2 ) G j ) + χ 3 G ¯ 2 i ( j = 1 N d i j ( 1 ) G j   ) }   ,
L Θ i ( Θ i ) = χ 2 Pr ( j = 1 N d i j ( 2 ) Θ j ) + χ 3 Pr H i ( j = 1 N d i j ( 1 ) Θ j ) ,
N G i ( G i , Θ i ) = χ κ H i { G i ( j = 1 N d i j ( 3 ) G j ) + χ G i H i ( j = 1 N d i j ( 2 ) G j ) χ H i ( j = 1 N d i j ( 1 ) G j ) ( j = 1 N d i j ( 1 ) G j ) ( j = 1 N d i j ( 1 ) G j ) ( j = 1 N d i j ( 2 ) G j ) χ 2 G i H i 2 ( j = 1 N d i j ( 1 ) G j )   }
N Θ i ( G i , Θ i ) = { χ 2 τ Pr ( j = 1 N d i j ( 1 ) Θ j ) ( j = 1 N d i j ( 1 ) Θ j ) + χ 3 τ Θ i Pr H i ( j = 1 N d i j ( 1 ) Θ j ) +   χ 2 τ Pr Θ i ( j = 1 N d i j ( 2 ) Θ j ) + χ 3 κ G i H i ( j = 1 N d i j ( 1 ) Θ j ) 2 χ 3 κ Θ i H i ( j = 1 N d i j ( 1 ) G j ) + E c ( j = 1 N d i j ( 2 ) G j ) ( j = 1 N d i j ( 2 ) G j ) 2 χ E c H i ( j = 1 N d i j ( 1 ) G j ) ( j = 1 N d i j ( 2 ) G j ) + χ 2 E c G ¯ 1 i ( j = 1 N d i j ( 1 ) G j ) ( j = 1 N d i j ( 1 ) G j ) }
Here:
H i = ( κ + χ 2 )   χ 2 c o s ( π i π N 1 ) .
G ¯ 1 i = 1   H i 2 + M ,
G ¯ 2 i = 1   H i 3 M H i .
As shown in Equation (44), the algebraic nonlinear system (S) is constituted by 2N equations. This nonlinear system can be solved accurately by means of Newton-Raphson iterative scheme (NRIS). Hence, thanks to this technique, the dimensionless quantities (Res)0.5 Cfs and (Res)−0.5Nus shown in Equations (31) and (32) can be computed numerically using the following expressions:
Re s 0.5 C f s = 1 χ 2 ( j = 1 N d 1 j ( 2 ) G j ) 1 χ κ ( j = 1 N d 1 j ( 1 ) G j ) ,
Re s 0.5 N u s = 1 χ ( j = 1 N d 1 j ( 1 ) Θ j ) τ χ 2 ( j = 1 N d 1 j ( 1 ) Θ j ) ( j = 1 N d 1 j ( 1 ) Θ j ) .
From the methodological point of view, we take χ = 10 and N = 70 as the best selected values during all subsequent analyses, in order to find out significant numerical results with an absolute accuracy of the order of 10 8 . Moreover, the average CPU time taken to compute the skin friction coefficient Res0.5Cfs and the rate of heat transfer Res−0.5Nus by GDQM is generally no more than 10 s.

5. Results and Discussion

The transformed set of differential equations that govern the flow are highly nonlinear and thereby solved numerical by applying generalized differential quadrature method (GDQM) to quantify the influences of different physical flow parameter. The impacts of dimensionless flow parameters such as magnetic parameter M, curvature parameter κ , Eckert number Ec, Prandtl number Pr, variable thermal conductivity parameter τ and temperature difference parameter λ on velocity g’(χ), temperature θ(χ) and entropy generation Ns are depicted in different graphs. In order to verify the accuracy of our numerical scheme, the local skin friction coefficient and Nusselt number are also computed using Runge-Kutta Fehlberg method (RKFM) as shown in Table 1. The comparison shows an excellent agreement and hence validates our numerical simulation. Table 1 also illustrates the influences of physical flow parameters on local skin friction coefficients and Nusselt number. It is inferred from Table 1 that skin friction coefficient increases with rising values of magnetic parameter and decreases with curvature parameter. The Eckert number, Prandtl number and variable thermal conductivity have no influence on skin friction coefficient. We also observed that the local Nusselt number decreases with rising values of magnetic parameter, Eckert number and variable thermal conductivity parameter. Further, it is noted that, local Nusselt number enhances with increasing values of curvature parameter and Prandtl number.
Figure 2a,b represent the effects of magnetic and curvature parameter on velocity profile g’(χ) respectively. It is observed that motion of fluid decelerates with increasing strength of applied magnetic field. This is because a resistive force knows as Lorentz force enhances with increasing strength of applied magnetic field. It is also observed that velocity decreases with rising values of curvature parameter κ . Further, it is found that for fixed value of χ the thickness of boundary layer is thick for flow over a curved boundary as compared to flat surface ( κ ). The effects of magnetic parameter(M), curvature parameter (κ), Eckert number (Ec), Prandtl number (Pr) and variable thermal conductivity parameter (τ) on temperature profile θ(χ) are shown in Figure 3a–e respectively. It is inferred from Figure 3a,b that temperature rises with increasing values of magnetic and curvature parameter respectively. This is because, the phenomenon of Ohmic heating increases with rising values of M and thus leads to rise the fluid temperature. We also observed that for fixed value of similarity variable χ the thickness of thermal boundary layer is thin for the fluid flow past over a flat stretching surface ( κ ) as compared to flow over a curved stretching surface. The Eckert number Ec is measure of the frictional forces between the fluid layers, therefore, with increasing Eckert number the frictional heating enhances and leads to rise the fluid temperature as presented in Figure 3c. The decreasing behavior of temperature with increasing values of Prandtl number is due to the fact that thermal diffusivity decreases with increasing Prandtl number and consequently leads to drop the fluid temperature as shown in Figure 3d. The thermal conductivity of fluid increase with rising values of variable thermal conductivity parameter (τ), therefore, the temperature of fluid rises with increasing τ as shown in Figure 3e. Figure 4a shows that entropy generation enhances with enhancing the strength of applied magnetic field. This is due to the dissipative nature of the Lorentz force. In addition, it is noticed that the rate of entropy generation is maximum at the curved boundary. Figure 4b illustrates that entropy generation reduce with increasing values of curvature parameter κ . Further, no significant effects are observed at the surface of curved boundary. In addition, the rate of entropy generation is less in the flow over a flat boundary ( κ ) as compared to the curved one. Figure 4c,d demonstrate the variations of entropy generation NS with Eckert number Ec and Prandtl number Pr, respectively. We found that entropy enhances with rising values of Ec and Pr. Significant effects are observed at the curved boundary and this is due to the presence of high thermal gradients at the surface of curved boundary. Figure 4e displays that, as the variable thermal conductivity parameter τ rises, the entropy generation NS enhance slightly at the surface of curved boundary and its vicinity. Furthermore, it is found that entropy generation NS decreases after certain vertical distance from the surface of stretching curved surface. Figure 4f demonstrates that, as temperature difference parameter λ increases, entropy generation NS reduces, therefore, to minimize the entropy generation inside the boundary layer, it is suggested to reduce the operating temperature Δ T (increase λ).

6. Closing Remarks

In the present investigation, we utilized the generalized differential quadrature method (GDQM) to get the numerical solutions of the reduced set of governing nonlinear differential equations. The impacts of different physical flow parameters are investigated by plotting various graphs. Following are the key outcomes of the present study.
  • The local skin friction coefficient enhances with magnetic parameter and reduces with increasing curvature parameter.
  • With an increase in magnetic parameter, Eckert number and variable thermal conductivity parameter, the local Nusselt number reduces but it enhances with rising values of curvature parameter and Prandtl number.
  • The fluid motion decelerates with increasing M and curvature parameter κ .
  • With rising values of magnetic parameter, Eckert number and variable thermal conductivity parameter, the temperature of fluid rises whereas decrement in temperature is observed with increasing values of Prandtl number and curvature parameter.
  • Less entropy is generated in the flow past over a flat stretching boundary as compared to the flow over a curved surface.
  • By increasing the curvature and temperature difference parameter, the entropy generation N s reduces.
  • With enhancing the values of magnetic parameter, Eckert number, Prandtl number and variable thermal conductivity parameter, Ns increases.

Author Contributions

M.I.A. and M.Q. formulated the problem. M.I.A. and A.W. solved the problem. M.Q., A.H. and A.W. computed and analyzed the results. All the authors equally contributed in writing and proof reading of the paper.

Funding

This research received no external funding. The APC was given by Ton Duc Thang University, Ho Chi Minh City, Vietnam. However, no grant number is available from source.

Acknowledgments

The corresponding author would like to acknowledge Ton Duc Thang University, Ho Chi Minh City, Vietnam for the financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bejan, A. The method of entropy generation minimization. In Energy and the Environment; Springer: Berlin/Heidelberg, Germany, 1999; pp. 11–22. [Google Scholar]
  2. Afridi, M.I.; Tlili, I.; Qasim, M.; Khan, I. Nonlinear Rosseland thermal radiation and energy dissipation effects on entropy generation in CNTs suspended nanofluids flow over a thin needle. Bound. Value Probl. 2018, 2018, 148. [Google Scholar] [CrossRef]
  3. Makinde, O.D. Second law analysis for variable viscosity hydromagnetic boundary layer flow with thermal radiation and Newtonian heating. Entropy 2011, 13, 1446–1464. [Google Scholar] [CrossRef]
  4. Butt, A.S.; Munawar, S.; Ali, A.; Mehmood, A. Entropy generation in the Blasius flow under thermal radiation. Phys. Scr. 2012, 85, 1–6. [Google Scholar] [CrossRef]
  5. Afridi, M.I.; Qasim, M. Entropy generation in three dimensional flow of dissipative fluid. Int. J. Appl. Comput. Math. 2018, 4, 16. [Google Scholar] [CrossRef]
  6. Vasu, B.; Ram, R.C.; Murthy, P.; Gorla, R.S.R. Entropy generation analysis in nonlinear convection flow of thermally stratified fluid in saturated porous medium with convective boundary condition. J. Heat Transf. 2017, 139, 1–10. [Google Scholar] [CrossRef]
  7. Afridi, M.I.; Qasim, M.; Shafie, S. Entropy generation in hydromagnetic boundary flow under the effects of frictional and Joule heating: Exact solutions. Eur. Phys. J. Plus 2017, 132, 1–11. [Google Scholar] [CrossRef]
  8. Butt, A.S.; Ali, A.; Mehmood, A. Entropy generation effects in Cu water nano fluid flow and heat transfer over a radially stretching surface. J. Nanofluids 2016, 5, 471–478. [Google Scholar] [CrossRef]
  9. Makinde, O.D.; Eegunjobi, A.S. Entropy generation in a couple stress fluid flow through a vertical channel filled with saturated porous media. Entropy 2013, 15, 4589–4606. [Google Scholar] [CrossRef]
  10. Rashidi, M.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]
  11. Schlichting, H.; Gersten, K. Boundary-Layer Theory; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  12. Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7, 26–28. [Google Scholar] [CrossRef]
  13. Crane, L.J. Flow past a stretching plate. J. Appl. Math. Phys. 1970, 21, 645–647. [Google Scholar] [CrossRef]
  14. 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]
  15. Wang, C.Y. Fluid flow due to a stretching cylinder. Phys. Fluids 1988, 31, 466–468. [Google Scholar] [CrossRef]
  16. Wang, C.Y. Stretching a surface in a rotating fluid. J. Appl. Math. Phys. 1988, 39, 177–185. [Google Scholar] [CrossRef]
  17. Dandapat, B.S.; Santra, B.; Vajravelu, K. The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet. Int. J. Heat Mass Transf. 2007, 50, 991–996. [Google Scholar] [CrossRef]
  18. Vajravelu, K.; Rollins, D. Hydromagnetic flow of a second grade fluid over a stretching sheet. Appl. Math. Comput. 2004, 148, 783–791. [Google Scholar] [CrossRef]
  19. Pal, D.; Mondal, H. Effects of temperature-dependent viscosity and variable thermal conductivity on MHD non-Darcy mixed convective diffusion of species over a stretching sheet. J. Egypt. Math. Soc. 2014, 22, 123–133. [Google Scholar] [CrossRef]
  20. Yazdi, M.H.; Hashim, I.; Sopian, K. Bin Slip boundary layer flow of a power-law fluid over moving permeable surface with viscous dissipation and prescribed surface temperature. Int. Rev. Mech. Eng. 2014, 8, 502–510. [Google Scholar]
  21. Hsiao, K.-L. MHD mixed convection of viscoelastic fluid over a stretching sheet with ohmic dissipation. J. Mech. 2008, 24, 29–34. [Google Scholar] [CrossRef]
  22. Vajravelu, K.; Prasad, K.V.; Vaidya, H. Influence of hall current on MHD flow and heat transfer over a slender stretching sheet in the presence of variable fluid properties. Commun. Numer. Anal. 2016, 2016, 17–36. [Google Scholar] [CrossRef]
  23. Vajravelu, K.; Prasad, K.V.; Vaidya, H.; Basha, N.Z.; Ng, C.-O. Mixed convective flow of a Casson fluid over a vertical stretching sheet. Int. J. Appl. Comput. Math. 2017, 3, 1619–1638. [Google Scholar] [CrossRef]
  24. Ali, F.; Sheikh, N.A.; Khan, I.; Saqib, M. Influence of a porous medium on the hydromagnetic free convection flow of micropolar fluid with radiative heat flux. J. Porous Media 2018, 21, 123–144. [Google Scholar] [CrossRef]
  25. Hsiao, K.-L. Conjugate heat transfer of magnetic mixed convection with radiative and viscous dissipation effects for second-grade viscoelastic fluid past a stretching sheet. Appl. Therm. Eng. 2007, 27, 1895–1903. [Google Scholar] [CrossRef]
  26. Pop, I.; Isa, S.S.P.M.; Arifin, N.M.; Nazar, R.; Bachok, N.; Ali, F.M. Unsteady viscous MHD flow over a permeable curved stretching/shrinking sheet. Int. J. Numer. Methods Heat Fluid Flow 2016, 26, 2370–2392. [Google Scholar] [CrossRef]
  27. Reddy, J.V.R.; Sugunamma, V.; Sandeep, N. Dual solutions for nanofluid flow past a curved surface with nonlinear radiation, Soret and Dufour effects. J. Phys. Conf. Ser. 2018, 1000, 12152. [Google Scholar] [CrossRef] [Green Version]
  28. Rudraswamy, N.G.; Kumar, K.G.; Gireesha, B.J.; Gorla, R.S.R. Soret and Dufour effects in three-dimensional flow of Jeffery nanofluid in the presence of nonlinear thermal radiation. J. Nanoeng. Nanomanuf. 2016, 6, 278–287. [Google Scholar] [CrossRef]
  29. Roşca, N.C.; Pop, I. Unsteady boundary layer flow over a permeable curved stretching/shrinking surface. Eur. J. Mech. 2015, 51, 61–67. [Google Scholar] [CrossRef]
  30. Rashad, A.M.; Rashidi, M.M.; Lorenzini, G.; Ahmed, S.E.; Aly, A.M. Magnetic field and internal heat generation effects on the free convection in a rectangular cavity filled with a porous medium saturated with Cu–water nanofluid. Int. J. Heat Mass Transf. 2017, 104, 878–889. [Google Scholar] [CrossRef]
  31. Sheikholeslami, M.; Vajravelu, K.; Rashidi, M.M. Forced convection heat transfer in a semi annulus under the influence of a variable magnetic field. Int. J. Heat Mass Transf. 2016, 92, 339–348. [Google Scholar] [CrossRef]
  32. Sheikholeslami, M.; Rashidi, M.M.; Ganji, D.D. Effect of non-uniform magnetic field on forced convection heat transfer of Fe3O4–water nanofluid. Comput. Methods Appl. Mech. Eng. 2015, 294, 299–312. [Google Scholar] [CrossRef]
  33. Shu, C. Differential Quadrature and Its Application in Engineering; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  34. Baskaya, E.; Komurgoz, G.; Ozkol, I. Investigation of oriented magnetic field effects on entropy generation in an inclined channel filled with ferrofluids. Entropy 2017, 19, 377. [Google Scholar] [CrossRef]
Figure 1. (a) Geometry of the flow domain. (b) Thermal and momentum boundary layer.
Figure 1. (a) Geometry of the flow domain. (b) Thermal and momentum boundary layer.
Entropy 20 00943 g001
Figure 2. Impacts on velocity profile g’(χ) with variations in (a) magnetic parameter M and (b) curvature parameter κ .
Figure 2. Impacts on velocity profile g’(χ) with variations in (a) magnetic parameter M and (b) curvature parameter κ .
Entropy 20 00943 g002
Figure 3. Impacts on temperature profile θ(χ)with variations in (a) magnetic parameter M (b) curvature parameter κ (c) Eckert number Ec (d) Prandtl number Pr and (e) variable thermal conductivity parameter τ .
Figure 3. Impacts on temperature profile θ(χ)with variations in (a) magnetic parameter M (b) curvature parameter κ (c) Eckert number Ec (d) Prandtl number Pr and (e) variable thermal conductivity parameter τ .
Entropy 20 00943 g003aEntropy 20 00943 g003b
Figure 4. Impacts on entropy generation number Ns with variations in (a) magnetic parameter M (b) curvature parameter κ (c) Eckert number EC (d) Prandtl number Pr (e) variable thermal conductivity parameter τ and (f) temperature difference parameter λ.
Figure 4. Impacts on entropy generation number Ns with variations in (a) magnetic parameter M (b) curvature parameter κ (c) Eckert number EC (d) Prandtl number Pr (e) variable thermal conductivity parameter τ and (f) temperature difference parameter λ.
Entropy 20 00943 g004aEntropy 20 00943 g004b
Table 1. Present numerical results for the skin friction coefficient Res0.5Cfs and the rate of heat transfer Res−0.5Nus at the curved surface by GDQM and RKFM, for various values of the physical parameters M, κ, Ec, Pr and τ.
Table 1. Present numerical results for the skin friction coefficient Res0.5Cfs and the rate of heat transfer Res−0.5Nus at the curved surface by GDQM and RKFM, for various values of the physical parameters M, κ, Ec, Pr and τ.
MκEcPrτ*GDQM*RKFM
Re s 0.5 C f s Re s 0.5 N u s Re s 0.5 C f s Re s 0.5 N u s
0.0100.32.00.21.07348861.09563461.07348861.0956346
0.51.32798491.01829021.32798491.0182902
1.01.53029130.94337631.53029130.9433763
2.01.86012860.79560161.86012860.7956016
3.02.13384600.64953662.13384600.6495366
0.250.32.00.21.28565251.05802251.28565261.0580225
101.18465731.06414281.18465731.0641428
201.13862921.06593531.13862921.0659353
301.12393411.06634821.12393411.0663482
10001.09632011.06689151.09632011.0668915
0.2100.12.00.21.18465731.11769211.18465731.1176921
0.41.18465731.03393801.18465731.0339380
0.71.18465730.92955821.18465730.9295582
1.01.18465730.80445341.18465730.8044534
1.51.18465730.54962361.18465730.5496236
0.2100.31.00.21.18465730.82214391.18465730.8221439
2.01.18465731.06414281.18465731.0641428
3.01.18465731.18013811.18465731.1801381
5.01.18465731.24997801.18465731.2499781
6.01.18465731.23918671.18465731.2391866
0.2100.32.00.01.18465731.73569481.18465731.7356948
0.31.18465730.82011821.18465730.8201182
0.71.18465730.17595901.18465730.1759590
1.01.1846573−0.11209861.1846573−0.1120984
1.51.1846573−0.41421811.1846573−0.4142179
*GDQM: Generalized differential quadrature method; *RKFM: Runge-Kutta-Fehlberg method

Share and Cite

MDPI and ACS Style

Afridi, M.I.; Wakif, A.; Qasim, M.; Hussanan, A. Irreversibility Analysis of Dissipative Fluid Flow Over A Curved Surface Stimulated by Variable Thermal Conductivity and Uniform Magnetic Field: Utilization of Generalized Differential Quadrature Method. Entropy 2018, 20, 943. https://doi.org/10.3390/e20120943

AMA Style

Afridi MI, Wakif A, Qasim M, Hussanan A. Irreversibility Analysis of Dissipative Fluid Flow Over A Curved Surface Stimulated by Variable Thermal Conductivity and Uniform Magnetic Field: Utilization of Generalized Differential Quadrature Method. Entropy. 2018; 20(12):943. https://doi.org/10.3390/e20120943

Chicago/Turabian Style

Afridi, Muhammad Idrees, Abderrahim Wakif, Muhammad Qasim, and Abid Hussanan. 2018. "Irreversibility Analysis of Dissipative Fluid Flow Over A Curved Surface Stimulated by Variable Thermal Conductivity and Uniform Magnetic Field: Utilization of Generalized Differential Quadrature Method" Entropy 20, no. 12: 943. https://doi.org/10.3390/e20120943

APA Style

Afridi, M. I., Wakif, A., Qasim, M., & Hussanan, A. (2018). Irreversibility Analysis of Dissipative Fluid Flow Over A Curved Surface Stimulated by Variable Thermal Conductivity and Uniform Magnetic Field: Utilization of Generalized Differential Quadrature Method. Entropy, 20(12), 943. https://doi.org/10.3390/e20120943

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