blob: 8983c9f6f1dfba987a4897716de1df2bd9443e9e (
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
|
*----
* TEST CASE TRDR22
* 17 X 17 PWR VOIDED ASSEMBLY
* 172-GROUP JEF2.2 DRAGLIB
* Author : Alain Hebert (2019)
*
* REF: none
*
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST
ASSMB DISCR LIBRARY CP CALC OUT SPHGEOM MTRACK MACRO MACRO2 SYS FLUX
EDIT ;
MODULE
LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: FLUD: MAC: DELETE: UTL: BIVACT:
BIVACA: SPH: OUT: ERROR: END: ;
PROCEDURE assertS ;
*----
* Microscopic cross sections from Jef 2.2 Draglib
*----
LIBRARY := LIB: ::
NMIX 8 CTRA WIMS ANIS 2
MIXS LIB: DRAGON FIL: DLIB_J2
MIX 1 579.9
H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2
B10 = B10 4.71444E-4 B11 = B11 1.90959E-3
MIX 2 579.9 O16 = O16 3.06711E-4
Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4
Zr91 = Zr0 4.18621E-2
MIX 3 579.9
H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2
Cr52 = Cr52 4.79927E-5
Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4
MoNat = Mo95 4.03755E-6
Al27 = Al27 2.35231E-6
Mn55 = Mn55 4.15901E-7
B10 = B10 4.60867E-6 B11 = B11 1.86674E-5
Zr91 = Zr0 8.92427E-4
MIX 4 579.9 O16 = O16 2.87335E-4
Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4
Zr91 = Zr0 3.92175E-2
MIX 5 579.9
H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2
Cr52 = Cr52 2.24991E-5
Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5
MoNat = Mo95 1.89281E-6
Al27 = Al27 1.10277E-6
Mn55 = Mn55 1.94976E-7
B10 = B10 4.66484E-6 B11 = B11 1.88950E-5
Zr91 = Zr0 4.18372E-4
MIX 6 579.9
H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2
Cr52 = Cr52 2.11122E-5
Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5
MoNat = Mo95 1.77614E-6
Al27 = Al27 1.03479E-6
Mn55 = Mn55 1.82957E-7
B10 = B10 4.66791E-6 B11 = B11 1.89074E-5
Zr91 = Zr0 3.92583E-4
MIX 7 579.9
H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2
Cr52 = Cr52 1.96591E-5
Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5
MoNat = Mo95 1.65389E-6
Al27 = Al27 9.63569E-7
Mn55 = Mn55 1.70365E-7
B10 = B10 4.67110E-6 B10 = B10 1.89203E-5
Zr91 = Zr0 3.65562E-4
MIX 8 933.6 O16 = O16 4.49355E-2
U235 = U235 7.39237E-4 1
U238 = U238 2.17285E-2 1
;
*----
* Geometry ASSMB : a 17 X 17 normal PWR assembly
* contains C1 : cell without fuel
* C2 : normal fuel cell
* C3 : peripheral cell
* C4 : corner cell
*----
ASSMB := GEO: :: CAR2D 9 9
X- DIAG X+ REFL Y- SYME Y+ DIAG
CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
C2 C2 C2 C2 C2 C2 C2 C3
C2 C2 C2 C2 C2 C2 C3
C1 C2 C2 C1 C2 C3
C2 C2 C2 C2 C3
C1 C2 C2 C3
C2 C2 C3
C2 C3
C4
MERGE 1 3 12 11 12 12 11 12 15
4 6 5 6 6 5 6 8
13 5 6 6 5 6 8
2 5 5 10 5 8
13 5 5 6 8
2 5 7 8
13 7 8
14 8
9
TURN H H B H H B H H A
H G G H G G H A
A E E F E E A
H H F H H A
H E G H A
H H A A
H A A
A A
A
::: C1 := GEO: CARCEL 2
MESHX 0.0 1.26472 MESHY 0.0 1.26472
RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ;
::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ;
::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ;
::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ;
;
*----
* Self-Shielding calculation SYBIL
* Transport calculation SYBIL
* Flux calculation for P1 homogeneous leakage
*----
DISCR := SYBILT: ASSMB ::
TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
MAXR 400 QUA2 6 3 ;
LIBRARY := SHI: LIBRARY DISCR :: LJ ;
ASSMB DISCR := DELETE: ASSMB DISCR ;
*----
* New lattice geometry with VOID boundary condition
*----
ASSMB := GEO: :: CAR2D 9 9
X- DIAG X+ VOID Y- SYME Y+ DIAG
CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
C2 C2 C2 C2 C2 C2 C2 C3
C2 C2 C2 C2 C2 C2 C3
C1 C2 C2 C1 C2 C3
C2 C2 C2 C2 C3
C1 C2 C2 C3
C2 C2 C3
C2 C3
C4
::: C1 := GEO: CARCEL 2
MESHX 0.0 1.26472 MESHY 0.0 1.26472
RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ;
::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ;
::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ;
::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ;
;
DISCR := SYBILT: ASSMB ::
TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
MAXR 400 QUA2 6 3 ;
CP := ASM: LIBRARY DISCR :: ARM ;
CALC := FLU: CP LIBRARY DISCR ::
TYPE K P1 HETE ;
assertS CALC :: 'K-EFFECTIVE' 1 0.4551677 ;
*----
* Editing using SPH model for transport-diffusion with flux-volume
* normalization
*----
OUT := EDI: LIBRARY DISCR CALC ASSMB ::
EDIT 4 UPS SAVE MICR RES MERGE CELL COND 4.0
;
*----
* TRANSPORT-DIFFUSION SPH EQUIVALENCE IN NON-FUNDAMENTAL MODE
*----
SPHGEOM := GEO: :: CAR2D 9 9
X- DIAG X+ ALBE 1 Y- SYME Y+ DIAG
MIX 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
MESHX 0.0 1.26471996 2.52943993 3.79415989 5.05887985
6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
MESHY 0.0 1.26471996 2.52943993 3.79415989 5.05887985
6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
;
MTRACK := BIVACT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 1 ;
MACRO := OUT :: STEP UP 'REF-CASE0001' STEP UP 'MACROLIB' ;
*----
* Set physical albedos
*----
MACRO := MAC: MACRO :: ALBP 1 0.0 0.0 ;
*----
* Perform SPH equivalence
*----
MACRO2 := SPH: MACRO MTRACK :: EDIT 2 ITER 1.0E-7 ;
*----
* Verification calculation
*----
SYS := BIVACA: MACRO2 MTRACK ;
FLUX := FLUD: SYS MTRACK :: EDIT 2 ADI 20 ACCE 5 3 ;
EDIT := OUT: FLUX SPHGEOM MACRO2 MTRACK ::
EDIT 0 INTG MIX ;
ECHO "Error statistics after correction" ;
ERROR: MACRO EDIT :: HREA 'NFTOT' ;
assertS FLUX :: K-EFFECTIVE 1 0.4551687 ;
ECHO "test TCDR22 completed" ;
END: ;
QUIT "LIST" .
|
