blob: 56dafcbae64483f62c370315fb06c234fcd5a4a2 (
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
|
*----
* TEST CASE FOR TINST:
* 17 X 17 PWR ASSEMBLY WITHOUT POISON
* DEPLETION WITH MULTICOMPO DATABASE
* 172-GROUP JEF2.2 DRAGLIB
*
* REF: TCDR14.x2m
* Geometry is the same as TCDR14 except that it is unfolded
* since DONJON can't deal with DIAG and SYME boundary conditions.
* Besides DONJON needs a CARD3D geometry with equally sized cells
* (RESINI: can't deal with bundle of different sizes) so the water
* gap is not included here
*
* A mixture number 3 has been included because we need at least
* one reflector material otherwise USPLIT: is unhappy
*
* AUTHOR : B. Toueg
* DATE : 2010/03/10
*
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST
MULTICOMPO ASSMB_HOM MATEX CORETRACK FUELMAP MACFL MACNOEVO
MACROEXT MATRIX COREFLUX POWER ;
MODULE
GEO: USPLIT: TRIVAT: RESINI: NCR: MACINI: TRIVAA: FLUD: FLPOW: TINST:
DELETE: END: ;
PROCEDURE assertS ;
* loading the L_MULTICOMPO generated with TCDR14.x2m
SEQ_ASCII saCOMPO :: FILE 'saCOMPO' ;
MULTICOMPO := saCOMPO ;
ASSMB_HOM := GEO: :: CAR3D 17 17 1
X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL
MIX 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 1 2 2 1 2 2 1 2 2 2 2 2
2 2 2 3 2 2 2 2 2 2 2 2 2 1 2 2 2 ! 1 dummy reflector material
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 1 2 2 2 2 2 2 2 2 2 1 2 2 2
2 2 2 2 2 1 2 2 1 2 2 1 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
MESHX -10.11776
-8.85304 -7.58832 -6.3236 -5.05888 -3.79416 -2.52944 -1.26472
0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304
10.11776 11.38248 ! and not 11.43248 because no water gap
MESHY -10.11776
-8.85304 -7.58832 -6.3236 -5.05888 -3.79416 -2.52944 -1.26472
0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304
10.11776 11.38248 ! and not 11.43248 because no water gap
MESHZ 0.0 1.26472
;
INTEGER MaxR := 17 17 * 1 * ;
ASSMB_HOM MATEX := USPLIT: ASSMB_HOM :: NGRP 2 MAXR <<MaxR>>
NREFL 1 RMIX 3
NFUEL 2 FMIX 1 2
;
INTEGER MCFD := 3 ;
CORETRACK := TRIVAT: ASSMB_HOM :: MAXR <<MaxR>> MCFD <<MCFD>> ;
FUELMAP MATEX := RESINI: MATEX ::
EDIT 3
:::
GEO: CAR3D 17 17 1
X- REFL X+ REFL Y- REFL Y+ REFL Z- REFL Z+ REFL
MIX 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 1 2 2 1 2 2 1 2 2 2 2 2
2 2 2 0 2 2 2 2 2 2 2 2 2 1 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 1 2 2 1 2 2 1 2 2 1 2 2 1 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 1 2 2 2 2 2 2 2 2 2 1 2 2 2
2 2 2 2 2 1 2 2 1 2 2 1 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
MESHX -10.11776
-8.85304 -7.58832 -6.3236 -5.05888 -3.79416 -2.52944 -1.26472
0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304
10.11776 11.38248 ! and not 11.43248 because no water gap
MESHY -10.11776
-8.85304 -7.58832 -6.3236 -5.05888 -3.79416 -2.52944 -1.26472
0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304
10.11776 11.38248 ! and not 11.43248 because no water gap
MESHZ 0.0 1.26472
;
NXNAME
'1' '2' '3' '4' '5' '6' '7' '8'
'9'
'10' '11' '12' '13' '14' '15' '16' '17'
NYNAME
'A' 'B' 'C' 'D' 'E' 'F' 'G' 'H'
'I'
'K' 'L' 'M' 'N' 'O' 'P' 'Q' 'R''
NCOMB ALL
FUEL WEIGHT 460. 460. ! arbitrary values needed for TINST:
;
FUELMAP := RESINI: FUELMAP ::
EDIT 0
BTYPE INST-BURN
INST-BVAL
SAME 0.
;
REAL Power := 31.9713 ; ! value taken from TCDR14
INTEGER iter max_iter := 0 5 ;
WHILE iter max_iter < DO
MACFL := NCR: MULTICOMPO FUELMAP ::
EDIT 0
MACRO
LINEAR
TABLE MULTICOMPO default 'BURN'
MIX 1 FROM 1
MICRO ALL
ENDMIX
MIX 2 FROM 2 ! can't be 1 !?
MICRO ALL
ENDMIX
;
MACNOEVO := NCR: MULTICOMPO ::
EDIT 0
MACRO
LINEAR
NMIX 3
COMPO MULTICOMPO default
MIX 3 FROM 1
SET 'BURN' 0.
MICRO ALL
ENDMIX
;
MACROEXT MATEX := MACINI: MATEX MACNOEVO MACFL :: EDIT 1 ;
MATRIX := TRIVAA: MACROEXT CORETRACK :: EDIT 1 ;
COREFLUX := FLUD: MATRIX CORETRACK ::
EDIT 2 ACCE 5 3 EXTE 1E-5 ADI 5 ;
POWER := FLPOW: FUELMAP COREFLUX CORETRACK MATEX ::
EDIT 1
PTOT <<Power>>
;
FUELMAP := TINST: FUELMAP POWER ::
EDIT 1
BURN-STEP 30.
;
EVALUATE iter := iter 1 + ;
IF iter max_iter < THEN
MACFL MACNOEVO MACROEXT MATRIX COREFLUX POWER := DELETE:
MACFL MACNOEVO MACROEXT MATRIX COREFLUX POWER ;
ENDIF ;
ENDWHILE ;
assertS COREFLUX :: 'K-INFINITY' 1 1.221152 ;
ECHO "test TINST completed" ;
END: ;
QUIT "LIST" .
|
