blob: 39755c3f35df5924070dbfe26fdd8a975190c468 (
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
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
|
*----
* Name : Mix_UOX_32.c2m
* Type : DRAGON procedure
* Use : Microlib generation with Draglibs for
* UOX calculation with 32 fuel regions
* Author : R. Vallerent
*
* Procedure called as:
*
* LIBRARY := Mix_UOX_32 ::
* <<dens_mod>> <<pbore>> <<N_U5_UOX>> <<densU_UOX>>
* <<temp_comb>> <<temp_mode>> ;
*
* Input data :
* dens_mod : Theoretical moderator density
* pbore : Natural boron concentration (ppm)
* N_U5_UOX : U235 enrichment in UOX-Fuel
* densU_UOX : Theoretical density of Uranium in UOX-Fuel
* temp_comb : Theoretical Fuel temperature (K)
* temp_mode : Theoretical moderator temperature (K)
*
* Output data :
* LIBRARY : Microlib with table of probabilities
PARAMETER LIBRARY ::
EDIT 0
::: LINKED_LIST LIBRARY ;
;
*----
* Input data recovery
*----
STRING Library ;
:: >>Library<< ;
REAL dens_mod pbore N_U5_UOX densU_UOX ;
:: >>dens_mod<< >>pbore<< >>N_U5_UOX<< >>densU_UOX<< ;
REAL temp_comb temp_mode ;
:: >>temp_comb<< >>temp_mode<< ;
*----
* Modules used in this procedure
*----
MODULE LIB: UTL: DELETE: END: ;
REAL TfuelK := temp_comb ;
REAL TmodeK := temp_mode ;
REAL TgainK := TmodeK ;
* Calculation of moderator data
* NOMBRE D AVOGADRO
REAL A := 6.022094E-1 ; ! Normalisation des concentrations par 10E-24
REAL M_H2O := 15.9994 2.0 1.00794 * + ;
REAL N_H2O := dens_mod A * M_H2O / ;
REAL b10 b11 := 0.199 0.801 ; ! Bore naturel
REAL coeff := 0.199 A 1.0E-06 * * 10.8110164 / ;
REAL N_B10 := coeff 0.659 pbore * * ; ! densite nominale
REAL N_2H2O N_B11 := 2.0 N_H2O * N_B10 b11 * b10 / ;
* Calculation of fuel data
REAL den_uo2 := densU_UOX ;
REAL M_U4 M_U5 M_U8 M_U6 :=
234.040952088 235.043929918 238.050788247 236.045568006 ;
REAL M_Pu8 M_Pu9 M_Pu0 M_Pu1 :=
238.049559894 239.052163381 240.053813545 241.056851456 ;
REAL M_Pu2 M_Am1 M_O16 :=
242.058742611 241.056829144 15.99491461956 ;
REAL N_U5 := N_U5_UOX ;
REAL N_U4 := N_U5 54.0E-04 * 0.72 / ;
REAL N_U6 N_Pu8 N_Pu9 N_Pu0 N_Pu1 N_Pu2 N_Am1 :=
1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 1.0E-14 ;
REAL N_U8 := 1.0 N_U5 - N_U4 - N_Pu8 7.0 * - ;
REAL N_O16F := N_U4 M_U4 / N_U5 M_U5 / + N_U8 M_U8 / + N_Pu8 M_Pu8 / + ;
EVALUATE N_O16F := N_O16F N_U6 M_U6 / + ;
EVALUATE N_O16F := N_O16F N_Pu9 M_Pu9 / + N_Pu0 M_Pu0 / + ;
EVALUATE N_O16F := N_O16F N_Pu1 M_Pu1 / + N_Pu2 M_Pu2 / + ;
EVALUATE N_O16F := N_O16F N_Am1 M_Am1 / + 2.0 * M_O16 * ;
REAL N_TOTF := N_U4 N_U5 + N_U6 + N_U8 + N_Pu8 + N_Pu9 + N_Pu0 + ;
EVALUATE N_TOTF := N_TOTF N_Pu1 + N_Pu2 + N_Am1 + N_O16F + ;
EVALUATE N_U4 := N_U4 N_TOTF / ;
EVALUATE N_U5 := N_U5 N_TOTF / ;
EVALUATE N_U6 := N_U6 N_TOTF / ;
EVALUATE N_U8 := N_U8 N_TOTF / ;
EVALUATE N_Pu8 := N_Pu8 N_TOTF / ;
EVALUATE N_Pu9 := N_Pu9 N_TOTF / ;
EVALUATE N_Pu0 := N_Pu0 N_TOTF / ;
EVALUATE N_Pu1 := N_Pu1 N_TOTF / ;
EVALUATE N_Pu2 := N_Pu2 N_TOTF / ;
EVALUATE N_Am1 := N_Am1 N_TOTF / ;
EVALUATE N_O16F := N_O16F N_TOTF / ;
* --------------------------------
* Cladding calculation
* --------------------------------
REAL P_Zr P_Fe P_Cr P_O16 P_Sn :=
98.115E-02 0.210E-02 0.100E-02 0.125E-02 1.450E-02 ;
REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 :=
51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ;
REAL P_Fe54 P_Fe56 P_Fe57 P_Fe58 :=
5.602E-02 91.865E-02 2.243E-02 0.290E-02 ;
*isot 5.8 % 91.72% 2.20% 0.28%
REAL P_Cr50 P_Cr52 P_Cr53 P_Cr54 :=
4.178E-02 83.701E-02 9.673E-02 2.448E-02 ;
REAL M_Zr := P_Zr90 89.904711 * P_Zr91 90.905643 * P_Zr92 91.905039 *
P_Zr94 93.906320 * P_Zr96 95.90829 * + + + + ;
REAL M_Fe := P_Fe54 53.9396148 * P_Fe56 55.9349421 *
P_Fe57 56.9353987 * P_Fe58 57.9332805 * + + + ;
REAL M_Cr := P_Cr50 49.9460496 * P_Cr52 51.9405119 *
P_Cr53 52.9406538 * P_Cr54 53.9388849 * + + + ;
REAL rho_clad := 6.55 ;
REAL N_Zr_clad := P_Zr rho_clad A * M_Zr / * ;
REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 :=
N_Zr_clad P_Zr90 * N_Zr_clad P_Zr91 * N_Zr_clad P_Zr92 *
N_Zr_clad P_Zr94 * N_Zr_clad P_Zr96 * ;
REAL N_Fe_clad := P_Fe rho_clad A * M_Fe / * ;
REAL N_Fe54 N_Fe56 N_Fe57 N_Fe58 :=
N_Fe_clad P_Fe54 * N_Fe_clad P_Fe56 * N_Fe_clad P_Fe57 *
N_Fe_clad P_Fe58 * ;
REAL N_Cr_clad := P_Cr rho_clad A * M_Cr / * ;
REAL N_Cr50 N_Cr52 N_Cr53 N_Cr54 :=
N_Cr_clad P_Cr50 * N_Cr_clad P_Cr52 * N_Cr_clad P_Cr53 *
N_Cr_clad P_Cr54 * ;
REAL N_O16G := P_O16 rho_clad A * 15.9949150 / * ;
* ;
LIBRARY := LIB: ::
EDIT 0
DEPL LIB: DRAGON FIL: <<Library>>
NMIX 40 CTRA APOL
SUBG
MIXS LIB: DRAGON FIL: <<Library>>
*----
* Cladding TG et TI
*----
MIX 1 <<TgainK>> NOEV
Zr90g = Zr90 <<N_Zr90>>
Zr91g = Zr91 <<N_Zr91>>
Zr92g = Zr92 <<N_Zr92>>
Zr94g = Zr94 <<N_Zr94>>
Zr96g = Zr96 <<N_Zr96>>
Fe54g = Fe54 <<N_Fe54>>
Fe56g = Fe56 <<N_Fe56>>
Fe57g = Fe57 <<N_Fe57>>
Fe58g = Fe58 <<N_Fe58>>
Cr50g = Cr50 <<N_Cr50>>
Cr52g = Cr52 <<N_Cr52>>
Cr53g = Cr53 <<N_Cr53>>
Cr54g = Cr54 <<N_Cr54>>
O16Gg = O16 <<N_O16G>>
*----
* Cladding fuel
*----
MIX 2 <<TgainK>> NOEV
Zr90f = Zr90 <<N_Zr90>>
Zr91f = Zr91 <<N_Zr91>>
Zr92f = Zr92 <<N_Zr92>>
Zr94f = Zr94 <<N_Zr94>>
Zr96f = Zr96 <<N_Zr96>>
Fe54f = Fe54 <<N_Fe54>>
Fe56f = Fe56 <<N_Fe56>>
Fe57f = Fe57 <<N_Fe57>>
Fe58f = Fe58 <<N_Fe58>>
Cr50f = Cr50 <<N_Cr50>>
Cr52f = Cr52 <<N_Cr52>>
Cr53f = Cr53 <<N_Cr53>>
Cr54f = Cr54 <<N_Cr54>>
O16Gf = O16 <<N_O16G>>
*----
* Gap
*----
MIX 3 <<TgainK>> NOEV
Al27 = Al27 1.00E-08
*----
* Moderator
*----
MIX 4 <<TmodeK>> NOEV ! pastille d'U
H1_H2O = H1_H2O <<N_2H2O>>
O16_H2O = O16 <<N_H2O>>
B10 = B10 <<N_B10>>
B11 = B11 <<N_B11>>
MIX 5
COMB 4 1.0 ! pastille d'U (2)
MIX 6
COMB 4 1.0 ! gaines T2
MIX 7
COMB 4 1.0 ! gaines T1
MIX 8
COMB 4 1.0 ! gaine TI
*----
* Fuel -> 8 * 4 evolving mixes
*----
* C1 Cell
MIX 9 <<TfuelK>> <<den_uo2>>
U234 = U234 <<N_U4>>
U235 = U235 <<N_U5>> 1
U236 = U236 <<N_U6>> 1
U238 = U238 <<N_U8>> 1
Pu238 = Pu238 <<N_Pu8>>
Pu239 = Pu239 <<N_Pu9>> 1
Pu240 = Pu240 <<N_Pu0>> 1
Pu241 = Pu241 <<N_Pu1>> 1
Pu242 = Pu242 <<N_Pu2>> 1
Am241 = Am241 <<N_Am1>>
O16F = O16 <<N_O16F>>
MIX 10
COMB 9 1.0
MIX 11
COMB 9 1.0
MIX 12
COMB 9 1.0
* C2 Cell
MIX 13
COMB 9 1.0
MIX 14
COMB 9 1.0
MIX 15
COMB 9 1.0
MIX 16
COMB 9 1.0
* C3 Cell
MIX 17
COMB 9 1.0
MIX 18
COMB 9 1.0
MIX 19
COMB 9 1.0
MIX 20
COMB 9 1.0
* C4 Cell
MIX 21
COMB 9 1.0
MIX 22
COMB 9 1.0
MIX 23
COMB 9 1.0
MIX 24
COMB 9 1.0
* C5 Cell
MIX 25
COMB 9 1.0
MIX 26
COMB 9 1.0
MIX 27
COMB 9 1.0
MIX 28
COMB 9 1.0
* C6 Cell
MIX 29
COMB 9 1.0
MIX 30
COMB 9 1.0
MIX 31
COMB 9 1.0
MIX 32
COMB 9 1.0
* C7 Cell
MIX 33
COMB 9 1.0
MIX 34
COMB 9 1.0
MIX 35
COMB 9 1.0
MIX 36
COMB 9 1.0
* C8 Cell
MIX 37
COMB 9 1.0
MIX 38
COMB 9 1.0
MIX 39
COMB 9 1.0
MIX 40
COMB 9 1.0
;
END: ;
QUIT .
|
