Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Temperature interpolation bug #1125

Closed
liangjg opened this issue Nov 26, 2018 · 1 comment
Closed

Temperature interpolation bug #1125

liangjg opened this issue Nov 26, 2018 · 1 comment

Comments

@liangjg
Copy link
Contributor

liangjg commented Nov 26, 2018

Users from NCSU reported that when using temperature interpolation method in OpenMC, keff shows abnormal behavior: for a pure u235 infinite case, the kinf drops to nearly zero at 610 K, with 600K and 900K library.
I found this bug does not show up in the last stable version "v0.10.0" so it should be introduced in recent updates. Preliminary debugging shows there are unexpected zeros in the "NU_FISSION" array within the lower temperature data (that's why the kinf is closer to the normal value when the target temperature is nearer to the higher temperature). I'll dig into it further. Just let me know if you have any idea about this @paulromano @openmc-dev/committers

@paulromano
Copy link
Contributor

Thanks for taking a preliminary look @liangjg. I'll take a look at the source and see if anything pops out at me. Part of the problem is that our test data currently does not have multiple temperatures, so temperature interpolation is not being tested. We should definitely remedy that so I'll look into setting up a new library for testing.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants