We computed by OpenMC two 1000 MWt Sodium Fast Reactor Benchmarks (MET-1000 and MOX-1000) with metallic and oxide fuel from NEA/NSC/R(2015)9. For a comparison we use the average value of the benchmark participants and standard deviations of the average value. We have reasonably well agreement for reactor multiplication factor, Doppler effect and control rod worth. The difference is below the standard deviation of the average value.
However for the sodium void effect we have considerable difference from the average benchmark results, the deviation is about 60 % for the metallic fuel MET-1000 and 35 % for the oxide fuel MOX-1000. This benchmark is well know and I think may be someone already computed these reactors and can say anything about the difference. If someone is interested about these problems and want to check it I am ready to share my OpenMC input files for the both problems.
@Denis_Solovyev Thanks for sharing these results. I’d be interested in taking a look at the inputs to see if I can reproduce the differences in sodium void effect you’re seeing. If you have equivalent Serpent/MCNP inputs, that would be helpful as well.
Hello, Paul Romano.
We send you our codes for MET-1000 calculation.
met_1000_core_up.py — Here is the code, where the core is with sodium.
met_1000_core_up_void.py — Here is the code, where the core is without sodium
We use standard OpenMC library ENDF/B-VII.1 for our calculation ( Official Data Libraries | OpenMC ) with linear interpolation of the cross section.
We do not have equivalent Serpent/MCNP results.
If you have any questions, please, let me know.
Denis Solovyev.
chanel_3d.py (10.5 KB)
control_rods.py (6.4 KB)
control_rods_up_calculation.py (9.2 KB)
coolant.py (4.0 KB)
coolant_cell_construction.py (910 Bytes)
empty_duct.py (3.2 KB)
fa_calculation.py (12.9 KB)
fuel_assembly.py (4.9 KB)
materials_met_1000.py (25.6 KB)
met_1000_core_up.py (7.0 KB)
met_1000_core_up_void.py (6.9 KB)
radial_reflector.py (4.8 KB)
radial_reflector_calculation.py (7.8 KB)
radial_shield.py (4.8 KB)
radial_shield_calculation.py (7.8 KB)
Hello, Paul Romano.
We found the reason of the differences in the sodium void effect between our calculation and the average benchmark results.
The problem was that for voided states we did not correctly removed sodium from the reactor core.
When we changed it we get the result close to benchmark report.
We add the revised version of our input:
control_rods_up_calculation.py (9.2 KB)
fa_calculation.py (12.8 KB)
met_1000_core_up.py file is used for calculation reactor core with sodium.
met_1000_core_up_void.py file is used for calculation reactor core without sodium.
1 Like
Glad you were able to figure it out!