summaryrefslogtreecommitdiff
path: root/Dragon/data/tdraglib_proc/TCDR15.c2m
blob: 747d4d23044d84a8ec10a9770420899c4671dba0 (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
*--------------------------------------------------------------------
* TEST CASE TCDR15
* HTTR Fuel
* Block ID:f393320   
* JAVIER ORTENSI
* Jef2.2 evaluation
* XMAS172 Group Structure
* Description:  models a fuel block in infinite medium with spatial
*               homogenization of the burnable poison region (BP x-section)
*               and the rest of the assembly (background x-section).
*--------------------------------------------------------------------
*  Define STRUCTURES and MODULES used
*--------------------------------------------------------------------
LINKED_LIST LATTIC LIBRARY LIBRARY2 OUT1G OUT2G FLUX TRACK SYS ;
SEQ_ASCII EDIT.out ;
SEQ_BINARY TRACK_FIL ;
PROCEDURE assertS ;

MODULE GEO: EXCELT: LIB: ASM: FLU: EDI: USS: END: DELETE: ;
*--------------------------------------------------------------------
*  Define variables used
*--------------------------------------------------------------------
LIBRARY := LIB: ::
  NMIX 14 CTRA WIMS
  SUBG
  ANIS 2
  MIXS LIB: DRAGON FIL: DLIB_J2
* Fuel Kernel
  MIX  1 300.0
     B10    = 'B10'   1.1010E-07
     B11    = 'B11'   4.4317E-07
     O16    = 'O16'   4.8150E-02
     O17    = 'O17'   1.8304E-05
    U235    = 'U235'  8.0367E-04 1 IRSET 0.0 81
    U238    = 'U238'  2.3280E-02 1 IRSET 0.0 81
* Porous C buffer
  MIX  2 300.0
    C12     = 'C0_GR'   5.8060E-02
* IPyC
  MIX  3 300.0
    C12     = 'C0_GR'   9.3909E-02
* SiC coating layer
  MIX  4 300.0
     C12    = 'C0'     4.8151E-02
     SI28   = 'Si0'    4.4410E-02
     SI29   = 'Si0'    2.2487E-03
     SI30   = 'Si0'    1.4927E-03
* OPyC
  MIX  5 300.0
     C12    = 'C0_GR'   9.3007E-02
* Graphite matrix (Compact)
  MIX 6 300.0
     C12    = 'C0_GR'   8.5085E-02
     B10    = 'B10'     2.4543E-08
     B11    = 'B11'     9.8790E-08
* Graphite sleeve
  MIX 7 300.0
     C12    = 'C0_GR'   8.8745E-02
     B10    = 'B10'     7.1811E-09
     B11    = 'B11'     2.8905E-08
* Helium Coolant
  MIX 8 300.0
     HE3    = 'He3'     1.3230E-09
     HE4    = 'He4'     9.6573E-04
* Block Graphite inside fuel cell
  MIX 9 300.0
     C12    = 'C0_GR'   8.8745E-02
     B10    = 'B10'     7.8482E-09
     B11    = 'B11'     3.1590E-08
* Block Graphite around fuel handling hole (rho adjusted)
  MIX 10 300.0
     C12    = 'C0_GR'   6.7791E-02
     B10    = 'B10'     5.9951E-09
     B11    = 'B11'     2.4131E-08
     HE3    = 'He3'     3.1239E-10
     HE4    = 'He4'     2.2802E-04
* Block Graphite outside fuel cell (rho adjusted)
  MIX 11 300.0
     C12    = 'C0_GR'   9.0845E-02
     B10    = 'B10'     8.0339E-09
     B11    = 'B11'     3.2338E-08
* Burnable Poison
  MIX 12 300.0
     C12    = 'C0_GR'   8.7821E-02
     B10    = 'B10'     3.5588E-04
     B11    = 'B11'     1.4071E-03
* Helium Coolant
  MIX 13 300.0
     HE3    = 'He3'     1.3230E-09
     HE4    = 'He4'     9.6573E-04
* Block Graphite around BP
  MIX 14 300.0
     C12    = 'C0_GR'   8.8745E-02
     B10    = 'B10'     7.8482E-09
     B11    = 'B11'     3.1590E-08
;
*--------------------------------------------------------------------
*  FUEL BLOCK SECTION
*--------------------------------------------------------------------
*  Define mixture for the fuel block
*--------------------------------------------------------------------
LATTIC := GEO: :: HEX 61
  HBC COMPLETE REFL
  CELL
* Ring 0 (center cell)
FH
* Ring 1  (6 cells)
F F F F F F
* Ring 2 (12 cells)
F F F F F F F F F F F F
* Ring 3 (18 cells)
BPE F F F F F BP F F F F F BP F F F F F
* Ring 4 (24 cells)
C C C C C C C C C C C C C C C C C C C C C C C C
* Spherical microstruture model
  BIHET SPHE  1  1
  (* NS     *) 5
  (* RS     *) 0.0 0.0297 0.0359 0.0389 0.0418 0.0464 
  (* milie  *)  15
  (* mixdil *)  6
  (* loop over NMISTR *)
  (* fract  *) 0.3040
  (* mixgr  *)  1 2 3 4 5
*  Side length based on the pin pitch
*    Fuel cell
 ::: F := GEO: HEXCEL 5
     SIDE 2.973354
     RADIUS 0.0 0.4995 1.3000 1.3130 1.7000 1.9800 
     MIX 8 15 8 7 8 9 ;
*  Fuel handling hole
 ::: FH := GEO: HEXCEL 0
     SIDE 2.973354
     MIX 10 ;
*  Burnable Poison Cell (Empty)
*  To be homogenized with fuel (use mat 8&9)
 ::: BPE := GEO: HEXCEL 1
     SIDE 2.973354
     RADIUS 0.0 0.6950
     MIX 8 9 ;
*  Burnable Poison Cell
 ::: BP := GEO: HEXCEL 2
     SIDE 2.973354
     RADIUS 0.0 0.6950 0.7500
     MIX 12 13 14 ;
* graphite cell
 ::: C := GEO: HEXCEL 0
     SIDE 2.973354
     MIX 11 ;
;
*--------------------------------------------------------------------
*  Use SYBILT, USS, ASM and FLU for flux solution
*--------------------------------------------------------------------
*----
*  Tracking calculation
*----
TRACK TRACK_FIL := EXCELT: LATTIC ::
    EDIT 0
    MAXR 400
    ANIS 2
    TRAK TISO 8 20.0
;
*---
*  Self Shielding calculation
*---
LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL ::
  EDIT 1
;
*---
*  Multigroup Flux Calculation
*---
 SYS := ASM: LIBRARY2 TRACK TRACK_FIL ::
  EDIT 1 ;
 FLUX := FLU: SYS LIBRARY2 TRACK TRACK_FIL ::
  EDIT 1 
  TYPE B B1 PNL
  EXTE 100 1E-5 THER 5 1E-5 ;
 OUT1G := EDI: FLUX LIBRARY2 TRACK LATTIC ::
  EDIT 3
 MERG MIX 1 1 1 1 1 1 1 1 1  1  1  0  0  0
 COND 7.4080E06 3.6790E06 6.7210E05 1.1110E05 1.9310E04 3.3550E03 
      1.5850E03 7.4850E02 2.7540E02 1.3010E02 6.1440E01 2.9000E01 
      1.3700E01 8.3200E00 5.0400E00 2.3800E00 1.2900E00 6.5000E-01 
      3.5000E-01 2.0000E-01 1.2000E-01 8.0000E-02 5.0000E-02 2.0000E-02 
      1.0000E-02 
 MICR ISOTXS ASCII ALL
 SAVE ;
 OUT1G := EDI: OUT1G FLUX LIBRARY2 TRACK ::
  EDIT 3
 MERG MIX 0 0 0 0 0 0 0 0 0  0  0  1  1  1
 COND 7.4080E06 3.6790E06 6.7210E05 1.1110E05 1.9310E04 3.3550E03 
      1.5850E03 7.4850E02 2.7540E02 1.3010E02 6.1440E01 2.9000E01 
      1.3700E01 8.3200E00 5.0400E00 2.3800E00 1.2900E00 6.5000E-01 
      3.5000E-01 2.0000E-01 1.2000E-01 8.0000E-02 5.0000E-02 2.0000E-02 
      1.0000E-02 
 MICR ISOTXS ASCII ALL
 SAVE ;
assertS FLUX :: K-INFINITY 1 1.144844 ;
ECHO "test TCDR15 completed" ;
EDIT.out := OUT1G ; 
END: ;