An algorithm for multi-group two-dimensional neutron diffusion kinetics in nuclear reactor cores

Abstract

The objective of this thesis is to introduce a new methodology for two{dimensional multi{ group neutron diffusion kinetics in a reactor core. The presented methodology uses a polyno- mial approximation in a rectangular homogeneous domain with non{homogeneous boundary conditions. As it consists on a truncated Taylor series, its error estimates varies with the size of the rectangle. The coefficients are obtained mainly by their relations with the independent term, which is determined by the differential equation. These relations are obtained by the boundary conditions only, and these relations are proven linear independent. A numerical scheme is made to assure faster convergence. The procedures done for one homogeneous rectangle are used to construct the solution of global orthogonal geometry with step{wise constant parameters steady state and time dependent problems by the iterative SOR algo- rithm. The dominant eigenvalue and its eigenfunction are obtained by the power method in the eigenvalue problem. The solution for the time dependent cases uses the modi ed Euler method in the time variable. Four classic test cases are considered for illustration.

O objetivo desta tese é introduzir uma nova metodologia para a cinética bidimensional multi- grupo de difusão de nêutrons em reatores nucleares. A metodologia apresentada usa uma aproximação polinomial em um domínio homogêneo retangular com condições de contornos não homogêneas. Como ela consiste em uma série de Taylor truncada, sua estimativa de erro varia de acordo com o tamanho do retângulo. Os coeficientes são obtidos principalmente pelas suas relações com o termo independente, que _e determinado pela equação diferencial. Estas relações são obtidas apenas pelas condições de contorno, e é demonstrado serem linearmente independentes. Um esquema numérico é feito para assegurar uma rápida convergência. Estes procedimentos feitos para um retângulo homogêneo são feitos para construir soluções para problemas de autovalor e dependentes do tempo de geometria ortogonal global com parâmetros seccionalmente constantes pelo método iterativo SOR. O autovalor dominante e sua autofunção são obtidos pelo método da potência no problema de autovalor. A solução para casos dependentes do tempo usam o método de Euler modificado na variável tempo. Quatro casos-teste clássicos são considerados para ilustração.