BRILLIANT

Grant Agreement: 662167

WP5. Deliverable D5.18

Thesis No. 2 offered by the project

Main author(s): A. Skwarcan-Bidakowski, B. Chmielarz, W. Gudowski (KTH)

Contributing author(s):

Date of issue of this report: 20 November, 2017

Start date of project: 01/07/2015 Duration: 36 Months

Project funded by the European Commission under the Horizon 2020 Euratom Framework

Programme for Nuclear Research &Training Activities (2014-2018)

Dissemination Level

PU Public X

RE Restricted to a group specified by the partners of the BRILLIANT project

CO Confidential, only for partners of the BRILLIANT project

Ref. Ares(2018)3480913 - 30/06/2018

BRILLIANT D-N°: 5.18 – Thesis nr. 2 offered by the Brilliant Project Dissemination level: PU Date of issue of this report : 20/11/2017

2

Master’s Thesis No. 2 offered by the project:

Nuclear reactor core model for the advanced

nuclear fuel cycle simulator FANCSEE. Advanced

use of Monte Carlo methods in nuclear reactor

calculations

Alexander Skwarcan-Bidakowski

Thesis supervisors: Blazej Chmielarz (blazejc@kth.se) and Wacław

Gudowski (waclaw@kth.se)

In co-operation with Uppsala University

Kungliga Tekniska högskolan – Nuclear Reactor Physics, Stockholm,

Sweden

BRILLIANT D-N°: 5.18 – Thesis nr. 2 offered by the Brilliant Project Dissemination level: PU Date of issue of this report : 20/11/2017

3

Table of content:

PART 1. EXECUTIVE SUMMARY .................................................................................................... 4

PART 2. THESIS: Nuclear reactor core model for the advanced nuclear fuel cycle

simulator FANCSEE. Advanced use of Monte Carlo methods in nuclear

reactor calculations ....................................................................................................... 5

BRILLIANT D-N°: 5.18 – Thesis nr. 2 offered by the Brilliant Project Dissemination level: PU Date of issue of this report : 20/11/2017

4

PART 1. EXECUTIVE SUMMARY

A detailed reactor core modelling of the LOVIISA-2 PWR and FORSMARK-3 BWR was performed in the Serpent 2 Continuous Energy Monte-Carlo code in order to incorporate these results to the advanced nuclear fuel cycle simulator – FANCSEE.

Both models of the reactors were completed but the approximations of the atomic densities of nuclides present in the core differed significantly.

In the LOVIISA-2 PWR, the predicted atomic density for the nuclides approximated by Chebyshev Rational Approximation method (CRAM) coincided with the corrected atomic density simulated by the Serpent 2 program. In the case of FORSMARK-3 BWR, the atomic density from CRAM poorly approximated the data returned by the simulation in Serpent 2. Due to boiling of the moderator in the core of FORSMARK-3, the model seemed to encounter problems of fission density, which yielded unusable results.

The results based on the models of the reactor cores are significant to the FANCSEE Nuclear fuel cycle simulator, which will be used as a dataset for the nuclear fuel cycle burnup in the reactors.

BRILLIANT D-N°: 5.18 – Thesis nr. 2 offered by the Brilliant Project Dissemination level: PU Date of issue of this report : 20/11/2017

5

PART 2. THESIS: Nuclear reactor core model for the advanced

nuclear fuel cycle simulator FANCSEE. Advanced use of Monte

Carlo methods in nuclear reactor calculations

TVE-F 17 009 juni

Examensarbete 15 hpJuni 2017

Nuclear reactor core model for the advanced

nuclear fuel cycle simulator FANCSEE. Advanced

use of Monte Carlo methods in nuclear reactor

calculations

Alexander Skwarcan-Bidakowski

Teknisk- naturvetenskaplig fakultet UTH-enheten Besöksadress: Ångströmlaboratoriet Lägerhyddsvägen 1 Hus 4, Plan 0 Postadress: Box 536 751 21 Uppsala Telefon: 018 – 471 30 03 Telefax: 018 – 471 30 00 Hemsida: http://www.teknat.uu.se/student

Abstract

Nuclear reactor core model for the advanced nuclearfuel cycle simulator FANCSEE. Advanced use ofMonte Carlo methods in nuclear reactor calculationsAlexander Skwarcan-Bidakowski

A detailed reactor core modeling of the LOVIISA-2 PWR and FORSMARK-3 BWR was performed in the Serpent 2 Continuous Energy Monte-Carlo code.

Both models of the reactors were completed but the approximations of the atomic densities of nuclides present in the core differed significantly.

In the LOVIISA-2 PWR, the predicted atomic density for the nuclides approximated by Chebyshev Rational Approximation method (CRAM) coincided with the corrected atomic density simulated by the Serpent 2 program. In the case of FORSMARK-3 BWR, the atomic density from CRAM poorly approximated the data returned by the simulation in Serpent 2. Due to boiling of the moderator in the core of FORSMARK-3, the model seemed to encounter problems of fission density, which yielded unusable results.

The results based on the models of the reactor cores are significant to the FANCSEE Nuclear fuel cycle simulator, which will be used as a dataset for the nuclear fuel cycle burnup in the reactors.

ISSN: 1401-5757, TVE-F 17 009 juniExaminator: Martin SjödinÄmnesgranskare: Changqing RuanHandledare: prof. Waclaw Gudowski, Blazej Chmielarz

1 Populärvetenskaplig sammanfattning av projektet

Kärnkraftverk kan vara svaret till en hållbar och pålitlig energiutvinning som ersätter förbrännin-gen av fossila bränslen. Idag står kolkraftverk för ca 40% av världens elektricitet men samtidigthar det varit en av anledningarna till de höga koldioxidutsläppen. Ersättning av kolkraftverk medexempelvis kärnkraft skulle kräva en del planering och uppskattning av bland annat, hur mycketkärnbränsle som krävs och hur mycket avfall som måste förvaras eller anrikas för vidare operation.

Dagen teknik tillåter användningen av dator-simuleringar för en bra uppskattning av dessa värden.FANCSEE är en sådan kärnbränsle-simulator som nuvarande utvecklas i Institutionen för Reak-torfysik på Kungliga Tekniska Högskolan och ska ha ett grafisk användargränssnitt till skillnad frånandra programvaror. Detta gör det enklare för beslutsfattare att göra informativa beslut. Pro-gramvaran använder sig av värden framtagna från simulationer av detaljerade modeller baseradepå verkliga reaktorer.

I detta projekt har två brett använda reaktormodeller undersökts och modelerats i ett programsom simulerar reaktionerna i en reaktorkärna. Resultaten från dessa simuleringar används sedan iFANSCEE för att uppskatta kärnavfall, bränslebehov, energinätkapacitet eller budget för en kärn-reaktor. Har man kunskap av hur mycket kärnavfall som produceras, kan man exempelvis i godtid planera hur mycket som ska förvaras eller berikas för fortsatt användning.

Just nu är FANCSEE fortfarande under utveckling. Fler reaktormodeller ska implementeras förbredare data-bibliotek. Med detta verktyg kommer kärnbränsleplaneringen gå fortare, och varaenklare än någonsin tidigare.

1

Contents

1 Populärvetenskaplig sammanfattning av projektet 1

2 Introduction 32.1 Principles of nuclear reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32.2 Objective . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32.3 FANCSEE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

3 Monte Carlo method 4

4 Serpent 2 5

5 Theory 65.1 Process of radioactive decay and transmutation through neutron induced reactions 75.2 Bateman Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85.3 Neutron cross-section . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95.4 Burnup calculation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105.5 Chebyshev Rational Approximation method . . . . . . . . . . . . . . . . . . . . . . 11

6 Method 126.1 LOVIISA-2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136.2 FORSMARK-3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 176.3 Implementation of core design to Serpent 2 syntax . . . . . . . . . . . . . . . . . . 216.4 Reactor profile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

7 Results 247.1 LOVIISA-2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 247.2 FORSMARK-3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26

8 Conclusions 28

9 Future work 28

10 Acknowledgements 29

2

2 Introduction

Nuclear energy provides about 11% of the electricity of the world, while burning of coal provides astaggering 40% [1]. The need for nuclear energy will likely increase as the burning of coal decreasesdue to the global problem of climate change. Though renewable energy sources are popular withthe public, they are not as reliable or stable as nuclear power, hence nuclear power is likely to carrythe burden of increased energy demand. For a number of reasons, including safety, politics andenvironmental concerns, nuclear power is still a controversial topic to the public, but only throughcontinuous research, development and communication these concerns can be addressed. Workingtowards sustainable, safe and economically viable nuclear power, simulations are crucial to findingscenarios of reactor cycles. These reactor cycles include, but are not limited to, simulations ofneutron reactions, core criticality or the burnup of nuclear fuel and nuclear transmutation, thelatter two whom will be the objective of this paper.

2.1 Principles of nuclear reactors

The energy generated by nuclear reactors come from the fission of certain heavy metals, ie. elementswith proton number between 90 and 100. Nuclear reactors mostly use uranium as nuclear fuel,which is extracted from natural uranium ore. Natural uranium is mostly composed of uranium-238with low enrichments of the fissile isotope uranium-235 (mass fraction 0.72%), which can be foundnaturally in the environment. The fuel however, is generally an enriched mixture of 238U and 235U .The measure of how much energy is extracted from the nuclear fuel is called Burnup, also knownas fuel utilization and is measured in fissions per initial metal atom, %FIMA or the total energyreleased per ton of initial mass of uranium [MWd/kgU ].

2.2 Objective

The objective of this paper is to simulate neutronics of two detailed nuclear reactor core designs inthe Monte-Carlo particle transport code [2] - Serpent 2, and to simulate burnup cycle calculationsof them. From these burnup cycles, depletion matrices will be derived and used in future work tomodel realistic fuel cycle scenarios in a fuel cycle simulator, FANCSEE being developed at KTH.Figure 1 shows the process of transforming the results of the burnup cycles to usable data in theFANCSEE software.

The nuclear reactors under consideration will be the PressurizedWater Reactor LOVIISA-2 (VVER-440/213) in Finland and the Boiling Water Reactor FORSMARK-3 (ABB-III) in Sweden. Bothreactor types were some of the first reactors developed in the 1950s and have since been improvedon. They are also the most common reactor types which make them important for the FANCSEEproject.

3

Figure 1: The Figure shows parts of the process, from simulating the Burnup cycle to the finalsoftware. Red - Developers side, green - User side, blue - focus of the project.

2.3 FANCSEE

The FANCSEE simulator was initially developed at the Royal Institute of Technology in Stock-holm, KTH, by Torbjörn Bäck and prof. Wacław Gudowski, and during the years progressivelydeveloped in cooperation with the University of Tartu. The FANCSEE software is in its final stagesof development by Błażej Chmielarz, KTH [3], with a ultimate goal to create a user-friendly, graph-ically controlled software, allowing advanced simulations of a nuclear fuel cycle scenario, even forvery complex and diversified scenarios and has the ability of tracking 1307 nuclides at any point inthe cycle for up to 1100 years. The target audience would not only include scientists in the field,but also for policymakers as a source of well informed decisions and for educational purposes forstudents. It could also serve as a way of planning nuclear waste repositories, fuel needs, energygrid capacity or budget in a nuclear park.

3 Monte Carlo method

The Monte Carlo method is a stochastic method of approximating solutions to mathematical prob-lems using computational algorithms that rely on repeated statistical sampling from a probabilitydistribution function. It is often used in problems with probabilistic structures where variablesfollow a pattern. The randomness of the Monte Carlo method is only applied to the degrees offreedom that are defined and limited by the user, and have great advantages over a deterministicmethod mainly because it does not require new analytical solutions for every change of the system,is a more flexible method and can be applied to a system of any complexity.

4

4 Serpent 2

Serpent 2 is a universe based stochastic Monte-Carlo code [2]. It utilizes a user defined geometry(2-3 dimensional) with materials and their densities to simulate neutron reactions and transmu-tations of isotopes in a reactor core. The components that assemble the reactor are defined inseparate universes, ie. a space which is filled with a specific geometrical structure with a definedmaterial. What is returned by the simulation are the reaction rates of the neutrons, containingatomic density vectors of all isotopes present in a defined space over time, core criticality and isnot limited to keeping track of only the fuel.

The Serpent 2 code has predefined geometrical templates, parameters of which can be modified tosuit the desired design. The advantage of the universe based approach is that the overall modelingcan be divided into smaller modules, or universes, which can be stacked inside of eachother. Theuniverses are a central part of the Serpent 2 syntax; geometries and designs are created in separatelevels, independent of eachother and have the option to be nested within one another. Creating afunctionable model requires the scale of the nests to match so they don’t overlap or leave undefinedspace.

When running the script on a Unix-based system (Ubuntu 14.04), any geometrical errors, whetherit be overlaps or undefined regions, the Serpent 2 program automatically throws an error, thereforeit is critical to check the geometry before proceeding with the simulation. The Serpent 2 programalso checks the atomic densities of the materials used, ensuring that the syntax is correctly defined.It does not however check if the definitions are reasonable.

An example of how universes are nested within eachother is the design of the core. Usually,the highest universe in a core is occupied by the smallest component; the fuel pin, defining theuniverse as a single pin. The pin is then nested inside of a lower universe, called a pin lattice, akind of assembly composed of pins. The assembly is then placed in a assembly lattice making upthe core of the reactor, as shown in Figure 2. Within these universes, geometrical structures canbe implemented and filled with a desired material. Each material is comprised of a list of nuclideswhose properties, namely neutron interaction cross-sections at different temperatures, are definedin the cross-section library JEFF-3.1.

Figure 2: The Figure shows the different universes and how they are nested within eachother.

5

5 Theory

The theory of the nuclear reactor is based on just a few basic concepts, one of them being nuclearfission. Nuclear fission occurs when fissionable nuclides collide with an incoming neutron splittinginto fission fragments and high energy neutrons. The probability of which fission fragments thenuclide is split into is described in a fission yield distribution curve, shown in Figure 3 for 235U .The fragments are shown with their respective mass number. The fission fragments are releasedat high speeds, and it is the kinetic energy of the fission fragments that stands for the majorityof the power generated by a nuclear reactor. The high energy neutrons collide with other nuclidesand continue the chain reaction. To control the reaction rate, control rods, moderators, neutronpoisons are implemented. The design of the reactor is also of great importance to the reactionrate, as well as many other variables. Neutron poisons are often passive means of reaction control,mostly used when inserting fresh fuel loads into the reactor to lower its high reactivity. In the caseof the PWR, neutron poisons are always used with fresh fuel batches, decreasing in concentrationthroughout the cycle. Some neutron poisons get depleted when acquiring neutrons which couldin other cases initiate new reactions of other uranium atoms. Control rods work under the sameprinciple as neutron poison in the sense of catching neutrons, but are used for more finer reactioncontrol or complete reactor shutdown, SCRAM, in which case they are called Safety rods. Reactorshave several control rods, where some are kept in the core for an even reaction profile, and otherare kept outside of the core or partly inserted, depending on the desired rate of reaction and power.

Figure 3: The Figure shows a fission yield distribution curve of 235U . Graph: Nuclear Data Center,Japan Atomic Energy Agency.

6

5.1 Process of radioactive decay and transmutation through neutroninduced reactions

The driving force of nuclear energy released in a reactor comes from the transmutation of isotopesthrough neutron induced chain reactions. What is necessary for a chain reaction to occur is thepresence of fissile elements, those being elements that are able to emit neutrons when splittinginto fission fragments through a neutron induced fission reaction. When a neutron collides with anucleus of an element or a isotope it either gets captured, scattered or it splits the target nucleusand releases energy (more detailed explanation and other reactions are discussed in section 5.3Neutron Cross-section). If the target nucleus is split (fission), it emits neutrons and creates fissionfragments which are isotopes of lesser mass number, shown in Figure 3. In the right circumstancesthe emitted neutrons may lead to a chain reaction. In the other case, if the neutron is capturedby the nucleus of the isotope, it can stabilize it, transmute it to another stable nuclide or becomeunstable and decay with half-lifes stretching from milliseconds to billions of years. It can decay inmultiple ways and transmute into other isotopes, some of which are shown in Table 1. If

AZX,

where A is the mass number, Z is the atomic number and X is the symbol of the isotope. When aisotope decays its isotopic composition is changed, leading to the transmutation of a new isotope.This new isotope is called the daughter nucleus.

Mode of decay ParticlesDaughterNucleus

Alpha decay (α) An Alpha particle is emitted (A - 4, Z - 2)Beta minus (β−) An Electron and an electron antineutrino is emitted (A, Z + 1)Beta plus (β+) A Positron and an electron neutrino is emitted (A, Z - 1)Gamma decay (γ) A Gamma ray is released (A, Z)Neutron emission A neutron is emitted from the nucleus (A - 1, Z)

Electron capture The nucleus captures a orbiting electron and emits aneutrino; daughter nucleus is in a excited state (A, Z - 1)

Spontaneous fission Nucleus shatters into two or more smaller nuclei andother particles

-

Table 1: Table of common radioactive decays.

The most common isotope used in reactor fuel is 238U , being about 99% of the mass in naturallyfound uranium. But it is not a fissile isotope and therefore nuclear fuel is usually enriched withthe fissile isotope 235U , ie. the ratio of 235U to 238U is increased. 235U is exchanged betweenbatches of uranium in the enrichment process, resulting in a depleted and an enriched batch ofuranium. When the nucleus of 235U absorbs a neutron, it either undergoes fission, splitting intofissile elements (Equation 1) while emitting 2-3 neutrons or changes its isotopic composition whileemitting a Gamma ray (Equation 2). These neutrons can continue the chain reaction by reactingwith other fissile elements.

An example of how 235U reacts when absorbing a neutron is shown below, as well as their proba-bilities of happening.

23592 U +1

0 n→ 23692 U

85%−−→ fission(eg.13956 Ba+9436 Kr + 310n) (1)

7

23592 U +1

0 n→ 23692 U

15%−−→ 23692 U + γ (2)

When the non-fissile isotope 238U captures a neutron, it transmutes into 239U , then, by β− decay,into neptunium-239 (239Np) and finally into plutonium-239 (239Pu), which is a fissile isotope(Equation 3). This process of decay takes approximately 2.3 days and yields plutonium with ahalf-life of 24,110 years.

23892 U +1

0 n→23992 U

β−

−−→ 23993 Np

β−

−−→ 23994 Pu (3)

Because 239Pu is a fissile isotope it has a high probability of undergoing fission when absorbinga neutron. If it undergoes fission, it splits into fission fragments, free neutrons, gamma rays anda neutrino, with a released energy of ∼ 200 MeV; one example of how it can split is shown inEquation 4. While there is a 73% probability that the neutron absorption results in fission, thereis also a 27% chance that the neutron will be captured by the nucleus and therefore transmutewhile emitting a γ − ray as shown in Equation 5.

23994 Pu+1

0 n→ 24094 Pu

73%−−→ fission(eg.13454 Xe+10340 Zr + 310n) (4)

23994 Pu+1

0 n→ 24094 Pu

27%−−→ 24094 Pu+ γ (5)

The above Equations are examples of fission fragments, but during fission the produced isotopescan vary as long as the mass number N is balanced. Usually, fission generates two fission fragmentsand a number of neutrons, depending on the energy spectrum and fissioned nuclide.

5.2 Bateman Equation

The Bateman Equation is a basis for all codes dealing with isotope concentration in a nuclearreactor. It is a mathematical model describing the decay chains of isotopes, or in other words, byhow much an isotope decays into its daughter isotope as a function of time based on its decay rateand the initial amount of the decaying isotope. The mathematical model is a first order differentialEquation, on the form

dn

dt= −λn ⇒ n(t) = e−λtn(0) (6)

where n is the number of atoms of an isotope per unit volume [atoms/cm3], also known as atomicdensity and λ is the reaction rate for the change of one isotope to another. Equation 6 describesthe decay of a isotope dependent on time. Equation 7 however, is a system of first order differentialEquations, which describe the rate at which isotope i decays and at the same time is created fromall parent isotopes that transmute into it through nuclear decay or neutron induced fission. If theparent isotopes decay into its daughter isotope at a quicker rate than the daughter isotope decays,the net rate will be positive, increasing the density of the daughter isotope i.

dni(t)

dt= −λini(t) +

k∑j 6=i

λj→inj(t) (7)

While the calculation of the chain can be solved explicitly with k = 2, the Equation quicklybecomes nigh impossible to solve analytically in consequence of infinite chains and numerically

8

identical reaction rate constants. Therefore, the Bateman Equation can be applied with a varietyof state-of-the-art numerical methods, one of them being the Chebyshev Rational ApproximationMethod (CRAM) for solving the burnup matrix exponential, as computed in Serpent 2 and theFANCSEE software.

5.3 Neutron cross-section

Every material in the reactor consists of different elements whether it is an alloy, the moderator, orthe fuel. Material cross-sections are unique for every element and define the ability of the elementto absorb a neutron. If the material is neutron poison or a control rod, their material compositionneeds to contain elements with high neutron capture cross-sections, since they need to capture ex-cess neutrons. The opposite applies to the cladding material around the fuel pellets, it consists ofelements with low neutron absorption cross-section, so that neutrons emitted from the fuel passesthrough the cladding with minimal losses.

For nuclear reactors to work, neutron interactions have to occur in the core. The neutrons caneither interact with each other, neutron-neutron interactions, or with the nuclides in matter, thelatter being more probable. That is why neutron-neutron interactions are disregarded in the sim-ulations and only neutron-matter interactions are considered. The probability of the reactionshappening is described by the neutron cross-section and has the dimension of area [m2]. Thestandard unit for measuring cross-sections is the barn which is equivalent to 10−28 m2 [4] and canbe conceptualized by area of the target nucleus - the larger the area, the higher the probability ofa reaction. The neutron cross-section is given by

σ =µ

n(8)

where n is the atomic density of the target nuclide [atoms/m−3] (number of nuclides per unitvolume), and µ is the attenuation coefficient [m−1], ie. is characterizes how easily a material canbe penetrated by a neutron or energy. If the coefficient is small, the neutron passes through thenuclide as if it was transparent. The cross-section depends on the nuclide as well as on the type ofreaction and the kinetic energy of the incident neutron.

There is a large variety of reactions that can occur, but they can be divided into three groups- fission, capture and scattering. Neutron capture is a reaction in which an incident neutron isabsorbed by the nucleus and puts the nucleus in a excited state, followed by the nucleus decayingback into its initial ground state, emitting gamma rays. When a Neutron scattering reaction takesplace, the target nucleus emits a neutron after a neutron-nucleus interaction. If the energy ofthe incident neutron is transferred to the nucleus and the emitted neutron, the reaction is calledelastic neutron scattering. If the neutron is first absorbed by the nucleus, changing its kineticenergy to the internal energy of the nucleus and remitting the neutron at a significantly lowerenergy, leaving the nucleus in a excited state, it is called inelastic neutron scattering. Even thoughthe fission reaction is what leads to power generation, it can only do so under proper conditions.These conditions are partially ensured by neutron scattering events, making them as importantas fission. Lastly, the fission reaction is what emits energy that is used to generate the majorityof power. The incident neutron collides and enters the fissionable nucleus, creating a compoundnucleus and exciting it to energy levels well above the critical energy of the nucleus. This leadsto the splitting of the nucleus into two large fission fragments, releasing large amounts of energyin the process and 2 or 3 free neutrons that then interact with other nuclei. A sequence of fission

9

events is called a fission chain reaction.

If the kinetic energy of the incident neutron is too high (fast neutron), the probability of a fis-sion reaction is typically lower. Neutron cross-sections are highly dependent on incident neutronenergies. By scattering reactions the fast neutron decreases its kinetic energy until it reaches thethermal spectrum of energy. The process leads to a higher probability of interaction with a fissilenuclide in a fission reaction.

Figure 4 shows the neutron cross-section, in barn, of 235U nuclides dependent on the kineticenergy of the incident neutron. If the kinetic energy of the neutron is high (1 MeV), the neutroncross-section is significantly lower, in relation to the neutron cross-section when the neutron has alower kinetic energy (∼ 10−10 MeV). From the figure it can be seen that the neutron cross-section(in barn) is orders of magnitude greater when the neutron has a lower kinetic energy, than if itwould have a high kinetic energy. The region in the middle of the figure is called a resonanceregion, where the kinetic energy of the neutron and the energy of the target nucleus equal to aexcited state of a compound nucleus, creating peaks in the cross-section.

Figure 4: The figure shows the neutron cross-section of 235U dependent on the energy of theincident neutron (log scale). Graph generated from JANIS NEA with the JEFF-3.1 library.

5.4 Burnup calculation

The main goal of the burnup calculation is to simulate behavior of the reactor during its operationcycle. It is strongly based on the geometry created in the Serpent 2 script and the isotopic com-positions of the materials used. The geometry of the reactor stays the same during the simulationbut the material composition changes every step of the cycle, because of the neutron reactions orthe spontaneous radioactive decays taking place. Generally, burnup can be defined as a process ofnuclear transmutation under neutron flux, which simply means the amount of isotopes (or nuclides)resulting from radioactive decay or neutron-induced reactions.

Before the simulation is run, burnup steps and the amount of neutrons simulated per step arechosen (in our case 107 neutrons per step). These burnup steps can be chosen as the unit time[days] or MegaWatt days per ton of uranium produced [MWd/tU]. The amount of neutrons presentin the system will only be a fraction of the amount of neutrons in a real reactor. This is due tothe fact that it would be impractical to simulate such a large number of neutrons and, if thechosen number is sufficiently high, it does not cause deviations in the final burnup since Serpent 2normalizes the reaction rates from thermal power levels of existing reactors, entered by the user.

10

What Serpent 2 does is it sequentially assigns each simulated neutron energy, direction and distancetraveled, using the Monte-Carlo method that stochastically chooses each value from an appropriatedistribution. If the neutron enters a material, the code looks at a probability of the neutron reactingwith nuclides inside the material and the reaction that would occur, depending on the energy ofthe neutron and the cross-section of each reaction for every nuclide in the material. This algorithmis called Neutron transport. Every time a reaction occurs, the program takes note of the reactionand how many neutrons are emitted depending on the reaction. From those values it calculates therate of reaction normalized to either the specific power density of fuel, total power, neutron fluxor fission rate in the reactor. In the next step, the material nuclides are corrected, ie the nuclidesthat have decayed or have been reacted with, are changed and the simulation starts again withthese new nuclides for the next burnup step. Each burnup step, Serpent 2 creates a burnup matrixfrom the Bateman equation and reaction rates yielded by Neutron transport. The burnup matrixis then solved by Chebyshev Rational Approximation Method (CRAM).

5.5 Chebyshev Rational Approximation method

The Chebyshev Rational Approximation method (CRAM) is a state-of-the-art method of burnupmatrix exponential calculation. It is the most accurate method of approximation using the burnupmatrix, A, from the Bateman Equation, shown in Equation 9. Generally, the eigenvalues of theburnup matrix are confined to a region near the negative real axis, which is a condition well-suitedfor the use of Chebyshev polynomials, used in CRAM, when approximating the atomic density ofthe nuclides [3] [5].

When the burnup matrix and the initial atomic density of the nuclide is inserted into CRAM(Equation 10) it returns a function dependent on time of the burnup in the system.

n(t) = eAtn0(t) (9)

n(t) = a0n(0) + 2Re

k/2∑i=1

ai(At+ θiI)−1

n(0) (10)

Where n(0) is the initial atomic density of the isotope, A is the burnup matrix, a0 is a limitingvalue of the function at −∞, ai is the residue at the poles θi and I is the identity matrix.

CRAM is used in both Serpent 2 and FANCSEE. In Serpent 2, burnup matrices are used andhave both a predictor and corrector step, where the program first predicts the final atomic den-sity, followed by a correction based on the predicted value. For FANCSEE, cycle-averaged burnupmatrices are used and are assumed to not need further corrections, thus obtaining solutions muchquicker. This is done because the atomic densities from FANCSEE should be quickly approximatedas opposed to the computationally expensive corrected values returned by Serpent 2.

11

6 Method

Both reactor types where developed at the early stages of the nuclear age, but have been improvedsignificantly since then. The LOVIISA-2 reactor in Finland is a Pressurized Water Reactor (PWR)which was first designed in the former Soviet Union. The reactors pressure vessel contains thecore surrounded by water, which acts as a moderator and a coolant, being at a temperatureof approximately 300 degrees centigrade at a pressure of 123 bar. Because of the pressure, themoderator stays in its liquid form throughout the cycle and the system has therefore a primaryand a secondary heat transfer circuit. The primary circuit cycles the moderator through the core,cooling it, and then pumps it into the steam generator where it transfers its heat to the water inthe secondary circuit, before it enters the condenser. The water in the secondary circuit reachesits saturation point and transforms into steam before entering the turbine, as shown in Figure 5.

Figure 5: The figure shows a diagram of the Pressurized Water Reactor. Photo: United StatesNuclear Regulatory Commission.

FORSMARK-3 is a Boiling Water Reactor (BWR) which was developed by ASEA - Atom inSweden with a different approach to design. FORSMARK-3 has its vessel at a pressure of 69 bar,allowing the moderator to boil when passing through the core, generating steam and drying it,all inside of the reactor vessel, instead of having pressure in the vessel high enough to preventboiling, like LOVIISA-2. This requires only one circuit, so the steam passes through the turbineafter leaving the vessel. A simplified diagram is shown in Figure 6.

Figure 6: The figure shows a diagram of the Boiling Water Reactor. Photo: United States NuclearRegulatory Commission.

12

6.1 LOVIISA-2

The design of the core starts at the smallest components of the reactor - the fuel pellets. Theseneed to be designed in detail, considering both its inner and outer diameter of their cylindricalgeometry. The inner diameter and the concave geometry of the pellet is there for when the fuelpellets swell up during operation where the core of the pellet would expand the most because ofheat. The pellets are kept in a tube consisting of a thin layer of cladding material, in the case ofthe PWR VVER-440/213 LOVIISA-2 NPP the cladding is made out of a Zirconium-based alloysuch as Zr-1%Nb. The usage of Zirconium alloy is central, because it has many properties desirablefor cladding in a PWR - it has a good corrosion and stress-corrosion cracking resistance, whichis critical to the tube because the fuel rods oscillate under operation and should not crack. Thealloy also has low neutron absorption, ie. it does not absorb as many neutrons as other materials,therefore letting them interact with other isotopes and improving neutron economy (a balance ofneutrons created and neutrons lost).

The Zr-1%Nb tubes, or pins as they are called, are filled with fuel pellets stacked on top ofeach other to a certain point. The rest of the tube contains an area called the gas plenum, whereexists a spring that applies a force on the pins so that they remain stacked during transport. Thegas plenum serves as a margin for when the pressure increases inside of the pin, because of gasesreleased in fission. The rods are then assembled together into fuel assemblies in a hexagonal geom-etry, each assembly containing 126 fuel rods [6] with the same fuel and one central unfuelled rod,serving as a tube for water flow, with a rod pitch of 1.22 cm [6], ie the distance from the center ofone fuel rod to an adjacent rod (Figure 7).

Surrounding the rods and the Zr-1%Nb shroud containing the rods is a liquid, that works asa moderator and a coolant; in the PWR - water mixed with Boric acid at an average concentrationof 650 ppm. The role of the coolant in the reactor is to cool the core and transport the heatgenerated by it to the steam generator. The role of a moderator is to control the rate of reactionof fission, which, in the absence of the moderator would cease instantly, and in the absence ofcoolant cause core meltdown. The way to decrease the rate of reaction is to absorb some of theneutrons so that there is only a necessary amount of neutrons present for chain reaction criticality.This is done partly by the diluted Boric acid in the moderator, because of its neutron absorbingproperties. The moderator also serves as a way to slow down the neutrons because they need tobe moving at the right speed for other nuclides to capture them at the right rate, this is calledneutron capture (Neutron cross-section 5.3).

The moderator works as a passive reaction control, but for more active control a control rodis implemented. It can be used as a means of increasing/decreasing the reactor activity or aspart of a routine reactor shutdown for inspection by lowering or raising them. The geometry of acontrol rod is nearly identical to that of a fuel assembly, containing rods of the same size, but filledwith Borated Stainless Steel consisting of 20% chromium, 16% nickel and 2% natural boron, ie.a combination of 10B and 11B isotopes. The control rods are approximately twice as long as thefuel rods, with an upper part being the control rod and the lower part being fuel. They are thisway so that the fuel can still be inside of the core when the control rods are raised. Throughoutthe length of the core, 11 spacer grids of INCONEL alloy 690 are placed equidistantly. These aremesh like components of the core that hold the fuel pins in place so they do not oscillate duringoperation, caused by the moderator flowing through the core.

13

The reactor core is assembled by arranging 313 assemblies, where 276 are fuel assemblies and37 are control rod assemblies [6]. The core is also partly surrounded by shielding assemblies alsocalled reflectors. The reflectors role in the core is to limit neutron leakage, by reflecting the neu-trons back into the core instead of them leaving the core. The core is placed in a 6 cm thick corebarrel and then in a pressure vessel of thickness 14 cm, lined with a stainless steel clad on theinside (Figure 8 and 9).

Figure 7: The figure shows the cross-section of the VVER-440 fuel assembly. The distance betweenthe pins is called the pin pitch.

14

Figure 8: The figure shows the horizontal cross-section of the VVER-440 core.

15

Figure 9: The figure shows the vertical cross-section of the VVER-440 core, where the top part ofthe core is the gas plenum and the lines crossing the assemblies are the spacers.

16

6.2 FORSMARK-3

The FORSMARK-3 NPP has been updated and modernized a number of times to this date. Mostimportantly, the fuel type has been changed to SVEA-100 which differs from its older counterpartsby the 10x10 lattice (100 fuel rods per assembly) and no central tubes, with a larger pitch betweenthe fuel rods (1.24 cm). The pellet formed fuel used in the FORSMARK-3 reactor is uraniumdioxide (UO2), largely consisting of 238U with an enrichment of 2.769% 235U [7].

Surrounding the fuel is the cladding, which consists of Zirconium alloy-2, a low neutron cross-sectional material protecting the fuel from corrosion. The fuel pellets are stacked on top of eachother inside of the cladding tube and collected together into a 10x10 square fuel assembly andplaced in a sheath of Zirconium alloy-4.

The BWR also distinguishes itself from other reactors by arranging 4 fuel assemblies in a 2x2manner and having a control blade in the form of a cross between the assemblies (Figure 10). Thecontrol blade has 4 symmetrical arms, each containing 21 control rods[8]. The rods are made outof Boron Carbide powder inside of a Type-304 stainless steel tube, then placed inside of the bladeframe. When 169 assembly bundles are constructed, they are placed beside each other with a pitchof 15.45 cm, and 24 single assemblies are added symmetrically on the periphery [6]. When the coreis assembled it is placed in a 16 cm thick vessel, with an interior cladding of Type-304L stainlesssteel (Figures 11 and 12). The interior cladding is of great importance because it is usually thevessel that dictates the life of a reactor and the cladding protects the vessel from corrosion. If thecladding looses its integrity, leading to sufficient vessel damage, the reactor has to be decommis-sioned.

The vessel is filled to the top of the core with a moderator composed of demineralized water(as opposed to borated water in the PWR) and kept under a constant pressure of about 69 bar [6].Just above the core (Figure 12), water is pumped with a temperature of 278C and is then pumpeddown below the core by 8 equidistant jet pumps. The water then travels through the core, heatingup and becoming saturated steam when rising to the top part of the vessel. As the water risesthrough the core its density decreases as it heats up. The water is defined as layers with differentwater densities as seen in Figure 12 and Table 2. The steam then reaches the steam separatorwhere the larger droplets get separated from the steam and then the steam dryer before it goes outto the outlet into the turbine at a temperature of 286C. Because the water in the BWR is boiling,its density varies significantly throughout the core.

17

Figure 10: The figure shows the cross-section of the FORSMARK-3 SVEA-100 fuel assemblybundle with a control blade between the assemblies.

18

Figure 11: The figure shows the horizontal cross-section of the FORSMARK-3 ABB-III core.

19

Figure 12: The figure shows the vertical cross-section of the FORSMARK-3 ABB-III core. 1, 2, 3,4, 5, 6, 7 and 8 are regions of water with decreasing density.

Water ID Axial position in core [cm]Water density[g/cm3]

8 221 - 368 0.23797 191 - 221 0.31956 162 - 191 0.36625 132 - 162 0.42504 103 - 132 0.49953 74 - 103 0.59102 44 - 74 0.68991 0 - 44 0.7454

Table 2: Table of moderator density depending on axial position, from top of core to bottom (Coreheight 368 cm). See Figure 12. [9]

20

6.3 Implementation of core design to Serpent 2 syntax

Converting the core design into Serpent 2 code needs detailed geometry and correct materialdefinitions to function correctly and provide correct data. All commands used in Serpent 2 arecalled cards. The Serpent 2 code has a predefined card that describes the fuel rods called PIN.The PIN card is a universe type that is defined by an string, followed by a material name and theradius of the cross-sectional area that is occupied by the material. Here the PIN is defined as ’pin1’, followed by the materials and their parameters.

When the necessary PINS are defined, they are assembled together using the LATTICE card,which input is its identification string, lattice type (1 - square, 2 - hex), (x,y) position of center,number of elements in both x and y direction and pitch, which defines the distance between eachelement. Filler pins, ’pin 3’, are added on the periphery. These are usually a moderator thatsurrounds the elements of the core. In this case ’lat 10’ is the LATTICE containing 126 ’pin 1’universes, surrounded by ’pin 3’ which is just defined as water. A single Central tube, ’pin 2’ isplaced in the middle of the hexagon.

Universes, either defined as pins or lattices, have the option to add geometrical structures dependingon their design. In the VVER-440, the fuel assemblies are placed inside of a hexagonal shroudtube, defined by the SURFACE card. The SURFACE card itself only defines a surface that can beused in multiple universes, but until the CELL card is used, the geometry is not generated. TheSURF card defines the boundary of the geometry, its (x,y) position and the radius of the surface,while the CELL fills the surfaces with a universe, in this case - ’lat 10’. It then fills the othersurfaces with materials by defining the area between surface 1 and surface 2 with the parameter1, defining the area outside of surface 1 and −2, the area inside of surface 2.

21

Every material used in the geometry has a defined nuclide composition. This is done by the MATE-RIAL card, containing parameters of nuclide names and their densities. The material identificationis on the form ZZAAA.XXc, where ’ZZ’ is the integer proton number (1-99), ’AAA’ is the massnumber and ’XX’ is the temperature estimate at which the isotope will be present in the material,expressed in hundreds of Kelvin. On the same line, the density of the isotope is expressed, either inatomic density (1024/cm3) or mass density (g/cm3). Combining the isotopes, the material ’fuel1’defines one of the fuels in the core.

With these basic cards, the whole geometry of the reactor core can be generated. For morecomplicated designs, surface cards are stacked, creating desired geometries.

22

6.4 Reactor profile

For accurate reaction rates and a balanced core activity, the fuel should not be approximated to beidentical in the core, but instead have the most fresh fuel at the periphery, and the most burned inthe center of the core, with intermediate burned fuel between the two. After a full burnup cycle,the central, most burned fuel is taken out and replaced by the second to oldest fuel and so on. Thepreviously fresh fuel on the periphery is moved towards the center of the core and its place it takenby new fresh fuel, as shown in Figure 13. If the core would have identical fuel, there would bemore reactions occurring in the center of the core, which is preferably avoided because a flattenedcore power profile is desired.

Figure 13: The blue arrows signify the movement of the nuclear fuel at the end of each fuel cycle.Different colors of the fuel indicate each zone (FORSMARK-3).

For the LOVIISA-2 reactor, the core is divided into 3 radial zones containing differently burnednuclear fuels. The fuels are moved at the end of each 12-month burnup cycle. In the case ofFORSMARK-3, the reactor has 5 radial zones and a burnup cycle of 12 months.

Determining the level at which the nuclear fuels are burned for the different zones is done bysimulating the cycle with identical fuel for all zones followed by a calculation of how much of eachzone has been depleted. The fuel is then redefined in the Serpent 2 code and placed in appropriatezones for a realistic power profile. The profile is essential for acquiring accurate reaction ratesthroughout the core.

The reactor core models presented were made to resemble the actual reactor cores active today.Components that were either confidential or protected by Patents were approximated, as well ascomponents the exact geometry of which did not affect the reactions directly.

23

7 Results

After modeling the core and running it through burnup simulation, the Serpent 2 code returnsseveral files. It returns files with depletion matrices for all fuel batches and control rods at eachburnup step, using a debugging function of the code. It also returns a set of mesh figures showingneutron reactions at every burnup step in the reactor core. The depletion file consists of a setof matrices with rows corresponding to different isotopes and columns to burnup steps. Thesematrices keep data of atomic density, mass density, fission rates, cross-sections and other usefuldata.

For each reactor the simulations have been done for 10 different power levels ranging from 60%to 105% with interval steps of 5%, where 100% power level corresponds to the Thermal Capacityof each reactor [MWt]. In some circumstances the power grid can only take a certain amount ofpower from the reactor or the reactivity in the core is too low to sustain the thermal capacity sothe power level is adjusted accordingly.

The results below show the atomic densities of some of the most interesting nuclides in one fuelbatch. The figures follow the life cycle of one fuel batch, from the point at which it entered thereactor core as fresh fuel (on the periphery) and was taken out as burned fuel, each year replacedcloser to the center of the core (6.4 Reactor profile). This explains the sudden discontinuities every365 days for the SERPENT data in the following figures.

7.1 LOVIISA-2

The Thermal Capacity of the LOVIISA-2 reactor is 1500MWt [10]. The Figures below show thechange in atomic density for some of the most important isotopes as a function of time in days forone full reactor cycle operating at 100% Thermal Capacity.

Figure 14: The figure shows a comparison of data calculated by Serpent 2 (corrected) and FANC-SEE (predicted) of how the isotope concentration for 235U and 239Pu changes through transmu-tation in the LOVIISA-2 reactor at 100% power.

Figure 14 depicts the atomic density of 235U and 239Pu as a function of time (in days) for a full lifecycle of one fuel batch generated by Serpent 2 (shown by the discontinuous line) and FANCSEE

24

(continuous line). It shows how the atomic density of 235U decreases throughout the cycle, whilstthe atomic density of 239Pu increases. The reason is that the 235U nuclides undergo fission throughneutron induced reactions, causing a nuclear chain reaction and creating nuclides of lower atomicmass. During a chain reaction, neutrons collide with 238U starting a process of transmutation andsince its transmutation chain leads to 239Pu, as shown in section 5.1 (Equation 3), the atomic den-sity of 239Pu increases substantially in the early days of the cycle. Eventually, the curve of 239Pu

balances out as a cause of it transmuting into other nuclides through both fission and capture andthe gradual depletion of 235U (5.1 Equation 4 and 5).

Figure 15 shows additional nuclides of 240Pu and 241Am on a semi-logarithmic graph. Here,240Pu is a nuclide created as a result of the transmutation of 239Pu among others. The only nu-clides that are diminishing are 235U and 238U (not shown), being the initial nuclides in the fuel. Itcan also be noted that there is a similarity between the curves of 240Pu, 241Am and their substrateisotope 239Pu, the transmutation rate of which they are highly dependent on.

Comparing the Serpent 2 generated data with the FANCSEE data reveals similarities betweenthe two. In Serpent 2, the CRAM uses the burnup matrix exponential, both when predicting andcorrecting the reaction rates of the isotopes. This yields a precise graph of the atomic densitiesbut has a long computational time (approximately 8 hours, depending on the number of neutronsin the system). In FANCSEE, the CRAM uses a cycle-averaged burnup matrix exponential andonly computes the predicted value. It assumes that the prediction does not need corrections,thus obtaining solutions in the matter of seconds. If the prediction through FANCSEE is wellapproximated to the corrected data from Serpent 2, the model of the reactor can be used in theFANCSEE program. As seen in Figures 14 and 15 for the LOVIISA-2 reactor, the prediction fromthe cycle-averaged burnup matrix exponential in FANCSEE is well approximated to the correctedgraph from Serpent 2, especially at both the point of origin and the end value.

Figure 15: The figure shows a comparison of data calculated by Serpent 2 (corrected) and FANC-SEE (predicted) of how the isotope concentration for 235U , 239Pu, 240Pu and 241Am changesthrough transmutation in the LOVIISA-2 reactor at 100% power.

Figure 16 is a mesh depicting the fission density of the LOVIISA-2 core at burnup step 9.80MWd/kgU, ie. where the process of fission occurs. As expected, the fission density is most densein the central part of the core, signified by the red and yellow colors. The blue scale is neutron

25

thermalization, caused by neutron scattering in the moderator, where the neutrons loose kineticenergy and enter the thermal spectrum, provided that they are not absorbed. These figures aregenerated for each burnup step and are convenient when analyzing the core model. If there are anydiscrepancies in the result, they can provide useful information on where the problem may reside.

Figure 16: The figure shows neutron thermalization (blue scale) and fission density (yellow scale)in a cross-section of the LOVIISA-2 reactor core at 100% power.

7.2 FORSMARK-3

The FORSMARK-3 has a Thermal Capacity of 3300MWt [11] and the following figures depict thechange in atomic density of some of the most important isotopes as a function of time in days inthe reactor operating at 100% Thermal Capacity.

Figure 17 and 18 show the atomic densities for 235U , 239Pu and 235U , 239Pu, 240Pu and 241Am

respectively. 238U is less relevant, so it is left out, since it mostly becomes 239Pu. The graph is alsolimited between 50 days and 1825 days, for clearer results. Here, it can be seen that the averageatomic density from FANCSEE is poorly approximated to the data generated from Serpent 2. Notonly is the approximation poor, but also the overall result of the atomic density. What the figuresseem to show, is that in the middle of the fuel cycle the 235U and 239Pu isotopes cross, increasingthe content of 239Pu to levels higher than 235U . The overall increase of fissile nuclides, observablein Figure 17, is typical to a breeder reactor, which breeds more fissile isotopes than it consumes(ie. having a conversion ratio > 1). The FORSMARK-3 reactor is not a breeder reactor and theatomic densities of 235U and 239Pu should approach each other at a later point, towards the endof the cycle.

26

Figure 17: The figure shows a comparison of data calculated by Serpent 2 and FANCSEE of howthe isotope concentration for 235U and 239Pu changes through transmutation in the FORSMARK-3reactor at 100% power.

Figure 18: The figure shows a comparison of data calculated by Serpent 2 and FANCSEE of howthe isotope concentration for 235U , 239Pu, 240Pu and 241Am changes through transmutation inthe FORSMARK-3 reactor at 100% power.

Looking at Figure 19, the fission density of a axial cross-section of the FORSMARK-3, it can benoted that the top part of the core has relatively low reactivity. This is caused mainly by therelatively low density of the moderator at the top part of the core. In a BWR, the moderator(water) boils, when passing through the core. This was implemented in the model by dividing themoderator into different levels with densities varying from high density at the lower part of the coreand low density at the top part of the core (see Table 2, 6.2 FORSMARK-3). The moderator notonly serves as a way to transport heat away from the core, but also as a medium for fast neutronsto scatter against, losing kinetic energy and entering the thermal spectrum. When the neutronsare in the thermal spectrum, the fission cross-section of 235U , which is dependent on the energy

27

of the incident neutron, is high enough to the point that a fission reaction can take place. Butbecause of the low density of the moderator in the top part of the core, these scattering reactionsoccur at significantly lower rates compared to the rest of the core. Low scattering reactions leadto few thermal neutrons created, which then results in even fewer fission reactions.The position of the control blades were also altered to increase the rate of reaction in the core,but so far it too yielded flawed results. The produced figures came from the simulation of theFORSMARK-3 reactor core with fully withdrawn control blades.

Figure 19: The figure shows neutron thermalization (blue scale) and fission density (yellow scale)in a cross-section of the FORSMARK-3 reactor core at 100% power. (Axial cross-section)

8 Conclusions

The process of modeling reactors and simulating them in the Monte-Carlo code Serpent 2, showedto be an appropriate way to approximate the atomic densities of the nuclides present in the core.

From the results, it could be seen that the LOVIISA-2 core model generated appropriate re-sults when it was used to simulate fuel burnup. The simple geometry of the reactor core anda homogeneous moderator throughout the core yielded an even reaction profile. In the case ofFORSMARK-3, the moderator had different densities depending on its axial position. This turnedout to be problematic because of the approximations made when defining the densities of themoderator, which resulted in incorrect reaction rates throughout the core. The problem may alsolay in the definition of the initial fuel composition in the Serpent 2 code - the enrichment of 235U

may have been too low, resulting in less reactions occurring in the core. A detailed revision of alldefined materials in the syntax will be made.

The fully functional models will prove to be important to the FANCSEE program when simu-lating the fuel cycle, since they can be used to estimate and calibrate the burnup of the nuclidesin the nuclear fuel.

Tracking nuclear fuel through its cycle should become easy with the help of a graphical userinterface in the FANCSEE software.

9 Future work

Future work for the FANCSEE project will include the modeling of more reactor cores in the Ser-pent 2 code, creating a database of burnup matrices for the software. Some of the major remainingreactor types to model are the Pressurized Heavy water Reactor, Gas Cooled Reactor, Light WaterGraphite Reactor and a Fast Breeder Reactor.

28

A correction of the FORSMARK-3 BWR code is also a priority, since it is a reactor of highinterest - a reactor of this type (BWR) is used widely throughout the word and is the dominantreactor type in Sweden.

10 Acknowledgements

I want to thank Prof. Wacław Gudowski for giving me the opportunity to be a part of the FANC-SEE project and for the feedback that I received throughout the duration of my work.

I also want to give huge thanks to Błażej Chmielarz for aiding me in any difficulties that came upduring the development of the models as well as the feedback and guidance that I received.

References

[1] http://www.tsp-data-portal.org/Breakdown-of-Electricity-Generation-by-Energy-Source#tspQvChart (2014)

[2] Leppänen, J., ’Serpent – a Continuous-energy Monte Carlo Reactor Physics Burnup CalculationCode.’ VTT Technical Research Centre of Finland. (June 18, 2015).

[3] Chmielarz, B., ’FANCSEE - Development of a fuel cycle code with Graphical User Interface’,PhD Thesis, Kungliga Tekniska högskolan – Nuclear Reactor Physics, 2016.

[4] http://www.nuclear-power.net/neutron-cross-section/

[5] M. PUSA and J. LEPPÄNEN, “Computing the Matrix Exponential in Burnup Calculations,”Nucl. Sci. Eng., 164, 140 (2010).

[6] J. Varley, ’World nuclear industry handbook, 1999’, Nuclear Engineering International., Wilm-ington Business Pub. 1999, p. 194, 205, 215.

[7] L. Agrenius, ’Clab – Data for fuel assemblies used in calorimetric and nuclear measurements’,Agrenius Ingenjörsbyrå AB, (December 2006)

[8] S. Palmtag, ’Initial Boiling Water Reactor (BWR) Input Specifications’, U.S. Department ofEnergy, Nuclear Energy, (February 28, 2015)

[9] W.Marshall, B. J. Ade, S. Bowman, J. S. Martinez-Gonzalez, ’Axial Moderator Density Dis-tributions, Control Blade Usage, and Axial Burnup Distributions for Extended BWR BurnupCredit’, Oak Ridge national Laboratory, p. 32 (modified) (August 2016)

[10] https://www.iaea.org/pris/CountryStatistics/ReactorDetails.aspx?current=158

[11] https://www.iaea.org/pris/CountryStatistics/ReactorDetails.aspx?current=532

[12] https://www.nrc.gov/docs/ML1125/ML11258A302.pdf

29