Issue with assigning fuel volume in depletion

Hello,

I am having some issues with creating a depletion for my VVER model. when I try to run it in jupyter notebook it gives me an error of ‘kernel died’. If I run it in a Linux shell it gives me an error saying ‘no boundary conditions were applied to any surface.’ i know there is no issue with my geometry as the code runs without the depletion settings. I think the issue is around assigning the fuel.volume function.

I am attaching screenshots of the error

and the python code i used

Any help is appreciated!

#Create depletion “operator”
au13.volume=2286853.171
au24.volume=2000996.525
au33.volume=1705977.806
au36.volume=40459.70995
au40.volume=571713.2928
au44.volume=1099448.64
fuel_mix_1.volume=56804.8464
fuel_mix_2.volume=43977.9456

chain_file = openmc.deplete.Chain.from_xml(“./chain_endfb71_pwr.xml”)
chain_file.nuclide_dict
model = openmc.Model(geometry=geometry, settings=settings_file, materials=materials_file)
op = openmc.deplete.CoupledOperator(model, “./chain_endfb71_pwr.xml”)

Perform simulation using the predictor algorithm

time_steps = [60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60] # Days
power = 1200e6 # W/cm, for 2D simulations only (use W for 3D)
integrator = openmc.deplete.PredictorIntegrator(op, time_steps, power, timestep_units=‘d’)
integrator.integrate()

settings.export_to_xml()

I had kernels die on me during depletion because my cross section files were missing some data and the program wasn’t skipping past them. Fixed it by using the TENDL library as this had all the ones I needed. Idk if this is helpful but might be worth checking.

As for the geometry have you specified boundaries, like a killzone around your geometry?

So I know it’s an issue with my depletion code written above because my geometry seems to be working without the depletion code.

I think my issue is linked to the fuel volume but I will try to use another library and see if it resolves my issue