APPLICATION OF RESPONSE MATRIX METHOD TO TRANSIENT SIMULATIONS OF NUCLEAR SYSTEMS

¹ KTH Royal Institute of Technology, Division of Nuclear Engineering AlbaNova University Center, 10691 Stockholm, Sweden
² Department of Mechanical and Nuclear Engineering, Kansas State University 3002 Rathbone Hall, Manhattan, KS 66506, USA

mickus@kth.se
jaroberts@ksu.edu
jandufek@kth.se

Published online: 22 February 2021

Abstract

Until recently, reactor transient problems were exclusively solved by approximate deterministic methods. The increase in available computing power made it feasible to approach the transient analyses with time-dependent Monte Carlo methods. These methods offer the first-principle solution to the space-time evolution of reactor power by explicitly tracking prompt neutrons, precursors of delayed neutrons and delayed neutrons in time and space. Nevertheless, a very significant computing cost is associated with such methods. The general benefits of the Monte Carlo approach may be retained at a reduced computing cost by applying a hybrid stochastic-deterministic computing scheme. Among such schemes are those based on the fission matrix and the response matrix formalisms. These schemes aim at estimating a variant of the Greens function during a Monte Carlo transport calculation, which is later used to formulate a deterministic approach to solving a space-time dependent problem. In this contribution, we provide an overview of the time-dependent response matrix method, which describes a system by a set of response functions. We have recently suggested an approach where the functions are determined during a Monte Carlo criticality calculation and are then used to deterministically solve the space-time behaviour of the system. Here, we compare the time-dependent response matrix solution with the transient fission matrix and the time-dependent Monte Carlo solutions for a control rod movement problem in a mini-core reactor geometry. The response matrix formalism results in a set of loosely connected equations which offers favourable scaling properties compared to the methods based on the fission matrix formalism.