Dinâmica da condensação de Bose-Einstein em gases fracamente interagentes

A presente dissertação estuda com detalhes a evolução temporal fora do equilíbrio de um condensado de Bose-Einstein homogêneo diluído imerso em um reservatório térmico. Nós modelamos o sistema através de um campo de Bose escalar complexo. É apropriado descrever o comportamento microscópico desse sistema por meio da teoria quântica de campos através do formalismo de Schwinger-Keldysh. Usando esse formalismo, de tempo real a dinâmica do condensado é solucionada por um grupo de equações integro-diferencial auto consistente, essas são solucionadas numericamente. Estudamos também o cenário quench, e como a densidade do gás e as interações entre as flutuações tem o efeito de provocar as instabilidades nesse sistema. Aplicamos esse desenvolvimento para estudar o comportamento de duas espécies homogêneas de um gás de Bose diluído imerso em um reservatório térmico.

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Otimização da programação das limpezas em redes de trocadores de calor utilizando técnicas de programação matemática

Redes de trocadores de calor são bastante utilizadas na indústria química para promover a integração energética do processo, recuperando calor de correntes quentes para aquecer correntes frias. Estas redes estão sujeitas à deposição, o que causa um aumento na resistência à transferência de calor, prejudicando-a. Uma das principais formas de diminuir o prejuízo causado por este fenômeno é a realização periódica de limpezas nos trocadores de calor. O presente trabalho tem como objetivo desenvolver um novo método para encontrar a programação ótima das limpezas em uma rede de trocadores de calor. O método desenvolvido utiliza o conceito de horizonte deslizante associado a um problema de programação linear inteira mista (MILP). Este problema MILP é capaz de definir o conjunto ótimo de trocadores de calor a serem limpos em um determinado instante de tempo (primeiro instante do horizonte deslizante), levando em conta sua influência nos instantes futuros (restante do horizonte deslizante). O problema MILP utiliza restrições referentes aos balanços de energia, equações de trocadores de calor e número máximo de limpezas simultâneas, com o objetivo de minimizar o consumo de energia da planta. A programação ótima das limpezas é composta pela combinação dos resultados obtidos em cada um dos instantes de tempo.O desempenho desta abordagem foi analisado através de sua aplicação em diversos exemplos típicos apresentados na literatura, inclusive um exemplo de grande porte de uma refinaria brasileira. Os resultados mostraram que a abordagem aplicada foi capaz de prover ganhos semelhantes e, algumas vezes, superiores aos da literatura, indicando que o método desenvolvido é capaz de fornecer bons resultados com um baixo esforço computacional

Construção de quadrados mágicos pelo método do passo uniforme

Lehmer (1929) analisa matematicamente o método do passo uniforme para construção de quadrados mágicos de ordem impar. Ele divide sua análise em várias etapas. Na primeira delas, envolvendo a discussão de condições necessárias e suficientes para o preenchimento do quadrado pelo método, o autor afirma que se dois números guardarem entre si uma certa relação, eles serão designados a ocupar a mesma célula do quadrado causando seu não preenchimento. A análise do preenchimento pelo método do passo uniforme envolve a resolução de um sistema linear módulo n. Nesse trabalho, discutimos o comportamento das soluções desse sistema quando o método falha no preenchimento. Como consequência, concluímos que números que guardam a relação mencionada nunca ocupam a mesma célula. A análise das condições necessárias e suficientes para obter quadrados mágicos segundo a definição de Lehmer (1929) envolve a resolução de equações de congruências lineares a duas variáveis. Nesse trabalho, detalhamos os resultados de Lehmer (1929). A análise das condições necessárias e suficientes para obtenção de quadrados mágicos, como são reconhecidos usualmente, também envolve a resolução de equações de congruências lineares a duas variáveis. Discutimos o comportamento das soluções dessas equações para obter diagonais principais mágicas. Como consequência, mostramos que diagonais principais mágicas são obtidas se e somente se as coordenadas iniciais guardarem certas relações

Átomos hidrogenoides em espaços com dimensionalidade D6=3: os casos não relativístico e relativístico

A possibilidade da existência de átomos de hidrogênio estáveis em dimensões superiores a três é abordada. O problema da dimensionalidade é visto como um problema de Física, no qual relacionam-se algumas leis físicas com a dimensão espacial. A base da análise deste trabalho faz uso das equações de Schrödinger (não relativística) e de Dirac (relativística). Nos dois casos, utiliza-se a generalização tanto do setor cinemático bem como o setor de interação coulombiana para variar o parâmetro topológico dimensão. Para o caso não relativístico, os auto-valores de energia e as auto-funções são obtidas através do método numérico de Numerov. Embora existam soluções em espaços com dimensões superiores, os resultados obtidos no presente trabalho indicam que a natureza deve, de alguma maneira, se manifestar em um espaço tridimensional.

Análise do mercado sucroalcooleiro e das elasticidades preço e renda da demanda por etanol hidratado

Este trabalho tem como objetivo estudar o mercado de etanol combustível e estimar as elasticidades preço e renda da demanda de etanol hidratado, no Brasil, no período de janeiro de 2003 a setembro de 2012. O método econométrico utilizado para analisar os dados e obter os resultados referentes à estimativa dos parâmetros das equações de demanda foi a metodologia de Johansen. Baseado na teoria da demanda incluiu-se, num primeiro momento, como variáveis explicativas: o preço do bem, o preço do bem substituto e a renda, num segundo momento incluiu-se também a frota de veículos, estimando assim o modelo VAR/VEC. Os resultados encontrados foram significativos e estão de acordo com a teoria econômica, nos levando a concluir que a demanda por etanol é bem elástica ao preço do etanol e ao preço da gasolina. Antes de desenvolvermos o modelo, analisamos o mercado de cana-de-açúcar, de etanol e da commodity concorrente ao etanol, o açúcar. Destacando a preocupação com o meio ambiente e a importância do etanol como energia renovável.

Caracterização microestrutural e análise de tensões residuais pelo método do furo cego em tubo de seção quadrada com costura.

Tensões residuais são uma das principais causas de falhas em componentes mecânicos submetidos a processos de fabricação. O objetivo do trabalho foi medir as tensões residuais presentes em um tubo quadrado soldado por resistência elétrica de alta frequência e caracterizar microestruturalmente o seu material. Para a caracterização, foram utilizadas técnicas de microscopia óptica (MO), microscopia eletrônica de varredura (MEV) e análise química quantitativa. Para a medição das tensões residuais, foi utilizado o método do furo cego, baseado na norma ASTM E837-08, onde rosetas (strain-gages) são coladas à peça para medir as deformações geradas devido à usinagem de um pequeno furo no local de medição. As deformações foram associadas às tensões residuais através de equações baseadas na Lei de Hooke. A caracterização revelou uma microestrutura composta basicamente de ferrita e perlita, típica de aços com baixo teor de carbono, corroborando com a especificação fornecida pelo fabricante. As tensões residuais encontradas foram trativas e mostraram-se elevadas, com alguns valores acima do limite de escoamento do material.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

O surgimento dos números irracionais

Este é um trabalho de pesquisa sobre um conjunto de números (irracionais) que é pouco trabalhado no ensino básico de matemática. Foi uma procura muito interessante e enriquecedora, pois encontrei matemáticos e historiadores com visões bem diferentes. Muitos deles não aceitavam este novo conjunto. Para Leopold Kronecker, só existia o conjunto dos números inteiros. Já para Cantor e Dedekind, o aparecimento dos irracionais foi extremamente importante para o desenvolvimento da matemática, abrindo novos horizontes. Menciono aqui um pouco da vida e da obra de alguns matemáticos que se envolveram com os números irracionais. Tratamos ainda da descoberta dos incomensuráveis, ou seja, como iniciou-se o problema da incomensurabilidade, e do retângulo áureo e sua importância em outras áreas. O trabalho mostra também dois grupos de números que não são mencionados quando ensinamos equações algébricas, que são os números algébricos e os números transcendentes, assim como teoremas essenciais para a prova da transcendência dos irracionais especiais e . Por fim, proponho uma aula para uma turma do 3 ano do Ensino Médio com o objetivo de mostrar a irracionalidade de alguns números, usando os teoremas pertinentes

Inflação em sistemas com fluidos de radiação com viscosidade e dissipação e análise de estabilidade

Este trabalho apresenta um estudo da estabilidade das equações da inflação morna com um fluido de radiação viscoso. A viscosidade do fluido é proveniente do constante decaimento de partículas neste, devido à dissipação do campo escalar da inflação, o ínflaton.Esta viscosidade, que pode ser volumar ou laminar, é tratada em termos de teorias termodinâmicas fora do equilíbrio. Este estudo se limita às equações de fundo da inflação morna, de modo que somente a viscosidade volumar tem um efeito significativo, sendo a viscosidade laminar importante somente no contexto de perturbações cosmológicas. A descrição da viscosidade em termos de uma termodinâmica fora do equilíbrio, porém, não pode ser realizada univocamente, pois a única informação que temos sobre processos irreversíveis é a segunda lei da termodinâmica. Portanto, parte-se em busca de teorias que estejam de acordo com esta lei e que, por argumentos plausíveis, sejam capazes de descrever o comportamento dos fluxos dissipativos próximo ao equilíbrio. O objetivo deste trabalho é estudar a estabilidade da inflação morna viscosa para teorias causais e não causais para o fluido de radiação com viscosidade, de forma que se possa observar o impacto da viscosidade no regime inflacionário e a relevância de se passar a considerar a causalidade. Para o fluido de radiação, as teorias consideradas são a teoria não causal de Eckart e as teorias causais de Israel-Stewart e de Denicol et al (hidrodinâmica dissipativa causal não linear). Obtém-se que as teorias causais, como era de se esperar, além de serem, por definição, consistentes no tocante à finitude da velocidade de propagação dos fluxos dissipativos, tornam o sistema dinâmico estável para valores de viscosidade mais distantes do equilíbrio. Observa-se também, nitidamente, que a teoria de Denicol et al é a mais robusta nesse sentido. Este trabalho, portanto, visa dar continuidade ao estudo dos efeitos não-isentrópicos na inflação, já que, além da dissipação do ínflaton na inflação morna, o impacto da viscosidade tem despertado bastante interesse.

Avaliação de ligações tubulares tipo T entre CHS com reforço tipo chapa.

Os perfis tubulares sem costura são utilizados em diversos países, principalmente devido às vantagens associadas à estética a sua elevada resistência à torção, cargas axiais e efeitos combinados. Canadá, Inglaterra, Alemanha e Holanda utilizam de forma veemente estas estruturas e possuem produção contínua e industrializada com alto nível de desenvolvimento tecnológico. O Brasil, porém, se limitava praticamente ao uso destes perfis nas coberturas espaciais. Devido ao aumento da utilização desses tipos de estruturas, fez-se necessário o aprofundamento dos estudos com métodos de análise coerentes para utilização de perfis tubulares, principalmente em relação às ligações, pois são consideradas regiões vulneráveis neste tipo de estrutura. Para atender a necessidade de normatização deste procedimento desenvolveu-se uma norma brasileira específica para o dimensionamento de estruturas em perfis tubulares. Considerando esta perspectiva, este trabalho apresenta uma análise de ligações tipo T com reforço tipo chapa com perfis tubulares circulares (CHS) para o banzo e para o montante efetuada com base na norma europeia, Eurocode 3, no CIDECT, na NBR 16239:2013 e ISO 14346. Desenvolveu-se no programa Ansys um modelo de elementos finitos para cada tipo de ligação analisada, calibrado e validado com resultados experimentais e numéricos existentes na literatura. Verificou-se a influência da compressão atuante no montante no comportamento global das ligações. As não-linearidades físicas e geométricas foram incorporadas aos modelos a fim de se mobilizar totalmente a capacidade resistente desta ligação. A nãolinearidade do material foi considerada através do critério de plastificação de von Mises através da lei constitutiva tensão versus deformação bilinear de forma a exibir um comportamento elasto-plástico com encruamento. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizada. A análise dos esforços resistentes obtidos em comparação com os resultados do modelo numérico, apresentaram valores a favor da segurança no cálculo utilizando as equações de dimensionamento. Por fim um estudo para fatores de correção das equações de dimensionamento foi também proposto.

Colapso gravitacional na teoria de Hořava-Lifshitz no limite infravermelho

O objetivo desta tese é verificar se algumas soluções de matéria da teoria da relatividade geral também satisfazem as equações da teoria de Horava-Lifshitz no limite infravermelho. Para isso, partimos das soluções mais simples possíveis, tais como é o caso de um fluido de radiação nula e de poeira, conhecidas na relatividade geral, e encontramos que estas não correspondem a quaisquer soluções na teoria de gravitação de Hořava-Lifshitz para o limite de baixas energias, infravermelho, no qual esta teoria deveria se reduzir à anterior. Este resultado nos remete a novos desafios na direção de ajustes teóricos que permitam que esta teoria descreva corretamente tanto o cenário cosmológico quanto o de formação de estruturas, estáveis ou colapsadas. Para tornar o trabalho mais claro, é feita uma introdução à teoria de Hořava- Lifshitz, tema central deste trabalho, e como ela se acopla com a matéria.

Equação de referência do Teste do Degrau medida em distância para indivíduos saudáveis e sedentários

O teste de caminhada em seis minutos (TC6M) avalia a capacidade respiratória durante o exercício. Recentemente, o teste do degrau em seis minutos (TD6M) está sendo estudado como uma proposta para essa mesma avaliação. Diante do exposto, o desenvolvimento de uma equação de referência se torna importante. O objetivo desse estudo foi desenvolver uma equação de referência padrão para o Teste do Degrau em Seis minutos. Esse estudo foi do tipo transversal, em que foram selecionados 452 indivíduos. Após a aplicação dos critérios de inclusão/exclusão, foram selecionados 326 sujeitos saudáveis e sedentários com idade entre 20 e 80 anos. Para serem considerados saudáveis, os participantes não podiam ter história de doenças (exceto hipertensão arterial sistêmica ou diabetes mellitus em tratamento) e foram submetidos à radiografia de tórax, espirometria e eletrocardiograma, que deveriam ser normais. O nível de sedentarismo foi avaliado através do International Physical Activity Questionnaire (IPAQ). Foram coletados os seguintes dados demográficos: idade, peso e altura. Todos os indivíduos realizaram TD6M na sua própria cadência (autocadenciado) em um degrau de 16,5 cm de altura, 65 de largura e 30 cm de comprimento. O número de subidas e descidas foi contado por um pedômetro digital. Foram mensuradas a pressão arterial, a frequência cardíaca e a saturação de oxigênio, antes e depois do TD6M. A análise estatística foi realizada pelo software STATA 12.0, e as equações foram desenvolvidas pelo modelo estatístico de regressão linear múltipla. Como resultado, observou-se que os participantes apresentaram exames normais, sendo 135 homens e 191 mulheres. O IPAQ demonstrou 157 ativos, 114 irregularmente ativos e 40 inativos, sendo que 14 indivíduos não responderam. A análise da distância alcançada em relação à idade e à diferença da frequência cardíaca tanto para homens quanto para mulheres mostrou significância estatística, demonstrando a sua importância para o desenvolvimento das equações para cada gênero. Para mulheres: Distância (m) = 88,83 - [(FC final- FC inicial)* 0,23] (Idade *0,37); para homens: Distância (m)= 110,20 - [(FC final- FC inicial)*0,18] (Idade * 0,59). O estudo concluiu que as equações de referência desenvolvidas nesse estudo foram realizadas em uma população de indivíduos saudáveis e sedentários e pode ser usada como padrão de referência para o TD6M.

Estudo de fluxos numéricos para modelos magneto-hidrodinâmicos utilizando o código Flash

A teoria magneto-hidrodinâmicos permite a estruturação de modelos computacionais, designados modelos MHDs, que são uma extensão da dinâmica dos fluidos para lidar com fluidos eletricamente carregados, tais como os plasmas, em que se precisa considerar os efeitos de forças eletromagnéticas. Tais modelos são especialmente úteis quando o movimento exato de uma partícula não é de interesse, sendo que as equações descrevem as evoluções de quantidades macroscópicas. Várias formas de modelos MHD têm sido amplamente utilizadas na Física Espacial para descrever muitos tipos diferentes de fenômenos de plasma, tais como reconexão magnética e interações de ventos estelares com diferentes objetos celestiais. Neste trabalho, o objetivo é analisar o comportamento de diversos fluxos numéricos em uma discretização de volumes finitos de um modelo numérico de MHD usando um esquema de malha entrelaçada sem separação direcional considerando alguns casos testes. Para as simulações, utiliza-se o código Flash, desenvolvido pela Universidade de Chicago, por ser um código de amplo interesse nas simulações astrofísicas e de fenômenos no espaço próximo à Terra. A metodologia consiste na inclusão de um fluxo numérico, permitindo melhoria com respeito ao esquema HLL.

Análise dos problemas direto e inverso em um modelo de evolução populacional através de transformações integrais e inferência bayesiana

Modelos de evolução populacional são há muito tempo assunto de grande relevância, principalmente quando a população de estudo é composta por vetores de doenças. Tal importância se deve ao fato de existirem milhares de doenças que são propagadas por espécies específicas e conhecer como tais populações se comportam é vital quando pretende-se criar políticas públicas para controlar a sua proliferação. Este trabalho descreve um problema de evolução populacional difusivo com armadilhas locais e tempo de reprodução atrasado, o problema direto descreve a densidade de uma população uma vez conhecidos os parâmetros do modelo onde sua solução é obtida por meio da técnica de transformada integral generalizada, uma técnica numérico-analítica. Porém a solução do problema direto, por si só, não permite a simulação computacional de uma população em uma aplicação prática, uma vez que os parâmetros do modelo variam de população para população e precisam, portanto, ter seus valores conhecidos. Com o objetivo de possibilitar esta caracterização, o presente trabalho propõe a formulação e solução do problema inverso, estimando os parâmetros do modelo a partir de dados da população utilizando para tal tarefa dois métodos Bayesianos.