Statistical Burnup Distribution Of Moving Pebbles In Htr-Pm Reactor

In the pebble-bed high temperature reactor under construction in China, called HTR-PM, the spherical fuel elements continuously flow downward in the cylindrical core. After the discharge, the burnup of each pebble is checked at the core outlet and, according to the achieved burnup level, the pebble might be disposed or reinserted into the upper section of the core, distributing randomly in the radial direction and defining a variable number of passes necessary to achieve the average maximum burnup of 90 MWd/kgU. Discrete Element Method (DEM) simulations have been carried out to achieve a clear understanding of the movement of 420,000 fuel pebbles in the HTR-PM core. At the same time, neutronic properties have been investigated for a single pebble and for the full core with Serpent 2 Monte Carlo code in order to perform a parametrization of the one-group microscopic cross sections at the core-level. The pebble movement has been coupled with the neutronic behavior of a single pebble in a dedicated burnup script called Moving Pebble Burnup (MPB), developed in Matlab. 3,000 single pebble burnup histories were simulated to obtain sufficient statistics and insight on the burnup process in the HTR-PM. The decrease of the average burnup gained per single pass implies that a miss-handling of recirculated fuel elements is unlikely to lead to exceedance of the maximum allowed burnup of 100 MWd/kgU. Furthermore, the core demonstrates a self-compensation effect of burnup, meaning that it always compensates burnup under- or over-runs in the successive passes. Finally, it is possible to conclude that the fuel cycle of the HTR-PM, as it has been laid out, is well-designed and feasible.

Statistical Burnup Distribution of Moving Pebbles in the High-Temperature Reactor HTR-PM

Jiri Krepel, Andreas Pautz, Horst-Michael Prasser, Fanny Vitullo

In the pebble-bed high-temperature reactor under construction in China, called the HTR-PM, the spherical fuel elements continuously flow downward in the cylindrical core. The burnup of each pebble is checked at the core outlet and, according to the achieved burnup level, the pebble might be disposed or reinserted into the upper section of the core. Upon reinsertion, each pebble is radially distributed in a random manner and, according to its downward path, faces different burnup conditions. Hence, the number of passes necessary to achieve the average discharge burnup of 90 MWd/kgU may vary. Discrete element method (DEM) simulations have been carried out to achieve a clear understanding of the movement of the 420000 fuel pebbles in the HTR-PM core. At the same time, neutronics properties have been investigated for a single pebble and for the full core with the Serpent 2 Monte Carlo code. As a result, one-group microscopic cross sections (XS) have been parametrized at the core level. The pebble movement has been loosely coupled with the depletion of a single pebble in a dedicated burnup script called moving pebble burnup (MPB), developed in MATLAB. 3000 single pebble burnup histories were simulated to obtain sufficient statistics and an insight into the HTR-PM burnup process. The decrease of the average burnup gained per single pass implies that a misshandling of recirculated fuel elements is unlikely to lead to an excess of the maximum allowed burnup of 100 MWd/kgU. The core demonstrates a self-compensation effect of burnup, meaning that it always compensates burnup under- or over-runs in the successive passes. In addition, gamma detection of Cs-137 has been studied as a practical method for monitoring the burnup of the discharged pebbles, turning out to be an applicable measurement technique. Finally, it is possible to conclude that the fuel cycle of the HTR-PM, as it has been laid out, is well designed and feasible.

ASME2020

Development of the control assembly pattern and dynamic analysis of the generation IV large gas-cooled fast reactor (GFR)

Gaëtan Girardin

During the past ten years, different independent factors, such as the rapidly increasing worldwide demand in energy, societal concerns about greenhouse gas emissions, and the high and volatile prices for fossil fuels, have contributed to the renewed interest in nuclear technology. It is in this context that the Generation IV international forum (GIF) launched the initiative, in 2000, to collaborate on the research and development (R&D) efforts needed for the next generation, i.e. Generation IV, of nuclear reactors. These advanced systems will be ideally deployed beyond the year 2030, following the Generation III or III+ nuclear power plants, which are mainly based on light water technology and are currently entering deployment. A particular goal set for Generation IV systems is closure of the nuclear fuel cycle. Thus, apart from improvements in safety, they are expected to offer a better utilization of natural resources, as also a minimization of long-lived radioactive wastes. Among the systems selected by the GIF, the Gas-cooled Fast Reactor (GFR) is a highly innovative system with advanced fuel geometry and materials (fuel pellets of mixed uranium-plutonium carbide within a plate-type, honeycomb structure made of SiC). It is in the context of the large, 2400 MWth reference GFR design that the present doctoral research has been conducted, the principal aim having been to develop and qualify the control assembly (CA) pattern and corresponding CA implementation scheme for this system. The work has been carried out in three successive and complementary phases: (1) validation of the neutronics tools, (2) the CA pattern development and related static analysis, and (3) dynamic core behavior studies for hypothetical CA driven transients. The deterministic code system ERANOS and its associated nuclear data libraries for fast reactors were developed and validated for the previous generation of sodium-cooled reactors. The validation of ERANOS for GFR applications was, therefore, the first task to be realized in the present research. This has entailed a systematic reanalysis of the GFR-relevant, integral data generated at PSI during the GCFR-PROTEUS experimental program of the 1970's. Thus, during the first phase of the thesis, the reference PROTEUS test lattice from these experiments has been analyzed with ERANOS-2.0 and its associated, adjusted nuclear data library ERALIB1, in order to derive a reference computational scheme to be used later for the GFR analysis. Additionally, benchmark calculations were performed with the Monte Carlo code MCNPX, allowing one to both check the deterministic results and to analyze the sensitivity to different modern data libraries. It has been found that, for the main reaction rate ratios, the new analysis of the GCFR-PROTEUS reference lattice generally yields good agreement – within 1σ measurement uncertainty – with experimental values and with the Monte Carlo simulations. As shown by the analysis, the predictions were in somewhat better agreement in the case of the adjusted ERALIB1 library. The applicability of ERANOS-2.0/ERALIB1 as the reference neutronics tool for the GFR analysis could thus be demonstrated. Furthermore, neutronics aspects related to the novel features of the GFR, for which new experimental investigations are needed, were highlighted. In the second phase of the research, the CA pattern was developed for the GFR, based on iterative neutronics and thermal-hydraulics calculations, 2D and 3D neutronics models for the reactor core having first been set up using the reference ERANOS-2.0/ERALIB1 computational scheme. For the thermal-hydraulics analysis, the CEA code COPERNIC was used. This design work was followed by the study of an appropriate CA implementation scheme (number of CAs and corresponding positions within the core). Detailed neutronics studies revealed the existence of large CA interaction effects, so-called shadowing/anti-shadowing effects leading to an amplification/reduction of the CA worth. The interactions between the absorber pins within a CA, and between the CAs themselves, were investigated in detail, with the goal to optimize the CA efficiency, in terms of the absorber fraction and minimization of the associated heterogeneity effects. The proposed CA pattern consists of 54 absorber pins placed in a triangular lattice. Each absorber pin is a stainless-steel tube filled with highly enriched 10B boron carbide pellets. As a result of the detailed investigations, the absorber pin diameter could be chosen such as to minimize the pin-to-pin influence within the assembly. In particular, a central part of the CA was designed without any absorber pins (zone filled with stagnant helium). A final reduction of the heterogeneity effect (difference between homogeneous and heterogeneous treatments) to 13% was achieved through this feature. Of special importance, the neutronics investigations performed for the reference GFR core ("2004-Core"), especially those related to the CA interactions, have directly contributed to a new core design ("2007-Core"), with the height-to-diameter ratio having been increased to 0.6, compared to 0.3 for the reference core. During the third phase, detailed coupled, 3D neutron-kinetics (NK) and 1D thermal-hydraulics (TH) models were developed for the GFR core, the aim being to arrive at an in-depth understanding of the 3D core behavior during CA driven transients, especially from the viewpoint of spatial effects. The coupled models were developed using the PARCS code for the 3D NK and the TRACE code for the 1D TH modeling. Particular attention was paid to have each individual fuel sub-assembly and CA represented, in order to allow the analysis of local deformations of the 3D distributions of power and safety related parameters, such as the coolant, cladding and fuel temperatures. The validation of the coupled full-core models was performed against reference ERANOS-VARIANT calculations. In particular, the CA worth and reactivity feedbacks were benchmarked, the discrepancies being shown to be relatively low (

EPFL2009

Plutonium management in the medium term

Rakesh Chawla

For many years various countries with access to commercial reprocessing services have been routinely recycling plutonium as UO2/PuO 2 mixed oxide (MOX) fuel in light water reactors (LWRs). This LWR MOX recycle strategy is still widely regarded as an interim step leading to the eventual establishment of sustainable fast reactor fuel cycles. The OECD/NEA Working Party on the Physics of Plutonium Fuels and Innovative Fuel Cycles (WPPR) has recently completed a review of the technical options for plutonium management in what it refers to as the "medium terra." For the purpose of the review, the WPPR considers the medium term to cover the period from now up to the point at which fast reactor fuel cycles are established on a commercial scale. The review identified a number of different designs of innovative plutonium fuel assemblies intended to be used in current LWR cores, in LWRs with significantly different moderation properties, as well as in high-temperature gas reactors. The full review report describes these various options and highlights their respective advantages and disadvantages. This paper briefly summarizes the main findings of the review.