## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Brazilian Journal of Chemical Engineering

##
*Print version* ISSN 0104-6632

### Braz. J. Chem. Eng. vol.29 no.3 São Paulo July/Sept. 2012

#### http://dx.doi.org/10.1590/S0104-66322012000300008

**FLUID DYNAMICS; HEAT AND MASS TRANSFER; AND OTHER TOPICS**

**Second law analysis for free convection in non-newtonian fluids over a horizontal plate embedded in a porous medium: (prescribed heat flux)**

**W. A. Khan ^{I, *}; R. S. R. Gorla^{II}**

^{I}Department of Engineering Sciences, PN Engineering College, Phone: + 092-3218929467, National University of Sciences and Technology, PNS Jauhar, Karachi 75350, Pakistan. E-mail: wkhan_2000@yahoo.com

^{II}Department of Mechanical Engineering, Cleveland State University, Cleveland, Ohio, USA

**ABSTRACT**

Second law characteristics of heat transfer and fluid flow due to free convection of non-Newtonian fluids over a horizontal plate with prescribed surface heat flux in a porous medium are analyzed. Velocity and temperature fields are obtained numerically using an implicit finite difference method under the similarity assumption and these results are used to compute the entropy generation rate Ns, irreversibility ratio φ and the Bejan number Be for both Newtonian and non-Newtonian fluids. The effects of power-law index n, heat flux variation parameter λ , and modified duty parameter, G on the dimensionless entropy generation rate Ns, and the Bejan number Be are investigated and presented graphically.

**Keywords:** Second law; Free convection; Non-Newtonian fluids; Porous medium; Plate.

**INTRODUCTION**

The second law analysis of a non-Newtonian fluid over a horizontal surface has many significant applications in thermal engineering and industries. Applications of horizontal surfaces can also be found in various industrial exchanger systems. One of the fundamental problems of the engineering processes is the improvement in thermal systems during the convection in any fluid. The second law analysis is one of the best tools for improving the performance of the engineering processes. It investigates the irreversibility due to fluid flow and heat transfer in terms of the entropy generation rate. Bejan (1979, 1984, 1996) employed this method in many convection problems and later on several other investigators used this method successfully in different heat transfer problems.

Saouli and Saouli (2004) performed the second law analysis for Newtonian fluids along an inclined heated plate, whereas Saouli and Saouli (2009) and Gorla and Pratt (2007) performed the same analysis for non-Newtonian fluids with and without viscous dissipation effects. They considered the upper surface of the liquid film to be free and adiabatic and the lower wall fixed with constant heat flux. Their results show that the entropy generation number and the irreversibility ratio decrease in the transverse direction and increase as the viscous dissipation parameter increases. Recently, Saouli *et al*. (2006) investigated the entropy generation in a laminar, gravity-driven conducting liquid film along an inclined heated plate in the presence of a transverse magnetic field.

In this study we investigated the second law analysis of heat transfer and fluid flow due to free convection of non-Newtonian fluids over a horizontal plate with prescribed surface heat flux. The velocity and temperature distributions are determined by solving the Darcy and energy equations subject to an appropriate set of boundary conditions. The dimensionless entropy generation rate, irreversibility ratio and the Bejan number were computed using numerical values of the velocity, temperature and their gradients. The dimensionless parameters governing the problem are the dimensionless temperature variation parameter, fluid behavior index, the dimensionless transverse and axial distances, and Rayleigh number.

**ANALYSIS**

Consider a steady free convection of a power-law fluid along an impermeable horizontal plate (Fig. 1) with prescribed wall heat flux Q(x) = Ax^{λ} embedded in a porous medium. It is assumed that the physical properties of the power-law fluids, except the density, are constant and that the Boussinesq approximation for the case of maximum density is valid.

The ambient temperature is as shown in Fig. 1, where are the Cartesian coordinates along the horizontal surface and normal to it, respectively. The velocity components are along and normal to the wall surface. The acceleration due to gravity g acts vertically downwards. Assuming no heat friction and the thermal equilibrium between the fluid and the porous matrix, the governing equations in dimensional form, for this case, can be written as:

with boundary conditions:

where the bars denote dimensional quantities. Using the following non-dimensional variables

where q_{0} is the characteristic heat flux, Ra is the modified Rayleigh number, which is defined by:

we get the following governing equations in dimensionless form:

with dimensionless boundary conditions:

Using the following transformations

we get:

and the governing equations reduce to:

with reduced boundary conditions:

where primes denote the derivatives w. r. t.η Equations (13) and (14) were solved numerically for the velocity and temperature profiles and gradients using an implicit finite difference method, where a step size of Δη was used in all the calculations. These gradients are used in the calculations of dimensionless entropy generation rates. It is important to note that, for λ, Equations (13) and (14) reduce to the uniform wall flux case.

**Entropy Generation Rate**

Following Bejan [1], the local rate of entropy generation in convective heat transfer from a power-law fluid can be written as:

where

Using Eqs. (17) and (18) in Eq. (16), we get dimensionless local rate of entropy generation:

where:

with:

where T^{0}is the reference temperature, N_{C} is the entropy generation in the axial direction, N_{Y} is the entropy generation in the normal direction to the horizontal surface, N_{F} is the entropy generation due to fluid friction and G is the modified duty parameter.

The irreversibility ratio Φ is defined by Bejan (1984) as the ratio of the entropy generation due to the fluid friction N_{F} to the total entropy generation due to heat transfer ( N_{C}+N_{Y}), i.e.

It is important to note that when φ=1, both the heat transfer and fluid friction irreversibilities make the same contribution to generating entropy. But, when φ>1, fluid friction dominates and, when φ>1, heat transfer dominates. Due to importance of the contribution of entropy generation due to heat transfer, Paoletti *et al*. (1989) defined the Bejan number (Be) as the ratio of the entropy generation due to heat transfer to the total entropy generation, i.e.

where the Bejan number ranges from 0 to 1. When the heat transfer irreversibility dominates, Be = 1, when the irreversibility is dominated by fluid friction, Be = 0, and when the heat transfer and fluid friction entropy generation rates are equal Be = 1/2.

**RESULTS AND DISCUSSION**

Entropy generation rates due to heat transfer and fluid flow in free convection of Newtonian and non-Newtonian fluids over a horizontal plate with prescribed surface heat flux are analyzed graphically. Velocity and temperature fields are obtained by solving Equations (14) and (15) numerically with boundary conditions (16) using an implicit finite difference method under the similarity assumption. In order to validate the accuracy of the numerical results for dimensionless entropy generation rates and Bejan numbers, dimensionless heat transfer coefficients for free convection along a horizontal plate with prescribed wall heat flux are compared for a Newtonian fluid (n=1) with those of Gorla and Kumari (2003), and Aldoss *et al*. (1994). The differences between their results and the present results is less than 0.0098. Therefore, our results are highly accurate. This comparison is presented in Table 1.

The dimensionless entropy generation rates for a Newtonian fluid (n=1), a pseudoplastic fluid (n=0.5), and a dilatant fluid (n=1.5) over a horizontal plate in porous media are plotted in Figures 2-4 respectively. In each case, the effects of the viscous dissipation parameter G and the Rayleigh number Ra on the dimensionless entropy generation rate Ns are shown for different values of the temperature variation parameter λ. As expected, the entropy generation rate for Newtonian fluids is maximum at the plate surface due to maximum heat transfer and friction. However, there is no contribution to the entropy generation rate at the free surface, due to zero velocity and temperature gradients (Figure 2). For a Newtonian fluid, the entropy generation rates increase with the viscous dissipation parameter due to the increase in friction and they decrease with the increase in the flux variation parameter λ(Fig. 2(a). Figure 2(b) shows that the entropy generation rates increase with the increase in the Rayleigh number Ra. The variation of the dimensionless entropy generation rate with the transverse distance is η shown in Figure 3 for a pseudoplastic fluid (n = 0.5). Like Newtonian fluids, the entropy generation rates increase with the increase in viscous dissipation G and the Rayleigh number Ra. This is due to the decrease in heat transfer along the axial distance. Again, the dimensionless entropy generation rates are lower for higher values of the flux variation parameter λ. For dilatant fluids (n = 1.5), the entropy generation rates decrease initially with the transverse distance and then increase asymptotically for all parameters as shown in Figures 4(a) and 4(b), due to the increase in heat transfer and fluid friction irreversibilities. For this reason, dilatant fluids are not recommended for free convection along a horizontal plate under uniform surface temperature.

The variation of the Bejan number with the transverse distance for Newtonian, pseudoplastic and dilatant fluids is shown in Figures 5-7, respectively. They show the effects of the viscous dissipation parameter and Rayleigh numbers on the Bejan number. It is clear from Figure 5 that, for the Newtonian fluids, the Bejan number is a decreasing function of the dimensionless transverse distance and it decreases with the increase in the viscous dissipation parameter and the Rayleigh number, as shown in Figures 5(a) and 5(b), whereas the Bejan number increases with the increase in the temperature variation parameter in both cases. The Bejan number is higher at the plate surface due to higher heat transfer irreversibility. Unlike Newtonian fluids, the Bejan number for pseudoplastic fluids is an increasing function of the dimensionless transverse distance and it decreases with the increase in the viscous dissipation number and the Rayleigh number, as shown in Figures 6(a) and 6(b). The Bejan number increases with the increase in the flux variation parameter in both cases. This is due to the fact that the heat transfer irreversibility increases with the increase in the flux variation parameter. In the case of dilatant fluids, the heat transfer irreversibility is higher at the plate surface and remains constant initially, then fluid friction irreversibility dominates and, as a result, the Bejan number decreases with the transverse distance, as shown in Figures 7(a) and 7(b). Due to the nature and rheological properties of dilatant fluids, the heat transfer irreversibilities decrease with the increase in viscous dissipation and Rayleigh numbers and, as a result, the Bejan number decreases. However, with the increase in the flux variation and axial distance, heat transfer increases and, as a result, the Bejan number increases, as shown in Figures 7(a) and 7(b), respectively.

**CONCLUSIONS**

Second law analysis of heat transfer and fluid flow due to free convection of Newtonian and non-Newtonian fluids over a horizontal plate is performed for prescribed surface flux boundary conditions. Velocity and temperature fields are obtained numerically using an implicit finite difference method under the similarity assumption. These results are then used to determine the dimensionless entropy generation rates and Bejan numbers. It is observed that the entropy generation rate and Bejan number depend upon the viscous dissipation parameter G, the temperature variation parameter λ , the Rayleigh number Ra, the fluid behaviour index n, and the transverse and axial distances.

**NOMENCLATURE**

**REFERENCES**

Aldoss, T. K., Chen, T. S. and Armaly, B. F., Mixed convection over nonisothermal horizontal surfaces in a porous medium: The entire regime. Numerical Heat Transfer, 25, 685 (1994). [ Links ]

Bejan, A., A study of entropy generation in fundamental convective heat transfer. J. Heat Transfer, 101, 718 (1979). [ Links ]

Bejan, A., Convection Heat Transfer. Wiley, New York (1984). [ Links ]

Bejan, A., Entropy Generation Minimization. CRC Press, Boca Raton, New York (1996). [ Links ]

Gorla, R. S. R. and Kumari, M., Free convection in Non-Newtonian fluids along a horizontal plate in a porous medium. Heat and Mass Transfer, 39, 101 (2003). [ Links ]

Gorla, R. S. R. and Pratt, D. M., Second law analysis of a Non-Newtonian laminar falling liquid film along an inclined heated plate. Entropy, 9, 30 (2007). [ Links ]

Paoletti, S., Rispoli, F., Sciubba, E., Calculation of exergetic losses in compact heat exchanger passages. ASME AES, 10, 21 (1989). [ Links ]

Saouli, S., Saouli-Aboud, S, Second law analysis of laminar falling liquid film along an inclined heated plate. International Communications in Heat and Mass Transfer, 31, 879 (2004). [ Links ]

Saouli, A., Saouli, S., Settou, N., Meza, N., Thermodynamic analysis of gravity-driven liquid film along an inclined heated plate with hydromagnetic and viscous dissipation effects. Entropy, 8, 188 (2006). [ Links ]

Saouli, S. and Saouli, A., Second-law analysis of laminar Non-Newtonian gravity-driven liquid film along an inclined heated plate with viscous dissipation. Brazilian Journal of Chemical Engineering, 26, No. 2, 407 (2009). [ Links ]

Submitted: May 27, 2010

Revised: August 30, 2011

Accepted: February 26, 2012

* To whom correspondence should be addressed