A User's Manual for DELSOL3 - prod.sandia.gov - Sandia National ...

More documents

Recommendations

Info

34 N' Figure 11-6. Code Defined Field (IUSERF = 1). RADMIN = 0.80; RADMAX = 7.15; NAZM = 12; NRAD = 11. Letting I = 1, NAZM, boundaries are defined as follows: I NRADMN(1) NRADMX(1) 1,2,12 1 3,11 1 4,lO 1 599 1 fh8 1 7 1 11 10 9 8 7 6
To directly specify such a field, the following variables must be defined: (1) RAD- MIN, RADMAX, NRAD, NAZM, INORTH, and AMAXN, to set up the zoning; and (2) NRADMN, NRADMX, DENSIT, AZMSEP, and FLAND, to character- ize the field zone by zone. In addition, for a performance calculation on a system optimized by DELSOL (where the results were saved on a file by specifying IO- TAPE=l in namelist OPT), all of the above variables can be read from the file by specifying ITAPE=3. In this case, IUSERF=2 will be the code default, and should not be set by the user to any other value. For the Lth azimuthal zone (L=l, NAZM) as numbered in Figure 11-4, all radial zones are occupied from the minimum radial zone number, NRADMN(L), through the maximum radial zone number, NRADMX(L). If no zones are occupied in the Lth azimuthal zone, then NRADMN(L)=NRADMX(L)=O. The radial/azimuthal zone boundaries may not exactly match the boundaries of the user’s field. The FLAND array can be used to trim the DELSOL zoning. If there is a land constraint, FLAND will be calculated automatically by DELSOL. In the absence of a land constraint, the user may specify FLAND. For the (K,L) zone, FLAND(K,L) is the fraction of the land area in the (K,L) zone that is occupied by the heliostat field. If the whole zone is occupied, FLAND(K,L)=l.O; if half the zone is occupied, FLAND(K,L)=0.5, etc. Field averaged performance is calculated and reported using this option. This option cannot be used for any initial performance run which will be used later for a field buildup (optimization) run. d) IUSERF-3 This option allows the user to specify the x and y (east and north) coordi- nates of the base of every heliostat relative to the tower base. For the performance at a single time (IPROB=2, namelist BASIC), an asymmetrical heliostat field can be used. However, in order to calculate daily or annual performances the field must be symmetric about the N-S axis. A special convention is used to group and number the heliostats. The heliostats are grouped into “rows” as illustrated in Figure 11-7. In a field that sur- rounds the receiver the rows will usually be completely or partially filled circles. In a north-only field the rows will be arcs or lines. The rows do not intersect. The rows are numbered starting with the row nearest the tower and proceeding outward. Within each row the heliostats are numbered starting with the heliostat on the N/S line or just east of the N/S line. The numbering increases in a clock- wise manner around the tower. Note that for a line or arc of heliostats (see row 4 in Figure 11-7) the number starts in the middle, proceeds to the eastern edge, goes to the western edge and then heads to the middle again. The code considers the shading and blocking by only those heliostats within f two rows of the row in which the heliostat of interest is located. For any of the options described above, the number of zones and hence the accuracy increases as NRAD and NAZM increase. The tradeoff is that computing time and cost will also increase. The variation in execution time is approximately linear with the number of zones, while the increase in accuracy with the number 35
Page 1 and 2: - SANDIA REPORT SAND86-8018 Unlimit
Page 3: SAND 86-8018 Unlimited Release Prin
Page 6 and 7: ' The computer code which this repo
Page 9 and 10: I. Introduction CONTENTS A. Differe
Page 11 and 12: 1. More Accurate Images from Canted
Page 13: VII. Comparison of DELSOL, MIRVAL,
Page 16 and 17: 111-1 Insolation Models 111-2 Atmos
Page 19 and 20: A USER’S MANUAL FOR DELSOLS: A CO
Page 21 and 22: owing, blocking, atmospheric attenu
Page 23: The fourth major enhancement in DEL
Page 26 and 27: W- 26 S ? Jr RECEIVER NORMAL COORDI
Page 28 and 29: ZONE FlEL 180" Figure 11-3. Surroun
Page 30 and 31: 30 ZONE 60 Figure 11-4. Method of Z
Page 32 and 33: Figure 11-5. North-only Field (INOR
Page 36 and 37: 36 Figure 11-7. Schematic Diagram o
Page 38 and 39: -SLEW (1)- Figure 11-8. Example of
Page 40 and 41: Figure 11-9. Radial Stagger Arrange
Page 42 and 43: Note that when IDENS=l or 2 the azi
Page 44 and 45: J1.D. Heliostats Either rectangular
Page 46 and 47: Given that heliostats will be mass
Page 48 and 49: (A) FLAT REFLECTED RAY FROM (B) FOC
Page 50 and 51: IROUND = 0 WM = 9.91 (m) HM = 9.93
Page 52 and 53: That is, the aperture is always loc
Page 54 and 55: controlled by the IAUTOP parameter
Page 56 and 57: only be used with rectangular cavit
Page 58 and 59: grid of points on the DELSOL assume
Page 60 and 61: W E Figure 11-17. Flux Points on a
Page 62 and 63: flux point grid defined during syst
Page 64 and 65: To see whether this problem is occu
Page 66 and 67: Cal cul ati onal Day 1 2 3 4 5 Tabl
Page 68 and 69: 1II.B-I. Position 01 the Sun-The su
Page 70 and 71: where SO is given by Equation (1II.
Page 72 and 73: c 0 w V c LL .r a 0 0 0 0 0 s /s 0
Page 74 and 75: 1II.C. Field Performance Calculatio
Page 76 and 77: 1II.F. Flux Density and SDillane Th
Page 78 and 79: III. G-1. Atmospheric Attenuation:
Page 80 and 81: eceiver) or with aperture area (cav
Page 82 and 83: In other words, PLOST,R is the same
Page 84 and 85:
Following the approach in the previ
Page 86 and 87:
86 - - - LLLLLL www CLee o m 4 I i
Page 88 and 89:
This approximate method of calculat
Page 90 and 91:
RECEIVER STARTUP WEATHER SUN PARASI
Page 92 and 93:
111.1. Relationship Between Perform
Page 95 and 96:
IV. System Optimization Calculation
Page 97 and 98:
power and annual energy. Section 11
Page 99 and 100:
, . IV.A-7. Calculation of Annual E
Page 101 and 102:
. production. When each design powe
Page 103 and 104:
large enough fraction of total capi
Page 105 and 106:
_____ Table IV-2. Design Parameters
Page 107 and 108:
DELSOL will examine NUMHTW (520) eq
Page 109 and 110:
The final step of cavity receiver o
Page 111 and 112:
. furthermore, that the assumed tur
Page 113 and 114:
point in the NMXFLX(1) array (I=l,
Page 115 and 116:
, flux level. DELSOL determines the
Page 117 and 118:
optimum for any power level is 13.
Page 119 and 120:
I V. System Costs and Economics Whi
Page 121 and 122:
for meteorological equipment. The d
Page 123 and 124:
0 co I I I 0 co J 0 c\! c, 0 3 Ln 0
Page 125 and 126:
w/2 x RWCAV W/2 x RWCAV IN (180. -
Page 127 and 128:
S M ~ , = ~ solar F multiple for re
Page 129 and 130:
nSTOR is determined from an assumed
Page 131 and 132:
CTK,REF (CSTREF) = $9.70 x loG VTK,
Page 133 and 134:
XHE,A = scaling exponent. ni is cal
Page 135 and 136:
. CCFIXED = 2.OE6 +'0.14 x DCC + 0.
Page 137 and 138:
(V.B - 6) The discount rate rDIS is
Page 139 and 140:
VI. Program Flow DELSOL can be used
Page 141 and 142:
.I c to reach a certain design powe
Page 143 and 144:
an optimum storage size based on en
Page 145 and 146:
I Read i n N + K t x ] ITAPE=O or 1
Page 147 and 148:
Read in System Definition from File
Page 149 and 150:
. TABLE VI-2 OUTPIJT FROM OPTIMIZAT
Page 151 and 152:
) Power level specific for detailed
Page 153 and 154:
V1.B. DescriDtion of Subroutines in
Page 155 and 156:
C cy c: C C
Page 157 and 158:
L C c C c C c: C c C C c C C C c: c
Page 159:
C C C C C C C C C C C C C C C C C C
Page 162 and 163:
At the time of the comparisons MIRV
Page 164 and 165:
164 n m W * ? 0 F F 0 ( z W /M W) X
Page 166 and 167:
166 Phase 1; CDRL Item 2, Pilot Pla
Page 168 and 169:
168 43. Vittitoe, C. N. et al., “
Page 170 and 171:
since the values have been already
Page 172 and 173:
Title Card $BASIC$ $FIELD$ $HSTAT$
Page 174 and 175:
IPROB NYEAR HRDEL UDAY UTIME NUAZ N
Page 176 and 177:
TDESP PLAT ALT INSOL SOLCON IWEATH
Page 178 and 179:
REFTIM REFSOL ASTART IATM ATMl ATM2
Page 180 and 181:
IHPR NLAND XTOWER YTOWER (I= 1 ,NLA
Page 182 and 183:
182 1 1-40 2 1-10 11-20 3 1-10 11-2
Page 184 and 185:
SIGTX SIGTY ICANT NCANTX NCANTY HCA
Page 186 and 187:
THT TOWL TOWD IREC W H RRECL IAUTOP
Page 188 and 189:
For cavities or flat plates, the fo
Page 190 and 191:
NXFLX FAZMIN FAZMAX 190 AZMF (North
Page 192 and 193:
NMXFLXY) identifies the exact point
Page 194 and 195:
REFPRL FSP FEP EFFSTR PF SMULT IP H
Page 196 and 197:
IHOPT NUMTHT THTST THTEND NUMREC WS
Page 198 and 199:
SMULT IPLFL(1) (I=l ,NUMOPT) IPROPT
Page 200 and 201:
CH CL CWR CWDR CWDA ITHT CTOWl CTOW
Page 202 and 203:
ICHE CHEREF PHEREF XHEP APHREF PPHR
Page 204 and 205:
CKNREF PKNREF XKN 204 Sodium-to-sal
Page 206 and 207:
TC TR FDEBT RDEBT ROE IDEP NDEP NYO
Page 208 and 209:
The next set of cards is read at th
Page 210 and 211:
value of insolation from the specif
Page 212 and 213:
For optimizing cavity depth, the IO
Page 214 and 215:
For the final detailed performance
Page 216 and 217:
from that specified during the opti
Page 218 and 219:
Sample Problem 4 - User Defined Fie
Page 221 and 222:
GLOSSARY Absorber: The portion of t
Page 223 and 224:
Heat Exchanger: A component in whic
Page 225 and 226:
Receiver System Efficiency: The rat
Page 227 and 228:
Variable AEVMAX AEVREF AFDC ALP(I)
Page 229 and 230:
Variable H H20 HCANT HM HPANL HRDEL
Page 231 and 232:
Variable PARL7,PARLI) PARL9,PARLlO
Page 233 and 234:
Variable W WEATH WEND WM WPANL WST
Page 235 and 236:
Analysis Review and Critique 6503 8
Page 237 and 238:
Martin Marietta Aerospace P.O. Box
Page 239:
Stone and Webster Engineering Corpo
show all

A User's Manual for DELSOL3 - prod.sandia.gov - Sandia National ...

Create successful ePaper yourself

Delete template?

Save as template?