Search

The development of fast reactors with complex designs and operation status requires more accurate and effective simulation. The Monte-Carlo method can generate multi-group cross-sections in arbitrary geometry without approximation on resonances treatment and leads to good results in combination with diffusion codes. However, in previous studies, the coupling of Monte-Carlo generated multi-group cross-sections (MC-MGXS) and transport solvers has shown relatively large biases in fast reactor problems. In this paper, the main contribution to the biases is proved to be the neglect of the angle-dependence of the total cross-sections. The flux-moment homogenization technique (MHT) is proposed to take into account this dependence. In this method, the angular dependence is attributed to the transfer cross-sections, keeping an independent form for the total sections. For the MET-1000 benchmark, the multi-group transport simulation results with MC-MGXS generated with MHT are improved by 700 pcm and an additional 120 pcm with higher order scattering. The factors that cause the residual bias are discussed. The core power distribution bias is also significantly reduced when MHT is used. It proves that the MC-MGXS with MHT can be applicable with transport solvers in fast reactor analysis.

With the rise of economic and safety standards for nuclear reactors, new concepts of Gen-IV reactors and modular reactors showed more complex designs that challenge current tools for reactor physics analysis. A Monte Carlo (MC) two-step method was proposed in this work. This calculation scheme uses the continuous-energy MC method to generate multi-group cross-sections from heterogeneous models. The multi-group MC method, which can adapt locally-heterogeneous models, is used in the core calculation step. This calculation scheme is verified using a Gen-IV modular lead-based fast reactor (LFR) benchmark case. The influence of homogenized patterns, scatter approximations, flux separable approximation, and local heterogeneity in core calculation on simulation results are investigated. Results showed that the cross-sections generated using the 3D assembly model with a locally heterogeneous representation of control rods lead to an accurate estimation with less than 270 pcm bias in core reactivity, 0.5% bias in control rod worth, and 1.5% bias on power distribution. The study verified the applicability of multi-group cross-sections generated with the MC method for LFR analysis. The study also proved the feasibility of multi-group MC in core calculation with local heterogeneity, which saves 85% time compared to the continuous-energy MC.

Fuel assemblies have a decisive impact on the performance and safety of nuclear reactors. Helical fuel has huge potential for application in small module reactors (SMRs) due to its advantages in volume power density and safety. Typical helical fuel elements are usually trilobal or cruciform in cross-section. The fuel rods are helically twisted in the axial direction, eliminating the need for spacer grids as the fuel rods are self-supporting. In this paper, a refined subchannel division approach is proposed based on the crossflow mechanism of helical fuel assemblies. Then, a refined helical fuel mixing model framework, including the helical fuel distributed resistance method and directed crossflow method, is developed and implemented in a helical fuel rod bundle to investigate the mixing characteristics. Validation is provided by a 5×5 helical fuel bundle mixing experiment. The refined model predicts about 92.7% of the data with the ±10% error range. Compared with existing helical fuel mixing models, the refined mixing model has higher axial accuracy and radial spatial resolution, and can accurately predict the twist angle-dependent crossflow rate and entry effect. Meanwhile, the refined helical fuel mixing model framework is universal and can be effectively used for the mixing prediction of arbitrary geometric helical fuel after the calibration of coefficients.

High-fidelity nuclear data libraries and neutronics simulation tools are essential for the development of fast reactors. The IAEA coordinated research project on “Neutronics Benchmark of CEFR Start-Up Tests” offers valuable data for the qualification of nuclear data libraries and neutronics codes. This paper focuses on the verification and validation of the CEFR start-up modelling using OpenMC Monte-Carlo code against the experimental measurements. The OpenMC simulation results agree well with the measurements in criticality, control rod worth, sodium void reactivity, temperature reactivity, subassembly swap reactivity, and reaction distribution. In feedback coefficient evaluations, an additional state method shows high consistency with lower uncertainty. Among 122 relative errors in the benchmark of the distribution of nuclear reaction, 104 errors are less than 10% and 84 errors are less than 5%. The results demonstrate the high reliability of OpenMC for its application in fast reactor simulations. In the companion paper, the influence of cross-section libraries is investigated using neutronics modelling in this paper.

Nuclear data is the basis of reactor physics analysis. This paper aim at studying the influence of major evaluated nuclear data libraries, CENDL-3.2, ENDF/B-VIII.0, JEFF-3.3, and JENDL-4.0u, on the neutronics modelling of CEFR start-up tests. Results show these four libraries have a good performance and consistency in the modelling CEFR start-up tests. The JEFF-3.3 results exhibit only an 8 pcm keffdifference with the measurement. The difference in criticality is decomposed by nuclide, which shows the large overestimation of CENDL-3.2 is mainly from the cross-section of 52Cr. Except for few cases, the calculation results are within 1σ of measurement uncertainty in control rod worth, sodium void reactivity, temperature reactivity, and subassembly swap reactivity. In the evaluation of axial and radial reaction distribution, there are about 65% of relative errors that are less than 5% and 82% of relative errors that are less than 10%.

The small modular liquid-metal fast reactor (SMFR) is an important component of advanced nuclear systems. SMFRs exhibit relatively low breeding capability and constraint space for control rod installation. Consequently, control rods are deeply inserted at beginning and are withdrawn gradually to compensate for large burnup reactivity loss in a long lifetime. This paper is committed to investigating the impact of control rod compensation operation on core neutronics characteristics. This paper presents a whole core fine depletion model of long lifetime SMFR using OpenMC and the influence of depletion chains is verified. Three control rod position schemes to simulate the compensation process are compared. The results show that the fine simulation of the control rod compensation process impacts significantly the fuel burnup distribution and absorber consumption. A control rod equivalent position scheme proposed in this work is an optimal option in the trade-off between computation time and accuracy. The control position is crucial for accurate power distribution and void feedback coefficients in SMFRs. The results in this paper also show that the pin level power distribution is important due to the heterogeneous distribution in SMFRs. The fuel burnup distribution at the end of core life impacts the worth of control rods.