Uso de sistema de informação geográfica para subsidiar o cálculo de estoque de carbono em florestas no âmbito do mecanismo de desenvolvimento limpo.

Com a entrada em vigor do Protocolo de Quioto, intensificam-se as expectativas pela regulamentação de um mercado de créditos de carbono. No caso de esses créditos terem sua origem em projetos de reflorestamento ou de florestamento, conforme previsto pelos chamados Mecanismos de Desenvolvimento Limpo (MDL), vem a ser fundamental o estabelecimento de uma metodologia para quantificação de estoques de carbono armazenados sob a forma de biomassa vegetal. Este trabalho propõe, como um método informatizado para cálculo de estoque de carbono em florestas, um conjunto de funcionalidades e um modelo de dados cadastrais totalmente integrados com um Sistema de Informações Geográficas de arquitetura aberta. A partir de mapas e imagens geo-referenciadas e com base em dados obtidos de pequenos transectos, o sistema calcula a quantidade total de carbono estocada em toda a floresta. Além de permitir apresentar esses resultados para diferentes agentes armazenadores de carbono, como, por exemplo, segmentos de floresta ou cada espécie vegetal, o sistema mantém registro histórico de dados dendrométricos, o que virá permitir a geração de gráficos de curvas de crescimento e, por conseguinte, estimativas futuras.

Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

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.