Monte Carlo Code Serpent Research Articles

Insights into calculating Reference Discontinuity Factors with Serpent Monte Carlo code

This study explores the calculation of Reference Discontinuity Factors (RDFs) using the Serpent Monte Carlo code, focusing on the methodology and potential pitfalls. In two-step reactor analyses, consistently generated RDFs are crucial for aligning homogeneous nodal diffusion results with the reference heterogeneous transport solution. However, the Serpent internal diffusion solver, based on the Analytic Function Expansion Nodal (AFEN) method, may not be compatible with other nodal methods such as the Nodal Expansion Method (NEM). Additionally, the solver can suffer from instabilities, particularly in multi-group calculations, leading to erroneous RDFs. Despite these challenges, Serpent can generate the necessary raw data for RDF calculation, which can be accurately processed using external diffusion solvers. Two numerical examples − a 1D fuel-reflector model and a 2D SMR core model − illustrate the effects of consistent and inconsistent RDFs on simulation accuracy. The study emphasizes the importance of using compatible diffusion solvers and thoroughly assessing RDFs to avoid errors in reactor simulations.

Benchmark Analysis of a Prismatic Type-high Temperature Gas-cooled Reactor with the ENDF/B-VII.0, ENDF/B-VII.1 and ENDF/B-VIII.0 Nuclear Data Libraries

In this work we performed a benchmark analysis of the High Temperature Engineering Test Reactor (HTTR) fully-loaded start-up critical core with a 30-bundle loading configuration using the Monte Carlo code Serpent 2 with the recently released nuclear data libraries ENDF/B-VII.0, ENDF/B-VII.1 and ENDF/B-VIII.0. The purpose of the work is to reveal the impacts of using different ENDF nuclear data libraries on neutronics calculations of a prismatic high temperature gas-cooled reactor (HTGR). The benchmark results obtained with Serpent 2 were compared against the available experimental values and those attained by different computer codes (MCNP5 and SCALE) using the nuclear data library ENDF/B-VII.0. The comparative results showed good agreement between Serpent 2, the experimental values, MCNP5 and SCALE. The results also exhibited a notable difference between using ENDF/B-VII.0, ENDF/B-VII.1 and ENDF/B-VIII.0 in predicting the neutronics parameters of the HTTR that suggests further investigation in future work.

Neutronics and Thermo-Fluids Simulation of Generic Pebble-Bed Fluoride-Salt-Cooled High-Temperature Reactor

The fluoride-salt-cooled high-temperature reactor (FHR) is one of the advanced reactors that has been attracting considerable interest from both the research community and the nuclear industry. To help facilitate the nuclear community’s familiarity with the FHR, Kairos Power has developed a generic FHR (gFHR) benchmark. In the research performed here, this benchmark was used to assess innovative modeling methods that combine stochastic and deterministic computer codes to perform the design and analysis of the gFHR. The Monte Carlo code Serpent 2 was used to generate few-group cross sections that were then used in the neutron diffusion and thermal-fluids code AGREE to perform full-core neutronics and thermal-fluids steady-state and transient core analysis. The Argonne National Laboratory code SAM was then used to model the gFHR system and to simulate the load-follow operation of the gFHR.

An OpenFOAM solver for multiphysics modeling of fusion reactor design: The nemoFoam code

This paper introduces a multiphysics code, named nemoFoam, designed within the OpenFOAM environment to address the coupling of neutron and photon transport with thermal-hydraulics in nuclear reactor simulations. The code, conceived both for fusion and fission applications, employs a modular approach where the neutronic module is currently based on multi-group neutron diffusion equations, including a mono-kinetic treatment for photons. The thermal-hydraulic module is built on the standard OpenFOAM solver. The paper focuses on presenting the key features of nemoFoam and discussing the structure of the code, the nuclear and thermal-hydraulic models, and the coupling strategy between them.In order to assess the performance of the neutronic module, the code is applied to a fusion case study, indeed a benchmark against the Serpent Monte Carlo code for neutron and photon transport is performed, applying nemoFoam to the Affordable, Robust and Compact (ARC) fusion reactor design. Simulation results demonstrate good agreement with Monte Carlo benchmarks, emphasizing the potential of the code for multiphysics simulations. The modular design of nemoFoam allows users to extend the implemented model equations to study additional phenomena, focus only on selected aspects independently and, potentially, to add new solvers in each module.

Neutronic assessment and optimization of ACP-100 reactor core models to achieve unit multiplication and radial power peaking factor

This study focuses on the neutronic and burnup assessment of the PWR-based Chinese ACP-100 SMR core. Four different core models with varying fuel compositions and assembly configurations were analyzed using the Monte Carlo code SERPENT 2. Three of these models shared uniformly enriched fuel assemblies of 3 wt.%, 4 wt.%, and 4.45 wt.% 235U respectively, while the fourth model (the proposed design for the ACP-100 in this study) featured a mix of differently enriched fuel assemblies along with the inclusion of burnable absorbers. The burnup analysis of the cores spanned 1570 EFPDs, evaluating several neutronic parameters (beta effective, reactivity coefficients) as well as fuel cycle characteristics (cycle burnup, discharge burnup and cycle length). Changes in the isotopic concentration of fissile, fertile, and poison materials were observed and the assembly-wise linear power distribution, neutron flux spectra, and radial power peaking factors of the models were determined. In the mixed core, a combination of borated water, IFBAs and control rods were employed to achieve exact criticality (unit keff) and uniform power distribution (unit radial power peaking factor). Insertion of control rods into 41 of the 57 assemblies, in conjunction with IFBAs placed strategically, attained a near unit radial power peaking factor (PPF) of 1.0586. Finally, the combination of IFBAs, sol-bor of 936 ppm concentration, and CRs with rod worth of 414.33 pcm per assembly brought down the keff to 1.000, but increased the PPF to 1.117. While the achievement of uniform power distribution came at the cost of fuel burnup, the proposed mixed core successfully remained critical for over 1.5 years and reached a discharge burnup of 23.37 GWd/t for a two-batch refueling scheme.

Analysis of Fission Matrix Databases with Nonuniform Fuel Temperature Profiles for Static and Transient Calculations

Monte Carlo codes are the most accurate way to solve the neutronics in a reactor core but can be computationally expensive, especially for when feedback effects are considered or for transient calculations. In this paper, we use the fission matrix (FM) method to perform static and transient calculations with point-kinetics equations and a quasi-static model for an adiabatic transient with feedback. This was applied to a three-dimensional model of the Transient Reactor Test Facility (TREAT) experimental reactor using the Monte Carlo code Serpent for reference calculations and to generate fission matrix databases (FMDBs). In previous works, FMDBs were generated with uniform fuel temperature profiles. Here, we analyze the use of FMDBs with nonuniform temperature profiles, for static and transient calculations. For static calculations, comparison between Serpent and the FM method using nonuniform and uniform FMDBs showed maximum differences in multiplication factors of 57.0 and 77.9 pcm, respectively. For the fission source distribution, comparisons showed maximum root-mean-square differences of 1.10% and 4.89% for nonuniform and uniform FMDBs, respectively. Similar results were obtained when using homogenized databases. Therefore, using nonuniform FMDBs produces a better approximation than uniform databases. For transient calculations, comparisons between both database sets showed differences of 0.9% and −1.8% for the peak and final total power.

Neutronic examination of the U–Mo accident tolerant fuel for VVER-1200 reactors

In this study, we investigated the possibility of employing accident tolerant fuel (ATF) in VVER-1200/V491 assembly without gadolinium-containing fuel rods using the Monte Carlo code Serpent 1.1.7 with ENDF/B-VII cross-section library. The analysis involves assembly design with reflective boundary conditions. To compare the neutronic performances, U–5Mo, U-7.5Mo, U–10Mo, and U–15Mo fuels were chosen in addition to ordinary UO2 fuel. The concentration of 135Xe, 149Sm, fissile and fertile isotopes with burnup, reactivity feedback with fuel temperature variation, and βeff values were calculated. The results indicate that the fuel cycle length increases by 54.27% for U–5Mo, 32.6% for U-7.5Mo, and 13.8% for U–10Mo, while it decreases by 16.4% for U–15Mo fuel. Additionally, the effect of 95Mo content in natural Mo was investigated by reducing the 95Mo concentration. According to the results, each proposed fuel's fuel cycle length extended when the depletion ratio of 95Mo increased. Additionally, the calculations for reactivity feedback guarantee safe operating conditions for all U-xMo fuels.

Burnup-Dependent Neutron Spectrum Behaviour of a Pressurised Water Reactor Fuel Assembly

Understanding the behaviour of a neutron spectrum with burnup is important for describing various phenomena associated with reactor operation. The quest to understand the neutron spectrum comes with a lot of questions. One question that is usually asked by students is: Does the neutron spectrum harden or soften with burnup? Most textbooks used by students do not provide a definite answer to this question. This paper seeks to answer this question using a 3D model of a standard 17 × 17 pressurised water reactor fuel assembly. Two cases were studied using the Serpent Monte Carlo code: the first considered the fuel assembly with constant boron concentration (traditionally found in many published papers), and the second considered boron iteration (where the boron concentration was reduced with burnup). Neutron spectra for the two cases at beginning of life and end of life were compared for spectral shifts. In addition, thermal spectral indices were used to assess spectrum hardening or softening with burnup. Spectral shifts to lower energies were observed in the thermal region of the neutron spectrum, whereas the fast region experienced no spectral shift. There was an increase in thermal spectral indices indicating that the spectrum became soft with burnup.

Sensitivity and Uncertainty Calculations in Support of the Noise Analysis of the TAPIRO Reactor

In the framework of an assessment of the application of neutron noise methodologies to the TAPIRO reactor, sensitivity and uncertainty analyses on impacting nuclear data for this experimental reactor have been performed using the Serpent Monte Carlo code. Integral parameters relevant to the noise analysis of TAPIRO have been investigated, namely the effective multiplication factor, the effective delayed neutron fraction and the effective generation time. The nuclear data uncertainty has been assessed using different data libraries (ENDF/B-VIII.0 and JEFF-3.3), different energy grids (ECCO-33 and XMAS- 172) and different techniques (Generalised Perturbation Theory and Unscented Transform). A ranking of the most impacting nuclides and reactions in terms of uncertainties has been performed, pointing out that the impact of nuclear data uncertainties on the responses is non-negligible and must be taken into account also in the noise analysis of the reactor.

Neutronic analysis of a heat-pipe cooled reactor

Heat-pipe cooled reactors belong to the group of nuclear reactors using heat pipes filled with liquid metals (such as sodium or potassium) to cool the core. Due to the passive heat removal system, there is no need to use closed loops with pumps, and the reactor can be operated with reduced operational requirements. Consequently, this system can be used in remote locations without access to an electrical grid, or it can be used for space applications.This paper deals with a neutronic study of the Special Purpose Reactor Design-B concept from the Idaho National Laboratory. Using Monte Carlo code Serpent 2 and ENDF/B VIII.0 library, a fixed-temperature model was created to calculate the safety characteristics of the system. This included reactivity coefficients, power distribution, neutron flux spectrum, and criticality safety. In a simplified depletion calculation, an effect of fuel depletion on safety systems was determined as well as a decay heat.

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). The Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

A Consistent One-Dimensional Multigroup Diffusion Model for Molten Salt Reactor Neutronics Calculations

Molten Salt Reactors (MSRs) have recently gained resurged research and development interest in the advanced reactor community. Several computational tools are being developed to capture the strong neutronics/thermal-hydraulics coupling effect in this special reactor configuration. This paper presents a consistent one-dimensional (1D) multigroup neutron diffusion model for MSR analysis, with the primary aim for fast and accurate calculations for long transients, as well as sensitivity and uncertainty analysis of the reactor. A fictitious radial leakage cross section is introduced in the model to properly account for the radial leakage effects of the reactor. The leakage cross section and other consistent neutronics parameters are generated with the Monte Carlo code Serpent using high-fidelity three-dimensional (3D) models. The accuracy of the 1D consistent model is verified by the reference solution from the Monte Carlo model on the Molten Salt Reactor Experiment (MSRE) configuration. The 1D consistent model successfully reproduced the integrated flux from the 3D model and the reactor multiplication factor keff with the error in the range of 95 to 397 pcm (per cent mille), depending on discretized energy group structures. The developed model is also extended to estimate the reactivity loss due to fuel circulation in MSRE. The estimate of reactivity loss in dynamics analysis is in great agreement with the experimental data. This model functions as the first step in the development of a 1D fully neutronics/thermal-hydraulics coupled model for short- and long-term MSRE transient analysis.

Calculation of radiation fields in the VVER-1000 concrete biological shield using Monte Carlo code Serpent

To calculate radiation fields in the concrete biological shield (CBS) of the VVER-1000 reactor in this work, we have developed and applied the Monte Carlo code Serpent simulation framework based on the variance reduction technique. We have quantified the radial, axial, and azimuthal variation of neutron and gamma-ray fluxes and the absorbed dose rate in the CBS. Using the calculation results, we estimate maximum neutron fluence and maximum absorbed dose in the VVER-1000 CBS over the period of 60 and 80 years of the reactor operation and localize the domains of highest radiation exposure. The obtained results are in good agreement with the available data on the VVER-1000 and other pressurized water reactors. We show that the fluence of neutrons with energy above 0.1 MeV decreases by half at a depth of 4 cm of concrete, and the gamma-ray absorbed dose decreases by 40 % at a depth of 13.5 cm. The outcomes from this research will help to assess the effects of prolonged irradiation of the VVER-1000 CBS, which is required for reliable risk assessment for extended operation of nuclear power plants.

CEFR control rod drop transient simulation using RAST-F code system

This study aimed to verify and validate the transient simulation capability of the hybrid code system RAST-F for fast reactor analysis. For this purpose, control rod (CR) drop experiments involving eight separate CRs and six CR groups in the China Experimental Fast Reactor (CEFR) start-up tests were utilized to simulate the CR drop transient. The RAST-F numerical solution, including the neutron population, time-dependent reactivity, and CR worth, was compared against the measurement values obtained from two out-of-core detectors. Moreover, the time-dependent reactivity and CR worth from RAST-F were verified against the results obtained by the Monte Carlo code Serpent using continuous energy nuclear data. A code-to-code comparison between Serpent and RAST-F showed good agreement in terms of time-dependent reactivity and CR worth. The discrepancy was less than 160 pcm for reactivity and less than 110 pcm for CR worth. RAST-F solution was almost identical to the measurement data in terms of neutron population and reactivity. All the calculated CR worth results agreed with experimental results within two standard deviations of experimental uncertainty for all CRs and CR groups. This work demonstrates that the RAST-F code system can be a potential tool for analyzing time-dependent phenomena in fast reactors.

Performance evaluation of a currently in-use dry storage cask design for spent accident tolerant fuel loading case under normal operation condition

Abstract The aim of this study is to assess the performance of an existing dry storage cask design for accident tolerant fuel loading case and to examine the compliance with the safety limits applied currently for dry storage. In the study, for a dry storage cask design currently in use, criticality calculations, dose rate evaluation and thermal analyses are performed in case of loading with the accident tolerant fuel discharged from a PWR. Firstly, for an accident tolerant fuel selected among the concepts proposed for use in light water reactors, burnup analyses are performed by utilizing the Serpent continuous energy code and spent fuel characteristics are determined. Then, criticality analyzes are carried out by using the Serpent Monte Carlo code for the case of loading the accident tolerant fuel into the selected dry storage cask design. Gamma and neutron dose rates at the outer surface and close distances of the storage cask are determined with the Serpent code. To evaluate the thermal performance of the storage cask, thermal analyzes are performed by using the ANSYS Fluent computational fluid dynamics code. The analysis results are compared with the nuclear safety criteria applied to dry storage casks. Results of the analysis show that the dry storage cask design currently in-use does not exceed the criticality, dose rate and maximum surface temperature limits when loaded with spent accident tolerant fuel.

Effects of sparging gas entrainment on void and temperature coefficients in the Transatomic Power molten salt reactor core

A model of the Transatomic Power Molten Salt Reactor core was simulated using the Monte Carlo code Serpent 2 to determine the effects of sparging gas entrainment on the void and temperature coefficients. The reactor core showed a positive void coefficient of around + 170 pcm per 1% increase in void and a shift in transition temperatures from negative to positive temperature coefficient of around –24 K per 1% increase in void. These observations were similar for the homogeneous model chosen in this work and for a heterogeneous (explicit bubble) model. The observed characteristics were caused by decreased fuel mixture self–shielding as gas entrainment levels increase and its synergistic effect with gas thermal expansion. Due to the important safety implications from the observed characteristics, similar reactor designs should be evaluated for their void and temperature coefficients with respect to gas entrainment levels.

Use of remix spent mixed fuel plutonium in the BN-1200 reactor

The VVER-1000 thermal neutron reactor can operate on mixed uranium-plutonium fuel with a content of reactor-grade plutonium up to 5% with a 100% loaded core. In this case, plutonium burns up to 56% of odd isotopes. The energy potential of such plutonium is very low, and its further use in thermal reactors is impractical. However, such plutonium can be used in fast neutron reactors. The paper presents the results of investigating the possibility for such isotopic plutonium composition to be used in the BN-1200 thermal neutron reactor and its value be increased for the plutonium recycle in the reactor. For this purpose, a precision model of the BN-1200 reactor has been developed using the Serpent Monte Carlo code. The model has been verified against the reference values of the nuclear fuel burnup and breeding ratios. The study has shown that such plutonium can be used in the BN-1200 reactor MOX fuel. Maintaining the operating cycle length requires the plutonium fraction in the MOX fuel to be increased up to 2%. In the BN-1200 reactor, the isotopic composition has been found to improve for the further recycle of plutonium in the thermal reactor, i.e. odd plutonium isotopes increase. The fewer odd plutonium isotopes at the beginning of the BN-1200 operating cycle, the greater their increase. It can be seen as the result of the calculation that plutonium from VVER-1000 spent mixed fuel must be loaded into the BN-1200 reactor at least twice to increase the fraction of odd isotopes to the level of reactor-grade plutonium.

Neutron-Physical and Radiation Characteristics of Different Low Enrichment Fuels (Thorium and Uranium) in a Proposed JRTR Research Reactor

Abstract: The Jordan Research and Training Reactor (JRTR) is a 5-MWth, open-pool type, light-water moderated and cooled reactor with a heavy-water reflector system. The core consists of standard Material Testing Reactor (MTR), the plate type of fuel assemblies with low enriched fuel of 19.75% U-235 enrichment. One of the most important purposes of JRTR is producing several radioisotopes e.g. (99Mo/99mTc, 131I and 192Ir isotopes) by using the neutron activation technology. A new computational model has been developed for different low enrichment fuels (Thorium and Uranium) in the JRTR by using the Serpent Monte Carlo code. The purpose of this paper is to validate Th-U233 fuel by comparison of the calculation results of important neutronics parameters, like keff, flux distribution, kinetics parameters, power peaking factors, heat decay and activity of spent fuel. Keywords: U3Si2, Th3Si2, JRTR, Light-water Reactors (LWRs), Extended cycle lengths.

Validation and application of the Dragon5 lattice code for neutronics and burnup analysis of VVER-1000 pin cell and assembly model

This work deals with the validation and assessment of neutronics and burnup analysis capacity of the deterministic lattice code Dragon5 for VVER-1000 hexagonal pin and assembly model configuration. Various VVER-1000 hexagonal pin-cell benchmark cases including the fresh uranium fuel and mixed oxide fuel of ORNL are simulated and compared with the HELIOS and OpenMC code. The data library differences of certain nuclides will lead to large deviations of the results. Accurate results can be presented by the Dragon5 lattice code when compared with OpenMC code. The OECD LEU and MOX assembly benchmarks are computed to further validate the prediction ability for VVER-1000 assembly, reasonable agreement between Dragon5 solutions and published values of other codes is observed for burnup parameters including infinite multiplication factors, concentrations of important actinides and fission product nuclides. Moreover, the post-irradiation of UO2 sample experiment of the Russian Kalinin-1 VVER-1000 reactor is simulated and compared with laboratory measured values to provide the experiment validation of the Dragon5 lattice code. The predicted concentrations of actinides up to Cm244 are compared with HELIOS and Serpent Monte Carlo code. The results show that the Dragon5 lattice code presents better prediction accuracy for Calculated-to-Experimental ratios (C/E's) of most actinide concentrations than the HELIOS and Serpent code. Therefore, the École Polytechnique de Montréal Dragon5 lattice physics code can be applied to the neutronics and burnup analysis of VVER-1000 hexagonal pin-cell and assembly model.

Improvement of Few-Group Homogenized Cross-Sections for RSG-GAS In-Core Fuel Management Code Batan-FUEL

Abstract This paper presents the generation and the verification of the few-group homogenized cross section for Batan-FUEL code. This code is routinely used for fuel management in the G. A. Siwabessy Multipurpose Reactor (RSG-GAS). The Monte Carlo code Serpent 2 code, in conjunction with the latest nuclear data library ENDF/B-VIII.0, was used. Calculations using the existing newly generated few-group cross section data were carried out for the 88th core. The calculated core parameters such as excess reactivity and control rod worth are compared to the available experimental data. On the other hand, the fuel burnup fraction and radial power peaking factor (PPF) are compared to the results of Serpent 2. It was shown that the new cross section data have more consistency and a better agreement excess reactivity and control rod worth compared to the experimental data. Similarly, the U-235 burnup fraction and radial power peaking factor by the new cross section data are also shown to concur well with Serpent 2. The newly generated few-group cross section is recommended to replace the existing ones for the fuel management of RSG-GAS.