Scribd
Upload
529 views

Lessons From..BURNCAL

BURNCAL is a Monte Carlo burn up code that: 1) Uses neutron absorption and reaction tallies from MCNP to directly determine changes in nuclide inventory, rather than using neutron flux as an intermediate input. 2) Solves the coupled differential equations that describe the production and depletion of isotopes over time. 3) Allows full use of MCNP's capabilities by integrating all neutron flux information into burnup calculations without post-processing or flux manipulation.

Uploaded by

Topan Setiadipura
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
529 views

Lessons From..BURNCAL

BURNCAL is a Monte Carlo burn up code that: 1) Uses neutron absorption and reaction tallies from MCNP to directly determine changes in nuclide inventory, rather than using neutron flux as an intermediate input. 2) Solves the coupled differential equations that describe the production and depletion of isotopes over time. 3) Allows full use of MCNP's capabilities by integrating all neutron flux information into burnup calculations without post-processing or flux manipulation.

Uploaded by

Topan Setiadipura
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
You are on page 1/ 23

Monte Carlo Burn Up Code

Computational Division
BATAN – INDONESIA
December 2007
(latest update 11 April 2007)
Motivation
“Compared to deterministic codes, the MC-based
depletion code has broader application, and
sometimes plays a unique role in research. As a
practical example, in the early 1990s, the Institute of
Nuclear Energy Research (INER) in Taiwan built a
research reactor which would require certain tools to
perform depletion calculations. MCNP, being the ideal
tool of modelling the research reactor, cannot do
depletion calculations. Therefore, they purchased a
linkaged code (similar to MCODE) from Germany at
the price of about $200,000”
[Zhiwen Xu,”Design Strategies for Optimizing High Burnup Fuel in PWR”, MIT
Dissertation, January 2003]
Lessons from
BURNCAL: A Nuclear Reactor Burnup Code Using MCNP Tallies
Edward J. Parma
Advanced Nuclear Concepts
Sandia National Laboratories

(latest update 11 April 2007)


Abstract

• BURNCAL is a Fortran computer code designed to aid in analysis,


prediction, and optimization of fuel burnup performance in a nuclear
reactor.
• The code uses output parameters generated by the MCNP to
determine the isotopic inventory as a function of time and power
density.
• BURNCAL is unique in comparison to other burnup codes in that it
does not use the calculated neutron flux as input to other computer
codes to generate the nuclide mixture for the next step, instead,
directly uses the neutron absorption tally/reaction information
generated by MCNP for each nuclide of interest to determine the
nuclide inventory for that region.
• The above uniqueness allows for the full capabilities of MCNP to be
incorporated into the calculation and a more accurate and robust
analysis to be perform.
Introduction

• Burnup and depletion codes have been developed and


used in the nuclear industry since the introduction of digital
computing.
• These codes solve the diffusion equation in one to three
dimensions using a few neutron-energy groups. Only a few
of the major fission products were included in the
calculation.
• As the digital computing capabilities and the neutron
cross-section libraries developed, more sophisticated
solution techniques, such as discrete-ordinates and Monte
Carlo methods, were employed to solve for keff and the
multigroup neutron flux.
Introduction

• Codes used today to calculate burnup reactivity


effects include coupled neutronics/burnup codes
and post-processing modules :
– neutronic : calculate keff and the energy-dependent
neutron flux for a reactor or a region within the reactor.
– burnup : uses the neutron flux, loading composition,
and power density to compute the nuclide composition
for the next time step.
The processes are repeated for the duration of the
burnup cycle.
Introduction

• The highes fidelity approach simply uses neutron


absoption and fission reaction information generated in
the neutronics code to determine the nuclide composition
at the next step.
• This type of model allows for all of the neutron flux
information to be integrated into the calculation without
post-processing and additional manipulation of the
neutron flux and cross-section set.
• Neutron absorption and fission reaction data for
individual nuclides are available as output for MCNP
through the use of ‘tallies’.
• The only requirement is that an energy-dependent cross-
section set is available for each nuclide of interest at the
required material temperature.
Introduction

• The Fortran code, BURNCAL, was developed to


incorporate the neutron absorption reaction output
available using the MCNP tally specification cards into a
burnup/production calculative model.
• The code uses the nuclide inventory, MCNP tally
information, power density, fission-product yield data,
and half-life data to determine the new nuclide inventory
for a given region at the new time. The new inventories
are then automatically placed back into the MCNP input
file and the case run for the next time step.
• At the completion of an analysis, the last-inventory can
be re-shuffle as a new loading scheme.
BURNCAL
• Code flow :
BURNCAL

• Design goals :
“to construct a burnup cade that would use
only the neutron absorption tally/reaction
information generated by MCNP in the
calculation of fissioned or neutron-
transmuted isotopes for different fuel-
loaded regions.”
BURNCAL
• Code packaged requirements:

– able to read and interpret an MCNP input file such that


isotopic concentrations can be identified for different
fueled and unfueled regions of a core.
– Able to read and interpret an MCNP output file such that
fission and transmutation reaction can be resolved for
each isotope and region of interest.
– Able to numerically calculate the production and removal
of a specific isotope through fission, transmutation, and
radioactive decay.
– Able to read, interpret, and write to an MCNP input file
such that the resultant isotopic cencentrations can be
reinserted into the MCNP input file at the proper locations.
PHYSICS_OFBURNUPPLUS_MATHEMATIC

• Production/removal equations :

Using the above rate information, a general coupled set of linear first-order differential
equations can be developed for all of the isotopes of interest. Due to the fact that a large
number of coupled equations can result, a numerical solution technique is the only
practical solution method.
PHYSICS_OFBURNUPPLUS_MATHEMATIC

• NEUTRON INTERACTION :
PHYSICS_OFBURNUPPLUS_MATHEMATIC

• Explicit diff.equation for each nuclide :

Time rate Production rate per unit volume of Removal rate per unit
of change isotope i volume of isotope i
in
concentrat
ion of from fission of from neutron from decay by by by
isotope i all fissionable transmutatio of all isotopes fissions neutron decay
nuclides n of all including capture
isotopes beta,alpha, (excludin
including etc. g fission)
(n,gamma),
(n,alpha), etc.
PHYSICS_OFBURNUPPLUS_MATHEMATIC
Fission production rate
Concentration of
isotope j , atoms/cm3
Neutron flux ,
All n/cm2-s
fissionable
nuclide

Microscopic fissions cross-


Fission fraction for the section of isotope j, cm2
production of isotope i from
fission of isotope j , n/cm2-s
PHYSICS_OFBURNUPPLUS_MATHEMATIC
Capture production rate
Concentration of
isotope k , atoms/cm3

Neutron flux ,
n/cm2-s

Integrated
over all Microscopic capture cross-
isotopes section (minus fissions) of
isotope k, cm2
PHYSICS_OFBURNUPPLUS_MATHEMATIC

Decay production rate Concentration of


isotope l , atoms/cm3

Integrate
over all Decay constant (equal to
isotopes ln(2)/half-life) of isotope l, 1/s
PHYSICS_OFBURNUPPLUS_MATHEMATIC
Concentration of
Fission removal rate isotope l , atoms/cm3

Neutron flux ,
n/cm2-s

Microscopic fissions cross-


section of isotope j, cm2
PHYSICS_OFBURNUPPLUS_MATHEMATIC

Capture removal rate Concentration of


isotope l , atoms/cm3

Neutron flux ,
n/cm2-s
Microscopic capture cross-
section of isotope j, cm2
PHYSICS_OFBURNUPPLUS_MATHEMATIC

Decay removal rate Concentration of


isotope l , atoms/cm3

Decay constant (equal to


ln(2)/half-life) of isotope l, 1/s
PHYSICS_OFBURNUPPLUS_MATHEMATIC
• Note that, although not explicitly defined as such, the neutron flux and cross sections are energy dependent. While, the
neutron flux and isotope concentration are both position and time dependent.

• Hence, the flux, cross section, and


concentration defined as :
φ ⇒φ( x, y , z , E , t )
σ ⇒σ( E )
N i ⇒N i ( x, y , z , t )
• Since only the flux and the cross sections are energy dependent, the
energy dependence can be integrated out of the differential equation as
follow
∞
σNφ ⇒ N ( x, y, z , t ) ∫σ ( E ) ⋅φ( x, y, z , t ) ⋅ dE
0
Computational Solution

• BURNCAL solve the above coupled set of linear first order differential
equation using 1st order Runge-Kutta Methode with sufficiently small
time steps to ensure adequate accuracy in the results.
Initial value problems of ODEs
Possible ‘better’ solution on the above problem :
1. Use Predictor-Corrector MCODE from MIT.
“ In fact the standard predictor-corrector algorithm is the proven
algorithm for burnup calculation for all licensing-level reactor
PARALLEL
physics codes…. It is conceptually simpler and numerically COMPUTATION
more accurate than the middle-timestep approach.” [Zhiwen ON
Xu,2003]
PUBLIC
2. Divide the complex chain to simple chain and solve it CLUSTER
analitically.[Zaki Su’ud,1996]
3. Matrix exponential method [Origen2.1].
Scientific Efforts

• Advance solution of initial value problem


of ODEs.
• Effective coupling between MC and burn
up calculation.
• Effectively run the above ideas using
parallel computing scheme on public
cluster.
• Application to find an optimum neutronic
design of an advance nuclear reactor such
as VHTR.

You might also like