Unexpected Behavior of Axial Power Distribution

Report any suspected bugs and unexpected behavior here

Unexpected Behavior of Axial Power Distribution

Postby Yousef1986 » Thu Feb 18, 2016 7:14 pm

Dear Jaakko,

I have encountered an unexpected behaviour of axial power distribution when I used the card "set cpd"

The syntax was:
set cpd 2 200 -100 100

What I found is the following:
- The 2D pin power distribution is realistic and makes sense. However, the 3D pin power distribution results in non-realistic values (I attached part of _core0.m file below for your reference; Sorry for that I could not upload the whole file).
For example, check-out pin no. 20, the power should be zero there and should not be printed but in the axial power distribution, it appearers with some power value. For the 2D power distribution case, pin no. 20 does not appear there which make since pin no. 20 has no fissile material.

Code: Select all
% ----- Distribution at level 0 (first lattice):

lvl0 = [
   1 8.33374E+03  0.00061 % a
   2 7.96612E+03  0.00058 % a
   3 7.56476E+03  0.00060 % a
   4 7.83344E+03  0.00059 % a
   5 7.84243E+03  0.00061 % a
   6 8.52551E+03  0.00055 % a
   7 8.48892E+03  0.00059 % a
   8 7.96790E+03  0.00057 % a
   9 8.17519E+03  0.00059 % a
  10 7.96723E+03  0.00060 % a
  11 8.48705E+03  0.00058 % a
  12 8.53867E+03  0.00056 % a
  13 7.83770E+03  0.00061 % a
  14 7.83747E+03  0.00059 % a
  15 7.56480E+03  0.00062 % a
  16 7.96484E+03  0.00060 % a
  17 8.33446E+03  0.00060 % a
  18 7.97284E+03  0.00064 % a
  19 7.15430E+03  0.00059 % a
  21 7.45144E+03  0.00060 % a
  23 8.99917E+03  0.00056 % a
  24 8.56027E+03  0.00055 % a
  26 8.30190E+03  0.00057 % a
  28 8.56402E+03  0.00056 % a
  29 9.00784E+03  0.00054 % a
  31 7.44293E+03  0.00058 % a
  33 7.17019E+03  0.00061 % a
  34 7.97026E+03  0.00058 % a
  35 7.55498E+03  0.00062 % a
  37 7.60476E+03  0.00061 % a
  38 9.03548E+03  0.00054 % a
  39 9.36174E+03  0.00056 % a
  41 9.75323E+03  0.00054 % a
  42 9.21348E+03  0.00055 % a
  44 9.21723E+03  0.00055 % a
  45 9.75545E+03  0.00055 % a
  47 9.35110E+03  0.00056 % a
  48 9.03238E+03  0.00055 % a
  49 7.61392E+03  0.00060 % a
  51 7.57139E+03  0.00061 % a
  52 7.83395E+03  0.00059 % a
  53 7.44341E+03  0.00060 % a
  54 9.03384E+03  0.00054 % a
  56 9.62463E+03  0.00053 % a
  57 1.00700E+04  0.00053 % a
  58 9.45090E+03  0.00053 % a
  59 8.86347E+03  0.00054 % a
  60 9.50140E+03  0.00052 % a
  61 8.85291E+03  0.00055 % a
  62 9.44924E+03  0.00052 % a
  63 1.00794E+04  0.00051 % a
  64 9.62678E+03  0.00056 % a
  66 9.03053E+03  0.00054 % a
  67 7.45600E+03  0.00061 % a
  68 7.83555E+03  0.00062 % a
  69 7.85006E+03  0.00061 % a
  71 9.35283E+03  0.00053 % a
  72 9.62096E+03  0.00053 % a
  73 1.27283E+03  0.00049 % a
  74 9.68089E+03  0.00052 % a
  75 9.07794E+03  0.00053 % a
  76 1.25990E+03  0.00052 % a
  77 8.68341E+03  0.00056 % a
  78 1.25940E+03  0.00051 % a
  79 9.06966E+03  0.00053 % a
  80 9.66774E+03  0.00052 % a
  81 1.27373E+03  0.00050 % a
  82 9.62320E+03  0.00052 % a
  83 9.36069E+03  0.00055 % a
  85 7.85234E+03  0.00061 % a
  86 8.52785E+03  0.00059 % a
  87 8.99773E+03  0.00054 % a
  89 1.00674E+04  0.00052 % a
  90 9.67528E+03  0.00049 % a
  92 9.96603E+03  0.00051 % a
  93 9.15103E+03  0.00053 % a
  95 9.14425E+03  0.00052 % a
  96 9.97824E+03  0.00052 % a
  98 9.67887E+03  0.00054 % a
  99 1.00703E+04  0.00051 % a
101 9.00715E+03  0.00056 % a
102 8.52614E+03  0.00060 % a
103 8.48682E+03  0.00058 % a
104 8.56715E+03  0.00056 % a
105 9.74360E+03  0.00053 % a
106 9.44533E+03  0.00054 % a
107 9.08930E+03  0.00054 % a
108 9.96983E+03  0.00052 % a
109 9.45875E+03  0.00055 % a
110 8.78466E+03  0.00054 % a
112 8.78122E+03  0.00057 % a
113 9.45401E+03  0.00053 % a
114 9.96958E+03  0.00052 % a
115 9.08076E+03  0.00057 % a
116 9.44567E+03  0.00054 % a
117 9.74953E+03  0.00053 % a
118 8.56572E+03  0.00060 % a
119 8.48675E+03  0.00058 % a
120 7.96253E+03  0.00057 % a
122 9.21067E+03  0.00057 % a
123 8.84948E+03  0.00055 % a
124 1.26028E+03  0.00049 % a
125 9.13079E+03  0.00052 % a
126 8.78311E+03  0.00057 % a
127 8.96488E+03  0.00053 % a
128 9.17231E+03  0.00054 % a
129 8.96703E+03  0.00055 % a
130 8.78610E+03  0.00055 % a
131 9.14305E+03  0.00053 % a
132 1.25923E+03  0.00049 % a
133 8.86027E+03  0.00057 % a
134 9.21368E+03  0.00054 % a
136 7.96790E+03  0.00056 % a
137 8.17286E+03  0.00059 % a
138 8.28763E+03  0.00055 % a
140 9.48809E+03  0.00054 % a
141 8.68074E+03  0.00056 % a
144 9.17039E+03  0.00056 % a
146 9.16713E+03  0.00055 % a
149 8.67342E+03  0.00055 % a
150 9.48969E+03  0.00053 % a
152 8.29956E+03  0.00059 % a
153 8.17509E+03  0.00060 % a
154 7.97191E+03  0.00056 % a
156 9.21668E+03  0.00057 % a
157 8.86096E+03  0.00051 % a
158 1.25983E+03  0.00049 % a
159 9.13609E+03  0.00055 % a
160 8.79538E+03  0.00054 % a
161 8.96751E+03  0.00053 % a
162 9.17400E+03  0.00055 % a
163 8.96938E+03  0.00055 % a
164 8.78055E+03  0.00056 % a
165 9.13883E+03  0.00055 % a
166 1.25964E+03  0.00048 % a
167 8.85203E+03  0.00053 % a
168 9.22009E+03  0.00056 % a
170 7.96139E+03  0.00061 % a
171 8.49301E+03  0.00059 % a
172 8.57132E+03  0.00058 % a
173 9.75946E+03  0.00054 % a
174 9.45188E+03  0.00052 % a
175 9.08499E+03  0.00053 % a
176 9.97218E+03  0.00050 % a
177 9.46596E+03  0.00053 % a
178 8.78532E+03  0.00054 % a
180 8.78665E+03  0.00055 % a
181 9.45879E+03  0.00053 % a
182 9.96775E+03  0.00052 % a
183 9.08258E+03  0.00054 % a
184 9.44686E+03  0.00054 % a
185 9.75068E+03  0.00054 % a
186 8.56847E+03  0.00056 % a
187 8.48785E+03  0.00057 % a
188 8.52104E+03  0.00057 % a
189 9.00753E+03  0.00054 % a
191 1.00809E+04  0.00055 % a
192 9.68090E+03  0.00053 % a
194 9.97324E+03  0.00052 % a
195 9.14171E+03  0.00054 % a
197 9.14392E+03  0.00054 % a
198 9.96267E+03  0.00054 % a
200 9.67064E+03  0.00054 % a
201 1.00702E+04  0.00053 % a
203 9.00073E+03  0.00054 % a
204 8.52278E+03  0.00060 % a
205 7.85156E+03  0.00060 % a
207 9.35855E+03  0.00053 % a
208 9.62149E+03  0.00054 % a
209 1.27424E+03  0.00050 % a
210 9.67025E+03  0.00053 % a
211 9.07705E+03  0.00054 % a
212 1.26010E+03  0.00048 % a
213 8.67707E+03  0.00056 % a
214 1.25909E+03  0.00048 % a
215 9.07070E+03  0.00054 % a
216 9.68214E+03  0.00052 % a
217 1.27412E+03  0.00048 % a
218 9.63594E+03  0.00054 % a
219 9.36162E+03  0.00055 % a
221 7.84881E+03  0.00060 % a
222 7.84131E+03  0.00060 % a
223 7.45149E+03  0.00058 % a
224 9.03898E+03  0.00056 % a
226 9.62909E+03  0.00054 % a
227 1.00691E+04  0.00050 % a
228 9.44940E+03  0.00052 % a
229 8.85751E+03  0.00056 % a
230 9.49702E+03  0.00054 % a
231 8.84847E+03  0.00058 % a
232 9.44646E+03  0.00052 % a
233 1.00778E+04  0.00052 % a
234 9.62997E+03  0.00054 % a
236 9.03285E+03  0.00054 % a
237 7.45260E+03  0.00060 % a
238 7.83316E+03  0.00060 % a
239 7.56680E+03  0.00062 % a
241 7.61153E+03  0.00061 % a
242 9.02161E+03  0.00056 % a
243 9.37149E+03  0.00056 % a
245 9.76453E+03  0.00054 % a
246 9.21674E+03  0.00054 % a
248 9.22011E+03  0.00054 % a
249 9.74733E+03  0.00053 % a
251 9.35787E+03  0.00056 % a
252 9.03710E+03  0.00053 % a
253 7.60844E+03  0.00058 % a
255 7.57372E+03  0.00059 % a
256 7.95493E+03  0.00059 % a
257 7.16441E+03  0.00063 % a
259 7.45437E+03  0.00059 % a
261 8.99968E+03  0.00056 % a
262 8.57893E+03  0.00056 % a
264 8.30680E+03  0.00057 % a
266 8.56829E+03  0.00057 % a
267 9.00185E+03  0.00056 % a
269 7.44499E+03  0.00061 % a
271 7.16796E+03  0.00062 % a
272 7.96977E+03  0.00060 % a
273 8.33072E+03  0.00062 % a
274 7.96880E+03  0.00058 % a
275 7.56261E+03  0.00059 % a
276 7.84239E+03  0.00060 % a
277 7.86130E+03  0.00060 % a
278 8.52496E+03  0.00057 % a
279 8.48705E+03  0.00057 % a
280 7.96628E+03  0.00059 % a
281 8.17860E+03  0.00058 % a
282 7.96604E+03  0.00058 % a
283 8.50059E+03  0.00057 % a
284 8.52789E+03  0.00059 % a
285 7.84849E+03  0.00058 % a
286 7.83034E+03  0.00059 % a
287 7.56903E+03  0.00062 % a
288 7.96545E+03  0.00059 % a
289 8.33510E+03  0.00059 % a
];

% Number and fraction of bins with statistical errors
% below 0.01, 0.02, 0.05 and 0.10:

n0 = [ 236 236 236 236 ];
f0 = [ 1.00000 1.00000 1.00000 1.00000 ];

% Distribution of statistical errors (1% bins):

prof0 = [ 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 ];

% Cumulative distribution of statistical errors (1% bins):

cum0 = [ 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 ];

% ----- Distribution at level 2 (axial):

lvl2 = [
   1   1   1 4.22133E+01  0.00784 % r
   1   1   2 4.16052E+01  0.00773 % r
   1   1   3 4.21411E+01  0.00785 % r
   1   1   4 4.22192E+01  0.00795 % r
   1   1   5 4.21435E+01  0.00769 % r
   1   1   6 4.21787E+01  0.00759 % r
   1   1   7 4.18686E+01  0.00794 % r
   1   1   8 4.18626E+01  0.00777 % r
   1   1   9 4.17206E+01  0.00762 % r
   1   1  10 4.23959E+01  0.00782 % r
   1   1  11 4.14863E+01  0.00766 % r
   1   1  12 4.26079E+01  0.00779 % r
   1   1  13 4.21108E+01  0.00778 % r
   1   1  14 4.25141E+01  0.00757 % r
   1   1  15 4.17339E+01  0.00788 % r
   1   1  16 4.22874E+01  0.00776 % r
   1   1  17 4.19611E+01  0.00778 % r
   1   1  18 4.20085E+01  0.00785 % r
   1   1  19 4.14147E+01  0.00777 % r
   1   1  20 4.21897E+01  0.00769 % r
   1   1  21 4.20175E+01  0.00780 % r
   1   1  22 4.14564E+01  0.00783 % r
   1   1  23 4.19055E+01  0.00802 % r
   1   1  24 4.14765E+01  0.00752 % r
   1   1  25 4.18087E+01  0.00762 % r
   1   1  26 4.24093E+01  0.00761 % r
   1   1  27 4.16309E+01  0.00779 % r
   1   1  28 4.19220E+01  0.00782 % r
   1   1  29 4.12901E+01  0.00777 % r
   1   1  30 4.17415E+01  0.00783 % r
   1   1  31 4.12562E+01  0.00805 % r
   1   1  32 4.19782E+01  0.00783 % r
   1   1  33 4.20081E+01  0.00768 % r
   1   1  34 4.18174E+01  0.00765 % r
   1   1  35 4.25728E+01  0.00767 % r
   1   1  36 4.14836E+01  0.00752 % r
   1   1  37 4.25099E+01  0.00780 % r
   1   1  38 4.18581E+01  0.00773 % r
   1   1  39 4.21708E+01  0.00770 % r
   1   1  40 4.15971E+01  0.00761 % r
   1   1  41 4.19467E+01  0.00778 % r
   1   1  42 4.14431E+01  0.00773 % r
   1   1  43 4.18067E+01  0.00778 % r
   1   1  44 4.23944E+01  0.00787 % r
   1   1  45 4.18076E+01  0.00757 % r
   1   1  46 4.16291E+01  0.00795 % r
   1   1  47 4.13774E+01  0.00793 % r
   1   1  48 4.17911E+01  0.00772 % r
   1   1  49 4.16352E+01  0.00751 % r
   1   1  50 4.24449E+01  0.00748 % r
   1   1  51 4.17182E+01  0.00750 % r
   1   1  52 4.15228E+01  0.00791 % r
   1   1  53 4.16941E+01  0.00824 % r
   1   1  54 4.16936E+01  0.00791 % r
   1   1  55 4.16488E+01  0.00790 % r
   1   1  56 4.16332E+01  0.00758 % r
   1   1  57 4.17645E+01  0.00767 % r
   1   1  58 4.19158E+01  0.00760 % r
   1   1  59 4.13296E+01  0.00771 % r
   1   1  60 4.18536E+01  0.00759 % r
   1   1  61 4.12952E+01  0.00780 % r
   1   1  62 4.21883E+01  0.00751 % r
   1   1  63 4.14594E+01  0.00781 % r
   1   1  64 4.19636E+01  0.00799 % r
   1   1  65 4.16835E+01  0.00756 % r
   1   1  66 4.17914E+01  0.00756 % r
   1   1  67 4.15777E+01  0.00777 % r
   1   1  68 4.21154E+01  0.00768 % r
   1   1  69 4.18054E+01  0.00741 % r
   1   1  70 4.19839E+01  0.00770 % r
   1   1  71 4.16147E+01  0.00807 % r
   1   1  72 4.21618E+01  0.00781 % r
   1   1  73 4.18884E+01  0.00773 % r
   1   1  74 4.14409E+01  0.00786 % r
   1   1  75 4.22875E+01  0.00739 % r
   1   1  76 4.24379E+01  0.00796 % r
   1   1  77 4.15298E+01  0.00766 % r
   1   1  78 4.14045E+01  0.00771 % r
   1   1  79 4.16623E+01  0.00814 % r
   1   1  80 4.18231E+01  0.00801 % r
   1   1  81 4.20321E+01  0.00769 % r
   1   1  82 4.11366E+01  0.00790 % r
   1   1  83 4.17560E+01  0.00800 % r
   1   1  84 4.16943E+01  0.00761 % r
   1   1  85 4.16636E+01  0.00756 % r
   1   1  86 4.19820E+01  0.00796 % r
   1   1  87 4.14670E+01  0.00792 % r
   1   1  88 4.18545E+01  0.00774 % r
   1   1  89 4.15513E+01  0.00743 % r
   1   1  90 4.12277E+01  0.00777 % r
   1   1  91 4.20770E+01  0.00769 % r
   1   1  92 4.24475E+01  0.00758 % r
   1   1  93 4.18947E+01  0.00761 % r
   1   1  94 4.17352E+01  0.00728 % r
   1   1  95 4.19022E+01  0.00762 % r
   1   1  96 4.12660E+01  0.00763 % r
   1   1  97 4.13522E+01  0.00778 % r
   1   1  98 4.17100E+01  0.00785 % r
   1   1  99 4.13702E+01  0.00765 % r
   1   1 100 4.24817E+01  0.00786 % r
   1   1 101 4.23687E+01  0.00755 % r
   1   1 102 4.15634E+01  0.00762 % r
   1   1 103 4.15200E+01  0.00808 % r
   1   1 104 4.19651E+01  0.00794 % r
   1   1 105 4.20251E+01  0.00804 % r
   1   1 106 4.17752E+01  0.00781 % r
   1   1 107 4.17397E+01  0.00767 % r
   1   1 108 4.15819E+01  0.00765 % r
   1   1 109 4.17386E+01  0.00736 % r
   1   1 110 4.24352E+01  0.00750 % r
   1   1 111 4.18662E+01  0.00764 % r
   1   1 112 4.21126E+01  0.00771 % r
   1   1 113 4.19519E+01  0.00776 % r
   1   1 114 4.20871E+01  0.00796 % r
   1   1 115 4.21011E+01  0.00774 % r
   1   1 116 4.20402E+01  0.00772 % r
   1   1 117 4.13132E+01  0.00762 % r
   1   1 118 4.13386E+01  0.00776 % r
   1   1 119 4.14585E+01  0.00761 % r
   1   1 120 4.15312E+01  0.00780 % r
   1   1 121 4.14686E+01  0.00798 % r
   1   1 122 4.16183E+01  0.00741 % r
   1   1 123 4.17677E+01  0.00787 % r
   1   1 124 4.18221E+01  0.00743 % r
   1   1 125 4.14885E+01  0.00808 % r
   1   1 126 4.18949E+01  0.00776 % r
   1   1 127 4.19414E+01  0.00770 % r
   1   1 128 4.15785E+01  0.00763 % r
   1   1 129 4.15202E+01  0.00776 % r
   1   1 130 4.17937E+01  0.00776 % r
   1   1 131 4.19238E+01  0.00752 % r
   1   1 132 4.16315E+01  0.00758 % r
   1   1 133 4.13724E+01  0.00779 % r
   1   1 134 4.17739E+01  0.00745 % r
   1   1 135 4.18645E+01  0.00786 % r
   1   1 136 4.09144E+01  0.00785 % r
   1   1 137 4.17325E+01  0.00788 % r
   1   1 138 4.18897E+01  0.00788 % r
   1   1 139 4.21760E+01  0.00782 % r
   1   1 140 4.14709E+01  0.00790 % r
   1   1 141 4.15321E+01  0.00743 % r
   1   1 142 4.11742E+01  0.00776 % r
   1   1 143 4.16211E+01  0.00769 % r
   1   1 144 4.11783E+01  0.00790 % r
   1   1 145 4.17908E+01  0.00805 % r
   1   1 146 4.11301E+01  0.00771 % r
   1   1 147 4.17832E+01  0.00770 % r
   1   1 148 4.10388E+01  0.00776 % r
   1   1 149 4.11442E+01  0.00760 % r
   1   1 150 4.15977E+01  0.00764 % r
   1   1 151 4.16817E+01  0.00800 % r
   1   1 152 4.12085E+01  0.00795 % r
   1   1 153 4.12661E+01  0.00760 % r
   1   1 154 4.17826E+01  0.00750 % r
   1   1 155 4.08805E+01  0.00796 % r
   1   1 156 4.12771E+01  0.00795 % r
   1   1 157 4.14629E+01  0.00789 % r
   1   1 158 4.12289E+01  0.00752 % r
   1   1 159 4.16310E+01  0.00759 % r
   1   1 160 4.11957E+01  0.00788 % r
   1   1 161 4.11406E+01  0.00785 % r
   1   1 162 4.16765E+01  0.00778 % r
   1   1 163 4.13636E+01  0.00760 % r
   1   1 164 4.13824E+01  0.00785 % r
   1   1 165 4.11074E+01  0.00779 % r
   1   1 166 4.08055E+01  0.00799 % r
   1   1 167 4.06018E+01  0.00757 % r
   1   1 168 4.17566E+01  0.00757 % r
   1   1 169 4.14616E+01  0.00732 % r
   1   1 170 4.13543E+01  0.00798 % r
   1   1 171 4.15353E+01  0.00786 % r
   1   1 172 4.06799E+01  0.00767 % r
   1   1 173 4.20393E+01  0.00776 % r
   1   1 174 4.12829E+01  0.00785 % r
   1   1 175 4.14133E+01  0.00779 % r
   1   1 176 4.10625E+01  0.00760 % r
   1   1 177 4.15138E+01  0.00744 % r
   1   1 178 4.11411E+01  0.00782 % r
   1   1 179 4.16467E+01  0.00802 % r
   1   1 180 4.16180E+01  0.00781 % r
   1   1 181 4.12934E+01  0.00807 % r
   1   1 182 4.09514E+01  0.00812 % r
   1   1 183 4.10150E+01  0.00795 % r
   1   1 184 4.12090E+01  0.00783 % r
   1   1 185 4.17860E+01  0.00757 % r
   1   1 186 4.16792E+01  0.00762 % r
   1   1 187 4.12389E+01  0.00771 % r
   1   1 188 4.10331E+01  0.00745 % r
   1   1 189 4.12844E+01  0.00790 % r
   1   1 190 4.15974E+01  0.00766 % r
   1   1 191 4.08339E+01  0.00781 % r
   1   1 192 4.13985E+01  0.00805 % r
   1   1 193 4.09881E+01  0.00798 % r
   1   1 194 4.14773E+01  0.00751 % r
   1   1 195 4.12239E+01  0.00779 % r
   1   1 196 4.15903E+01  0.00788 % r
   1   1 197 4.13709E+01  0.00752 % r
   1   1 198 4.18615E+01  0.00778 % r
   1   1 199 4.08262E+01  0.00765 % r
   1   1 200 4.13919E+01  0.00786 % r
   2   1   1 3.98321E+01  0.00783 % r
   2   1   2 4.00725E+01  0.00782 % r
   2   1   3 4.02663E+01  0.00792 % r
   2   1   4 4.00962E+01  0.00807 % r
   2   1   5 4.04263E+01  0.00787 % r
   2   1   6 4.01721E+01  0.00793 % r


By the way, I used set pop 1000000 1000 200
And, I am using Serpent2.1.24
Yousef Alzaben
Ph.D. Student
Karlsruhe Institute of Technology (KIT)
Yousef1986
 
Posts: 73
Joined: Fri Sep 04, 2015 3:33 pm
Location: Karlsruhe, Germany

Re: Unexpected Behavior of Axial Power Distribution

Postby jesse.johns » Fri Feb 19, 2016 7:55 pm

Yousef,

If I recall correctly, the first index of the axial distribution is the pin number.

In the attached file, you only show pin 1 and partial of pin 2. If I add up all of the axial nodes, then I get a power equal to the power in lvl0 for pin 1.

It is interesting that the distribution is flat, however. Do you have leakage at the top and bottom of your assembly?

Can you put the axial output for pin 20?

Jesse
jesse.johns
 
Posts: 114
Joined: Wed Apr 20, 2011 5:19 am

Re: Unexpected Behavior of Axial Power Distribution

Postby Yousef1986 » Mon Feb 22, 2016 12:40 pm

Hi Jesse,

What you said is exactly correct, I misunderstood the indices. I thought the outer loop for printing power distribution corresponds to axial layers, but it appeared to be for pin numbers.

I do not have axial leakages since I am using a reflective boundary condition through all faces (i.e. set bc 2 2 2), which explains why I do have almost flat power distribution.

After I had understood the meaning of indices, pin 20 does not exist in the output file since it is a non-fissile rod (no power produced) and "set cpd" card only print those rods with power > 0

Thank you Jesse for clarifying this issue.
Yousef Alzaben
Ph.D. Student
Karlsruhe Institute of Technology (KIT)
Yousef1986
 
Posts: 73
Joined: Fri Sep 04, 2015 3:33 pm
Location: Karlsruhe, Germany


Return to Bug reports

Who is online

Users browsing this forum: No registered users and 1 guest