Simulação Computacional da Estabilização de Sistemas Subcríticos segundo o Modelo Unidimensional de Difusão de Nêutrons Monoenergético

Abstract

Com o objetivo de descrever o perfil do fluxo escalar de nêutrons monoenergéticos em geometria unidimensional usando o modelo de difusão, desenvolvemos um método numérico espectronodal de difusão (END) em duas versões. A primeira versão considera problemas de autovalor, cujas soluções numéricas geradas pelo método END, sendo livres de erros de truncamento espacial, coincidem com a solução analítica dominante, afora erros da aritmética finita computacional. A segunda versão considera problemas de fonte fixa, onde a fonte é obtida analiticamente a partir da reconstrução espacial do fluxo escalar gerado pelo método END para problemas de autovalor. Neste caso também são obtidas soluções livres de erros de truncamento espacial. Um código computacional que implementa o método END foi desenvolvido em linguagem de programação C# (C Sharp), a fim de gerar as soluções numéricas. A titulo de ilustração, são apresentadas soluções numéricas, geradas através deste código.