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

More documents

Recommendations

Info

The next set of cards is read at the start of the optimization calculation. At this point in the calculation, the design will be affected by flux limits and aim strategies. IAUTOP=l on the $REC$ card specifies the 1-dimensional smart aim- ing option necessary for a flux limited external receiver design. In the $OPT$ namelist, tower height and receiver dimensions are varied over a range compatible with the range of power levels. The variation in receiver height to width is specified to allow spreading of the flux along the height. IFLX=l in $NLFLUX$ turns on the flux calculation, and the flux limit is set by FLXLIM. Since the maximum flux should occur at the center of the north panel for the default design point, only that one point is tested for flux levels. The default values in $NLFLUX$ will result in the flux being calculated at the correct point and compared to the limit FLXLIM. The absence of storage is indicated by setting SMULT=1.0 in both $OPT$ and $NLEFF$. It is not necessary to set the reference costs to zero because the code will calculate zero storage costs if no energy is available to charge storage. The optimum systems for all power levels specified by the non-zero values of IPLFL will be saved on a file on Unit 30 because the value of IOTAPE=l is specified. - Comments on Output The first part of the output consists of the printing of the performance namelists, followed by summaries of the heliostat design and receiver used in the initial performance calculation, zone by zone density and heliostat counts, an insolation table, and the zone by zone yearly average performance. A printing of the optimization namelists follows, then a summary of the optimization variables and design constraints. The default choice of IPROPT in $OPT$ produces the abbreviated list indicating the combinations of optimization variables searched, along with pertinent information related to each set considered. The optimization results are then presented in a series of tables. Note in the system design summary that all design variables fall within the ranges searched except for the receiver height in the three smallest power sizes, which is at the maximum allowed height to width ratio. Although energy costs would probably not change significantly, a thorough user would want to rerun this case to allow consideration of larger height to width ratios. Note also in the field layout summary that the field extends to RADMAX in all cases for a number of azimuthal zones. The user could rerun this case with a larger value for RADMAX to allow the code more flexibility in zone selection for field build-up. In no case was the flux limit met for this series of calculations. It is likely that a finer grid of choices for receiver dimensions and tower height would have resulted in optimum systems in which the peak flux was much closer to the allowed flux on the receiver. Sample Problem lb looks at the flux distribution over the receiver for one of these systems in which the flux limit was not met. 208
Sample Problem lb - Performance/Flux Run for an Optimized Design Problem Statement For the 125 MWe system optimized in Sample Problem la, produce a detailed performance summary of the system on an annual basis, and show that the flux at any point on the receiver does not exceed the flux limit of 0.85 MW/m2. InDut Cards SAMPLE PROBLEM Ib $BASIC IPROB=O. ITAPE=3. TDESP-125.0 $ $FIELD B BHSTAT $ OHEC 3 BNLFLUX IFtX-1, IFXOUT(3. 1)=1, IFLAUT=1. DIAMF.12.0, NXFLX=13. NYFLX=13. FAZMIN=O.O. FAZMAX=360.0. FZMIN=12.0. FZMAX=-12.0 $ BNLEFF $ BREC W=-I. 8 Analysis of InDut For an annual performance calculation, IPROB=O is specified. ITAPE=3 indicates that the optimized system is stored on a file on Unit 30, and TDESP specifies the power level of the system to be so analyzed. All values relating to the heliostats, field, and receiver will be defined from the file on Unit 30, and should not be redefined in these namelist inputs, even though the namelist cards must still remain in the input in the correct order. A flux map over the whole receiver is specified by the appropriate choices of the values in Namelist $NLFLUX$. The values of FZMIN, FZMAX, FAZMIN, and FAZMAX are distances from the point specified by the values of XFC, YFC, and ZFC, which in this case were the default values. IFXOUT(3,1)=1 specifies that the flux map will be for the third of NYEAR days (in this case, day 81) at the first time for which calculations are done, that is, solar noon (time=O.O). The value of DIAMF would have defaulted to the receiver diameter anyway, but was included as a precautionary/informational action. Comments on OutDut As for the first part of Problem la, printouts of namelists and then summaries of heliostats and the receiver are printed. Positions of flux points are next printed, followed by field and insolation summaries. Next, yearly average and design point performance summaries are printed, followed by an annual power production table and summary. Finally, a table of fluxes at the requested locations is printed. Note that the insolation value used in the design point summary and calculation of fluxes is different from the value used in the annual power production table. The former uses the value REFSOL, while the latter uses the calculated 209
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 34 and 35:
34 N' Figure 11-6. Code Defined Fie
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 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 ...

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?