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

Added JENDL nuclear data #3

Merged
merged 13 commits into from
Feb 23, 2019
Merged

Added JENDL nuclear data #3

merged 13 commits into from
Feb 23, 2019

Conversation

shimwell
Copy link
Member

Hi all

Perhaps this is not quite ready to merge, I wanted to share this code for downloading and converting JENDL 4.0 nuclear data and ask a question about metastables elements.

The other convert scripts have a section where the script opens the ACE file and changes the ZAID number to ensure it is identified as metastable.

I wanted to know if this is necessary when converting an NJOY file ? It appears that the metastables are identified when the h5 files are produced so I guess not. I just wanted to double check because this is the main difference between this script and the others.

The other difference is the use of openmc.data.IncidentNeutron.from_njoy(filename) which is slower than openmc.data.IncidentNeutron.from_ace(filename)

Best,

Jon

Copy link
Contributor

@paulromano paulromano left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for contributing this @shimwell. I have a few requests for changes if you don't mind. Also, I'd prefer the name of this script be generate_jendl.py. My naming convention is:

  • convert_*.py = Convert existing ACE files
  • generate_*.py = Generate ACE files using NJOY and then convert

convert_jendl.py Outdated Show resolved Hide resolved
convert_jendl.py Outdated Show resolved Hide resolved
convert_jendl.py Outdated Show resolved Hide resolved
convert_jendl.py Outdated Show resolved Hide resolved
@paulromano
Copy link
Contributor

One more comment I forgot to put before. I know you fusion guys don't care, but it would be good to include the JENDL thermal scattering sublibrary too. I might make some updates to this script though once merged, so feel free to ignore that and I can take care of it.

@shimwell
Copy link
Member Author

I also spotted a few convert scripts where we are downloading files without using the handy download function in the openmc _utils.py script. So I was able to remove a decent amount of lines in these scripts.

@shimwell
Copy link
Member Author

shimwell commented Feb 23, 2019

At the moment I am running all the download scripts to check that they work but it takes a while with the njoy parts. I guess this repo is not well suited to CI testing due to the large amounts of data downloading and time required for processing. Perhaps a topic for another PR but in the future we could carry out some automated testing with the different libraries. Perhaps something like this ...


# Run all download scripts within a Circle CI container with OpenMC installed
python3 data/convert_fendl.py
python3 data/convert_jeff32.py
python3 data/convert_nndc71.py
python3 data/convert_tendl.py
python3 data/generate_endf71.py
python3 data/generate_jendl.py

# set OPENMC_CROSS_SECTIONS path to each library 
# run simple sphere model for each isotope in the library with a neutron spectra tally
# flag up inconsistencies with the tally results

@paulromano
Copy link
Contributor

Thanks for cleaning up some of the other scripts -- I had been meaning to get rid of the manual downloading!

@paulromano paulromano merged commit 9f29cd7 into openmc-dev:master Feb 23, 2019
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

Successfully merging this pull request may close these issues.

2 participants