Continuing to investigate convergence

Related to my previous question I’d like to dig a little deeper. I’ve taken some time to read the OpenMC documentation in more detail.

For my work, I think it could be interesting to inspect the convergence rates of k_{eff} and the Shannon entropy of the source distribution.

In addition to this, the section on nonlinear diffusion acceleration (NDA) sounds interesting too.

I have the following questions:

  • When I run some of the example code, I can see that the outputs contain k_{eff}, but what should the values converge to? Is there a built-in facility to measure the convergence rate? How can I output the Shannon entropy of the source distribution as described in Chapter 8?
  • How can I output the result of the neutron balance equation described in Chapter 12?
  • What code would be most interesting to run convergence experiments on? Is any one of the example codes more suitable than the others?

I know I’m asking a lot of questions here, I’m happy about any pointers you can give me.