Solução LTSN da equação de transporte em geometria cartesiana unidimensional para c=1

Abstract

Nos últimos anos, o método LTSN - que resolve diversos problemas de transporte em uma placa plana - emergiu de forma contínua em nossa literatura. Porém, não é de nosso conhecimento que se tenha aplicado este método a problemas isotrópicos de transporte de partículas neutras em uma placa plana, quando o parâmetro albedo vale c = 1. Sabemos que para esta situação a equação de transporte unidimensional apresenta dois autovalores que se encontram no infinito. Consequentemente, a formulação LTSN não pode ser aplicada, pois a solução LTSN é utilizada para problemas onde a matriz LTSN é diagonalizável, e isto ocorre quando c 6= 1. Para a resolução destes tipos de problemas, nós modificamos a solução LTSN que aproxima a solução de Case quando o albedo é unitário, combinando de forma adequada a decomposição de Schur e a expansão de Heavside. A convergência provada do Método LTSN permite que determinemos a solução com precisão prescrita. Apresentamos simulações e comparações numéricas com resultados disponíveis na literatura. Por este procedimento, esperamos terminar o estudo da praticabilidade do LTSN para resolver problemas do transporte em uma placa plana.

In the last years, the LTSN method - for transport problems in slabgeometry appears frequently in the literature. However, to our knowledge, this method has not been applied to the solution of neutral particle transport problems in a slab with isotropic scattering for c = 1. We know that in this situation the neutron transport equation presents two eigenvalues that coalesce to infinity. Therefore, the LTSN formulation can not be applied to this type of problem, because the LTSN solution is derived for problems in which the LTSN matrix is non-defective, that is for problems with c 6= 1. To solve these types of problems we modify the LTSN solution that approaches the Case solution when the albedo is unitary, combining in an adequate form the Schur decomposition and the Heavside expansion. The convergence of the LTSN method allows us to determine the solution with prescribed accuracy. We present numerical simulations and comparisons with results available in the literature. By this procedure we hope to complete the study of the LTSN formulation to solve transport problems in slab-geometry.