Entropy generation and pumping power to heat transfer ratio (PPR) of a laminar flow, for a circular tube immersed in an isothermal fluid, are studied analytically in this paper. Two different fluids, namely, water and ethylene glycol, are chosen to study the influence of fluid properties on entropy generation and PPR. The expressions for dimensionless entropy generation, Bejan number and PPR are derived in a detailed way and their variations with Reynolds number, external Biot number, and the dimensionless temperature difference are illustrated. The results of the analysis are compared with those for a laminar flow in a circular tube with uniform wall temperature boundary condition. Finally, a criterion is established to determine which type of thermal boundary conditions is more suitable for a particular fluid, with respect to its influence on entropy generation.
1. Introduction
Heat transfer is a fundamental source of thermodynamic irreversibility in all real engineering devices. When heat is transferred across a finite temperature difference, some capacity to do work is lost. In convection, apart from heat transfer, fluid friction is the other source of loss of available work. Both heat transfer and fluid friction generate entropy. This entropy generation must be minimized to reduce the loss of available work. Entropy generation minimization is no longer an avant-garde philosophy but a mainstream one in the design of thermal systems.
In the past, many researchers have studied the problem of entropy generation minimization in fluid flow with heat transfer. Bejan [1, 2] outlined the method for evaluating the entropy generation in fluid flow with heat transfer. He found out the entropy generated in fluid flow with heat transfer over a flat plate, in a duct, for cylinders in cross flow and in various other geometrical configurations. Şahin [3] studied the entropy generated in a circular duct with uniform wall temperature for two fluids, namely, water and glycerol. The effect of temperature on viscosity was taken into account in that study. He found out that the total energy loss (due to pumping process and entropy generation) can be minimized with respect to the duct length for viscous fluids with temperature dependent viscosity. Şahin [4, 5] also considered the effect of duct geometries on the entropy generation, both for uniform wall temperature and for uniform heat flux boundary conditions. The various duct geometries like circular, square, rectangular, equilateral triangle, and sinusoidal were considered. In general, circular geometry was found to be preferable, especially for high Reynolds number. A similar study for hexagonal cross-section was carried out by Jarungthammachote [6] and for a helical coil by Ko and Ting [7]. Dağtekin et al. [8] analyzed the entropy generation in a circular tube with internal fins of different shapes. The study of entropy generation in a pipe was extended to non-Newtonian fluids by Mahmud and Fraser [9, 10] and by Hung [11].
As far as the authors are aware, the studies carried out so far have been confined to tubes/ducts with uniform wall temperature or with uniform wall heat flux. The present study is concerned with entropy generation in a circular duct immersed in an isothermal external fluid, which is a more general thermal boundary condition. According to Sparrow and Patankar [12], both the uniform heat flux and uniform wall temperature boundary conditions are the limiting cases of this boundary condition. To study the effect of fluid properties, two different fluids, water and ethylene glycol, have been considered in this paper. The effects of dimensionless temperature difference and external Biot number on the dimensionless entropy generation, Bejan number, and PPR have been investigated. A comparison has been made with uniform wall temperature boundary condition.
2. Analysis
Consider a circular tube immersed in an isothermal fluid (at temperature Te) as shown in Figure 1. Another fluid flows inside the tube. The flow inside the tube is assumed to be laminar and enters the tube at a uniform temperature Ti. The energy balance for a control volume (CV) of infinitesimal length dx (neglecting axial conduction, changes in kinetic and potential energies, and viscous dissipation) is given by
(1)m˙CpdT=q˙′dx,
where q˙′, the rate of heat transfer per unit length from the external fluid to the internal fluid, is given by
(2)q˙′=h^2πr(Te-T)
(h^ in (2) is the overall heat transfer coefficient between the external fluid and the internal fluid).
An infinitesimal control volume in a circular duct immersed in an isothermal fluid.
Solving the two equations above for T(x), we get
(3)T(x)=Te-(Te-Ti)e(-4h^x/ρUDCp).
Using nondimensional parameters, we can write this equation as
(4)T(x)=Te-(Te-Ti)e(-4Stλ),where St is the Stanton number given by
(5)St=h^ρUCpand λ is the non-dimensional duct length, given by
(6)λ=xD.
The temperature at the outlet of the duct, Tout, is given by
(7)Tout=T(L)=Te-(Te-Ti)e(-4StλL),
where
(8)λL=LD.
Writing (7) in terms of mass flow rate, we get
(9)Tout=T(L)=Te-(Te-Ti)e(-N^ukπL/m˙Cp).
The rate of total heat transferred to the duct is given by
(10)Q˙=m˙Cp(Tout-Ti)=m˙Cp(Te-Ti-(Te-Ti)e-4StλL).
The entropy generated in this infinitesimal CV is given by
(11)dS˙gen=-q˙′dxTwall+m˙ds.
Now,
(12)q˙′=h^(Te-T)πD=h(Twall-T)πD⟹Twall=h^(Te-T)h+T.
Substituting the values of q˙′ and Twall from (12) into (11), we get the following:
(13)dS˙gen=-h^(Te-T)πDdxh^(Te-T)/h+T+m˙ds.
For an incompressible fluid,
(14)m˙ds=m˙CpdTT-m˙dPρT,
where the pressure drop per unit length is related to the friction factor as
(15)-dPdx=fρU22D.
Thus we get
(16)m˙ds=m˙CpdTT+m˙fU22DTdx.
Substituting (16) to (13) we get
(17)dS˙gen=-h^(Te-T)πDdxh^(Te-T)/h+T+m˙CpdTT+m˙fU22DTdx.
Substituting the value of T(x) from (4) into (17) and integrating (17) from x=0 to x=L, we get the total entropy generation rate in the duct as
(18)S˙gen=m˙Cp[-1αln|1-ταe-4StλL1-τα|+ln|1-τe-4StλL1-τ|+fEcτ8Stln|e4StλL-τ1-τ|],
where St is the Stanton number (defined in (5)),
(19)α=1-N^uNu,(20)τ=Te-TiTe.
Ec is the Eckert number given as
(21)Ec=U2Cp(Te-Ti).
In the literature, dimensionless entropy generation has also been used to quantify the entropy generated. It is given by
(22)ψ=S˙genm˙Cp=-1αln|1-ταe-4StλL1-τα|+ln|1-τe-4StλL1-τ|+fEcτ8Stln|e4StλL-τ1-τ|.
If we take the limit α→0, that is, (N^u→Nu), (18) reduces to
(23)limα→0S˙gen=m˙Cp{limα→0(-1αln|1-ταe-4StλL1-τα|)limα→0S˙gen=+limα→0(ln|1-τe-4StλL1-τ|)limα→0S˙gen=+limα→0(fEcτ8Stln|e4StλL-τ1-τ|)}.
Using the identity lima→0(ln(1+a)/a)=1, we get in the simplified form the following:
(24)limα→0S˙gen=m˙Cp{-τ(1-e-4StλL)+ln|1-τe-4StλL1-τ|limα→0S˙gen=m˙Cp+fτEc8Stln|e4StλL-τ1-τ|},
which is the same as the expression for entropy generation for fully developed laminar flow in a circular pipe with uniform wall temperature as given in [3]. This corroborates the conclusion of Sparrow and Patankar [12] that, as the external Nusselt number tends to duct side Nusselt number, the isothermal external fluid boundary condition tends to the uniform wall temperature boundary condition.
Equation (22) can also be written as function of Re instead of St (St=N^u/RePr):
(25)ψ=-1αln|1-ταe-π1/Re1-τα|+ln|1-τe-π1/Re1-τ|+π2Re2ln|eπ1/Re-τ1-τ|,
where
(26)π1=4N^uλLkμCp,(27)π2=(fRe)μ38N^uTe(Dρ)2k.
2.1. Bejan Number
The dimensionless entropy generation gives a good idea of the rate of total entropy generation. However, it does not convey which out of the two entropy generation mechanisms, namely, heat transfer and fluid friction, dominates. To resolve this, Paoletti et al. [13], defined a dimensionless number, called Bejan number, as
(28)Be=EntropyGeneratedDueToHeat TransferTotalEntropyGenerated.
At Be=1, all the entropies are generated due to heat transfer and at Be=0, all the entropies generated are due to fluid friction.
In the current analysis, if we take limTe→Tiψ, it gives us the dimensionless entropy generated due to fluid friction only (as Te→Ti, there is no heat transfer, and so entropy generated will be due to fluid friction only):
(29)FF=limTe→Tiψ=π1×π2×Re,
where FF is dimensionless entropy generated due to fluid friction only.
So,
(30)Be=ψ-FFψ=1-FFψ.
2.2. Pumping Power to Heat Transfer Ratio
Another parameter which is a good indicator of the performance of a heat exchanger pipe is the ratio of pumping power to heat transfer (PPR), which can be expressed as
(31)PPR=AΔPUQ˙.
Here the pressure drop ΔP along the length of the duct can be obtained by integrating the infinitesimal pressure drop dP given by
(32)dP=-fρU22Ddx.
Using (32) to calculate pressure drop and (10) to calculate heat transfer rate and inserting the same into (31), we get
(33)PPR=μ2λL(fRe)(Re/τ)2ρ2D2CpTe[1-e(-π1/Re)].
2.3. Fluid Properties
To study the effect of fluid properties, two different fluids, water and ethylene glycol, are considered.
Water. The thermophysical properties of water are calculated using the following equations, given in [6]:
(34)ρw(kg/m3)=1000(1-(T+288.9414)(T-3.9863)2508929.2(T+68.12963)),(35)kw(W/mk)=-0.3838+0.005254×(T+273)kw(W/mk)=-0.000006369×(T+273)2,(36)μw(Ns/m2)=(exp[-2.10+(-4.45)(273(T+273))(273(T+273))2μw(Ns/m2)=+6.55(273(T+273))2])(0.00179),(37)Cpw(J/kgK)=(7.25575005×101-6.62445402×10-1(T+273)+2.56198746×10-33(T+27)2-4.3659123×10-6(T+273)3+2.78178981×10-93(T+27)4)×R0.018.
In the above equations, T is the temperature of water in degree C and R is the universal gas constant.
Ethylene Glycol. The density of ethylene glycol is taken as constant, since it varies roughly by only 1% in the temperature range considered. For determining the viscosity of ethylene glycol, the following empirical correlation given by Sherman [14] is used:
(38)μ(T)=μ(Tref)(TTref)nexp[B(1T-1Tref)].
In (38), Tref is the reference temperature whose value is taken as 273 K. B and n are fluid dependent parameters, whose numerical values for ethylene glycol are found out using regression. To estimate thermal conductivity and specific heat, the fourth degree polynomial equations are used. The coefficients of these polynomial equations are also determined using regression. The data for regression is taken from [15]. The regression is carried out using the commercially available software Microsoft Excel.
The fluid properties are evaluated at Tavg=(Tout+Ti)/2 using an iterative procedure. A value of Tout is guessed. Using this guessed value, the fluid properties are found out at Tavg. These fluid properties are then used to find out the new Tout using (7). This process is repeated till the difference in the successive values of Tout is less than 1 K. The Tout for the last iteration is taken as the final value. The Tavg, which is calculated using the final value of Tout, is used to estimate the fluid properties.
3. Results and Discussion
The Nu and N^u are functions of external Bi. Their values are given in tabular form in [12] and are presented in Table 1. The value of friction factor for laminar flow is taken from [16] and is represented as:
(39)f=64Re.
Nu and N^u as a function of Bi taken from [12].
Bi
Nu
N^u
0
4.364
0
0.1
4.330
0.1909
0.25
4.284
0.4471
0.5
4.221
0.8075
1
4.122
1.345
2
3.997
1.998
5
3.840
2.773
10
3.758
3.163
100
3.663
3.597
∞
3.657
3.657
The numerical values of the constant parameters are given in Table 2. For these parameters, the effect of Bi and τ on the dimensionless entropy generation ψ, the Bejan number Be, and the PPR has been studied. The results are detailed below.
Values of different constant parameters used in the analysis.
Parameters
Numerical values
Te
310 K
D
0.08 m
L
1.0 m
3.1. Effect on ψ
Figure 2 shows the variation of ψ with Re for different values of (a) Bi and (b) τ for water. As can be seen from the figure, ψ decreases with Re. To explain this, we consider (4), where we find that if we let Re tend to infinity, then the T(x) tends to Ti (since St=Nu/(Re·Pr)); that is,
(40)limRe→∞T(x)=Ti.
Variation of dimensionless entropy generation ψ of water with Re for different values of (a) Bi and (b) τ.
Therefore, at higher Re the temperature gradients inside the fluid decrease and there is less entropy generated due to heat transfer. Even though as Re increases, the entropy generated due to pressure increases, but as we will see later, the contribution of entropy generated due to pressure to total entropy generated is very less (Be is close to 1), so ψ decreases when Re increases.
It is also seen from Figure 2(a) that as Re increases, the curves for different Bi tend to merge. Similarly, the curves for different τ also tend to merge at higher Re (see Figure 2(b)). This is because, as Re increases, the contribution of heat transfer to entropy generation decreases and total entropy generated is roughly due to pressure only. Since entropy generation due to pressure is independent of Bi or τ, the curves for different Bi (and for different τ) tend to merge at higher Re.
We see from Figure 2(a) that ψ increases with increase in Bi. This is because as Bi increases, N^u increases (see Table 1), which means that heat transfer rate to the fluid increases. This leads to increase in temperature gradients inside the fluid. This increase in temperature gradient inside the fluid leads to increase in entropy generation.
Figure 2(b) shows that, as the dimensionless temperature difference τ increases, ψ increases. This is because an increase in τ causes the heat transfer to the fluid to increase, which leads to an increase in temperature gradients and hence in ψ.
Figure 3 shows the variation of ψ with Re for different values of (a) Bi and (b) τ for ethylene glycol. Similar trends, as observed for water, can be observed here.
Variation of dimensionless entropy generation ψ of ethylene glycol with Re for different values of (a) Bi (b) τ.
3.2. Effect on Be
Figure 4 shows the variation of Be with Re for different values of (a) Bi and (b) τ for water. As can be seen from the figure, the value of Bejan number is very high (close to 1). This means that heat transfer irreversibility contributes almost 100% to the total entropy generation. This is corroborated by Ben-Mansour and Sahin [17] for circular pipes and by Jarungthammachote [6] for hexagonal pipes. We can also see from the figures that Be decreases with increase in Re. This is because, as pointed out in the previous section, the temperature gradients inside the fluid decrease with increase in Re. This leads to a decrease in heat transfer irreversibility. Moreover, the irreversibility due to pressure difference increases with increase in Re. Due to these two reasons, Be decreases with increase in Re.
Variation of Bejan number of water with Re for different values of (a) Bi and (b) τ.
We can also see from Figure 4(a) that, as Bi increases, Be increases for a fixed Re. From Figure 4(b), we see that, as τ increases, Be increases. The result of these two figures can be explained on the basis of heat transfer irreversibility. As Bi or τ increases, the heat transfer to the fluid increases, due to which temperature gradients inside the fluid increase. This leads to increase in heat transfer irreversibility, which further leads to increase in Be.
Figure 5 shows the variation of Be with Re for different values of (a) Bi and (b) τ for ethylene glycol. Trends, similar to the case of water, can also be observed for ethylene glycol.
Variation of Bejan number of ethylene glycol with Re for different values of (a) Bi and (b) τ.
3.3. Effect on PPR
Figure 6 shows the variation of PPR with Re for different values of (a) Bi and (b) τ for water. As can be seen from the figure, for very small values of Re, PPR tends to 0. This is because, as Re→0, there is no flow and so no pumping power is required. As Re increases, PPR increases, because, to maintain a flow at higher Re, more pumping power is required. Figure 6(a) shows that, as Bi increases, the PPR decreases. This is because, as Bi increases, N^u increases (see Table 1), which means more heat transfer to the fluid for the same Re. This leads to decrease in PPR. For the same reason, PPR decreases with increase in τ, as shown in Figure 6(b).
Variation of PPR of water with Re for different values of (a) Bi and (b) τ.
Figure 7 shows the variation of PPR with Re for different values of (a) Bi and (b) τ for ethylene glycol. Trends, similar to the case of water, can be observed.
Variation of PPR of ethylene glycol with Re for different values of (a) Bi and (b) τ.
3.4. Comparison between Two Fluids
The comparison of total entropy generation, entropy generation due to heat transfer, and entropy generation due to pressure between water and ethylene glycol for the same mass flow rate, same Biot number, and same dimensionless temperature difference is shown in Table 3.
Comparison of entropy generation due to temperature difference and entropy generation due to pressure for water and ethylene glycol.
m˙(kg/s)
Water
Ethylene glycol
S˙gen,Total×103(W/K)
S˙gen,ΔT×103 (W/K)
S˙gen,ΔP×106 (W/K)
S˙gen,Total×103 (W/K)
S˙gen,ΔT×103 (W/K)
S˙gen,ΔP×106 (W/K)
0.05
0.839260096
0.839254054
0.006042
0.159089483
0.159012754
0.076729
0.06
0.847344381
0.847335678
0.008703
0.159603359
0.159492849
0.110510
0.07
0.853184312
0.853172464
0.011848
0.159987558
0.159837123
0.196506
0.08
0.857601029
0.857585551
0.015478
0.160292650
0.160096145
0.248721
0.09
0.861058756
0.861039165
0.019591
0.160546880
0.160298159
0.307081
As can be seen from Table 3, the total entropy generated and entropy generated due to heat transfer are higher for water while the entropy generated due to pressure is higher for ethylene glycol.
To investigate the reason for this trend, (9) must be paid attention to. For the same mass flow rate, the temperature at the outlet depends on k/Cp. Even though Cp of water is higher than that of ethylene glycol, the k for water is still higher, with the result that k/Cp is higher for water. So Tout for water is higher than Tout for ethylene glycol. Due to this, temperature gradients inside the water are higher which leads to higher entropy generation for water than for ethylene glycol.
Entropy generation due to pressure is higher for ethylene glycol than that of water. This is because the viscosity of ethylene glycol is much higher than that of water.
3.5. Comparison with Uniform Wall Temperature
Figure 8 shows the comparison of (a) ψ and (b) PPR between the circular duct with uniform wall temperature and circular duct immersed in an isothermal fluid. For comparison purpose, the temperature of isothermal external fluid is kept same as uniform wall temperature. As can be seen from the figure, ψ is higher for uniform wall temperature condition than for isothermal external fluid condition, while PPR is lower for uniform wall temperature condition than for isothermal external fluid condition; other parameters are kept the same.
Comparison of (a) ψ and (b) PPR for different boundary conditions.
While comparing ψ and PPR for different types of boundary conditions, the question arises: which type of boundary conditions is more suitable? To answer this, it should be noted that the total exergy loss in the pipe is the sum of ψ and PPR. As can be seen from Table 4, for less viscous fluids like water, ψ is much higher than PPR, a trend which becomes more pronounced for smaller Re. But more viscous fluids like ethylene glycol, on the other hand, have higher PPR than ψ. From this, it can be concluded that, for fluids whose viscosity is low, isothermal external fluid boundary condition should be preferred, while, for more viscous fluids, like ethylene glycol, uniform wall temperature boundary condition is to be preferred.
Comparison of dimensionless entropy generation ψ with PPR for water and ethylene glycol.
Re
Water
Ethylene glycol
ψ
PPR
ψ
PPR
100
2.62E-05
5.63E-10
8.78E-07
7.76E-06
500
8.08E-06
1.15E-08
1.82E-07
0.000193
1000
4.28E-06
4.48E-08
9.8E-08
0.000774
1500
2.91E-06
9.99E-08
7.33E-08
0.00174
2000
2.21E-06
1.77E-07
6.33E-08
0.003094
2500
1.78E-06
2.76E-07
5.91E-08
0.004834
3000
1.49E-06
3.96E-07
5.8E-08
0.00696
4. Conclusions
Dimensionless entropy generation, Be, and PPR of laminar flow in a circular tube immersed in an isothermal fluid are studied in this paper. The conclusions obtained can be summarised as follows.
With increase in Re, both the dimensionless entropy generation and Be decrease, while PPR increases.
With increase in Bi, the dimensionless entropy generation and Be increase, while PPR decreases.
With increase in τ, the dimensionless entropy generation and Be increase, while PPR decreases.
The entropy generated due to heat transfer is more for water than for ethylene glycol, while the entropy generated due to pressure is more for ethylene glycol than for water, for the same mass flow rate.
PPR for ethylene glycol is higher than that for water, since ethylene glycol is more viscous than water.
For more viscous fluids like ethylene glycol, uniform wall temperature boundary condition should be preferred, since it leads to less PPR. For less viscous fluid like water, a tube immersed in an isothermal fluid boundary condition should be preferred, since it leads to less entropy generation.
The values of N^u and f used to derive the expression for dimensionless entropy generation have been taken from [12] and [16], respectively. In [12], the value of N^u has been derived numerically, while in [16] the value of f has been derived analytically and verified experimentally. Thus, the margin of probable error in the values of N^u and f is very less. In the highly unlikely case of small errors creeping in the values of N^u and f, the trends of dimensionless entropy generation presented here will remain unchanged qualitatively.
NomenclatureA:
Cross-sectional area (m2)
Be:
Bejan number
Bi:
Biot number her/k
Cp:
Specific heat capacity (J/kg K)
CV:
Control volume
D:
Cross-sectional diameter (m)
Ec:
Eckert number (defined in (21))
f:
Friction factor
FF:
Dimensionless entropy generated due to fluid friction
h:
Internal heat transfer coefficient q˙′′/(Twall-T) (W/m2 K)
h^:
Overall heat transfer coefficient q˙′′/(Te-T) (W/m2 K)
he:
External heat transfer coefficient q˙′′/(Te-Twall) (W/m2 K)
k:
Thermal conductivity (W/m K)
L:
Duct length (m)
m˙:
Mass flow rate (kg/s)
Nu:
Internal Nusselt number hD/k
N^u:
Overall Nusselt number h^D/k
P:
Pressure (Pa)
PPR:
Pumping power to heat transfer ratio (defined in (31))
Pr:
Prandtl number
q˙′′:
Local heat flux (W/m2)
q˙′:
Heat rate per unit length (W/m)
Q˙:
Heat rate (W)
r:
Cross-section radius (m)
Re:
Reynolds number
s:
Specific entropy (J/kg·K)
S˙:
Entropy rate (W/K)
St:
Stanton number (defined in (5))
T:
Temperature (K)
U:
Average fluid velocity (m/s).
Greek Lettersα:
Dimensionless parameter (defined in (19))
τ:
Dimensionless temperature difference (defined in (20))
λ:
Dimensionless length (defined in (6))
λL:
Dimensionless length (defined in (8))
ψ:
Dimensionless entropy generation (defined in (22))
ΔP:
Pressure drop (Pa)
ρ:
Fluid density (kg/m3)
μ:
Fluid viscosity (Ns/m2)
π1:
Dimensionless parameter (defined in (26))
π2:
Dimensionless parameter (defined in (27)).
Subscriptsavg:
Average
e:
External
gen:
Generated
i:
Inlet
out:
Outlet
ref:
Reference
Total:
Total
w:
Water
wall:
Wall
ΔP:
Pressure difference
ΔT:
Temperature difference.
Acknowledgment
Sincere thanks are due to Dr. T. K. Vashist of Infotech Enterprises Limited, Bangalore, for his contribution and guidance towards revising the paper.
BejanA.Second law analysis in heat transfer198058-97217322-s2.0-0018506296BejanA.1982New York, NY, USAJohn Wiley & SonsSahinA. Z.Second law analysis of laminar viscous flow through a duct subjected to constant wall temperature19981207683ŞahinA. Z.Irreversibilities in various duct geometries with constant wall heat flux and laminar flow19982364654732-s2.0-003210425410.1016/S0360-5442(98)00010-3ŞahinA. Z.A second law comparison for optimum shape of duct subjected to constant wall temperature and laminar flow1998335-64254302-s2.0-0032045250JarungthammachoteS.Entropy generation analysis for fully developed laminar convection in hexagonal duct subjected to constant heat flux20103512537453792-s2.0-7864987462510.1016/j.energy.2010.07.020KoT. H.TingK.Entropy generation and thermodynamic optimization of fully developed laminar convection in a helical coil2005321-22142232-s2.0-1464439289410.1016/j.icheatmasstransfer.2004.04.039DaǧtekinI.ÖztopH. F.ŞahinA. Z.An analysis of entropy generation through a circular duct with different shaped longitudinal fins for laminar flow20054811711812-s2.0-884426383610.1016/j.ijheatmasstransfer.2004.08.013MahmudS.FraserR. A.Inherent irreversibility of channel and pipe flows for non-Newtonian fluids20022955775872-s2.0-003666007110.1016/S0735-1933(02)00377-9MahmudS.FraserR. A.Second law analysis of forced convection in a circular duct for non-Newtonian fluids20063112222622442-s2.0-3374584675710.1016/j.energy.2005.09.003HungY.-M.Viscous dissipation effect on entropy generation for non-Newtonian fluids in microchannels2008359112511292-s2.0-5084910424610.1016/j.icheatmasstransfer.2008.06.005SparrowE. M.PatankarS. V.Relationships amoung boundary conditions and nusselt numbers for thermally developed duct flows19779934834852-s2.0-0017526421PaolettiS.RispoliF.SciubbaE.Calculation of exergetic losses in compact heat exchanger passages1989ASME AES- 10ShermanF. S.1990New York, NY, USAMcGraw HillIncroperaF. P.DeWittD. P.20065thSingaporeJohn Wiley & SonsFoxR. W.McDonaldA. T.20045thSingaporeJohn Wiley & SonsBen-MansourR.SahinA. Z.Entropy generation in developing laminar fluid flow through a circular pipe with variable properties20054211112-s2.0-2684445985310.1007/s00231-005-0637-6