# Difficulties obtaining a flux or current tally score

**URL:** <https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247>\
**Category:** User Support\
**Created:** [July 31, 2023, 10:52am UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247 "2023-07-31T10:52:00Z")\
**Posts on this page:** 10\
**Page:** 1

<div class="post-metadata">

**Author:** ![KSartakov](https://avatars.discourse-cdn.com/v4/letter/k/76d3ee/32.png) [@KSartakov](https://openmc.discourse.group/u/KSartakov)\
**Post date:** [July 31, 2023, 10:52am UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/1 "2023-07-31T10:52:00Z")

</div>

Hello

I’m running OpenMC on a human phantom geometry that consists of an outer volume and multiple inner volumes. I will add a picture for reference.

 ![output](https://global.discourse-cdn.com/free1/uploads/openmc/original/2X/8/83f7c06afe52f8b86fe306409f7fa55523e1772e.png)

I’m trying to set the model up so that I can place the source in any inner volume (for example the brain) and have the model calculate the flux/current for any other chosen volume (for example the heart). But I’m running into a problem where I only get a numbered tally result for the most outer volume (skin) and any other volume inside gives me the tally result of 0. I have tried both with SurfaceFilter and MaterialFilter but the result is the same. I also tried adding SurfaceFromFilter to counter this result but with that added even the most outer volume gives me the 0 result.

At this point I can’t seem to figure out where this error could be resulting from and how to move on from this. Any help and insight would be greatly appreciated.

---

<div class="post-metadata">

**Author:** ![yrrepy](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/yrrepy/32/2365_2.png) [@yrrepy](https://openmc.discourse.group/u/yrrepy)\
**Post date:** [August 1, 2023, 9:47pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/2 "2023-08-01T21:47:15Z")

</div>

Hi Krislin,

You could overlay a mesh tally and confirm that particles are transporting within the voxel geometry.

I suspect there’s an issue with how the tallies are defined, perhaps too many filters.  
When I tally voxels I typically use a single particle filter, and a single _distribcell_ filter.  
Can you provide your XMLs?

Perry

---

<div class="post-metadata">

**Author:** ![KSartakov](https://avatars.discourse-cdn.com/v4/letter/k/76d3ee/32.png) [@KSartakov](https://openmc.discourse.group/u/KSartakov)\
**Post date:** [August 2, 2023, 2:36pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/3 "2023-08-02T14:36:37Z")

</div>

Hi Perry,

thank you for the response.  
So far I have only been using the surface filter, photon particle filter and energy function filter in my tallies. Occasionally adding cell from filter as well for testing purposes.

Here are my XML files:  
[geometry.xml](https://openmc.discourse.group/uploads/short-url/jpU98qFn5tsGGrqK0UC4JMXfX4v.xml) (719 Bytes)  
[materials.xml](https://openmc.discourse.group/uploads/short-url/yBRbXPY1vZTcgTOtYSMlxoOCwms.xml) (69.3 KB)  
[settings.xml](https://openmc.discourse.group/uploads/short-url/mJ7gr0R6wPR7W9SsX73FTIiG0s6.xml) (366 Bytes)  
[tallies.xml](https://openmc.discourse.group/uploads/short-url/4bklTAOjVHSAC9DLHW5syLzRrPm.xml) (1.3 KB)

Krislin

---

<div class="post-metadata">

**Author:** ![yrrepy](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/yrrepy/32/2365_2.png) [@yrrepy](https://openmc.discourse.group/u/yrrepy)\
**Post date:** [August 2, 2023, 10:20pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/4 "2023-08-02T22:20:49Z")

</div>

Hmm, I’m not familiar with CAD/mesh geometry.

I would try:

- score flux
- photon particle filter
- material filter, but only one material filter specifying one material, per tally  
(have individual tallies with individual filters for each material)

If this works, then add the energy-function filter.

---

<div class="post-metadata">

**Author:** ![KSartakov](https://avatars.discourse-cdn.com/v4/letter/k/76d3ee/32.png) [@KSartakov](https://openmc.discourse.group/u/KSartakov)\
**Post date:** [August 4, 2023, 2:54pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/5 "2023-08-04T14:54:52Z")

</div>

Thank you for the suggestion but using material filter and removing energy function filter still gives the same results- 0 flux tally for inner materials and numbered result for the outer material.

---

<div class="post-metadata">

**Author:** ![Daedalus](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/daedalus/32/2928_2.png) [@Daedalus](https://openmc.discourse.group/u/Daedalus)\
**Post date:** [August 4, 2023, 11:39pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/6 "2023-08-04T23:39:47Z")

</div>

Make sure it isn’t a tally sensitiity issue and make your source really hot to see if it returns anything. I’ve had that issue in the past.

---

<div class="post-metadata">

**Author:** ![KSartakov](https://avatars.discourse-cdn.com/v4/letter/k/76d3ee/32.png) [@KSartakov](https://openmc.discourse.group/u/KSartakov)\
**Post date:** [August 9, 2023, 2:08pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/7 "2023-08-09T14:08:21Z")

</div>

Increasing the source energy doesn’t bring any changes. But how would I go about checking the flux tally sensitivity?

---

<div class="post-metadata">

**Author:** ![Daedalus](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/daedalus/32/2928_2.png) [@Daedalus](https://openmc.discourse.group/u/Daedalus)\
**Post date:** [August 9, 2023, 2:22pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/8 "2023-08-09T14:22:43Z")

</div>

Check openmc.trigger. There are ways to adjust the tally threshhold (if that is the problem).

---

<div class="post-metadata">

**Author:** ![KSartakov](https://avatars.discourse-cdn.com/v4/letter/k/76d3ee/32.png) [@KSartakov](https://openmc.discourse.group/u/KSartakov)\
**Post date:** [August 10, 2023, 8:24pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/9 "2023-08-10T20:24:16Z")

</div>

It seems like tally sensitivity isn’t the problem here. I have been playing around with different energy levels and triggers but the results remain the same. Is there any way to conduct some sort of ‘diagnostics’? Without any error messages I feel like I’m just stumbling around in the dark.

---

<div class="post-metadata">

**Author:** ![Daedalus](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/daedalus/32/2928_2.png) [@Daedalus](https://openmc.discourse.group/u/Daedalus)\
**Post date:** [August 10, 2023, 9:11pm UTC](https://openmc.discourse.group/t/difficulties-obtaining-a-flux-or-current-tally-score/3247/10 "2023-08-10T21:11:56Z")

</div>

I have had a similar problem in the past twice with reactor simulations, but both times I was using local geometry, not a dagmc import. The first time I had my regions defined wrong and it ended up overriding all the cells into water (even though the voxel slice showed the correct geometry), which sounds similar to your issue where you are only getting skin registrations. If there is a way to check whether all of the source particles are interacting with the skin when that shouldn’t be the case, that might be something worth checking. The second time I had the issue I was generating multi-group cross sections for a full reactor simulation, but I wasn’t running enough particles (not particle energy, particle count) for certain pieces of aluminum to register any values. These are the only possible causes that come to mind, as I have no experience with imported geometries. Hope this helps.
