The doctoral dissertations of the former Helsinki University of Technology (TKK) and Aalto University Schools of Technology (CHEM, ELEC, ENG, SCI) published in electronic format are available in the electronic publications archive of Aalto University - Aaltodoc.
Aalto

Development of a New Monte Carlo Reactor Physics Code

Jaakko Leppänen

Dissertation for the degree of Doctor of Science in Technology to be presented with due permission of the Department of Engineering Physics and Mathematics for public examination and debate in Auditorium F1 at Helsinki University of Technology (Espoo, Finland) on the 18th of June, 2007, at 12 o'clock noon.

Dissertation in PDF format (ISBN 978-951-38-7019-5)   [2247 KB]
VTT Publications 640, ISSN 1455-0849

Dissertation is also available in print (ISBN 978-951-38-7018-8)
Copyright © 2007 VTT Technical Research Centre of Finland
VTT Publications 640, ISSN 1235-0621
VTT-PUBS-640
TKK-DISS-2315

Abstract

Monte Carlo neutron transport codes are widely used in various reactor physics applications, traditionally related to criticality safety analyses, radiation shielding problems, detector modelling and validation of deterministic transport codes. The main advantage of the method is the capability to model geometry and interaction physics without major approximations. The disadvantage is that the modelling of complicated systems is very computing-intensive, which restricts the applications to some extent. The importance of Monte Carlo calculation is likely to increase in the future, along with the development in computer capacities and parallel calculation.

An interesting near-future application for the Monte Carlo method is the generation of input parameters for deterministic reactor simulator codes. These codes are used in coupled LWR full-core analyses and typically based on few-group nodal diffusion methods. The input data consists of homogenised few-group constants, presently generated using deterministic lattice transport codes. The task is becoming increasingly challenging, along with the development in nuclear technology. Calculations involving high-burnup fuels, advanced MOX technology and next-generation reactor systems are likely to cause problems in the future, if code development cannot keep up with the applications. A potential solution is the use of Monte Carlo based lattice transport codes, which brings all the advantages of the calculation method.

So far there has been only a handful of studies on group constant generation using the Monte Carlo method, although the interest has clearly increased during the past few years. The homogenisation of reaction cross sections is simple and straightforward, and it can be carried out using any Monte Carlo code. Some of the parameters, however, require the use of special techniques that are usually not available in general-purpose codes. The main problem is the calculation of neutron diffusion coefficients, which have no continuous-energy counterparts in the Monte Carlo calculation.

This study is focused on the development of an entirely new Monte Carlo neutron transport code, specifically intended for reactor physics calculations at the fuel assembly level. The PSG code is developed at VTT Technical Research Centre of Finland and one of the main applications is the generation of homogenised group constants for deterministic reactor simulator codes. The theoretical background on general transport theory, nodal diffusion calculation and the Monte Carlo method are discussed. The basic methodology used in the PSG code is introduced and previous studies related to the topic are briefly reviewed. PSG is validated by comparison to reference results produced by MCNP4C and CASMO-4E in infinite two-dimensional LWR lattice calculations. Group constants generated by PSG are used in ARES reactor simulator calculations and the results compared to reference calculations using CASMO-4E data.

Keywords: reactor physics, Monte Carlo method, neutron transport codes, PSG code, deterministic reactor simulator codes, homogenisation, homogenised group constants, nodal diffusion method, LWR lattice calculations

This publication is copyrighted. You may download, display and print it for Your own personal use. Commercial use is prohibited.

© 2007 Helsinki University of Technology


Last update 2011-05-26