Power distribution HTR prismatic fuel element

Hello everyone,

I’m currently running into problems to get the power distribution on a HTR prismatic fuel element with TRISO explicitly modeled (in a random distribution not a lattice inside the fuel compact).

To obtain the power distribution, i built a fission tally with a DistribCellFilter with the fuel compact cell as an argument, it runs and all and seems to give good results since in the dataframe i have coherent numbers. Sadly, the openmc.lib.find_cell method seems to return me a cell filled with void instead of the fuel compact cells or the cells of my model so the image is blank.

Another way i imagined was to make a very fine MeshFilter (approx 4000*4000) on a fission tally, but when i try to run the model it crashes after 50 cycles (i have 50 inactive cycles).

Could some of you help me obtain the power distribution of my prismatic fuel element via one of those methods of via another one i didn’t think of ?

Thank in advance,
Sincerely,
Pierre