1 s2.0 S0038092X97000054 Main
1 s2.0 S0038092X97000054 Main
33-42, 1997
Pergamon 0 1997 Elsevier Science Ltd
PLI: SOO3%092X( 97)00005-4 All rights reserved. Printed in Great Britain
0038-092X/97 %17.00+0.00
Abstract-Artificial roughness has been found to enhance the heat transfer from the absorber plate to
the air in a solar air heater duct. However, this improvement is invariably accompanied by increased
pumping power. In this work, the effect of roughness and operating parameters on the thermal as well
as the hydraulic performance of roughened solar air heaters is discussed and the thermohydraulic perfor-
mance of roughened solar air heaters is compared with that of conventional smooth solar air heaters.
The optimum design and operating conditions have been determined. On the basis of thermohydraulic
considerations it has been found that the systems operating in a specified range of Reynolds number
show better thermohydraulic performance depending upon the insolation. A relationship between
the system and operating parameters that combine to yield optimum performance has been developed.
0 1997 Elsevier Science Ltd.
them are smaller relative roughness height of plant, and qth is the efficiency of thermal conver-
the elements with lower flow velocities sion of the power plant.
(Reynolds numbers less than 35,000) and a The rate of useful thermal energy may be
minimum insolation of 600 W m-‘. Further, obtained from the equation:
the range of Reynolds numbers investigated by
qll = F’[Z(za) - UL(I~- ti)/2]A* (2)
them was very high (Reynolds numbers from
about 15,000 to 60,000) and not normally appli- where
cable to solar air heaters which generally oper- F’=h/(h+ U,).
ate at Reynolds numbers below 15,000 so as to
allow a sufficient rise in temperature. An extens- For solar air heaters with absorber plates having
ive investigation of heat transfer and fluid flow artificial rib roughness, the heat transfer coeffi-
in solar air heater ducts has been conducted by cient, h, is given by the equation (Gupta, 1994):
one of the present authors (Gupta, 1994). It
was observed that non-transverse ribs are For e+ ~35:
advantageous when compared with transverse h =0.0024(e/D)0~oo’(W/H)-0~06(Re)‘~084
ribs for the enhancement of the heat transfer
coefficient. Correlations for the heat transfer exp[ - 0.04( 1 - a/60)2](k/D) (3)
coefficient and friction factor for non-transverse and for e’ 235:
ribs developed by the present authors have
been used in this work to evaluate the thermo- h=0.0071(e/D)-0~24(W/H)-0~028(Re)0~88
hydraulic performance of solar air heaters. exp[ -0.475( 1 -a/6O)‘](k/D). (4)
The rate of useful energy gain in a roughened
2. THERMOHYDRAULIC PERFORMANCE solar air heater may also be calculated using
OF SOLAR AIR HEATERS the equations:
Absorber Plate
m
0 1
2 6 10 14 18
Re x 10e3
Fig. 1. Schematic layout of roughened solar air heater test Fig. 2. Energy balance for solar collectors having roughened
section. absorber plates.
D. Gupta et al.
a Curve e/D
Rote of energy gain Rate of energy gain 1 O-023 1 ??lOOOW/m2
2 O-036 Ti = 300K
II lnrolotion
3 0.043 1-O m/s
0.7: 2- “W =
4 o-050
5 Smooth
0.613-
c-
efficiency maximum
)-
FLOW RATE, m
Fig. 3. Energy gain and losses for a solar air heater. 0.5t i-
0.5:2-
0.75 -
I = 1000W/m2
0-U9,
4 8 12 16 20 24
Re x 10’3
b Curve e/D
I = 500W/m2
1 O-023 = 300K
Ti
2 O-036
VW = l-Om/s
3 o-050
4 Smooth
0.68
2 6 10 14 18 22
Re x 10m3
Fig. 4. Efficiency versus Reynolds number.
1501 ’ I I 1 I I I
Fig. 5. Effect of relative roughness height on effective effi-
ciency: (a) insolation 1000 W m-‘; (b) insolation 4 8 12 16 20 24 28
500 W m-‘. (c) Effect of relative roughness height on useful
Re x 10’3
energy collected and power consumption.
Thermohydraulic peforrnance of solar air heaters with roughened absorber plates 31
3.3. Efect of relative roughness height ness height. The effective efficiency increases
Figure 5(a,b) show that the effect of relative with an increase in insolation.
roughness height on the effective efficiency is Figure 6(c) shows a plot of the difference in
insignificant at lower Reynolds numbers (4000 the rate of the useful energy collected (qJ and
to about 7000) in comparison to that at higher the rate of thermal energy equivalent to the
Reynolds numbers. At higher Reynolds num- power consumption (Z’,,,/c) for different values
bers, the effective efficiency increases with a of insolation. It reveals that the maximum for
decrease in relative roughness height for a given the rate of net energy gain shifts to a lower
Reynolds number. It is observed that the Reynolds number as the insolation decreases
Reynolds number at which the effective effi- and hence the point of maximum effective effi-
ciency attains a maximum value depends on ciency also shifts to a lower Reynolds number,
the relative roughness height. The value of the as shown in Fig. 6(a,b). This appears to be due
Reynolds number corresponding to this condi- to a reduction in the plate temperature as the
insolation decreases, lowering the rate of heat
tion decreases as the relative roughness height
transfer to the air, whereas the power expendi-
increases. These figures show that the thermo-
ture remains the same.
hydraulic performance of roughened surfaces is
Figure 7 shows the effect of insolation on the
best if the relative roughness height is the lowest
effective efficiency for a given relative roughness
(e/D=O.O33) and the Reynolds number corre-
height and Reynolds number. The figure shows
sponding to the best thermohydraulic perfor-
that the effective efficiency of roughened solar
mance under a given set of conditions is
air heaters increases with the insolation at
moderate (Re= 14,000) in the range investi-
higher flow rates (Re > lO,OOO), whereas the
gated. At higher Reynolds numbers the perfor-
effective efficiency is constant for all insolation
mance of a smooth air heater is seen to be
values at Reynolds numbers of around 10,000.
better than that of a roughened solar air heater.
At lower flow rates (Re < 10,000) the effective
The reason for the fact that the Reynolds
efficiency decreases as the insolation increases.
number corresponding to the maximum effec-
Figure 7 also reveals that the effect of insolation
tive efficiency shifts to a lower value as the is more prominent at higher Reynolds numbers
relative roughness height increases, may be and the effective efficiency increases at a much
attributed to the decrease of the rate of useful faster rate with insolation as compared to that
energy collected as the relative roughness height at lower Reynolds number.
increases. The rate of useful energy collected
decreases, whereas the friction losses rise with 3.5. Efect ofjlow rate (Reynolds number)
increasing relative roughness height causing The variation of the effective efficiency with
increased energy consumption. Figure 5(c) Reynolds number for different values of the
shows the effect of the increase in relative rough- relative roughness height is shown for insolation
ness height on the rate of net energy gain, i.e. of 1000 and 500 W mm2 in Fig. 5(a,b), respec-
the difference of the rate of useful heat gain tively. The variation of the effective efficiency
and the thermal energy rate equivalent to for smooth solar air heaters is also shown.
mechanical power expenditure. Figure 5(c) These figures show that the effective efficiency
clearly shows that the maximum in the rate of for all the roughened and smooth solar air
the net energy gain shifts to a lower Reynolds heaters attains a maximum value at a particular
number as the roughness height increases for a Reynolds number and reveals the important
given set of environmental conditions. Thus, it fact that, beyond a particular Reynolds number,
can be concluded that an increase in roughness the smooth solar air heater becomes thermohy-
height will result in the shift of the point of draulically more efficient than the roughened
maximum effective efficiency towards a lower solar air heater even though the thermal effi-
Reynolds number. ciency of the roughened solar air heater will
always be higher. Figure 8(a,b) show the ratio
3.4. Efiect of insolation of the effective efficiency of a roughened solar
Figure 6(a,b) show the effect of insolation on air heater to that of a smooth one plotted
the effective efficiency. It is observed that the against Reynolds number for different relative
effective efficiency attains a maximum value for roughness heights at a given insolation. These
each value of insolation and these maxima figures show that the values of the ratio are
correspond to lower Reynolds numbers as the nearly the same for all relative roughness heights
insolation decreases for a given relative rough- for Reynolds numbers between 4000 and about
38 D. Gupta et al.
a b
Curve I (W/m*)
1 1300 T; =300K
2 1100 “;, = 1-O m/s Curve I (W/m*)
I
3 1000
e/D = 0.023
4 800 1 1300 Ti = 300K
5 600 23 1100
1000 VW = 1.Om/s
6 400 0.70 e/D= O-043
O-7( 4 800
5 600
O-6! 0.65 -
E
t
P
O-6(
2 0.60 -
0.55 -
0.50 1 \ 6
I I 1 I I I
8 12 16 20 i 4 8 12 16 20 24
Re x 10’3 Re x 10m3
C = 300K
Ti
VW = l*Om/s
320
e/D = DO36 For I =1300W/m*
280
240
40- pm/c
o-. I I
4 8 12 16 20 24 28 32
Re x 10B3
Fig. 6. Effect of intensity of solar radiation on effective efficiency: (a) relative roughness height 0.023; (b) relative roughness
height 0.043. (c) Effect of intensity of solar radiation on energy balance for a roughened solar air heater.
9000 and thereafter the values start to differ for always be thermohydraulically less efficient than
different relative roughness heights. The limiting a smooth air heater.
value of Reynolds number at which the rough- Thus, it is seen that roughened solar air
ened solar air heater becomes thermohydrauli- heaters are thermohydraulically more efficient
tally less efficient than the smooth solar air if they operate in a given range of Reynolds
heater is indicated by the point where the ratio numbers (say 3500 to 19,000 for a relative
of effective efficiencies attains a value less than roughness height of 0.023 and insolation of
unity. A roughened solar air heater operating 1000 W m-3. However, this range of enhanced
at Reynolds numbers higher than this value will operation depends on the system and operating
Thermohydraulic peformance of solar air heaters with roughened absorber plates 39
a
I Ti =300K l-Re= 4000 Curve e/D I =1000W/m2
Vw =l.Om/s 2-Re = 7000 1 0.023
e/D = O-023 3 -Re = 10000 2 0.036
1.20 4-Re = 14000 3 0.043
5-Re = 18000 4 0.050
I
III I
I
II
8 12 16 20 24
Re x 10e3
1
0.60 -
b Curve e/D I =500 W/m2
1 0.023
2 0.036
*
0.404 -k 3 O-050
500 700 900 1100 1300 c” 1.1
i
INSOLATION, W/m2
Reynolds number
which yield optimum conditions to achieve max- C, specific heat of air, J kg-’ K-’
c constant = rhhtlm~tr~F
imum effective efficiency. D equivalent diameter of duct, m
e height of roughness element, m
F plate efficiency factor
5. CONCLUSIONS f friction factor
friction factor for smooth uassages
On the basis of the above discussion the height of the duct, m _ -
following conclusions can be drawn: convective heat transfer coefficient. W m-’ K-’
intensity of solar radiation, W mei
(1) For each value of the relative roughness thermal conductivity of air, W m-l K-’
height, the effective efficiency attains a max- length of collector, m
imum as the operating flow rate is varied, mass flow rate, kg s-i
mechanical energy consumed for propelling air
and the numerical value of the maximum through collector, W
effective efficiency is found to decrease with pressure drop across collector length, N m-’
an increase in relative roughness height. useful heat gain, W
ambient temperature, K
The optimum operating flow rate shifts to inlet air temperature, K
a lower value as the relative roughness outlet air temperature, K
height increases. plate temperature, K
mean fluid temperature, K
(2) The effective efficiency of a roughened solar mean plate temperature, K
air heater increases as the insolation overall heat loss coefficient, W mm2 K-’
increases for Reynolds numbers higher than velocity of air in solar air heater duct, m s-i
wind velocity, m s-’
10,000. However, at lower Reynolds num- W width of the solar air heater duct, m
bers (Re c 10,000) the thermohydraulic
efficiency decreases with increasing Greek letters
insolation. a angle of attack of roughness elements
(3) An optimum operating condition exists for ncrr effective efficiency
a given insolation and this optimum condi- Y kinematic viscosity of fluid, m* s-i
p density of fluid, kg md3
tion shifts to a higher value of Reynolds (ra) transmittance-absorptance product for absorber cover
number as insolation increases. combination
(4) Roughened solar air heaters are thermohy-
draulically advantageous for lower Dimensionless parameters
Reynolds numbers, whereas a smooth solar e/D relative roughness height
air heater will perform better thermo- P/e relative roughness pitch
e+ roughness Reynolds number [ = (e/D)(m)( Re)]
hydraulically, although the thermal effi- Nu Nusselt number = hD/k
ciency of a roughened solar air heater may Pr Prandtl number = pC,lk
be more than that of a smooth heater, Re Reynolds number = p VD/p
beyond a certain limiting value of Reynolds
number; this limiting Reynolds number has
REFERENCES
been found to lie in the range of
13,000-l 9,000, the actual value depending Altfeld K., Leiner L. and Fiebig M. (1988) Second law
optimization of flat plate solar air heaters-Part I. Solar
upon the actual relative roughness height Energy 41, 2, 127-132.
and insolation. Aftfeld K., Leiner L. and Fiebig M. (1988) Second law
(5) The optimum design conditions for rough- optimization of flat plate solar air heaters-Part II. Solar
Energy 41, 4, 309317.
ened solar air heaters have been determined Cortes A. and Piacentini R. (1990) Improvement of the
for varying relative roughness height and efficiency of a bare solar collector by means of turbulence
insolation and for a relative roughness pitch promoters. Appl.Energy 36,253-256.
Gupta D. (1994) Investigations on fluid flow and heat
of 10 and for an angle of attack of rough- transfer in solar air heaters with roughened absorbers.
ness elements of 60”. The following empiri- Ph.D. thesis, University of Roorkee, India.
cal relationship relates the system and Han J. C., Glicksman L. R. and Rohsenow W. M. (1978)
Heat transfer and friction for rib roughened surfaces.
operating parameters that combine to yield ht. J. Heat Mass Transfer 21, 1143-1156.
the optimum conditions: Han J. C. (1984) Heat transfer and friction in channels with
two opposite rib roughened walls. Trans. ASME, J. Heat
(Re) = 1311.2(Z)“~281(e/D)-o~121. Transfer 106, 774-78 1.
Lewis M. J. (1975) Optimizing the thermo-hydraulic perfor-
mance of roughened surfaces. Ini. .I. Heat Mass Transfer
NOMENCLATURE 18, 1243-1248.
Prasad B. N. and Saini J. S. (1988) Thermo-hydraulic opti-
A area of cross section, mz mization of artificially roughened solar air heaters. Proc.
A, area of absorber plate, mz NSEC, Solar Energy Society of India, Hyderabad.
42 D. Gupta er al.
Prasad K. and Mullick S. C. (1983) Heat transfer character- Vilemas J. V. and Simonis V. M. (1985) Heat transfer and
istics of a solar air heater used for drying purposes. Appl. friction of rough ducts carrying gas flow with variable
Energy 13, 83-93. physical properties. ht. .I. Heat Mass Transfer 28, I,
Sparrow E. M. and Tao W. Q. (1983) Enhanced heat 59-68.
transfer in a flat rectangular duct with streamwise. peri- Webb R. L., Eckert E. R. G. and Goldstein R. J. (1971)
odic disturbances at one principal wall. Trans. ASME, Heat transfer and friction in tubes with repeated rib
J. Heat Transfer 105, 851-861. roughness. ht. J. Heat Mass Transfer 14, 601-617.