blob: 1e1af235f9ece81a9e177648698d01ef78e06481 (
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
|
* DRAGON CROSS-SECTION LIBRARY PROCEDURE DECK
***********************************************************************
* Purpose:
* define properties values of the mixtures.
*
* author:
* MAJID FASSI FEHRI 09/08/31
* CNSC
*
* Modified:
* Richard CHAMBON 11/02/24
* Ecole Polytechnique
*
***********************************************************************
PARAMETER LIBR ::
::: LINKED_LIST LIBR ; ;
MODULE END: INFO: LIB: ;
INTEGER CALCU ;
STRING LIBNAM := "E6MLIB" ;
PROCEDURE PERTRB ;
REAL TFuel TCool TMode DCool DMode
PCool PMode Bor Xe Sm Np Gad XeRef SmRef NpRef ;
REAL FWGTU5 FWGTU8 FWGTO16 CWGTH1 CWGTD2 CWGTO16
MWGTH1 MWGTD2 MWGTO16 ;
REAL PTtemp PTdens Gaptemp Gapdens CTtemp CTdens
Shtemp Shdens Fdens ;
REAL BWGT0 BWGT1 GWGT5 GWGT7 ;
REAL XWGT5 SWGT9 NWGT9 ;
ECHO "LIBRY PROCEDURE" ;
* Get current calculation value
:: >>CALCU<< ;
* Get current properties values
PERTRB :: <<CALCU>>
>>TFuel<< >>TCool<< >>TMode<< >>DCool<< >>DMode<<
>>PCool<< >>PMode<< >>Bor<< >>Xe<< >>Sm<< >>Np<<
>>Gad<< ;
* Set non-changing properties values
EVALUATE Fdens := 10.39475 ;
EVALUATE Shdens := 6.3918 ;
EVALUATE PTdens := 6.57 ;
EVALUATE CTdens := 6.44 ;
EVALUATE Gapdens := 0.00118 ;
EVALUATE Shtemp := TCool 5.0 + ;
EVALUATE PTtemp := TCool 2.0 - ;
EVALUATE CTtemp := TMode 2.0 + ;
EVALUATE Gaptemp := PTtemp CTtemp + 2.0 / ;
* Change units fron fraction to percent
EVALUATE BWGT0 := 0.199 ;
! EVALUATE BWGT0 := BWGT0 0.000109 * Bor * ;
EVALUATE BWGT0 := BWGT0 Bor * 100.0 * ;
EVALUATE BWGT1 := 0.801 ;
! EVALUATE BWGT1 := BWGT1 0.000109 * Bor * ;
EVALUATE BWGT1 := BWGT1 Bor * 100.0 * ;
EVALUATE GWGT5 := 0.482837 ;
! EVALUATE GWGT5 := GWGT5 0.000033 * Gad * ;
EVALUATE GWGT5 := GWGT5 Gad * 100.0 * ;
EVALUATE GWGT7 := 0.517163 ;
! EVALUATE GWGT7 := GWGT7 0.000033 * Gad * ;
EVALUATE GWGT7 := GWGT7 Gad * 100.0 * ;
* Arbitrary coefficients
EVALUATE XWGT5 := Xe 0.00001 * ;
EVALUATE SWGT9 := Sm 0.00001 * ;
EVALUATE NWGT9 := Np 0.00001 * ;
*----
* FUEL WGT%
*----
INFO: :: EDIT 0
TMP: <<TFuel>> K
ENR: 0.711 WGT%
LIB: WIMSAECL FIL: <<LIBNAM>>
ISO: 3 U235 U238 O16
CALC WGT% UO2 >>FWGTU5<< >>FWGTU8<< >>FWGTO16<< ;
*----
* Coolant WGT%
*----
INFO: :: EDIT 0
TMP: <<TCool>> K
PUR: <<PCool>> WGT%
LIB: WIMSAECL FIL: <<LIBNAM>>
ISO: 3 H1H2O DD2O O16
CALC WGT% D2O >>CWGTH1<< >>CWGTD2<< >>CWGTO16<< ;
*----
* Moderator WGT%
*----
INFO: :: EDIT 0
TMP: <<TMode>> K
PUR: <<PMode>> WGT%
LIB: WIMSAECL FIL: <<LIBNAM>>
ISO: 3 H1H2O DD2O O16
CALC WGT% D2O >>MWGTH1<< >>MWGTD2<< >>MWGTO16<< ;
ECHO "PARAMETERS USED" ;
ECHO FWGTU5 FWGTU8 FWGTO16 ;
ECHO CWGTH1 CWGTD2 CWGTO16 ;
ECHO MWGTH1 MWGTD2 MWGTO16 ;
ECHO XWGT5 SWGT9 NWGT9 ;
ECHO BWGT0 BWGT1 GWGT5 GWGT7 ;
ECHO "PARAMETERS USED SUITE" ;
ECHO PCool PMode Bor Xe Sm Np Gad ;
ECHO DCool DMode Fdens Shdens PTdens CTdens Gapdens ;
ECHO TCool TMode TFuel Shtemp PTtemp CTtemp Gaptemp ;
ECHO "LIBRARY CREATION STARTED" ;
*----
* Create library
*----
LIBR := LIB: ::
EDIT 0
NMIX 10
CTRA WIMS
DEPL LIB: WIMSAECL FIL: <<LIBNAM>>
MIXS LIB: WIMSAECL FIL: <<LIBNAM>>
* Coolant
MIX 1 <<TCool>> <<DCool>>
CH1 = H1H2O <<CWGTH1>>
CD2 = DD2O <<CWGTD2>>
CO16 = O16 <<CWGTO16>>
* Pressure tube
MIX 2 <<PTtemp>> <<PTdens>>
ZR90PT 49.3443000
ZR91PT 10.8805220
ZR92PT 16.8136710
ZR94PT 17.4101970
ZR96PT 2.8645910
NB93 2.5800000
FE54 0.0026659
FE56 0.0429768
FE57 0.0010016
FE58 0.0001359
CR50 0.0003376
CR52 0.0067696
CR53 0.0007823
CR54 0.0001984
NI58 0.0024819
NI60 0.0009815
NI64 0.0000365
B10 0.0000243
* Gap
MIX 3 <<Gaptemp>> <<Gapdens>>
C 27.11
O16 72.89
* Calendria tube
MIX 4 <<CTtemp>> <<CTdens>>
ZR90CT 49.7984827
ZR91CT 10.9806700
ZR92CT 16.9684300
ZR94CT 17.5704469
ZR96CT 2.8909577
FE54 0.0076934
FE56 0.1240245
FE57 0.0028905
FE58 0.0003921
NI58 0.0390021
NI60 0.0154242
NI64 0.0005736
CR50 0.0041737
CR52 0.0836990
CR53 0.0096726
CR54 0.0024534
B10 0.0000596
* Moderator
MIX 5 <<TMode>> <<DMode>>
MH1 = H1H2O <<MWGTH1>>
MD2 = DD2O <<MWGTD2>>
MO16 = O16 <<MWGTO16>>
MB10 = B10 <<BWGT0>>
MB11 = B11 <<BWGT1>>
MG5 = GD155 <<GWGT5>>
MG7 = GD157 <<GWGT7>>
* Central Ring Fuel
MIX 6 <<TFuel>> <<Fdens>>
FU3 = U233 0.0 1
FU4 = U234 0.0 1
FU5 = U235 <<FWGTU5>> 1
FU6 = U236 0.0 1
FU8 = U238 <<FWGTU8>> 1
FO16 = O16 <<FWGTO16>>
FP9 = PU239 0.0 1
FP0 = PU240 0.0 1
FP1 = PU241 0.0 1
FP2 = PU242 0.0 1
AM241 0.0 1
AM242M 0.0 1
AM243 0.0 1
NP9 = NP239 <<NWGT9>>
XE5 = XE135 <<XWGT5>>
SM9 = SM149 <<SWGT9>>
* Inner Ring Fuel
MIX 7 COMB 6 1.0
* Intermediate Ring Fuel
MIX 8 COMB 6 1.0
* Outer Ring Fuel
MIX 9 COMB 6 1.0
* Sheath
MIX 10 <<Shtemp>> <<Shdens>>
ZR90CL 49.7847918
ZR91CL 10.9776511
ZR92CL 16.9637650
ZR94CL 17.5656164
ZR96CL 2.8901629
FE54 0.0119675
FE56 0.1929270
FE57 0.0044963
FE58 0.0006100
NI58 0.0049639
NI60 0.0019631
NI64 0.0000730
CR50 0.0041737
CR52 0.0836990
CR53 0.0096726
CR54 0.0024534
B10 0.0000596
;
ECHO "LIBRARY CREATION FINISHED" ;
END: ;
|
