Skip to content

Instantly share code, notes, and snippets.

@paulromano

paulromano/Depletion changing material.ipynb

Created January 9, 2022 19:45
Show Gist options
  • Save paulromano/3fe537f80cc974661d766e74ed7fe8ee to your computer and use it in GitHub Desktop.
Save paulromano/3fe537f80cc974661d766e74ed7fe8ee to your computer and use it in GitHub Desktop.
Download ZIP
Example of running a depletion calculation while changing a material in the middle
Display the source blob
Display the rendered blob
Raw
Depletion changing material.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "4b368e3d",
"metadata": {},
"outputs": [],
"source": [
"import openmc.deplete\n",
"from math import pi\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"id": "f1a911ba",
"metadata": {},
"source": [
"First create a simple one material, one region problem and deplete it over three 1-day timesteps."
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "83e7599d",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" %%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" ############### %%%%%%%%%%%%%%%%%%%%%%%%\n",
" ################## %%%%%%%%%%%%%%%%%%%%%%%\n",
" ################### %%%%%%%%%%%%%%%%%%%%%%%\n",
" #################### %%%%%%%%%%%%%%%%%%%%%%\n",
" ##################### %%%%%%%%%%%%%%%%%%%%%\n",
" ###################### %%%%%%%%%%%%%%%%%%%%\n",
" ####################### %%%%%%%%%%%%%%%%%%\n",
" ####################### %%%%%%%%%%%%%%%%%\n",
" ###################### %%%%%%%%%%%%%%%%%\n",
" #################### %%%%%%%%%%%%%%%%%\n",
" ################# %%%%%%%%%%%%%%%%%\n",
" ############### %%%%%%%%%%%%%%%%\n",
" ############ %%%%%%%%%%%%%%%\n",
" ######## %%%%%%%%%%%%%%\n",
" %%%%%%%%%%%\n",
"\n",
" | The OpenMC Monte Carlo Code\n",
" Copyright | 2011-2021 MIT, UChicago Argonne LLC, and contributors\n",
" License | https://docs.openmc.org/en/latest/license.html\n",
" Version | 0.13.0-dev\n",
" Git SHA1 | 65b02c57b9c09cf0c4c914a31034ed076688165e\n",
" Date/Time | 2022-01-09 13:43:52\n",
" OpenMP Threads | 12\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading geometry XML file...\n",
" Reading U235 from /opt/data/hdf5/nndc_hdf5_v15/U235.h5\n",
" Minimum neutron data temperature: 294 K\n",
" Maximum neutron data temperature: 294 K\n",
" Preparing distributed cell instances...\n",
" Reading plot XML file...\n",
" Writing summary.h5 file...\n",
" Reading I135 from /opt/data/hdf5/nndc_hdf5_v15/I135.h5\n",
" Reading Xe135 from /opt/data/hdf5/nndc_hdf5_v15/Xe135.h5\n",
" Reading Xe136 from /opt/data/hdf5/nndc_hdf5_v15/Xe136.h5\n",
" Reading Cs135 from /opt/data/hdf5/nndc_hdf5_v15/Cs135.h5\n",
" Reading Gd157 from /opt/data/hdf5/nndc_hdf5_v15/Gd157.h5\n",
" Reading Gd156 from /opt/data/hdf5/nndc_hdf5_v15/Gd156.h5\n",
" Reading U234 from /opt/data/hdf5/nndc_hdf5_v15/U234.h5\n",
" Reading U238 from /opt/data/hdf5/nndc_hdf5_v15/U238.h5\n",
" Maximum neutron transport energy: 20000000 eV for U235\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.44518\n",
" 2/1 0.37906\n",
" 3/1 0.35856\n",
" 4/1 0.33546\n",
" 5/1 0.33864\n",
" 6/1 0.34190\n",
" 7/1 0.32486\n",
" 8/1 0.35235\n",
" 9/1 0.33251\n",
" 10/1 0.33518\n",
" 11/1 0.32198\n",
" 12/1 0.32571 0.32384 +/- 0.00187\n",
" 13/1 0.33900 0.32889 +/- 0.00517\n",
" 14/1 0.32974 0.32911 +/- 0.00366\n",
" 15/1 0.31343 0.32597 +/- 0.00423\n",
" 16/1 0.34141 0.32854 +/- 0.00430\n",
" 17/1 0.33988 0.33016 +/- 0.00398\n",
" 18/1 0.32317 0.32929 +/- 0.00356\n",
" 19/1 0.33678 0.33012 +/- 0.00325\n",
" 20/1 0.33961 0.33107 +/- 0.00305\n",
" 21/1 0.32598 0.33061 +/- 0.00280\n",
" 22/1 0.34779 0.33204 +/- 0.00293\n",
" 23/1 0.34242 0.33284 +/- 0.00281\n",
" 24/1 0.34736 0.33387 +/- 0.00280\n",
" 25/1 0.32826 0.33350 +/- 0.00264\n",
" 26/1 0.33177 0.33339 +/- 0.00247\n",
" 27/1 0.33458 0.33346 +/- 0.00232\n",
" 28/1 0.33354 0.33347 +/- 0.00219\n",
" 29/1 0.33125 0.33335 +/- 0.00207\n",
" 30/1 0.33253 0.33331 +/- 0.00197\n",
" 31/1 0.33983 0.33362 +/- 0.00190\n",
" 32/1 0.32817 0.33337 +/- 0.00182\n",
" 33/1 0.33153 0.33329 +/- 0.00174\n",
" 34/1 0.33852 0.33351 +/- 0.00168\n",
" 35/1 0.33025 0.33338 +/- 0.00162\n",
" 36/1 0.33850 0.33358 +/- 0.00157\n",
" 37/1 0.33694 0.33370 +/- 0.00152\n",
" 38/1 0.33765 0.33384 +/- 0.00147\n",
" 39/1 0.33440 0.33386 +/- 0.00142\n",
" 40/1 0.32965 0.33372 +/- 0.00138\n",
" 41/1 0.33707 0.33383 +/- 0.00133\n",
" 42/1 0.34793 0.33427 +/- 0.00137\n",
" 43/1 0.32330 0.33394 +/- 0.00136\n",
" 44/1 0.33243 0.33389 +/- 0.00132\n",
" 45/1 0.32191 0.33355 +/- 0.00133\n",
" 46/1 0.32493 0.33331 +/- 0.00132\n",
" 47/1 0.33701 0.33341 +/- 0.00128\n",
" 48/1 0.33396 0.33343 +/- 0.00125\n",
" 49/1 0.34328 0.33368 +/- 0.00124\n",
" 50/1 0.33067 0.33360 +/- 0.00121\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.7577e-02 seconds\n",
" Reading cross sections = 3.9985e-02 seconds\n",
" Total time in simulation = 4.8107e-02 seconds\n",
" Time in transport only = 3.7646e-02 seconds\n",
" Time in inactive batches = 8.3598e-03 seconds\n",
" Time in active batches = 3.9747e-02 seconds\n",
" Time synchronizing fission bank = 2.6843e-03 seconds\n",
" Sampling source sites = 2.4146e-03 seconds\n",
" SEND/RECV source sites = 2.5521e-04 seconds\n",
" Time accumulating tallies = 2.9472e-03 seconds\n",
" Time writing statepoints = 2.2034e-03 seconds\n",
" Total time for finalization = 5.9802e-05 seconds\n",
" Total time elapsed = 9.9555e-02 seconds\n",
" Calculation Rate (inactive) = 1.19621e+06 particles/second\n",
" Calculation Rate (active) = 1.00637e+06 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33654 +/- 0.00174\n",
" k-effective (Track-length) = 0.33360 +/- 0.00121\n",
" k-effective (Absorption) = 0.34157 +/- 0.00442\n",
" Combined k-effective = 0.33457 +/- 0.00111\n",
" Leakage Fraction = 0.85930 +/- 0.00182\n",
"\n",
" Creating state point openmc_simulation_n0.h5...\n",
" Maximum neutron transport energy: 20000000 eV for U235\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.44339\n",
" 2/1 0.37826\n",
" 3/1 0.35221\n",
" 4/1 0.34922\n",
" 5/1 0.34501\n",
" 6/1 0.33781\n",
" 7/1 0.33154\n",
" 8/1 0.33461\n",
" 9/1 0.33471\n",
" 10/1 0.33565\n",
" 11/1 0.34487\n",
" 12/1 0.33343 0.33915 +/- 0.00572\n",
" 13/1 0.33519 0.33783 +/- 0.00356\n",
" 14/1 0.33699 0.33762 +/- 0.00252\n",
" 15/1 0.33378 0.33685 +/- 0.00210\n",
" 16/1 0.33426 0.33642 +/- 0.00177\n",
" 17/1 0.33301 0.33593 +/- 0.00157\n",
" 18/1 0.34614 0.33721 +/- 0.00187\n",
" 19/1 0.32589 0.33595 +/- 0.00207\n",
" 20/1 0.32563 0.33492 +/- 0.00212\n",
" 21/1 0.33598 0.33502 +/- 0.00192\n",
" 22/1 0.34082 0.33550 +/- 0.00182\n",
" 23/1 0.33245 0.33527 +/- 0.00169\n",
" 24/1 0.32735 0.33470 +/- 0.00166\n",
" 25/1 0.32265 0.33390 +/- 0.00174\n",
" 26/1 0.33420 0.33392 +/- 0.00163\n",
" 27/1 0.34306 0.33445 +/- 0.00162\n",
" 28/1 0.33313 0.33438 +/- 0.00153\n",
" 29/1 0.34539 0.33496 +/- 0.00156\n",
" 30/1 0.33775 0.33510 +/- 0.00149\n",
" 31/1 0.32327 0.33454 +/- 0.00152\n",
" 32/1 0.33729 0.33466 +/- 0.00146\n",
" 33/1 0.33310 0.33459 +/- 0.00139\n",
" 34/1 0.32289 0.33411 +/- 0.00142\n",
" 35/1 0.32959 0.33392 +/- 0.00138\n",
" 36/1 0.34504 0.33435 +/- 0.00139\n",
" 37/1 0.31932 0.33380 +/- 0.00145\n",
" 38/1 0.31966 0.33329 +/- 0.00148\n",
" 39/1 0.32981 0.33317 +/- 0.00144\n",
" 40/1 0.33276 0.33316 +/- 0.00139\n",
" 41/1 0.33519 0.33322 +/- 0.00134\n",
" 42/1 0.33414 0.33325 +/- 0.00130\n",
" 43/1 0.34086 0.33348 +/- 0.00128\n",
" 44/1 0.33009 0.33338 +/- 0.00125\n",
" 45/1 0.32816 0.33323 +/- 0.00122\n",
" 46/1 0.32060 0.33288 +/- 0.00124\n",
" 47/1 0.34275 0.33315 +/- 0.00123\n",
" 48/1 0.33844 0.33329 +/- 0.00121\n",
" 49/1 0.32663 0.33312 +/- 0.00119\n",
" 50/1 0.33314 0.33312 +/- 0.00116\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 5.2577e-02 seconds\n",
" Time in transport only = 3.6513e-02 seconds\n",
" Time in inactive batches = 5.9433e-03 seconds\n",
" Time in active batches = 4.6634e-02 seconds\n",
" Time synchronizing fission bank = 2.8639e-03 seconds\n",
" Sampling source sites = 2.5966e-03 seconds\n",
" SEND/RECV source sites = 2.5190e-04 seconds\n",
" Time accumulating tallies = 6.8783e-03 seconds\n",
" Time writing statepoints = 5.3778e-03 seconds\n",
" Total time for finalization = 6.4615e-05 seconds\n",
" Total time elapsed = 5.5925e-02 seconds\n",
" Calculation Rate (inactive) = 1.68258e+06 particles/second\n",
" Calculation Rate (active) = 857749 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33492 +/- 0.00219\n",
" k-effective (Track-length) = 0.33312 +/- 0.00116\n",
" k-effective (Absorption) = 0.33368 +/- 0.00373\n",
" Combined k-effective = 0.33312 +/- 0.00115\n",
" Leakage Fraction = 0.86337 +/- 0.00156\n",
"\n",
" Creating state point openmc_simulation_n1.h5...\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" Maximum neutron transport energy: 20000000 eV for U235\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.43764\n",
" 2/1 0.37151\n",
" 3/1 0.35618\n",
" 4/1 0.34880\n",
" 5/1 0.33780\n",
" 6/1 0.33063\n",
" 7/1 0.34847\n",
" 8/1 0.32855\n",
" 9/1 0.32569\n",
" 10/1 0.33595\n",
" 11/1 0.33918\n",
" 12/1 0.32217 0.33067 +/- 0.00850\n",
" 13/1 0.32066 0.32734 +/- 0.00594\n",
" 14/1 0.34358 0.33140 +/- 0.00584\n",
" 15/1 0.32657 0.33043 +/- 0.00463\n",
" 16/1 0.34450 0.33278 +/- 0.00445\n",
" 17/1 0.34323 0.33427 +/- 0.00404\n",
" 18/1 0.33571 0.33445 +/- 0.00351\n",
" 19/1 0.33439 0.33444 +/- 0.00309\n",
" 20/1 0.33000 0.33400 +/- 0.00280\n",
" 21/1 0.32163 0.33287 +/- 0.00277\n",
" 22/1 0.32287 0.33204 +/- 0.00266\n",
" 23/1 0.33314 0.33213 +/- 0.00245\n",
" 24/1 0.33211 0.33212 +/- 0.00227\n",
" 25/1 0.33497 0.33231 +/- 0.00212\n",
" 26/1 0.33975 0.33278 +/- 0.00204\n",
" 27/1 0.33596 0.33297 +/- 0.00192\n",
" 28/1 0.35064 0.33395 +/- 0.00206\n",
" 29/1 0.33884 0.33421 +/- 0.00197\n",
" 30/1 0.33380 0.33418 +/- 0.00187\n",
" 31/1 0.34804 0.33484 +/- 0.00189\n",
" 32/1 0.32764 0.33452 +/- 0.00184\n",
" 33/1 0.33207 0.33441 +/- 0.00176\n",
" 34/1 0.32265 0.33392 +/- 0.00175\n",
" 35/1 0.34831 0.33450 +/- 0.00178\n",
" 36/1 0.34693 0.33497 +/- 0.00177\n",
" 37/1 0.32072 0.33445 +/- 0.00179\n",
" 38/1 0.34476 0.33481 +/- 0.00176\n",
" 39/1 0.33968 0.33498 +/- 0.00171\n",
" 40/1 0.33298 0.33492 +/- 0.00165\n",
" 41/1 0.34382 0.33520 +/- 0.00162\n",
" 42/1 0.32851 0.33499 +/- 0.00158\n",
" 43/1 0.31258 0.33431 +/- 0.00168\n",
" 44/1 0.33554 0.33435 +/- 0.00163\n",
" 45/1 0.32995 0.33423 +/- 0.00159\n",
" 46/1 0.31969 0.33382 +/- 0.00159\n",
" 47/1 0.32416 0.33356 +/- 0.00157\n",
" 48/1 0.32749 0.33340 +/- 0.00154\n",
" 49/1 0.31304 0.33288 +/- 0.00159\n",
" 50/1 0.32534 0.33269 +/- 0.00156\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 4.9787e-02 seconds\n",
" Time in transport only = 3.5364e-02 seconds\n",
" Time in inactive batches = 5.6371e-03 seconds\n",
" Time in active batches = 4.4150e-02 seconds\n",
" Time synchronizing fission bank = 2.7776e-03 seconds\n",
" Sampling source sites = 2.5122e-03 seconds\n",
" SEND/RECV source sites = 2.5014e-04 seconds\n",
" Time accumulating tallies = 5.2078e-03 seconds\n",
" Time writing statepoints = 5.8251e-03 seconds\n",
" Total time for finalization = 6.1818e-05 seconds\n",
" Total time elapsed = 5.3166e-02 seconds\n",
" Calculation Rate (inactive) = 1.77395e+06 particles/second\n",
" Calculation Rate (active) = 905996 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.32995 +/- 0.00282\n",
" k-effective (Track-length) = 0.33269 +/- 0.00156\n",
" k-effective (Absorption) = 0.32968 +/- 0.00409\n",
" Combined k-effective = 0.33273 +/- 0.00161\n",
" Leakage Fraction = 0.86487 +/- 0.00158\n",
"\n",
" Creating state point openmc_simulation_n2.h5...\n",
" Maximum neutron transport energy: 20000000 eV for U235\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.43897\n",
" 2/1 0.36752\n",
" 3/1 0.33475\n",
" 4/1 0.34311\n",
" 5/1 0.32238\n",
" 6/1 0.32624\n",
" 7/1 0.32126\n",
" 8/1 0.32441\n",
" 9/1 0.33653\n",
" 10/1 0.31893\n",
" 11/1 0.31984\n",
" 12/1 0.33756 0.32870 +/- 0.00886\n",
" 13/1 0.33054 0.32932 +/- 0.00515\n",
" 14/1 0.33394 0.33047 +/- 0.00382\n",
" 15/1 0.32740 0.32986 +/- 0.00302\n",
" 16/1 0.31232 0.32693 +/- 0.00383\n",
" 17/1 0.33812 0.32853 +/- 0.00361\n",
" 18/1 0.33259 0.32904 +/- 0.00316\n",
" 19/1 0.32548 0.32864 +/- 0.00282\n",
" 20/1 0.34703 0.33048 +/- 0.00312\n",
" 21/1 0.33728 0.33110 +/- 0.00289\n",
" 22/1 0.33859 0.33172 +/- 0.00271\n",
" 23/1 0.33579 0.33204 +/- 0.00251\n",
" 24/1 0.32914 0.33183 +/- 0.00234\n",
" 25/1 0.33862 0.33228 +/- 0.00222\n",
" 26/1 0.33775 0.33263 +/- 0.00211\n",
" 27/1 0.32328 0.33208 +/- 0.00205\n",
" 28/1 0.33470 0.33222 +/- 0.00194\n",
" 29/1 0.33050 0.33213 +/- 0.00184\n",
" 30/1 0.32773 0.33191 +/- 0.00176\n",
" 31/1 0.33071 0.33185 +/- 0.00167\n",
" 32/1 0.33949 0.33220 +/- 0.00163\n",
" 33/1 0.33229 0.33220 +/- 0.00156\n",
" 34/1 0.32804 0.33203 +/- 0.00150\n",
" 35/1 0.32284 0.33166 +/- 0.00149\n",
" 36/1 0.34618 0.33222 +/- 0.00153\n",
" 37/1 0.33959 0.33250 +/- 0.00150\n",
" 38/1 0.34136 0.33281 +/- 0.00148\n",
" 39/1 0.33620 0.33293 +/- 0.00143\n",
" 40/1 0.32989 0.33283 +/- 0.00139\n",
" 41/1 0.33808 0.33300 +/- 0.00135\n",
" 42/1 0.32021 0.33260 +/- 0.00137\n",
" 43/1 0.33591 0.33270 +/- 0.00133\n",
" 44/1 0.33040 0.33263 +/- 0.00129\n",
" 45/1 0.33325 0.33265 +/- 0.00126\n",
" 46/1 0.33018 0.33258 +/- 0.00122\n",
" 47/1 0.32765 0.33245 +/- 0.00120\n",
" 48/1 0.33027 0.33239 +/- 0.00117\n",
" 49/1 0.32830 0.33228 +/- 0.00114\n",
" 50/1 0.32343 0.33206 +/- 0.00113\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 5.4866e-02 seconds\n",
" Time in transport only = 4.0356e-02 seconds\n",
" Time in inactive batches = 7.1365e-03 seconds\n",
" Time in active batches = 4.7730e-02 seconds\n",
" Time synchronizing fission bank = 2.9047e-03 seconds\n",
" Sampling source sites = 2.6269e-03 seconds\n",
" SEND/RECV source sites = 2.6190e-04 seconds\n",
" Time accumulating tallies = 4.7966e-03 seconds\n",
" Time writing statepoints = 5.9634e-03 seconds\n",
" Total time for finalization = 6.4340e-05 seconds\n",
" Total time elapsed = 5.8373e-02 seconds\n",
" Calculation Rate (inactive) = 1.40124e+06 particles/second\n",
" Calculation Rate (active) = 838054 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33116 +/- 0.00215\n",
" k-effective (Track-length) = 0.33206 +/- 0.00113\n",
" k-effective (Absorption) = 0.33397 +/- 0.00448\n",
" Combined k-effective = 0.33231 +/- 0.00112\n",
" Leakage Fraction = 0.86295 +/- 0.00184\n",
"\n",
" Creating state point openmc_simulation_n3.h5...\n"
]
}
],
"source": [
"material = openmc.Material()\n",
"material.add_nuclide('U235', 1.0)\n",
"material.set_density('g/cm3', 5.0)\n",
"material.volume = 4/3 * pi * 10.0**3\n",
"\n",
"model = openmc.Model()\n",
"sph = openmc.Sphere(r=10.0, boundary_type='vacuum')\n",
"cell = openmc.Cell(fill=material, region=-sph)\n",
"model.geometry = openmc.Geometry([cell])\n",
"\n",
"model.settings = openmc.Settings()\n",
"model.settings.particles = 1000\n",
"model.settings.batches = 50\n",
"model.settings.inactive = 10\n",
"\n",
"chain_file = 'chain_simple.xml'\n",
"op = openmc.deplete.Operator(model, chain_file)\n",
"power = 40.0e3\n",
"time_steps = [1.0]*3\n",
"integrator = openmc.deplete.PredictorIntegrator(op, time_steps, timestep_units='d', power=power)\n",
"integrator.integrate()"
]
},
{
"cell_type": "markdown",
"id": "d1b6aab3",
"metadata": {},
"source": [
"Now let's look at the density of $^{135}$Xe from the results:"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "b07d2c6f",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Xe135 [atoms]')"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAERCAYAAABhKjCtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAsCUlEQVR4nO3deXxU5dn/8c/FEgj7krCvAgKKIBDBrYrWBbWVurX4gEtrHyoVl+fRWuvTXxerrVtt1YpCW6qWilatFRdA675rwiqbBARJEAhb2BLIcv3+mIMOYZIMkMmZTL7v12tezFnnOgyca86573Pd5u6IiIhU1CDsAEREJDkpQYiISExKECIiEpMShIiIxKQEISIiMSlBiIhITCmXIMxsmpltNLNP41j3FDOba2alZnZxhWV3m9liM1tqZg+YmSUuahGR5JNyCQJ4FBgd57pfAFcCT0TPNLMTgZOAwcAg4Djg1BqLUESkDki5BOHubwNboueZWR8zm21mOWb2jpkNCNZd7e4LgfKKuwGaAmlAE6AxsCHx0YuIJI+USxCVmApc6+7DgZuAyVWt7O4fAG8AXwavOe6+NOFRiogkkUZhB5BoZtYCOBF4OqoZoUk12/QFBgLdglmvmtk33P2dhAUqIpJkUj5BELlK2ubuxx7ENhcAH7r7TgAzmwWcAChBiEi9kfK3mNx9O/C5mV0CYBFDqtnsC+BUM2tkZo2JNFDrFpOI1CsplyDMbAbwAdDfzPLM7CpgHHCVmS0AFgNjgnWPM7M84BJgipktDnbzDLASWAQsABa4+wu1fCgiIqEylfsWEZFYEnYFYWbdzewNM1sSPHB2fYx1LHgILdfMFprZsKhlV5jZiuB1RaLiFBGR2BJ2BWFmnYHO7j7XzFoCOcB33H1J1DrnAtcC5wIjgfvdfaSZtQOygSwizyTkAMPdfWtVn5mRkeG9evVKyPGIiKSinJycTe6eGWtZwnoxufu+Zwhw9x1mthToCiyJWm0M8LhHstSHZtYmSCyjgFfdfQuAmb1K5OnoGVV9Zq9evcjOzq7xYxERSVVmtqayZbXSSG1mvYChwEcVFnUF1kZN5wXzKpsvIiK1JOEJInhQ7VnghqDLaU3vf4KZZZtZdkFBQU3vXkSk3kpoggieIXgW+Ie7/yvGKvlA96jpbsG8yuYfwN2nunuWu2dlZsa8jSYiIocgkb2YDPgrsNTd76tktZnA5UFvpuOBwqDtYg5wlpm1NbO2wFnBPBERqSWJLLVxEnAZsMjM5gfzbgV6ALj7I8DLRHow5QK7ge8Hy7aY2W+AT4LtbtvXYC0iIrUjkb2Y3gWqHGQn6L10TSXLpgHTEhCaiIjEIeVKbYiISM2oD9VcRURShruzZdde1m0rJn/bbvK3FVNSVs7Vp/ap8c9SghARSSJ7S8tZX1hM3rbdrNtWzLptReRvLWJdYRH524pYt62I4pL9B8Hs0LKJEoSISF3m7mwvKt3/5L/t6xN//tYiCnbuoWIFpMyWTejSJp2BnVrxzQEd6NImna5t0r/6s02zxgmJVwlCRKSGlJaVs357caUn/3Xbiti1t2y/bdIaNQhO9k0Z1T+TLsGJv1vwZ+c2TWnSqGEox6MEISISpx3FJV+d/POCE3/0yX/99mLKK/z6b9c8ja5t0jkiszkn98uga9Sv/y5t0slokUbUcMhJRQlCRAQoK3cKduzZ71f/vpP/vuntxaX7bdOogdG5TVO6tknn+D7t9zv5d22bTpfW6aSnhfPrvyYoQYhIvbB7b2nQ8yfGyb+wiPWFxZSU7f/zv1XTRnRt24xubdMZ0bvdfr/8u7VNJ6NFExo2SM5f/zVBCUJE6jx3Z9POvfud/POC2z7rCiOJYOvukv22adjA6NSqKV3aNGVYj7b7Nfp2bZtO59ZNadk0MY2/dYUShIgkveKSMr4sjGr4rXDyX1dYzN7S/bt+Nk9rGLnN0yadId3a7Hfy79ImnY4tm9CooZ4VrooShIiEyt3Zuruk0pN//rZiNu3cs982ZpG+/13apDOoa2vOPrrTAV0/W6U3StrG37pCCUJEEmpvaTkbthcfcPL/6hbQtmKKSvbv+tm0cYOvTvQDO7c64OTfqXVT0hrp13+iKUGIyCFzd7YXlx7wtG9+1Ml/w47iAx78ymgR6fp5ZMeWjOrf4YD7/22bNdav/ySgBCEicSnaW8bTOWtZvn7HVyf//G1F7Nyzf9fPtIYNvur6eXK/jP0e+trX+Nu0cd3t+lmfKEGISJXKy51/zcvn3jnLWb+9mDbNGtO1TTo92jfjhKDvf5fgSeCubdPJaN6EBinc9bM+UYIQkUq9v3ITd7y0lMXrtjOkW2seuHQoI3q3CzssqSVKECJygNyNO7lz1lL+s3QjXdukc//YY/n24C66MqhnEpYgzGwa8C1go7sPirH8J8C4qDgGApnBcKOrgR1AGVDq7lmJilNEvrZ55x7uf20F//joC9IbN+Tm0f35wUm91WZQTyXyCuJR4E/A47EWuvs9wD0AZvZt4H8qjDt9mrtvSmB8IhIoLinj0fdX89DruewuKePSEd254YwjyWjRJOzQJESJHJP6bTPrFefqlwIzEhWLiMTm7ryw8EvumrWM/G1FnD6gA7eeO4C+HVqGHZokgdDbIMysGTAamBQ124FXzMyBKe4+tYrtJwATAHr06JHIUEVSSs6aLfzmxaXMX7uNgZ1bcffFgzmpb0bYYUkSCT1BAN8G3qtwe+lkd883sw7Aq2a2zN3fjrVxkDymAmRlZXmsdUTka2s27+Ku2ct4edF6OrRswt0XD+aiYd1SuiqpHJpkSBBjqXB7yd3zgz83mtlzwAggZoIQkfgU7i7hwddX8NgHq2nUoAE3nNGPCaccQbO0ZDgNSDIK9V+GmbUGTgXGR81rDjRw9x3B+7OA20IKUaTO21tazvQP1/DA6ysoLCrhkuHduPGs/nRs1TTs0CTJJbKb6wxgFJBhZnnAL4HGAO7+SLDaBcAr7r4ratOOwHNBHZZGwBPuPjtRcYqkKndnzuIN3DlrKas37+bkvhnceu5AjurSKuzQpI5IZC+mS+NY51Ei3WGj560ChiQmKpH6YcHabdzx0lI+Xr2Fvh1a8Lcrj2NU/0wVwJODopuPIikkf1sRd89exvPz19G+eRq3f2cQY4/rroFx5JAoQYikgB3FJUx+cyV/ffdzDPjxqD5MHNWn3g+ZKYdHCUKkDistK2fGJ2v546ufsXnXXi4Y2pWbzu5P1zbpYYcmKUAJQqQOcnfeWL6R3768jNyNOxnRux1/O28gg7u1CTs0SSFKECJ1zJJ127nj5SW8l7uZXu2bMeWy4Zx1VEc1QEuNU4IQqSM2bC/m3jnLeWZuHq3TG/PLbx/FuJE9NTazJIwShEiS2723lClvrWLq26soLS/nhyf3ZtJp/WjdTA3QklhKECJJqqzceTYnj3tfWc7GHXs475jO/HT0AHq0bxZ2aFJPKEGIJKF3V2zi9peWsGz9Dob2aMPD44cxvKeG+pTapQQhkkRWbNjBb19eyhvLC+jWNp0HLx3KtwZ3VgO0hEIJQiQJbNq5hz+8+hlPfrKWZmkN+dk5A7jixF4a6lNCpQQhEqLikjL++u7nPPzmSopKyhg/sgfXn3Ek7ZqnhR2aiBKESBjKy52ZC9Zxz5zl5G8r4oyBHfnZuQPok9ki7NBEvqIEIVLLPv58C7e/tISFeYUM6tqKey8Zwgl92ocdlsgBlCBEasnnm3Zx56ylzFm8gU6tmvL7S4ZwwdCuNNBQn5KklCBEEmzrrr3c/9oKpn+4hrRGDbjxzCP54TeOID1NDdCS3JQgRBJkT2kZj7+/hgdfX8HOPaV877ju/M+ZR9KhpYb6lLohkUOOTgO+BWx090Exlo8Cngc+D2b9y91vC5aNBu4HGgJ/cfc7ExWnSE1zd15etJ67Zi/jiy27OfXITG49dyD9O7UMOzSRg5LIK4hHgT8Bj1exzjvu/q3oGWbWEHgIOBPIAz4xs5nuviRRgYrUlLlfbOWOl5aSs2Yr/Tu25PEfjOCUIzPDDkvkkCRyTOq3zazXIWw6AsgNxqbGzJ4ExgBKEJK01m7ZzV2zl/Hiwi/JaNGE3114DN/N6k5DNUBLHRZ2G8QJZrYAWAfc5O6Lga7A2qh18oCRle3AzCYAEwB69OiRwFBFDlRYVMLkN3L523uradAArju9LxNO7UOLJmH/1xI5fGH+K54L9HT3nWZ2LvBvoN/B7sTdpwJTAbKysrxGIxSpRElZOU989AV//M9nbCsq4cKh3bjp7CPp3FpDfUrqCC1BuPv2qPcvm9lkM8sA8oHuUat2C+aJhM7d+c/Sjfxu1lJWFezihCPa83/nDWRQ19ZhhyZS40JLEGbWCdjg7m5mI4AGwGZgG9DPzHoTSQxjgf8KK06RfT7NL+T2l5bw4aotHJHZnL9cnsU3B3ZQpVVJWYns5joDGAVkmFke8EugMYC7PwJcDEw0s1KgCBjr7g6UmtkkYA6Rbq7TgrYJkVB8WVjEPXOW89y8fNo2S+O2MUdz6YgeNG6ooT4ltVnknJwasrKyPDs7O+wwJEXs3FPKlLdW8ud3VlFeDt8/uRfXnNaXVk011KekDjPLcfesWMvU1UKkgtKycp7OyeP3r3zGpp17+PaQLtx8dn+6t9NQn1K/KEGIRHnrswJ++9JSlm/YwfCebfnz5cMZ2qNt2GGJhEIJQgRYtn47v315GW9/VkCPds2YPG4Y5wzqpAZoqdeUIKRe27ijmPte+Yx/Zq+lRZNG/Py8gVx2Qk+aNFKlVRElCKmXivaW8Zd3VvHwWyvZW1rOFSf24rrT+9FWQ32KfEUJQuqV8nLnX/PyuXfOctZvL+bsoztyyzkD6Z3RPOzQRJKOEoTUG++v3MQdLy1l8brtDO7WmvvHHsvIIzTUp0hllCAk5eVu3Mmds5byn6Ub6dK6KX/83rGcP6SLhvoUqYYShKSszTv3cP9rK/jHR1+Q3rghN4/uzw9O6k3TxmqAFomHEoSknOKSMh59fzUPvZ7L7pIyLh3RnRvOOJKMFk3CDk2kTlGCkJTh7ryw8Evunr2MvK1FnD6gAz87ZwD9OmqoT5FDoQQhKSFnzRZ+8+JS5q/dxsDOrZh+1WBO7pcRdlgidZoShNRpazbv4q7Zy3h50Xo6tGzC3RcP5qJh3TTUp0gNUIKQOqlwdwkPvr6Cxz5YTaMGDbjhjH5MOOUImqXpn7RITdH/JqlT9paWM/3DNTzw+goKi0q4ZHg3bjyrPx1bNQ07NJGUU2WCMLMH4tjHdnf/eQ3FIxKTuzNn8QbunLWU1Zt3c3LfDG49dyBHdWkVdmgiKau6K4gxwC+qWecW4IAEYWbTgG8BG919UIzl44CfAgbsACa6+4Jg2epgXhlQWtlgFlI/LMzbxu0vLuXj1Vvo26EFf7vyOEb1z1SlVZEEqy5B/MHdH6tqBTOrrFj+o8CfgMcrWf45cKq7bzWzc4CpwMio5ae5+6Zq4pMUlr+tiHtmL+Pf89fRvnkat39nEGOP604jDfUpUiuqTBDu/sfqdlDZOu7+tpn1qmK796MmPwS6VfdZUn88Pz+fm59ZCMCPR/Vh4qg+tNRQnyK1Kq6fYmZ2t5m1MrPGZvaamRWY2fgajOMqYFbUtAOvmFmOmU2owc+ROmBRXiE/eWYhg7u15vWbRnHz6AFKDiIhiPda/Sx3306kTWE10Bf4SU0EYGanEUkQP42afbK7DwPOAa4xs1Oq2H6CmWWbWXZBQUFNhCQh2rJrL1dPzyGjeRqPjB9O1zbpYYckUm/FmyD23Yo6D3ja3Qtr4sPNbDDwF2CMu2/eN9/d84M/NwLPASMq24e7T3X3LHfPyszMrImwJCRl5c51M+ZRsGMPD48fTnvVThIJVbwJ4kUzWwYMB14zs0yg+HA+2Mx6AP8CLnP3z6LmNzezlvveA2cBnx7OZ0nd8PtXlvNu7iZuG3M0Q7q3CTsckXovrgfl3P0WM7sbKHT3MjPbRaQLbKXMbAYwCsgwszzgl0DjYH+PEOk+2x6YHHRX3NedtSPwXDCvEfCEu88+hGOTOmT2p+uZ/OZKLh3RnbEjeoQdjogQZ4Iws4bAyUAvM4ve5r7KtnH3S6vap7v/EPhhjPmrgCHxxCWpIXfjTm56egFDurXmV+cfHXY4IhKIt9TGC0RuKS0CyhMXjtQ3O/eUcvX0HJo0asDD44fTpJEG8xFJFvEmiG7uPjihkUi94+7c/MwCVhXsZPpVI+miHksiSSXeRupZZnZWQiORemfq26t4edF6fjp6ACf21dgNIskm3iuID4k0HDcASojUT3J3V6U0OSTv527irtnLOPeYTkw45YiwwxGRGOJNEPcBJwCL3N0TGI/UA+u2FTFpxjyOyGzB3RcPUdE9kSQV7y2mtcCnSg5yuIpLypg4PYe9peVMuWw4LZpoSBKRZBXv/85VwJtmNgvYs2+mu1fazVUkll+/sJgFeYU8Mn44fTJbhB2OiFQh3gTxefBKC14iB+3Jj79gxsdr+fGoPowe1CnscESkGvE+Sf1rADNrEUzvTGRQknoWrN3GL55fzDf6ZXDjWf3DDkdE4hBvue9BZjYPWAwsDspw65FXicvmnXuYOD2HzJZNuH/sUBo2UKO0SF0QbyP1VOB/3b2nu/cEbgT+nLiwJFWUlpVz3ZPz2LRrL1MuG0675rpDKVJXxJsgmrv7G/sm3P1NoHlCIpKUcs8ry3kvdzO3f2cQg7q2DjscETkIcfdiMrP/B/w9mB5PpGeTSKVmLfqSKW+tYtzIHnw3q3vY4YjIQYr3CuIHQCaR8RueBTKA7ycqKKn7cjfu4KanF3Bs9zb84ttHhR2OiByCeK8gznD366JnmNklwNM1H5LUdTuKS5jw9xzS0xry8PhhqtAqUkfFewXxszjnST3n7tz09ALWbN7Ng5cOo3NrVWgVqauqvIIws3OAc4GuZvZA1KJWQGkiA5O66eG3VjJn8QZ+ft5ATujTPuxwROQwVHcFsQ7IJjJYUE7UayZwdnU7N7NpZrbRzGKOKW0RD5hZrpktNLNhUcuuMLMVweuKeA9IwvPuik3cO2c53xrcmatO7h12OCJymKq8gnD3BcACM3vC3UsOYf+PAn8CHq9k+TlAv+A1EngYGGlm7YiMYZ0FOJBjZjPdfeshxCC1IG/rbq6dMZe+HVpw10WDVaFVJAXE2wbRy8yeMbMlZrZq36u6jdz9bWBLFauMAR73iA+BNmbWmcjVyavuviVICq8Co+OMVWpZpELrXErLnCmXZdFcFVpFUkK8CeJvRH7dlwKnEbkimF4Dn9+VSCnxffKCeZXNP4CZTTCzbDPLLigoqIGQ5GC4O794/lMW5Rdy3/eOpXeGnp8USRXxJoh0d38NMHdf4+6/As5LXFjxc/ep7p7l7lmZmZlhh1PvzPh4Lf/MzuPa0/ty5lEdww5HRGpQvAliTzDc6Aozm2RmFwA1Ucw/H4h+xLZbMK+y+ZJE5n2xlV/O/JRTjszkhjOODDscEalh8SaI64FmwHXAcCKlNmqiZ9FM4PKgN9PxQKG7fwnMAc4ys7Zm1hY4K5gnSWLTzj1MnD6XTq2b8sDYY1WhVSQFxTsexCfB250cRIkNM5sBjAIyzCyPSM+kxsE+HwFeJvKcRS6we9++3X2Lmf0G2Pe5t7l7VY3dUotKy8qZ9MRctu7ey7MTT6RNM1VoFUlF1T0o96ugveGQ1nH3S6vaNhjj+ppKlk0DplW1vYTjrtnL+HDVFn5/yRBVaBVJYdVdQfzQzLZXsdyAscCvaiwiSWovLlzHn9/5nMtP6MlFw7uFHY6IJFB1CeLPQMs41pF64LMNO7j5mYUM79mWn5+nCq0iqa66J6l/XVuBSHLbXlzCj/6eQ7O0RkweN4y0RvH2bxCRukr/y6Va5eXOjf9cwNotu5k8bhgdWzUNOyQRqQVKEFKth99ayatLNnDruQMZ0btd2OGISC1RgpAqvfVZAfe+spwxx3bh+yf1CjscEalFVSYIM8uoMD0+KM89wVSuM+Wt3bKb65+cR/+OLfndhceoQqtIPVPdFcQr+96Y2c+By4iMB3EmcF8C45KQFZeUcfX0HMrKnUfGD6dZmiq0itQ31f2vj/7JeCHwDXffZWZPAHMTF5aEyd35v+c+ZfG67fz1iix6qUKrSL1UXYJIN7OhRK40Grr7LgB3LzGzsoRHJ6GY/tEXPDs3j+u/2Y9vDlSFVpH6qroE8SVf30raYmad3f1LM2uPxqROSTlrtnLbC4s5rX8m13+zX9jhiEiIqntQ7rRKFm0DTqnxaCRUG3cU8+N/5NC5dTp//N5QGqhCq0i9Vm03VzPrZGadgveZZnYhMMDddyc8Oqk1JWXlTHpiHoVFJUy5bDitmzUOOyQRCVl13Vx/BHwAfGhmE4EXiYwk9y8zu6oW4pNa8ruXl/Hx51u488LBDOzcKuxwRCQJVNcGMQk4GkgH1gB93X19MIjPG8BfExyf1ILn5+cz7b3PufLEXnxnaMyhv0WkHqouQZQEt5J2m9lKd18P4O5bzcwTH54k2rL127nl2UUc16st/3fewLDDEZEkUl0bhJvZvpvR5+2baWZN49hWklxhUaRCa4umjXjov4bRuKG+UhH5WnVnhAsAB3D3vKj57YEbq9u5mY02s+Vmlmtmt8RY/gczmx+8PjOzbVHLyqKWzYzraCRu5eXO/z41n/ytRTw8bhgdVKFVRCqorpvrF5XMzwfyq9rWzBoCDxEpy5EHfGJmM919SdR+/idq/WuBoVG7KHL3Y6s7ADk0f3ojl9eWbeTX5x9NVi9VaBWRAx3yPQUzW1TNKiOAXHdf5e57gSeBMVWsfykw41Djkfi9sXwjf/jPZ1wwtCuXn9Az7HBEJElVeQURPPMQcxHQqZp9dwXWRk3nASMr+ZyeQG/g9ajZTc0sm8gT23e6+78r2XYCMAGgR48e1YQkX2zezfUz5jGgUyt+e4EqtIpI5arrxfQU8A+CdogKavKm9VjgGXePru/U093zzewI4HUzW+TuKytu6O5TgakAWVlZ6llVhaK9Zfxoeg5mxpTxw0lPaxh2SCKSxKpLEAuBe93904oLzOyMarbNB7pHTXej8naLscA10TOCdg7cfZWZvUmkfeKABCHxcXdufW4Ry9ZvZ9qVx9GjfbOwQxKRJFddG8QNwPZKll1QzbafAP3MrLeZpRFJAgf0RjKzAUBbIk9s75vX1syaBO8zgJOAJRW3lfg9/sEanpuXz/+ccSSn9e8QdjgiUgdU14vpnSqWZVezbamZTQLmAA2Bae6+2MxuA7LdfV+yGAs86e7Rt4cGAlPMrJxIErszuveTHJzs1Vv4zYtLOGNgByad1jfscESkjrD9z8uVrGR2JPAw0NHdB5nZYOB8d7890QEejKysLM/OrjJv1Tsbtxdz3oPv0jytIc9POpnW6SrCJyJfM7Mcd8+KtSzebq5/Bn4GlAC4+0Iiv/wliZWUlXPNE3PZWVzKI5cNV3IQkYMSb4Jo5u4fV5inAYOS3B0vLeWT1Vu56+LBDOikCq0icnDiTRCbzKwPQXdXM7uYyGhzkqSem5fHo++v5qqTe3P+kC5hhyMidVB13Vz3uYbIswYDzCwf+BwYl7Co5LAsWbedn/1rESN6t+OWcwaEHY6I1FHVDRjUAyLPIrj7GUAmkdHkTgb02HISKtxdwtXTc2id3lgVWkXksFR39njTzG4OCu/h7ruAZmY2HfhDwqOTg1Je7tzw1Dy+LCxi8rjhZLZsEnZIIlKHVZcghgN9gPlmdrqZXQ98TOShthGJDk4Ozv2vreCN5QX84ttHM7xn27DDEZE6rroH5bYCPwoSw3+AdcDxFcaGkCTw2tIN3P/aCi4a1o3xI3X3T0QOX3VtEG3MbArwfWA08Awwy8xOr43gJD6rN+3ihqfmc3SXVtxxwSBVaBWRGlFdL6a5wGTgGncvBV4xs2OByWa2xt0vTXSAUrXde0u5enoODRsYj4wfTtPGqtAqIjWjugRxSsXbSe4+HzjRzP47YVFJXNydW55dxPINO3js+yPo3k4VWkWk5lR5i6mqtgZ3/3PNhyMH42/vrWbmgnXcdFZ/TjkyM+xwRCTFqJN8HfXRqs389uWlnHlURyae2ifscEQkBSlB1EEbthdzzRPz6NGuGb//7hAaNFCjtIjUvHhLbUiS2FtazsTpOezeW8oT/z2SVk1VoVVEEkMJoo65/aUlzP1iGw/91zCO7Ngy7HBEJIUl9BaTmY02s+Vmlmtmt8RYfqWZFZjZ/OD1w6hlV5jZiuB1RSLjrCuezcnj8Q/W8N/f6M15gzuHHY6IpLiEXUEE9ZseAs4E8oBPzGxmjKFDn3L3SRW2bQf8EsgiUmI8J9h2a6LiTXaf5hdy63OLOP6Idvx0tCq0ikjiJfIKYgSQG1SC3Qs8CYyJc9uzgVfdfUuQFF4l8iR3vbRt916unp5Du+Zp/Om/htFIFVpFpBYk8kzTFVgbNZ0XzKvoIjNbaGbPmFn3g9w25ZWVO9c9OZ+N2/cwedwwMlqoQquI1I6wf4q+APRy98FErhIeO9gdmNkEM8s2s+yCgoIaDzBsf/zPZ7z9WQG/Ov9ohvZQhVYRqT2JTBD5QPeo6W7BvK+4+2Z33xNM/oVIefG4to3ax1R3z3L3rMzM1Hqa+NUlG3jw9Vy+m9WNS0d0r34DEZEalMgE8QnQz8x6m1kaMBaYGb2CmUV3xTkfWBq8nwOcZWZtzawtcFYwr95YVbCT/31qPoO7tea2MarQKiK1L2G9mNy91MwmETmxNwSmuftiM7sNyHb3mcB1ZnY+UApsAa4Mtt1iZr8hkmQAbnP3LYmKNdns2hOp0NqooTF53DBVaBWRUJi7hx1DjcnKyvLs7Oywwzgs7s61M+bx8qIvefwHIzm5X0bYIYlICjOzHHfPirUs7EZqqeCv737Oiwu/5CdnD1ByEJFQKUEkkQ9WbuZ3s5Yx+uhOXH3qEWGHIyL1nBJEkviysIhrZ8ylV/tm3HPJYDVKi0joVKwvCewpLWPi9LkU7S3jyQnH01IVWkUkCShBJIHbXljC/LXbeHjcMPp2UIVWEUkOusUUsn9mr+UfH33Bj049gnOOUYVWEUkeShAhWpRXyM///Skn9W3PT87qH3Y4IiL7UYIIyZZdkQqtGc3TeGDsUFVoFZGkozaIEJSVO9fNmEfBjj08ffUJtFeFVhFJQkoQIfj9K8t5N3cTd110DEO6twk7HBGRmHRfo5bN/nQ9k99cyaUjuvO943qEHY6ISKWUIGpR7sad3PT0AoZ0a82vzj867HBERKqkBFFLdgYVWps0asDD44fTpJEqtIpIclMbRC1wd25+ZgGrCnYy/aqRdGmTHnZIIiLV0hVELZj69ipeXrSeW84ZwIl9VaFVROoGJYgEez93E3fNXsZ5x3Tmv7+hCq0iUncoQSTQum1FTJoxjyMyW3DXxarQKiJ1S0IThJmNNrPlZpZrZrfEWP6/ZrbEzBaa2Wtm1jNqWZmZzQ9eMytum+yKS8qYOD2HvaXlTLlsOC2aqLlHROqWhJ21zKwh8BBwJpAHfGJmM919SdRq84Asd99tZhOBu4HvBcuK3P3YRMWXaL9+YTEL8gqZctlw+mS2CDscEZGDlsgriBFArruvcve9wJPAmOgV3P0Nd98dTH4IdEtgPLXmyY+/YMbHa/nxqD6cfXSnsMMRETkkiUwQXYG1UdN5wbzKXAXMippuambZZvahmX2nso3MbEKwXnZBQcFhBVwTFqzdxi+eX8w3+mVwoyq0ikgdlhQ3xs1sPJAFnBo1u6e755vZEcDrZrbI3VdW3NbdpwJTAbKysrxWAq7E5p17mDg9h8yWTXhg7FAaNlCjtIjUXYm8gsgHukdNdwvm7cfMzgD+Dzjf3ffsm+/u+cGfq4A3gaEJjPWwlZaVc92T89i0ay9TLhtO2+ZpYYckInJYEpkgPgH6mVlvM0sDxgL79UYys6HAFCLJYWPU/LZm1iR4nwGcBEQ3biede15Zznu5m7njO4MY1LV12OGIiBy2hN1icvdSM5sEzAEaAtPcfbGZ3QZku/tM4B6gBfB08IzAF+5+PjAQmGJm5USS2J0Vej8llVmLvmTKW6sYN7IHl2R1r34DEZE6wNxDvW1fo7Kysjw7O7tWPzN34w7G/Ok9+nVsyVM/Ol5F+ESkTjGzHHfPirVMT1Ifhh3FJUz4ew7paQ15ePwwJQcRSSlJ0YupLnJ3bnp6AWs27+YfPxxJ59aq0CoiqUVXEIfo4bdWMmfxBn52zgCOP6J92OGIiNQ4JYhD8O6KTdw7ZznfGtyZq07uHXY4IiIJoQRxkPK27ubaGXPp26EFd12kCq0ikrqUIA5CpELrXErLnCmXZdFcFVpFJIXpDBcnd+cXz3/KovxC/nJ5Fr0zmocdkohIQukKIk4zPl7LP7PzuPb0vpxxVMewwxERSTgliDjM+2Irv5z5KacemckNZxwZdjgiIrVCCaIam3buYeL0uXRq3ZT7xx6rCq0iUm+oDaIKpWXlTHpiLlt37+XZiSfSppkqtIpI/aEEUYW7Zi/jw1VbuO+7Q1ShVUTqHd1iqsSLC9fx53c+5/ITenLhsJQYCVVE5KAoQcTw2YYd3PzMQob3bMvPzzsq7HBEREKhBFHB9uISfvT3HJqlNWLyuGGkNdJfkYjUTzr7RSkvd2785wLWbtnN5HHD6NiqadghiYiERgkiysNvreTVJRu49dyBjOjdLuxwRERCldAEYWajzWy5meWa2S0xljcxs6eC5R+ZWa+oZT8L5i83s7MTGSfAW58VcO8ryxlzbBe+f1KvatcXEUl1CUsQZtYQeAg4BzgKuNTMKrb4XgVsdfe+wB+Au4JtjwLGAkcDo4HJwf4SYu2W3Vz/5Dz6d2zJ7y48RhVaRURI7BXECCDX3Ve5+17gSWBMhXXGAI8F758BvmmRs/MY4El33+PunwO5wf5qXHFJGVdPz6Gs3Hlk/HCapenREBERSGyC6AqsjZrOC+bFXMfdS4FCoH2c2wJgZhPMLNvMsgsKCg46SHfo37Elf/zesfRShVYRka/U+Z/L7j4VmAqQlZXlB7t9elpD7vvesTUdlohInZfIK4h8oHvUdLdgXsx1zKwR0BrYHOe2IiKSQIlMEJ8A/cyst5mlEWl0nllhnZnAFcH7i4HX3d2D+WODXk69gX7AxwmMVUREKkjYLSZ3LzWzScAcoCEwzd0Xm9ltQLa7zwT+CvzdzHKBLUSSCMF6/wSWAKXANe5elqhYRUTkQBb5wZ4asrKyPDs7O+wwRETqDDPLcfesWMv0JLWIiMSkBCEiIjEpQYiISExKECIiElNKNVKbWQGw5hA3zwA21WA4YUqVY0mV4wAdSzJKleOAwzuWnu6eGWtBSiWIw2Fm2ZW15Nc1qXIsqXIcoGNJRqlyHJC4Y9EtJhERiUkJQkREYlKC+NrUsAOoQalyLKlyHKBjSUapchyQoGNRG4SIiMSkKwgREYlJCUJERGKqdwnCzEab2XIzyzWzW2Isb2JmTwXLPzKzXiGEWa04juNKMysws/nB64dhxFkdM5tmZhvN7NNKlpuZPRAc50IzG1bbMcYrjmMZZWaFUd/JL2o7xniZWXcze8PMlpjZYjO7PsY6Sf/dxHkcdeJ7MbOmZvaxmS0IjuXXMdap2fOXu9ebF5Gy4yuBI4A0YAFwVIV1fgw8ErwfCzwVdtyHeBxXAn8KO9Y4juUUYBjwaSXLzwVmAQYcD3wUdsyHcSyjgBfDjjPOY+kMDAvetwQ+i/FvLOm/mziPo058L8Hfc4vgfWPgI+D4CuvU6Pmrvl1BjABy3X2Vu+8FngTGVFhnDPBY8P4Z4JtmZrUYYzziOY46wd3fJjIWSGXGAI97xIdAGzPrXDvRHZw4jqXOcPcv3X1u8H4HsJQDx4VP+u8mzuOoE4K/553BZOPgVbGXUY2ev+pbgugKrI2azuPAfyxfrePupUAh0L5WootfPMcBcFFw6f+MmXWPsbwuiPdY64oTglsEs8zs6LCDiUdwm2IokV+s0erUd1PFcUAd+V7MrKGZzQc2Aq+6e6XfSU2cv+pbgqhPXgB6uftg4FW+/lUh4ZlLpO7NEOBB4N/hhlM9M2sBPAvc4O7bw47nUFVzHHXme3H3Mnc/FugGjDCzQYn8vPqWIPKB6F/S3YJ5Mdcxs0ZAa2BzrUQXv2qPw903u/ueYPIvwPBaiq2mxfOd1Qnuvn3fLQJ3fxlobGYZIYdVKTNrTOSk+g93/1eMVerEd1PdcdS17wXA3bcBbwCjKyyq0fNXfUsQnwD9zKy3maURacSZWWGdmcAVwfuLgdc9aPFJItUeR4V7wecTufdaF80ELg96zBwPFLr7l2EHdSjMrNO++8FmNoLI/79k+/EBRHooERkzfqm731fJakn/3cRzHHXlezGzTDNrE7xPB84EllVYrUbPX40OdcO6yN1LzWwSMIdIT6Bp7r7YzG4Dst19JpF/TH83s1wiDY5jw4s4tjiP4zozOx8oJXIcV4YWcBXMbAaRXiQZZpYH/JJI4xvu/gjwMpHeMrnAbuD74URavTiO5WJgopmVAkXA2CT88bHPScBlwKLgnjfArUAPqFPfTTzHUVe+l87AY2bWkEgS+6e7v5jI85dKbYiISEz17RaTiIjESQlCRERiUoIQEZGYlCBERCQmJQgREYlJCUJERGJSghCpwMzaR5V+Xm9m+cH7nWY2OQGf96iZfW5mV1eyfGfwZ599cdR0DCKx6DkIkSqY2a+Ane5+bwI/41Ei5aafqWT5TndvUdm0SKLoCkIkTsHAMi8G739lZo+Z2TtmtsbMLjSzu81skZnNDur/YGbDzewtM8sxsznxlMMOSqh8EOzr9kQfl0hllCBEDl0f4HQita6mA2+4+zFEyjWcFySJB4GL3X04MA24I4793g88HOwrqWobSf1Sr2oxidSwWe5eYmaLiNTEmh3MXwT0AvoDg4BXg1pwDYnvhH8ScFHw/u/AXTUYs0jclCBEDt0eAHcvN7OSqAJv5UT+bxmw2N1POIR9q3FQQqdbTCKJsxzINLMTIDIuQZyjlb3H11U4xyUqOJHqKEGIJEgwXvjFwF1mtgCYD5wYx6bXA9cEt66SdghPSX3q5ioSsuq6ucZYX91cpVboCkIkfIXAbyp7UG6ffQ/KARtqJSqp93QFISIiMekKQkREYlKCEBGRmJQgREQkJiUIERGJ6f8DrhVdI0mfnDYAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"results = openmc.deplete.ResultsList.from_hdf5('depletion_results.h5')\n",
"time, n_Xe135 = results.get_atoms('1', 'Xe135')\n",
"days = 24*60*60\n",
"plt.plot(time/days, n_Xe135)\n",
"plt.xlabel('Time [d]')\n",
"plt.ylabel('Xe135 [atoms]')"
]
},
{
"cell_type": "markdown",
"id": "b76b5435",
"metadata": {},
"source": [
"For the next depletion run, we need to get a new set of materials that is equivalent to the last depletion step. For this, we can use the `ResultsList.export_to_materials` method. The change we'll make to the material is to remove ${135}$Xe so that it starts with a near-zero concentration."
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "0b835890",
"metadata": {},
"outputs": [],
"source": [
"# Get materials at the end of the last simulation\n",
"model.materials = results.export_to_materials(len(time_steps))\n",
"\n",
"# Now change the material by getting rid of Xe135\n",
"model.materials[0].remove_nuclide('Xe135')"
]
},
{
"cell_type": "markdown",
"id": "e3a11a2c",
"metadata": {},
"source": [
"Now run the depletion integrator again --- note that we don't supply previous results (which would override the material change we made above)."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "1e15dfb0",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" %%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" %%%%%%%%%%%%%%%%%%%%%%%%\n",
" ############### %%%%%%%%%%%%%%%%%%%%%%%%\n",
" ################## %%%%%%%%%%%%%%%%%%%%%%%\n",
" ################### %%%%%%%%%%%%%%%%%%%%%%%\n",
" #################### %%%%%%%%%%%%%%%%%%%%%%\n",
" ##################### %%%%%%%%%%%%%%%%%%%%%\n",
" ###################### %%%%%%%%%%%%%%%%%%%%\n",
" ####################### %%%%%%%%%%%%%%%%%%\n",
" ####################### %%%%%%%%%%%%%%%%%\n",
" ###################### %%%%%%%%%%%%%%%%%\n",
" #################### %%%%%%%%%%%%%%%%%\n",
" ################# %%%%%%%%%%%%%%%%%\n",
" ############### %%%%%%%%%%%%%%%%\n",
" ############ %%%%%%%%%%%%%%%\n",
" ######## %%%%%%%%%%%%%%\n",
" %%%%%%%%%%%\n",
"\n",
" | The OpenMC Monte Carlo Code\n",
" Copyright | 2011-2021 MIT, UChicago Argonne LLC, and contributors\n",
" License | https://docs.openmc.org/en/latest/license.html\n",
" Version | 0.13.0-dev\n",
" Git SHA1 | 65b02c57b9c09cf0c4c914a31034ed076688165e\n",
" Date/Time | 2022-01-09 13:43:53\n",
" OpenMP Threads | 12\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading geometry XML file...\n",
" Reading I135 from /opt/data/hdf5/nndc_hdf5_v15/I135.h5\n",
" Reading Xe136 from /opt/data/hdf5/nndc_hdf5_v15/Xe136.h5\n",
" Reading Cs135 from /opt/data/hdf5/nndc_hdf5_v15/Cs135.h5\n",
" Reading Gd157 from /opt/data/hdf5/nndc_hdf5_v15/Gd157.h5\n",
" Reading Gd156 from /opt/data/hdf5/nndc_hdf5_v15/Gd156.h5\n",
" Reading U234 from /opt/data/hdf5/nndc_hdf5_v15/U234.h5\n",
" Reading U235 from /opt/data/hdf5/nndc_hdf5_v15/U235.h5\n",
" Reading U238 from /opt/data/hdf5/nndc_hdf5_v15/U238.h5\n",
" Minimum neutron data temperature: 0 K\n",
" Maximum neutron data temperature: 1.7976931348623157e+308 K\n",
" Preparing distributed cell instances...\n",
" Reading plot XML file...\n",
" Writing summary.h5 file...\n",
" Reading Xe135 from /opt/data/hdf5/nndc_hdf5_v15/Xe135.h5\n",
" Maximum neutron transport energy: 20000000 eV for I135\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.44499\n",
" 2/1 0.37975\n",
" 3/1 0.35298\n",
" 4/1 0.34118\n",
" 5/1 0.33717\n",
" 6/1 0.34054\n",
" 7/1 0.32571\n",
" 8/1 0.33094\n",
" 9/1 0.31767\n",
" 10/1 0.33711\n",
" 11/1 0.32746\n",
" 12/1 0.33361 0.33053 +/- 0.00307\n",
" 13/1 0.33729 0.33278 +/- 0.00287\n",
" 14/1 0.31204 0.32760 +/- 0.00557\n",
" 15/1 0.35363 0.33281 +/- 0.00676\n",
" 16/1 0.32952 0.33226 +/- 0.00555\n",
" 17/1 0.33687 0.33292 +/- 0.00473\n",
" 18/1 0.32477 0.33190 +/- 0.00423\n",
" 19/1 0.32252 0.33086 +/- 0.00387\n",
" 20/1 0.33911 0.33168 +/- 0.00356\n",
" 21/1 0.33617 0.33209 +/- 0.00324\n",
" 22/1 0.32190 0.33124 +/- 0.00308\n",
" 23/1 0.33006 0.33115 +/- 0.00284\n",
" 24/1 0.32309 0.33057 +/- 0.00269\n",
" 25/1 0.32976 0.33052 +/- 0.00250\n",
" 26/1 0.33036 0.33051 +/- 0.00234\n",
" 27/1 0.32291 0.33006 +/- 0.00224\n",
" 28/1 0.33655 0.33042 +/- 0.00215\n",
" 29/1 0.31822 0.32978 +/- 0.00213\n",
" 30/1 0.34582 0.33058 +/- 0.00217\n",
" 31/1 0.32458 0.33030 +/- 0.00209\n",
" 32/1 0.32461 0.33004 +/- 0.00201\n",
" 33/1 0.32295 0.32973 +/- 0.00194\n",
" 34/1 0.34932 0.33055 +/- 0.00203\n",
" 35/1 0.34864 0.33127 +/- 0.00208\n",
" 36/1 0.34792 0.33191 +/- 0.00210\n",
" 37/1 0.33887 0.33217 +/- 0.00203\n",
" 38/1 0.31751 0.33165 +/- 0.00203\n",
" 39/1 0.34412 0.33208 +/- 0.00200\n",
" 40/1 0.32791 0.33194 +/- 0.00194\n",
" 41/1 0.33789 0.33213 +/- 0.00189\n",
" 42/1 0.33964 0.33236 +/- 0.00184\n",
" 43/1 0.32730 0.33221 +/- 0.00179\n",
" 44/1 0.33543 0.33230 +/- 0.00174\n",
" 45/1 0.33573 0.33240 +/- 0.00169\n",
" 46/1 0.32221 0.33212 +/- 0.00167\n",
" 47/1 0.34005 0.33233 +/- 0.00164\n",
" 48/1 0.33499 0.33240 +/- 0.00160\n",
" 49/1 0.34347 0.33269 +/- 0.00158\n",
" 50/1 0.33390 0.33272 +/- 0.00154\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.8436e-01 seconds\n",
" Reading cross sections = 4.7836e-01 seconds\n",
" Total time in simulation = 4.4880e-02 seconds\n",
" Time in transport only = 3.4566e-02 seconds\n",
" Time in inactive batches = 6.1937e-03 seconds\n",
" Time in active batches = 3.8687e-02 seconds\n",
" Time synchronizing fission bank = 2.5904e-03 seconds\n",
" Sampling source sites = 2.3432e-03 seconds\n",
" SEND/RECV source sites = 2.3338e-04 seconds\n",
" Time accumulating tallies = 3.0039e-03 seconds\n",
" Time writing statepoints = 2.1722e-03 seconds\n",
" Total time for finalization = 5.6857e-05 seconds\n",
" Total time elapsed = 5.3259e-01 seconds\n",
" Calculation Rate (inactive) = 1.61456e+06 particles/second\n",
" Calculation Rate (active) = 1.03395e+06 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33326 +/- 0.00231\n",
" k-effective (Track-length) = 0.33272 +/- 0.00154\n",
" k-effective (Absorption) = 0.33551 +/- 0.00539\n",
" Combined k-effective = 0.33285 +/- 0.00157\n",
" Leakage Fraction = 0.86215 +/- 0.00216\n",
"\n",
" Creating state point openmc_simulation_n0.h5...\n",
" Maximum neutron transport energy: 20000000 eV for I135\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.44304\n",
" 2/1 0.37678\n",
" 3/1 0.35418\n",
" 4/1 0.34945\n",
" 5/1 0.33601\n",
" 6/1 0.32262\n",
" 7/1 0.32733\n",
" 8/1 0.32685\n",
" 9/1 0.33419\n",
" 10/1 0.33031\n",
" 11/1 0.32766\n",
" 12/1 0.33874 0.33320 +/- 0.00554\n",
" 13/1 0.33716 0.33452 +/- 0.00346\n",
" 14/1 0.33721 0.33519 +/- 0.00254\n",
" 15/1 0.33802 0.33576 +/- 0.00204\n",
" 16/1 0.32795 0.33446 +/- 0.00212\n",
" 17/1 0.33964 0.33520 +/- 0.00194\n",
" 18/1 0.32893 0.33441 +/- 0.00185\n",
" 19/1 0.31602 0.33237 +/- 0.00262\n",
" 20/1 0.33439 0.33257 +/- 0.00235\n",
" 21/1 0.34165 0.33340 +/- 0.00228\n",
" 22/1 0.34085 0.33402 +/- 0.00217\n",
" 23/1 0.32933 0.33366 +/- 0.00203\n",
" 24/1 0.33219 0.33355 +/- 0.00188\n",
" 25/1 0.33482 0.33364 +/- 0.00175\n",
" 26/1 0.32928 0.33336 +/- 0.00166\n",
" 27/1 0.34823 0.33424 +/- 0.00179\n",
" 28/1 0.35674 0.33549 +/- 0.00210\n",
" 29/1 0.33788 0.33562 +/- 0.00199\n",
" 30/1 0.32127 0.33490 +/- 0.00202\n",
" 31/1 0.33130 0.33473 +/- 0.00193\n",
" 32/1 0.34438 0.33517 +/- 0.00189\n",
" 33/1 0.33792 0.33529 +/- 0.00181\n",
" 34/1 0.32681 0.33493 +/- 0.00177\n",
" 35/1 0.33487 0.33493 +/- 0.00170\n",
" 36/1 0.34786 0.33543 +/- 0.00170\n",
" 37/1 0.33966 0.33558 +/- 0.00165\n",
" 38/1 0.34130 0.33579 +/- 0.00160\n",
" 39/1 0.32394 0.33538 +/- 0.00160\n",
" 40/1 0.32715 0.33511 +/- 0.00157\n",
" 41/1 0.34174 0.33532 +/- 0.00153\n",
" 42/1 0.33771 0.33539 +/- 0.00148\n",
" 43/1 0.33322 0.33533 +/- 0.00144\n",
" 44/1 0.33487 0.33531 +/- 0.00140\n",
" 45/1 0.32783 0.33510 +/- 0.00137\n",
" 46/1 0.32464 0.33481 +/- 0.00137\n",
" 47/1 0.33517 0.33482 +/- 0.00133\n",
" 48/1 0.33256 0.33476 +/- 0.00129\n",
" 49/1 0.32880 0.33461 +/- 0.00127\n",
" 50/1 0.33417 0.33460 +/- 0.00124\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 5.4805e-02 seconds\n",
" Time in transport only = 4.0073e-02 seconds\n",
" Time in inactive batches = 5.6291e-03 seconds\n",
" Time in active batches = 4.9176e-02 seconds\n",
" Time synchronizing fission bank = 2.7622e-03 seconds\n",
" Sampling source sites = 2.4806e-03 seconds\n",
" SEND/RECV source sites = 2.5121e-04 seconds\n",
" Time accumulating tallies = 5.3500e-03 seconds\n",
" Time writing statepoints = 5.2063e-03 seconds\n",
" Total time for finalization = 6.2409e-05 seconds\n",
" Total time elapsed = 5.8172e-02 seconds\n",
" Calculation Rate (inactive) = 1.77647e+06 particles/second\n",
" Calculation Rate (active) = 813405 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33615 +/- 0.00292\n",
" k-effective (Track-length) = 0.33460 +/- 0.00124\n",
" k-effective (Absorption) = 0.33953 +/- 0.00351\n",
" Combined k-effective = 0.33477 +/- 0.00129\n",
" Leakage Fraction = 0.86047 +/- 0.00140\n",
"\n",
" Creating state point openmc_simulation_n1.h5...\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" Maximum neutron transport energy: 20000000 eV for I135\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.43749\n",
" 2/1 0.37037\n",
" 3/1 0.35902\n",
" 4/1 0.33888\n",
" 5/1 0.32159\n",
" 6/1 0.33254\n",
" 7/1 0.33918\n",
" 8/1 0.32945\n",
" 9/1 0.32723\n",
" 10/1 0.34119\n",
" 11/1 0.33134\n",
" 12/1 0.31736 0.32435 +/- 0.00699\n",
" 13/1 0.31267 0.32046 +/- 0.00561\n",
" 14/1 0.32451 0.32147 +/- 0.00409\n",
" 15/1 0.34099 0.32537 +/- 0.00503\n",
" 16/1 0.34485 0.32862 +/- 0.00523\n",
" 17/1 0.34059 0.33033 +/- 0.00474\n",
" 18/1 0.33071 0.33038 +/- 0.00411\n",
" 19/1 0.35215 0.33280 +/- 0.00436\n",
" 20/1 0.34348 0.33387 +/- 0.00404\n",
" 21/1 0.32824 0.33336 +/- 0.00369\n",
" 22/1 0.32528 0.33268 +/- 0.00344\n",
" 23/1 0.33203 0.33263 +/- 0.00316\n",
" 24/1 0.32846 0.33233 +/- 0.00294\n",
" 25/1 0.33623 0.33259 +/- 0.00275\n",
" 26/1 0.33216 0.33257 +/- 0.00257\n",
" 27/1 0.33786 0.33288 +/- 0.00244\n",
" 28/1 0.35845 0.33430 +/- 0.00270\n",
" 29/1 0.34805 0.33502 +/- 0.00266\n",
" 30/1 0.33342 0.33494 +/- 0.00252\n",
" 31/1 0.32471 0.33445 +/- 0.00245\n",
" 32/1 0.34187 0.33479 +/- 0.00236\n",
" 33/1 0.33047 0.33460 +/- 0.00226\n",
" 34/1 0.33385 0.33457 +/- 0.00216\n",
" 35/1 0.34163 0.33485 +/- 0.00209\n",
" 36/1 0.33033 0.33468 +/- 0.00202\n",
" 37/1 0.33540 0.33471 +/- 0.00194\n",
" 38/1 0.33524 0.33473 +/- 0.00187\n",
" 39/1 0.33022 0.33457 +/- 0.00181\n",
" 40/1 0.33805 0.33469 +/- 0.00176\n",
" 41/1 0.33910 0.33483 +/- 0.00170\n",
" 42/1 0.33385 0.33480 +/- 0.00165\n",
" 43/1 0.32752 0.33458 +/- 0.00162\n",
" 44/1 0.34962 0.33502 +/- 0.00163\n",
" 45/1 0.33518 0.33502 +/- 0.00158\n",
" 46/1 0.33601 0.33505 +/- 0.00154\n",
" 47/1 0.34462 0.33531 +/- 0.00152\n",
" 48/1 0.33196 0.33522 +/- 0.00148\n",
" 49/1 0.32787 0.33503 +/- 0.00145\n",
" 50/1 0.33053 0.33492 +/- 0.00142\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 6.2339e-02 seconds\n",
" Time in transport only = 4.2191e-02 seconds\n",
" Time in inactive batches = 6.5212e-03 seconds\n",
" Time in active batches = 5.5818e-02 seconds\n",
" Time synchronizing fission bank = 2.9475e-03 seconds\n",
" Sampling source sites = 2.6677e-03 seconds\n",
" SEND/RECV source sites = 2.6321e-04 seconds\n",
" Time accumulating tallies = 1.0428e-02 seconds\n",
" Time writing statepoints = 6.1328e-03 seconds\n",
" Total time for finalization = 7.4346e-05 seconds\n",
" Total time elapsed = 6.5828e-02 seconds\n",
" Calculation Rate (inactive) = 1.53346e+06 particles/second\n",
" Calculation Rate (active) = 716620 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33584 +/- 0.00213\n",
" k-effective (Track-length) = 0.33492 +/- 0.00142\n",
" k-effective (Absorption) = 0.33614 +/- 0.00455\n",
" Combined k-effective = 0.33514 +/- 0.00124\n",
" Leakage Fraction = 0.86230 +/- 0.00184\n",
"\n",
" Creating state point openmc_simulation_n2.h5...\n",
" Maximum neutron transport energy: 20000000 eV for I135\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.43652\n",
" 2/1 0.36965\n",
" 3/1 0.34291\n",
" 4/1 0.33751\n",
" 5/1 0.33269\n",
" 6/1 0.33107\n",
" 7/1 0.32994\n",
" 8/1 0.33638\n",
" 9/1 0.34974\n",
" 10/1 0.34152\n",
" 11/1 0.31659\n",
" 12/1 0.34465 0.33062 +/- 0.01403\n",
" 13/1 0.32762 0.32962 +/- 0.00816\n",
" 14/1 0.32287 0.32794 +/- 0.00601\n",
" 15/1 0.33523 0.32939 +/- 0.00488\n",
" 16/1 0.34528 0.33204 +/- 0.00478\n",
" 17/1 0.32589 0.33116 +/- 0.00414\n",
" 18/1 0.34534 0.33294 +/- 0.00400\n",
" 19/1 0.34441 0.33421 +/- 0.00375\n",
" 20/1 0.32839 0.33363 +/- 0.00340\n",
" 21/1 0.32874 0.33318 +/- 0.00311\n",
" 22/1 0.33208 0.33309 +/- 0.00284\n",
" 23/1 0.32915 0.33279 +/- 0.00263\n",
" 24/1 0.33391 0.33287 +/- 0.00244\n",
" 25/1 0.33268 0.33286 +/- 0.00227\n",
" 26/1 0.33532 0.33301 +/- 0.00213\n",
" 27/1 0.33315 0.33302 +/- 0.00200\n",
" 28/1 0.33462 0.33311 +/- 0.00189\n",
" 29/1 0.33976 0.33346 +/- 0.00182\n",
" 30/1 0.32486 0.33303 +/- 0.00178\n",
" 31/1 0.33359 0.33305 +/- 0.00169\n",
" 32/1 0.34359 0.33353 +/- 0.00168\n",
" 33/1 0.33078 0.33341 +/- 0.00161\n",
" 34/1 0.32398 0.33302 +/- 0.00159\n",
" 35/1 0.32786 0.33281 +/- 0.00154\n",
" 36/1 0.33146 0.33276 +/- 0.00148\n",
" 37/1 0.33186 0.33273 +/- 0.00143\n",
" 38/1 0.33356 0.33276 +/- 0.00137\n",
" 39/1 0.33102 0.33270 +/- 0.00133\n",
" 40/1 0.33215 0.33268 +/- 0.00128\n",
" 41/1 0.33335 0.33270 +/- 0.00124\n",
" 42/1 0.34316 0.33303 +/- 0.00125\n",
" 43/1 0.34486 0.33339 +/- 0.00126\n",
" 44/1 0.33854 0.33354 +/- 0.00123\n",
" 45/1 0.33576 0.33360 +/- 0.00120\n",
" 46/1 0.32339 0.33332 +/- 0.00120\n",
" 47/1 0.33626 0.33340 +/- 0.00117\n",
" 48/1 0.34710 0.33376 +/- 0.00119\n",
" 49/1 0.33917 0.33390 +/- 0.00117\n",
" 50/1 0.33579 0.33394 +/- 0.00114\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 0.0000e+00 seconds\n",
" Reading cross sections = 0.0000e+00 seconds\n",
" Total time in simulation = 6.8805e-02 seconds\n",
" Time in transport only = 4.7870e-02 seconds\n",
" Time in inactive batches = 6.3043e-03 seconds\n",
" Time in active batches = 6.2501e-02 seconds\n",
" Time synchronizing fission bank = 2.7223e-03 seconds\n",
" Sampling source sites = 2.4579e-03 seconds\n",
" SEND/RECV source sites = 2.4881e-04 seconds\n",
" Time accumulating tallies = 9.6293e-03 seconds\n",
" Time writing statepoints = 6.1461e-03 seconds\n",
" Total time for finalization = 6.4295e-05 seconds\n",
" Total time elapsed = 7.2150e-02 seconds\n",
" Calculation Rate (inactive) = 1.58623e+06 particles/second\n",
" Calculation Rate (active) = 639989 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.33470 +/- 0.00236\n",
" k-effective (Track-length) = 0.33394 +/- 0.00114\n",
" k-effective (Absorption) = 0.33728 +/- 0.00486\n",
" Combined k-effective = 0.33432 +/- 0.00099\n",
" Leakage Fraction = 0.86203 +/- 0.00198\n",
"\n",
" Creating state point openmc_simulation_n3.h5...\n"
]
}
],
"source": [
"# We need to create a new operator; otherwise, it will still use the number densities\n",
"# from the end of the previous simulation\n",
"op = openmc.deplete.Operator(model, chain_file)\n",
"\n",
"integrator = openmc.deplete.PredictorIntegrator(op, time_steps, timestep_units='d', power=power)\n",
"integrator.integrate()"
]
},
{
"cell_type": "markdown",
"id": "66f5229e",
"metadata": {},
"source": [
"And finally let's look at the $^{135}$Xe density again:"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "4b205c9c",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Xe135 [atoms]')"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAERCAYAAABhKjCtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqpklEQVR4nO3deZwU9Z3/8deHubjPGRC5YRBBvEdEIB5REU9MYn6L0UQTXZOs5txNYrIqBtlNYnZjYoxRElk1iUdioqJBERONJ8pIVBgOGe5BhIHhPub8/P7oAtuxZ7qBrqnumffz8egH3d+q6n4XDfXpqm9Vfc3dERERaaxd1AFERCQzqUCIiEhCKhAiIpKQCoSIiCSkAiEiIgmpQIiISEKtrkCY2Uwz22Rmi1KY93QzW2BmdWZ2WaNpt5tZmZktMbM7zczCSy0iknlaXYEA7gcmpTjvWuBq4KH4RjMbB4wHjgNGA6cAZ6QtoYhIFmh1BcLdXwKq4tvMbJiZPWtmb5nZy2Z2dDDvand/F2ho/DZAeyAfKADygI3hpxcRyRytrkA0YQbwNXc/GfgP4O7mZnb314EXgA3BY467Lwk9pYhIBsmNOkDYzKwzMA74U1w3QkGSZYqBkUD/oGmumX3C3V8OLaiISIZp9QWC2F7SNnc/4SCW+RQwz913AZjZM8BpgAqEiLQZrf4Qk7vvAFaZ2WcBLOb4JIutBc4ws1wzyyPWQa1DTCLSprS6AmFmDwOvAyPMrMLMrgGuAK4xs3eAMmByMO8pZlYBfBa418zKgrd5DFgBLATeAd5x96daeFVERCJlut23iIgk0ur2IEREJD1aVSd1YWGhDx48OOoYIiJZ46233trs7kWJpoVWIMxsAPAg0IfYhWcz3P0XjeYx4BfABcAe4Gp3XxBMuwq4KZh1urs/kOwzBw8eTGlpafpWQkSklTOzNU1NC3MPog74d3dfYGZdgLfMbK67L46b53xgePA4Ffg1cKqZ9QSmAiXEistbZjbL3beGmFdEROKE1gfh7hv27w24+05ip4n2azTbZOBBj5kHdDezvsB5wFx3rwqKwlxSv7+SiIikQYt0UpvZYOBE4I1Gk/oB6+JeVwRtTbUneu/rzKzUzEorKyvTlllEpK0LvUAEt7r4M/DN4KK1tHL3Ge5e4u4lRUUJ+1lEROQQhFoggquQ/wz8wd3/kmCW9cCAuNf9g7am2kVEpIWEViCCM5TuA5a4+8+amG0W8IXg9hdjge3uvgGYA0w0sx5m1gOYGLSJiEgLCfMspvHA54GFZvZ20PYDYCCAu98DzCZ2ims5sdNcvxhMqzKz24D5wXLT3P0jYzyIiEi4QisQ7v4K0OwwnR67z8f1TUybCcwMIZqISEarb3D21NSxt6aePcFjb23dh88PtMfmyc1px1fPHJb2HK3qSmoRkZZSW98Qt7GuCzbi+zfgCTbmtR9u8BMvE2vbXVNPTV3jQS6bV9SlQAVCRCRV7k5NfcNHf4Xv3zDXxp7vrq47sIGO37B/uFFP0FYTW6a2/uBudJqf246O+Tl0zMuhQ34OHfNz6ZCfQ89O+fTvkUOHvNzY9Pz903PokJ9Lx7z4ttyPTO+YF3uP/NxwupNVIEQkMu7OvtqGj/2ajj+80nij/rGNeVMb+Np66hsObiPePq9dbMMdbJT3b4z7dG0f2yjnxW24G2/MP7aBz6VT8LxDXg65Odl3b1QVCBFJ2aad+6jaXZPwOPiHG+vGv7hjx893V+/fcH90/oMZccCMAxvvDsEv6I4Fsdc9OuZ/bKMdv7Fv7hd4x2Aj3q5ds92mbY4KhIg0afveWuat3MKr5Zt5pXwzKyt3J10mp53FHUb58Nd2l/a59O5SkPgXeN6Hh1w+tjHP+3Bj3z6vHXFjy0vIVCBE5IDqunoWrNl2oCC8W7GNBoeO+TmcOqQnl58ykH49OiTcqO/fsOfnaCPeWqhAiLRhDQ3Okg92BAVhC2+u2sK+2gZy2hknDOjODZ8czoTiQk4Y0D20jlDJXCoQIm3Muqo9B/YQXluxhardNQAM792ZKacMZEJxIacO7UmX9nkRJ5WoqUCItHLb9tTw+ootvFy+mVfLN7Nmyx4Aencp4MwRRUwoLmR8cSF9uraPOKlkGhUIkVZmX209b63ZyitBQVi4fjvu0Lkgl7FDe3L1uMFMKC6kuHdn9RVIs1QgRLJcfYOz+P0dBwrC/NVVVNc1kNvOOGlgD7559lFMGN6L4/p3Jy8Lz8WX6KhAiGQZd2dt1Z4DBeG1FVvYtqcWgKOP6MKVYwcxobiQMUN60qlA/8Xl0Olfj0gW2LKrmtdWfHg9QsXWvQD07daec0f2YcLwQk4b1oveXdSPIOmjAiGSgfbW1PPm6qpYQVi+mcUbYoMxdmmfy7hhvfjy6UMZX1zIkMJO6keQ0KhAiGSA+gZn4frtvLK8klfKN7NgzTZq6hvIz2nHyYN68J3zRjC+uJDRR3bNynv6SHZSgRCJgLuzavPuj1yPsHNfHQCj+nbl6vGDGV9cyJjBPemQnxNxWmmrVCBEWkjlzmpeWxE7ZPRq+Wbe374PgH7dO3DhsX0ZX1zIuGG96NW5IOKkIjGhFQgzmwlcBGxy99EJpn8HuCIux0igKBhudDWwE6gH6ty9JKycImHZXV3Hm6uqDpxttPSDnQB075jHuGG9uL64kAnFhQzs2VH9CJKRwtyDuB+4C3gw0UR3/ynwUwAzuxj4VqNxp89y980h5hNJq9r6Bt6t2MYry2NnGy1Yu5W6Bic/tx1jBvfke5P6MaG4kFFHdiVHt5WWLBDmmNQvmdngFGe/HHg4rCwiYXB3yjftOrCHMG9lFbuq6zCDY/t1419PH8qE4kJOHtSD9nnqR5DsE3kfhJl1BCYBN8Q1O/CcmTlwr7vPaGb564DrAAYOHBhmVBE+2L6PV4OC8Er5ZjbtrAZgcK+OTD7hSCYUx65H6N4xP+KkIocv8gIBXAy82ujw0gR3X29mvYG5ZrbU3V9KtHBQPGYAlJSUHNz4giJJ7NxXyxsrP+xHWL5pFwA9O+UzblgvPjG8kHHDChnQs2PESUXSLxMKxBQaHV5y9/XBn5vM7HFgDJCwQIikU01dA2+v23agILy9bhv1DU77vHaMGdKLz5b0Z3xxISOP6KrhKaXVi7RAmFk34Azgyri2TkA7d98ZPJ8ITIsoorRy7s6yjTsPnHr6xqoq9tTU087guP7d+eoZwxhfXMhJg7pTkKt+BGlbwjzN9WHgTKDQzCqAqUAegLvfE8z2KeA5d48f6LYP8Hhw2l8u8JC7PxtWTml73t+298AewqvlW9i8K9aPMLSoE5edHNtDGDu0F906aMAcadvCPIvp8hTmuZ/Y6bDxbSuB48NJJW3R9r21vB7c6O7V8s2s3Bz7PVLYuYAJxb0YHwyYc2T3DhEnFcksmdAHIZJW1XWxAXP2j7O8sGIbDQ4d83MYO7QXVwS3wz6qjwbMEWmOCoRkvYYGZ8kHOw4UhDdXbWFfbQM57YwTB3Tna58czoThhRzfvzv5ubrRnUiqVCAkK62r2vORG91V7a4B4Kg+nbl8zMADA+Z0aa9+BJFDpQIhWWHr7hpeX7nlQOfymi17AOjTtYAzRxQxIehH6NNVA+aIpIsKhGSkhgbntRUfFoRF72/HHToX5DJ2aC++OG4wE4YXMqxI/QgiYVGBkIw07enF3P/aavJyjBMH9uBb5xzF+OJCju/fTQPmiLQQFQjJOO9WbOOB11dz+ZgB3HThKDoV6J+pSBT0U0wySn2Dc9MTiyjqXMAPLhip4iASIRUIySgPvbmWdyu2c9NFo3QGkkjEVCAkY1TurOb2Z5cyobiQi4/rG3UckTZPBUIyxo+eWUJ1bQPTJh+jM5NEMoAKhGSEeSu38JcF6/nyGUMZWtQ56jgiggqEZICaugZufmIRA3p24PqziqOOIyIBnSIikZv56iqWb9rFfVeVaOxmkQyiPQiJ1Ppte/nF88uZOKoPZ4/sE3UcEYmjAiGRmvZUGQBTLzkm4iQi0lhoBcLMZprZJjNb1MT0M81su5m9HTxuiZs2ycyWmVm5md0YVkaJ1t+XbmRO2Ua+cc5w+mmwHpGME+YexP3ApCTzvOzuJwSPaQBmlgP8CjgfGAVcbmajQswpEdhbU8/UWWUM792ZL40fEnUcEUkgtALh7i8BVYew6Big3N1XunsN8AgwOa3hJHJ3v1jOuqq93HbpaA3iI5Khov6feZqZvWNmz5jZ/oPQ/YB1cfNUBG0Jmdl1ZlZqZqWVlZVhZpU0WVG5i3v/sZJPn9iPsUN7RR1HRJoQZYFYAAxy9+OBXwJPHMqbuPsMdy9x95KioqJ05pMQuDtTnyyjIK8d379gZNRxRKQZkRUId9/h7ruC57OBPDMrBNYDA+Jm7R+0SSvw9LsbeKV8M989bwRFXQqijiMizYisQJjZERbccMfMxgRZtgDzgeFmNsTM8oEpwKyockr67NxXy21PL+a4/t343KmDoo4jIkmEdiW1mT0MnAkUmlkFMBXIA3D3e4DLgK+aWR2wF5ji7g7UmdkNwBwgB5jp7mVh5ZSWc8fc5VTuqua3V5WQ00434xPJdKEVCHe/PMn0u4C7mpg2G5gdRi6JRtn727n/tVVceeogjuvfPeo4IpKCqM9ikjagIRglrmenfP5j4oio44hIilQgJHR/LF3HP9du4wcXjKRbR40SJ5ItVCAkVFW7a/jxs0sZM6QnnzqxyctZRCQDqUBIqH7yzFJ27atj+qWjNUqcSJZRgZDQvLWmikdL13HNJ4ZwVJ8uUccRkYOkAiGhqKtv4D8fX8SR3drzjbOHRx1HRA6BRpSTUNz/2mqWfrCTez9/Mh3z9c9MJBtpD0LS7oPt+7hj7nt88ujeTBylUeJEspUKhKTdbX9dTF2Dc+vFx6hjWiSLqUBIWr30XiV/fXcDN5xVzMBeHaOOIyKHQQVC0mZfbT23PLmIIYWduO6MoVHHEZHDpN5DSZsZL61k9ZY9/P6aUynIzYk6jogcJu1BSFqs2bKbu14o5+Ljj2TC8MKo44hIGqhAyGFzd6bOKiM/px03XahR4kRai2YPMZnZnSm8xw53vylNeSQLzSn7gBeXVXLLRaPo07V91HFEJE2S9UFMBm5JMs+NgApEG7W7uo4fPrWYkX278oXTNEqcSGuSrEDc4e4PNDeDmfVIYx7JMnf+bTkbtu/jrs+dRG6OjliKtCbN/o92958ne4Om5jGzmWa2ycwWNTH9CjN718wWmtlrZnZ83LTVQfvbZlaaLINEY9kHO7nvlVVMOWUAJw/S7wSR1ialn3xmdruZdTWzPDP7m5lVmtmVSRa7H5jUzPRVwBnufixwGzCj0fSz3P0Edy9JJaO0LHfn5icW0aV9Lt+bdHTUcUQkBKkeE5jo7juAi4DVQDHwneYWcPeXgKpmpr/m7luDl/OA/ilmkQzw5wXreXN1Fd8/fyQ9OuVHHUdEQpBqgdjfV3Eh8Cd3357mHNcAz8S9duA5M3vLzK5rbkEzu87MSs2stLKyMs2xJJFte2r40ewlnDyoB5edrLou0lqleiX102a2FNgLfNXMioB96QhgZmcRKxAT4ponuPt6M+sNzDWzpcEeyce4+wyCw1MlJSWejkzSvJ/OWca2vbVMv3Q07drpZnwirVVKexDufiMwDihx91pgN7FTYA+LmR0H/BaY7O5b4j5vffDnJuBxYMzhfpakx9vrtvHQm2u5etxgRvbtGnUcEQlRSnsQZpZD7Bf+YDOLX+Znh/rBZjYQ+AvweXd/L669E9DO3XcGzycC0w71cyR96hucm55YSO8uBXzzHI0SJ9LapXqI6Slih5QWAg2pLGBmDwNnAoVmVgFMBfIA3P0eYhfg9QLuDsYMqAvOWOoDPB605QIPufuzKeaUEP3hjTUsWr+Duz53Il3a50UdR0RClmqB6O/uxx3MG7v75UmmXwtcm6B9JXD8x5eQKG3auY+fzlnGJ4YXcuGxfaOOIyItINWzmJ4xs4mhJpGM9t9/XUJ1bQPTJo/WKHEibUSqexDziB32aQfUAga4u6uXsg14bcVmnnj7fb5+9nCGFHaKOo6ItJBUC8TPgNOAhe6uU0nbkJq6Bm5+YhEDe3bk384cFnUcEWlBqRaIdcAiFYe257evrGRF5W7+74un0D5Po8SJtCWpFoiVwItm9gxQvb/R3Q/5NFfJfBVb93Dn35Yz6ZgjOGtE76jjiEgLS7VArAoe+cFD2oAfPrWYdmbccvGoqKOISARSKhDu/kMAM+scvN4VZiiJ3vOLNzJ38Ua+f/7RHNm9Q9RxRCQCqd7ue7SZ/RMoA8qCm+gdE240icremnpufaqMo/p05ksThkQdR0Qikup1EDOAb7v7IHcfBPw78JvwYkmU7nphORVb9zL90mPJ0yhxIm1Wqv/7O7n7C/tfuPuLgE6Ib4XKN+1ixksr+cxJ/RkzpGfUcUQkQimfxWRmNwO/C15fSezMJmlF3J1bnlxEh7wcvn+BRokTaetS3YP4ElBE7O6rfwYKgS+GFUqiMeud93ltxRa+O+loCjsXRB1HRCKW6h7EOe7+9fgGM/ss8Kf0R5Io7NhXy/S/LuH4/t24fMzAqOOISAZIdQ/i+ym2SZb62XPvsWVXNdMvPZYcjRInIiTZgzCz84ELgH5mdmfcpK5AXZjBpOUsWr+dB19fzefHDuLY/t2ijiMiGSLZIab3gVLgEuCtuPadwLfCCiUtp6HB+c8nFtGzUwHfnjgi6jgikkGaLRDu/g7wjpk9FIxFfVDMbCZwEbDJ3UcnmG7AL4jtpewBrnb3BcG0q4Cbglmnu/sDB/v5ktwj89fxzrpt/PxfTqBbB40SJyIfSrUPYrCZPWZmi81s5f5HCsvdD0xqZvr5wPDgcR3wawAz60lsiNJTgTHAVDPrkWJWSdGWXdX85NmljB3ak8knHBl1HBHJMKkWiP8jtvGuA84CHgR+n2whd38JqGpmlsnAgx4zD+huZn2B84C57l7l7luBuTRfaOQQ/PiZpeyurmP6pRolTkQ+LtUC0cHd/waYu69x91uBC9Pw+f2IjTWxX0XQ1lT7x5jZdWZWamallZWVaYjUNsxfXcWf3qrgX08fSnHvLlHHEZEMlGqBqA6GG11uZjeY2aeAziHmSpm7z3D3EncvKSoqijpOVqitb+CmxxfRr3sHvvbJ4qjjiEiGSrVAfAPoCHwdOJnYrTauSsPnrwcGxL3uH7Q11S5pcP+rq1m2cSe3XnIMHfNTvVZSRNqalAqEu893913uXuHuX3T3zwR9BodrFvAFixkLbHf3DcAcYKKZ9Qg6pycGbXKYNmzfyx3Pv8c5I3tz7qg+UccRkQzWbIEws1uTvUFz85jZw8DrwAgzqzCza8zsK2b2lWCW2cRu+ldO7Pbh/wbg7lXAbcD84DEtaJPDdNvTi2lwZ+rFGs5DRJqX7PjCtWa2o5npBkwBbk000d0vb+7N3d2B65uYNhOYmSSfHIQXl21i9sIP+M55IxjQs2PUcUQkwyUrEL8Bkp3iooGDssC+2nqmzipjaFEnrv2ERokTkeSSXUn9w5YKIuG65x8rWLNlD3+49lQKcnOijiMiWUDjSbYBqzfv5u4XV3DJ8Ucyvrgw6jgikiVUIFo5d+eWWWUU5LTjpgtHRh1HRLKICkQr98yiD3jpvUr+feJR9O7aPuo4IpJFkp3mWtjo9ZVmdmdwewvdvCfD7aquY9pTiznmyK5cOXZQ1HFEJMsk24N4bv8TM7sJ+DyxcSHOBX4WYi5Jg188/x4bd+5j+qWjyc3RzqKIHJxkp7nG7yV8GviEu+82s4eABeHFksO19IMdzHx1NVNOGciJA3WndBE5eMkKRAczO5HYnkaOu+8GcPdaM6sPPZ0ckoYG56bHF9GtQx7fPU+jxInIoUlWIDbw4aGkKjPr6+4bzKwXGpM6Yz22oILSNVu5/bLj6NEpP+o4IpKlkl0od1YTk7YBp6c9jRy2rbtr+NHsJZwyuAeXndQ/6jgiksWS9lya2RFmdkTwvMjMPg0c7e57Qk8nB+32OcvYsa+O2y4dTbt2OtFMRA5dstNcv0zsbqzzzOyrwNPERpL7i5ld0wL55CAsWLuVR+av5UvjB3P0EV2jjiMiWS5ZH8QNwDFAB2ANUOzuHwRjNLwA3BdyPklRXX0DNz+xiD5d2vONc46KOo6ItALJCkRtcChpj5mtcPcPANx9q5l5+PEkVb+ft4ay93dw9xUn0blAo8SJyOFL1gfhZpYXPL9wf6OZtU9hWWkhm3bs43+fe4/Tjyri/NFHRB1HRFqJZBv5TwEO4O4Vce29gH9P9uZmNsnMlplZuZndmGD6HWb2dvB4z8y2xU2rj5s2K6W1aaOm/3UJ1fUNTLvkGHQHFBFJl2Snua5ton09sL65Zc0sB/gVsdtyVADzzWyWuy+Oe59vxc3/NeDEuLfY6+4nJFuBtu7V8s3Meud9vnnOcAYXdoo6joi0Iod8mMjMFiaZZQxQ7u4r3b0GeASY3Mz8lwMPH2qetqi6rp6bn1zEoF4d+coZw6KOIyKtTLN7EME1DwknAckOdvcD1sW9rgBObeJzBgFDgL/HNbc3s1JiV2z/2N2fSPJ5bc5vX17Fysrd3P/FU2ifp1HiRCS9kp3u8ijwB4J+iEbSObjAFOAxd4+/v9Mgd19vZkOBv5vZQndf0XhBM7sOuA5g4MCBaYyU2dZV7eHOvy3ngmOP4MwRvaOOIyKtULIC8S7wP+6+qPEEMzsnybLrgQFxr/vTdL/FFOD6+IagnwN3X2lmLxLrn/hYgXD3GcAMgJKSkjZz6u0Pnyojp51x80Wjoo4iIq1Usj6IbwI7mpj2qSTLzgeGm9kQM8snVgQ+djaSmR0N9CB2xfb+th5mVhA8LwTGA4sbL9tWzV28keeXbOJb5xxF324doo4jIq1UsrOYXm5mWmmSZevM7AZgDpADzHT3MjObBpS6+/5iMQV4xN3jf/2PBO41swZiRezH8Wc/tWV7auq4dVYZI/p04erxg6OOIyKtWEqX3JrZUcCvgT7uPtrMjgMucffpzS3n7rOB2Y3abmn0+tYEy70GHJtKtrbml38vZ/22vfzpK6eRp1HiRCREqW5hfgN8H6gFcPd3if3ylxa0fONOfvPSSj57cn9OGdwz6jgi0sqlWiA6uvubjdo0YFALcndufnIRnQpyufH8o6OOIyJtQKoFYrOZDSM43dXMLiM22py0kCfffp95K6v43qSj6dW5IOo4ItIGpHrbz+uJnUp6tJmtB1YBV4SWSj5i+95apv91CccP6M6UUwYkX0BEJA2SDRg0EGLXIrj7OUARsdHkJgBt56q0iP3vc8uo2l3Nf2mUOBFpQckOMb1oZt8NbryHu+8GOprZ74E7Qk8nLKzYzu/mreELpw1mdL9uUccRkTYkWYE4GRgGvG1mnzSzbwBvEruobUzY4dq6+gbnpicWUti5gG9P1ChxItKykl0otxX4clAYngfeB8Y2GhtCQvLwm2t5p2I7v5hyAl3b5yVfQEQkjZL1QXQ3s3uBLwKTgMeAZ8zsky0Rri3bvKua259dyrhhvbjk+COjjiMibVCys5gWAHcD17t7HfCcmZ0A3G1ma9z98rADtlU/mr2UvbX1TJs8WqPEiUgkkhWI0xsfTnL3t4FxZvavoaVq495YuYU/L6jg+rOGUdy7c9RxRKSNavYQU3N9De7+m/THkdr6Bm5+chH9unfghrOGRx1HRNqwVC+UkxYy85VVvLdxF7/9Qgkd8jVKnIhER7cDzSDvb9vLz59fzrmj+nDOqD5RxxGRNk4FIoNMe2oxjjP1Yo0SJyLRU4HIEC8s3cSzZR/w9bOH079Hx6jjiIioQGSCfbX1TJ1VxrCiTlw7YWjUcUREgJALhJlNMrNlZlZuZjcmmH61mVWa2dvB49q4aVeZ2fLgcVWYOaN294srWFu1h9suHU1+rmq2iGSG0M5iCm7w9yvgXKACmG9msxKMLf2ou9/QaNmewFSghNgYFG8Fy24NK29UVlbu4p4XV3DpCUcyblhh1HFERA4I8+fqGKA8uFV4DfAIMDnFZc8D5rp7VVAU5hK71Uer4u7c8mQZBXnt+MGFI6OOIyLyEWEWiH7AurjXFUFbY58xs3fN7DEz2z8aTqrLYmbXmVmpmZVWVlamI3eL+evCDbxSvpnvnDeC3l3aRx1HROQjoj7g/RQw2N2PI7aX8MDBvoG7z3D3EncvKSoqSnvAsOzcV8u0pxYzul9Xrjh1UNRxREQ+JswCsR6IHx+zf9B2gLtvcffq4OVviY0/kdKy2e7nzy+nclc10y89lhyNEiciGSjMAjEfGG5mQ8wsH5gCzIqfwcz6xr28BFgSPJ8DTDSzHmbWA5gYtLUKi9/fwf2vreZzYwZywoDuUccREUkotLOY3L3OzG4gtmHPAWa6e5mZTQNK3X0W8HUzuwSoA6qAq4Nlq8zsNmJFBmCau1eFlbUlNTQ4Nz+5iO4d8vjueUdHHUdEpEnm7lFnSJuSkhIvLS2NOkazHp2/lu/9eSH/89njuezk/lHHEZE2zszecveSRNOi7qRuU6p21/CjZ5YyZnBPPnNSwpOyREQyhgpEC7r92aXs2lfHbZdqlDgRyXwqEC3krTVbeWT+Oq6ZMIQRR3SJOo6ISFIqEC2grr6Bm55YRN9u7fn62RolTkSygwpEC3jw9TUs2bCDqRePolOBBvETkeygAhGyjTv28bO573HmiCLOO+aIqOOIiKRMBSJktz29mJr6Bn54yTHqmBaRrKICEaKXl1fy9LsbuP7MYgb16hR1HBGRg6ICEZLqunpuebKMIYWd+PIZGiVORLKPekxDMuMfK1m1eTcPfmkM7fNyoo4jInLQtAcRgrVb9nDXC+VceFxfTj8qe25BLiISTwUizdydqbMWkdvOuPnCUVHHERE5ZCoQaTanbCMvLKvkW+cexRHdNEqciGQvFYg02l1dx7Snyjj6iC5cPW5w1HFERA6LCkQa3fn35by/fR/TLx1Nbo7+akUku2krlibvbdzJfS+v4l9KBlAyuGfUcUREDluoBcLMJpnZMjMrN7MbE0z/tpktNrN3zexvZjYoblq9mb0dPGY1XjaTuDs3PbGIzu1z+d75GiVORFqH0K6DMLMc4FfAuUAFMN/MZrn74rjZ/gmUuPseM/sqcDvwL8G0ve5+Qlj50unxf67nzVVV/PjTx9KzU37UcURE0iLMPYgxQLm7r3T3GuARYHL8DO7+grvvCV7OA7JuDM7te2r579lLOHFgd/5fyYCo44iIpE2YBaIfsC7udUXQ1pRrgGfiXrc3s1Izm2dmlza1kJldF8xXWllZeViBD8VPn1tK1e4apl86mnbtdDM+EWk9MuJWG2Z2JVACnBHXPMjd15vZUODvZrbQ3Vc0XtbdZwAzAEpKSrxFAgfeWbeNP7yxlqvHDeaYI7u15EeLiIQuzD2I9UD8MZf+QdtHmNk5wH8Cl7h79f52d18f/LkSeBE4McSsB62+IdYxXdS5gG+fe1TUcURE0i7MAjEfGG5mQ8wsH5gCfORsJDM7EbiXWHHYFNfew8wKgueFwHggvnM7cg+9sYaF67dz80Wj6NI+L+o4IiJpF9ohJnevM7MbgDlADjDT3cvMbBpQ6u6zgJ8CnYE/BYPprHX3S4CRwL1m1kCsiP240dlPkarcWc3tc5YxobiQi47rG3UcEZFQhNoH4e6zgdmN2m6Je35OE8u9BhwbZrbD8aPZS6iubWDaZI0SJyKtl66kPkjzVm7hL/9cz5fPGMrQos5RxxERCY0KxEGoqWvgpicWMaBnB64/qzjqOCIiocqI01yzxX2vrKJ80y7uu6pEo8SJSKunPYgUVWzdw51/W87EUX04e2SfqOOIiIROBSJF056KnUQ19ZJjIk4iItIyVCBS8LclG3lu8Ua+cc5w+nXvEHUcEZEWoQKRxN6aeqbOKmN47858afyQqOOIiLQYdVIn8asXyqnYupdHrhtLfq7qqYi0HdriNWNF5S7ufWkFnz6xH2OH9oo6johIi1KBaIK7c8uTi2ifl8P3LxgZdRwRkRanAtGEp97dwKvlW/jueSMo6lIQdRwRkRanApHAjn213Pb0Yo7r343PnToo+QIiIq2QOqkTuGPue2zeVc19V5WQo1HiRKSN0h5EI2Xvb+eB11Zz5amDOK5/96jjiIhERgUiTkMwSlzPTvn8x8QRUccREYmUCkScR0vX8c+12/jBBSPp1lGjxIlI26YCEdiyq5ofP7OUMUN68qkT+0UdR0QkcqEWCDObZGbLzKzczG5MML3AzB4Npr9hZoPjpn0/aF9mZueFmRPgJ88uZXd1HdMvHa1R4kRECLFAmFkO8CvgfGAUcLmZjWo02zXAVncvBu4AfhIsOwqYAhwDTALuDt4vFKWrq/hjaQXXfmIoR/XpEtbHiIhklTD3IMYA5e6+0t1rgEeAyY3mmQw8EDx/DDjbYj/fJwOPuHu1u68CyoP3S7u6+tgocUd2a8/Xz9YocSIi+4V5HUQ/YF3c6wrg1Kbmcfc6M9sO9Ara5zVaNmHHgJldB1wHMHDgwIMOua+ugWP7deOcUX3omK/LQkRE9sv6LaK7zwBmAJSUlPjBLt+5IJeffvb4tOcSEcl2YR5iWg8MiHvdP2hLOI+Z5QLdgC0pLisiIiEKs0DMB4ab2RAzyyfW6Tyr0TyzgKuC55cBf3d3D9qnBGc5DQGGA2+GmFVERBoJ7RBT0KdwAzAHyAFmunuZmU0DSt19FnAf8DszKweqiBURgvn+CCwG6oDr3b0+rKwiIvJxFvvB3jqUlJR4aWlp1DFERLKGmb3l7iWJpulKahERSUgFQkREElKBEBGRhFQgREQkoVbVSW1mlcCaQ1y8ENicxjhRai3r0lrWA7Qumai1rAcc3roMcveiRBNaVYE4HGZW2lRPfrZpLevSWtYDtC6ZqLWsB4S3LjrEJCIiCalAiIhIQioQH5oRdYA0ai3r0lrWA7Qumai1rAeEtC7qgxARkYS0ByEiIgmpQIiISEJtrkCY2SQzW2Zm5WZ2Y4LpBWb2aDD9DTMbHEHMpFJYj6vNrNLM3g4e10aRMxkzm2lmm8xsURPTzczuDNbzXTM7qaUzpiqFdTnTzLbHfSe3tHTGVJnZADN7wcwWm1mZmX0jwTwZ/92kuB5Z8b2YWXsze9PM3gnW5YcJ5knv9svd28yD2G3HVwBDgXzgHWBUo3n+DbgneD4FeDTq3Ie4HlcDd0WdNYV1OR04CVjUxPQLgGcAA8YCb0Sd+TDW5Uzg6ahzprgufYGTguddgPcS/BvL+O8mxfXIiu8l+HvuHDzPA94AxjaaJ63br7a2BzEGKHf3le5eAzwCTG40z2TggeD5Y8DZZmYtmDEVqaxHVnD3l4iNBdKUycCDHjMP6G5mfVsm3cFJYV2yhrtvcPcFwfOdwBI+Pi58xn83Ka5HVgj+nncFL/OCR+OzjNK6/WprBaIfsC7udQUf/8dyYB53rwO2A71aJF3qUlkPgM8Eu/6PmdmABNOzQarrmi1OCw4RPGNmx0QdJhXBYYoTif1ijZdV300z6wFZ8r2YWY6ZvQ1sAua6e5PfSTq2X22tQLQlTwGD3f04YC4f/qqQ6Cwgdt+b44FfAk9EGyc5M+sM/Bn4prvviDrPoUqyHlnzvbh7vbufAPQHxpjZ6DA/r60ViPVA/C/p/kFbwnnMLBfoBmxpkXSpS7oe7r7F3auDl78FTm6hbOmWyneWFdx9x/5DBO4+G8gzs8KIYzXJzPKIbVT/4O5/STBLVnw3ydYj274XAHffBrwATGo0Ka3br7ZWIOYDw81siJnlE+vEmdVonlnAVcHzy4C/e9Djk0GSrkejY8GXEDv2mo1mAV8IzpgZC2x39w1RhzoUZnbE/uPBZjaG2P+/TPvxAcTOUCI2ZvwSd/9ZE7Nl/HeTynpky/diZkVm1j143gE4F1jaaLa0br9yD3XBbOTudWZ2AzCH2JlAM929zMymAaXuPovYP6bfmVk5sQ7HKdElTizF9fi6mV0C1BFbj6sjC9wMM3uY2FkkhWZWAUwl1vmGu98DzCZ2tkw5sAf4YjRJk0thXS4DvmpmdcBeYEoG/vjYbzzweWBhcMwb4AfAQMiq7yaV9ciW76Uv8ICZ5RArYn9096fD3H7pVhsiIpJQWzvEJCIiKVKBEBGRhFQgREQkIRUIERFJSAVCREQSUoEQEZGEVCBEGjGzXnG3fv7AzNYHz3eZ2d0hfN79ZrbKzL7SxPRdwZ/D9udIdwaRRHQdhEgzzOxWYJe7/0+In3E/sdtNP9bE9F3u3rmp1yJh0R6ESIqCgWWeDp7famYPmNnLZrbGzD5tZreb2UIzeza4/w9mdrKZ/cPM3jKzOancDju4hcrrwXtND3u9RJqiAiFy6IYBnyR2r6vfAy+4+7HEbtdwYVAkfglc5u4nAzOB/0rhfX8B/Dp4r4y6t5G0LW3qXkwiafaMu9ea2UJi98R6NmhfCAwGRgCjgbnBveBySG2DPx74TPD8d8BP0phZJGUqECKHrhrA3RvMrDbuBm8NxP5vGVDm7qcdwnurc1Aip0NMIuFZBhSZ2WkQG5cgxdHKXuXDu3BeEVY4kWRUIERCEowXfhnwEzN7B3gbGJfCot8Arg8OXWXsEJ7S+uk0V5GIJTvNNcH8Os1VWoT2IESitx24rakL5fbbf6EcsLFFUkmbpz0IERFJSHsQIiKSkAqEiIgkpAIhIiIJqUCIiEhC/x8YmTbUSw1O7gAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"results = openmc.deplete.ResultsList.from_hdf5('depletion_results.h5')\n",
"time, n_Xe135 = results.get_atoms('1', 'Xe135')\n",
"days = 24*60*60\n",
"plt.plot(time/days, n_Xe135)\n",
"plt.xlabel('Time [d]')\n",
"plt.ylabel('Xe135 [atoms]')"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "a7acf699",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.7"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment