In this work, two approaches, based on the certified Reduced Basis method, have been developed for simulating the movement of nuclear reactor control rods, in time-dependent non-coercive settings featuring a 3D geometrical framework. In particular, in a first approach, a piece-wise affine transformation based on subdomains division has been implemented for modelling the movement of one control rod. In the second approach, a “staircase” strategy has been adopted for simulating the movement of all the three rods featured by the nuclear reactor chosen as case study. The neutron kinetics has been modelled according to the so-called multi-group neutron diffusion, which, in the present case, is a set of ten coupled parametrized parabolic equations (two energy groups for the neutron flux, and eight for the precursors). Both the reduced order models, developed according to the two approaches, provided a very good accuracy compared with high-fidelity results, assumed as “truth” solutions. At the same time, the computational speed-up in the Online phase, with respect to the fine “truth” finite element discretization, achievable by both the proposed approaches is at least of three orders of magnitude, allowing a real-time simulation of the rod movement and control.

PB - SISSA UR - http://urania.sissa.it/xmlui/handle/1963/34963 IS - in press U1 - 35188 U2 - Mathematics ER - TY - JOUR T1 - Deal2lkit: a Toolkit Library for High Performance Programming in deal.II Y1 - 2015 A1 - Alberto Sartori A1 - Nicola Giuliani A1 - Mauro Bardelloni A1 - Luca Heltai AB - We present version 1.0.0 of the deal2lkit (deal.II ToolKit) library. deal2lkit is a collection of modules and classes for the general purpose finite element library deal.II. Its principal aim is to provide a high level interface, controlled via parameter files, for those steps that are common in all finite element programs: mesh generation, selection of the finite element type, application of boundary conditions and many others. Each module can be used as a building block independently on the others, and can be integrated in existing finite element codes based on deal.II, drastically reducing the size of programs, rendering their use automatically parametrised, and reducing the overall time-to-market of finite element programming. Moreover, deal2lkit features interfaces with the SUNDIALS (SUite of Nonlinear and DIfferential/ALgebraic equation Solvers) and ASSIMP (Open Asset Import Library) libraries. Some examples are provided which show the aim and scopes of deal2lkit. The deal2lkit library is released under the GNU Lesser General Public License (LGPL) and can be retrieved from the deal2lkit repository https://github.com/mathLab/deal2lkit. PB - SISSA UR - http://urania.sissa.it/xmlui/handle/1963/35006 U1 - 35235 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - Comparison of a Modal Method and a Proper Orthogonal Decomposition approach for multi-group time-dependent reactor spatial kinetics JF - Annals of Nuclear Energy Y1 - 2014 A1 - Alberto Sartori A1 - Davide Baroli A1 - Antonio Cammi A1 - Davide Chiesa A1 - Lelio Luzzi A1 - Roberto R. Ponciroli A1 - Ezio Previtali A1 - Marco E. Ricotti A1 - Gianluigi Rozza A1 - Monica Sisti AB -In this paper, two modelling approaches based on a Modal Method (MM) and on the Proper Orthogonal Decomposition (POD) technique, for developing a control-oriented model of nuclear reactor spatial kinetics, are presented and compared. Both these methods allow developing neutronics description by means of a set of ordinary differential equations. The comparison of the outcomes provided by the two approaches focuses on the capability of evaluating the reactivity and the neutron flux shape in different reactor configurations, with reference to a TRIGA Mark II reactor. The results given by the POD-based approach are higher-fidelity with respect to the reference solution than those computed according to the MM-based approach, in particular when the perturbation concerns a reduced region of the core. If the perturbation is homogeneous throughout the core, the two approaches allow obtaining comparable accuracy results on the quantities of interest. As far as the computational burden is concerned, the POD approach ensures a better efficiency rather than direct Modal Method, thanks to the ability of performing a longer computation in the preprocessing that leads to a faster evaluation during the on-line phase.

PB - Elsevier VL - 71 UR - http://urania.sissa.it/xmlui/handle/1963/35039 U1 - 35270 U2 - Physics U4 - 1 ER - TY - Generic T1 - A reduced order model for multi-group time-dependent parametrized reactor spatial kinetics T2 - 22nd International Conference on Nuclear Engineering ICONE22 Y1 - 2014 A1 - Alberto Sartori A1 - Davide Baroli A1 - Antonio Cammi A1 - Lelio Luzzi A1 - Gianluigi Rozza AB -In this work, a Reduced Order Model (ROM) for multigroup time-dependent parametrized reactor spatial kinetics is presented. The Reduced Basis method (built upon a high-fidelity "truth" finite element approximation) has been applied to model the neutronics behavior of a parametrized system composed by a control rod surrounded by fissile material. The neutron kinetics has been described by means of a parametrized multi-group diffusion equation where the height of the control rod (i.e., how much the rod is inserted) plays the role of the varying parameter. In order to model a continuous movement of the rod, a piecewise affine transformation based on subdomain division has been implemented. The proposed ROM is capable to efficiently reproduce the neutron flux distribution allowing to take into account the spatial effects induced by the movement of the control rod with a computational speed-up of 30000 times, with respect to the "truth" model.

JF - 22nd International Conference on Nuclear Engineering ICONE22 PB - American Society of Mechanical Engineers (ASME) CY - Prague, Czech Republic SN - 978-079184595-0 UR - http://urania.sissa.it/xmlui/handle/1963/35123 N1 - 2014 22nd International Conference on Nuclear Engineering, ICONE 2014; Prague; Czech Republic; 7 July 2014 through 11 July 2014; Code 109131; U1 - 35360 U2 - Mathematics U4 - 1 ER -