%0 Journal Article %J SOFTWAREX %D 2018 %T deal2lkit: A toolkit library for high performance programming in deal.II %A Alberto Sartori %A Nicola Giuliani %A Mauro Bardelloni %A Luca Heltai %B SOFTWAREX %V 7 %P 318–327 %G eng %R 10.1016/j.softx.2018.09.004 %0 Journal Article %J Annals of Nuclear Energy, 87, 2 (2016): pp. 198-208 %D 2016 %T A multi-physics reduced order model for the analysis of Lead Fast Reactor single channel %A Alberto Sartori %A Antonio Cammi %A Lelio Luzzi %A Gianluigi Rozza %X In this work, a Reduced Basis method, with basis functions sampled by a Proper Orthogonal Decomposition technique, has been employed to develop a reduced order model of a multi-physics parametrized Lead-cooled Fast Reactor single-channel. Being the first time that a reduced order model is developed in this context, the work focused on a methodological approach and the coupling between the neutronics and the heat transfer, where the thermal feedbacks on neutronics are explicitly taken into account, in time-invariant settings. In order to address the potential of such approach, two different kinds of varying parameters have been considered, namely one related to a geometric quantity (i.e., the inner radius of the fuel pellet) and one related to a physical quantity (i.e., the inlet lead velocity). The capabilities of the presented reduced order model (ROM) have been tested and compared with a high-fidelity finite element model (upon which the ROM has been constructed) on different aspects. In particular, the comparison focused on the system reactivity prediction (with and without thermal feedbacks on neutronics), the neutron flux and temperature field reconstruction, and on the computational time. The outcomes provided by the reduced order model are in good agreement with the high-fidelity finite element ones, and a computational speed-up of at least three orders of magnitude is achieved as well. %B Annals of Nuclear Energy, 87, 2 (2016): pp. 198-208 %I Elsevier %V 87 %P 208 %G en %U http://urania.sissa.it/xmlui/handle/1963/35191 %1 35471 %2 Mathematics %4 1 %# MAT/08 %$ Approved for entry into archive by Maria Pia Calandra (calapia@sissa.it) on 2016-06-16T11:08:05Z (GMT) No. of bitstreams: 0 %& 198 %R doi:10.1016/j.anucene.2015.09.002 %0 Journal Article %J NERS-14-1062; ASME J of Nuclear Rad Sci, 2, 2 (2016) 021019 %D 2016 %T A Reduced Basis Approach for Modeling the Movement of Nuclear Reactor Control Rods %A Alberto Sartori %A Antonio Cammi %A Lelio Luzzi %A Gianluigi Rozza %X This work presents a reduced order model (ROM) aimed at simulating nuclear reactor control rods movement and featuring fast-running prediction of reactivity and neutron flux distribution as well. In particular, the reduced basis (RB) method (built upon a high-fidelity finite element (FE) approximation) has been employed. The neutronics has been modeled according to a parametrized stationary version of the multigroup neutron diffusion equation, which can be formulated as a generalized eigenvalue problem. Within the RB framework, the centroidal Voronoi tessellation is employed as a sampling technique due to the possibility of a hierarchical parameter space exploration, without relying on a “classical” a posteriori error estimation, and saving an important amount of computational time in the offline phase. Here, the proposed ROM is capable of correctly predicting, with respect to the high-fidelity FE approximation, both the reactivity and neutron flux shape. In this way, a computational speedup of at least three orders of magnitude is achieved. If a higher precision is required, the number of employed basis functions (BFs) must be increased. %B NERS-14-1062; ASME J of Nuclear Rad Sci, 2, 2 (2016) 021019 %I ASME %V 2 %P 8 %8 02/2016 %G en %U http://urania.sissa.it/xmlui/handle/1963/35192 %N 2 %1 35473 %2 Mathematics %4 1 %$ Submitted by Gianluigi Rozza (grozza@sissa.it) on 2016-05-30T00:18:10Z No. of bitstreams: 1 sartori_post-icone.pdf: 1377810 bytes, checksum: 5c09a1296518bbbd4b61b709ced8dbca (MD5) %R 10.1115/1.4031945 %0 Journal Article %J Communications in Computational Physics %D 2016 %T Reduced basis approaches in time-dependent noncoercive settings for modelling the movement of nuclear reactor control rods %A Alberto Sartori %A Antonio Cammi %A Lelio Luzzi %A Gianluigi Rozza %X

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.

%B Communications in Computational Physics %I SISSA %8 2016 %G en %U http://urania.sissa.it/xmlui/handle/1963/34963 %N in press %1 35188 %2 Mathematics %$ Submitted by Gianluigi Rozza (grozza@sissa.it) on 2015-11-05T15:06:51Z No. of bitstreams: 1 manuscript.pdf: 2056729 bytes, checksum: 2a2fad22dc931e8c4094642e21c479cd (MD5) %0 Journal Article %D 2015 %T Deal2lkit: a Toolkit Library for High Performance Programming in deal.II %A Alberto Sartori %A Nicola Giuliani %A Mauro Bardelloni %A Luca Heltai %X 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. %I SISSA %G en %U http://urania.sissa.it/xmlui/handle/1963/35006 %1 35235 %2 Mathematics %4 1 %# MAT/08 %$ Submitted by asartori@sissa.it (asartori@sissa.it) on 2015-11-13T13:12:29Z No. of bitstreams: 1 paper_deal2lkit.pdf: 1085486 bytes, checksum: 1dc7e9b9790fe38d2ced7b5328633cd8 (MD5) %0 Journal Article %J Annals of Nuclear Energy %D 2014 %T Comparison of a Modal Method and a Proper Orthogonal Decomposition approach for multi-group time-dependent reactor spatial kinetics %A Alberto Sartori %A Davide Baroli %A Antonio Cammi %A Davide Chiesa %A Lelio Luzzi %A Roberto R. Ponciroli %A Ezio Previtali %A Marco E. Ricotti %A Gianluigi Rozza %A Monica Sisti %X

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.

%B Annals of Nuclear Energy %I Elsevier %V 71 %P 229 %8 09/2014 %G en %U http://urania.sissa.it/xmlui/handle/1963/35039 %1 35270 %2 Physics %4 1 %$ Approved for entry into archive by Maria Pia Calandra (calapia@sissa.it) on 2015-11-18T12:07:02Z (GMT) No. of bitstreams: 0 %& 217 %R 10.1016/j.anucene.2014.03.043 %0 Conference Proceedings %B 22nd International Conference on Nuclear Engineering ICONE22 %D 2014 %T A reduced order model for multi-group time-dependent parametrized reactor spatial kinetics %A Alberto Sartori %A Davide Baroli %A Antonio Cammi %A Lelio Luzzi %A Gianluigi Rozza %X

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.

%B 22nd International Conference on Nuclear Engineering ICONE22 %7 American Society Mechanical Engineering %I American Society of Mechanical Engineers (ASME) %C Prague, Czech Republic %P V005T17A048-V005T17A048 %8 07/2014 %@ 978-079184595-0 %G en %U http://urania.sissa.it/xmlui/handle/1963/35123 %1 35360 %2 Mathematics %4 1 %$ Approved for entry into archive by Maria Pia Calandra (calapia@sissa.it) on 2015-12-03T10:06:02Z (GMT) No. of bitstreams: 0 %R 10.1115/ICONE22-30707