The Reactor Antineutrino Anomaly (RAA) refers to the deficit observed between the average event rate measured in reactor antineutrino experiments with respect to the theoretical prediction. This anomaly was first identified in 2011 ($2.5\,σ$) as a consequence of the Huber-Muller reactor antineutrino flux calculation. It was thought to be resolved in 2021 as a result of previous reactor antineutrino flux calculations, with a reduction to about $1\,σ$. In this work, we present the RAA obtained with the latest reactor antineutrino flux calculation published in 2023 by a French research group, which was never used before for the calculation of the RAA. It is the first summation flux model which includes a comprehensive uncertainty budget. The result indicates a revival of the RAA at the level of $2.2\,σ$. We also consider the usual simplest explanation of the RAA by active-sterile neutrino oscillations. We present the constraints on the oscillation parameters and we derive a tension of $3.8\,σ$ with the results of gallium source experiments (Gallium Anomaly) taking into account also the solar neutrino and KATRIN bounds, that of the combined short-baseline reactor spectral ratio measure
Since decades, the CEA has been involved in the development of 4th generation reactors cooled by sodium. The dedicated experimental platform PLATEAU was erected to study the different issues and built a results database for the code validation. As experiments with sodium are complex, the tests are led in mock-ups using water as a simulant fluid. The MICAS mock-up, model of the upper plenum, aims at studying the gas entrainment at the free surface. Gas is prohibited in the reactor core due to neutron effects, which may induce reactivity variations. Since vortices are difficult to handle by numerical codes in complex geometries such as reactor ones, criteria have been developed based on the Rankine vortex model to forecast the air entrainment. However, in some case, numerical results do not correlate the experimental ones and those discrepancies are may be due to the vortex model. This article aims at comparing the vortices characteristics calculated using the Rankine vortex model with experimental results in order to improve the gas entrainment criteria. The parameters of the Rankine model are determined using the velocity fields measured by the Particle Image Velocity (PIV) method
Nuclear reactors are often modeled as inflexible baseload generators with fixed downtimes and restrictive ramping constraints. In practice, however, a reactor's operational flexibility is closely tied to its fuel cycle and associated reactivity margin. A key physical constraint for power maneuverability is xenon poisoning, caused from the transient buildup of neutron-absorbing xenon following a power reduction. This transient can delay or prevent subsequent power ramp-up due to suppressed core reactivity. Additionally, if a reactor is shutdown during periods of low reactivity, restart times can vary significantly, leading to prolonged downtimes. This work introduces a physics-informed modeling framework that embeds fuel cycle dynamics within a unit commitment (UC) formulation. The framework tracks reactivity margin, dynamically enforces xenon induced constraints, and endogenously schedules refueling outages based on core conditions. By capturing intracycle reactivity evolution, the model enables operation dependent nuclear dispatch that reflects both techno-economic requirements and irreducible nuclear physics limits. Application to a representative reactor fleet shows that flexibl
This study presents a cross-disciplinary reactor-to-repository framework to compare different advanced reactors with respect to their spent nuclear fuel (SNF). The framework consists of (1) OpenMC for simulating neutronics, fuel depletion, and radioactive decays; (2) NWPY for computing the repository footprint for SNF disposal given the thermal constraints; and (3) PFLOTRAN for simulating radionuclide transport in the geosphere to compute the peak dose rate, which is used to quantify the repository performance and environmental impact. We first perform the meta-analysis of past comparative analyses to identify the factors led previously to inconsistent conclusions. We then demonstrate the new framework by comparing five reactor types. Significant findings are that (1) the repository footprint is neither linearly related to SNF volume nor to decay heat, due to the repository's thermal constraint, (2) fast reactors have significantly higher I-129 inventory, which is often the primarily dose contributor from repositories, and (3) the repository performance primarily depends on the waste forms. The TRISO-based reactors, in particular, have significantly higher SNF volumes, but result i
A simple and rapid method is proposed for assessing the reduction in the lifetime of steel walls of the reactor vessel under neutron irradiation. The method is based on modeling the number of radiation defects by the behavior of a general time-dependent random process of death and birth and queuing theory. Necessary data for assessments: the estimated operating time of the reactor (in years), the actual operating time of the reactor, the accumulated fluence depending on time, the temperature on the walls of the reactor vessel, the neutron absorption cross section of the steel of the reactor walls, the energy of fast neutrons striking the walls. The main problem: getting this accurate data.
Canada has a rich history of nuclear technology development. Since the 1940s, nuclear research infrastructure and facilities, such as National Research Universal (NRU) reactor at Canadian Nuclear Laboratories in Chalk River, Ontario, have played a key role for R&D and for building Canadian expertise and competency in nuclear technology. The NRU reactor retired in 2018. Since the needs of stakeholders vary with time, consideration of a new multipurpose research reactor for Canada must contemplate current user input for their requirements. As an early step in such consideration, a systematic approach was employed to engage various national stakeholders from academia, industry, and research organizations, including the government. Several virtual workshops were held, each with a specific theme around utilizations. To provide international context, our workshops also included presentations from leading nuclear laboratories around the world. We provide a summary of Canada's approach and the main findings of this early stakeholder engagement.
Neutronic calculations for reactors are a daunting task when using Monte Carlo (MC) methods. As high-performance computing has advanced, the simulation of a reactor is nowadays more readily done, but design and optimization with multiple parameters is still a computational challenge. MC transport simulations, coupled with machine learning techniques, offer promising avenues for enhancing the efficiency and effectiveness of nuclear reactor optimization. This paper introduces a novel benchmark problem within the OpenNeoMC framework designed specifically for reinforcement learning. The benchmark involves optimizing a unit cell of a research reactor with two varying parameters (fuel density and water spacing) to maximize neutron flux while maintaining reactor criticality. The test case features distinct local optima, representing different physical regimes, thus posing a challenge for learning algorithms. Through extensive simulations utilizing evolutionary and neuroevolutionary algorithms, we demonstrate the effectiveness of reinforcement learning in navigating complex optimization landscapes with strict constraints. Furthermore, we propose acceleration techniques within the OpenNeoMC
The trajectories of neutrons in the reactor, the points of their fission of uranium nuclei, the points of neutron absorption, fission chains and chain reactions are considered from the standpoint of fractal geometry and percolation theory. In the study of the stationary critical operating mode of a nuclear reactor, models of Cayley trees and Laplacian fractals are used. This approach allows us to obtain the neutron multiplication equation and an expression for the critical size of the reactor. Models of irreversible growth and various fractal dimensions are also considered as applied to the evolution of neutrons in a reactor. Prospects for the development of the proposed approach to describing reactors, primarily the kinetics and processes of neutron transfer, are indicated.
Antineutrinos from nuclear reactors have the potential to be used for reactor monitoring in the mid- to far-field under certain conditions. Antineutrinos are an unshieldable signal and carry information about the reactor core and the distance they travel. Using gadolinium-doped water Cherenkov detectors for this purpose has been previously proposed alongside rate-only analyses. As antineutrinos carry information about their distance of travel in their energy spectrum, the analyses can be extended to a spectral analysis to gain more knowledge about the detected core. Two complementary analyses are used to evaluate the distance between a proposed gadolinium-doped water-based liquid scintillator detector and a detected nuclear reactor. Example cases are shown for a detector in Boulby Mine, near the Boulby Underground Laboratory in the UK, and six reactor sites in the UK and France. The analyses both show strong potential to range reactors, but are limited by the detector design.
In this paper a study on a fusion reactor core is presented in heavy ion inertial fusion (HIF), including the heavy ion beam (HIB) transport in a fusion reactor, a HIB interaction with a background gas, reactor cavity gas dynamics, the reactor gas backflow to the beam lines, and a HIB fusion reactor design. The HIB has remarkable preferable features to release the fusion energy in inertial fusion: in particle accelerators HIBs are generated with a high driver efficiency of ~30-40%, and the HIB ions deposit their energy inside of materials. Therefore, a requirement for the fusion target energy gain is relatively low, that would be ~50 to operate a HIF fusion reactor with a standard energy output of 1GW of electricity. In a fusion reactor the HIB charge neutralization is needed for a ballistic HIB transport. Multiple mechanical shutters would be installed at each HIB port at the reactor wall to stop the blast waves and the chamber gas backflow, so that the accelerator final elements would be protected from the reactor gas contaminant. The essential fusion reactor components are discussed in this paper.
The Molten Salt Fast Reactor (MSFR) is one of the six GEN-IV reactor designs. In the MSFR, the liquid fuel is the coolant, which moves throughout the primary circuit. This complex phenomenology requires multiphysics modeling. In the present paper, a model of the MSFR is developed in the multiphysics code Cardinal, considering neutronic-thermal hydraulic feedback and the transport of delayed neutron precursors (DNPs) and decay heat precursors (DHPs). OpenMC is used to solve neutronic equations, and NekRS is used to solve mass, momentum, energy, DNPs, and DHPs distribution. A RANS k-t turbulence model is used in NekRS. DNPs and DHPs are modeled using a convective-diffusion equation with modified source terms considering radioactive decay. Cardinal results showed a reasonable behavior for temperature, heat source, velocity, DNPs, and DHPs. However, the current limitations in OpenMC do not allow the modification of delayed neutron source locations. Ongoing efforts look to include this feature in future work to introduce DNP feedback in OpenMC.
The period of a nuclear reactor is usually represented as a deterministic quantity. In this article, the reactor period is described as a random variable. An approach based on the first passage time is used. In this case, the effect of the change in entropy on the value of the period of the reactor is determined. Relations are obtained that describe the effect of entropy on the period of the reactor. Expressions are obtained for the dispersion of the reactor period. This value cannot be found from a deterministic approach. Knowing the effect of total entropy change on the period of the reactor makes it possible to control the behavior of the reactor.
The large quantities of antineutrinos produced through the decay of fission fragments in nuclear reactors provide an opportunity to study the properties of these particles and investigate their use in reactor monitoring. The reactor antineutrino spectra are measured using specialized, large area detectors that detect antineutrinos through inverse beta decay, electron elastic scattering, or coherent elastic neutrino nucleus scattering; although, inverse beta decay is the only demonstrated method so far. Reactor monitoring takes advantage of the differences in the antineutrino yield and spectra resulting from uranium and plutonium fission providing an opportunity to estimate the fissile material composition in the reactor. Recent experiments reveal a deviation between the measured and calculated antineutrino flux and spectra indicating either the existence of yet undiscovered neutrino physics, uncertainties in the reactor source term calculation, incorrect nuclear data, or a combination of all three. To address the nuclear data that impact the antineutrino spectrum calculations and measurements, an international group of over 180 experts in antineutrino physics, reactor analysis, det
Increasing the distance from which an antineutrino detector is capable of monitoring the operation of a registered reactor, or discovering a clandestine reactor, strengthens the Non-Proliferation of Nuclear Weapons Treaty. This report presents calculations of reactor antineutrino interactions, from quasi-elastic neutrino-proton scattering and elastic neutrino-electron scattering, in a water-based detector operated >10 km from a commercial power reactor. It separately calculates signal from the proximal reactor and background from all other registered reactors. The main results are interaction rates and kinetic energy distributions of charged leptons scattered from quasi-elastic and elastic processes. Comparing signal and background distributions evaluates reactor monitoring capability. Scaling the results to detectors of different sizes, target media, and standoff distances is straightforward. Calculations are for two examples of a commercial reactor (P_th~3 GW) operating nearby (L~20 km) an underground facility capable of hosting a detector (~1 kT H2O) project. These reactor-site combinations are Perry-Morton on the southern shore of Lake Erie in the U.S. and Hartlepool-Boulby
Increasing the distance from which an antineutrino detector is capable of monitoring the operation of a registered reactor, or discovering a clandestine reactor, strengthens the Non-Proliferation of Nuclear Weapons Treaty. This paper presents calculations of reactor antineutrino interactions from quasi-elastic neutrino-proton scattering and elastic neutrino-electron scattering in a water-based detector operated $\gtrsim10$ km from a commercial power reactor. It separately calculates signal from the proximal reactor and background from all other registered reactors. The main results are differential and integral interaction rates from the quasi-elastic and elastic processes. There are two underground facilities capable of hosting a detector ($\sim1$ kT H$_2$O) project nearby ($L\sim20$ km) an operating commercial reactor ($P_{th}\sim3$ GW). These reactor-site combinations, which are under consideration for project WATCHMAN, are Perry-Morton on the southern shore of Lake Erie in the United States and Hartlepool-Boulby on the western shore of the North Sea in England. The signal rate from the proximal reactor is about five times greater at the Morton site than at the Boulby site due t
Antineutrinos from a nuclear reactor comprise an unshieldable signal which carries information about the core. A gadolinium-doped, water-based Cherenkov detector could detect reactor antineutrinos for mid- to far-field remote reactor monitoring for non-proliferation applications. Two novel and independent reconstruction and analysis pathways have been developed and applied to a number of representative reactor signals to evaluate the sensitivity of a kiloton-scale, gadolinium-doped Cherenkov detector as a remote monitor. The sensitivity of four detector configurations to nine reactor signal combinations was evaluated for a detector situated in Boulby Mine, close to the Boulby Underground Laboratory in the UK. It was found that a 22~m detector with a gadolinium-doped, water-based liquid scintillator fill is sensitive to a $\rm \sim 3~GW_{th}$ reactor at a standoff of $\rm \sim150~km$ within two years in the current reactor landscape. A larger detector would be required to achieve a more timely detection or to monitor smaller or more distant reactors.
The sensitivity to the sterile neutrino mixing at very short baseline reactor neutrino experiments is investigated. If the reactor core is relatively large as in the case of commercial reactors, then the sensitivity is lost for $Δm^2 \gtrsim$ 1 eV$^2$ due to smearing of the reactor core size. If the reactor core is small as in the case of the experimental fast neutron reactor Joyo, the ILL research reactor or the Osiris reactor, on the other hand, then sensitivity to $\sin^22θ_{14}$ can be as good as 0.03 for $Δm^2 \sim$ several eV$^2$ because of its small size.
Reactor antineutrinos have played a significant role in establishing the standard model of particle physics and the theory of neutrino oscillations. In this article, we review the reactor antineutrino flux and in particular the reactor antineutrino anomaly (RAA) coined over a decade ago. RAA refers to a deficit of the measured antineutrino inverse beta decay rates at very short-baseline reactor experiments compared to the theoretically improved predictions (i.e. the Huber-Mueller model). Since the resolution of several previous experimental anomalies have led to the discovery of non-zero neutrino mass and mixing, many efforts have been invested to study the origin of RAA both experimentally and theoretically. The progress includes the observation of discrepancies in antineutrino energy spectrum between data and the Huber-Mueller model, the re-evaluation of the Huber-Mueller model uncertainties, the potential isotope-dependent rate deficits, and the better agreement between data and new model predictions using the improved summation method. These developments disfavor the hypothesis of a light sterile neutrino as the explanation of RAA and supports the deficiencies of Huber-Mueller
Reactor antineutrinos are used to study neutrino oscillation, search for signatures of non-standard neutrino interactions, and to monitor reactor operation for safeguard applications. The flux and energy spectrum of reactor antineutrinos can be predicted from the decays of the nuclear fission products. A comparison of recent reactor calculations with past measurements at baselines of 10-100m suggests a 5.7% deficit. Precision measurements of reactor antineutrinos at very short baselines O(1-10 m) can be used to probe this anomaly and search for possible oscillations into sterile neutrino species. This paper studies the experimental requirements for a new reactor antineutrino measurement at very short baselines and calculates the sensitivity of various scenarios. We conclude that an experiment at a typical research reactor provides 5σ discovery potential for the favored oscillation parameter space with 3 years of data collection.
KamLAND is a reactor neutrino oscillation experiment with a very long baseline. This experiment successfully measured oscillation phenomena of reactor antineutrinos coming mainly from 53 reactors in Japan. In order to extract the results, it is necessary to accurately calculate time-dependent antineutrino spectra from all the reactors. A simple model of reactor cores and code implementing it were developed for this purpose. This paper describes the model of the reactor cores used in the KamLAND reactor analysis.