crp/MSRE
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

minor

Browse source
This commit is contained in:
church89 2022-11-18 10:32:12 +01:00
parent f46c6c718e
commit 568851e60a
1 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
msre.py
View file

@ -398,7 +398,7 @@ def depletion(model, mass, power):
def control_rod_worth(model):
rho = []
cell = model.geometry.get_cells_by_name('CR1')[0]
for inch in np.arange(0,51,10):
for inch in np.arange(0,51,5):
cm = inch*2.54
setattr(cell, 'translation', [0,0,19.2+cm])
res=model.run()
@ -406,10 +406,10 @@ def control_rod_worth(model):
keff=sp.keff.n
pcm = (keff-1)/keff*1e5
rho.append(pcm)
dx = 2*2.54
dx = 5*2.54
d_rho = np.diff(rho)/dx
plt.figure()
plt.scatter(np.arange(0,51,10)[1:], d_rho)
plt.scatter(np.arange(0,51,5)[1:], d_rho)
plt.xlabel('Withdrawn of control rod n. 1 [inch]')
plt.ylabel('Reactivity worth')
plt.savefig('reac_rod_worth')
@ -480,5 +480,5 @@ if __name__ == '__main__':
power = 8e6 #total thermal power [W]
#run(build(make_tally=True, plot_geom=True), mass, power)
control_rod_worth(build(make_tally=False, plot_geom=False))
depletion(build(make_tally=False, plot_geom=False), mass, power)
#depletion(build(make_tally=False, plot_geom=False), mass, power)
#triton_adder(mass)