Abstract:The Generation IV Sodium-cooled Fast Reactor (SFR) R&D program of the CEA (French Alternative Energies and Atomic Energy Commission), embodied in France by the ASTRID program (600 MWe, 2010–2019), must meet more stringent safety requirements than previous programs. For this reason, the behavior of the reactor in accident phenomenology is studied. We can distinguish three main initiator families that can lead to severe accident conditions in an SFR core: reactivity insertions, local subassembly failures and loss of primary flow. They can evolve in different ways depending on the transient and the reactor characteristics, especially the latter (boiling stabilization, primary power excursion). In order to assess the accidental behavior of innovative SFRs, in particular in case of sodium boiling, dedicated tools have been developed at CEA in recent years. In general, these tools are based on "fine" 3D codes that have highlighted the possibility of stabilized boiling. However, it is difficult to analyze stability with such tools in a dynamic scenario and it is also very complicated to validate complex numerical tools on a dynamic nonlinear phenomenology. In this paper, we propose a methodology that allows to study the stability of SFR nuclear cores under boiling conditions resulting from an unprotected loss-of-coolant accident. This methodology relies on a semi-analytical approach directly applied on a discretized model with possibly thousands of unknowns, rather than on a handful of unknowns derived from a reduced-order model as usually done in the literature when it comes to stability analyses. The nuclear reactor is modeled by a 1D drift thermal-hydraulic model, a 1D neutron diffusion model, and a OD heat transfer model. The semi-analytical methodology developed in this article is based on a continuation method, the Numerical Asymptotic Method, which allows us to compute all the stationary solutions of interest (single-phase liquid, two-phase) with good numerical efficiency (accuracy and computational time). We then perform a linear stability analysis of some of these stationary solutions, which amounts to solving generalized eigenvalue problems with a singular mass matrix. This allows us to study the effect of neutron coupling on the stability of steady-state solutions in the framework of Dynamical Systems Theory through bifurcation diagrams. We can see that neutron coupling extends the stability zone in the studied situation.

Numerical simulation of two-dimensional fractional neutron diffusion model describing dynamical behaviour of sodium-cooled fast reactor

Numerical approximation of a two-dimensional fractional neutron diffusion model describing dynamics of neutron flux in a nuclear reactor

Fractional Bernstein polynomial method for solving time-fractional neutron diffusion systems

An efficient computational technique for solving a time-fractional reaction-subdiffusion model in 2D space

Neutron Generation Based Method for Monte Carlo Three-Dimensional Reactor Time-Dependent Simulation

Research of 3-D Hexagonal Nodal Transport Method for Fast Reactor

Subdiffusive processes in BWRs

Three Dimensional Numerical Simulation of Sodium Flowon Small Water Leak in Steam Generator of Fast Reactor

Stability and bifurcation analysis of a Sodium-cooled Fast Reactor under two-phase flow using 1D multiphysics model

Steady and transient solutions of neutron diffusion equations via computational methods based on Hartley series and higher order finite difference schemes

Numerical simulation of three dimensional unsteady process of hydrogen diffusion in nuclear reactor containment

Three dimensional numerical simulation of sodium flow on small water leak in steam generator of fast reactor

Numerical simulation of non-steady state neutron kinetics of the TRIGA Mark II reactor Vienna

Neutron Point Kinetics Model With A Distributed-Order Fractional Derivative

Numerical simulation of fractional-order reaction–diffusion equations with the Riesz and Caputo derivatives

A novel simulation methodology of fractional order nuclear science model

A high order numerical scheme for solving a class of non‐homogeneous time‐fractional reaction diffusion equation

Research on the High-Performance Computing Method for the Neutron Diffusion Spatiotemporal Kinetics Equation Based on the Convolutional Neural Network

An iterative scheme to include turbulent diffusion in advective-dominated transport of delayed neutron precursors

Comprehensive Numerical Analysis of Time-Fractional Reaction–Diffusion Models with Applications to Chemical and Biological Phenomena

A Consistent One-Dimensional Multigroup Diffusion Model for Molten Salt Reactor Neutronics Calculations