Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Nete trabalho é apresentada uma solução analílica para o problema de ordenada discreta unidimensional e multigrupo de transporle de neutrons em simetria planar. A idéia básica da formulação proposta consiste na aplicação da transformada de Laplace na equação de ordenada discreta. Para a solução do sistema linear resultante, uma solução explícila para a matriz lnversa é estabelecida. Dessa forma, o fluxo angular é obtido, por inversão analítica, em termos do fluxo angular em x=O. Essa formulação é aplicada a problemas de domínio finito e semi-infinito. No primeiro caso, os valores de fluxo angular desconhecidos na fronteira em x=O, são determinados a partir dos valores conhecidos do fluxo angular em x=a; no segundo caso é usada a condição de que o fluxo angular é limilado no infinito. Foram tratados problemas homogêneos e heterogêneos para a placa plana com um grupo de neutrons e multigrupo.O problema inverso, que consiste na determinação do fluxo incidente na fronteira a partir de valores do fluxo escalar no interior do domínio, também foi resolvido. Os resullados obtidos para os problemas acima descritos, apresentaram uma boa comparação com os resultados disponíveis na literatura.

Neste trabalho, estendemos, de forma analítica, a formulação LTSN à problemas de transporte unidimensionais sem simetria azimutal. Para este problema, também apresentamos a solução com dependência contínua na variável angular, a partir da qual é estabelecido um método iterativo de solução da equação de transporte unidimensional. Também discutimos como a formulação LTSN é aplicada na resolução de problemas de transporte unidimensionais dependentes do tempo, tanto de forma aproximada pela inversão numérica do fluxo transformado na variável tempo, bem como analiticamente, pela aplicação do método LTSNnas equações nodais. Simulações numéricas e comparações com resultados disponíveis na literatura são apresentadas. / In this work, we extend, in an analytical way, the LTSNformulation to one-dimensional transport problems without azimuthal symmetry. For this problem, we also present the solution with continuous dependence on the angular variable, from which an iterative method of the one-dimensional transport equation solution is established. We also discuss how the LTSN formulation is applied to the resolution of the one-dimensional time dependent transport problem, both in an approximate way, by the numerical inversion of the transformed time flux, and analytically, by the application of the LTSN method into the nodal equations. Numerical simulations and comparisons with results available in the literature are presented.

O principal objetivo dessa tese consiste em determinar uma solução numéricada equação bidimensional do transporte de nêutrons para elevadas ordens de quadratura angular. Diagonalizando a matriz de transporte LTSN bidimensional , construímos dois algoritmos que se diferenciam pela forma de representar os termos de fuga transversal, que surgem nas equações LTSN integradas transversalmente. Esses termos no método LTSN2D − Diag são expressos como combinação linear dos autovetores multiplicados por exponenciais dos respectivos autovalores. No método LTSN2D − DiagExp os termos de fuga transversal são representados por uma função exponencial com constante de decaimento heuristicamente identificada com parâmetros materiais característicos do meio. A análise epectral desenvolvida permite realizar a diagonalização. Um estudo sobre o condicionamento é feito e também associamos um número de condicionamento ao termo de fuga transversal. Definimos os erros no fluxo aproximado e na fórmula da quadratura, e estabelecemos uma relação entre eles. A convergência ocorre com condições de fronteira e quadratura angular adequadas. Apresentamos os resultados numéricos gerados pelos novos métodos LTSN2D − Diag e LTSN2D − DiagExp para elevadas ordens de quadratura angular para um problema ilustrativo e comparamos com resultados disponíveis na literatura.

Neste trabalho é apresentada uma solução analítica, expressa na forma integral, para a equação SN de transporte unidimensional dependente do tempo em geometria cartesiana, válida para domínios limitado e ilimitado (0 < x < ∞), usando a técnica da dupla transformada de Laplace. A principal ideia consiste na aplicação da transformada de Laplace na variável tempo e na solução da equação resultante pelo método LTSN usando condiçõs de contorno apropriadas para problemas com domínio limitado e ilimitado. Também é apresentada uma nova solução analítica para o conjunto de equações para o fluxo angular, em um retângulo, usando a técnica da transformada dupla de Laplace. Sua principal ideia engloba os passos seguintes: aplicação da transformada de Laplace em uma variável espacial, solução da equação resultante pelo método LTSN e reconstrução do fluxo angular duplamente transformado usando o teorema de inversão da transformada de Laplace. Aqui é assumido que o fluxo angular na fronteira é aproximado por uma função exponencial. Os resultados obti- dos por esses métodos, para os problemas considerados, são comparados com os resultados disponíveis na literatura. / In this work a general analytical solution is reported, expressed in integral form for the time-dependent, one-dimensional SN transport equation in cartesian geometry valid for bounded and unbounded domain (0 < x < ∞), using the double Laplace transform technique. The main idea consists in the application of the Laplace transform technique in time variable and solution of the resulting equation by the LTSN method, using appropriated boundary conditions for bounded and unbounded domain problems. It is also presented a new analytical solution for the set of SN equations for the angular flux, in a retangle, using the double Laplace transform technique. Its main idea comprehends the following steps: application of the Laplace transform in one space variable, solution of the resulting equation by the LTSN method and reconstruction of the double Laplace transformed angular flux using the inversion theorem of the Laplace transform. It is also assumed that the angular flux at boundary is approximated by an exponential function. The results attained by these methods for the problems considered are compared against results available in the literature.

Nesse trabalho é apresentada uma solução analítica para a equação de transporte bidimensional em um domínio retangular considerando o espalhamento Compton, utilizando o método LTSN e a aproximação PN na variável angular. Esse procedimento leva a uma formulação para a taxa de dose total devido µa radiação gama. Nessa derivação, a energia depositada pelo elétron livre é considerada. Para alcançarmos esse objetivo dois problemas são resolvidos: o primeiro para o transporte de fótons, assumindo como seção de choque diferencial o núcleo de espalhamento de Klein-Nishina bem como modelo multigrupo em energia; e o segundo para o transporte de elétrons livres, considerando a seção de choque de espalhamento diferencial de Rutherford. Simulações numéricas e validações com os resultados obtidos pelo método de Monte Carlo são apresentadas para ambos os problemas. / In this work we present the analytical solution for the two-dimensional transport equation in a rectangle considering Compton scattering, using the LTSN method and the PN approach for the angular variable. This procedure leads to a formulation for the total dosis due the gamma radiation. In this derivation, the energy deposited by the free electron is taken into account. To reach this objective we solve two model problems: the ¯rst one for the transport of photons, assuming the Klein-Nishina scattering kernel as the scattering di®erential cross section as well as the multigroup model in the energy variable; and the second one the transport of free electrons considering the screened Rutherford scattering kernel. Numerical simulations and validations with Monte Carlo results are reported for both problems.

No decurso deste trabalho, são apresentadas soluções analíticas para problemas de transporte de nêutrons em geometrias cilíndrica e cartesiana. Para a geometria cilíndrica, usa-se a transformada de Hankel de ordem zero juntamente com o método SN para um problema cilíndrico unidimensional, considerando simetria azimutal e espalhamento isotrópico. Este método é aqui chamado HTSN. O problema cilíndrico com espalhamento anisotrópico é tratado usando o método da decomposição, que possibilita construir um processo recursivo em que a solução HTSN entra como condição uma inicial. Para a geometria cartesiana, a equação de transporte em uma e duas dimensões é derivada em relação à variável angular tantas vezes quantas for a ordem da expansão do espalhamento anisotrópico. Este processo leva a um conjunto de equações íntegro-diferenciais mais a equação puramente diferencial que pode ser resolvida pelo método de separação de variáveis. Seguindo este processo, foi possível chegar às soluções de Case para o problema de transporte em uma dimensão. Também são apresentadas simulações numéricas para um problema de transporte em geometria cilíndrica isotrópico e com anisotropia quadrática. / In this work, we are reported analytical solutions for the transport equation for neutral particles in cylindrical and cartesian geometry. For the cylindrical geometry, it is applied the Hankel transform of order zero in the SN approximation of the one-dimensional cylindrical transport equation, assuming azimuthal symmetry and isotropic scattering. This procedure is coined HTSN method. The anisotropic problem is handled using the decomposition method, generating a recursive approach, which the HTSN solution is used as initial condition. For cartesian geometry, the one and two dimensional transport equation is derived in the angular variable as many time as the degree of the anisotropic scattering. This procedure leads to set of integro-differential plus one differential equation that can be really solved by the variable separation method. Following this procedure, it was possible to come out with the Case solution for the one-dimensional problem. Numerical simulations are reported for the cylindrical transport problem both isotropic and anisotropic case of quadratic degree.

No presente trabalho foi desenvolvido um modelo para o decaimento da energia cinética turbulenta em uma camada limite planetária convectiva. O ponto de partida foi a equação de balanço para a função espectro de energia cinética turbulenta. Para parametrizar o termo de tranferência inercial foi considerado que a energia é transferida de grandes para pequenos turbilhões até a escala de Kolmogorov, onde a energia é dissipada na forma de calor. O termo de fonte ou perda de energia cinética por efeito térmico foi parametrizado considerando que se pode escrevê-lo como o produto de duas funções, onde o decaimento no tempo pode ser expresso por uma função cosseno. O termo de produção de turbulência por efeito mecânico foi parametrizado a partir de análise dimensional e considerando-se a teoria de similaridade de Monin-Obukov. Foi observado que a função espectro de energia decai de forma mais lenta quando se considera a fonte de turbulência mecânica do que quando predomina a turbulência gerada por convecção. A energia cinética turbulenta obtida pelo modelo proposto foi comparada com resultados de LES (Large Eddy Simulation). / In this present work a model for decline of the turbulent kinetic energy in the Convective Boundary Layer was developed. The starting point was the budget equation of the turbulent kinetic energy spectrum. To parametrize the term of inertial transfer, it was considered that the energy is transferred of great to small eddies until the scale of Kolmogorov, where the energy is dissipated in the heat form. In the buoyancy effect term was considered that it can written as the commodity of two functions, where the decline in the time can be express as a cosine function. The turbulence production term generated by shear of the wind (mechanics turbulence) was parametrized from the dimensional analysis and considering it the theory of the similarity of Monin-Obukhov. It was observed that the energy spectrum function to decays slower when is considered the mechanic source of what when predominates the turbulence generated by convection. The turbulent kinetic energy gotten by considered model was compared with results of LES (Large Eddy Simulation).

O trabalho proposto apresenta um novo método analítico para a resolução da equação de Schrödinger dependente do tempo. Esse método é empregado na simulação de reações químicas envolvendo compostos orgânicos, a fim de estimar o impacto ambiental decorrente do lançamento de substâncias químicas em corpos hídricos. O método é baseado no emprego de transformações de Bäcklund, que fornecem soluções em forma fechada expressas como funções arbitrárias de argumentos conhecidos. Essa formulação proporciona uma redução significativa no tempo de processamento requerido para a obtenção de resultados numéricos, permitindo a simulação de sistemas químicos contendo até 200 átomos. Simulações simbólicas são apresentadas. / This work presents a new analytical method for solving the time dependent Schrödinger equation. This method is employed to simulate chemical reactions involving organic compounds, in order to estimate environmental damages due to the disposal of chemicals in water bodies. The method is based on Bäcklund transformations, which supply closed form solutions as arbitrary functions of known arguments. This formulation provides a significant reduction on the time processing required to obtain numerical results, allowing to carry out the simulation of chemical systems containing up to 200 atoms. Symbolic simulations are reported.

Neste trabalho é apresentada uma solução analítica para a aproximação da equação de transporte linear unidimennal em geometria plana, considerando modelo de multigrupo e espalhamento anisotrópico. A idéia principal desse método consiste em aplicar a transformada de Laplace ao sistema de equações diferenciais ordinárias Este procedimento gera um Sistema linear para o fluxo angular transformado. Resolvendo esse sistema pelo algoritmo de Trzaska, o fluxo angular é obtido em termos do fluxo angular na fronteira x = O pela técnica de inversão de Heaviside. Os resultados obtidos por este método para problemas de placa plana, homogênea e heterogênea, para um e dois grupos de energia, em dom1nio finito e semi-inf1nito, bem como os problemas inversos determinação do parâmetro c, da espessura critica de uma placa e do fluxo angular incidente na fronteira de uma placa plana, para um e dois grupos de energia, foram comparados com os resultados disponíveis na literatura e apresentaram boa concordância. / In this work is presented an approximated analytical solution for the one-dimensional slab-geometry linear transport equation by considering multigroup model and an1sotropic scattering. The main idea of this approach is based on the application of the Laplace transform into the set of the PN ordlnary d1fferential equations. This procedure leads to a linear system to be solved for the transformed angular flux by the Trzaska's algorithm. Once this system is solved, the angular flux is then obtained as a function of the angular flux at the boundary x = O by using the Heaviside's expansion tecn nique. The results achieved by th1s method for the homogenecus and heterogeneous slab-geometry problems in a finite and semi-infinite domain, considering the multigroup model and anisotropic scattering as well for the inverse problems: determination of the c parameter, criticai thickness of a slab and the incoming angular flux at the boundary, were compared with the ones available in the literature showing a very good agreement.