How can we add the value of "source intensity’’ in setting in an OpenMC program?
The intensity of a source is given as the
strength attribute of the
Source class. For example, if you want a source with an intensity of 10:
src = openmc.Source(..., strength=10.0)
src = openmc.Source(...) src.strength = 10.0
You can think of the strength as a global tally multiplier.
Thanks a lot, i appreciate your help.
Is this value in source particles/s?
Yeah I think so cuz the source intensity means the emmision rate, so the unit is particle/s.
The units are actually arbitrary. Normally tallies are given “per source particle”. If you treat your source strength in source particles/sec, then this results in tallies with units of “per sec” since the tallies are just multiplied by the source strength. However, if you consider your source strength to be something else, for example source particles/year, then you would get tally results “per year”.
Hi @paulromano. My question is because I would like to understand how to use the strength attribute to convert the flux units. The documentation says:
“As described in Scores, all tally scores are normalized per source particle simulated. However, for analysis of a given system, we usually want tally scores in a more natural unit. For example, neutron flux is often reported in units of particles/cm^2-sec. For a fixed source simulation, it is usually straightforward to convert units if the source rate is known. For example, if the system being modeled includes a source that is emitting 10^4 neutrons per second, the tally results just need to be multiplied by 10^4. This can either be done manually or using the openmc.Source.strength attribute”.
I have always used the first option. I mean, I get the flux in ‘particle-cm/source particle’ (flux obtained with strength = 1 – default value), multiply this flux by the source intensity in ‘source particle/sec’ and divide by the volume (cm^3) corresponding to the tally result to obtain the flux in ‘particle/cm^2-sec’ (see post with the equation).
If I properly understood the statement in the documentation, in the case that I want to use the second approach, I have to do the following:
openmc.Source.strength = 10^4
which will give us the flux in particle-cm/sec. Then dividing the flux tally by the volume [cm^3], we obtain the units ‘particle/cm^2-sec’. In other words, the Source intensity variable in the mentioned equation is not required, is this correct?
Maria – yes, that’s correct. If you are running a fixed source simulation and you specify
source = openmc.Source() source.strength = 1e4 ...
Then all tallies will be multiplied by 1e4.
Thank you @paulromano.