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.

Solução da equação de transporte multidimensional em geometria cartesiana e meio infinito usando derivada fracionária

Neste trabalho, foi construída uma forma integral para a solução das equações de transporte em uma, duas e três dimensões, considerando o núcleo de espalhamento de Klein-Nishina, espalhamento isotrópico e o núcleo de espalhamento de Rutherford, respectivamente, seguindo a mesma idéia proposta em trabalhos recentes, nos quais foi construída uma solução para a equação de transporte de nêutrons em geometria cartesiana, usando derivada fracionária. A metodologia consiste em igualar a derivada fracionária do fluxo angular à equação integral, determinar a ordem da derivada fracionária comparando o núcleo da equação integral com o da definição de Riemann-Liouville. Essa formulação foi aplicada ao cálculo de dose absorvida. São apresentadas soluções geradas a partir do emprego do método da derivada fracionária e comparadas a resultados disponíveis na literatura.

Competições magnéticas no sistema Fe xCo 1-xTa 2 O 6

São apresentados resultados obtidos a partir de difração de raios-X (DRX), difração de nêutrons (DN), susceptibilidade magnética (X(T)), magnetização (M(H)), espectroscopia Mõssbauer (EM) e calor específico (Cp) de amostras do sistema FexCo1-x Ta206. Difratogramas de DRX e de DN e as curvas M(H) indicam que as amostras estão bem cristalizadas e homogêneas, e que o sistema é uma solução sólida para toda faixa de substituição Fe -> Co. Os ajustes de DN revelam fases magnéticas com dois vetores de propagação (::I::~~ ~) para o CoTa206 e (~ O ~) e (O ~ ~) para o FeTa206' A segunda configuração permanece a mesma para as amostras ricas em Fe (0,46 :S x < 1,00), enquanto que as amostras ricas em Co (0,09 :S x < 0,46) apresentam configuração magnética indexada pelos vetores de propagação (::I::~~ O). O diagrama de fase temperatura vs. x exibe um ponto bicrítico em torno de T = 4,9 K e x = 0,46. A temperatura de Néel máxima das região rica em Fe é 9,5 K, e 7,1 K para a região rica em Co. No ponto bicrítico, o sistema mostra coexistência de ambas estruturas magnéticas. Esse comportamento bicrítico é interpretado como sendo induzido pelas competições entre as diferentes fases magnéticas e pela variação das propriedades cristalográficas.

Estrutura nuclear de estrelas compactas

Este trabalho tem como objetivo o estudo da matéria nuclear a altas densidades considerando-se as fases hadrônica e de quarks à temperatura nula e finita, com vistas a aplicações no estudo de propriedades estáticas globais de estrelas compactas. Parte dos cálculos apresentados nesta dissertação foram realizados por diferentes autores. Entretanto, em geral, estes trabalhos limitaram-se ao estudo da matéria nuclear em regiões de densidades e temperaturas específicas. Este estudo visa, por sua vez, o desenvolvimento de um tratamento amplo e consistente para estes sistemas, considerando-se diferentes regimes de densidade e temperatura para ambas as fases, hadrônica e de quarks. Buscamos com isso adquirir conhecimento suficiente que possibilite, não somente a ampliação do escopo dos modelos considerados, como também o desenvolvimento, no futuro, de um modelo mais apropriado à descrição de propriedades estáticas e dinâmicas de estrelas compactas. Ainda assim, este trabalho apresenta novos aspectos e resultados inéditos referentes ao estudo da matéria nuclear, como descrevemos a seguir. No estudo da matéria nuclear na fase hadrônica, consideramos os modelos da teoria quântica de campos nucleares desenvolvidos por J. D. Walecka, J. Zimanyi e S. A. Moszkowski, e por J. Boguta e A. R. Bodmer, e conhecidos, respectivamente, como Hadrodinâmica Quântica, ZM e Não-Linear. Nestes modelos a matéria nuclear é descrita a partir de uma formulação lagrangeana com os campos efetivos dos bárions acoplados aos campos dos mésons, responsáveis pela interação nuclear Neste estudo consideramos inicialmente a descrição de propriedades estáticas globais de sistemas nucleares de muitos corpos à temperatura nula, como por exemplo, a massa efetiva do núcleon na matéria nuclear simétrica e de nêutrons. A equação de estado da matéria de nêutrons possibilita a descrição de propriedades estáticas globais de estrelas compactas, como sua massa e raio, através da sua incorporação nas equações de Tolman, Oppenheimer e Volkoff (TOV). Os resultados obtidos nestes cálculos estão em plena concordância com os resultados apresentados por outros autores. Consideramos posteriormente o estudo da matéria nuclear com graus de liberdade de bárions e mésons à temperatura finita, com particular atenção na região de transição de fase. Para este estudo, incorporamos aos modelos considerados, o formalismo da mecânica estatística à temperatura finita. Os resultados obtidos, para as propriedades da matéria nuclear à temperatura finita, concordam também com os resultados obtidos por outros autores. Um aspecto inédito apresentado neste trabalho refere-se à incorporação de valores para os pontos críticos da transição de fase, ainda não determinados por outros autores. O comportamento do calor específico também é analisado de forma inédita nesta dissertação no tratamento utilizado com os modelos Não-Linear e ZM. Utilizamos a equação de estado da matéria de nêutrons à temperatura finita nas equações TOV, determinando propriedades globais de uma estrela protoneutrônica Observamos neste trabalho que ocorre um aumento da massa máxima da estrela com o aumento da temperatura, comportamento este já previsto por outros autores em diferentes modelos. Posteriormente incorporamos ao formalismo à temperatura finita, o equilíbrio químico, a presença de graus de liberdade leptônicos para elétrons e múons e a neutralidade de carga. Apresentamos nesta etapa do trabalho, uma forma alternativa para a incorporação destes ingredientes, baseada na determinação de uma fração relativa entre os potenciais químicos de prótons e nêutrons, à temperatura nula, extendendo este resultado à temperatura finita. Este procedimento permite a determinação da distribuição de núcleons e léptons no interior de uma estrela protoneutrônica, onde incluímos ainda a presença de neutrinos confinados. No estudo da matéria de quarks, consideramos o modelo de sacola do Massachussets Institute of Technology (MIT). Incorporando as equações TOV neste estudo, determinamos propriedades globais de estrelas de quarks, bem como a distribuição dos diferentes sabores de quarks no interior estelar. Como principal resultado, obtivemos uma equação de estado geral para a matéria hadrônica e de quarks, introduzida nas equações TOV, e analisamos a existência de estrelas híbridas. Os resultados obtidos nesta etapa do trabalho são totalmente coerentes com aqueles obtidos por outros autores.

Estudo de difusão de impurezas introduzidas por implantação iônica em polímeros

No presente trabalho estudamos de forma sistemática a difusão de impurezas em filmes poliméricos usando as técnicas de implantação iônica e análise por feixe de íons, retroespalhamento Rutherford e de perfil de profundidade por nêutrons. Com este propósito foram implantadas e realizadas medidas em diferentes intervalos de temperatura para diferentes sistemas como (Au, Ag) no fotoresiste AZ1350, (Bi, Eu, Er, B) no fotoresiste S1813 e (Xe, Kr) no termoplástico Poliestireno. Como resultado mostrou-se que para implantações em baixas energias e fluências Au, Ag seguem uma difusão regular. Os valores obtidos para os parâmetros de difusão são semelhantes indicando assim um mecanismo de difusão verdadeiramente atômico. É mostrado também que com o aumento da fluência, devido aos danos gerados pelo processo de implantação, átomos são aprisionados na região implantada levando a um mecanismo de difusão por aprisionamento e liberação. Contudo, mostrou-se que a energia de ativação indica que este processo de difusão ainda é de caráter atômico. Da análise dos valores de D observamos um efeito de massa associado onde D(T)Au < D(T)Ag, pois a massa de Ag é duas vezes menor que a de Au. Para os elementos como Bi, Eu e Er, considerados quimicamente mais ativos que Au, não foram observados efeitos de possíveis ligações químicas nestes sistemas. Valores de energia de ativação de Bi apresentaram-se próximos aos de Au para as fluências de implantação aplicadas, o mesmo ocorrendo para a difusão de Er e Eu. O B mostrou que depois de implantado difunde durante ou imediatamente após a implantação. Difusão esta dada na presença de armadilhas saturáveis induzidas por radiação, indo além da difusão térmica, por ordem de magnitude de ≈10-12, contra ≈10-13 cm2s-1, respectivamente. Quanto a Xe e Kr, observou-se que estes também difundem durante ou imediatamente após a implantação, e que a fração do gás retido no pico depende da fluência implantada. Implantações em baixa temperatura (80 K) e posteriores análises foram determinados in situ por RBS na faixa de 80 a 300 K. Verificou-se que a difusão segue um perfil regular. Em cada caso mostrou-se que a dependência dos valores de D como função da temperatura segue um comportamento tipo Arrhenius, com valores de energia de ativação para os metais (Au, Ag, Bi) entre 580 e 680 meV, para os lantanídeos (Er, Eu) valores entre 525 a 530 meV, para o semi-metal (B) 100 meV, e finalmente para os gases nobres Kr e Xe entre 40 e 67 meV.

Estudo computacional de líquidos iônicos do tipo dialquilimidazólio

Neste trabalho foram desenvolvidos parâmetros de campo de força dos ânions tetrafenilborato (BPh4 -) e hexafluorfosfato (PF6 -) dentro da metodologia AMBER para a simulação computacional por dinâmica molecular de Líquidos Iônicos formados por estes ânions e cátions do tipo dialquilimidazólio, o 1-n-butil-3-metilimidazólio (BMI+), o 1-etil-3-metilimidazólio (EMI+) e o 1,3-dimetilimidazólio (MMI+). A validação destes parâmetros foi realizada por comparação entre as freqüências dos modos normais obtido por cálculo ab initio com aquelas obtidas por mecânica molecular, juntamente com uma comparação entre as estruturas moleculares e momentos multipolares obtidos pelas duas metodologias. Seguiu-se então a validação por comparação dos resultados dos cálculos de dinâmica molecular com dados experimentais, como densidades, entalpias de vaporização, condutividade elétrica, estrutura radial e espacial e também dados de difração de nêutrons. Foi atingida uma concordância bastante grande entre dados experimentais e cálculo teórico principalmente no que diz respeito à estrutura dos Líquidos Iônicos e foi possível racionalizar em termos de tamanho do grupo alquila do cátion (n-butil, metil ou etil) e do tamanho do ânion tanto propriedades estruturais deste líquidos quanto características dinâmicas dos mesmos.

Solução analítica das equações difusivas da teoria geral de perturbação pelo método da transformada de Laplace

Neste trabalho, apresentamos uma solução analítica para as equações difusivas unidimensionais da Teoria Geral de Perturbação em uma placa heterogênea, isto é, apresentamos as soluções analíticas para os problemas de autovalor para o fluxo de nêutrons e para o fluxo adjunto de nêutrons, para o cálculo do fator de multiplicação efetivo (keff), para o problema de fonte fixa e para o problema de função auxiliar. Resolvemos todos os problemas mencionados aplicando a Transformada de Laplace em uma placa heterogênea considerando um modelo de dois grupos de energia e realizamos a inversão de Laplace do fluxo transformado analiticamente através da técnica da expansão de Heaviside. Conhecendo o fluxo de nêutrons, exceto pelas constantes de integração, aplicamos as condições de contorno e de interface e resolvemos as equações algébricas homogêneas para o fator de multiplicação efetivo pelo método da bissecção. Obtemos o fluxo de nêutrons através da avaliação das constantes de integração para uma potência prescrita. Exemplificamos a metodologia proposta para uma placa com duas regiões e comparamos os resultados obtidos com os existentes na literatura.

Um modelo de sacola difusa para a matéria nuclear

Neste trabalho desenvolvemos um modelo efetivo para a descrição da matéria nuclear, que incorpora os resultados obtidos, para a descrição de um núcleon, pelo modelo de sacola difusa. O sistema nuclear será descrito via uma função de energia interna, que compreende um termo livre e outro que leva em conta a interação entre os núcleons. A parte livre, por se tratar de um sistema de férmions, corresponderá à energia de um gásde Fermi livre. Além disso, para evitar a superposição de dois ou mais núcleons, introduzimos um volume de exclusão a la Van der Waals. Na parte integrante, a troca de píons entre os núcleons será levada em conta via um potêncial efetivo. A função energia interna dependerá da densidade da matéria nuclear e também de um parâmetro que determinará o volume esperado de cada núcleon na matéria nuclear. O valor deste parâmetro será um pouco diferente do valor encontrado para um núcleons isolado, devido à interação entre eles. Obtém-se então resultados para a energia de ligação por núcleon para a matéria nuclear simétrica e para a matéria de nêutrons, bem como para a equação de estado da matéria de nêutrons.

Measurement systems of soil water matric potential and evaluation of soil moisture under different irrigation depths

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Efeito de campo aleatório no modelo Sherrington-Kirkpatrick

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior

TL, OSL and C-14 dating results of the sediments and bricks from mummified nuns' grave

Este artigo apresenta os resultados de datação por TL e OSL de solos, e fragmentos de tijolos de um túmulo, que foi ocupado por duas freiras mumificadas encontradas no Mosteiro da Luz, localizado no Estado de São Paulo, Brasil. As idades encontradas por TL e OSL foram comparadas às obtidas a partir de C-14 dos colágenos contidos em amostras de osso das múmias. A maioria das idades obtidas são do século XVIII. A espectroscopia de radiação-gama foi utilizada para avaliar concentrações de radioisótopos naturais nas amostras e para calcular as taxas de dose anual que resultaram em 3,0 a 5,3 Gy/kano. As concentrações radioativas são próximas daquelas obtidas através de Análise por Ativação de Nêutrons. Os conteúdos de elementos U, Th e Ce são superiores aos encontrados na maioria dos sedimentos.

Variação entre a lâmina de irrigação aplicada e armazenada no solo sob condições de irrigação com baixa uniformidade de distribuição de água

In sprinkler irrigation is important to have a good uniformity of application so that the water in the root zone does not show areas with little available water while others have percolating beyond the plant roots, even though the soil allows a lateral redistribution of water. A usual way to obtain the uniformity of irrigation is by measures of dispersion. The aim of this work was to evaluate the uniformity of water stored in a sand soil before and after non uniformity depth irrigation applied from one sprinkler, as well the variation of storage in a 0-0,45 m layer soil using a neutron probe. The statistical design was completely randomized, and for the variable Christiansen uniformity coefficient modified (CUCHH) were 10 treatments (the irrigation depth, 8 intervals of readings storage, more the existing one before irrigation added to the irrigation depth). For the variable soil water, the treatments were the same, excepting the irrigation depth. Despite the low surface uniformity (16.3%), there was not significantly difference between the storage uniformity before, after irrigation and the potential, however there was from these to the surface uniformity. From de irrigation depth, 15.3 mm, only 6.1 mm remained in the layer from 0 to 0.45 m. There was not significantly difference between the water stored in the soil before and after irrigation within a period of up to 134 hours, being the increase in storage due to irrigation was just 11.7%.

Estudo de N-metilformamida em meio não aquoso. O caso NMF-acetona

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo de N-metilformamida em meio aquoso e não aquoso (THF)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Férmions em referenciais acelerados: desintegração de prótons e outras aplicações

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)