blob: 868ef26ebffe559893527f55f968c8e6fad47dfc (
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
|
*-------------------------------------------------------------------------------------------------
* Name : TEST_GEO_latt_tspc_S30.x2m
* Author : G. Cirillo, D. Manzione
* HEXAGONAL LATTICE MODEL (7 HEXAGONS) IN HEXAGONAL BOX TEST SURFACE GEO
* TYPGEO 12 NBFOLD 0
*-------------------------------------------------------------------------------------------------
* Define STRUCTURES and MODULES used
*-------------------------------------------------------------------------------------------------
SEQ_ASCII ALFR_LATT_b :: FILE './lattice_hex_tspc_S30_7cells.dat' ;
SEQ_ASCII LatTEST_b :: FILE './LatTEST_b_S30.ps' ;
SEQ_ASCII Lines_HEX.m :: FILE './Lines_HEX.m' ;
LINKED_LIST
LATT MICRO MICRO2 TRACK CP FLUX OUT ;
MODULE
G2S: LIB: SALT: MCCGT: TLM: TONE: ASM: FLU: EDI: DELETE: END: ;
SEQ_BINARY
TRACK_f ;
PROCEDURE assertS ;
*
*-------------------------------------------------------------------------------------------------
* Input Data
*-------------------------------------------------------------------------------------------------
REAL temp_fuel := 1200.0 ; (*in K*)
REAL temp_fuel2 := 900.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 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 XMAS 172 energy discr
*-------------------------------------------------------------------------------------------------
MICRO := LIB: ::
EDIT 2
NMIX 8 (*MAXIMUM OF MATERIAL MIXTURES*)
CTRA NONE
ANIS 2 (*LINEAR ANISOTROPY SCATTERING*)
ADED 2 NFTOT NINEL
MIXS LIB: DRAGON FIL: DLIB_8R0_315
MIX 1 <<temp_clad>> <<densHe>> GAS
He4 = He4 1.0
MIX 2 <<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 3 COMB 1 1.0
MIX 4 <<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 5 <<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 + Self-shielding:
*-------------------------------------------------------------------------------------------------
*-----------------------------------------------------
* Tracking of the imported geometry
*-----------------------------------------------------
LatTEST_b := G2S: ALFR_LATT_b ::
EDIT 2
DRAWMIX ;
TRACK TRACK_f := SALT: ALFR_LATT_b ::
TITL 'TEST_GEO imported'
EDIT 2
ANIS 1
TSPC 12 50.0 LONG
MERGMIX
;
Lines_HEX.m := TLM: TRACK TRACK_f ::
EDIT 99
NTPO 1
DIRECTIONS NoPause DIR 0 ;
*-----------------------------------------------------
* Resonance self-shielding calculation
*-----------------------------------------------------
MICRO2 := TONE: MICRO TRACK TRACK_f :: EDIT 1 MXIT 3 TRAN PIJ ;
TRACK TRACK_f := DELETE: TRACK TRACK_f ;
*-----------------------------------------------------
TRACK TRACK_f := SALT: ALFR_LATT_b ::
TITL 'TEST_GEO imported'
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
;
*-------------------------------------------------------------------------------------------------
* Flux calculation
*-------------------------------------------------------------------------------------------------
CP := ASM: MICRO2 TRACK TRACK_f :: ARM ;
FLUX := FLU: CP MICRO2 TRACK TRACK_f :: TYPE K ;
*-------------------------------------------------------------------------------------------------
* Editing
*-------------------------------------------------------------------------------------------------
OUT := EDI: MICRO2 TRACK FLUX ALFR_LATT_b CP ::
EDIT 4
MERG COMP
COND
;
assertS FLUX :: 'K-EFFECTIVE' 1 0.1612311 ;
ECHO "test TEST_GEO_latt_tspc_S30 completed" ;
END: ;
QUIT "LIST" .
|
