Journal of Applied Mathematics and Physics
Vol.04 No.07(2016), Article ID:69242,8 pages

Analysis of Entropy Generation in a Rectangular Porous Duct

Shyam Lal Yadav, Ashok Kumar Singh

Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi, India

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 3 May 2016; accepted 25 July 2016; published 28 July 2016


In this paper, we have considered a fully developed flow of a viscous incompressible fluid in a rectangular porous duct saturated with the same fluid. The duct is heated from the bottom for forced and mixed convection. The Brinkman model is used to simulate the momentum transfer in the porous duct. Using the momentum and thermal energy equations, the entropy generation has been obtained due to the heat transfer, viscous and Darcy dissipations. It is found from the mathematical analysis that the entropy generation is double when the viscous as well as the Darcy dissipations terms are taken in the thermal energy equation in comparison when the viscous as well as the Darcy dissipations terms are not taken in the thermal energy equation. This result clearly shows that there is no need of taking the viscous and Darcy dissipations terms in the thermal energy equation to obtain the entropy generation.


Brinkman Model, Heat Flux, Fully Developed Flow, Viscosity Ratio, Entropy Generation

1. Introduction

The study of transport phenomena in the channel filled with saturated porous media has attracted considerable attention of scientists, engineers and experimentalists in present time. This attention is mainly due to the applications of this phenomenon in the field of electronics cooling system, geothermal system, storage of nuclear waste materials, microelectronics heat transfer equipment, coal and grain storage, crude oil production, catalytic converters, ground water pollution, fiber and granular insulations, solidification of castings, etc. The advancement in the thermal systems as well as the energy utilization during the convection in any fluid is one of the fundamental problems of the technological processes, because the improved thermal systems will provide better material processing, energy conservation and environmental effects. Also, because of applications to the cooling of electronic equipment there has been an increased interest in the forced convections in the channels and ducts filled with the porous media. One of the important viscous fluid flow situations in the porous media is the Poiseuille flow in a rectangular duct. In addition, convection through the porous medium may be found in the fiber and granular insulation, including structures for high power density, electric machines and nuclear reactors.

Viscosity is the measure of a fluid’s resistance to the flow and it describes the internal friction of a moving fluid. Due to the added resistance of the porous structure, the effective viscosity of the porous medium is generally found different from the fluid viscosity and it depends on the properties of porous media (e.g., porosity and permeability). The effective viscosity of the porous matrix is generally larger than the fluid viscosity, however, the available literatures, unfortunately, exhibit a remarkable disagreement with respect to the viscosity ratio. Brinkman [1] has taken both viscosities as same. Lundgren [2] has shown that, for the flow phenomenon near a random bed of spheres, the effective viscosity is larger than the fluid viscosity for high porosity and smaller for low porosity. Further, the viscosity ratio has been taken as greater than one in the studies of Adler and Mills [3] and Kim and Russell [4] while in Koplik, J., Levine, H., Zee, A. [5] , Kolodziej [6] , Nield and Bejan [7] and Breugem [8] , it has been taken as less than one.

In addition to the analysis based on the basic conservation laws, the analysis of second-law of thermodynamics is important in understanding the entropy generation, which is attributed to the thermodynamic irreversibility. This kind of thermodynamic analysis is significant for studying the optimum operating conditions, which helps in designing a system with less entropy and destruction of available work (energy). According to the Gouy-Stodola theorem, the lost available work is directly proportional to the entropy generation. The utilization of the second law of thermodynamics in convective heat transfer is very well presented by Bejan [9] [10] . By comparing the results of various cross sectional duct geometries, Sahin [11] has shown that the circular duct geometry is favourable one when the frictional contribution of the entropy generation becomes dominant. Narusawa [12] has presented the theoretical and numerical studies to the cross-sectional rate of the entropy generation for a fully developed flow, consisting of contributions from heat transfer and viscous dissipation. By applying the first and second laws of thermodynamics to forced convection in a cylindrical annular space with isoflux boundary conditions, Mahmud and Fraser [13] have obtained the general expressions for the velocity and temperature distributions, entropy generation number and Bejan number. Mohammad [14] has demonstrated the flow field and heat transfer with laminar forced convection in conduits field with a porous material to different degrees. He calculated that the effect of the Darcy number on the heat transfer in a fully developed flow region may largely be neglected for in the case of homogeneously filled channel.

Sauoli and Sauoli [15] have obtained the effect of the second law on laminar falling liquid film along an inclined heated plate by considering the upper surface of the liquid film free and adiabatic and the lower wall having constant heat flux. Their results show that the entropy generation increases for all values of the group parameters. Makinde and Osalusi [16] have studied the impact of the second law in the case of laminar flow in a channel field with saturated porous media. Yadav and Singh [17] have shown that the Darcy dissipation dominates the viscous dissipation everywhere of the channel except near the walls in a fully developed horizontal porous channel and further Yadav and Singh [18] have presented in details the analysis of entropy generation in an annular porous duct.

Heat transfer from solid walls to flowing fluids is an area of extreme scientific interest as well as of immense practical importance. Certain flows, passing through the bodies with high porosity, do not follow the Darcy’s law and Brinkman’s model is applicable for that type of flows. Neale and Nader [19] have also given the practical significance of the Brinkman extension of Darcy law. Accurate description of fluid flow behavior in the porous media is essential for successful design and operation of systems in this area. In the present paper, our main focus is to see the effect of viscous and Darcy dissipations in the energy equation of rectangular porous duct by employing the Brinkman model to calculate the non-dimensional entropy generation.

2. Mathematical Formulation

Consider the steady, laminar, two dimensional incompressible fluid flow in a saturated porous medium bounded by a rectangular duct of width, 2W and height, H with the origin of a coordinate system located at the corner of the rectangular duct as sketched in Figure 1, and the direction of gravity, g, in the negative y direction. The temperature at the bottom and the top are kept at constant temperatures and, respectively with small temperature difference; while the side boundaries are taken as adiabatic. For steady, low-Rayleigh number mixed convection, pressure variation may be decomposed into two parts, a variation in the

Figure 1. Schematic Diagram.

main flow (z) direction and a cross sectional variation in the x-y plane; i.e. with


We also assume that the fluid motion can adequately be described by Boussinesq approximation. Under the assumptions of constant thermo-physical properties and linear Boussinesq approximations, the governing conservation equations, namely the equation of continuity, momentum and thermal energy for the isotropic and homogeneous porous medium may be written as:






For the fully developed flow, the above equations may be observed that, these are identical to the governing equations of the Rayleigh-Benard convection in the porous medium (i.e. natural convection in rectangular porous cavity). The entropy is generated due to the presence of irreversibility, and entropy generation is adopted as a quantative measure of the irreversibility associated with a thermal process. Therefore, under supercritical conditions the longitudinal convective rolls are superposed to the forced convective flow. The local volumetric rate of entropy generation, based on second law of thermodynamics, is given as


The first term on the right hand side of above equation is due to the heat transfer while the second and third terms are due to the viscous and Darcy dissipations. It is well known that the viscous and Darcy dissipations are negligible in the first law of thermodynamics in many applications in clear fluid flow but it cannot be neglected in porous medium. The ratio, becomes, where is the ratio of viscous dissipation to conduction in the first law with and indicating the characteristics scales, which, in the present reference, mention for an average temperature, , and the temperature difference, , respectively. Hence, even if the dissipations are not negligible in the first law (which is the case of present study); it may be important in the second law due to the presence of the factor. Our objective is to evaluate the cross sectional entropy generation, , for fully developed mixed convection by integrating over a cross section; i.e.



The subscript, , stands for the fully developed mixed convection. Since the temperature variations are small compared to the absolute temperature T, the factors, and, in Equation (6) are assumed to be constant in performing the integrations.

For the rate of entropy generation due to the heat transfer, we obtain, after applying integration by parts as well as the adiabatic side wall conditions,


The second term on the right hand side is zero as we apply conservation of energy, Equation (5), the hydrodynamic boundary conditions and integration by parts; i.e.


Also by integrating conservation of energy, Equation (5), over a cross section it may be shown that heat trans- fer at the top boundary is balanced by heat transfer at the bottom boundary; i.e.




with In terms of a Rayleigh number, a Nusselt number and the mixing cup temperature, , Equation (11) may be expressed as




with. Given and for a specified fluid the entropy generation rate due to heat transfer is proportional to the Nusselt number.

The rate of entropy generation due to viscous and Darcy dissipations, upon application of the hydrodynamic conditions as well as continuity, becomes


where. The first term on the right hand side of Equation (14), we may substitute conservation of momentum from Equation (2) to Equation (4), then repeated applications of integration by parts and the hydrodynamic boundary conditions yield,






The first term on the right hand side is due to the forced convection (buoyancy induced flow). To represent the first term in Equation (17) in terms of non-dimensional parameters relevant to forced convection, we introduce


where are mean z-direction velocity and hydraulic diameter respectively, yielding


The rate of entropy generation due to the viscous dissipation associated with the buoyancy-induced flow in porous medium, on the other hand, may be expressed as


where with



The rate of non-dimensional entropy generation over a cross section for fully developed mixed convection in porous medium, is


The right hand side of Equation (23) consists of contributions from heat transfer, cross-sectional (buoyancy-induced) flow and longitudinal main flow.

Now if Equation (5) includes the viscous as well as Darcy dissipations, i.e.


then Equation. (23) becomes


3. Result and Discussions

Entropy is a thermo-dynamical property that is a measure of the energy unavailable for useful work in a thermodynamic process, such as in energy conversion devices, engines, or machines, which can only be driven by convertible energy. When a substance is heated or cooled, there is a change in the entropy and has a theoretical entropy minimization (maximum efficiency) while converting the energy to useful work. From Equation (23) and Equation (25), we see that, if the thermal energy equation includes the viscous as well as Darcy dissipations then the second term (due to viscous dissipation) on the right hand side of the entropy generation rate is double and there is no change in first term (due to heat transfer).

4. Conclusion

This paper presents the analytical calculation for the non-dimensional entropy generation and the obtained result shows that there is no need to include the viscous and Darcy dissipations in the energy equation.


The author (SLY) would like to thank the Council of Scientific and Industrial Research, New Delhi, India, for financial support in the form of a Junior Research Fellowship.

Cite this paper

Shyam Lal Yadav,Ashok Kumar Singh, (2016) Analysis of Entropy Generation in a Rectangular Porous Duct. Journal of Applied Mathematics and Physics,04,1336-1343. doi: 10.4236/jamp.2016.47143


  1. 1. Brinkman, H.C. (1947) On the Permeability of Media Consisting of Closely Packed Porous Particles. Applied Scientific Research, A1, 81-86.

  2. 2. Lundgren, T.S. (1972) Slow Flow through Stationary Random Beds and Suspensions of Spheres. Journal of Fluid Mechanics, 51, 273-299.

  3. 3. Adler, P.M. and Mills, P.M. (1979) Motion and Rupture of a Porous Sphere in a Linear Flow Field. Journal of Rheology, 23, 25-37.

  4. 4. Kim, S. and Russell, W.B. (1985) Modelling in Porous Media by Renormalization of Stokes Equations. Journal of Fluid Mechanics, 154, 269-286.

  5. 5. Koplik, J., Levine, H. and Zee, A. (1983) Viscosity Renormalization in the Brinkman Equation. Physics of Fluids, 26, 2864-2870.

  6. 6. Kolodziej, J.A. (1988) Influence of the Porosity of a Porous Medium on the Effective Viscosity in Brinkman’s Filtration Equation. Acta Mechanica, 75, 241-254.

  7. 7. Nield, D.A. and Bejan, A. (2013) Convection in Porous Media. Springer, New York.

  8. 8. Breugem, W.P. (2007) The Effective Viscosity of Channel-Type Porous Medium. Physics of Fluids, 19, 1-16.

  9. 9. Bejan, A. (2013) Convection Heat Transfer. 4th Edition, Wiley, New York.

  10. 10. Bejan, A. (1982) Second Law Analysis in Heat Transfer and Thermal Design. Advances in Heat Transfer, 15, 1-58.

  11. 11. Sahin, A.Z. (1998) A Second Law Comparison for Optimum Shape of Duct Subjected to Constant Wall Temperature and Laminar Flow. Heat and Mass Transfer, 33, 425-430.

  12. 12. Narusawa, U. (2001) The Second Law Analysis of Mixed Convection in Rectangular Ducts. Heat and Mass Transfer, 37, 197-203.

  13. 13. Mahmud, S. and Fraser, R.A. (2003) The Second Law Analysis in Fundamental Convective Heat Transfer Problems. International Journal of Thermal Sciences, 42, 177-186.

  14. 14. Mohammad, A.A. (2003) Heat Transfer Enhancement in Heat Exchanger Fitted with Porous Media. Part, I.: Constant Wall Temperature. International Journal of Thermal Sciences, 42, 385-395.

  15. 15. Sauoli, S. and Aiboud-Sauoli, S. (2004) Second Law Analysis of Laminar Falling Liquid Film along an Inclined Heated Plate. International Communication in Heat and Mass Transfer, 31, 879-886.

  16. 16. Makinde, O.D. and Osalusi, E. (2005) Second Law Analysis of Laminar Flow in a Channel Field with Saturated Porous Media. Entropy, 7, 148-160.

  17. 17. Yadav, S.L. and Singh, A.K. (2014) Effects of Viscous and Darcy Dissipations on Entropy Generation Rate of Flow through a Horizontal Porous Channel. International Journal of Energy and Technology, 6, 1-7.

  18. 18. Yadav, S.L. and Singh, A.K. (2016) Analysis of Entropy Generation in Annular Porous Duct. Transport in Porous Media, 111, 425-440.

  19. 19. Neale, G. and Nader, W. (1974) Practical Significance of Brinkman’s Extension of Darcy’s Law-Coupled Parallel Flows within a Channel and a Bounding Porous-Medium. The Canadian Journal of Chemical Engineering, 52, 475-478.

List of Symbol

aspect ratio,

specific heat at constant pressure

hydraulic diameter

Eckert number,

friction factor,

gravitational acceleration

height of the duct

thermal conductivity

average Nusselt number at the bottom based on H, with

being the temperature gradient at the bottom averaged over the width


Peclet number,

Prandtl number

Rayleigh number,

cross sectional entropy generation

local volumetric rate of entropy generation

mixing-cup temperature

temperature of the bottom boundary

temperature of the top boundary

velocity in x (lateral) direction

velocity in y (vertical) direction

. velocity in z (main flow) direction

mean velocity in z-direction

half width of a duct

transverse coordinate

vertical coordinate

longitudinal coordinate

Greek Symbols

non-dimensional temperature

non-dimensional mixing-cup temperature

absolute viscosity

viscosity in presence of porous media

kinematic viscosity




forced convection

due to heat transfer

mixed convection

due to viscous and Darcy dissipations