Neutron diffusion equation solution using finite element method

Este trabalho apresenta a implementação da difusão neutrônica no seio do programa MEF usando como dados de entrada, as seções de choque homogenizadas produzidas pelo programa HAMMER. Esta implementação resgata a possibilidade de se fazer o cálculo de difusão neutrônica utilizando-se o método dos elementos finitos, tarefa já realizada anteriormente no IEN[2]. A esta implementação, acrescenta-se o cálculo de erro do fluxo neutrônico e se redefine a discretização, baseada neste erro, utilizando-se a remalhagem adaptativa[6] Este acréscimo permite homogeneizar o erro do cálculo do fluxo em todo o domínio tornando-se um programa moderno e um bom substituto ao programa CITATION. Apesar da utilização do método de potência, considerado um método de convergência lenta, para a solução do problema especial de autovalor, os resultados tem mostrado que isto não é uma limitação para a obtenção do fluxo neutrônico mesmo quando se utilizam malhas com mais de 10.000 elementos finitos. Neste trabalho, apresenta-se inicialmente a formulação variacional do problema de difusão neutrônica, a discretização em elementos finitos e finalmente a solução de vários exemplos de difusão neutrônica limitados a dois grupos de energia com várias regiões