Um estudo sobre equações auxiliares em formulações nodais explícitas para solução da equação de transporte de nêutrons bidimensional

Abstract

Neste trabalho, uma abordagem analítica é utilizada juntamente com esquemas nodais para resolução de um problema de transporte de nêutrons de fonte fixa, definido em meios homogêneos e heterogêneos, em geometria cartesiana bidimensional, com espalhamento isotrópico. A metodologia é desenvolvida a partir da versão em ordenadas discretas da equação bidimensional de transporte. Utiliza-se esquemas nodais para obtenção de equações unidimensionais integradas transversalmente, as quais são resolvidas via método ADO (Analytical Discrete Ordinates). Expressões explícitas em termos das variáveis espaciais são determinadas para os fluxos angulares médios em cada região em que o domínio foi subdividido. Para obtenção da solução em todo o domínio, acopla-se as soluções de cada região às regiões vizinhas, através de um sistema linear. Neste contexto, o objetivo principal é o estudo relacionado às equações auxiliares necessárias para aproximação dos termos desconhecidos nos contornos do domínio ou das interfaces, que surgem devido ao processo de integração. Três abordagens distintas para aproximação dos termos de fuga transversais, oriundos da integração transversal, são estudadas: aproximações por constantes, aproximações lineares e aproximações exponenciais, as quais são incorporadas ao termo fonte. Adicionalmente, neste trabalho quatro esquemas de quadraturas numéricas são utilizados para aproximar o termo integral de espalhamento: quadratura Simétrica de Nível (LQN), Legendre-Chebyshev quadrangular (PNTN), Legendre-Chebyshev triangular (PNTNSN) e Quadruple Range (QR). Resultados numéricos são obtidos para o fluxo escalar médio em regiões do domínio e comparados com resultados disponíveis na literatura bem como gerados pelo código AHOT. A análise dos resultados confirma a viabilidade da proposta das equações auxiliares alternativas, mantendo a eficiência computacional já verificada em outras abordagens do método ADO, no entanto indica que estudos complementares necessitam ser realizados para caracterizar vantagens adicionais no uso de tais propostas.In this work, an analytical approach is used along with nodal schemes for solving xed source neutron transport problems de ned in two-dimensional homogeneous and heterogeneous medium with isotropic scattering. The methodology is developed from the discrete ordinates version of the two-dimensional transport equation. Nodal procedures are performed to derive one-dimensional transverse integrated equations, which are solved by the ADO method. Explicit expressions in terms of the spatial variables are obtained for averaged angular uxes in each region in which the domain is subdivided. The solution for the whole domain is obtained trough the coupling of the local solutions in a general linear system. In this context, the main goal is to analyse the use of different auxiliary equations to describe the unknown transverse leakage terms on the boundaries and interfaces. Three di erent approaches are used to approximate the unknown transverse leakage terms: constant, linear and exponential approximations, which are incorporated into the source term. Four numerical quadrature schemes are used to approximate the integral scattering term: Level Symmetric quadrature scheme LQN, Legendre-Chebyshev quadrangular PNTN, Legendre-Chebyshev triangular PNTNSN and Quadruple Range (QR). The numerical results obtained for region-averaged scalar uxes are compared with results available in the literature as well as numerical results provided by the AHOT code. The analysis of the results con rm the feasibility with computational e ciency as usual for the ADO method. However subsequent work is indicated to evidence advantages of the proposed auxiliary equations