SERPENT validation and optimization with mesh adaptive search on stereolithography geometry models
Not Available
Talamo, Alberto (ORCID:0000000156850483)↗
SEARCH · Engineering Papers
Results for “Serpent”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Serpent and MCNP Calculations of the Energy Deposition in the Transformational Challenge Reactor
This paper focuses on the calculation of the energy deposition in the Transformational Challenge Reactor by two major Monte Carlo codes: Serpent and MCNP. The first software computation relies on Kinetic Energy Released per unit Mass (KERMA) factors while the second one relies on Q-values. The results from these two independent computation methodologies are in very good agreement; however, Serpent runs much faster than MCNP (for the same computational model) and allows for a detailed energy deposition distribution from a 1-mm-side square mesh with a relative statistical error between 0.5% and 1%. This detailed energy deposition is suitable for multiphysics analyses aimed at design optimizations. In order to calculate the energy deposition, Serpent needs enhanced ACE files (distributed by the software developers). Unlike other Monte Carlo software that uses inputs based on Python or Java languages, the Serpent input syntax is very similar to that of MCNP; a Python script can convert a MCNP input to a Serpent input in seconds. For simulations not requiring the calculation of the energy deposition, Serpent can also read nuclear data from MCNP ACE files, which eventually improves the comparison of the results of the two codes.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Verification of the AGREE and Serpent for the steady state HTR-10 benchmark problems
This paper provides results from steady state HTR-10 benchmark calculations performed using AGREE and Serpent that are compared to experimental results as well as calculations performed by INET. The purpose of completing this benchmark is to validate AGREE and Serpent for the prediction of HTGR operation so that they may ultimately be used to support licensing and deployment efforts of advanced reactors. The benchmark consists of several problems ranging from control rod worth calculations to k{sub eff} calculations at various temperatures. Overall, AGREE and Serpent show good agreement with the reference solutions and are effectively able to predict HTGR operation for a variety of steady state cases. The largest difference from the reference was for the initial core single rod worth, which is possibly due to the larger core helium cavity causing inaccuracies in the diffusion calculation whereas the control rod worth determined by Serpent and AGREE is much closer to the reference result for the full core loading. Moreover, in every benchmark problem, increasing the number of energy groups in the cross sections results in improved agreement of the AGREE result with the reference, although it also causes an increase in computation time. For problems with experimental data available, the accuracy of results generated by Serpent and AGREE is comparable to results obtained by other benchmark participants. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
MONTE CARLO ANALYSES OF THE YALINA THERMAL FACILITY WITH SERPENT STEREOLITHOGRAPHY GEOMETRY MODEL
This paper analyzes the YALINA Thermal subcritical assembly of Belarus using two different Monte Carlo transport programs, SERPENT and MCNP. The MCNP model is based on combinatorial geometry and universes hierarchy, while the SERPENT model is based on Stereolithography geometry. The latter consists of unstructured triangulated surfaces defined by the normal and vertices. This geometry format is used by 3D printers and it has been created by: the CUBIT software, MATLAB scripts, and C coding. All the Monte Carlo simulations have been performed using the ENDF/B-VII.0 nuclear data library. Both MCNP and SERPENT share the same geometry specifications, which describe the facility details without using any material homogenization. Three different configurations have been studied with different number of fuel rods. The three fuel configurations use 216, 245, or 280 fuel rods, respectively. The numerical simulations show that the agreement between SERPENT and MCNP results is within few tens of pcms.
CUBIT↗
Criticality calculations of the Very High Temperature Reactor Critical Assembly benchmark with Serpent and SCALE/KENO-VI
Within the framework of the IAEA Coordinated Research Project on HTGR Uncertainty Analysis in Modeling, criticality calculations of the Very High Temperature Critical Assembly experiment were performed as the validation reference to the prismatic MHTGR-350 lattice calculations. Criticality measurements performed at several temperature points at this Japanese graphite-moderated facility were recently included in the International Handbook of Evaluated Reactor Physics Benchmark Experiments and represent one of the few data sets available for the validation of HTGR lattice physics. Here, this work compares VHTRC criticality simulations utilizing the Monte Carlo codes Serpent and SCALE/KENO-VI. Reasonable agreement was found between Serpent and KENO-VI, but only the use of the latest ENDF cross section library release, namely the ENDF/B-VII.1 library, led to an improved match with the measured data. Furthermore, the fourth beta release of SCALE 6.2/KENO-VI showed significant improvements from the current SCALE 6.1.2 version, compared to the experimental values and Serpent.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Verification of the Serpent-Griffin Workflow using the SNAP 8 Experimental Reactor
The Systems for Nuclear Auxiliary Power (SNAP) program accumulated extensive experimental measurements over a span of 15 years. This work builds upon previous studies which validated Serpent against experimental data for various criticality configurations of the SNAP 8 Experimental Reactor (S8ER). A 2-stage sequence is applied here with Serpent used for the generation of few-group cross sections and Griffin as the transport eigenvalue solver. Sensitivity studies are performed to qualify the effect and uncertainty of various parameters, comparing against reference models for the S8ER. The results from this work will then be expanded to create a generalized methodology for the Serpent-Griffin 2-stage approach for microreactor applications.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Serpent transient analyses of GIACINT geometrical change experiments
This paper simulates seven transient experiments with geometrical movements performed at the GIACINT facility of Belarus. These experiments include three slow control rods insertion in more than 25 s, one fast control rod insertion in less than 1 s, and three water moderator draining transients. In the fast control rod insertion experiment, the control rods are first pushed by a mechanical spring and then let fall by gravity in the assembly. In this experiment, the control rods speed varies because of the initial spring force, the gravitational acceleration, and the water buoyancy force. In all other transient experiments, the geometry change occurs with a constant speed. The Serpent computer program was utilized to simulate these transient experiments with geometrical changes by performing two separate computations. One time-independent computation is performed in criticality mode and the other time-dependent computation is performed in dynamic-source mode. The first computation writes the neutron population and delayed neutron precursors for starting the transient simulation. The second computation divides the transient time into several time bins, reads the two files from the first computation, and updates them at the end of each time bin. The Serpent simulation results obtained for this set of experiments are in good agreement with the time dependent experimental measurements. In addition, Serpent and MCNP simulations results show an excellent agreement for the water moderator transient
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Evaluating the X2 initial core zero power physics tests with Serpent-Ants
The validation of the Ants nodal neutronics code for VVER applications is started by modelling the zero power physics tests for the initial core of the Khmelnitsky 2 nuclear power plant as described in the X2 benchmark using the Serpent-Ants two step neutronics calculation chain. The Ants prediction compare favorably against an earlier continuous energy Monte-Carlo reference solution as well as the measured data, except for discrepancies in the SCRAM worth between predicted and measured values. Similar discrepancies have been previously reported for VVER-1000 reactors, and the Ants predicted SCRAM worths match well with the Serpent results. Additionally, the assembly power distribution and axial power distribution predicted by Ants for the critical HZP state of the reactor are also compared to a reference Serpent solution with very good results. Lastly, an investigation into the effects of the few-group structure used in the nodal calculations and the use or lack of leakage correction for group constants on the results shows that the best accuracy is reached with seven or eight energy groups without leakage correction, although a reasonable accuracy can also be obtained with two or three energy groups using fundamental mode leakage correction. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Comparative analysis of energy deposition modes available in Serpent 2 within the framework of the supercritical water reactor - Fuel qualification test reactor physics benchmark
A joint European Canadian Chinese development of a supercritical water-cooled small modular reactor (SCW-SMR) technology is in progress since September 2020 in the framework of a Horizon 2020 project called ECC-SMART. As a main purpose of the project, proper estimates of energy deposition and its spatial distribution are prerequisites for the accurate analysis of safety related parameters of the SCW-SMR concept under development. A supercritical water reactor fuel computational benchmark model, provided by Canadian Nuclear Laboratories, was applied for detailed comparison of different energy deposition calculation options available in the Serpent 2 Monte Carlo code. The effect of energy deposition options on the normalization of the results as well as on the spatial distribution of the energy deposition are discussed. Consistent energy deposition calculation methods are presented between three Monte Carlo codes, viz., Serpent 2, MCNP6 and OpenMC. Although resource-intensive, the use of the coupled neutron-photon transport mode of Serpent 2 is recommended for accurate spatial and quantitative characterization of energy deposition in the SCW-SMR fuel assemblies, accounting for both neutron and photon heating of all the materials. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
IBR-2M Reactor Modeling with MCNP, SERPENT, CUBIT, and GIMP Computer Programs
IBR-2M is a fast research reactor located in Dubna (Russia), which operates in supercritical condition for ~800 μs every 200 ms. Two reflector parts are rotating in opposite directions generating 1.8 GWth peak power when the two parts align with the fuel zone. The reactor core uses plutonium fuel and sodium coolant. The control rods are located in stationary reflector zones. This report describes the computational modeling of the IBR-2M reactor. Two independent computational models were generated using MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry based on macro bodies, e.g. boxes, cylinders, hexagons, spheres, cones, etc., and surfaces. The SERPENT model is based on the STereo Lithographic (STL) geometry representation that is used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The water cooling system of the stationary and the rotating reflectors has a very complicated geometry, involving serpentines and irregular surfaces profiles. The materials in these areas have been homogenized using the GIMP image editor according to the material distributions obtained from the colors of the design drawings. The MCNP and the SERPENT inputs consist of ~1500 and ~1750 lines, respectively. The calculated effective multiplication factors from these two Monte Carlo computer programs with different models agree within 50 pcm.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
HTTR 3 D Cross Section Generation with Serpent and MAMMOTH
The High Temperature Test Reactor (HTTR) is a graphite moderated thermal reactor with a nominal power output of 30 MWth. The purpose of this work is to prepare a set of neutron cross sections for use in the transient analyses of three Loss of Forced Cooling (LOFC) events planned for this reactor. The cross sections are generated with the Serpent 2 3 D Monte Carlo code to include the axial heterogeneity and the strong axial coupling of this core. The approach yields good agreement of results when compared to previously used 2 D lattice models, which require a fine spatial discretization and, surprisingly, comes at a higher computational cost. In addition, the temperature distribution must be considered to obtain reasonable cross sections and the axial leakage component dominates the spectral effects in the core. Consequently, a temperature model as a function of the fuel and moderator temperatures is adopted, and the core is depleted to 390 effective full-power days (EFPDs). A full tabulation of cross sections is prepared at this burnup point for a variety of fuel and moderator temperatures. Next, the multiphysics reactor application MAMMOTH is used to evaluate the quality of these cross sections using a diffusion solver. As expected, the homogenization error in the cross sections is significant and the Super Homogenization (SPH) correction from MAMMOTH is necessary to preserve key reaction rates. The SPH corrected MAMMOTH results are in excellent agreement with the Serpent results and reproduce both the reference power profile and temperature coefficients at each tabulation point. This fact confirms the accuracy of the SPH corrected cross sections and demonstrates the maturity of using Serpent 2 and MAMMOTH for 3 D cross section generation, even for reactors as complex as the HTTR. Finally, the authors recommend that the cross sections should be parametrized with the local burnup and that a full core burnup calculation with coupled thermal fluids would provide a better estimate of the initial condition for any subsequent transient analysis. Unfortunately, this was not within the scope of this work, and direct analysis with these tabulations should be considered as a first order approximation, at best.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Implementation of Weighted Delta-Tracking with Scattering in the Serpent 2 Monte Carlo Code
Serpent 2, a continuous-energy Monte Carlo Reactor physics burnup code developed at the VTT Technical Research Centre of Finland, uses a common rejection sampling technique, known as Woodcock delta-tracking. This method improves the efficiency of Monte Carlo in the complicated geometric layout of nuclear reactors and is widely used despite drawbacks that make it inefficient in specific cases, such as in the presence of strong absorbers. At present, Serpent actively uses delta-tracking, falling back to a surface tracking option in situations where delta-tracking is inefficient to partially offset these weaknesses. Morgan and Kotlyar introduced a Weighted Delta Tracking (WDT) routine to address these issues, replacing the rejection sampling of absorption events with an implicit event. In this research, we examine the theoretical basis for rejection sampling and implicit events, and how they relate to WDT. We then discuss implementation of a routine that extends the WDT routine to include scattering events. The routine is implemented and tested within Serpent. We examine various test cases including a Boiling Water Reactor (BWR) pin cell, and homogenous fuel media based on the content of the fuel of the Transient Test Reactor (TREAT) at Idaho National Lab (INL)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Integration of OpenMC methods into MAMMOTH and Serpent
OpenMC, a Monte Carlo particle transport simulation code focused on neutron criticality calculations, contains several methods we wish to emulate in MAMMOTH and Serpent. First, research coupling OpenMC and the Multiphysics Object-Oriented Simulation Environment (MOOSE) has shown promising results. Second, the utilization of Functional Expansion Tallies (FETs) allows for a more efficient passing of multiphysics data between OpenMC and MOOSE. Both of these capabilities have been preliminarily implemented into Serpent. Results are discussed and future work recommended.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A collision history-based approach to Sensitivity/Perturbation calculations in the continuous energy Monte Carlo code SERPENT
In this work, the implementation of a collision history-based approach to sensitivity/perturbation calculations in the Monte Carlo code SERPENT is discussed. The proposed methods allow the calculation of the eects of nuclear data perturbation on several response functions: the eective multiplication factor, reaction rate ratios and bilinear ratios (e.g., eective kinetics parameters). SERPENT results are compared to ERANOS and TSUNAMI Generalized Perturbation Theory calculations for two fast metallic systems and for a PWR pin-cell benchmark. New methods for the calculation of sensitivities to angular scattering distributions are also presented, which adopts fully continuous (in energy and angle) Monte Carlo estimators.
97 MATHEMATICS AND COMPUTING↗
Reactivity Effects at Mayak Production Association January 2, 1958 Criticality Accident Using Serpent 2 and OpenFOAM
The process nuclear criticality accident that occurred at the Mayak Production Association (Chelyabinsk 40) on January 2, 1958 involving a vessel of uranyl nitrate solution claimed the lives of three workers and left a fourth worker with continuing health problems. There are a myriad of uncertain parameters involved with this accident: What was the molarity of the solution? How much solution was in the vessel at the time of the accident? In what position was the vessel and the solution when it went critical? How important was the impact of reflection due to the workers and/or the floor? These uncertain parameters have made this accident particularly difficult to analyze in the past. This work aims to lower the uncertainty on some of these parameters. A most-probable solution composition is determined by comparing literature on the physical properties of uranyl nitrate solutions to those presented in LA-13638, which describes the accident in question. Using this most-probable solution, the main contributions to the reactivity of the system and hence the eventual accident, are identified through Serpent 2 and OpenFOAM analyses. Serpent 2, a Monte Carlo software tool, is used to perform calculations of the reactivity effects of lowering the vessel toward the floor and the reactivity added by the close proximity of workers. OpenFOAM, a C++ partial differential equation solver toolkit, is used to simulate the fluid inside the vessel as the vessel is tipped. This is done by treating the solution and air inside the vessel as two incompressible, isothermal, and immiscible fluids using a volume of fluid (VoF) approach. The goal of this approach is simply to track the interface between the two fluids, and hence give an accurate description of the geometrical structure of the solution as the vessel is tipped. These two unique tools are then coupled to provide a time-dependent flow simulation to study the effect that the changing geometrical structure had on the criticality of the system, which is novel to the criticality safety field. This work provides a more accurate picture of the accident going forward.
97 MATHEMATICS AND COMPUTING↗
Modeling Continuous Online Refueling with SCALE 6.3.1 and Serpent-2 in the EIRENE Novel Molten Salt Reactor Design
The Molten Salt Reactor (MSR) is a Generation IV advanced fission reactor design in which the coolant, and in some cases the fuel itself, is in the form of liquid molten alkali-halide salts with a fluoride or chloride ionic base. In liquid-fueled MSRs, reactor refueling may be performed online, where fresh fuel salt is added to the core during operation without the need for shutdown periods. Refueling for these reactor designs is typically modeled using either a batch or continuous refueling approach, both of which can be simulated with the SCALE 6.3.1 and Serpent-2 code systems. However, the specific manner in which they are implemented can vary depending upon the desired rate of refueling, whether the refueling rate is constant or variable, and consideration of salt drainage for systems where the in-core salt volume is kept constant. In a novel MSR fuel cycle concept termed the “Sourdough” fuel cycle, fuel salt is allowed to “grow” within the core, with excess salt either being transferred to an external holding tank to maintain a constant core volume or diverted to an upper plenum within the core to allow for volume growth. In this work, continuous refueling was modeled in a thermal-spectrum, LEU-fueled, small MSR design operating with the Sourdough fuel cycle using the SCALE 6.3.1 and Serpent-2 codes. Two different continuous refueling approaches were simulated, with the performance of each being compared to determine which is most suitable for use with the Sourdough fuel cycle concept.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Direct NeTS sampling of nuclear graphite $S(α, β, T)$ in Serpent
For advanced reactor applications, Neural Thermal Scattering (NeTS) modules were developed to predict the thermal scattering law (TSL or $S(α, β, T)$) of a nuclear graphite neutron moderator. NeTS are multi-layer, feedforward artificial neural networks, which act as universal function approximators designed for TSL datasets. In this case, a 4-layer neural network with 164 neurons per layer is trained using FLASSH evaluated data in PyTorch and serialized as a torchscript dictionary to predict $S(α, β, T)$ on-the-fly. Relative, absolute and maximum percent deviations of NeTS from File 7 data generated using the FLASSH code are on the order of 0.01%, 0.1% and 1%, respectively, with low inference latencies of 0.000172 s per $S(α, β, T)$ at a given temperature. Capturing the full dimensionality of possible inelastic neutron-lattice interactions, NeTS functionality is embedded in the Serpent Monte Carlo code, where $S(α, β, T)_{NeTS}$ sampling is conducted on-the-fly and compared to ACE look-up-tables for predicting TREAT criticality. k-eff differences between sampling algorithms of 6 pcm are observed and are within the order of Monte Carlo uncertainty. Compared to discrete and continuous-energy ACE files (30 MB and 131 MB per temperature), the NeTS format is on the order of 200–300 kB for a continuous-temperature, interpolation-free representation of $S(α, β, T)$ and cross sections. NeTS-in-Serpent runtimes comparable with ACE look-up tables are achieved by scaling NeTS for high performance computing architectures with hybrid OpenMP + MPI parallelization. This work validates a novel, self-contained reactor physics framework for predictive cross sections, and demonstrates a general methodology for embedding modern machine learning libraries within existing neutronic analysis frameworks.
Nuclear Criticality Safety Program (NCSP)↗
Activation Analysis for the Inboard Region of FNSF Using SERPENT
This study describes an application of the SERPENT 2 code with the TENDL-2017 nuclear data library and the latest available model features of the Fusion Energy System Studies–Fusion Nuclear Science Facility (FNSF), to evaluate the activation of components after shutdown at 1, 10, and 100 years, assuming a plant lifetime of 8.5 full-power years. The primary parameters evaluated include the specific activity, decay heat, and waste disposal rating (WDR). The specific activity and decay heat are calculated with SERPENT 2 using a 360-deg model of the FNSF, while the WDR is calculated and classified based on the waste disposal limits established by the U.S. Nuclear Regulatory Commission under 10 CFR 61.55 as well as by using the Fetter approach. A python-based script developed for a previous high-level waste classification and analysis study was implemented and adapted to this research to calculate the WDR by comparing nuclide concentrations to the values established in 10 CFR 61.55 to generate a waste classification for each component surveyed. As only three short-lived isotopes have limitations for classifications beyond Class A, of which only 63Ni is present in appreciable quantities, there is a limit to the amount that short-lived isotopes contribute to the most significant waste analyzed here. In most cases, a handful of long-lived isotopes can be problematic, such as 59Ni and 94Nb, for example, which are solely responsible for multiple Class C classifications. The results herein reported heavily depend on the specific materials and mass/volume fractions in the specific model used in this study, which has changed and evolved since the inception of the FNSF concept and past studies. Furthermore, the more significant contributions of this study may be the development of a modeling and simulation toolkit and a strategy to perform these calculations, so to help evaluate and optimize future fusion facilities.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗