Fast reactor depletion caculate

User Support
1

When I calculate the burn-up of the reactor, the power is set to about 1e7, and the Keff will gradually decrease. However, when the set power value is less than 1e7, the keff change is not regular. How should I handle the situation, please? I would appreciate it if anyone answered the question.
Here’s my code:
model = openmc.Model(geometry=geom, settings=settings)
operator = openmc.deplete.CoupledOperator(model, “./chain_casl_sfr.xml”)
power = 1e6
time_steps = [30] * 12
integrator = openmc.deplete.LEQIIntegrator(operator, time_steps, power, timestep_units=‘d’)
integrator.integrate()

computational results
image

2

I have the same problem with you
:pensive: have you solved this problem?

3

Hi, my power does not match the core design, and when my power is scaled up, the problem is solved.