A novel collision source method has been developed to solve the Linear Boltzmann Equation (LBE) more efficiently by adaptation of the angular quadrature order. The angular adaptation method is unique in that the flux from each scattering source iteration is obtained, with potentially a different quadrature order used for each. Traditionally, the flux from every iteration is combined, with the same quadrature applied to the combined flux. Since the scattering process tends to distribute the radiation more evenly over angles (i.e., make it more isotropic), the quadrature requirements generally decrease with each iteration. This method allows for an optimal use of processing power, by using a high order quadrature for the first few iterations that need it, before shifting to lower order quadratures for the remaining iterations. This is essentially an extension of the first collision source method, and is referred to as the adaptive collision source (ACS) method. The ACS methodology has been implemented in the 3-D, parallel, multigroup discrete ordinates code TITAN. This code was tested on a variety of test problems including fixed-source and eigenvalue problems. The ACS implementation in TITAN has shown a reduction in computation time by a factor of 1.5-4 on the fixed-source test problems, for the same desired level of accuracy, as compared to the standard TITAN code. / Ph. D.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto. / In this masters dissertation we describe a hybrid direct method for calculating the disadvantage factor and the neutron flux distribution in fuel-moderator lattices. For the mathematical model, we used the discrete ordinates (SN) transport equation, considering linearly anisotropic scattering in the monoenergetic model and isotropic scattering in the energy multigroup model in slab geometry. We describe a spectral analysis of the monoenergetic and two-group SN equations, in a way which is very similar to the Case method. The basic idea is to use higher order angular quadrature set in the highly absorbing fuel region (SNF)and lower order angular quadrature set in the diffusive moderator region (SNM) i.e., NF > NM. Therefore, we apply special continuity conditions for the fuel existing fluxes that constitute the incoming fluxes for the moderator region, and conversely for the moderator existing fluxes that constitute the incoming fluxes for the fuel region, based on the equivalence of the SN and PN-1 equations, which characterize the hybrid model. As a direct hybrid method, we use NF + NM linear algebraic equations composed of (NF + NM)/2 reflexive boundary conditions and (NF + NM)/2 continuity conditions to solve for the NF + NM expansion coefficients. With these coefficients we can calculate the numerical values for the angular fluxes and for the scalar fluxes at any location of domain. We present numerical results to illustrate the efficiency and the accuracy of the offered method.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto. / In this masters dissertation we describe a hybrid direct method for calculating the disadvantage factor and the neutron flux distribution in fuel-moderator lattices. For the mathematical model, we used the discrete ordinates (SN) transport equation, considering linearly anisotropic scattering in the monoenergetic model and isotropic scattering in the energy multigroup model in slab geometry. We describe a spectral analysis of the monoenergetic and two-group SN equations, in a way which is very similar to the Case method. The basic idea is to use higher order angular quadrature set in the highly absorbing fuel region (SNF)and lower order angular quadrature set in the diffusive moderator region (SNM) i.e., NF > NM. Therefore, we apply special continuity conditions for the fuel existing fluxes that constitute the incoming fluxes for the moderator region, and conversely for the moderator existing fluxes that constitute the incoming fluxes for the fuel region, based on the equivalence of the SN and PN-1 equations, which characterize the hybrid model. As a direct hybrid method, we use NF + NM linear algebraic equations composed of (NF + NM)/2 reflexive boundary conditions and (NF + NM)/2 continuity conditions to solve for the NF + NM expansion coefficients. With these coefficients we can calculate the numerical values for the angular fluxes and for the scalar fluxes at any location of domain. We present numerical results to illustrate the efficiency and the accuracy of the offered method.