# Depletion calculations do not deplete **URL:** \ **Category:** User Support\ **Created:** [January 29, 2024, 10:52am UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873 "2024-01-29T10:52:34Z")\ **Posts on this page:** 12\ **Page:** 1
**Author:** ![Kladdy](https://avatars.discourse-cdn.com/v4/letter/k/839c29/32.png) [@Kladdy](https://openmc.discourse.group/u/Kladdy)\ **Post date:** [January 29, 2024, 10:52am UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/1 "2024-01-29T10:52:34Z")
I installed OpenMC 0.14.0 and downloaded the depletion example jupyter notebook from [https://nbviewer.org/github/openmc-dev/openmc-notebooks/blob/main/depletion.ipynb](https://nbviewer.org/github/openmc-dev/openmc-notebooks/blob/main/depletion.ipynb). When I run the calculations, however, I do not get anywhere near the same results as in the notebook. For me, the reactivity just jumps around the initial value within the uncertainties: ![output](https://global.discourse-cdn.com/free1/uploads/openmc/original/2X/3/32834a947b24736e3f71d07cfdc0951eb3527a21.png) While in the example linked above, the reactivity clearly goes down with time. This means that there actually is no depletion going on in my case, while in the example some form of depletion is absolutely happening. The only difference between my jupyter notebook and the one above is the part where the rectangular prism is defined, as this was changed in 0.14.0: ```python bound_box = openmc.model.RectangularPrism(1.24, 1.24, boundary_type="reflective") root_cell = openmc.Cell(fill=pin_univ, region=-bound_box) ``` I have built OpenMC from source on a M1 MacBook Pro, using OpenMPI. The CMake command I used was ```bash HDF5_ROOT=/opt/homebrew CXX=/opt/homebrew CC=/opt/homebrew/opt/llvm/bin/clang CXX=/opt/homebrew/opt/llvm/bin/clang++ LDFLAGS=-L/usr/local/lib cmake -DCMAKE_INSTALL_PREFIX=$HOME/openmc-mpi -DOPENMC_USE_MPI=ON -DOPENMC_USE_OPENMP=OFF .. ``` One weird thing is that even though i explicitly built OpenMC without OpenMP (I use MPI with OpenMPI instead, see command above), I get the following output when OpenMC is ran which seems to indicate that I am using OpenMP: ```txt | The OpenMC Monte Carlo Code Copyright | 2011-2023 MIT, UChicago Argonne LLC, and contributors License | https://docs.openmc.org/en/latest/license.html Version | 0.14.0 Git SHA1 | fa2330103de61a864c958d1a7250f11e5dd91468 Date/Time | 2024-01-29 11:44:49 MPI Processes | 1 OpenMP Threads | 10 Reading settings XML file... Reading cross sections XML file... Reading materials XML file... Reading geometry XML file... ``` Does anyone have a clue why my depletion is not running? I did run the OpenMC test suite, and It seems like most of the tests involving depletion has failed. I use the _endfb-vii.1-hdf5_ cross section library from [Official Data Libraries | OpenMC](https://openmc.org/official-data-libraries/) with the simplified depletion chain from [Depletion Chains | OpenMC](https://openmc.org/depletion-chains/). ---
**Author:** ![chiara.genoni](https://avatars.discourse-cdn.com/v4/letter/c/ec9cab/32.png) [@chiara.genoni](https://openmc.discourse.group/u/chiara.genoni)\ **Post date:** [January 29, 2024, 7:41pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/2 "2024-01-29T19:41:34Z")
I am getting the same issue when I run the depletion example jupyter notebook. The output plot is exactly the same, there is not any depletion occurring. ---
**Author:** ![Kladdy](https://avatars.discourse-cdn.com/v4/letter/k/839c29/32.png) [@Kladdy](https://openmc.discourse.group/u/Kladdy)\ **Post date:** [January 30, 2024, 7:45pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/4 "2024-01-30T19:45:21Z")
Interesting, so it seems like this is not an isolated incident. Which version of OpenMC are you running, and which cross section library/depletion chain are you using? ---
**Author:** ![chiara.genoni](https://avatars.discourse-cdn.com/v4/letter/c/ec9cab/32.png) [@chiara.genoni](https://openmc.discourse.group/u/chiara.genoni)\ **Post date:** [January 30, 2024, 8:07pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/5 "2024-01-30T20:07:54Z")
I am using version 14.0 of OpenMC, chain\_casl\_pwr.xml as depletion chain, endfb80 as cross section library. The plot I got is very similar: ![image](https://global.discourse-cdn.com/free1/uploads/openmc/original/2X/1/14346c97644b9422a7ac6ebf3e9f852d4940f236.png) ---
**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 30, 2024, 10:09pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/6 "2024-01-30T22:09:03Z")
One thing I’ll mention – there was recently a new version of SciPy released (1.12) and it caused a lot of our depletion tests to fail, and I have a sneaking suspicion that this may be the cause of the problems you all are seeing here. I haven’t had time to dig into the root cause yet and find a fix, but in the meantime I would suggest downgrading scipy to see if that helps (`pip install 'scipy==1.11.4'` or `conda install 'scipy=1.11.4'`). ---
**Author:** ![chiara.genoni](https://avatars.discourse-cdn.com/v4/letter/c/ec9cab/32.png) [@chiara.genoni](https://openmc.discourse.group/u/chiara.genoni)\ **Post date:** [January 31, 2024, 4:54pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/7 "2024-01-31T16:54:35Z")
Okay now it seems to work. Results are much more similar to the ones obtained in the example, the difference could be due to the different cross section library, but at least there is depletion. ![image](https://global.discourse-cdn.com/free1/uploads/openmc/original/2X/b/b59df3a76f544200f1c4dbf236109b47555b92f8.png) Thank you for your help!!! ---
**Author:** ![Kladdy](https://avatars.discourse-cdn.com/v4/letter/k/839c29/32.png) [@Kladdy](https://openmc.discourse.group/u/Kladdy)\ **Post date:** [February 4, 2024, 8:19pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/8 "2024-02-04T20:19:06Z")
This worked for me as well. Thanks a lot Paul! 🙂 ---
**Author:** ![chiara.genoni](https://avatars.discourse-cdn.com/v4/letter/c/ec9cab/32.png) [@chiara.genoni](https://openmc.discourse.group/u/chiara.genoni)\ **Post date:** [February 5, 2024, 4:31pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/9 "2024-02-05T16:31:07Z")
Can I ask what type of cross section library you used? Because now there is depletion but I feel like it is even too much compared to the results I should obtain and I do not think that difference could only be due to the cross section library. ---
**Author:** ![Kladdy](https://avatars.discourse-cdn.com/v4/letter/k/839c29/32.png) [@Kladdy](https://openmc.discourse.group/u/Kladdy)\ **Post date:** [February 7, 2024, 9:40pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/11 "2024-02-07T21:40:06Z")
I get the same result as you - i.e. way too much depletion compared to the example linked in the original post. My graph looks very much like yours, and the reactivity goes down a lot more than the example. I use ENDF/B-VII.1 as the nuclear data library and the Simplified Chain (Thermal Spectrum) from [Depletion Chains | OpenMC](https://openmc.org/depletion-chains/). I also tested the ENDF/B-VII.1 Chain (Thermal Spectrum), but it still gives me way too much depletion. ![output1](https://global.discourse-cdn.com/free1/uploads/openmc/original/2X/d/d8fc27914cec7e447537249347048f7dae9f73c2.png) ---
**Author:** ![dongbing](https://avatars.discourse-cdn.com/v4/letter/d/a87d85/32.png) [@dongbing](https://openmc.discourse.group/u/dongbing)\ **Post date:** [March 20, 2024, 3:35am UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/12 "2024-03-20T03:35:53Z")
It really works, thank you! ---
**Author:** ![KUNG](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/kung/32/3725_2.png) [@KUNG](https://openmc.discourse.group/u/KUNG)\ **Post date:** [May 22, 2024, 9:43am UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/13 "2024-05-22T09:43:18Z")
This solution worked for me too. Thanks! In my specific case the mpiexec runs of a adapted pincell depletion example just did not deplete (kinf looked constant, no pu239 fission rate increase)… it looked like it would redo the n0 step over and over again. After downgrading concentrations are changing as expected… Cheers. ---
**Author:** ![thilagam](https://avatars.discourse-cdn.com/v4/letter/t/e274bd/32.png) [@thilagam](https://openmc.discourse.group/u/thilagam)\ **Post date:** [April 23, 2025, 12:35pm UTC](https://openmc.discourse.group/t/depletion-calculations-do-not-deplete/3873/14 "2025-04-23T12:35:31Z")
Dear Dr. Paul Romano, Thank you for your kind guidance on this topic. I encountered the same issue after a recent update. It has now been resolved by reverting back to SciPy 1.11.4. With best regards, L Thilagam