Normalizing Heating Tally Data [eV/src-particle] to Power Density [W/cm^3]

User Support
1

Hello Everyone,

I am currently collecting heating tallies in a cylindrical mesh and I am interested in normalizing these heating tallies [eV/src-particle] to power density [W/cm^3].

I have already looked at the 8.3 Section, “Normalization of Tally Results” but I am having trouble getting to power density [W/cm^3] from the normalization factor in that section.

The end goal is to determine the power peaking factor in a localized region of the core. Are there other ways of normalization that will allow me to determine the power peaking factor?

Any help would be greatly appreciated

Regards, Matthew

2

Its trivial really. I wish you the best of luck neutronics lead <3

1 Like
3

k_tally = (1.60218e-19)*k_tally # convert ev/source to joules/source
k_tally_summed = np.sum(k_tally)

norm = power_level/k_tally_summed # source/s normalization factor

vol_power = normk_tally_summed/(4040fuel_lengthnp.pi*fuel_radius**2) # units of joules/s/cm3 or watts/cm3 for each rod

Here k_tally is the kappa_fission tally which starts as a vector of values where each value represents on fuel cell. I convert it to joules from eV to start. k_tally is then summed to get the total joules/source-particle for the core. norm is the normalization factor described by the 8.3 section where power level is in watts and is provided by the user. The volumetric power is then the summed kappa-fission tally times the normalization factor divided by the total fissionable volume over the entire core (in my case a 40x40 lattice of cylinders).

1 Like
4

@Daedalus, when I am multiplying norm with kappa_fission, aren’t I just calculating the value of norm cancelling out kappa_fission out of the equation?

5 
linear_power = k_tally*norm/(fuel_length/num_segs)   # joules/s/segment_length   (W/cm)

In the volumetric power check yes. The way I did it for individual rod linear powers was like this where you use the vector k_tally instead of the summed version to get a vector of linear powers. k_tally_summed isn’t cancelled in this method.

1 Like