Soluςão da Equaςão de Difusςão de Nêutrons em Geometria Retangular pela Técnica da Diagonalizaςão de Matrizes

Abstract

Recentemente (Petersen et.al., Kerntechnik, 74:60-64 (2009)) resolveu de forma numérica as equações de difusão de nêutrons com dois grupos de energia para uma placa heterogênea utilizando o método de diferenças finitas com condições de contorno do tipo albedo. Neste trabalho nós propomos a resolução do mesmo problema, porém de forma analítica, utilizando a técnica de diagonalização de matrizes. Partindo de um fator de multiplicação prescrito (keff), com a técnica da diagonalização nós desacoplamos o sistema de equações diferenciais e resolvemos separadamente cada equação, tanto para o fluxo rápido de nêutrons quanto para o térmico, e recalculamos o keff pelo método da potência até um critério de parada escolhido. Essas equações foram resolvidas utilizando resultados analíticos conhecidos para equações diferenciais ordinárias de segunda ordem. Para a determinação das constantes geradas foram utilizadas as condições de contorno e de continuidade de corrente e fluxo de nêutrons nas interfaces. Os resultados obtidos foram comparados com a solução numérica proposta por (Petersen et.al., Kerntechnik, 74:60-64 (2009)) para duas regiões.