Evaluation of PBMR-400 Core Design Steady State Condition with Serpent and AGREE

The significant recent advances in computer speed and memory have made possible an increasing fidelity and accuracy in reactor core simulation with minimal increase in the computational burden. This has been important for modeling some of the smaller advanced reactor designs for which simplified approximations such as few groups homogenized diffusion theory are not as accurate as they were for large light water reactor cores. For narrow cylindrical cores with large surface to volume ratios such the Ped Bed Modular Reactor (PBMR), neutron leakage from the core can be significant, particularly with the harder neutron spectrum and longer mean free path than a light water reactor. In this paper the core from the OECD PBMR-400 benchmark was analyzed using multigroup Monte Carlo cross sections in the HTR reactor core simulation code AGREE. Homogenized cross sections were generated for each of the discrete regions of the AGREE model using a full core SERPENT Monte Carlo model. The cross sections were generated for a variety of group structures in AGREE to assess the importance of finer group discretization on the accuracy of the core eigenvalue and flux predictions compared to the SERPENT full core Monte Carlo solution. A significant increase in the accuracy was observed by increasing the number of energy groups, with as much as a 530 pcm improvement in the eigenvalue calculation when increasing the number of energy groups from 2 to 14. Significant improvements were also observed in the AGREE neutron flux distributions compared to the SERPENT full core calculation.

One-step coupled calculations (Serpent-OpenFOAM) for a fuel rod of the IPR-R1 triga reactor

In this work, a single step of coupled calculations for a fuel rod of IPR-R1 TRIGA was performed. The used me-thodology allowed to simulate the fuel pin behavior in steady-state mode for different power levels. The aim of this paper is to present a practical approach to perform coupled calculations between neutronic (Monte Carlo) and thermal-hydraulic (CFD) codes. For this purpose, is necessary to evaluate the influence of the water thermal-physical properties temperature variations on keff parameter. Besides that, Serpent Nuclear Code was used for the neutronics evaluation, while OpenFOAM was used for thermal-hydraulics. OpenFOAM si- mula-tions were made by using a modified chtMultiRegionFoam solver, developed to read Serpent output correctly. The neutronic code was used without any modifications. The results shows that this coupled calculations were consistent and that leads to encouraging further methodology development and its use for full core simulation. Also, the results shows good agreement with calculations performed using other version of OpenFOAM and Milonga as neutronic code.

Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations

Decay calculations play an important role in reactor physics simulations. In particular, the isotopic decay of the burned fuel during refueling is important for predicting the startup reactivity of the following burn-up cycle. In addition, there is a growing interest in high-fidelity simulations where the mesh in the burn-up region can involve millions of regions. However, existing models repeatedly solve the same Bateman equations for each region, which is a waste of calculational resources. RMC is a Monte Carlo neutron transport code developed for advanced reactor physics analysis including criticality calculations and burn-up calculations. This paper presents a decay calculation method named the Decay Chain Method (DCM) to optimize the RMC code for large-scale decay calculations. Unlike traditional methods, the Decay Chain Method solves the Bateman equations one decay chain at a time rather than one region at a time. The decay calculation in the burn-up mode then treats the decay steps as zero power burn-up steps with some optimized calculational methods to further reduce the calculational time. These methods were evaluated for a single pin example and for a Virtual Environment for Reactor Applications (VERA) full-core example. The calculations for the single pin example verify the accuracy of the decay step treatment in the burn-up mode and show the improved efficiency. The single pin is divided into 1–1,000,000 decay regions to study the efficiency differences between the Transmutation Trajectory Analysis (TTA) and DCM methods. Both methods have a linear complexity with respect to the number of regions but DCM costs just one-sixtieth of the TTA time. In the simplified VERA full core example, the DCM method reduces the decay calculation time to 0.32 min from 75.26 min while the accuracy remains unchanged.