Gamma noise to non-invasively monitor nuclear research reactors

Autonomous nuclear reactor monitoring is a key aspect of the International Atomic Energy Agency's strategy to ensure nonproliferation treaty compliance. From the rise of small modular reactor technology, decentralized nuclear reactor fleets may strain the capacities of such monitoring and requires new approaches. We demonstrate the superior capabilities of a gamma detection system to monitor the criticality of a zero power nuclear reactor from beyond typical vessel boundaries, offering a powerful alternative to neutron-based systems by providing direct information on fission chain propagation. Using the case example of the research reactor CROCUS, we demonstrate how two bismuth germanate scintillators placed outside the reactor vessel can precisely observe reactor criticality using so called noise methods and provide core status information in seconds. Our results indicate a wide range of applications due to the newly gained geometric flexibility that could find use in fields beyond nuclear safety.

First in-core gamma spectroscopy experiments in a zero power reactor

Gamma rays in nuclear reactors, arising either from nuclear reactions or decay processes, significantly contribute to the heating and dose of the reactor components. Zero power research reactors offer the possibility to measure gamma rays in a purely neutronic environment, allowing for validation experiments of dose estimates, computed spectra, and prompt to delayed gamma ratios. The resulting data can contribute to models, code validation and photo atomic/nuclear data evaluation. To date, most experiments have relied on flux measurements using TLDs, ionization chambers, or spectrometers set in low flux areas. The CROCUS reactor allows for flexible detector placement in and around the core, and has recently been outfitted with gamma detection capabilities to fulfill the need for in-core gamma spectroscopy, as opposed to flux. In this paper we report on the experiments and accompanying simulations of gamma spectrum measurements inside a zero power reactor core, CROCUS. It is a two-zone, uranium-fueled light water moderated facility operated by the Laboratory for Reactor Physics and Systems Behaviour (LRS) at the Swiss Federal Institute of Technology Lausanne (EPFL). Herein we also introduce, in detail, the new LEAF system: A Large Energy-resolving detection Array for Fission gammas. It consists of an array of four detectors – two large ø 127 254 mm Bismuth Germanate (BGO) and two smaller ø 12 50 mm Cerium Bromide (CeBr3) scintillators. We describe the calibration and characterization of LEAF followed by first in-core measurements of gamma ray spectra in a zero power reactor at different sub-critical and critical states, and different locations. The spectra are then compared to code results, namely MCNP6.2 pulse height tallies. We were able to distinguish prompt processes and delayed peaks from decay databases. We present thus experimental data from hitherto inaccessible core regions. We provide the data as validation means for codes that attempt to model these processes for energies up to 10 MeV. We finally draw conclusions and discuss the future uses of LEAF. The results indicate the possibility of isotope tracking and burn-up validation.

PISTIL, a reactivity modulation device to probe the transfer function of the nuclear reactor CROCUS

The present article summarizes the development and testing of a reactivity modulation device developed by the French Atomic and Alternative Energies Commission (CEA). It was installed in the CROCUS reactor of the Swiss Federal Institute of Technology in Lausanne (EPFL). Experimental tests were performed in the framework of a collaboration between CEA and EPFL.

The so-called PISTIL device aims at measuring the nuclear reactor transfer function in the frequency range of interest between 1 mHz and 200 Hz, in order to probe the in-core kinetic behavior of prompt and delayed neutrons. The reactivity modulation is obtained from the rotation of cadmium foils.

The design of the system was driven with the objective of installing PISTIL at the center of the CROCUS reactor. Neutronic simulations with TRIPOLI-4 Monte Carlo code were performed to select the suitable design parameters and meet the safety requirements of the reactor operation.

The total reactivity worth of the device, as estimated by TRIPOLI-4 Monte Carlo calculation, was approximately 0.16 $ and the maximum amplitude of the reactivity modulation was about 0.013 $. In-core reactivity calibration was then performed and were consistent as compared to TRIPOLI-4 estimations.

Hybrid halide perovskite neutron detectors

Interest in fast and easy detection of high-energy radiation (x-, γ-rays and neutrons) is closely related to numerous practical applications ranging from biomedicine and industry to homeland security issues. In this regard, crystals of hybrid halide perovskite have proven to be excellent detectors of x- and γ-rays, offering exceptionally high sensitivities in parallel to the ease of design and handling. Here, we demonstrate that by assembling a methylammonium lead tri-bromide perovskite single crystal (CH₃NH₃PbBr₃ SC) with a Gadolinium (Gd) foil, one can very efficiently detect a flux of thermal neutrons. The neutrons absorbed by the Gd foil turn into γ-rays, which photo-generate charge carriers in the CH₃NH₃PbBr₃ SC. The induced photo-carriers contribute to the electric current, which can easily be measured, providing information on the radiation intensity of thermal neutrons. The dependence on the beam size, bias voltage and the converting distance is investigated. To ensure stable and efficient charge extraction, the perovskite SCs were equipped with carbon electrodes. Furthermore, other types of conversion layers were also tested, including borated polyethylene sheets as well as Gd grains and Gd₂O₃ pellets directly engulfed into the SCs. Monte Carlo N-Particle (MCNP) radiation transport code calculations quantitatively confirmed the detection mechanism herein proposed.

TOWARDS THE VALIDATION OF NOISE EXPERIMENTS IN THE CROCUS REACTOR USING THE TRIPOLI-4 MONTE CARLO CODE IN ANALOG MODE

Intrinsic neutron noise experiments offer a non-invasive manner to measure the prompt decay constant or reactivity of fissile systems. Using the fluctuations in the density of fission chains, one can infer the kinetics parameters via correlation analysis such as the Rossi-alpha method. The models allowing for the interpretation of these measurements typically rely on the assumption of the system behaving according to point kinetics. When dealing with systems where point kinetics fail to predict the true time correlation – such as heterogeneous or large cores – the direct simulation of fission chains using Monte Carlo methods appears as the only reliable candidate to provide reference predictions for the correlation functions. Monte Carlo methods using explicit fission model libraries are thus being examined as tools for prediction in noise analysis. In this work we illustrate the developments and simulation results of the analog transport capabilities of the TRIPOLI-4 Monte Carlo code coupled with the LLNL fission library FREYA, as applied to a set of neutron noise experiments carried out in the CROCUS zero-power reactor with emphasis on the identification of spatial effects. To validate the general capability of the code to predict noise correlations, we examine time distributions of the whole core fission and explicit detection reactions. We present the methodology to achieve a good agreement between experiments and simulations. We reproduced experimental results for relative α, within typical biases, and conclude on the general feasibility of the analog method. We further explore a decoupled core model and analyze it using the noise method. The results indicate an effective method to treat decoupled systems.

HIGHLY LOCALIZED AZIMUTHAL MEASUREMENTS IN THE CROCUS REACTOR TOWARDS THE VALIDATION OF HIGH-FIDELITY NEUTRONICS CODES

Highly localized in-core measurements are necessary for the validation of neutron transport calculations with high spatial resolution. In the present work, a miniature neutron detector developed at EPFL in collaboration with PSI was used to carry out a set of thermal neutrons counting measurements in the zero-power CROCUS reactor core within a spatial range in order of mm. The miniature detector, positioned close to the core reflector, shows a gradient of +(4.29 ± 0.10)% in the count rate profile in the radial direction within 1.3 cm, with higher values pointing towards the core reflector because of the higher share of neutrons in the thermal range. On the contrary, in a control rod guide tube the count rate gradient is -(4.37 ± 0.10)% and it is directed towards the core center. The measured values are compared with the azimuthal trend of the normalized 6Li reaction rate calculated with an iterative three-steps method performed with the Monte Carlo code Serpent 2. These measurements proved the feasibility of resolving spatial effects in the mm-range and they represent a basis for further investigating highly spatially-resolved phenomena in the CROCUS core.

VALIDATION OF AXIAL VOID PROFILE MEASURED BY NEUTRON NOISE TECHNIQUES IN CROCUS

Recently a joint project has been carried out between the Paul Scherrer Institut, the Ecole Polytechnique Federale de Lausanne and swissnuclear, an industrial partner, in order to determine the axial void distribution in a channel installed in the reflector of the zero power research reactor CROCUS, using neutron noise techniques. The main objective of the present paper is to report on the validation of the results against an alternative measurement technique using gamma-ray attenuation and simulations with the TRACE code. For the gamma-ray attenuation experiments, the channel used in CROCUS is installed out of the core in a Plexiglass water tank. The source and detector are fixed and the channel is moved axially to keep the geometry of the source/detector arrangement untouched. This is key to measure the void effect by gamma attenuation due to the low contrast of this technique. The paper compares the experimental results obtained with both techniques, with the outcomes of simulations carried out with the TRACE code. Even though the quantitative void fraction estimations are not consistent, the trends obtained with the simulation and experimental techniques are the same. The discrepancies between the various experimental techniques and the simulation outcomes are related to the heterogeneous distribution of the water-air mixture in the radial sections of the channel.

AN INITIATIVE FOR THE DEVELOPMENT AND APPLICATION OF OPEN-SOURCE MULTI-PHYSICS SIMULATION IN SUPPORT OF R&D AND E&T IN NUCLEAR SCIENCE AND TECHNOLOGY

Modelling and simulation (M&S) have gradually become irreplaceable tools in the field of Nuclear Science and Technology (NS&T), including nuclear energy systems. This is partly due to growing computational resources and advances in computational science and partly to the difficulties to finance, build and license new experimental facilities. However, the utilization of M&S for research and development (R&D) and education and training (E&T) applications is somewhat hampered by limited accessibility to controlled and sensitive nuclear M&S tools as well as by the desires of the developers of these tools to retain their intellectual property (IP). Open-source software and open-access data are growingly perceived as means to accelerate innovation by promoting synergistic collaborative developments while lowering the barriers associated to code distribution, modification, and sharing. Open-source software development is ideal for R&D and E&T purposes because it permits the enhancement of understanding, the use of advanced computational methods and it promotes the cooperation among researchers and scientists, without rigorous constraints on quality assurance or reliance on proprietary data for technology-specific validation. As a fundamental research tool, this helps to mitigate constraints related to dual use of such technology. It is in this context that an initiative is being launched under the aegis of the International Atomic Energy Agency (IAEA) to promote the development and application of open-source multi-physics simulation in support of R&D and E&T in NS&T. This paper presents scope and objectives of this initiative.

DEVELOPMENT OF 3D PIN-BY-PIN CORE SOLVER TORTIN AND COUPLING WITH THERMAL-HYDRAULICS

Currently, safety analyses mostly rely on codes which solve both the neutronics and the thermal-hydraulics with assembly-wise nodes resolution as multiphysics heterogeneous transport solvers are still too time and memory expensive. The pin-by-pin homogenized codes can be seen as a bridge between the heterogeneous codes and the traditional nodal assembly-wise calculations. In this work, the pin-by-pin simplified transport solver Tortin has been coupled with a sub-channel code COBRA-TF. The verification of the 3D solver of Tortin is presented at first, showing very good agreement in terms of axial and radial power profile with the Monte Carlo code SERPENT for a small minicore and with the state-of-the-art nodal code SIMULATE5 for a quarter core without feedback. Then the results of Tortin+COBRA-TF are compared with SIMULATE5 for one assembly problem with feedback. The axial profiles of power and moderator temperature show good agreement, while the fuel temperature differ by up to 40 K. This is caused mainly by different gap and fuel conductance parameters used in COBRA-TF and in SIMULATE5.

Delayed gamma fraction determination in the zero power reactor CROCUS

Gamma rays are an inextricable part of a nuclear reactor’s radiation field, and as such require characterization for dose rate estimations required for the radiation protection of personnel, material choices, and the design of nuclear facilities. Most commonplace radiation transport codes used for shielding calculations only included the prompt neutron induced component of the emitted gamma rays. The relative amount of gamma rays that are emitted from delayed processes – the delayed gamma fraction – amount to a significant contribution, e.g. in a typical zero power reactor at steady state is estimated to be roughly a third. Accurate predictions of gamma fields thus require an estimation of the delayed content in order to meaningfully contribute. As a consequence, recent code developments also include delayed gamma sources and require validation data. The CROCUS zero power research reactor at EPFL is part of the NEA IRPhE and has therefore been characterized for benchmark quality experiments. In order to provide the means for delayed gamma validation, a dedicated experimental campaign was conducted in the CROCUS reactor using its newly developed gamma detection capabilities based on scintillators. In this paper we present the experimental determination of the delayed gamma fraction in CROCUS using in-core neutron and gamma detectors in a benchmark reactor configuration. A consistent and flexibly applicable methodology on how to estimate the delayed gamma fraction in zero power reactors has hitherto not existed – we herein present a general experimental setup and analysis technique that can be applied to other facilities. We found that the build-up time of relevant short lived delayed gamma emitters is likely attributed to the activation of the aluminium cladding of the fuel. Using a CeBr3 scintillator in the control rod position of the CROCUS core, we determined a delayed gamma fraction of (30.6±0.6)%.

High resolution measurements with miniature neutron scintillators in the SUR-100 zero power reactor

Three 1-mm3 miniature fiber-coupled scintillators have been used to perform cm-wise resolution measurements of the thermal neutron flux within experimental channels of the SUR-100 facility, a zero power thermal reactor operated by the Institute of Nuclear Technology and Energy Systems at the University of Stuttgart. The detection system is developed at the École Polytechnique Fédérale de Lausanne in collaboration with the Paul Scherrer Institut. Thermal neutrons count rates were measured along the experimental channels I and II, which cross the reactor at the center and tangentially to the core, respectively. The reactor was modelled with the Monte Carlo neutron transport code Serpent-2.1.31. The comparison of experimental and computed reaction rate distributions showed a good agreement within the core region, with discrepancies within 2σ. An unexpected discrepancy, probably caused by a geometric inconsistency in the computational model of the reactor, was observed in the reflector region of the experimental channel I, where a 20% difference (i.e. 8σ) was found between experimental and simulated results. Significant discrepancies, respectively worth 10σ and 15σ, were noticed at distance, in the lead shielding region, for both experimental channels I and II. In addition, reaction rate gradients across the 2.6 cm and 5.4 cm diameters of both channels were measured. A horizontal reaction rate gradient of (9.09 ± 0.20) % was measured within 2.4 cm across the diameter of the experimental channel II, with a difference from computed results of 2%. The absence of a vertical reaction rate gradient inside the experimental channel I was confirmed by measurements.

Design of a 150-miniature detectors 3D core-mapping system for the CROCUS reactor

The present article provides an overview of the design of a three-dimensional (3D) full-core mapping system for the CROCUS reactor, operated at the École polytechnique fédérale de Lausanne (EPFL), Switzerland. The system is composed of 149 miniature neutron detectors distributed within the core double lattice at three main axial levels. The miniature detector technology is based on the optimization of the well-proven coupling of a miniature ZnS:6LiF(Ag) scintillator to a state-of-the-art silicon photomultiplier (SiPM) via jacketed optical fibers. The challenges in the mechanical design, the detector optimization, the core criticality, and the development of the acquisition electronics are strongly interconnected and their combination is addressed in this article. The 3D full-core mapping system is foreseen to be installed in CROCUS in autumn 2021 and it will pave the way for the investigation of 3D dynamic phenomena in nuclear reactor cores.

CORE SIM+ SIMULATIONS OF COLIBRI FUEL RODS OSCILLATION EXPERIMENTS AND COMPARISON WITH MEASUREMENTS

At EPFL, the CROCUS reactor has been used to carry out experiments with vibrating fuel rods. The paper presents a first attempt to employ the measured data to validate CORE SIM+, a neutron noise solver developed at Chalmers University of Technology. For this purpose, the original experimental data are processed in order to extract the necessary information. In particular, detector recordings are scrutinized and detrended, and used to estimate CPSDs of detector pairs. These values are then compared with the ones derived from the CORE SIM+ simulations of the experiments. The main trend of the experimental data along with the values for some detectors are successfully reproduced by CORE SIM+. Further work is necessary on both the experimental and computational sides in order to improve the validation process.

Kilogram-Scale Crystallogenesis of Halide Perovskites for Gamma-Rays Dose Rate Measurements

Gamma-rays (γ-rays), wherever present, e.g., in medicine, nuclear environment, or homeland security, due to their strong impact on biological matter, should be closely monitored. There is a need for simple, sensitive γ-ray detectors at affordable prices. Here, it is shown that γ-ray detectors based on crystals of methylammonium lead tribromide (MAPbBr3) ideally meet these requirements. Specifically, the γ-rays incident on a MAPbBr3 crystal generates photocarriers with a high mobility-lifetime product, allowing radiation detection by photocurrent measurements at room temperatures. Moreover, the MAPbBr3 crystal-based detectors, equipped with improved carbon electrodes, can operate at low bias (≈1.0 V), hence being suitable for applications in energy-sparse environments, including space. The γ-ray detectors reported herein are exposed to radiation from a 60Co source at dose rates up to 2.3 Gy h-1 under ambient conditions for over 100 h, without any sign of degradation. The excellent radiation tolerance stems from the intrinsic structural plasticity of the organic-inorganic halide perovskites, which can be attributed to a defect-healing process by fast ion migration at the nanoscale level. The sensitivity of the γ-ray detection upon volume is tested for MAPbBr3 crystals reaching up to 1000 cm3 (3.3 kg in weight) grown by a unique crystal growth technique.

TRANSIENT CALCULATIONS OF SPERT III EXPERIMENTS

Cold-startup and hot-standby reactivity accident tests conducted at the SPERT III E-core research reactor are analysed with the coupled neutron-kinetic/thermal-hydraulic code system DYN3D-ATHLET. Homogenised 2-group cross sections for DYN3D are thereby generated with the Monte Carlo neutron transport code Serpent 2 in combination with the ENDF/B-VII.1 cross section library. Results in terms of maximum power, energy release, and reactivity compensation are in good agreement with the experimental values. The time-dependent contributions to the reactivity feedback are investigated for both a cold-startup test and a hot-standby test. These findings prove the suitability of the combined application of the simulation codes to predict the reactor dynamic behaviour in the event of prompt-critical and super-prompt critical transients even for small reactor cores. Furthermore, static core characteristics of the SPERT III E-core reactor at cold-startup condition are analysed with using a static DYN3D model, a detailed Serpent reference model, and a simplified Serpent model consistent with the DYN3D model. The critical control rod position and the excess reactivities of both the control rods and the transient rod obtained with the Serpent reference model are consistent with the experimental values. For the same parameters, the DYN3D model is in good agreement with the Serpent simplified model.

COUPLING OF nTRACER TO COBRA-TF FOR HIGH-FIDELITY ANALYSIS OF VVERs

Paul Scherrer Institut is developing a high-resolution multi-physics core solver for VVER analysis. This work presents the preliminary stages of the development, specifically the coupling of the 3D pin-by-pin neutronic solver nTRACER to the sub-channel thermal-hydraulic code COBRA-TF for single assembly multi-physics steady state calculations. The coupling scheme and the modifications performed in the codes are described in details. The results of the coupled nTRACER/COBRA-TF calculations are compared to the ones of a standalone nTRACER calculation where the feedbacks are provided by a simplified 1D thermal-hydraulic solver. The agreement is very good with fuel temperature differences around 10 K which can be attributed to the different correlations used in the various solvers. The cross-comparison of the two multi-physics computational routes serves as a preliminary verification of the coupling scheme developed between nTRACER and COBRA-TF.

Local and high distance neutron and gamma measurements of fuel rods oscillation experiments

We report in the present article on the successful observation using noise analysis of the lateral oscillation of one fuel rod by ±2.5 mm around nominal at 0.1 Hz frequency, using an mm3 miniature neutron scintillator at the rod level, and a BGO gamma detector seven meters away from the reactor core center. The experiment was conducted as part of the COLIBRI program in the CROCUS reactor, which is dedicated to the investigation of reactor noise induced by fuel vibrations. It consists in experiments on rod lateral displacement (static) and oscillation (dynamic) with different rods’ numbers at various relevant amplitudes and frequencies. Its main motivation is the increased amplitudes in the neutron noise distributions recorded in ex- and in-core detectors that have been observed in recent years in Siemens pre-Konvoi type of PWR reactors. The obtained experimental data are used for the purpose of code validation, especially within the framework of the European project CORTEX on reactor noise applications. During the first phase of COLIBRI, the observation of a spatial dependence of the perturbation noise, also called neutron modulation, was demonstrated. In the second phase of COLIBRI starting 2021, it is planned to use a core mapping array of neutron detectors to record its propagation. It consists in about 150 miniature scintillators coupled to optical fibers and SiPM readouts, to be distributed in the reactor core. As a feasibility test, experiments were performed using a miniature scintillator prototype placed on a fuel rod, and oscillating the instrumented rod or the one directly adjacent to the detector. In addition, it is theoretically possible to measure branching or perturbation reactor noise using gamma radiation. Following recent developments on gamma measurements in CROCUS, the fuel oscillation was simultaneously recorded with a gamma detection array, LEAF. Its large BGO detectors were used by placing them at the maximum distance to the core, i.e. seven meters away with a clear line of sight using an experimental channel through the shielding of the reactor cavity.

DOSIMETRY MODELING AND EXPERIMENTAL VALIDATION FOR THE PETALE PROGRAM IN THE CROCUS REACTOR

The PETALE experimental program in the CROCUS reactor intends to provide integral measurements on reactivity worth and dosimetry measurement to constrain nuclear data relative to stainless steel heavy reflectors. The experimental setup consists in eight successive plates of pure iron, pure nickel, pure chromium, or nuclear-grade stainless steel set at the close periphery of the core. The plates are interleaved with up to nine dosimeters that consist of thin activation foils with different possible materials to be sensitive to different ranges of the neutron spectrum. A precise measurement with a good estimation of the uncertainties and correlations is required, especially when comparing reaction rates, e.g. transmission measurement and/or spectral indices.

The present work focuses on the validation of the dosimetry technics developed in preparation of this experimental program. Different aspects are discussed: monitors, efficiency calibration, self-absorption correction, self-shielding and nuclear data uncertainties. The different sources of uncertainties for the experiment-calculation comparisons are characterized, taking into account all the correlation between the different dosimeters. These correlations are a mandatory element for the aimed Bayesian assimilation in order to avoir overfitting when considering dosimeter providing a similar information.

IMPACT OF VARIOUS SOURCE OF COVARIANCE INFORMATION ON INTEGRAL PARAMETERS UNCERTAINTY DURING DEPLETION CALCULATIONS WITH CASMO-5

This paper describes the effect of input uncertainties on a set of integral parameters (kinf, nuclide compositions) associated with the validation of CASMO-5 against PIE data. The nuclear data under consideration are the cross-sections, fission spectrum and neutron multiplicities and fission yields. Various sources of covariance information are considered, novel ones (ENDFB-VIII.0, JEFF-3.3) as well as more widely distributed ones (JENDL-4.0, ENDF/B-VII.1, Scale 6.1 and Scale 6.2). All possible nuclide reaction pairs (cross sections, fission spectrum and averaged number of neutron per fission) have been perturbed, e.g. all isotopes available in both the respective covariance libraries and the CASMO-5 library. The evolution of the uncertainty estimates with exposure is complemented with sensitivity analysis to determine the main contributors to the uncertainty. The Pearson coefficient defined between the model output and a given input is used in this work to assess the part of the variance in the model output coming from the considered input uncertainty. It is a very promising measure of sensitivity as it is computationally cheap even though it assumes linearity of the output with respect to input perturbations. The evolution of the uncertainty with exposure, both in terms of trends and magnitude are however very different. Sensitivity analysis allows determining why the trends and magnitudes are different.

ANALYSIS OF THE FIRST COLIBRI FUEL RODS OSCILLATION CAMPAIGN IN THE CROCUS REACTOR FOR THE EUROPEAN PROJECT CORTEX

The Horizon2020 European project CORTEX aims at developing an innovative core monitoring technique that allows detecting anomalies in nuclear reactors, such as excessive vibrations of core internals, flow blockage, or coolant inlet perturbations. The technique will be mainly based on using the fluctuations in neutron flux recorded by in-core and ex-core instrumentation, from which the anomalies will be differentiated depending on their type, location and characteristics. The project will result in a deepened understanding of the physical processes involved, allowing utilities to detect operational problems at a very early stage. In this framework, neutron noise computational methods and models are developed. In parallel, mechanical noise experimental campaigns are carried out in two zero-power reactors: AKR-2 and CROCUS. The aim is to produce high quality neutron noise-specific experimental data for the validation of the models. In CROCUS, the COLIBRI experimental program was developed to investigate experimentally the radiation noise induced by fuel rods vibrations. In this way, the 2018 first CORTEX campaign in CROCUS consisted in experiments with a perturbation induced by a fuel rods oscillator. Eighteen fuel rods located at the periphery of the core fuel lattice were oscillated between ±0.5 mm and ±2.0 mm around their central position at a frequency ranging from 0.1 Hz to 2 Hz. Signals from 11 neutron detectors which were set at positions in-core and ex-core in the water reflector, were recorded. The present article documents the results in noise level of the experimental campaign. Neutron noise levels are compared for several oscillation frequencies and amplitudes, and at the various detector locations concluding to the observation of a spatial dependency of the noise in amplitude.

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

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 miss-handling 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 137Cs 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.

The COLIBRI experimental program in the CROCUS reactor: characterization of the fuel rods oscillator

The present article presents the mechanical characterization of the fuel rods oscillator developed for the purposes of the COLIBRI experimental program in CROCUS. COLIBRI aims at investigating the radiation noise related to fuel vibrations. The main motivation is the increased amplitudes in the neutron noise distributions recorded in ex- and in-core detectors that have been observed in recent years in Siemens pre-Konvoi type of pressurized water reactors. Several potential explanations have been put forward, but no definitive conclusions could yet be drawn. Among others, changes in fuel assembly or pin vibration patterns, due to recent modifications of assembly structural designs, were pointed out as a possible cause. Computational dynamic tools are currently developed within the Horizon 2020 European project CORTEX, to help with understanding the additional noise amplitude. The COLIBRI program is used for their validation. An in-core device was designed, tested, and licensed between 2015 and 2019 for fuel rods oscillation in CROCUS, in successive steps from out-of-pile tests with dummy fuel rods to critical in-core tests. The characterization of its mechanical behavior is presented, in air and in water, and as a function of the load, for safety and experimental purposes. The device allows simultaneously oscillating up to 18 fuel rods. The maximum oscillation amplitude is 5 mm, while the maximum allowed frequency is 2 Hz, i.e. in the frequency range in which the induced neutron flux fluctuations are most pronounced in nuclear power plants.

Developing and testing a miniature fiber-coupled scintillator for in-core neutron counting in CROCUS

An advanced neutron detection system for highly localized measurements in nuclear reactor cores was developed and tested in the Laboratory for Reactor Physics and System Behaviour (LRS) at the École polytechnique fédérale de Lausanne (EPFL), Switzerland, in close collaboration with the Detector group of the Laboratory for Particle Physics (LTP) at the Paul Scherrer Institute (PSI), Switzerland. The miniature-size detector is based on the coupling of a ZnS:6LiF scintillator/converter screen of 1 mm2 and 0.2 mm thickness with a 10-m optical fiber, the latter being connected to a silicon photomultiplier (SiPM). In this development version, the output signal is processed via analog read-out electronics. The present work documents the characterization of a detection system prototype in the mixedradiation fields o f t he C ARROUSEL f acility a nd i ts t esting in the CROCUS zero-power reactor operated at LRS. The fibercoupled scintillator shows a linear response with the reactor power increase up to 6.5 W (i.e. around 108 cm-2s-1 total neutron flux), with a s ubsequent l oss o f l inearity d ue t o e lectronic dead time of the analog system. Nevertheless, the detector shows excellent neutron counting capabilities whether compared to other localized detection systems available at LRS, e.g. miniature fission chambers and an sCVD diamond detector.

Neutron noise experiments in the AKR-2 and CROCUS reactors for the European project CORTEX

The present article gives an overview of the first experimental campaigns carried out in the AKR-2 and CROCUS reactors within the framework of the Horizon 2020 European project CORTEX. CORTEX aims at developing innovative core monitoring techniques that allow detecting anomalies in nuclear reactors, e.g. excessive vibrations of core internals. The technique will be mainly based on using the fluctuations in neutron flux, i.e. noise analysis. The project will result in a deepened understanding of the physical processes involved. This will allow utilities to detect operational problems at a very early stage, and to take proper actions before such problems have any adverse effect on plant safety and reliability. The purpose of the experimental campaigns in the AKR-2 and CROCUS reactors is to produce noise-specific experimental data for the validation of the neutron noise computational models developed within this framework. The first campaigns at both facilities consisted in measurements at reference static states, and with the addition of mechanical perturbations. In the AKR-2 reactor, perturbations were induced by two devices: a rotating absorber and a vibrating absorber, both sets in experimental channels close to the core. In CROCUS, the project benefited from the COLIBRI experimental program: 18 periphery fuel rods were oscillated at a maximum of ±2 mm around their central position in the Hz range. The present article documents the experimental setups and measurements for each facility and perturbation type.

Design and Simulation of Gamma Spectrometry Experiments in the CROCUS Reactor

Gamma rays in nuclear reactors, arising either from fission or decay processes, significantly contribute to the heating and dose of the reactor components. Zero power research reactors offer the possibility to measure gamma rays in a purely neutronic environment, allowing for validation experiments of computed spectra, dose estimates, reactor noise and prompt to delayed gamma ratios. This data then contributes to models, code validation and photo atomic nuclear data evaluation. In order to contribute to aforementioned experimental data, gamma detection capabilities are being added to the CROCUS reactor facility. The CROCUS reactor is a two-zone, uranium-fueled light water moderated facility operated by the Laboratory for Reactor Physics and Systems Behaviour (LRS) at the Swiss Federal Institute of Technology Lausanne (EPFL). With a maximum power of 100W, it is a zero power reactor used for teaching and research, most recently for intrinsic and induced neutron noise studies. For future gamma detection applications in the CROCUS reactor, an array of four detectors - two large 5”x10” Bismuth Germanate (BGO) and two smaller Cerium Bromide (CeBr3) scintillators - was acquired. The BGO detectors are to be arbitrarily positioned in the core reflector and out of the vessel for measurements at arbitrary distances. The CeBr3 detectors on the other hand are small enough to be set in the guide tubes of the control rods for in-core measurements. We present a study of the neutron and gamma flux in the core and reflector using the MCNP 6.2 and Serpent 2 Monte Carlo codes for coupled neutron and photon transport criticality calculations. More specifically, we investigate and compare predicted spectra as well as reactivity worth of different envisioned experimental setups. We further predict pulse height spectra as well as doses to the crystals with and without cadmium shielding to estimate allowable reactor powers with respect to detector radiation hardness. The results serve as basis for calibration and aid in the design and regulatory approval of the experiments.

In-core dosimetry for the validation of neutron spectra in the CROCUS reactor

The present article describes the preliminary validation study of simulated in-core and reflector n eutron spectra in preparation of oncoming experimental programs in the zeropower reactor CROCUS at EPFL. For this purpose, a set of activation foils were irradiated at three characteristic positions in the CROCUS reactor, and the subsequent activities were analyzed via γ spectrometry. The experimental setup was then modeled with the Monte Carlo neutron transport code Serpent2 and associated with an analysis tool to include the effect of the reactor power history during experiments.

The comparison of calculated and measured reaction rates (C/E) indicates a general consistency (at 2σ) between calculated and measured spectra. However, offsets of C/E values were observed in (n, γ) reactions, up to 18% for 115In and 8% for 63Cu dosimeters. This could be caused by an unexpected isotopic composition, uncertainties in nuclear data, or the spectrometry analysis.

In addition, a 100-groups spectrum unfolding was performed using the experimentally determined reaction rates and the Serpent2 spectra as the prior knowledge. The unfolded spectra were mainly adjusted in the thermal and fast ranges, while few modifications w ere m ade i n t he e pithermal r egion d ue t o the low contribution of epithermal neutrons in activation processes. Moreover, within energy groups where the capture reactions show resonant behavior, flux depletion (up to 38% as compared to the prior spectra) is observed due to the absence of self-shielding effect in the unfolding process. For this purpose, an unfolding method based on energy groups weighting is developed and tested.

Data assimilation of post-irradiation examination data for fission yields from GEF

Nuclear data, especially fission yields, create uncertainties in the predicted concentrations of fission products in spent fuel which can exceed engineering target accuracies. Herein, we present a new framework that extends data assimilation methods to burnup simulations by using post-irradiation examination experiments. The adjusted fission yields lowered the bias and reduced the uncertainty of the simulations. Our approach adjusts the model parameters of the code GEF. We compare the BFMC and MOCABA approaches to data assimilation, focusing especially on the effects of the non-normality of GEF’s fission yields. In the application that we present, the best data assimilation framework decreased the average bias of the simulations from 26% to 14%. The average relative standard deviation decreased from 21% to 14%. The GEF fission yields after data assimilation agreed better with those in JEFF3.3. For Pu-239 thermal fission, the average relative difference from JEFF3.3 was 16% before data assimilation and after it was 12%. For the standard deviations of the fission yields, GEF’s were 100% larger than JEFF3.3’s before data assimilation and after were only 4% larger. The inconsistency of the integral data had an important effect on MOCABA, as shown with the Marginal Likelihood Optimization method. When the method was not applied, MOCABA’s adjusted fission yields worsened the bias of the simulations by 30%. BFMC showed that it inherently accounted for this inconsistency. Applying Marginal Likelihood Optimization with BFMC gave a 2% lower bias compared to not applying it, but the results were more poorly converged.

Power calibration methodology at the CROCUS reactor

In the present article, we detail the method used to experimentally determine the power of the CROCUS zero-power reactor, and to subsequently calibrate its ex-core monitor fission chambers. Knowledge of the reactor power is a mandatory quantity for a safe operation. Furthermore, most experimental research programs rely on absolute fission rates in design and interpretation – for instance, tally normalization of reaction rate studies in dosimetry, or normalization of power spectral density in neutron noise measurements. The minimization of associated uncertainties is only achieved by an accurate power determination method. The main experiment consists in the irradiation, and therefore, the activation of several axially distributed Au-197 foils in the central axis of the core, which activities are measured with a High-Purity Germanium (HPGe) gamma spectrometer. The effective cross sections are determined by MCNP and Serpent Monte Carlo simulations. We quantify the reaction rate of each gold foil, and derive the corresponding fission rate in the reactor. The variance weighted average over the distributed foils then provides a calibration factor for the count rates measured in the fission chambers during the irradiation. We detail the calibration process with minimization of respective uncertainties arising from each sub-step, from power control after reactivity insertion, to the calibration of the HPGe gamma spectrometer. Biases arising from different nuclear data choices are also discussed.

Uncertainty propagation for the design study of the PETALE experimental programme in the CROCUS reactor

The PETALE experimental programme in the CROCUS reactor intends to provide integral measurements to constrain stainless steel nuclear data. This article presents the tools and the methodology developed to design and optimize the experiments, and its operating principle. Two acceleration techniques have been implemented in the Serpent2 code to perform a Total Monte Carlo uncertainty propagation using variance reduction and correlated sampling technique. Their application to the estimation of the expected reaction rates in dosimeters is also discussed, together with the estimation of the impact of the nuisance parameters of aluminium used in the experiment structures.

Bayesian Monte Carlo assimilation for the PETALE experimental programme using inter-dosimeter correlation

This article presents the methodology developed to generate and use dosimeter covariances and to estimate nuisance parameters for the PETALE experimental programme. In anticipation of the final experimental results, this work investigates the consideration of these experimental correlations in the Bayesian assimilation process on nuclear data. Results show that the assimilation of a given set of dosimeters provides a strong constraint on some of the posterior reaction rate predictions of the other dosimeters. It confirms that, regarding the assimilation process, the different sets of dosimeters are correlated.

Assessment of representativity of the PETALE experiments for validation of Swiss LWRs ex-core dosimetry calculations

The international experimental program PETALE will be carried out at the CROCUS research reactor of EPFL. The program aims at measuring neutron penetration in slabs made of materials composing typical LWR reactor pressure vessel. The measurements will be used for code and nuclear data validation and for the evaluation of the reflecting properties of these materials. In this paper the representativity of the PETALE experiments is assessed with respect to operational LWR reactors dosimetry and activation evaluations using the Paul Scherrer Institute (PSI) in-house tool NUSS. The NUSS tool allows the stochastic sampling of nuclear data using covariance matrices available in modern nuclear data libraries and the subsequent running of a Monte Carlo code with the modified data files. The representativity can then be assessed based on the Pearson correlation coefficients. The ultimate goal of the work is first of all to assess if the planned PETALE measurements could be applicable beyond their primary purpose and serve for extending the PSI validation database for LWR reactor dosimetry evaluations. Secondly, provided that the PETALE measurements are found useful for the task above, the information on the correlations between the PETALE neutron detectors’ responses and the reactor dosimetry quantities of interest, shall be presented and discussed.

Data assimilation of post irradiation examination experiments to adjust fission yields

Nuclear data, especially fission yields, create uncertainties in the predicted concentrations of fission products in spent fuel. Herein, we present a new framework that extends data assimilation methods to burnup simulations by using data from post-irradiation examination experiments. The adjusted fission yields improve the bias and reduce the uncertainty of predicted fission product concentrations in spent fuel. Our approach modifies fission yields by adjusting the model parameters of the code GEF with post-irradiation examination experiments. We used the BFMC data assimilation method to account for the non-normality of GEF's fission yields. In the application that we present, the assimilation decreased the average bias of the predicted fission product concentrations from 26% to 15%. The average relative standard deviation decreased from 21% to 14%. The GEF fission yields after data assimilation agreed better with those in ENDF/B-VIII.O. For Pu-239 thermal fission, the average relative difference from ENDF/B-VIII.O was 16% before data assimilation and 11% after. For the standard deviations of the fission yields, GEF's were, on average, 16% larger than those from ENDF/B-VIII.O before data assimilation and 15% smaller after.

Uncertainty propagation based on correlated sampling technique for nuclear data applications

A correlated sampling technique has been implemented to estimate the impact of cross section modifications on the neutron transport and in Monte Carlo simulations in one single calculation. This implementation has been coupled to a Total Monte Carlo approach which consists in propagating nuclear data uncertainties with random cross section files. The TMC-CS (Total Monte Carlo with Correlated Sampling) approach offers an interesting speed-up of the associated computation time. This methodology is detailed in this paper, together with two application cases to validate and illustrate the gain provided by this technique: the highly enriched uranium/iron metal core reflected by a stainless-steel reflector HMI-001 benchmark, and the PETALE experimental programme in the CROCUS zero-power light water reactor.