Solução do problema de Cinética de Difusão de Nêutrons

Neste trabalho, resolve-se a Equação de Cinética Espacial da teoria de Difusão de Nêutrons em geometria cartesiana unidimensional, domínio homogêneo e heterogêneo, com dois grupos de energia e seis grupos de precursores de nêutrons atrasados. A metodologia proposta consiste em expandir os fluxos escalares de nêutrons (rápido e térmico) e as concentrações de precursores de nêutrons atrasados em séries de Taylor para a variável espacial, trazendo a dependência temporal para os coeficientes dessas séries. Truncando as séries de Taylor na ordem quadrática, obtém-se um conjunto de sistemas recursivos de equações diferenciais ordinárias, onde aplica-se o Método da Decomposição modificado, dividindo a matriz dos coeficientes em duas, uma diagonal constante e outra, inserida no termo fonte, com os termos restantes e a dependência temporal. Através desse procedimento, elimina-se o caráter de rigidez das equações e dispensa-se a utilização da continuação analítica. O trabalho tem por objetivo obter uma solução com representação analítica livre de rigidez, com controle de erro, análise de estabilidade e convergência, utilizando a própria equação diferencial como estimativa.