Um método sintético de difusão para aceleração do esquema de fonte de espalhamento em cálculos SN unidimensionais de fonte fixa

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.

Um problema inverso na modelagem da difusão do calor

O presente trabalho aborda um problema inverso associado a difus~ao de calor em uma barra unidimensional. Esse fen^omeno e modelado por meio da equac~ao diferencial par- cial parabolica ut = uxx, conhecida como equac~ao de difus~ao do calor. O problema classico (problema direto) envolve essa equac~ao e um conjunto de restric~oes { as condic~oes inicial e de contorno {, o que permite garantir a exist^encia de uma soluc~ao unica. No problema inverso que estudamos, o valor da temperatura em um dos extremos da barra n~ao esta disponvel. Entretanto, conhecemos o valor da temperatura em um ponto x0 xo no interior da barra. Para aproximar o valor da temperatura no intervalo a direita de x0, propomos e testamos tr^es algoritmos de diferencas nitas: diferencas regressivas, leap-frog e diferencas regressivas maquiadas.

Modelagem computacional da estabilização de sistemas subcríticos segundo o modelo unidimensional de difusão de nêutrons monoenergéticos

Um método numérico espectronodal (END) livre de erros de truncamento espacial é desenvolvido para problemas unidimensionais de difusão de nêutrons monoenergéticos em duas versões. Na versão de problemas de autovalor, o método gera soluções numéricas para o perfil do fluxo escalar e para o fator de multiplicação efetivo (k), que coincidem com a solução analítica dominante, afora os erros da aritmética finita computacional. Na versão de fonte fxa, o método também gera soluções numéricas analíticas para o problema de fonte fixa correspondente, onde a fonte de fissão, com dependência espacial, é obtida analiticamente, a partir da reconstrução espacial do fluxo escalar gerado pelo método END para problemas de autovalor. Alguns experimentos numéricos são apresentados para dois problemas modelos a fim de ilustrar a precisão do método.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Problemas inversos em processos difusivos com retenção

Um Estudo para a solução numérica do modelo de difusão com retenção, proposta por Bevilacqua et al. (2011), é apresentado, bem como uma formulação implícita para o problema inverso para a estimativa dos parâmetros envolvidos na formulação matemática do modelo. Através de um estudo minucioso da análise de sensibilidade e do cálculo do coeficiente de correlação de Pearson, são identificadas as chances de se obter sucesso na solução do problema inverso através do método determinístico de Levenberg-Marquardt e dos métodos estocásticos Algoritmo de Colisão de Partículas (Particle Collision Algorithm - PCA) e Evolução Diferencial (Differential Evolution - DE). São apresentados os resultados obtidos através destes três métodos de otimização para três casos de conjunto de parâmetros. Foi observada uma forte correlação entre dois destes três parâmetros, o que dificultou a estimativa simultânea dos mesmos. Porém, foi obtido sucesso nas estimativas individuais de cada parâmetro. Foram obtidos bons resultados para os fatores que multiplicam os termos diferenciais da equação que modela o fenômeno de difusão com retenção.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.