blob: edd342591c1981311801e2658808357eb9cb43ee (
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
|
*DECK Mix_UOX
*----
* Name : Mix_UOX.c2m
* Type : DRAGON procedure
* Use : Microlib generation with Draglibs for
* UOX calculation with 4 fuel regions
* Author : A. Hebert
*
* Procedure called as:
*
* LIBRARY := Mix_UOX ::
* <<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 0.6022 * 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 dens_mod pbore * * ;
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 ;
*isot 4.35% 83.79% 9.50% 2.36%
REAL P_Sn115 P_Sn116 P_Sn117
P_Sn118 P_Sn119 P_Sn120 P_Sn122 P_Sn124 :=
0.33E-02 14.21E-02 7.56E-02
24.06E-02 8.59E-02 32.92E-02 4.75E-02 6.04E-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 M_Sn := P_Sn115 114.90335 *
P_Sn116 115.901748 * P_Sn117 116.902961 *
P_Sn118 117.901613 * P_Sn119 118.903316 *
P_Sn120 119.902207 * P_Sn122 121.903451 *
P_Sn124 123.905283 * + + + + + + + ;
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 / * ;
REAL N_Sn_clad := P_Sn rho_clad A * M_Sn / * ;
REAL N_Sn115 N_Sn116
N_Sn117 N_Sn118 N_Sn119 N_Sn120
N_Sn122 N_Sn124 :=
N_Sn_clad P_Sn115 *
N_Sn_clad P_Sn116 * N_Sn_clad P_Sn117 * N_Sn_clad P_Sn118 *
N_Sn_clad P_Sn119 * N_Sn_clad P_Sn120 * N_Sn_clad P_Sn122 *
N_Sn_clad P_Sn124 * ;
* ;
LIBRARY := LIB: ::
EDIT 1
DEPL LIB: DRAGON FIL: <<Library>>
NMIX 7 CTRA APOL
SUBG
MIXS LIB: DRAGON FIL: <<Library>>
*----
* Cladding fuel
*----
MIX 1 <<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>>
Sn115f = Sn115 <<N_Sn115>>
Sn116f = Sn116 <<N_Sn116>>
Sn117f = Sn117 <<N_Sn117>>
Sn118f = Sn118 <<N_Sn118>>
Sn119f = Sn119 <<N_Sn119>>
Sn120f = Sn120 <<N_Sn120>>
Sn122f = Sn122 <<N_Sn122>>
Sn124f = Sn124 <<N_Sn124>>
*----
* Gap
*----
MIX 2 <<TgainK>> NOEV
Al27 = Al27 1.00E-08
*----
* Moderator
*----
MIX 3 <<TmodeK>> NOEV
H1_H2O = H1_H2O <<N_2H2O>>
O16_H2O = O16 <<N_H2O>>
B10m = B10 <<N_B10>>
B11m = B11 <<N_B11>>
*----
* Fuel
*----
MIX 4 <<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 5
COMB 4 1.0
MIX 6
COMB 4 1.0
MIX 7
COMB 4 1.0
;
END: ;
QUIT .
|
