Optimization of the TiO2 nanofluid as a coolant in the VVER-1000 nuclear reactor based on the thermal reactivity feedback coefficients via the genetic algorithm

The purpose of this study is to display the neutronic simulation of nanofluid application to reactor core. The variations of VVER-1000 nuclear reactor primary neutronic parameters are investigated by using different volume fraction of nanofluid as coolant. The effect of using nanofluid as coolant on reactor dynamical parameters which play an important role in the dynamical analysis of the reactor and safety core is calculated. In this paper coolant and fuel temperature reactivity coefficients in a VVER-1000 nuclear reactor with nanofluid as a coolant are calculated by using various volume fractions and different sizes of TiO2 (Titania) nanoparticle. For do this, firstly the equivalent cell of the hexagonal fuel rod and the surrounding coolant nanofluid is simulated. Then the thermal hydraulic calculations are performed at different volume fractions and sizes of the nanoparticle. Then, using WIMS and CITATION codes, the reactor core is simulated and the effect of coolant and fuel temperature changes on the effective multiplication factor is calculated. For doing optimization, an artificial neural network is trained in MATLAB using the observed data. The different sizes and various volume fractions are inputs, fuel and coolant temperature reactivity coefficients are outputs. The optimal size and volume fraction is determined using the neural network by implementing the genetic algorithms. In the optimization, volume fraction of 7% and size 77 nm are optimal values.

Investigation of the safety features of advanced PWR assembly using SiC, Zr, FeCrAl and SS-310 as cladding materials

In this work, SiC (Silicon carbide), FeCrAl (ferritic), SS-310 (stainless steel 310) and Zirconium are simulated by MCNPX (Monte Carlo N‐Particle eXtended) code as cladding materials in advanced PWR (Pressurized Water Reactor) assembly. A number of reactor safety parameters are evaluated for the candidate cladding materials as reactivity, cycle length, radial power distribution of fuel pellet, reactivity coefficients, spectral hardening, peaking factor, thermal neutron fraction and delayed neutron fraction. The neutron economy presented by Zr and SiC models is analyzed through the burnup calculations on the unit cell and assembly levels. The study also provided the geometric conditions of all cladding materials under consideration in terms of the relation between fuel enrichment and cladding thickness from the viewpoint to achieve the same discharge burnup as the Zircaloy cladding. It was found that the SiC model participated in extending the life cycle by 2.23% compared to Zr. The materials other than SiC largely decreased discharge burnup in comparison with Zircaloy. Furthermore, the claddings with lower capture cross-sections (SiC and Zr) exhibit higher relative fission power at the pellet periphery. The simulation also showed that using SiC with a thickness of 571.15 μm and 4.83% U-235 can satisfy the EOL irradiation value as Zr. For reactivity coefficient, the higher absorbing materials (SS-310 and FeCrAl) exhibit more negative FTCs, MTCs and VRCs at the BOL But, at the intermediate stages of burnup Zr and SiC have a strong trend of negative reactivity coefficients. Finally, the delayed neutron fraction of SiC and Zr models is the highest among all the four models.

Evaluation of the Esfr End of Equilibrium Cycle State: Spatial Distributions of Reactivity Coefficients

The Horizon 2020 ESFR-SMART project investigates the behaviour of the commercial-size European Sodium-cooled Fast Reactor (ESFR) throughout its lifetime. This paper reports work focused on the End of Equilibrium Cycle (EOEC) loading of the ESFR, including neutronic analysis, core- and zone-wise reactivity coefficients, and more detailed local mapping of important safety-relevant parameters. Sensitivity and uncertainty analysis on these parameters have also been performed and a detailed investigation into decay heat mapping carried out. Due to the scope of this work the results have been split into three papers. The nominal operating conditions and both zone-wise and local mapping of reactivity coefficients are considered in this paper; the sensitivity and uncertainty analysis are detailed in Margulis et al. [1]; and the decay heat mapping calculations are reported in Jimenez-Carrascosa et al. [2]. The work was performed across four institutions using both continuous-energy Monte Carlo and deterministic reactor physics codes. A good agreement is observed between the methods, verifying the suitability of these codes for simulation of large, complicated reactor configurations; and giving confidence in the results for the most limiting ESFR EOEC core state for safety analysis. The results from this work will serve as basis for the transient calculations planned for the next stage of work on the ESFR, allowing for more in-depth studies to be performed on the multiphysics behaviour of the reactor.

SPX Benchmark Part II: Transient Results

The paper presents a transient simulation phase of the new benchmark on a large sodium fast reactor (SFR). This phase of the benchmark is devoted to the modelling of selected operational transients performed during start-up tests of the French SFR Superphénix. Six operational transients were selected for the analysis. The specifications of a simplified thermal hydraulic model equipped with point kinetics reactivity data and boundary conditions for the selected transients are given in the paper. The developed model contains necessary thermal hydraulic description of the primary system components, assumptions to account for thermal expansion reactivity feedbacks from out-of-core structures, neutron kinetics parameters, power distribution, and reactivity coefficients. The neutronic input parameters were obtained with the help of the Monte Carlo code Serpent during the first phase of the benchmark related to static neutronic characterization of the core. In this study, the solution of the transient benchmark was obtained with three thermal hydraulic system codes, namely TRACE, SIM-SFR, and ATHLET. The numerical results, compared to the available experimental data, exhibit a reasonable mutual agreement. Particular discrepancies between calculations and experiments could not be fully resolved. Therefore, a set of recommendations for achieving an improved agreement was proposed. In general, the proposed transient benchmark can be seen as an effective tool for validation and cross comparisons of system codes applied for safety analyses of SFRs, including approbation and comparison of different modelling features for thermal expansion of the out-of-core structures.

Modeling of Reactivity Effects and Transient Behaviour of Large Sodium Fast Reactor

In the paper the reactivity characteristics of the core of the large sodium fast reactor Superphenix (SPX) were evaluated and compared with available experimental data. The analysis was performed using the TRACE system code modified for the fast reactor applications. The simplified core model was developed aiming to overcome the lack of detailed information on design and realistic core conditions. Point Kinetics neutronic model with all relevant reactivity feedbacks was used to calculate transient power. The paper focuses on challenging issue of modelling of the transient thermal responses of primary system structural elements resulting in reactivity feedbacks specific to such large fast reactor which cannot be neglected. For these effects, the model was equipped with dedicated heat structures to reproduce important feedbacks due to vessel wall, diagrid, strongback, control rod drive lines thermal expansion. Peculiarly, application of the model was considered for a whole range of core conditions from zero power to 100% nominal. The developed core model allowed reproducing satisfactorily the core reactivity balance between zero power at 180?C and full power conditions. Additionally, the reactivity coefficients k, g, h at three power levels were calculated and satisfactory agreement with experimental measurements was also observed. The study demonstrated feasibility of application of relatively simple model with adjusted parameters for analysis of different conditions of very complex system.

SPX Benchmark Part I: Results of Static Neutronics

In the paper, the specification of a new neutronics benchmark for a large Sodium cooled Fast Reactor core and results of modelling by different participants are presented. The neutronics benchmark describes the core of the French sodium cooled reactor Superphénix at its startup configuration, which in particular was used for experimental measurement of reactivity characteristics. The benchmark consists of the detailed heterogeneous core specification for neutronic analysis and results of the reference solution. Different core geometries and thermal conditions from cold “as fabricated” up to full power were considered. The reference Monte Carlo solution of Serpent 2 includes data on multiplication factor, power distribution, axial and radial reaction rates distribution, reactivity coefficients and safety characteristics, control rods worth, kinetic data. The results of modelling with seven other solutions using deterministic and Monte Carlo methods are also presented and compared to the reference solution. The comparisons results demonstrate appropriate agreement of evaluated characteristics. The neutronics results will be used in the second phase of the benchmark for evaluation of transient behaviour of the core.

Neutronic Performance of the VVER-1000 Reactor Using Thorium Fuel with ENDF Library

In this paper, neutronic calculations and the core analysis of the VVER-1000 reactor were performed using MCNP6 code together with both ENDF/B-VII.1 and ENDF/B-VIII libraries. The effect of thorium introduction on the neutronic parameters of the VVER-1000 reactor was discussed. The reference core was initially filled with enriched uranium oxide fuel and then fueled with uranium-thorium fuel. The calculations determine the delayed neutron fraction βeff, the temperature reactivity coefficients, the fuel consumption, and the production of the transuranic elements during reactor operation. βeff and the Doppler coefficient (DC) are found to be in agreement with the design values. It is found that the core loaded with uranium and thorium has lower delayed neutron fraction than the uranium oxide core. The moderator temperature coefficients of the uranium-thorium core are found to be higher than those of the uranium core. Results indicated that thorium has lower production of minor actinides (MAs) and transuranic elements (mainly plutonium isotopes) compared with the relatively large amounts produced from the uranium-based fuel UO2.

Neutronic analysis of fuel pin design for the long-life core in a pressurized water reactor

This work presents the neutronic analysis of fuel design for a long-life core in a pressurized water reactor (PWR). In order to achieve a high burnup, a high enrichment U-235 is traditionally considered without special constraints against proliferation. To counter the excess reactivity, Erbium was selected as a burnable poison due to its good depletion performance. Calculations based on a standard fuel model were carried out for the PWR type core using SRAC code system. A parametric study was performed to quantify the neutronically achievable burnup at a number of enrichment levels and for a numerous geometries covering a wide design space of lattice pitch. The fuel temperature and coolant temperature reactivity coefficients as well as the small and large void reactivity coefficients are also investigated. It was found that it is possible to achieve sufficient criticality up to 100 GWd/tHM burnup without compromising the safety parameters.

Model of Reactivity Accident of the RBMK-1000 of the Chornobyl NPP 4th Power Unit

The article presents the results of modeling of the reactivity accident, which resulted in the destruction of reactor RBMK-1000 of the 4th power unit of the Chornobyl NPP on April 26, 1986. The RBMK-1000 reactivity accident model was developed on the basis of the kinetics of the nuclear reactor, taking into account the change in the reactivity of the reactor. Reactivity changes as a result of both external influence (movement of control rods; change in the reactor inlet coolant temperature (density)) and due to the action of reactivity feedback by the parameters of the reactor core (change in the fuel temperature, coolant temperature, concentration of 135Хе, graphite stack temperature, etc.). A similar approach was applied by the authors of the article for the study of transient processes with the operation of accelerated unit unloading mode on VVER-1000, and the validity of such model is confirmed. The study of the reactivity accident on RBMK-1000 was carried out for various combinations of values of the effectiveness of control rods; reactivity coefficients of the coolant temperature and fuel temperature; changes in the temperature of the coolant at the inlet to the reactor. In most of the studied RBMK-1000 reactor accident scenarios, the critical values of fuel enthalpy, at which the process of fuel destruction begins, are reached first. An important result of the research is the conclusion that it is not necessary to reach supercriticality on instantaneous neutrons, supercriticality on delayed neutrons is also sufficient to initiate fuel destruction.