# Using high temperature cross-sections in OpenMC

**URL:** https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249
**Category:** User Support
**Created:** [November 29, 2016, 8:42pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249 "2016-11-29T20:42:59Z")
**Posts on this page:** 9
**Page:** 1

<div class="post-metadata">

### Author: ![Julien](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/julien/32/1082_2.png) [@Julien](https://openmc.discourse.group/u/Julien)
#### Post date: [November 29, 2016, 8:42pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/1 "2016-11-29T20:42:59Z")

</div>

Hi all,

I am trying to use cross-sections at 800 K. I am not sure I am using the right way to set the temperature though. In my Python input file (I am using the Python API), I just set the “temperature” variable in my openmc.material to 800 K. For example, graphite = openmc.Material(name=‘Graphite’, temperature = ‘800’).

The problem is that when OpenMC runs, it does not find the corresponding cross-sections for the nuclides. This is the output on the command line:

_| The OpenMC Monte Carlo Code_  
_Copyright | 2011-2016 Massachusetts Institute of Technology_  
_License | [http://openmc.readthedocs.io/en/latest/license.html](http://openmc.readthedocs.io/en/latest/license.html)_  
_Version | 0.8.0_  
_Git SHA1 | 81c3bb593d30296efc531c15aba4b268cee687cf_  
_Date/Time | 2016-11-29 15:22:36_

---

<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: [November 29, 2016, 9:04pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/2 "2016-11-29T21:04:20Z")

</div>

Hi Julien,

The NNDC data unfortunately only includes cross sections at 293 K. You might want to look into [using JEFF data](http://openmc.readthedocs.io/en/latest/usersguide/install.html#using-jeff-cross-sections-from-oecd-nea) instead which contains temperatures up to 1800 K.

Best,  
Paul

---

<div class="post-metadata">

### Author: ![Julien](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/julien/32/1082_2.png) [@Julien](https://openmc.discourse.group/u/Julien)
#### Post date: [November 30, 2016, 4:29am UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/3 "2016-11-30T04:29:30Z")

</div>

Hi Paul,

Thanks for your help. I installed the JEFF data. However, I am not sure how I am supposed to use JEFF library. If I set the cross-section path to ‘/home/julien/openmc/data/jeff-3.2-hdf5/cross\_sections.xml’, OpenMC can’t find high temperature cross sections. There is another folder with sub-folders of ACE library at different temperatures but I can’t find any ‘cross\_sections.xml’ file to point to in this folder.  
How am I suppose to point to the JEFF library path so that I can use high temperature cross sections?

Thanks a lot.

Julien

---

<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: [December 5, 2016, 1:55pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/4 "2016-12-05T13:55:10Z")

</div>

Hi Julien,

With your version of OpenMC you should be using the HDF5 libraries that were converted from ACE files (the openmc-get-jeff-data script does this for you), not the ACE files directly. What is the error message you get when you try to run using …/jeff-3.2-hdf5/cross\_sections.xml?

Best,  
Paul

---

<div class="post-metadata">

### Author: ![Julien](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/julien/32/1082_2.png) [@Julien](https://openmc.discourse.group/u/Julien)
#### Post date: [January 5, 2017, 8:27pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/5 "2017-01-05T20:27:21Z")

</div>

Hi Paul,

Apologies for replying so late.

When I point my cross section path to ‘/home/julien/openmc/data/  
jeff-3.2-hdf5/cross\_sections.xml’ OpenMC is unable to find the cross sections of the isotopes even if the hdf5 corresponding files actually exist. For example when I start a pretty simple run case with uranium dioxyde, when OpenMC reads the cross sections it can’t find the cross section for the first nuclide (and I modified the composition to see if the first nuclides was the problem but it appears that OpenMC can’t find the cross sections of any nuclides). The error message I have is:

_Reading settings XML file... Reading geometry XML file... Reading cross sections XML file... Reading materials XML file..._

_Reading U235 from /home/julien/openmc/data/jeff-3.2-hdf5/U235.h5_  
_ERROR: Nuclear data library does not contain cross sections for U235 at or near_  
_294 K._  
_ERROR STOP_

Now, I also tried to use the MCNP cross sections (I converted them in HDF5 using the script openmc-convert-mcnp71-data that I found on the develop branch on Github) and I have the exact same error. I also changed the temperature in Material to see if that would work with other temperature but that did not fix the error.

I would be very grateful if someone could help me fix this error, preferably for the MCNP data case since I am planning on using MCNP cross sections.

Thanks a lot.

Julien

---

<div class="post-metadata">

### Author: ![Julien](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/julien/32/1082_2.png) [@Julien](https://openmc.discourse.group/u/Julien)
#### Post date: [January 6, 2017, 4:53pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/6 "2017-01-06T16:53:10Z")

</div>

Update on the problem:

I checked the source code where the HDF5 data are read. I discovered that the temperatures read and provided by nuclide\_header.F90 were in micro Kelvin for some reason. Therefore, the variable temps\_available in nuclide\_header.F90 is in mirco K while the temperature desired (temp\_desired) is in Kelvin. Obviously, OpenMC thinks that my desired temperature does not match any available temperatures. I divided the elements of temps\_available by 1e6 and OpenMC can read the cross section but I encounter segfault later in the process. Does someone has an idea about how to deal with this problem?

Thanks a lot.

Julien

---

<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: [January 6, 2017, 7:51pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/7 "2017-01-06T19:51:22Z")

</div>

Hi Julien,

We recently switched the units of energy in OpenMC from MeV to eV. You need to make sure that the HDF5 data you are using was produced by the same version of OpenMC that you are running with. It appears that your data was produced by an older version of OpenMC (and hence is in MeV) while the executable you are running expects data in eV. I would suggest regenerating your HDF5 cross sections and see if that fixes the problem.

Best regards,  
Paul

---

<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: [January 6, 2017, 8:06pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/8 "2017-01-06T20:06:56Z")

</div>

The way that temperature is treated with the master branch is fundamentally different than how it is treated in the latest develop branch. With the master branch, everything was based on ACE data. To select data at a particular temperature, you would indicate a cross section “suffix”, e.g. 71c or 80c, rather than specifying a temperature directly. The error you are seeing about temperatures not matching does not exist on the master branch. The fact that the code is able to open the HDF5 nuclear data at all means that it must be a later version. So, you may want to double check that you have installed the version that you think you did.

---

<div class="post-metadata">

### Author: ![Julien](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/julien/32/1082_2.png) [@Julien](https://openmc.discourse.group/u/Julien)
#### Post date: [January 7, 2017, 8:44pm UTC](https://openmc.discourse.group/t/using-high-temperature-cross-sections-in-openmc/249/9 "2017-01-07T20:44:55Z")

</div>

Paul,

You were right. For some reason I mixed up the versions provoking the errors I posted. I started a new develop-version OpenMC from scratch, converted the MCNP data and I can now use hight temperature cross sections. I am really sorry for the confusion. Thanks a lot anyway.

Julien
