Pin power distribution for VVER assembly

User Support
1

Hi there, I am trying to plot pin power distribution for hexagonal assembly.
There is an example for square assembly [here] (Multigroup Mode Part II: MGXS Library Generation with OpenMC — OpenMC Documentation) but no example for hexagonal assembly.
Can anyone describe how to plot pin power distribution for hexagonal assembly in openmc?

2

@TOWHID Take a look at Distributed-Cell-Tally-Visualization for Hexagonal Lattice example.

Hope this will help you

3

Thank you for your suggestion @Pranto

4

Hi @Pranto, i tried according to your last suggestion and i was able to plot pin power distribution. But now i am having problem to plot pin power distribution, when i am using burnable absorber pin cell with normal pin cell in my assembly. As it seems for some reason it only showing the burnable absorber mixed one’s but not the normal pins. Can you kindly tell me what i have to do if i want to plot in this case.
i am adding my script here. So kindly take a look at it and let me know what i am doing wrong or what i have to do?
again test (copy).py (76.5 KB)

5 
# Normal fuel pin
fuel = openmc.Cell(cell_id=2, name='Fuel')
fuel.fill = uo2

# IFBA pin
fuel = openmc.Cell(cell_id=10, name='Fuel')
fuel.fill = iba

You have two variables with the same name. Use a different one.

1 Like
6

@Pranto thanks for the help, I missed this one. Now it’s working.
But i am having trouble normalizing two different tally.
here is the code line

"ce_fission_rates1 /= np.mean(ce_fission_rates1[ce_fission_rates1 > 0.])
ce_fission_rates2 /= np.mean(ce_fission_rates2[ce_fission_rates2 > 0.]) "

The problem is that the 2nd one is normalized using only the 2nd one’s data
but i want to normalize the 2nd one but using the mean value of 1st one.
Kindly tell me how should i do it?

7

@TOWHID You can use mean() method to calculate the mean of the first data set, ce_fission_rates1.mean()

8

thanks a lot for your help

9

Hello, Pranto, I have some questions about this example, that is, when distribcell is changed from fuel_cell to main_cell, of course, the following fuel_cell.id is also changed to main_cell.id, the drawing becomes blank, why is this? I hope you can answer my questions, thank you very much!