Critical boron concentration

Hello everyone. I’m trying to find critical boron concentrate in SMR core but I’m new to OpenMC. I have followed to use openmc.search_for_keff to find the critical boron concentration but when I change fuel enrichment to 2.6, it shows f(a) and f(b) must have different signs. What should I do if I would like to find it in SMR core which I have modeled similar to PWR?

This is my script.
6.py (16.3 KB)

Thank you.
Yos.

If f(a) and f(b) have the same sign, it’s possible that the function you’re finding a root for is not monotonic, or that you have not successfully bracketed the root you need to find.

I would bet on the latter being the case. You should check what keff is at 0 ppm boron. If it is supercritical, you should just change your bracket argument to be [0, 2500]. Otherwise, you won’t be able to find criticality using boron if the reactor is subcritical without boration to begin with.