TY - JOUR

T1 - Development of time-dependent reaction rates to optimise predictor–corrector algorithm in ALEPH burn-up code

AU - Fiorito, Luca

AU - Stankovskiy, Alexey

AU - Van den Eynde, Gert

AU - Labeau, Pierre-Etienne

A2 - Kalcheva, Silva

N1 - Score = 10

PY - 2013/5/31

Y1 - 2013/5/31

N2 - Shells coupling Monte-Carlo transport and deterministic depletion codes are extensively used in the nuclear field to simulate material changes throughout irradiation. The dynamic behaviour of the phenomenon is described by the system of coupled ordinary differential equations, with generally a stiff matrix of coefficients that current codes keep constant in time along each burn-up interval. The matrix coefficients represent decay constants and microscopic reaction rates of the numerous nuclides involved in the calculations. For a typical burn-up problem, their determination consumes most of the required computational time while only a small fraction is spent by the depletion solver.
This work presents a unique and innovative feature of the ALEPH Monte-Carlo burn-up code which optimises the depletion algorithm by using time-dependent matrix coefficients. Linear polynomials interpolate the evolution of the matrix coefficients along a few consecutive time steps. Then, trend curves are
constructed and used to extrapolate the effective reaction rates in the following intervals, thus reducing the total required computational time . This technique has been implemented in the version 2 of the ALEPH Monte-Carlo burn-up code and validated against the REBUS experimental benchmark. The results revealed a considerable computational time saving without any drawbacks in the accuracy.

AB - Shells coupling Monte-Carlo transport and deterministic depletion codes are extensively used in the nuclear field to simulate material changes throughout irradiation. The dynamic behaviour of the phenomenon is described by the system of coupled ordinary differential equations, with generally a stiff matrix of coefficients that current codes keep constant in time along each burn-up interval. The matrix coefficients represent decay constants and microscopic reaction rates of the numerous nuclides involved in the calculations. For a typical burn-up problem, their determination consumes most of the required computational time while only a small fraction is spent by the depletion solver.
This work presents a unique and innovative feature of the ALEPH Monte-Carlo burn-up code which optimises the depletion algorithm by using time-dependent matrix coefficients. Linear polynomials interpolate the evolution of the matrix coefficients along a few consecutive time steps. Then, trend curves are
constructed and used to extrapolate the effective reaction rates in the following intervals, thus reducing the total required computational time . This technique has been implemented in the version 2 of the ALEPH Monte-Carlo burn-up code and validated against the REBUS experimental benchmark. The results revealed a considerable computational time saving without any drawbacks in the accuracy.

KW - reaction rate

KW - burn-up code

KW - linear interpolation

UR - http://ecm.sckcen.be/OTCS/llisapi.dll/open/ezp_133085

UR - http://knowledgecentre.sckcen.be/so2/bibref/11013

U2 - 10.1016/j.anucene.2013.05.046

DO - 10.1016/j.anucene.2013.05.046

M3 - Article

SN - 0306-4549

VL - 62

SP - 307

EP - 315

JO - Annals of nuclear energy

JF - Annals of nuclear energy

IS - 1

ER -