blob: 60c56adfcf180d0461d036161b1b9b6c51833fed (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
|
*-------------------------------------------------------------------------------------------------
* Name : TEST_GEO_hex_sect_tspc.x2m
* Author : G. Cirillo
* HEXAGONAL PIN MODEL TO TEST GLOW GEO
*-------------------------------------------------------------------------------------------------
* Define STRUCTURES and MODULES used
*-------------------------------------------------------------------------------------------------
SEQ_ASCII TDT_LAT :: FILE './ALFRED_1cell.dat' ;
SEQ_ASCII TDT_LAT_SECT :: FILE './ALFRED_1cell_sect.dat' ;
SEQ_ASCII HexTEST.ps :: FILE './HexTEST.ps' ;
SEQ_ASCII HexTEST2.ps :: FILE './HexTEST2.ps' ;
SEQ_ASCII Lines_HEX.m :: FILE './Lines_HEX.m' ;
LINKED_LIST
MICRO MICRO2 TRACK TRACK2 CP FLUX OUT ;
MODULE
GEO: G2S: LIB: SALT: MCCGT: TLM: TONE: ASM: FLU: EDI: END: ;
SEQ_BINARY
TRACK_f TRACK_f2 ;
PROCEDURE assertS ;
*
*-------------------------------------------------------------------------------------------------
* Input Data
*-------------------------------------------------------------------------------------------------
REAL temp_fuel := 1200.0 ; (*in K*)
REAL temp_clad := 600.0 ; (*in K*)
REAL temp_cool := 600.0 ; (*in K*)
REAL densUO2 := 10.3235 ; (*in g/cm^3*)
REAL densPuO2 := 10.7771 ; (*in g/cm^3*)
(*REAL densMOX := 10.41561 ; in g/cm^3*)
REAL densAIM1 := 7.972 ; (*in g/cm^3*)
REAL densLEAD := 10.503 ; (*in g/cm^3*)
REAL densHe := 0.0001663 ; (*in g/cm^3*)
*-------------------------------------------------------------------------------------------------
* Microlib: DRAGLIB ENDFB_VII 315 energy discr
*-------------------------------------------------------------------------------------------------
MICRO := LIB: ::
EDIT 2
NMIX 8 (*MAXIMUM OF MATERIAL MIXTURES*)
CTRA NONE
ANIS 2 (*LINEAR ANISOTROPY SCATTERING*)
MIXS LIB: DRAGON FIL: DLIB_315
MIX 1 <<temp_fuel>> <<densUO2>>
U235 = U235 1.762963E-03 1 O16 = O16 1.18235352E-01 1
U238 = U238 8.79681435E-01 1 (*O17 = O17 4.77255E-05*)
U234 = U234 5.81718E-06 1 (*O18 = O18 2.66717E-04*)
O17 = O17 3.144425E-04
MIX 2 <<temp_fuel>> <<densPuO2>>
Pu238 = Pu238 2.030129548E-02 1 Pu241 = Pu241 5.313709197E-02 1
Pu239 = Pu239 4.951133413E-01 1 Pu242 = Pu242 6.696932147E-02 1
Pu240 = Pu240 2.3501943E-01 1 Am241 = Am241 1.16169418E-02 1
O16 = O16 1.175300115E-01 1
O17 = O17 3.1256696E-04
(*O17 = O17 4.74407E-05*)
(*O18 = O18 2.6512626E-04*)
MIX 3 COMB 1 0.80204511 2 0.19795489
MIX 4 COMB 3 1.0
MIX 5 COMB 3 1.0
MIX 6 <<temp_clad>> <<densHe>> GAS
He4 = He4 1.0
MIX 7 <<temp_clad>> <<densAIM1>>
B10 = B10 1.1013297525043955E-05 2 Ni58 = Ni58 0.10415644320971433 2
B11 = B11 4.898670247495604E-05 2 Ni60 = Ni60 0.04150274249380582 2
C0 = C0 9E-04 2
(*C12 = C12 0.0008892061605771179*) Ni61 = Ni61 0.0018342065822099006 2
(*C13 = C13 1.0793839422882127E-05*) Ni62 = Ni62 0.005943974570405113 2
N14 = N14 0.0001494115819963943 2 Ni64 = Ni64 0.0015626331438648383 2
N15 = N15 5.884180036056514E-07 2 Cu63 = Cu63 0.00030815637842357554
O16 = O16 9.973468728094326E-05 2 Cu65 = Cu65 0.00014184362157642448
O17 = O17 2.653127191E-07
(*O17 = O17 4.02684437425691E-08*) Zr90 = Zr90 0.00015211835252769947 2
(*O18 = O18 2.2504427531418012E-07*) Zr91 = Zr91 3.3542662359302346E-05 2
Al27 = Al27 0.00015 2 Zr92 = Zr92 5.183430439088447E-05 2
Si28 = Si28 0.007809248969041125 2 Zr94 = Zr94 5.367331302703024E-05 2
Si29 = Si29 0.0004106987748868437 2 Zr96 = Zr96 8.831367695083468E-06 2
Si30 = Si30 0.0002800522560720316 2 Nb93 = Nb93 0.00015 2
P31 = P31 0.00045 Mo92 = Mo92 0.0021047837070591823 2
S32 = S32 9.476861080827586E-05 Mo94 = Mo94 0.001348697757094504 2
S33 = S33 7.698819377540353E-07 Mo95 = Mo95 0.0023550696583545722 2
S34 = S34 4.445152400876023E-06 Mo96 = Mo96 0.00249980045270409 2
S36 = S36 1.6354853094087033E-08 Mo97 = Mo97 0.0014516393261295652 2
Ca40 = Ca40 0.00019332358625973815 Mo98 = Mo98 0.00371810585345211 2
Ca42 = Ca42 1.3547188215775078E-06 2 Mo100 = Mo100 0.001521903245205976 2
Ca43 = Ca43 2.8940716633277043E-07 2 Sn112 = Sn112 4.571964310655134E-07
Ca44 = Ca44 4.575631664190023E-06 2 Sn114 = Sn114 3.1663618573474405E-07
Ca46 = Ca46 9.172845337918173E-09 Sn115 = Sn115 1.6454847928465057E-07
Ca48 = Ca48 4.474832428236633E-07 2 Sn116 = Sn116 7.098010841901756E-06
Ti46 = Ti46 0.0003168038086310422 2 Sn117 = Sn117 3.7815425451416758E-06
Ti47 = Ti47 0.00029191129292780875 2 Sn118 = Sn118 1.2027521650336542E-05
Ti48 = Ti48 0.0029538020902329385 2 Sn119 = Sn119 4.3019900178835224E-06
Ti49 = Ti49 0.00022128762152676788 2 Sn120 = Sn120 1.6453584513264257E-05
Ti50 = Ti50 0.00021619518668144252 2 Sn122 = Sn122 2.377274090454982E-06
V50 = V50 7.353610074676385E-07 Sn124 = Sn124 3.021695244932352E-06
V51 = V51 0.0002992646389925323 Sb121 = Sb121 2.840390648290279E-05
Cr50 = Cr50 0.006051845602941554 2 Sb123 = Sb123 2.159609351709721E-05
Cr52 = Cr52 0.12136407540945154 2 (*Ta180 = Ta180 1.7915400047357467E-08*)
Cr53 = Cr53 0.014026700887079016 2 (*Ta181 = Ta181 0.00014998208459995264*)
(*Ta181 = Ta181 1.5E-04*)
Cr54 = Cr54 0.00355737810052789 2 (*W180 = W180 3.5237265941784446E-07*)
Mn55 = Mn55 0.015 2 (*W182 = W182 7.868114842886196E-05*)
W182 = W182 7.903352109E-05
Fe54 = Fe54 0.03688469009223297 W183 = W183 4.2721807585340764E-05
Fe56 = Fe56 0.6004294481283934 2 W184 = W184 9.19745762621493E-05
Fe57 = Fe57 0.01411455230058187 W186 = W186 8.627009506423013E-05
Fe58 = Fe58 0.0019113094787918548 2
Co59 = Co59 0.0003
MIX 8 <<temp_cool>> <<densLEAD>>
Na23 = Na23 1E-05 As75 = As75 5E-06
Mg24 = Mg24 7.7909641699443E-06 Ag107 = Ag107 5.137617589720804E-06
Mg25 = Mg25 1.0300414880798231E-06 Ag109 = Ag109 4.862382410279197E-06
Mg26 = Mg26 1.1789943419758772E-06 Sn112 = Sn112 4.571964310655134E-08
Ca40 = Ca40 9.666179312986907E-08 Sn114 = Sn114 3.166361857347441E-08
Ca42 = Ca42 6.773594107887539E-08 3 Sn115 = Sn115 1.645484792846506E-08
Ca43 = Ca43 1.4470358316638523E-08 3 Sn116 = Sn116 7.098010841901756E-07
Ca44 = Ca44 2.2878158320950115E-07 3 Sn117 = Sn117 3.781542545141676E-07
Ca46 = Ca46 4.586422668959087E-10 Sn118 = Sn118 1.2027521650336543E-06
Ca48 = Ca48 2.2374162141183165E-08 3 Sn119 = Sn119 4.3019900178835226E-07
Fe54 = Fe54 5.645558222706856E-07 Sn120 = Sn120 1.6453584513264259E-06
Fe56 = Fe56 9.190152877956246E-06 3 Sn122 = Sn122 2.3772740904549819E-07
Fe57 = Fe57 2.1603686136746366E-07 Sn124 = Sn124 3.021695244932352E-07
Fe58 = Fe58 2.9254438405605885E-08 3 Sb121 = Sb121 5.680781296580559E-06
Cu63 = Cu63 6.8479195205239015E-06 Sb123 = Sb123 4.319218703419443E-06
Cu65 = Cu65 3.1520804794761E-06 Pb204 = Pb204 0.01377877075349251 3
(*Zn64 = Zn64 4.8080491426635065E-06*) Pb206 = Pb206 0.2395190658912157 3
Zn66 = Zn66 2.796254092194871E-06 Pb207 = Pb207 0.22070984690476395 3
Zn67 = Zn67 4.135741306808834E-07 Pb208 = Pb208 0.5258423164505278 3
Zn68 = Zn68 1.916879606631653E-06 Bi209 = Bi209 6E-05
Zn70 = Zn70 6.524302782908711E-08
;
*-------------------------------------------------------------------------------------------------
* Tracking + PostScript
*-------------------------------------------------------------------------------------------------
*---------------------------------------------------------
* PostScript
*---------------------------------------------------------
HexTEST.ps := G2S: TDT_LAT_SECT ::
EDIT 2
DRAWMIX ;
HexTEST2.ps := G2S: TDT_LAT ::
EDIT 2
DRAWMIX ;
*---------------------------------------------------------
* Tracking sectorized geometry for self-shielding calc
*---------------------------------------------------------
TRACK2 TRACK_f2 := SALT: TDT_LAT ::
TITL 'Tracking for self-shielding'
EDIT 2
ANIS 1
TSPC 12 50.0
;
*-------------------------------------------------------------------------------------------------
* Resonance self-shielding calculation
*-------------------------------------------------------------------------------------------------
MICRO2 := TONE: MICRO TRACK2 TRACK_f2 :: EDIT 1 MXIT 3 TRAN PIJ ;
*---------------------------------------------------------
* Tracking sectorized geometry for the flux calc
*---------------------------------------------------------
TRACK TRACK_f := SALT: TDT_LAT_SECT ::
TITL 'Tracking for flux calculation'
EDIT 2
ANIS 2
TSPC 12 50.0 LONG
;
TRACK := MCCGT: TRACK TRACK_f ::
EDIT 2
GAUS 4 AAC 80 TMT STIS 1 EPSI 1E-5 MCU 10000
;
Lines_HEX.m := TLM: TRACK TRACK_f ::
EDIT 99
NTPO 1
DIRECTIONS NoPause DIR 0 ;
*-------------------------------------------------------------------------------------------------
* Flux calculation
*-------------------------------------------------------------------------------------------------
CP := ASM: MICRO2 TRACK TRACK_f :: ARM ;
FLUX := FLU: CP MICRO2 TRACK TRACK_f :: TYPE K P1 ;
*-------------------------------------------------------------------------------------------------
* Editing
*-------------------------------------------------------------------------------------------------
OUT := EDI: MICRO2 TRACK FLUX TDT_LAT_SECT CP ::
EDIT 4
MERG COMP
COND
;
assertS FLUX :: 'K-EFFECTIVE' 1 1.341990 ;
ECHO "test TEST_GEO_hex_sect_tspc completed" ;
END: ;
QUIT "LIST" .
|
