summaryrefslogtreecommitdiff
path: root/Donjon/data/test_channBWR.x2m
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: ;