Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O esquema iterativo de fonte de espalhamento (SI) é tradicionalmente aplicado para a convergência da solução numérica de malha fina para problemas de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas com fonte fixa. O esquema SI é muito simples de se implementar sob o ponto de vista computacional; porém, o esquema SI pode apresentar taxa de convergência muito lenta, principalmente para meios difusivos (baixa absorção) com vários livres caminhos médios de extensão. Nesta dissertação descrevemos uma técnica de aceleração baseada na melhoria da estimativa inicial para a distribuição da fonte de espalhamento no interior do domínio de solução. Em outras palavras, usamos como estimativa inicial para o fluxo escalar médio na grade de discretização de malha fina, presentes nos termos da fonte de espalhamento das equações discretizadas SN usadas nas varreduras de transporte, a solução numérica da equação da difusão de nêutrons em grade espacial de malha grossa com condições de contorno especiais, que aproximam as condições de contorno prescritas que são clássicas em cálculos SN, incluindo condições de contorno do tipo vácuo. Para aplicarmos esta solução gerada pela equação da difusão em grade de discretização de malha grossa nas equações discretizadas SN de transporte na grade de discretização de malha fina, primeiro implementamos uma reconstrução espacial dentro de cada nodo de discretização, e então determinamos o fluxo escalar médio em grade de discretização de malha fina para usá-lo nos termos da fonte de espalhamento. Consideramos um número de experimentos numéricos para ilustrar a eficiência oferecida pela presente técnica (DSA) de aceleração sintética de difusão. / The scattering source iterative (SI) scheme is traditionally applied to converge finemesh numerical solutions to fixed-source discrete ordinates neutron transport problems. The SI scheme is very simple to implement under a computational viewpoint. However, the SI scheme may show very slow convergence rate, mainly for diffusive media (low absorption) with several mean free paths in extent. In this work we describe an acceleration technique based on an improved initial guess for the scattering source distribution within the slab. In other words, we use as initial guess for the fine-mesh average scalar flux in the scattering source terms of the SN discretized equations used in the transport sweeps, the coarse-mesh solution of the neutron diffusion equation with special boundary conditions to account for the classical SN prescribed boundary conditions, including vacuum boundary conditions. To apply this coarse-mesh diffusion solution into the fine-mesh SN transport sweep discretized equations, we first perform within-node spatial reconstruction, and then we determine the fine-mesh average scalar flux for use in the scattering source terms. We consider a number of numerical experiments to illustrate the efficiency of the offered diffusion synthetic acceleration (DSA) technique.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto. / A numerical nodal method that is free from all spatial truncation errors is developed for one-speed slab-geometry discrete ordinates (SN) fixed-source adjoint neutral particle transport problems. The unknown in the method are the node-edge and the node-average adjoint angular fluxes, and the numerical values obtained for these quantities are those of the analytic solution of the adjoint SN equations. The method is based on the use of the standard spatially discretized SN balance adjoint equation, which holds in each spatial node and for each discrete ordinates direction, and a nonstandard adjoint auxiliary equation that contains a Greens function for the node-average adjoint angular fluxes in terms of the exiting adjoint angular fluxes from the node edges and the adjoint interior source. Numerical results are given to illustrate the methods accuracy.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Um método espectronodal é desenvolvido para problemas de transporte de partículas neutras de fonte fixa, multigrupo de energia em geometria cartesiana na formulação de ordenadas discretas (SN). Para geometria unidimensional o método espectronodal multigrupo denomina-se método spectral Greens function (SGF) com o esquema de inversão nodal (NBI) que converge solução numérica para problemas SN multigrupo em geometria unidimensional, que são completamente livre de erros de truncamento espacial para ordem L de anisotropia de espalhamento desde que L < N. Para geometria X; Y o método espectronodal multigrupo baseia-se em integrações transversais das equações SN no interior dos nodos de discretização espacial, separadamente nas direções coordenadas x e y. Já que os termos de fuga transversal são aproximados por constantes, o método nodal resultante denomina-se SGF-constant nodal (SGF-CN), que é aplicado a problemas SN multigrupo de fonte fixa em geometria X; Y com espalhamento isotrópico. Resultados numéricos são apresentados para ilustrar a eficiência dos códigos SGF e SGF-CN e a precisão das soluções numéricas convergidas em cálculos de malha grossa. / A spectral nodal method is described for neutral particle energy multigroup fixed-source transport problems in cartesian geometry in the discrete ordinates (SN) formulation. For slab geometry the offered multigroup spectral nodal method is referred to as the spectral Greens function (SGF) method with the one-node block inversion (NBI) iterative scheme, which converges numerical solutions to multigroup slab-geometry SN problems, that are completely free from spatial truncation errors for scattering anisotropy of order L, provided L < N. For X; Y-geometry, the offered multigroup spectral nodal method is based on transverse integrations of the SN equations inside the discretization nodes, separately in x- and y- coordinate directions. Since the transverse-leakage terms are approximated by constants, the resulting nodal method is referred to as the multigroup SGF-contant nodal (SGF-CN) method, which is applied for multigroup X; Y-geometry fixed-source SN problems with isotropic scattering. Numerical results are presented to illustrate the efficiency of the SGF and SGF-CN codes and the accuracy of the converged numerical solutions in coarse-mesh calculations.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O esquema iterativo de fonte de espalhamento (SI) é tradicionalmente aplicado para a convergência da solução numérica de malha fina para problemas de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas com fonte fixa. O esquema SI é muito simples de se implementar sob o ponto de vista computacional; porém, o esquema SI pode apresentar taxa de convergência muito lenta, principalmente para meios difusivos (baixa absorção) com vários livres caminhos médios de extensão. Nesta dissertação descrevemos uma técnica de aceleração baseada na melhoria da estimativa inicial para a distribuição da fonte de espalhamento no interior do domínio de solução. Em outras palavras, usamos como estimativa inicial para o fluxo escalar médio na grade de discretização de malha fina, presentes nos termos da fonte de espalhamento das equações discretizadas SN usadas nas varreduras de transporte, a solução numérica da equação da difusão de nêutrons em grade espacial de malha grossa com condições de contorno especiais, que aproximam as condições de contorno prescritas que são clássicas em cálculos SN, incluindo condições de contorno do tipo vácuo. Para aplicarmos esta solução gerada pela equação da difusão em grade de discretização de malha grossa nas equações discretizadas SN de transporte na grade de discretização de malha fina, primeiro implementamos uma reconstrução espacial dentro de cada nodo de discretização, e então determinamos o fluxo escalar médio em grade de discretização de malha fina para usá-lo nos termos da fonte de espalhamento. Consideramos um número de experimentos numéricos para ilustrar a eficiência oferecida pela presente técnica (DSA) de aceleração sintética de difusão. / The scattering source iterative (SI) scheme is traditionally applied to converge finemesh numerical solutions to fixed-source discrete ordinates neutron transport problems. The SI scheme is very simple to implement under a computational viewpoint. However, the SI scheme may show very slow convergence rate, mainly for diffusive media (low absorption) with several mean free paths in extent. In this work we describe an acceleration technique based on an improved initial guess for the scattering source distribution within the slab. In other words, we use as initial guess for the fine-mesh average scalar flux in the scattering source terms of the SN discretized equations used in the transport sweeps, the coarse-mesh solution of the neutron diffusion equation with special boundary conditions to account for the classical SN prescribed boundary conditions, including vacuum boundary conditions. To apply this coarse-mesh diffusion solution into the fine-mesh SN transport sweep discretized equations, we first perform within-node spatial reconstruction, and then we determine the fine-mesh average scalar flux for use in the scattering source terms. We consider a number of numerical experiments to illustrate the efficiency of the offered diffusion synthetic acceleration (DSA) technique.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto. / A numerical nodal method that is free from all spatial truncation errors is developed for one-speed slab-geometry discrete ordinates (SN) fixed-source adjoint neutral particle transport problems. The unknown in the method are the node-edge and the node-average adjoint angular fluxes, and the numerical values obtained for these quantities are those of the analytic solution of the adjoint SN equations. The method is based on the use of the standard spatially discretized SN balance adjoint equation, which holds in each spatial node and for each discrete ordinates direction, and a nonstandard adjoint auxiliary equation that contains a Greens function for the node-average adjoint angular fluxes in terms of the exiting adjoint angular fluxes from the node edges and the adjoint interior source. Numerical results are given to illustrate the methods accuracy.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Um método espectronodal é desenvolvido para problemas de transporte de partículas neutras de fonte fixa, multigrupo de energia em geometria cartesiana na formulação de ordenadas discretas (SN). Para geometria unidimensional o método espectronodal multigrupo denomina-se método spectral Greens function (SGF) com o esquema de inversão nodal (NBI) que converge solução numérica para problemas SN multigrupo em geometria unidimensional, que são completamente livre de erros de truncamento espacial para ordem L de anisotropia de espalhamento desde que L < N. Para geometria X; Y o método espectronodal multigrupo baseia-se em integrações transversais das equações SN no interior dos nodos de discretização espacial, separadamente nas direções coordenadas x e y. Já que os termos de fuga transversal são aproximados por constantes, o método nodal resultante denomina-se SGF-constant nodal (SGF-CN), que é aplicado a problemas SN multigrupo de fonte fixa em geometria X; Y com espalhamento isotrópico. Resultados numéricos são apresentados para ilustrar a eficiência dos códigos SGF e SGF-CN e a precisão das soluções numéricas convergidas em cálculos de malha grossa. / A spectral nodal method is described for neutral particle energy multigroup fixed-source transport problems in cartesian geometry in the discrete ordinates (SN) formulation. For slab geometry the offered multigroup spectral nodal method is referred to as the spectral Greens function (SGF) method with the one-node block inversion (NBI) iterative scheme, which converges numerical solutions to multigroup slab-geometry SN problems, that are completely free from spatial truncation errors for scattering anisotropy of order L, provided L < N. For X; Y-geometry, the offered multigroup spectral nodal method is based on transverse integrations of the SN equations inside the discretization nodes, separately in x- and y- coordinate directions. Since the transverse-leakage terms are approximated by constants, the resulting nodal method is referred to as the multigroup SGF-contant nodal (SGF-CN) method, which is applied for multigroup X; Y-geometry fixed-source SN problems with isotropic scattering. Numerical results are presented to illustrate the efficiency of the SGF and SGF-CN codes and the accuracy of the converged numerical solutions in coarse-mesh calculations.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência. / As neutron fission events do not take place in the non-multiplying regions of nuclear reactors, e.g., moderator, reflector, and structural core, these regions do not generate power and the computational efficiency of nuclear reactor global calculations can hence be improved by eliminating the explicit numerical calculations within the non-multiplying regions around the active domain. Discussed here is the computational efficiency of approximate discrete ordinates (SN) albedo boundary conditions for two-energy group eigenvalue problems in X,Y geometry. Albedo, the Latin word for whiteness, was originally defined as the fraction of incident light reflected diffusely by a surface. This Latin word has remained the usual scientific term in astronomy and in this dissertation this concept is extended for the reflection of neutrons. The non-standard SN albedo substitutes approximately the reflector region around the active domain, as we neglect the transverse leakage terms within the nonmultiplying reflector. Should the problem have no transverse leakage terms, i.e., onedimensional slab geometry, then the offered albedo boundary conditions are exact. By computational efficiency we mean analyzing the accuracy of the numerical results versus the CPU execution time of each run for a given model problem. Numerical results to two symmetric test problems are shown to illustrate this efficiency analysis.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência. / As neutron fission events do not take place in the non-multiplying regions of nuclear reactors, e.g., moderator, reflector, and structural core, these regions do not generate power and the computational efficiency of nuclear reactor global calculations can hence be improved by eliminating the explicit numerical calculations within the non-multiplying regions around the active domain. Discussed here is the computational efficiency of approximate discrete ordinates (SN) albedo boundary conditions for two-energy group eigenvalue problems in X,Y geometry. Albedo, the Latin word for whiteness, was originally defined as the fraction of incident light reflected diffusely by a surface. This Latin word has remained the usual scientific term in astronomy and in this dissertation this concept is extended for the reflection of neutrons. The non-standard SN albedo substitutes approximately the reflector region around the active domain, as we neglect the transverse leakage terms within the nonmultiplying reflector. Should the problem have no transverse leakage terms, i.e., onedimensional slab geometry, then the offered albedo boundary conditions are exact. By computational efficiency we mean analyzing the accuracy of the numerical results versus the CPU execution time of each run for a given model problem. Numerical results to two symmetric test problems are shown to illustrate this efficiency analysis.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese. / In this work, three techniques for numerically solving one-speed neutral particle inverse transport problems for nuclear engineering applications are developed. It is well known that direct steady-state monoenergetic transport problems are characterized by estimating the flux of particles as a distribution function of space and direction-of-motion independent variables, when the material parameters (cross sections), the geometry, and the incoming flux at the boundaries of the domain (boundary conditions), as well as the interior source distribution are known. Conversely, inverse problems, in this work, seek for estimates to the incident boundary flux, or interior source, or void fractions in homogeneous slabs. The mathematical model used for direct and inverse problems is the time-independent one-speed slab-geometry transport equation with linearly anisotropic scattering in the discrete ordinates (SN) formulation. In the boundary-value inverse problems, given the existing flux at one boundary of the slab, as measured by a neutron detector, for example, we seek for accurate estimate for the incident flux at the opposite boundary. On the other hand, in the interior source inverse SN problems, we seek for accurate estimate for the interior source stored within the slab for shielding purpose, given the exiting flux at the boundary of the slab. Furthermore, as with the void fraction inverse SN problems, given the exiting flux at one boundary of the slab due to prescribed incident flux at the opposite boundary, we seek for accurate estimate of the void fraction within the slab in the context of non-destructive testing applications in industry. The computer code developed in this work presents the coarse-mesh spectral Greens function (SGF) nodal method for direct SN problems in slab geometry to generate accurate numerical solutions to the three inverse SN problems described above. For the boundary-value and interior source inverse SN problems, we use the proportionality property of the leakage of particles; moreover, for the void fraction inverse SN problems, we offer the technique that we refer to as the physical bisection method. We present numerical results to illustrate the accuracy of the three techniques, as described in this dissertation.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese. / In this work, three techniques for numerically solving one-speed neutral particle inverse transport problems for nuclear engineering applications are developed. It is well known that direct steady-state monoenergetic transport problems are characterized by estimating the flux of particles as a distribution function of space and direction-of-motion independent variables, when the material parameters (cross sections), the geometry, and the incoming flux at the boundaries of the domain (boundary conditions), as well as the interior source distribution are known. Conversely, inverse problems, in this work, seek for estimates to the incident boundary flux, or interior source, or void fractions in homogeneous slabs. The mathematical model used for direct and inverse problems is the time-independent one-speed slab-geometry transport equation with linearly anisotropic scattering in the discrete ordinates (SN) formulation. In the boundary-value inverse problems, given the existing flux at one boundary of the slab, as measured by a neutron detector, for example, we seek for accurate estimate for the incident flux at the opposite boundary. On the other hand, in the interior source inverse SN problems, we seek for accurate estimate for the interior source stored within the slab for shielding purpose, given the exiting flux at the boundary of the slab. Furthermore, as with the void fraction inverse SN problems, given the exiting flux at one boundary of the slab due to prescribed incident flux at the opposite boundary, we seek for accurate estimate of the void fraction within the slab in the context of non-destructive testing applications in industry. The computer code developed in this work presents the coarse-mesh spectral Greens function (SGF) nodal method for direct SN problems in slab geometry to generate accurate numerical solutions to the three inverse SN problems described above. For the boundary-value and interior source inverse SN problems, we use the proportionality property of the leakage of particles; moreover, for the void fraction inverse SN problems, we offer the technique that we refer to as the physical bisection method. We present numerical results to illustrate the accuracy of the three techniques, as described in this dissertation.