Burnup calculation problem

Hi,
I want to calculate up to a burnup of 20 MWd/kgHM. 50steps of 0.2MWd/kgHM & 10steps of 1.0MWd/kgHM. How can i do this successive depletion calculation? Will one integration after another do?

Thanks
-Imtiaj

Thank you julius.

openmc.deplete.SICELIIntegrator works well.

-Imtiaj

Please,disregard my previous comment.I suggest you do it this way

Depletion by burnup

Yeah thanks again. I was also facing problem cz it was a array not successive burnup. so i tried spending more time on openmc documentation. Thank u.

It sounds like you’ve already figured it out, but for posterity, the easiest way to do this would be:

timesteps = 50*[0.2] + 10*[1.0]
integrator = Integrator(operator, timesteps=timesteps, timestep_units='MWd/kg')

where Integrator is really one of the various integrator classes. As you’ve noted, the values are not cumulative.