Fuel Assembly Burn-up calculation K-eff value

[Yiming Zhong]

Hi all,

I am trying to investigate a 17*17 PWR fuel assembly depletion problem. The depletion time is one refuelling cycle which is 18 months. The depletion step size was divided have a very short step (1 hour) at the beginning of the calculation, and a bit longer steps (one month)*18 at the later part of cycle time.  The results of K-eff values decrease but have some point waving. The k-eff will increase at some steps. 

I can't figure out why it looks like that. In my mind, the k-eff value should smoothly decrease.  The k-eff values in each step and plots would be attached in this mail. 

Please contact me if you have any idea why that could happen and how to optimise that.

Thanks,

Yiming

[Khurrum]

Hye Yiming

Have you tried any other integrator beside predictor integrator?

[Andrew Johnson]

Hi Yiming,

Looking at the uncertainties on your k value, I don't think this is purely a statistical artifact. I do find it interesting that your fuel assembly is pretty subcritical, but that shouldn't be causing this significant of an issue.

Without more information on your model, my initial reaction would be to check the convergence of your fission source for the fresh model using the Shannon entropy setting. I would also like to repeat Khurrum's comment on using other time-integration schemes (CE/CM, CE/LI, etc.) The predictor is conditionally stable, meaning sufficiently large time steps can cause instabilities and inaccuracies. If you are able to use one of the other integrators, please reach out with an update.

Regards,

Andrew

[Zhong Yiming]

Hi Andrew, Khurrum,

Thanks for your advice. I will try the other integrators and update the results here.

Thanks,

Yiming

Andrew Johnson <1drew.e...@gmail.com> 于2020年7月7日周二 下午12:23写道：

--
You received this message because you are subscribed to a topic in the Google Groups "OpenMC Users Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/openmc-users/mIMxbLR7zRs/unsubscribe.
To unsubscribe from this group and all its topics, send an email to openmc-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/openmc-users/eaa940e0-c324-4c52-9f40-065d627e8f0fn%40googlegroups.com.

[Zhong Yiming]

Hi,

After the above discussion, I rebuild the 17*17 nuclear assembly model geometry with 3D pitch universes Rectangular Lattices function. Now the results of eigenvalue calculation look good, even with the predictor integrator. I also used the CE/CM and CE/LI integrator as well. The results look very close. I put the calculation results in this email. Hope that could be helpful.

Meanwhile, I have another question. Due to the application of 'diff_burnable_mats' argument, there are many fuel materials in the fuel assembly model. But it brings some challenge to the analysis of the results. The first required parameter in the openmc.deplete.ReactionRates [Python API is the local_mats(list of str) Material ID. But there are maybe some materials listed in the material.XML file. Is there any method to get the reaction rate or atom number in all the burnable materials at the same time? Thanks for your advice.

Best,

Yiming

Zhong Yiming <zhongy...@gmail.com> 于2020年7月7日周二 下午3:25写道：