Second Law Analysis of Laminar Flow in a Channel Filled with Saturated Porous Media: a Numerical Solution

This paper investigates entropy generation due to forced convection in a porous medium sandwiched between two parallel plates one of them being subjected to a uniform heat flux and the other one insulated. Our results showed that viscous dissipation will affect the entropy generation rate at the centerline of the channel since viscous dissipation is a quadratic function of velocity [1-3]. Neglecting the Darcy dissipation term in comparison with the terms added by Al-Hadrami et al. [4], will lead to the misunderstanding that fluid friction has no effect on the entropy generation rate at the tube centerline where the velocity derivative vanishes due to symmetry. Though the term added by [4] is O(Da) compared to the Darcy term one should not drop it unless the clear flow solution is sought [5-7]. Moreover, as stated by Nield [1], one should not use just the term involving velocity derivatives, as some authors have done in the past, for example [8-11]. Though in this paper the viscous dissipation effects in the energy equation are neglected, we have take them into account when it came to the entropy generation analysis.


Problem statement
Figure 1 shows the problem under consideration.When the flow in hydrodynamically fully developed one knows that the Brinkman flow model leads to the following momentum equation .0 with the porous medium shape factor s being defined as S=a/K 1/2 wherein K is the permeability and a is the channel width.The above will lead to two ordinary differential equations that one of them needs numerical integration to be solved while the other one leads to a simple exponential answer.

Numerical details
A numerical solution is presented for the momentum and the thermal energy equation, i.e.Eqns.(1-2) subjected to the aforementioned boundary conditions.To achieve this goal, we applied an implicit finite difference scheme, backward in x and central in y.Details of the method may be found in [12] and so we neglect repeating them.Length to height ratio of the channel is set at 10.A 0.05x0.005grid system is applied and it is observed that moving from this grid to a more intense one-being 0.001x0.0001-willchange the results within less than 1%.Changing the thermal boundary condition to a constant heat flux wall, we validated our work by defining a Nusselt number Nu=ha/(2k).It is observed that our predicted Nu differs from that of [3] in the fourth figure.One may recover the known analytical data for this case as reported by [3]; however, the details and exact figures are not presented in this report for the sake of brevity.
We present a first law of thermodynamics approach to verify our solution.Considering a slice of the channel with the dimensions a.dx and applying the first law of thermodynamics, one obtains wherein U and T m are mean velocity and bulk temperature.Rearranging the above equation in terms of non-dimensional parameters one finds Since a uniform temperature is applied at the channel inlet, the solution to Eq. ( 5) reads Our work resulted in Fig. 2.This figure may also be used for validation purpose.

Results and discussion
Dimensionless temperature distribution is illustrated in Figs. 3 through 5. Showing the results, xPe is called x for short.Figures 3 and 4 are designed to illustrate the dimensionless temperature profiles in the duct cross-section for some streamwise locations.The slopes at upper and bottom walls show that the former is adiabatic while the latter is at a constant heat flux in such a way that moving toward the upper wall both the transverse temperature gradient and the temperature itself decrease.Figure 5 shows the dimensionless temperature contours through the channel.Moving downstream, the isotherm values increase as fluid becomes heated near the bottom wall.Figures 6 through 10 illustrate the dimensionless entropy generation rate, Ns, for some values of s and Br while in all figures the values of Pe and Ω are fixed to 10 and 1, respectively.As a common trend in all figures one comprehends that the entropy generation rate decreases from the channel inlet to the outlet due to a fall in temperature gradients and this plunge is associated with a decrease in the heat transfer irreversibility term.As seen, for small values of s the plots of Ns experience a minimum while for large values of s one observes two local minima.For large values of s, the velocity changes in a thin region adjacent to the walls and this is accompanied by a steep temperature gradient near the bottom wall leading to a local maximum for Ns that happens near the bottom wall.Though near the upper wall the temperature gradients are small, Ns reaches a local maximum for a decrease in the temperature appearing in the denominator of Ns.Ns is defined here similar to Bejan [13], without the assumption that temperature difference in a cross-section is small, as Some authors have dropped the denominators to define an entropy generation number but this definition may work out to be a representative of entropy generation rate only if Ω 1 , i.e. a q kT in ′ ′ is negligible, as stated by Bejan [13] and applied by some authors [14][15][16][17].This can be mathematically interpreted as θ >> Ω ; which is not the case here.Far from the walls, in the channel centerline regions, the velocity is distributed somehow uniformly so that the velocity gradients are small in this region.On the other hand, the velocity magnitude contributes to fluid friction irreversibility (FFI).For this reason, unlike small s cases, the entropy generation plot will not pass through one minimum in the channel centerline.Actually, when s is large Ns recovers in the channel centerline for the velocity attains its maximum in this region and this fact is much more voiced for higher values of Br.Comparing Fig. 6 with Fig. 7, one realizes that moving from Br=0.1 to Br=1 with s=10 stretches the Ns plot by a factor, say, approximately 8. The two minima for large s will happen near the walls and for small s the so-called minimum takes place at the centerline.However, the maximum entropy generation rate depends on the interaction between heat transfer and fluid friction irreversibility that, in turn, depends on Pe, s, Br, and Ω .
energy equation in the absence of flow work, viscous dissipation, thermal dispersion and axial conduction effect is u ˆis the normalized velocity (local velocity divided by the average velocity, i.e. u ˆ=u/U)the effective thermal conductivity, ρ is the fluid density, c p is the specific at constant pressure, and q ′ ′ is the wall heat flux.The dimensionless coordinates are (x,y)=( y x, )/a and the thermal boundary conditions are 0

Figure 3 Figure 4 Figure 5
Figure 3 Dimensionless temperature distribution at different longitudinal locations (s=0.1)

Figure 6
Figure 6 Ns versus y for some longitudinal locations (s=10 and Br=1).

Figure 9 Figure 10
Figure 9  Ns versus y for some x with s=0.1 and Br=0.1 (legends are identified on Fig.6).