# Source intensity

**URL:** <https://openmc.discourse.group/t/source-intensity/1454>\
**Category:** User Support\
**Created:** [September 22, 2021, 12:55pm UTC](https://openmc.discourse.group/t/source-intensity/1454 "2021-09-22T12:55:56Z")\
**Posts on this page:** 11\
**Page:** 1

<div class="post-metadata">

**Author:** ![younes](https://avatars.discourse-cdn.com/v4/letter/y/d6d6ee/32.png) [@younes](https://openmc.discourse.group/u/younes)\
**Post date:** [September 22, 2021, 12:55pm UTC](https://openmc.discourse.group/t/source-intensity/1454/1 "2021-09-22T12:55:56Z")

</div>

Hi;  
How can we add the value of "source intensity’’ in setting in an OpenMC program?  
Thank you.

---

<div class="post-metadata">

**Author:** ![paulromano](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/paulromano/32/486_2.png) [@paulromano](https://openmc.discourse.group/u/paulromano)\
**Post date:** [September 22, 2021, 2:01pm UTC](https://openmc.discourse.group/t/source-intensity/1454/2 "2021-09-22T14:01:24Z")

</div>

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:

```auto
src = openmc.Source(..., strength=10.0)

```

or

```auto
src = openmc.Source(...)
src.strength = 10.0

```

You can think of the strength as a global tally multiplier.

---

<div class="post-metadata">

**Author:** ![younes](https://avatars.discourse-cdn.com/v4/letter/y/d6d6ee/32.png) [@younes](https://openmc.discourse.group/u/younes)\
**Post date:** [September 22, 2021, 4:19pm UTC](https://openmc.discourse.group/t/source-intensity/1454/3 "2021-09-22T16:19:38Z")

</div>

Hi Paul;  
Thanks a lot, i appreciate your help.  
Younes.

---

<div class="post-metadata">

**Author:** ![MaryAAZ](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/maryaaz/32/1244_2.png) [@MaryAAZ](https://openmc.discourse.group/u/MaryAAZ)\
**Post date:** [September 22, 2021, 6:48pm UTC](https://openmc.discourse.group/t/source-intensity/1454/4 "2021-09-22T18:48:36Z")

</div>

Hi @paulromano,

Is this value in source particles/s?

Maria

---

<div class="post-metadata">

**Author:** ![younes](https://avatars.discourse-cdn.com/v4/letter/y/d6d6ee/32.png) [@younes](https://openmc.discourse.group/u/younes)\
**Post date:** [September 22, 2021, 7:09pm UTC](https://openmc.discourse.group/t/source-intensity/1454/5 "2021-09-22T19:09:01Z")

</div>

Hi @Mary;  
Yeah I think so cuz the source intensity means the emmision rate, so the unit is particle/s.  
Younes.

---

<div class="post-metadata">

**Author:** ![paulromano](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/paulromano/32/486_2.png) [@paulromano](https://openmc.discourse.group/u/paulromano)\
**Post date:** [September 22, 2021, 8:55pm UTC](https://openmc.discourse.group/t/source-intensity/1454/6 "2021-09-22T20:55:43Z")

</div>

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”.

---

<div class="post-metadata">

**Author:** ![MaryAAZ](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/maryaaz/32/1244_2.png) [@MaryAAZ](https://openmc.discourse.group/u/MaryAAZ)\
**Post date:** [September 22, 2021, 10:57pm UTC](https://openmc.discourse.group/t/source-intensity/1454/7 "2021-09-22T22:57:47Z")

</div>

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](https://docs.openmc.org/en/stable/usersguide/tallies.html#usersguide-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](https://openmc.discourse.group/t/unit-transformation-in-a-fixed-source-problem/1135) 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?

Thanks,  
Maria

---

<div class="post-metadata">

**Author:** ![paulromano](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/paulromano/32/486_2.png) [@paulromano](https://openmc.discourse.group/u/paulromano)\
**Post date:** [September 29, 2021, 3:24am UTC](https://openmc.discourse.group/t/source-intensity/1454/8 "2021-09-29T03:24:54Z")

</div>

Maria – yes, that’s correct. If you are running a fixed source simulation and you specify

```auto
source = openmc.Source()
source.strength = 1e4
...

```

Then all tallies will be multiplied by 1e4.

---

<div class="post-metadata">

**Author:** ![MaryAAZ](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/maryaaz/32/1244_2.png) [@MaryAAZ](https://openmc.discourse.group/u/MaryAAZ)\
**Post date:** [September 29, 2021, 7:02pm UTC](https://openmc.discourse.group/t/source-intensity/1454/9 "2021-09-29T19:02:03Z")

</div>

Thank you @paulromano.

---

<div class="post-metadata">

**Author:** ![kevinm387](https://avatars.discourse-cdn.com/v4/letter/k/bcef8e/32.png) [@kevinm387](https://openmc.discourse.group/u/kevinm387)\
**Post date:** [April 19, 2022, 11:30pm UTC](https://openmc.discourse.group/t/source-intensity/1454/10 "2022-04-19T23:30:24Z")

</div>

Hi all,

Just want to ask a follow up question about how this applies to problems with multiple sources. For instance, if you had two identical sources, each with `source_strength = 1`, would you get tallies that are multiplied by two? Or would OpenMC do some sort of normalization? And if there is normalization, how would that extend to problems with many sources of varying strengths?

For a problem with multiple sources, would I essentially need to have all my source\_strengths add up to a total of 1 to get tallies comparable to a problem with only a single source with `source_strength = 1`?

---

<div class="post-metadata">

**Author:** ![MaryAAZ](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/maryaaz/32/1244_2.png) [@MaryAAZ](https://openmc.discourse.group/u/MaryAAZ)\
**Post date:** [April 20, 2022, 3:18am UTC](https://openmc.discourse.group/t/source-intensity/1454/11 "2022-04-20T03:18:07Z")

</div>

Hello,

I think the post below can help you with your question.

> [@Photon flux with NaN value](https://openmc.discourse.group/t/photon-flux-with-nan-value/1122):
>
> Hello all, I am having an issue in a shielding simulation. My model consists of different cylinders, each one filled with a different material (see figure). In the center, I have positioned two point sources, one neutron source and the other is a photon source. As I have both energy spectra, I am using openmc.stats.Tabular to inform them (following the recommendations in this [post](https://openmc.discourse.group/t/source-energy-distribution-for-neutron-of-different-energy-group/477)). I wish to tally the fluxes, then I created a mesh (101 x 101 elements) and particle filters. [image] The mesh h…

Maria
