blob: 22c013b6389c48548ba356571e5e5ce9b89685d4 (
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
|
***********************************************************
* Input file : test_channBWR.x2m *
* AUTHOR: Augusto Hernandez Solis (KTH) *
* auhs@kth.se *
* *
* Test of a BWR fuel assembly with real power conditions *
* *
* The aim is to test the THM module of DONJON *
***********************************************************
LINKED_LIST GEOM GEOM2 MATEX FMAP THM ;
MODULE GEO: USPLIT: RESINI: THM: END: ;
PROCEDURE assertS ;
GEOM := GEO: :: CAR3D 1 1 27
EDIT 5
X- VOID X+ VOID
Y- VOID Y+ VOID
Z- VOID Z+ VOID
MESHX 0.000 15.375
MESHY 0.000 15.375
MESHZ 0.000 14.720 29.440 44.160 58.880 73.600
88.320 103.040 117.760 132.480 147.200 161.920
176.640 191.360 206.080 220.800 235.520 250.240
264.960 279.680 294.400 309.120 323.840 338.560
353.280 368.000 382.720 397.440
MIX
PLANE 1 1
PLANE 2 2
PLANE 3 SAME 2
PLANE 4 SAME 2
PLANE 5 SAME 2
PLANE 6 SAME 2
PLANE 7 SAME 2
PLANE 8 SAME 2
PLANE 9 SAME 2
PLANE 10 SAME 2
PLANE 11 SAME 2
PLANE 12 SAME 2
PLANE 13 SAME 2
PLANE 14 SAME 2
PLANE 15 SAME 2
PLANE 16 SAME 2
PLANE 17 SAME 2
PLANE 18 SAME 2
PLANE 19 SAME 2
PLANE 20 SAME 2
PLANE 21 SAME 2
PLANE 22 SAME 2
PLANE 23 SAME 2
PLANE 24 SAME 2
PLANE 25 SAME 2
PLANE 26 SAME 2
PLANE 27 SAME 1
;
GEOM2 MATEX := USPLIT: GEOM :: EDIT 5 NGRP 2 MAXR 27
NREFL 1 RMIX 1
NFUEL 1 FMIX 2
;
FMAP MATEX := RESINI: MATEX ::
::: GEO: CAR3D 1 1 27
EDIT 5
X- VOID X+ VOID
Y- VOID Y+ VOID
Z- VOID Z+ VOID
MESHX 0.000 15.375
MESHY 0.000 15.375
MESHZ 0.000 14.720 29.440 44.160 58.880 73.600
88.320 103.040 117.760 132.480 147.200 161.920
176.640 191.360 206.080 220.800 235.520 250.240
264.960 279.680 294.400 309.120 323.840 338.560
353.280 368.000 382.720 397.440
MIX
PLANE 1 0
PLANE 2 2
PLANE 3 SAME 2
PLANE 4 SAME 2
PLANE 5 SAME 2
PLANE 6 SAME 2
PLANE 7 SAME 2
PLANE 8 SAME 2
PLANE 9 SAME 2
PLANE 10 SAME 2
PLANE 11 SAME 2
PLANE 12 SAME 2
PLANE 13 SAME 2
PLANE 14 SAME 2
PLANE 15 SAME 2
PLANE 16 SAME 2
PLANE 17 SAME 2
PLANE 18 SAME 2
PLANE 19 SAME 2
PLANE 20 SAME 2
PLANE 21 SAME 2
PLANE 22 SAME 2
PLANE 23 SAME 2
PLANE 24 SAME 2
PLANE 25 SAME 2
PLANE 26 SAME 2
PLANE 27 SAME 1
;
NXNAME '01'
NYNAME 'A'
NCOMB 1
B-ZONE 1
BTYPE INST-BURN
INST-BVAL BUND
1600.000 7000.000 9800.000 10800.000 11300.000
11400.000 11300.000 10900.000 10900.000 10800.000
10600.000 10300.000 10300.000 10100.000 10400.000
10200.000 10400.000 10200.000 9900.000 9400.000
8900.000 8100.000 6700.000 4900.000 1700.000
BUNDLE-POW BUND
43.236 211.550 287.214 305.744 311.920
307.288 301.111 290.302 288.758 282.581
274.861 264.051 259.419 250.154 231.624
222.359 219.271 210.006 197.653 183.755
168.313 146.695 117.356 78.752 20.074
ADD-PARAM PNAME 'C-BORE' PARKEY 'ppmBore' GLOBAL
ADD-PARAM PNAME 'T-FUEL' PARKEY 'TF' GLOBAL
ADD-PARAM PNAME 'T-COOL' PARKEY 'TCA' GLOBAL
ADD-PARAM PNAME 'D-COOL' PARKEY 'DCA' GLOBAL
;
INTEGER iter := 1 ;
WHILE iter 10 < DO
IF iter 1 = THEN
THM FMAP := THM: FMAP ::
EDIT 10
CRITFL 1.80E+6 (* W/m^2 *)
ASSMB 91 0
RADIUS 0.0041 0.0044 0.0050 0.0 (* m *)
RODMESH 5 10
SPEED 1.5939
SAHA
INLET 7090000.0000 (* Pa *) 550.150 (* K *)
RELAX 0.5
SET-PARAM 'T-COOL' 500.0
SET-PARAM 'T-FUEL' 1000.0
SET-PARAM 'D-COOL' 0.75
;
ELSE
THM FMAP := THM: THM FMAP :: EDIT 5 ;
ENDIF ;
EVALUATE iter := iter 1 + ;
ENDWHILE ;
assertS THM :: "MIN-T-FUEL" 1 583.3722 ;
assertS THM :: "MAX-T-FUEL" 1 928.3032 ;
ECHO "test test_channBWR completed" ;
END: ;
|