i wanna get the relationship between temperatures and keff, when i create statepoint.h5 related to temperature range, i met the error and i can’t solve it
error:FileNotFoundError: [Errno 2] Unable to open file (unable to open file: name = ‘/home/ww/work/examples/molten_salt_reactor/statepoint.973.h5’, errno = 2, error message = ‘No such file or directory’, flags = 0, o_flags = 0)
here is my last code:
def simulate_keff(temperatures):
keff_values={}
materials_dict= create_materials(temperatures)
cells = create_cells(materials_dict, temperatures)
# universes = create_universes(cells)
# print(universes)
for temp, cells_list in cells.items():
lattice_outer=openmc.Universe(cells=[cells_list[15]])
lattice_inner=openmc.Universe(cells=cells_list[:3])
final_universe=create_lattice_and_universe(lattice_outer,lattice_inner,cells_list)
geometry = openmc.Geometry(root=final_universe)
geometry.export_to_xml()
batches=100
inactive=10
particles=1000
settings=openmc.Settings()
settings.batches=batches
settings.inactive=inactive
settings.particles=particles
settings.run_mode='eigenvalue'
settings.temperature['method']='interpolation'
uniform_dist=openmc.stats.Box([-24,-24,-50],[24,24,50],only_fissionable=True)
settings.source=openmc.Source(space=uniform_dist)
settings.export_to_xml()
statepoint_file = '/home/ww/work/examples/molten_salt_reactor/statepoint.{}.h5'.format(temp)
openmc.run(output={'statepoint': statepoint_file})
with openmc.StatePoint(statepoint_file) as sp:
keff_values[temp] = sp.k_combined
return keff_values
temperatures = [973,983,993,1003,1013,1023,1033,1043,1053,1063,1073,1083,1093,1100]
keff_values = simulate_keff(temperatures)
print(keff_values)