Document Type : Full Lenght Research Article
Authors
1 University of Gour Banga, Malda 732 103, WB, India
2 Department of Applied Mathematics, Vidyasagar University, Midnapore 721 1102, India
3 Faculty of Military Science, Stellenbosch University, Private Bag X2, Saldanha 7395, South Africa
Abstract
Keywords
Main Subjects
Entropy generation in hydromagnetic and thermal boundary layer flow due to radial stretching sheet with Newtonian heating
Sanatan Das1* , Rabindra Nath Jana2 and Oluwole Daniel Makinde3
I*Department of Mathematics, University of Gour Banga, Malda 732 103, India 2Department of Applied Mathematics, Vidyasagar University, Midnapore 721 102, India 3Faculty of Military Science, Stellenbosch University, Private Bag X2, Saldanha 7395, South Africa
Corresponding author e-mail: tutusanasd@yahoo.co.in
|
Journal of Heat and Mass Transfer Research
Journal homepage: http://jhmtr.journals.semnan.ac.ir |
|
A B S T R A C T The entropy generation during hydromagnetic boundary layer flow of a viscous incompressible electrically conducting fluid due to radial stretching sheet with Newtonian heating in the presence of a transverse magnetic field and the thermal radiation has been analyzed. The governing equations are then solved numerically by using the fourth order Runge-Kutta method with shooting technique. The effects of the pertinent parameters on the fluid velocity, temperature, entropy generation number, Bejan number, as well as the shear stress at the surface of the sheet are discussed graphically and quantitatively. It is examined that because the presence of a magnetic field, the entropy generation in a thermal system can be controlled and reduced.
© 2016 Published by Semnan University Press. All rights reserved.
|
History: Received 3 January 2014 Received in revised form 4 June 2014 Accepted 29 September 2015
Keywords: Hydromagnetic flow, Newtonian heating, Entropy generation, radial stretching sheet |
|
1. Introduction
The boundary layer flow of a continuously stretching sheet has attracted considerable attention in recent years due to its numerous applications in industry. It occurs frequently in manufacturing involving the hot metal rolling, wire drawing, glass-fiber production, paper production, drawing of plastic films, and metal spinning, as well as metal and polymer extrusion processes. The boundary layer flow and the heat transfer over a continuous stretching surface have been discussed by a significant number of researchers [1-7] due to its wide applications in industrial and engineering processes. However, less attention has been paid to the boundary layer flow over a radial stretching sheet. Hayat et al. [8] have studied an axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet. Sahoo and Sharma [9] have investigated the MHD flow and heat transfer from a continuous surface in a uniform free stream of a non-Newtonian fluid. Sajid et al. [10] have presented an unsteady axisymmetric flow and heat transfer over a radial stretching sheet. Ariel [11] has presented the extended homotopy perturbation method and computation of flow past a stretching sheet. Sahoo [12] has examined the effects of the partial slip on axisymmetric flow of an electrically conducting viscoelastic fluid past a stretching sheet. Effects of slip, viscous dissipation and Joule heating on the MHD flow and heat transfer of a second grade fluid past a radial stretching sheet has been investigated by Sahoo [13].
The second law of thermodynamics is more reliable than the first law of thermodynamics due to the limitation of the efficiency of the first law in heat transfer engineering systems (Oztop and Al-Salem [14]). In order to access the best design of the thermal systems, one can employ the second law of thermodynamics by minimizing the irreversibility. Entropy generation is a criterion of the destruction of the available system work. The evaluation of the entropy generation is carried out to improve system performance. Heat transfer, mass transfer, viscous dissipation, etc. can be used as sources of entropy generation. Entropy generation can be used as a quantitative measure of irreversibilities that are associated with a process, because of this fact that the greater the entropy generation indicates the greater the extent of irreversibilities. In many engineering and industrial processes, entropy production, destroys the available energy in the system. It is therefore imperative to determine the rate of entropy generation in a system, in order to optimize energy in the system for efficient operation in the system. According to the second law of thermodynamics, all the flows and heat transfer processes undergo changes that are irreversible. These irreversible changes are mostly caused by the energy losses during the processes. Although measures can be taken to reduce these irreversible effects, it is impossible to recover all of the lost energy. This process causes the entropy of the system to increase. Due to this, the entropy generation rate is used as a standard metric to study the irreversibile effects. This method was proposed by Bejan [15,16].
In recent years, many papers have been published on the applications and entropy generation rates of the second law of thermodynamics. Odat et al. [17] have explored the entropy generation effects in the laminar flow past a flat plate under the influence of the magnetic field and found that the entropy generation rate increases with the magnetic field intensity. This study has revealed that the magnetic field is one of the causes that are responsible for the entropy production in the systems. Saouli and Aiboud-Saouli [18] investigated the second law analysis of laminar falling liquid film along an inclined heated plate. Esfahani and Jafarian [19] have presented the entropy generation analysis of a flat plate boundary layer with different solution techniques. An irreversibility analysis for gravity driven non-Newtonian liquid film along an inclined isothermal plate has been presented by Makinde [20]. Arikoglu et al. [21] have examined the effect of slip on entropy generation in a single rotating disk in MHD flow. Aiboud and Saouli [22] have illustrated the application of the second law analysis of thermodynamics to viscoelastic magneto-hydrodynamic flow over a stretching surface analytically by using Kummer's functions. Makinde [23] has conducted a thermodynamic analysis on a gravity-driven liquid film along an inclined heated plate. He assumed viscosity to be a variable quantity and considered the convective cooling effect. Makinde [24] has studied a variable viscosity boundary layer flow over a flat plate under the effects of thermal radiation and Newtonian heating, as well as explored entropy generation effects in this flow. Makinde [25] has examined the entropy generation on magnetohydrodynamic flow and heat transfer over a flat plate with a convective boundary condition. The effect of viscous dissipation and thermal radiation on entropy generation in Blasius flow has been displayed numerically by Butt et al. [26]. Rashidi et al. [27] have analyzed the entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. The results showed that as the thermal radiation parameter increases, the generated entropy decreases. Butt and Ali [28] have illustrated the effects of magnetic field on entropy generation in the flow and heat transfer due to a radially stretching surface. Butt and Ali [29,30] have carried out the entropy analysis of flow and heat transfer caused by a moving surface. Recently, Butt and Ali [31] have presented the entropy analysis of magnetohydrodynamic flow and heat transfer over a convectively heated radially stretching surface.
The aim of the present paper is to explore the entropy generation in MHD boundary layer flow of a viscous incompressible electrically conducting fluid due to a radial stretching sheet in the presence of a transverse magnetic field by taking the convective boundary condition into account. The partial differential equations governing the flow are reduced to nonlinear ordinary differential equations, which are solved numerically by shooting technique by using fourth order Runge-Kutta method. A systematic study of the effects of the various pertinent parameters on the flow and heat transfer characteristic is carried out with the help of graphs and table.
2. Mathematical formulation
Consider a steady two-dimensional boundary layer flow due to the stretching of the sheet along the radial direction with the velocity , where is a constant, as shown in Fig. 1. The sheet is located in the plane and the fluid is confined to the region . The flow is due to the stretching of the sheet. A uniform magnetic field of strength is applied perpendicularly to the sheet, i.e. in the -direction. The lower surface of the sheet is heated by convection from a hot fluid at temperature which provides a heat transfer coefficient while the temperature of the ambient cold fluid is . The cold fluid on the upper side of the sheet is assumed to be an electrically conducting Newtonian fluid with constant fluid property. There is no external electric field. Also, it is assumed that the magnetic Reynolds number is small enough so that the induced magnetic field can be neglected.
An order-of-magnitude analysis of the momentum equation (normal to the sheet) by using the usual boundary layer approximations shows that the pressure gradient is constant. Thus neglecting the pressure gradient, the continuity, momentum and energy equations in a viscous MHD incompressible boundary layer flow can be written, respectively, as follows
(1) |
|
(2) |
|
(3) |
Where and are the velocity components along the and -directions, respectively, the temperature of the fluid, the electric conductivity of the fluid, the kinematic viscosity, the fluid density, the thermal conductivity, the specific heat at constant pressure and the radiative heat flux.
The appropriate boundary conditions are
(4) |
Where is the thermal conductivity.
For an optically thick fluid, in addition to emission there is also self absorbed and usually the absorption coefficient is large and dependent on the wavelength so that we can apply Rosseland approximation for radiative flux. The Rosseland approximation [32] applies to optically thick medium and taking into account the model proposed by Magyari and Pantokratoras [33], the net radiation heat flux is given by
(5) |
Where is the Stefan-Boltzmann constant and the Rosseland mean absorption coefficient. The Rosseland mean absorption coefficients for the optically thick media are expressed as function of the thermodynamic properties of the media. It is assumed that the temperature difference between the fluid temperature and the free stream temperature is small, so that the term may be expressed as a linear function of temperature. This is done by expanding in a Taylor series about a free stream temperature as follows:
(6) |
Neglecting higher-order terms in the equation (6) beyond the first order in we get
(7) |
The use of the equation (5) and (7), the equation (3) becomes
(8) |
The following similarity variables are introduced (Butt and Ali [28]):
(9) |
Fig.1: Geometry of the problem |
where is the independent similarity variable, the dimensionless stream function and the dimensionless temperature.
The use of equation (9), equations (2) and (8) reduce to
(10) |
|
(11) |
Where is the magnetic parameter, the radiation parameter, and the Prandtl number which measures the ratio of momentum diffusivity to the thermal diffusivity.
The thermal radiation is quite significant and the quality of the final product can be controlled by the control of cooling rate via the radiation parameter. In polymer industry, the thermal radiation effect may play an important role in the control of heat transfer process if the process is directed in a thermally controlled environment.
The corresponding boundary conditions are
(12) |
Where is the surface convection parameter or so-called Biot number. When , the convective boundary condition reduces to a uniform surface temperature boundary condition. For uniformly heated ( ) radial stretching sheet and in the absence of thermal radiation ( ) the present problem reduces to the Butt and Ali [28].
3. Numerical method for solution
The governing non-linear ordinary differential equations (10) and (11) cannot be solved analytically. This set of equations with the boundary conditions (12) must be solved numerically by applying most efficient fourth-order Runge-Kutta integration scheme with shooting algorithm. Equations (10) and (11) and boundary conditions (12) are reduced to a set of simultaneous first order differential equations by setting , , , and as follows
(13) |
With the boundary conditions
(14) |
Where and are determined such that and . The essence of this method is to reduce the boundary value problem to an initial value problem and then to use the shooting numerical technique to guess and until the boundary conditions and are satisfied. It is important to note that the infinity ( ) in the above equations represents the boundary layer thickness. We compare the calculated values for and at with the given boundary condition and and adjust the estimated values, and , to give a better approximation for the solution. The accuracy of the assumed missing initial condition is checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. The step size is used to obtain the numerical solution with seven-decimal place accuracy as the criterion of convergence. The numerical computations are done by MATLAB built-in routine. The method is adequately explained in literature and it has second order convergence, unconditionally stable. It gives accurate result for boundary layer equations. In the present study, a uniform grid is used which is concentrated toward the wall.
4. Results and discussion
In order to gain a clear physical insight of the problem, we have discussed the effects of different values of magnetic parameter , radiation parameter , Prandtl number and the Biot number on the velocity, temperature and shear stress at the surface of the sheet. The default values of the other parameters are mentioned in the description of the respected figures. Fig.2 shows that the fluid velocity decreases with an increase in magnetic parameter . A drag-like Lorentz force is created by the application of the transverse magnetic field of the electrically conducting fluid. This force has the tendency to slow down the fluid flow. Fig.3 shows that the fluid temperature increases by increasing the magnetic parameter in the boundary layer region and the thermal boundary layer thickness increases as the magnetic field becomes stronger. It is important to notice that large resistances on the fluid particles, which cause heat to be generated in the fluid, as the transverse applied magnetic field increases. It can be observed from Fig.4 that the fluid temperature increases by increasing values of radiation parameter . The increase in radiation parameter means the release of heat energy from the flow region and so the fluid temperature decreases. A decrease in the values of for given and means a decrease in the Rosseland radiation absorptivity . Since the divergence of the radiative heat flux increases, decreases which in turn causes to increase the rate of radiative heat transfer to the fluid and hence the fluid temperature increases. In the presence of thermal radiation, the thermal boundary layer always found to be increased. This means that the thermal boundary layer increases and there is more uniform temperature distribution across the boundary layer.
The effect of Prandtl number Pr on the heat transfer process is shown in the Fig.5. This figure shows that an increase in Prandtl number number Pr results in a decrease in the temperature distribution because, the thermal boundary layer thickness decreases by increasing the Prandtl number Pr. In summary, an increase in the Prandtl number means slow rate of the thermal diffusion. The graph also shows that the wall temperature decreases, as the values of Prandtl number Pr increase. It is also found that the thermal boundary layer thickness reduces as Pr increases. Fig.6 demonstrates the effects of the Biot number on fluid temperature. It can be seen that the temperature profiles within the boundary layer increase by increasing in Biot number, leading to an increase in the thermal boundary layer thickness. Biot number is the ratio of the hot fluid side convection resistance to the cold fluid side convection resistance of a surface. For fixed cold fluid properties, Biot number is directly proportional to the heat transfer coefficient associated with the hot fluid. The thermal resistance on the hot fluid side is inversely proportional to . Thus, as increases, the hot fluid side convection resistance decreases and consequently, the surface temperature increases. It is also noticed that for large values of , i.e. , the temperature profile attains its maximum value 1; thus the convective boundary condition becomes the prescribed surface temperature case. The thermal boundary layer thickness increases as increases.
Fig.2: Velocity profiles for different
|
Fig.3: Temperature profiles for different when , and
|
Fig.4: Temperature profiles for different when , and
|
Fig.5: Temperature profiles for different when , and
|
Fig.6: Temperature profiles for different when , and |
At the engineering point of view, the shear stress at the surface of the sheet is an important characteristic in the heat transfer studies, since it is directly related to the heat transfer coefficients. The increased shear stress is generally a disadvantage in the technical applications, while the increased heat transfer can be exploited in some applications such as heat exchangers, but should be avoided in other applications such as a gas turbine, for instance. The numerical values of the surface temperature and the shear stress at the surface of the sheet are entered in the Table 1 for several values of , , Pr and . It can be seen from the Table 1 that the surface temperature increases by increasing values of either or or while it decreases by increasing in . This is because of the fact that the Biot number is directly proportional to the heat transfer coefficient associated with hot fluid. The thermal resistance on the hot fluid side is inversely proportional to . Thus, as increases, the hot fluid side convection resistance decrease; as a result, the surface temperature increases. The shear stress increases by increasing values of .
4.1. Entropy generation
Non-equilibrium conditions happen due to the exchange of the energy and momentum within the fluid and at solid boundaries that result in entropy generation. According to Woods [34], the local volumetric rate of entropy generation for a viscous incompressible conducting fluid in the presence of a magnetic field is given by
(15) |
The first term in equation (15) is the irreversibility due to the heat transfer, the second term is the entropy generation due to viscous dissipation and the third term is a local entropy generation due to the effect of magnetic field (Joule
The dimensionless entropy generation number may be defined by the following relationship:
(16) |
The use of (9), the entropy generation number in dimensionless form is
(17) |
Where is the Reynolds number, is the Brinkmann number and the non-dimensional temperature difference.
The entropy generation number can be written as a summation of the entropy generation due to heat transfer denoted by and the entropy generation due to fluid friction with magnetic field denoted by given as
(18) |
In order to obtain an idea of whether entropy generation due to the heat transfer dominates over entropy generation due to the fluid friction and the magnetic field, or vice versa, the Bejan number is defined to be the ratio of entropy generation due to the heat transfer to the entropy generation number [35].
(19) |
Where is the irreversibility ratio. Heat transfer dominates for and fluid friction with magnetic effects dominates when . The contribution of both heat transfer and fluid friction to entropy generation is equal when . The Bejan number takes the values between 0 and 1 (see Cimpean et al. [36]). The value of is the limit at which the heat transfer irreversibility dominates, is the opposite limit at which the irreversibility is dominated by the combined effects of fluid friction and magnetic field and is the case in which the heat transfer and fluid friction with magnetic field entropy production rates are equal. Further, the behavior of the Bejan number is studied for the optimum values of the parameters at which the entropy generation takes its minimum.
The influences of the different governing parameters on entropy generation and Bejan number are presented in Figs. 7-14. It can be seen from Fig.7 that the entropy generation number increases by increasing in magnetic parameter . An increase in the magnetic field intensity causes an increase in the entropy generation. It reveals that the magnetic field is a source of entropy generation in addition to the fluid friction and the heat transfer. Also, it can be seen that the entropy generation effects are prominent on the surface of the stretching sheet and in the region close to it. However, in the free stream region the entropy effects are negligible. This implies that in order to control the entropy which is generated in boundary layer flow, the value of the magnetic parameter should be reduced, which is an issue of interest in nuclear-MHD propulsion. Fig.8 shows that the entropy generation number increases by increasing Reynolds number due to higher heat transfer rates at the sheet surface. An increase in Reynolds number , the entropy effects due to heat transfer entropy effects become prominent and fluid friction and magnetic field are lessen near the stretching sheet surface. However, as the distance increases from the sheet surface, these effects are negligible. Fig.9 illustrates the effects of the Biot number on the entropy generation. Increasing enhances the entropy generation. It means that the convective surface boundary condition acts as a strong source of irreversibility. Therefore, the entropy can be minimized by reducing the convection through the boundary. The closeness of the curves in Fig.9 can be attributed that the entropy effects are dominated by the entropy effects due to fluid friction and magnetic field. It can be observed from Fig.10 that entropy generation number increases by increasing the group parameter due to the viscous heating effects. An increase in the values of the group parameter due to the combined effects of viscous heating and temperature difference yields a higher entropy generation number.
Fig.11 shows that the Bejan number increases by increasing the values of magnetic parameter . For large values of , the entropy effects due to fluid friction and magnetic field are fully dominated by heat transfer entropy effects near the sheet surface. Fig.12 illustrates the effects of the radiation parameter on the Bejan number . By increasing values of , the entropy effects due to heat transfer become strong and hence Bejan number increases. This is due to the fact that the effect of the radiation parameter is to enhance the temperature significantly in the flow region. The increase in radiation parameter means the generation of heat energy in the flow region. An increase in Biot number leads to increase the Bejan number as shown in Fig.13. Also, an increase in the values of the Biot number results in an increase in the dominant effect of heat transfer irreversibility at the sheet surface. This means that the sheet surface acts as a strong source of irreversibility. Fig.14 reveals that the Bejan number decreases by increasing group parameter . This is quite true as higher values of , which increase the magnitude of fluid friction with magnetic field irreversibility but has no effect on the heat transfer irreversibility , increases the values of leading to lower Bejan number. The group parameter is an important dimensionless number for irreversibility analysis. It determines the relative importance of viscous effects to that of temperature gradient entropy generation. The graphs of the Bejan number are useful to obtain an idea on whether heat transfer irreversibility dominates fluid friction irreversibility or vice versa.
5. Conclusion
Entropy generation analysis in hydromagnetic boundary layer flow of a viscous incompressible electrically conducting fluid due to a radially stretching sheet with Newtonian heating in the presence of a transverse magnetic field has been carried out. The velocity and temperature profiles are obtained numerically and used to compute the entropy generation number. The effects of the pertinent parameters on velocity and temperature profiles are presented graphically.
The influences of the same parameters and the dimensionless group parameter on the entropy generation rate and Bejan number are also discussed.
From the results the following conclusions could be drawn:
Table 1. The shear stress and the temperature at the sheet .
|
|
||||||||||||
0.5 |
|
0.72 |
0.1 |
1.57758 |
0.08345 |
||||||||
1 |
|
0.72 |
0.1 |
1.70391 |
0.08358 |
||||||||
2 |
|
0.72 |
0.1 |
1.94118 |
0.08384 |
||||||||
|
|
||||||||||||
1 |
0.1 |
|
|
|
0.08018 |
||||||||
1 |
0.5 |
|
|
|
0.08283 |
||||||||
1 |
1 |
|
|
|
0.08474 |
||||||||
|
|
||||||||||||
|
|
0.72 |
|
|
0.08358 |
||||||||
|
|
2 |
|
|
0.07279 |
||||||||
|
|
3 |
|
|
0.06597 |
||||||||
|
|
||||||||||||
|
|
|
0.1 |
|
0.08358 |
||||||||
|
|
|
1 |
|
0.47701 |
||||||||
|
|
|
5 |
|
0.82016 |
||||||||
Fig.7: for different when , and |
Fig.8: for different when , and |
Fig.9: for different when , and |
Fig.10: for different when , and |
Fig.11: Bejan number for different when , and |
Fig.12: Bejan number for different when , and |
Fig.13: Bejan number for different when , and |
Fig.14: Bejan number for different when and |
Acknowledgement
The authors are thankful to the honourable reviewers for constructive suggestions which helped a lot in improving the manuscript.
Nomenclature
|
a constant |
|
applied magnetic field |
|
Bejan number |
|
Biot number |
|
Brinkmann number |
|
specific heat at constant pressure |
|
volumetric rate of entropy generation |
|
non-dimensional steam function |
|
first order derivative with respect to |
|
second order derivative with respect to |
|
third order derivative with respect to |
|
heat transfer coefficient |
|
thermal conductivity |
|
Rosseland mean absorption coefficient |
|
magnetic parameter |
|
entropy generation due to heat transfer |
|
entropy generation due to fluid friction and magnetic field |
|
entropy generation number |
|
Prandtl number |
|
radiative heat flux |
|
radiation parameter |
|
Reynolds number |
|
fluid temperature |
|
hot fluid temperature |
|
free stream temperature |
, |
velocity components in and -directions |
|
Stretching velocity of the sheet in radial direction |
|
Greek symbols |
|
similarity variable |
|
kinematic viscosity |
|
irreversibility distribution ratio |
|
the fluid density |
|
Stefan-Boltzman constant |
|
non-dimensional temperature difference |
|
non-dimensional temperature |
|
first order derivative with respect to |
|
second order derivative with respect to |
References
[1] Carragher, P., Crane, L.J.(1982). Heat transfer on a continuous stretching sheet. Z. Angew Math. Mech. 62, 564-565.
[2] Wang, C.Y. (1990). Liquid film on an unsteady stretching sheet. Q. Appl. Math. 48, 601-610.
[3] Ariel, P. D. (2001). Axisymmetric flow of a second grade fluid past a stretching sheet. Int. J. Eng. Sci. 39(5), 529-553.
[4] Chiam, T.C. (1994). Stagnation point flow towards a stretching plate. J. Phys. Soc. Jpn. 63, 2443-2444.
[5] Liu, I. C. (2004). Flow and heat transfer of an electrically conducting fluid of second grade over a stretching sheet subject to a transverse magnetic field. Int. J. Heat Mass Transfer 47(19–20), 4427-4437.
[6] Wang, C.Y . (2007). Natural convection on a vertical radially stretching sheet. Math. Anal. App. 332, 877- 883.
[7] Ariel, P.D. (2007). Axisymmetric flow due to a stretching sheet with partial slip. Comput. Math. Appl. 54, 1169-1183.
[8] Hayat, T., Sajid, M. (2007). Analytic solution for axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet. Int. J. Heat Mass Transf 50, 75-84.
[9] Sahoo, B., Sharma, H. G. (2007). MHD flow and heat transfer from a continuous surface in a uniform free stream of a non-Newtonian fluid. Appl. Math. Mech.-Eng. Ed. 28(11),1467-1477.
[10] Sajid, M., Ahmad, I., Hayat, T., Ayub, M. (2008). Series solution for unsteady axisymmetric flow and heat transfer over a radially stretching sheet. Commun. Nonlinear Sci. Numer. Simul. 13(10), 2193-2202.
[11] Ariel, P.D. (2009). Extended homotopy perturbation method and computation of flow past a stretching sheet. Comput. Math. Appl. 58, 2402-2409.
[12] Sahoo, B. (2009). Effects of partial slip on axisymmetric flow of an electrically conducting viscoelastic fluid past a stretching sheet. Cent. Eur. J. Phys., DOI; 10.2478/s11534-009-0105-x.
[13] Sahoo, B. (2010). Effects of slip, viscous dissipation and joule heating on the MHD flow and heat transfer of a second grade fluid past a radially stretching sheet. Appl. Math. Mech.- Eng. Ed. 31, 159-173.
[14] Oztop, H. F., Al-Salem, K. A. (2012). review on entropy generation in natural and mixed convection heat transfer for energy systems, Renewable Sustainable Energy Rev. 16 (1), 911-920.
[15] Bejan, A. (1980). Second law analysis in heat transfer. Energy Int. J. 5, 721-732.
[16] Bejan, A . (1994). Entropy Generation Through Heat and Fluid Flow, Wiley, Canada.
[17] Odat, M. Q. A., Damseh, R. A., Nimr, M. A. A. (2004). Effect of magnetic field on entropy generation due to laminar forced convection past a horizontal flat plate. Entropy 4, 293-303.
[18] Saouli, S., Aiboud-Saouli, S. (2004). Second law analysis of laminar falling liquid film along an inclined heated plate. Int. Comm. Heat Mass Transfer 31(6), 879-886.
[19] Esfahani, J. A., Jafarian, M. M. (2005). Entropy generation analysis of a flat plate boundary layer with various solution methods. Scientia Iranica 12, 233-240.
[20] Makinde, O. D. (2006). Irreversibility analysis for gravitry driven non-Newtonian liquid film along an inclined isothermal plate. Physica Scripta 74, 642-645.
[21] Arikoglu, A., Ozkol. I., Komurgoz, G. (2008).Effect of slip on entropy generation in a single rotating disk in MHD flow. Appl. Energy 85(12), 1225-1236.
[22] Adboud, S., Saouli, S. (2010).Entropy analysis for viscoelastic magneto hydrodynamic flow over a stretching surface. Int. J. Non Linear Mech. 45 (5), 482-489.
[23] Makinde, O. D. (2010).Thermodynamic second law analysis for a gravity driven variable viscosity liquid film along an inclined heated plate with convective cooling. J. Mech. Sci. Tech. 24 (4), 899-908.
[24] Makinde, O. D. (2011).Second law analysis for variable viscosity hydromagnetic boundary layer flow with thermal radiation and Newtonian heating, Entropy 13, 1446-1464.
[25] Makinde, O. D. (2012).Entropy analysis for MHD boundary layer flow and heat transfer over a flat plate with a convective surface boundary condition. Int. J. Exergy 10 (2), 142 - 154.
[26] Butt, A.S., Munawar, S., Ali, A., Mehmood, A. (2012).Entropy generation in the Blasius flow under thermal radiation. Phys. Scr. 85 (3), 035008.
[27] Rashidi, S., Abelman, N., Freidooni, M. (2013). Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid, Int. J. Heat and Mass Transfer 62, 515-525.
[28] Butt, A.S., Ali, A. (2013). Effects of magnetic field on entropy generation in flow and heat transfer due to a radially stretching surface. Chin. Phys. Lett. 30(2), 024701.
[29] Butt, A.S., Ali, A.( 2014). Entropy analysis of flow and heat transfer caused by a moving plate with thermal radiation. J. Mech. Sci. Tech. 28 (1), 343-348.
[30] 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 129(3), 1-13.
[31] Butt, A.S., Ali, A.( 2014). Entropy analysis of magnetohydrodynamic flow and heat transfer over a convectively heated radially stretching surface. J. Taiwan Institute of Chem. Eng. 45, 1197-1203.
[32] Rosseland, S. (1931). Astrophysik und atom-theoretische Grundlagen, Springer-Verlag, Berlin.
[33] Magyari, E., Pantokratoras, A. (2011). Note on the effect of thermal radiation in the linearized Rosseland approximation on the heat transfer characteristics of various boundary layer flows. Int. Commun. Heat and Mass Trans. 38, 554-556.
[34] Woods, L. C. (1975). Thermodynamics of Fluid Systems. Oxford University Press, Oxford, UK.
[35] Paoletti, S, Rispoli, F, Sciubba, E.( 1989). Calculation of exergetic losses in compact heat exchanger passages. ASME AES 10, 21-29.
[36] Cimpean, D., Lungu, N., Pop, I. (2008). A problem of entropy generation in a channel filled with a porous medium. Creative Math. Inf. 17, 357-362.