A Rano Vitch 1981
A Rano Vitch 1981
A Rano Vitch 1981
E. A R A N O V I T C H
Joint Research Centre, lspra (Italy)
Wlm2.)jm O TRANSPAREN
I C O V E ~
FLUIDTUBE
2000
FRONTLOSSES
(RADtATION+NATURAL
CONVEOT.O.
1500
\
-g}00 ~ - ABSORBERP L
y/
A ~
'"SULA"ON
I BACK LOSSES
500~
UV • V IR
I
0.5 1 1.S 2 9.5 jam EDGE LOSSES
Fig. 1. Spectral distribution o f solar radiation at sea Fig. 2. Schema o f a fiat plate collector. The green-
level (very similar to a blackbody spectrum at 6000 house effect is very important for the efficiency o f
K). the system.
32
o[..[Wl
'.
~'
I J I [] [I I
~N/.~[
'~1 /I~ J
t , , / , . . ~ . ,..., ='.. :.. L.L d. .....
1.0
70 °C, a glass cover at a temperature of 40 °C
3 4 5 6 ? $ 9 10 11 12 13 1,; 15~
with a spacing o f 2.5 cm.
(b) WAVELENGTH (MICRONS)
g = 981 cm/s 2
Fig. 4. Absorption spectrum for "Mylar" W-2 film
(a) low range, (b) high range. d = 2.5 cm
34
25 -- 3,5 9
cc
4
N 2O
2
0
I I HORZIONTAL
16
}
1S
4
~3 J I
5Q 100 TEHPERAllJREoc 2
45 °
Fig. 5. Kinematic viscosity and thermal conductivity I I I
"0
o f air.
16--
/~ T = 70 - - 4 0 = 30 °C
70 + 40
T- + 273 = 328 K
2
v = 17.5 X 1 0 - 2 cm2/s I I I
10 4 10 5 10 6 10 7 _'~r
(981)(2.5)s(30)
Gr = = 45670. Fig. 6. Natural c o n v e c t i o n in inclined cells.
(17.5 × 10-2)~(328)
Nu = 0.062(Gr) °'s27 (14)
The heat transfer coefficient U~n between the
absorber plate and the cover plate can be cal- for s = 9 0 ° and 1.5 × 105 < Gr < 107
culated if the function which relates the
Nusselt number, Nu, to the Grashof number, Dropkin-Somerscales recommend:
Gr, and the slope, s, is known. A great
Nu = [0.060 - - 0.017 (s/90)] (Gr) 11s (15)
n u m b e r o f experimental correlations have
been proposed for the calculation o f heat with Gr > 2.105
transfer by natural convection in closed
inclined cells. It can be noted t h a t there
exists a dispersion o f the order of 20% It can be observed that the Dropkin-Somer-
between the different formulas recommended scales formula gives somewhat lower values
by authors such as Tabor, de Graaf, van der than the Tabor formulas which might be con-
Held, Dropkin, Somerscales, etc. sidered more reliable because t h e y were ex-
In Fig. 6 are represented correlations perimentally established with air, whereas the
recommended by Tabor and Dropkin-Somer- former was established with liquids.
scales for different inclinations. Moreover, the Dropkin-Somerscales for-
Tabor recommends: mula is recommended for Grashof numbers
Nu = 0.152(Gr) °'2sl (11) superior to 1.5 × 10 5 where, as in solar col-
lectors, the Grashof number between absorber
fors=0and104 <~Gr<107 plate and glass cover is generally inferior to
10 5.
Nu = 0.093(Gr) °'sl° (12) On the other hand, the Dropkin-Somerscales
formula has a great merit for simplicity and
for s = 45 ° and 104 < Gr < 1 0 7
can easily be used in design applications. So,
after changing the values of the constants in
Nu = 0.033(Gr) °'ssl (13)
order to be nearer to the values recommended
for s = 90 ° and 1.5 X 104 < Gr < (1.5)105 by Tabor, one obtains:
35
Nu = [0.078 -- 0.026 (s/90)] (Gr) l/a (16) This example shows that for non-selective sur-
faces, heat losses by radiation are predom-
with 104 < Gr < 105 inant over heat losses by natural convection.
= 96 W/m 2 . Fig. 7. H e a t t r a n s f e r b y n a t u r a l c o n v e c t i o n .
36
BACK LOSSES AND EDGE LOSSES design purposes, to use a single formula which
accounts globally for radiative and convective
Back losses and edge losses, for conve- losses.
nience, are generally characterized b y a single Two correlations proposed by Tabor and
heat transfer coefficient, Ub. In a poorly MacAdams are represented in Fig. 8. The
designed solar collector, the edge losses can linear equation for MacAdams:
become o f significant importance (in the case,
for instance, o f a metallic contact between Ua = 5.7 + 3.8V (25)
the absorber plate and the external casing). where
In a properly designed collector the edge
losses should be kept inferior to 5% o f the Ua = heat transfer coefficient from cover to
total losses and the heat transfer coefficient, the ambient air (W/m 2 °C)
Ub, can be represented with an acceptable V = wind velocity {m/s)
approximation b y :
fits well with the values recommended b y
Vb -~ w (21) Tabor, when the wind velocity is in the range
eb o f 2 - 5 m/s.
where
~a(w/m2oC)
)'b = thermal conductivity of the insulating 50
material
e b = thickness o f the insulating material.
~ ~AC AOA~
The corresponding losses are then:
I
qb "" hb(Tp -- Ta) (22)
fl
I
.I I
~, -- Ul(Tm -- 7'.)
tanh mL
X x+dx = , ~ (32)
Xm
m L [1 + ( ~ ) tanh mL ]
~kb
Ub ~ -- (41)
eb HEAT LOSSES w i r e 2
(j.t)2 I
UI = I 1 + Ub (46)
u, + uo~ ua
In Fig. 11 the overall heat transfer coefficient,
derived from the preceding example has been
represented as a function o f (Tp -- Ta).
Although it varies slightly with temperature,
( ,.~12
for the characterization of a collector, a con-
stant average value will be assumed.
Ut W/m2°C
Fig. 12. Transmission through one cover.
U| ...-
where n is the refractive index of the trans-
parent medium (the refractive index of air is
1).
At the second interface, the quantity (1 --
-- TG~ 5 o c ; s.45°;V. 3m/s p)p is reflected and the quantity (1 - - p ) 2 is
£ p , 0.96; £ c . 0 . 8 6
transmitted. Summing up the transmitted
terms after the multiple reflections, the
transmission coefficient, rF, becomes:
I i I I I" I I
10 20 30 40 50 60 70 ( Tp-To ]
oC TF = (1 __p)2 + (1 __p.)2op2 + (1 _ _ p ) 2 . p 4 +
Fig. 11. Overall heat transfer coefficient: single-glass
black paint collector. 1--p
+ .... - - (48)
l+p
(n - - 1) 2 1 - - 0.041
p(0) - (n + 1) 2 (54) (T,)2 - - 0.921
1 + 0.041
K is t h e e x t i n c t i o n c o e f f i c i e n t w h i c h can v a r y
Example 6. f r o m 0 . 0 4 / c m f o r an e x c e l l e n t glass t o 0 . 3 2 /
Calculate t h e Fresnel t r a n s m i s s i o n c o e f - c m f o r a p o o r glass; L is t h e a c t u a l p a t h o f
ficients in t h e case o f 1, 2 a n d 3 glass covers, radiation through the medium.
at n o r m a l i n c i d e n c e , a s s u m i n g t h a t t h e refrac-
tive i n d e x o f glass is 1.53. F r o m e q n . ( 5 6 ) : E x a m p l e 8.
Calculate t h e a b s o r p t i o n in t h e case o f 1, 2
2 × 1.53
a n d 3 glass covers, at n o r m a l incidence, as-
N = 1 r F a ( 0 ) = 2 × 1.53 + 1 ( 1 . 5 3 - - 1 ) 2 - s u m i n g t h a t t h e e x t i n c t i o n c o e f f i c i e n t is e q u a l
to 0.18/cm and that the thickness of each
= 0.92
glass c o v e r is 4 m m . A t n o r m a l i n c i d e n c e t h e
2 × 1.53 a c t u a l pass is e q u a l t o t h e t h i c k n e s s :
N=2 TF.2(0) = N = 1 ; L = 0.4 c m ; Ta = e -(°'~sX°'4) = 0 . 9 3
2 × 1.53 + 2(1.53 - - 1) ~
E q n . (48) gives:
E x a m p l e 10.
1 --0.312 Calculate t h e c o m b i n e d t r a n s m i s s i o n coef-
(Tr)1 = - 0.524
1 + 0.312 ficient f o r 1, 2 and 3 glass covers o f 4 m m
41
l OUTLET
9: EFFICIENCY -1
n
Y
T
' - I ~MyABB COVER
(BLACK PA,"T)2 ~LASSiICOVER ~S~" -.
.2--
I
I I II I I I
0.02 0.04 0.06 0.08 0.1 0.12 [ Tp-To )/1
~ INLET
Fig. 16. Instantaneous efficiency curves. Fig. 17. Temperature distribution in flow direction.
43
It can be seen from equation (65) that the and there is little to be gained in that direc-
temperature profile in the flow direction is tion. The increase of efficiency can be ex-
not linear except for small values o f the coef- pected essentially by the reduction of the
ficient ~. Introducing the average fluid tem- heat losses and various solutions are presented
perature, TIn, and the outlet fluid tempera- here. It should be noted that these solutions
ture, Tfe, one obtains: generally imply a certain deterioration in the
(Ta).I -- UI(Tm - - T.) 1 absorption process and one should be careful
- -- [1 - - e -F'~] ( 6 7 ) to verify that the overall effect is positive.
(rcx)~I -- U,(Tti -- T.) F'{3
Tti + T,e
Tm--
2 qf = ha(Tcl -- T a) (70)
The mean fluid temperature. Tin. is very near q! = egl(Tc42 -- T41) + [1 -- 0.33(s/90)] ×
the average temperature derived from the
inlet and outlet temperature, except for very X ( T c 2 - - T e l ) 4/s (71)
small mass flows. (~f = e g 2 ( T 4 - - T 4 2 ) +
In the preceding sections, the heat losses
and efficiency of a " n o r m a l " collector have + [1 - - 0 . 3 3 ( s / 9 0 ) ] (Tp - - T o 2 ) 4/s (72)
been determined. The next sections will be qb = Ub (Tp -- Ta) (73)
devoted to various ameliorations which can
be considered to increase the efficiency of a ql = q~ + qb (74)
solar collector.
As it was said earlier, two distinct processes egl - (75)
2/ec -- 1
are involved: the absorption of energy and
the inevitable losses of energy which increase 1
with temperature. As far as the absorption of eg2 - (76)
1 1
energy is concerned, the only improvements --+ ----1
which can be expected concern the trans- ep ec
mittance of the transparent covers. Some new
types o f glasses, still very expensive, with ~b
Ub = - (77)
better transmission properties, are becoming eb
avaflabld. The absorption coefficient of black
paint absorber plates is o f the order of 0.96 Us = 5.7 + 3.8V (78)
44
The electromagnetic theory of light shows but it seems that the two most promising
that metals are good reflectors o f long wave solutions from an economic and industrial
radiation, the reflectivity being a function of point of view are the so-called "black nickel"
the wavelength and the electrical conductivity: and "black chrome" solutions.
In Fig. 22 the heat losses and the efficiency
r~. = 1 -- 0.365(),x) -1/2 (82)
curve of a one-glass selective collector have
k = wavelength in pm been calculated, assuming:
x = conductivity in ohms per mm 2 per m ep = 0.15 and a = 0.90.
length
For a black paint collector, the radiative
rx = reflectivity. losses are predominant, here it is the convec-
tive losses which are the most important with
The reflectivity of good conductors with the following repartition:
polished surfaces such as copper, for instance, - -radiativelosses: 2 0 %
is found to be between 0.95 and 0.98 at a - -convective losses: 5 4 %
wavelength of 2 p m and near 1 for very long - -back losses: 2 6 %
wavelengths. W h e n the convective losses are dominating it
For a given wavelength ),, the relations: can be of a certain interest to consider evacu-
ated collectors, where by reduction of the
r~ + ~ = 1 (83)
pressure, the air is maintained stagnant inside
the collector (Gr < 1500). S o m e attempts
an = e~ (84)
have been made to manufacture fiat plate
show that to a reflectivity of 0.95 corresponds evacuated collectors. F r o m a technological
an emissivity, en, o f 0.05. point of view it is more realisticto consider
Electrolytical, chemical or electrochemical cylindrical geometries as in the Phillips' col-
s u r f a c e treatments can deposit films which lector for instance.
are visibly black and at the same time trans-
parent to wavelengths above 2 pm. F o r Reflective transparent covers
instance black deposits are produced indus- A n o t h e r means of reducing the global
trially b y electroplating in particular of nickel emissivity coefficient b e t w e e n the absorber
in the presence of zinc and sulphides. These plate and the transparent cover, is to increase
deposits on a bright metal base will show in the reflectivity of the transparent cover in the
the visible spectrum an absorption of 80 - infrared spectrum. It has been observed for
90% and retain the emissivity of the metal instance that deposits of indium-oxide, In2Oa,
base in the infrared spectrum. on glass will increase the reflectivity of the
An ideal selective surface will have the
profile shown in Fig. 21.
Many different combinations of coatings
and metal substrates have been investigated
HEAT LOSSES ( w/m 2 ) (wire 2 OC)
500
SELECTIVE S U R F A C E / 1 , UL
~
1
~,.0.95
300 m
it
200
100
j ~ q r
o 0 60 o
6
HONEYCOMB STRUCTURES 4 NO HONEY COHB
~ /..~'/~. ~o o ° r i d - 2 oo~L/d.5
2
By inserting "honeycomb" structures, also f ~o v -
oo
o
I~ n ° g ° 0 ° ° ° ¢ ° 90°°°°°
called Francia structures, between the absorber o o /.SO
plate and the transparent cover, both the heat 6
of heat losses. 4 °o
some cases becomes more important than the signed collectors where heat losses are mini-
heat transfer in the absence o f honeycombs. mired will lead to acceptable performances.
The relative transmission of various honey- The choice of the basic materials for the
comb materials as a function of incident absorber, the transparent cover and the insu-
angle is presented in Fig. 24. It can be seen lating material will play an important role in
that for paper and aluminium honeycombs the determination of the quality of the col-
the relative transmission falls too rapidly to lectors. However, it should not be forgotten
permit practical applications. Even in the case that the quality of a collector is not only
of Mylar honeycombs, for incident angles based on its thermal performance but also
greater than 30 °, the reduction in transmission on its cost and durability. The final product
is not negligible. should be the result of a careful balance
The transmission-absorptance product for a between these three aspects.
collector with and without a honeycomb is
shown in Fig. 25. By comparisons with curves
BIBLIOGRAPHY
of Fig. 13, the influence of the honeycombs
can be evaluated. Arnold, J. N., D. K. Edwards, and I. Carton, Effect
of tilt and horizontal aspect ratio on natural con-
LO
-- I ~ I . vection in a rectangular honeycomb, Trans. ASME,
O.G
J. Heat Transfer, 99 (1977) 120.
Buchberg, H., I. Catton, and D. K. Edwards, Natural
ie 0.8
convection in enclosed spaces: A review of appli-
cation to solar energy collection, "Plans. ASME, J.
0.7
Heat Transfer, 98 (2) (1976) 182.
Charters, W. W. S. and L. J. Peterson, Free convection
suppression using honeycomb cellular materials,
NORMAL R L ECTI V Y GLASS Solar Energy, 13 (1972) 4.
~ 0.5 Randall, K. R., J. W. Mitchell, and M. M. El Wakil,
Natural convection characteristics of flat-plate
i 0.4 collectors, In J. R. Howell and T. Min (eds.),
Heat Transfer in Solar Energy Systems, American
~"0,3 Society of Mechanical Engineers, New York, 1977.
Beckman, W. A., S. A. Klein, and J. A. Duffle, Solar
Heating Design, Wiley, New York, 1977.
Hottel, H. C. and B. B. Woertz, Performance of fiat-
0.1
plate solar-heat exchangers, Trans. ASME, 14
(1942) 91.
I I I Lin, R. J. H., Optimization of coatings for fiat plate
15 $0 45 QO solar collectors, Report C00/2930-4 to ERDA
ANGLE OF INCIDENCE, DEGREF..S (Jan. 1977).
Fig. 25. Transmission-absorptance product for a McDonald G. E., Spectral reflectance properties of
single-glazed flat-plate solar collector incorporating a black chrome for use as a solar selective coating,
honeycomb convection suppressor (from NBS 899). NASA Tech. Memo NASA TMX 0171596 (1974).
899). Tabor, H., Selective surfaces for solar collectors, In
Low Temperature Engineering Appl|cations of
Solar Energy, ASHRAE, New York, 1967.
CONCLUSIONS Hottel, H. C. and B. B. Woertz, The performance of
fiat-plate solar-heat collectors, Trans. ASME, 64
The complex nature of the photothermal (1942) 91.
conversion solar energy which involves, on Whillier, A., Solar energy collection and its utilization
for house heating, ScD. Thesis, MIT, 1953.
one hand, the collection and the absorption Klein, S. A., Calculation of fiat-plate loss coefficients,
o f sunlight and, on the other hand, its trans- Solar Energy, 17 (1975) 79.
formation into useful energy, has been illu- Tabor, H., Radiation, convection, and conduction
strated by the description o f the basic thermal coefficients in solar collectors, Bull. Res. Council
of Israel, 6C (1958) 155.
and optical phenomena taking place in flat- Barley, C. D. and C. B. Winn, Optimal sizing of solar
plate solar collectors. Because o f the low collectors by the method of relative areas, Solar
densities o f solar fluxes, only carefully de- Energy, 21 (1978) 279.