crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Update pincell depletion example plots

Browse source
This commit is contained in:
Paul Romano 2022-01-09 15:47:08 -06:00
parent 141a32487a
commit 2ccbf0669b
2 changed files with 24 additions and 26 deletions
Download patch file Download diff file Expand all files Collapse all files

25
examples/pincell_depletion/restart_depletion.py
View file

@ -72,21 +72,20 @@ time, Xe_capture = results.get_reaction_rate('1', 'Xe135', '(n,gamma)')
###############################################################################
days = 24*60*60
plt.figure()
plt.plot(time/days, keff[:, 0], label="K-effective")
plt.xlabel("Time (days)")
plt.ylabel("Keff")
fig, ax = plt.subplots()
ax.errorbar(time/days, keff[:, 0], keff[:, 1], label="K-effective")
ax.set_xlabel("Time [d]")
ax.set_ylabel("Keff")
plt.show()
plt.figure()
plt.plot(time/days, n_U235, label="U 235")
plt.xlabel("Time (days)")
plt.ylabel("n U5 (-)")
fig, ax = plt.subplots()
ax.plot(time/days, n_U235, label="U235")
ax.set_xlabel("Time [d]")
ax.set_ylabel("U235 atoms")
plt.show()
plt.figure()
plt.plot(time/days, Xe_capture, label="Xe135 capture")
plt.xlabel("Time (days)")
plt.ylabel("RR (-)")
fig, ax = plt.subplots()
ax.plot(time/days, Xe_capture, label="Xe135 capture")
ax.set_xlabel("Time [d]")
ax.set_ylabel("Xe135 capture rate")
plt.show()
plt.close('all')

25
examples/pincell_depletion/run_depletion.py
View file

@ -117,21 +117,20 @@ time, Xe_capture = results.get_reaction_rate('1', 'Xe135', '(n,gamma)')
###############################################################################
days = 24*60*60
plt.figure()
plt.plot(time/days, keff[:, 0], label="K-effective")
plt.xlabel("Time (days)")
plt.ylabel("Keff")
fig, ax = plt.subplots()
ax.errorbar(time/days, keff[:, 0], keff[:, 1], label="K-effective")
ax.set_xlabel("Time [d]")
ax.set_ylabel("Keff")
plt.show()
plt.figure()
plt.plot(time/days, n_U235, label="U235")
plt.xlabel("Time (days)")
plt.ylabel("n U5 (-)")
fig, ax = plt.subplots()
ax.plot(time/days, n_U235, label="U235")
ax.set_xlabel("Time [d]")
ax.set_ylabel("U235 atoms")
plt.show()
plt.figure()
plt.plot(time/days, Xe_capture, label="Xe135 capture")
plt.xlabel("Time (days)")
plt.ylabel("RR (-)")
fig, ax = plt.subplots()
ax.plot(time/days, Xe_capture, label="Xe135 capture")
ax.set_xlabel("Time [d]")
ax.set_ylabel("Xe135 capture rate")
plt.show()
plt.close('all')