Avanços no método LTSn para cálculo de criticalidade e desenvolvimento da primeira versão do código LTSn

O objetivo deste trabalho consiste no desenvolvimento de alguns avanços, teóricos e numéricos, no método LTSN visando implementar a primeira versão de um código computacional, para resolver a equação de transporte utilizando formulação LTSN na forma de multigrupos em geometria plana. Os avanços para o método LTSN estão fundamentados na solução iterativa para o sistema de equações que constituem as condições de contorno, um novo método para a busca do valor de keff baseado no método da bissecção. O desenvolvimento desta metodologia permitiu realizar o calculo muito rápido com altas ordens de quadratura e com esforço computacional muito reduzido. Juntos os avanços matemáticos e numéricos, implementados nesta primeira versão de um código fortran, tal como nos códigos já conhecidos permite solucionar a equação de transporte na forma de multigrupos, tanto para o cálculo direto como para o adjunto, com fontes arbitrárias. Este código utiliza de recursos computacionais da linguagem FORTRAN e as bibliotecas LAPACK, para a otimização de seus algoritmos, facilitando o desenvolvimento futuro. A validação deste trabalho foi feita utilizando dois problemas: um relativo ao fluxo angular e escalar, tanto para o fluxo direto como para o adjunto, cuja importância está relacionada com busca de convergência, relação de reciprocidade e comprovação da solução adjunta, e; um problema de criticalidade, para comprovar a eficácia do algoritmo de busca iterativa de keff e espessura crítica. Com este trabalho se abrem muitas possibilidades tanto teóricas como numéricas a investigar com o método LTSN.

Determinação de parâmetros de criticalidade, fluxo e potência em placas planas pelo método LTSn

O objetivo deste trabalho consiste em aplicar o método LTSn em cálculos de parâmetros críticos como Keff, espessura e concentração atômica e obtenção do fiuxo escalar, da potência específica e do enriquecimento do combustível em placa plana homogenea e heterogênea, considerando modelo multigrupo e em diversas ordens de quadraturas. O método LTSn consiste na aplicação da transformada de Laplace em um conjunto de equações~de ordenadas discretas gerado pela aproximação SN, resultando em um sistema de equações algébricas simbólicas dependentes do parâmetro complexo s e reconstrução dos fluxos angulares pela técnica de expansão de Heaviside. A aplicação do método LTSn reduz a soluçào de um problema de autovalor, a solução de uma equação transcedental, possibilitando a obtenção de parâmetros críticos. Simulações numéricas são apresentadas.

Fases e criticalidade no modelo ashkin - teller de três cores

The usual Ashkin-Teller (AT) model is obtained as a superposition of two Ising models coupled through a four-spin interaction term. In two dimension the AT model displays a line of fixed points along which the exponents vary continuously. On this line the model becomes soluble via a mapping onto the Baxter model. Such richness of multicritical behavior led Grest and Widom to introduce the N-color Ashkin-Teller model (N-AT). Those authors made an extensive analysis of the model thus introduced both in the isotropic as well as in the anisotropic cases by several analytical and computational methods. In the present work we define a more general version of the 3-color Ashkin-Teller model by introducing a 6-spin interaction term. We investigate the corresponding symmetry structure presented by our model in conjunction with an analysis of possible phase diagrams obtained by real space renormalization group techniques. The phase diagram are obtained at finite temperature in the region where the ferromagnetic behavior is predominant. Through the use of the transmissivities concepts we obtain the recursion relations in some periodical as well as aperiodic hierarchical lattices. In a first analysis we initially consider the two-color Ashkin-Teller model in order to obtain some results with could be used as a guide to our main purpose. In the anisotropic case the model was previously studied on the Wheatstone bridge by Claudionor Bezerra in his Master Degree dissertation. By using more appropriated computational resources we obtained isomorphic critical surfaces described in Bezerra's work but not properly identified. Besides, we also analyzed the isotropic version in an aperiodic hierarchical lattice, and we showed how the geometric fluctuations are affected by such aperiodicity and its consequences in the corresponding critical behavior. Those analysis were carried out by the use of appropriated definitions of transmissivities. Finally, we considered the modified 3-AT model with a 6-spin couplings. With the inclusion of such term the model becomes more attractive from the symmetry point of view. For some hierarchical lattices we derived general recursion relations in the anisotropic version of the model (3-AAT), from which case we can obtain the corresponding equations for the isotropic version (3-IAT). The 3-IAT was studied extensively in the whole region where the ferromagnetic couplings are dominant. The fixed points and the respective critical exponents were determined. By analyzing the attraction basins of such fixed points we were able to find the three-parameter phase diagram (temperature £ 4-spin coupling £ 6-spin coupling). We could identify fixed points corresponding to the universality class of Ising and 4- and 8-state Potts model. We also obtained a fixed point which seems to be a sort of reminiscence of a 6-state Potts fixed point as well as a possible indication of the existence of a Baxter line. Some unstable fixed points which do not belong to any aforementioned q-state Potts universality class was also found

Estudo sobre a detecção de invariância de escala discreta em sistemas com criticalidade auto-organizada

Recent studies have shown evidence of log-periodic behavior in non-hierarchical systems. An interesting fact is the emergence of such properties on rupture and breakdown of complex materials and financial failures. These may be examples of systems with self-organized criticality (SOC). In this work we study the detection of discrete scale invariance or log-periodicity. Theoretically showing the effectiveness of methods based on the Fourier Transform of the log-periodicity detection not only with prior knowledge of the critical point before this point as well. Specifically, we studied the Brazilian financial market with the objective of detecting discrete scale invariance in Bovespa (Bolsa de Valores de S˜ao Paulo) index. Some historical series were selected periods in 1999, 2001 and 2008. We report evidence for the detection of possible log-periodicity before breakage, shown its applicability to the study of systems with discrete scale invariance likely in the case of financial crashes, it shows an additional evidence of the possibility of forecasting breakage

Condições de contorno albedo para cálculos globais de reatores nucleares térmicos com o modelo de ordenadas discretas a dois grupos de energia

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.

Modelo de Ginzburg-Landau a partir da teoria de campos a temperatura finita

Neste trabalho, utilizamos o formalismo de teorias quânticas de campos a temperatura finita, tal como desenvolvidas por Matsubara, aplicado a uma hamiltoniana de N campos escalares com autointeração quártica a N grande. Obtém-se uma expressão, na primeira aproximação quântica, para o coeficiente do termo quadrático da hamiltoniana ("massa quadrada"), renormalizado, como função da temperatura. A partir dela, estudamos o processo de quebra espontânea de simetria. Por outro lado, a mesma hamiltoniana é conhecida como modelo de Ginzburg-Landau na literatura de matéria condensada, e que permite o estudo de transições de fase em materiais ferromagnéticos. A temperatura é introduzida através do termo quadrático na hamiltoniana, de forma linear: é proporcional à diferença entre a variável de temperatura e a temperatura crítica. Tal modelo, porém, possui validade apenas na regi~ao de temperaturas próximas à criticalidade. Como resultado de nossos cálculos na teoria de campos a temperatura finita, observamos que, numa faixa de valores em torno da temperatura crítica, a massa quadrática pode ser aproximada por uma relação linear em relação à variável de temperatura. Isso evidencia a compatibilidade da abordagem de Ginzburg-Landau, na vizinhança da criticalidade, com respeito ao formalismo de campos a temperatura finita. Discutimos também os efeitos causados pela presença de um potencial químico no sistema.

Condições de contorno tipo albedo para cálculos globais de reatores nucleares com o modelo de dois grupos de energia na formulação de transporte SN em geometria bidimensional cartesiana

Os eventos de fissão nuclear, resultados da interação dos nêutrons com os núcleos dos átomos do meio hospedeiro multiplicativo, não estão presentes em algumas regiões dos reatores nucleares, e.g., moderador, refletor, e meios estruturais. Nesses domínios espaciais não há geração de potência nuclear térmica e, além disso, comprometem a eficiência computacional dos cálculos globais de reatores nucleares. Propomos nesta tese uma estratégia visando a aumentar a eficiência computacional dessas simulações eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas (baffle e refletor) em torno do núcleo ativo. Apresentamos algumas modelagens e discutimos a eficiência da aplicação dessas condições de contorno aproximadas tipo albedo para uma e duas regiões nãomultiplicativas, na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. A denominação Albedo, palavra de origem latina para alvura, foi originalmente definida como a fração da luz incidente que é refletida difusamente por uma superfície. Esta denominação latina permaneceu como o termo científico usual em astronomia e, nesta tese, este conceito é estendido para reflexão de nêutrons. Estas condições de contorno tipo albedo SN não-convencional substituem aproximadamente as regiões de baffle e refletor no em torno do núcleo ativo do reator, desprezando os termos de fuga transversal no interior dessas regiões. 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 tese, são exatas. Por eficiência computacional entende-se a análise da 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 considerando dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

Estimativa de fontes externas uniformes de nêutrons em regiões combustíveis para uma distribuição prescrita de potência

Projetos de reatores nucleares foram classificados em quatro gerações (Gen) pelo Departamento de Energia dos Estados Unidos da América (DOE), quando o DOE introduziu o conceito de reatores de geração IV (Gen IV). Reatores Gen IV são um conjunto de projetos de reator nuclear, em sua maioria teóricos, atualmente sendo pesquisados. Entre os projetos Gen IV, incluem-se os projetos dos ADS (Accelerator Driven Systems), que são sistemas subcríticos estabilizados por fontes externas estacionárias de nêutrons. Estas fontes externas de nêutrons são normalmente geradas a partir da colisão de prótons com alta energia contra os núcleos de metais pesados presentes no núcleo do reator, fenômeno que é conhecido na literatura como spallation, e os prótons são acelerados num acelerador de partículas que é alimentado com parte da energia gerada pelo reator. A criticalidade de um sistema mantido por reações de fissão em cadeia depende do balanço entre a produção de nêutrons por fissão e a remoção por fuga pelos contornos e absorção de nêutrons. Um sistema está subcrítico quando a remoção por fuga e absorção ultrapassa a produção por fissão e, portanto, tende ao desligamento. Entretanto, qualquer sistema subcrítico pode ser estabilizado pela inclusão de fontes estacionárias de nêutrons em seu interior. O objetivo central deste trabalho é determinar as intensidades dessas fontes uniformes e isotrópicas de nêutrons, que se deve inserir em todas as regiões combustíveis do sistema, para que o mesmo estabilize-se gerando uma distribuição prescrita de potência elétrica. Diante do exposto, foi desenvolvido neste trabalho um aplicativo computacional em linguagem Java que estima as intensidades dessas fontes estacionárias de nêutrons, que devem ser inseridas em cada região combustível para que estabilizem o sistema subcrítico com uma dada distribuição de potência definida pelo usuário. Para atingir este objetivo, o modelo matemático adotado foi a equação unidimensional de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas (SN) e o convencional método de malha fina diamond difference (DD) foi utilizado para resolver numericamente os problemas SN físicos e adjuntos. Resultados numéricos para dois problemas-modelos típicos são apresentados para ilustrar a acurácia e eficiência da metodologia proposta.

Anomalia na densidade em um gás de rede com interações competitivas

Apesar de ser um líquido comum na natureza, muitas dúvidas ainda pairam sobre várias características da água. A existência de uma relação comum entre o tipo de potencial intermolecular, criticalidade e as várias formas de anomalia existentes nessa substância ainda é uma questão em aberto, apesar da intensa pesquisa que têm-se feito ao longo dos anos sobre esse assunto. Nesta dissertação, propomos a hipótese de que a anomalia na densidade esta correlacionada à presença de multicriticalidade e que ambos os fenômenos surgem de um potencial de duas escalas. Para dar suporte a esta hipótese, além de trabalhos anteriores, usamos e estudamos um gás de rede com interações que competem (primeiros vizinhos Vi atrativos e segundos vizinhos V2repulsivos). Construímos para este sistema um diagrama de fases J.lV8. T usando dois métodos: aproximação de campo médio e simulações. Encontramos na aproximação de campo médio duas linhas críticas, uma das quais encontra a linha de 1~ ordem separando duas fases líquidas, e um ponto tricrítico, dado que V2/V1 < -0.5. Se V2/V1 > -0.5, a transição líquido-líquido desaparece, dando lugar a apenas duas fases, uma líquido e uma gás, separadas por uma linha de coexistência terminada em um ponto crítico Com a aproximação de campo médio não encontramos anomalia na densidade. Os resultados obtidos com as simulações alteram qualitativamente o diagrama de fases. Tanto as linhas críticas quanto os pontos tricríticos tem suas posições modificadas com relação ao campo médio. Neste caso encontramos um comportamento anômalo na densidade se V2/Vl < -0.5. Concluímos que o potencial de duas escalas competitivas é um ingrediente necessário ao aparecimento de anomalia na densidade e coexistência entre duas fases líquidas. Ainda, mostramos que essa anomalia pode estar associada não apenas a dois pontos críticos, como se espera para a água, mas a uma multicriticalidade em geral, tal como linhas críticas.

Anomalia na densidade em um gás de rede com potencial repulsivo atenuado

Apesar de ser um líquido comum na natureza, muitas dúvidas ainda pairam sobre várias características da água. A existência de uma relação comum entre o tipo de potencial intermolecular, criticalidade e as várias formas de anomalia existentes nessa substância ainda é uma questão em aberto, apesar da intensa pesquisa que têm-se feito ao longo dos anos sobre esse assunto . Nesta dissertação, propomos a hipótese de que a anomalia na densidade está correlacionada à presença de multicriticalidade, e que ambos os fenômenos surgem de um potencial de duas escalas. Para dar suporte a esta hipótese, além de trabalhos anteriores, usamos e estudamos um gás de rede com interações que competem (primeiros vizinhos V1 repulsivos e segundos vizinhos V2 atrativos). Construímos para este sistema um diagrama de fases µ vs. T usando dois métodos: aproximação de campo médio e simulações. Encontramos na aproximação de campo médio duas linhas críticas, uma das quais encontra a linha de 1a ordem separando duas fases líquidas, e um ponto tricrític Com a aproximação de campo médio não encontramos anomalia na densidade. Os resultados obtidos com as simulações alteram qualitativamente o diagrama de fases. Tanto as linhas críticas quanto os pontos tricríticos tem suas posições modificadas com relação ao campo médio. Neste caso encontramos um comportamento anômalo na densidade. Concluímos que o potencial de duas escalas competitivas é um ingrediente necessário ao aparecimento de anomalia na densidade e coexistência entre duas fases líquidas. Ainda, mostramos que essa anomalia pode estar associada não apenas a dois pontos críticos, como se espera para a água, mas a uma multicriticalidade em geral, tal como linhas críticas.

Solução LTSn para problemas de transferência radiativa com polarização em geometria plana

O método LTSN tem sido utilizado na resolução de uma classe abrangente de problemas de transporte de partículas neutras que são reduzidos a um sistema linear algébrico depois da aplicação da transformada de Laplace. Na maioria dos casos estudados os autovalores associados são reais e simétricos. Para o problema de criticalidade os autovalores associados são reais ou imaginários puros e simétricos, e para o o problema de multigrupo podem aparecer autovalores complexos. O objetivo deste trabalho consiste na generalização da formulação LTSN para problemas de transporte com autovalores complexos. Por esse motivo é focada a solução de um problema radiativo de transporte com polarização em uma placa plana. A solução apresentada fundamenta-se na aplicação da transformada de Laplace ao conjunto de equações SN dos problemas resultantes da decomposição da equação de transferência radiativa com polarização em série de Fourier, seguindo o procedimento de Chandrasekhar. Esse procedimento gera 2L + 2 sistemas lineares de ordem 4N dependentes do parâmetro complexo "s". Aqui, L é o grau de anisotropia e N a ordem de quadratura. A solução desse sistema simbólico é obtida através da aplicação da transformada inversa de Laplace depois da inversão da matriz simbólica pelo método da diagonalização. Para a obtenção das constantes de integração é assumido que os componentes do vetor de Stokes são reais e as matrizes dos autovalores e autovetores são separadas em suas partes real e imaginária. A solução LTSN para autovalores complexos é validada através da comparação da solução para uma placa com espessura unitária, grau de anisotropia L = 13, albedo de espalhamento simples $ = 0:99, coe ciente de re exão de Lambert ¸0 = 0:1 e N = 150, segundo dados da literatura consultada.

Utilizando padrões organizacionais e avaliação de risco para adaptar a metodologia de desenvolvimento de software

Este trabalho descreve PMT – Pattern-based Methodology Tailoring, uma abordagem para a adaptação de metodologias de desenvolvimento de software, baseada em padrões e em critérios de risco. Seu principal objetivo é estabelecer meios de se adaptar uma linguagem de padrões organizacionais ao contexto de um projeto específico, o que é obtido através da seleção sistemática dos padrões organizacionais mais adequados aos requisitos do projeto. O trabalho é motivado pelo levantamento de que os arcabouços de processos de software existentes pouco fazem para compreender as necessidades de um projeto antes de definir a metodologia a ser aplicada. PMT utiliza uma análise dos riscos e do contexto de criticalidade para guiar o processo de adaptação. Padrões organizacionais que descrevem técnicas preventivas para os riscos identificados são selecionados por um mecanismo sistemático de seleção, o qual é suportado por uma ferramenta, chamada PMT-Tool.

Estudo da difusão translacional em um modelo para água

Apesar da água ser o líquido mais comum na natureza, suas características ainda não estão totalmente explicadas. A relação entre a forma do potencial intermolecular efetivo que representa as interações presentes na água, as várias anomalias existentes na água e a possível existência de dupla criticalidade ainda é uma questão em aberto. Recentemente descobriu-se que a água apresenta, além de anomalias termodinâmicas, anomalia na difusão translacional e rotacional. Mostrou-se que para um modelo SPC/E para água estas anomalias dinâmicas estão conectadas com a temperatura de máxima densidade (TMD) e que as anomalias dinâmicas ocupam uma região maior que a TMD no diagrama de fases p vs. T. Nesta dissertação investigamos a relação entre a anomalia na densidade e a anomalia na difusão translacional da água em um modelo onde as moléculas interagem através de um potencial que contém duas fontes de interação: uma isotrópica atrativa do tipo van der Waals, v, e uma direcional, u, que pode ser nula ou repulsiva. Reproduzimos o diagrama de fases pressão vs. temperatura aonde a anomalia na densidade e a dupla criticalidade se fazem presente. Mostramos que este sistema apresenta anomalia na difusão translacional e que esta se encontra na mesma região do diagrama de fases pressão vs. temperatura em que a TMD está presente. Diferentemente do que ocorre para água SPC/E, a região de anomalia na difusão translacional ocupa um espaço menor, encontrando-se na parte interna da TMD. Sugerimos que a discretização do potencial afeta mais fortemente a anomalia dinâmica do que a termodinâmica e que leva a um encolhimento na região de anomalia na difusão.

Undersampled Critical Branching Processes on Small-World and Random Networks Fail to Reproduce the Statistics of Spike Avalanches

The power-law size distributions obtained experimentally for neuronal avalanches are an important evidence of criticality in the brain. This evidence is supported by the fact that a critical branching process exhibits the same exponent t~3=2. Models at criticality have been employed to mimic avalanche propagation and explain the statistics observed experimentally. However, a crucial aspect of neuronal recordings has been almost completely neglected in the models: undersampling. While in a typical multielectrode array hundreds of neurons are recorded, in the same area of neuronal tissue tens of thousands of neurons can be found. Here we investigate the consequences of undersampling in models with three different topologies (two-dimensional, small-world and random network) and three different dynamical regimes (subcritical, critical and supercritical). We found that undersampling modifies avalanche size distributions, extinguishing the power laws observed in critical systems. Distributions from subcritical systems are also modified, but the shape of the undersampled distributions is more similar to that of a fully sampled system. Undersampled supercritical systems can recover the general characteristics of the fully sampled version, provided that enough neurons are measured. Undersampling in two-dimensional and small-world networks leads to similar effects, while the random network is insensitive to sampling density due to the lack of a well-defined neighborhood. We conjecture that neuronal avalanches recorded from local field potentials avoid undersampling effects due to the nature of this signal, but the same does not hold for spike avalanches. We conclude that undersampled branching-process-like models in these topologies fail to reproduce the statistics of spike avalanches.