Honeycomb solar pond collector and storage system

Honeycomb solar pond collector and storage system

Efwgy Printed Vol. 8. No. II. pp. 883490. io Great Britain. 1983 HONEYCOMB 9 SOLAR POND COLLECTOR STORAGE SYSTEM 0360-SW/83 33.00 + .I0 1983 Fm-...

460KB Sizes 9 Downloads 138 Views

Efwgy Printed

Vol. 8. No. II. pp. 883490. io Great Britain.

1983

HONEYCOMB

9

SOLAR POND COLLECTOR STORAGE SYSTEM

0360-SW/83 33.00 + .I0 1983 Fm-gamon Press Ltd.

AND

N. D. KAUSHIKA,~ M. B. BANJZRJEE, and YOJNA KArrt Centre of Energy Studies, Indian Institute of Technology, Hauz Khas, ND-I 10016, India (Received 2 1 September 1982)

Abstract-An analysis of a honeycomb-stabilized, saltless solar pond as a solar energy collector and long term (spanning seasons) storage system is presented. The solar pond is considered with a nonconvective zone made up of an oil layer and air honeycomb configuration. A heat flow model is developed using the two loss mechanisms (conduction and radiation). The efficiency of heat collection and the storage characteristics of the system are excellent for hot water production and process heat applications. NOTATION pond area, mr fraction of the pond area occupied by honeycomb material geometrical shape factor for radiation loss depth of the bottom convective zone, m upward heat loss coefficient, W/m*-K bottom heat loss coeBicient, W/m*-K sides heat loss coefficient, W/mZ-K radiation heat loss coefficient, W/m*-K thermal conductivity of air, W/m-K thermal conductivity of silicon oil, W/m-K thermal conductivity of side and bottom insulation, W/m-K thermal conductivity of the core material used for side and bottom insulation, W/m-K thermal conductivity of the honeycomb material, W/m-k thickness of the side and bottom insulation covers, m thickness of the core material used for side and bottom insulation, m highest number of harmonic m considered retrieved heat flux per unit of pond area, W/m* solar insolation at time 1, W/m* average value of the solar insolation W/m* amplitude of the mth harmonic of S(r) W/m* ambient air temperature at time f, K average value of T”(t), K amplitude of the mth harmonic T,,(r), K transmittivity of the oil layer for the solar incident angle 0 transmittivity of the honeycomb structure for the solar incident angle 0 average value of the sky temperature, K water temperature in the bottom convective zone at time I, K amplitude of the mth harmonic of T,(r) depth of the honeycomb structure, m total heat loss coefficient W/m*-K absorptivity of the convective zone for the solar incident angle 0 effective absorptivity-transmissitity product of nonconvective zone and including the oil-water interface heat capacity per unit volume of water, J/m3-K Stefan-Botlzamm constant emissivity of the oil surface INTRODUCTION

ponds have received worldwide attention in the last decade as an alternative source of low temperature heat for industrial process heating, space heating and electrical power generation. The most common version of the solar pond is the saline pond,12 which has many attractive features such as relatively low cost, large as well as long term storage capacity and simple technology. However, it possesses some operational problems such as rather complicated maintenance of the salt gradient, the disturbance of the salt gradient due to wave motion caused by winds, corrosive properties of salt water, etc. Recently, Ortabasi er al.3 have proposed that these problems may be circumvented in a honeycomb-stabilized, saltless solar pond. These authors have demonstrated the technical feasibility of the concept tsupported

by the Solar Energy Research Centre, University of Queensland, Brisbane 4067, Australia. 883

N. D. IC\L'SHKi et al.

884

in a laboratory model pond and have indicated that an air honeycomb structure floating on a thin layer (1 cm) of silicone oil can indeed serve as the nonconvective zone of the solar pond. In the present paper, we evaluate the solar collector and storage characteristics of a honeycomb solar pond and investigate its applications. SYSTEM

CONFIGURATION

The honeycomb solar pond is envisioned to have its nonconvective zone about 11 cm thick and a bottom convective zone of about 0.5 m depth. In a pond of such shallow depth, the side and bottom heat losses through the surrounding earth region would be very large unless insulation is provided to inhibit the loss. The configuration of a practical system would resemble a rectangle box, as is illustrated schematically in Fig, 1. The side walls and the floor are packed at the cores with insulation. The layers of thickness I3 on either side of the side and bottom insulation medium are made of material such as ferrocement to protect insulation from erosion by the ground soil and water. The bottom convective zone (water pool) is either blackened at the sides and bottom walls or else the water mass itself is blackened. HEAT

FLOW

MODEL

The incident solar radiation is intercepted and filtered out by the honeycomb before it falls on the oil layer. A part of it is reflected at the air-oil interface. The remainder is attenuated throughout the pond depth and is completely absorbed in the bottom convective zone (the zone used for heat collection, storage and extraction). The solar radiation reaching is the effective the bottom convective zone is (arM(t), where (ar)eff absorptivity-transmissivity product corresponding to the nonconvective zone and including the oil-water interface. It may be expressed as

(1)

SOIL CONVECTIVE

ZONE

Fig. I. Schematic diagram of the honeycomb solar pond.

Honeycomb solar pond collector and storage system

885

The thermal energy balance in the convective zone may be written as

The only approximation in writing Eq. (2) is the assumption that there is negligible heat storage in the insulation. The solar intensity and atmospheric air temperature variations with time are periodic and are represented by the Fourier series S(t)=Re

TA(t) = Re

(

>

(3)

TAO+ i TAmeimwr,

(4)

So+~.Smehw’

m

m

,

>

the real part of the following means quantity and where Re w = 2n/(l x 365.25 x 24 x 60 x 60) set-’ for the annual cycle variation. Consequently, T,(t) will also be periodic and T,( t ) = Re

Td + i

T,, eimt .

m

Equation (2) in combination with Eqs. (3)-(5), leads to the following expression for the heat flux: Q(t) = (zr),&, - U,(T,, - T,,,,) + (a~).,; -i

S,,, eimo’

(imop,h + U,)T,, e’““‘+ Uyi m

TAmeimor.

(6)

m

The average efficiency of heat extraction from the pond may be expressed as tl

=[SQ(f,dr]i[lS(i)dr].

(7a)

Evaluating the integrals over the annual cycle, we have

0)

rl = (ar )eir- u,(A T/S,),

where AT is the difference between the storage zone temperature and the atmospheric air temperature. If it is assumed that the honeycomb structure suppresses convection totally, the heat loss from the hot reservoir (bottom convective zone) can be evaluated by using as the two loss mechanisms only conduction and radiation. The conductive heat loss may be found from first principles; since edge effects will not be negligible in most cases, we have incorporated appropriate geometrical shape factors4y’ to allow for these effects. Following Ref. 4, the conductive heat-loss coefficients are given by h, = f [(Ak,/t,) + 2.16/k: + 0.60tJc,],

h, = f [(Ak,/t,) + 2.16& + 0.601,&],

(9)

h, = f [(4fhk,/t,) + 2.16hk, + 0.60&k,],

(10)

where k, = k,kdt, + tJ/(t,kl + tzk,), k, = A,kh + (1 - A,)k,, EGY Vol. 8 No. II-E

t, = t, + t,, r, = t, = f, + l,,

886

N. D. KAUSHIKAer al.

+ I&,), k: = 2/[( I/k,) + (l/k,)], 1 = J(A ) (for a pond with square cross-section). The radiation losses from the oil surface to the sky are calculated from the Stefan-Boltzmann law, again with the use of a shape factor to account for radiation trapping in the honeycomb structure. Following Holland6 and Tien and Yuen,’ the radiation loss coefficient is given by k, = k, = k,k.,(f, + f,)/(&,

h, = F’a

(T,” - T,4)/(T, - T,),

(11)

where l/F’= l/F+ (1 - LJ/E~. The shape factor F is taken from Ref. 8; T, is the temperature at the oil surface and may be expressed as

7’1, = Tco- (T, - T,o>(~,lk,>l([(k,l~,> + M - ’ + (OJ}.

(12)

Thus h, may be calculated from Eqs. (11, 12) by an iterative procedure similar to that of Ref. 9. Finally, the total heat loss coefficient is given by UL=h,+hB+h,+h,.

(13)

The ratio R of heat losses through the side to those through the surface is R = h,/(h, + h,). CALCULATIONS

AND

(14) RESULTS

The solar collection and storage behaviour of the pond may be judged from its year-round thermal performance. We have, therefore, made calculations for the retrieved heat flux and the efficiency of heat retrieval at constant temperature (T,, = 0) using Eqs. (6) and (7), respectively. The data on annual variation of solar intensity and atmospheric air temperature have been taken from Kaushika ef al.‘O and correspond to New Delhi (latitude 28”N) during 1974. The values of T,(8), T,(8), a(e), and the resultant (ar),g are approximately constant and correspond to a solar angle of incidence on the horizontal surface at 2 p.m. at equinox. T,(8) has been calculated by using expressions given by Holland,6 T,(8) has been derived from the solar absorption data of silicone oil supplied by Dow Corning Australia Pty. Ltd. The resultant value of (r~r)~~is 0.72. The heat-loss coefficients have been evaluated using the Eqs. (8)-(13). The values of parameters

Table 1. Variations of U, (W/m?-“K) with A and h.

50

2.18

2.23

2.30

2.37

2.43

100

2.17

2.21

2.25

2.30

2.35

200

2.16

2.19

2.22

2.25

2.29

500

2.15

2.17

2.19

2.21

2.23

1000

2.15

2.16

2.18

2.19

2.21

Honeycomb

solar

pond

collector

and storage

system

887

characterizing the pond geometry, as well as other thennophysical parameters used in these calculations, are as follows: I, = 0.10 m, l4 = 0.075 or 0.15 m, k, = 0.034 W/m-K, ES= 0.94, kz = 0.16 W/m-K tz = 0.01 m, k, = 0.035 W/m-K, (silicone oil), G = 5.669 x lo-* W/m-KJ, A, = 0.1363, A = 0.1 to 1000 m’, k, = 0.21 W/m-K, h = 0.10 to l.Om, l,=O.Ol m, r,,,,= 313-363 K, k, =0.42 W/m-K, F = 0.175 (for a cell size of 0.017 x 0.017 x 0.10 m). The variations of R and UL with geometric parameters (area, depth of the bottom convective zone) of the pond are illustrated in Figs. 2 (a, b) and Table 1, respectively. For smaller and deeper ponds, thicker side insulation is necessary to render the side losses insignificant coinpared to the heat losses through the surface. The variation of heat flux that can be retrieved from a honeycomb solar pond during year-round operation is shown in Figs. 3 (a, b). Corresponding annual average efficiencies of the honeycomb solar pond are compared with those of a salt gradient solar pond in Table 2. The honeycomb pond efficiency as a function of collector parameter (AT/So) is depicted in Fig. 4; for comparison, typical curves for the salt gradient solar pond collector’.” and the flat plate collector’* with two covers and selective coating on the absorber are also shown. The honeycomb solar pond is superior to the other collector/storage systems for heat extraction at temperatures lower than 65°C. The honeycomb solar pond is, however, less attractive for power production because organic Rankine cycle turbines are often used in power production systems. These turbines require source temperatures higher than 70°C and a source to sink temperature difference of 50°C or more.

0

10

20

30

40 POND

50

60

AREA

(r-n21

Fig. 2(a).

70

80

go

100

WI0

N. D. ~~J~HIICA

888

et al.

,

0.G I4 I 0.15 m

0.

R

0.

DEPTH

OF CONVECTIVE ZONE

0.

PONO

AREA

(rn’)

(b)

Fig. 2. Curves for the variation of R: (a) I, = 0.075 m, (b) l4 = 0.15 m.

Table 2. Heat extraction efficiencies (annual average) for a solar pond with honeycomb and salt gradient. *at extraction temperature, %

Efficiency of honeycomb solar pond,%

with edge losses

without

T

C)ptimm

efficiency salt gradient pond 5mlar ( Ref. 9))

cIf

I

k

edge 1ssses

40

55

57

39

so

45

49

35

60

35

41

32

70

25

33

30

80

14

25

28

90

4

17

27

Honeycomb

889

solar pond collector and storage system

-801

I (cl)

-25 J

F

M

A

M

J

J

A

S

0

N

MONTHS (b)

Fig. 3. Annual variation of Q(f): (a) the pond has side losses corresponding the side losses in the pond are set equal to zero.

to I, = 0.075 m, (b)

890

N. D.

et al.

khISHlKA

AT (So:220 0.91

0

W/m')

22 I

LL, I

66

0.1

0.2

0.3

88

110 I

I

0.4

0.5

J

0.6

AT /So Fig. 4. The efficiency as a function of the collector parameter AT/S,,. Curve A : flat plate collector with 2 glass covers and selective coating on the absorber. Curve B: salt gradient pond with a nonconvective zone of I m. Curve C: honeycomb solar pond with side losses corresponding to 1, = 0.075 m. Curve D: honeycomb solar pond with side losses set equal to zero.

REFERENCES I. 2. 3. 4. 5.

6. 7. 8. 9.

10. II. 12.

H. Tabor, Solar Energy 7, I89 (1963). H. Weinberger, Solar Energy 8, 45 (1964). U. Grtabasi, F. H. Dyksterhuis, and N. D. Kaushika, Solar Energy 31, 229 (1983). I. Langmuir, E. Q. Adams, and F. S. Meikle, Trans. Am. Efectrochem. Sot. 24, 53 (1913). L. M. K. Boelter, V. M. Cherry, and H. A. Johnson, Heat Transfer. University of California Press, Berkeley (1942). K. G. T. Hollands, Solar Energy 9, 1.59 (1965). C. L. Tien and W. W. Yuen, ht. J. Heat Mass Transfer 18, 1409 (1975). H. C. Hottel and J. D. Keller, Trans ASME 55, 39 (1933). M. S. Sodha, N. D. Kaushika, and S. K. Rao, Energy Res. 5, 321 (1981). N. D. Kaushika, P. K. Bansal, and M. S. Sodha, Appl. Energy 7, 169 (1980). C. F. Kooi, Solar Energy 23, 37 (1979). J. A. Duffie and W. A. Beckman, Solar Engineering of Thermal Processes. Wiley, New York (1980).