Modelo OLAM (ocean-land-atmosphere-model): descrição, aplicações, e perspectivas

O modelo OLAM foi desenvolvido com objetivo de estender a capacidade de representar os fenômenos de escala global e regional simultaneamente. Este modelo apresenta inovações quanto aos processos dinâmicos, configuração de grade, estrutura de memória e técnicas de solução numérica das equações prognósticas. As equações de Navier-Stokes são resolvidas através da técnica de volumes finitos que conservam massa, momento e energia. No presente trabalho, apresenta-se uma descrição sucinta do OLAM e alguns resultados de sua aplicação em simulações climáticas da precipitação mensal para a região norte da América do Sul, bem como em rodadas para previsão numérica de tempo regional. Os resultados mostram que o modelo consegue representar bem os aspectos meteorológicos de grande escala. Em geral, seu desempenho melhora quando são adotadas grades de maior resolução espacial, nas quais se verificam melhorias significativas tanto na estimativa da precipitação mensal regional, quanto na previsão numérica de tempo.

Análise computacional dos campos de velocidade e temperatura do ambiente interno da usina termelétrica Santana - Amapá

A simulação numérica do escoamento de ar em ambientes internos é na atualidade o método mais apropriado para análise de conforto térmico em ambientes internos. O escoamento de ar nesses ambientes configura-se como um escoamento complexo, pois, em regra geral, é uma combinação de escoamentos cisalhantes livres (jatos) e cisalhantes de parede, além disso, esses escoamentos são governados por forças de inércia e forças de empuxo, caracterizando-o como de convecção mista. A combinação desses mecanismos cria um escoamento com características complexas, como zonas de recirculação, vórtices, descolamento e recolamento de camada-limite dentre outras. Portanto, a precisão da solução estará diretamente ligada, principalmente, na habilidade do modelo de turbulência adotado de reproduzir as características turbulentas do escoamento de ar e da transferência térmica. O objetivo principal do presente trabalho foi a simulação computacional do ambiente térmico interno do galpão que abriga os geradores e motores Wärtzilä da Usina Termelétrica Santana no estado do Amapá. A formulação matemática baseada na solução das equações gerais de conservação inclui uma análise dos principais modelos de turbulência aplicados ao escoamento de ar em ambientes internos, assim como os processos de transferência de calor associados. Na modelagem numérica o método de volumes finitos é usado na discretização das equações de conservação, através do código comercial Fluent-Airpak, que foi usado nas simulações computacionais para a análise dos campos de velocidade e temperatura do ar. A utilização correta do programa computacional foi testada e validada para o problema através da simulação precisa de casos retirados da literatura. Os resultados numéricos foram comparados a dados obtidos de medições experimentais realizados no galpão e apresentou boa concordância, considerando a complexidade do problema simulado, o objetivo da simulação em face da diminuição da temperatura no interior do galpão e, também, em função das limitações encontradas quando da tomada das medições experimentais. Além disso, foram feitas simulações de estratégias de melhoria do ambiente térmico da Usina, baseadas na realidade levantada e nos resultados da simulação numérica. Finalmente, foram realizadas simulações do protótipo de solução proposto para a diminuição da temperatura interna do galpão o que possibilitará um aumento, na faixa de 20 a 30%, do tempo de permanência no interior do galpão.

Influência da convecção termossolutal nos espaçamentos dendríticos primários durante a solidificação direcional horizontal de ligas Sn-Pb

Este trabalho apresenta um estudo teórico-experimental sobre os efeitos da convecção termossolutal sobre os espaçamentos dendríticos primários de ligas hipoeutéticas Sn-Pb. Assim, um dispositivo experimental de configuração horizontal refrigerado a água foi construído e amostras do referido sistema de ligas foram solidificadas direcionalmente sob condições transientes de extração de calor. Um método teórico-experimental é aplicado para determinar as velocidades de deslocamento da isoterma liquidus, os gradientes térmicos e as taxas de resfriamento. São também apresentados resultados para os respectivos coeficientes de transferência de calor na interface metal/molde os quais foram calculados a partir de uma análise comparativa realizada entre os perfis experimentais de temperatura e valores teóricos fornecidos por um método numérico baseado em volumes finitos. Alguns modelos teóricos para a previsão dos espaçamentos dendríticos primários são comparados com os resultados experimentais obtidos. Finalmente, um estudo comparativo é realizado entre resultados encontrados neste trabalho e aqueles apresentados na literatura para os espaçamentos dendríticos primários das ligas Sn-Pb investigadas quando solidificadas direcionalmente nos sistemas verticais ascendente e descendente, sob as mesmas condições. A análise das microestruturas indica que os espaçamentos dendríticos primários são bastante influenciados pela direção de crescimento do sólido.

Um modelo de turbulência baseado no conceito de vórtice

O fenômeno da turbulência está presente na maioria dos escoamentos observados na indústria e na natureza. Muitas são as considerações a respeito das dificuldades relacionadas à caracterização dos escoamentos turbulentos. Uma das muitas questões trata do procedimento de análise do problema através da descrição estatística dos campos por grandezas “médias”, o que leva ao problema de fechamento e à modelagem do tensor de Reynolds, normalmente com modelos baseados no conceito de viscosidade turbulenta. Os modelos de turbulência já existentes apresentam algumas deficiências na previsão do escoamento, além de outras limitações, o que justifica a busca por novas abordagens para o tratamento da turbulência. Neste trabalho, o problema de fechamento é tratado segundo a modelagem turbulenta baseada no conceito de viscosidade turbulenta. Um novo modelo de turbulência é proposto, que admite a existência de vórtices imersos no escoamento e aplica conceitos e definições relacionados à identificação de vórtices, com o uso do critério de identificação Q , que caracteriza a região do escoamento ocupada pelo vórtice. Propõe-se a investigação da aplicabilidade do critério Q em conjunto com o modelo k − ε , para o desenvolvimento de um novo modelo de turbulência chamado k − ε −Q . Validou-se a aplicabilidade do modelo através de um código numérico computacional para tratamento de escoamentos turbulentos. A solução numérica foi obtida através da discretização do domínio fluido, utilizando o método de volumes finitos e o método multigrid foi utilizado para resolver o sistema linear resultante. Como verificação, foi utilizado este modelo de turbulência para simular o escoamento em uma cavidade quadrada com tampa deslizante e o escoamento turbulento sobre um degrau. Os resultados obtidos foram confrontados com dados experimentais e demonstraram que o modelo aqui proposto se apresenta mais eficiente que o clássico modelo k − ε , no tratamento da turbulência nesses dois problemas clássicos.

Análise termoelástica experimental e computacional por elementos finitos da parede externa em concreto armado de um forno de cozimento de anodo

Nesta dissertação foi desenvolvido um estudo sobre o comportamento de uma estrutura em concreto armado submetida a elevadas variações térmicas. A estrutura analisada consiste na parede externa de um forno de cozimento de anodo utilizado na produção de Alumínio, e está submetida a variações térmicas de cerca de 125C devido à própria operação do forno. As motivações principais para o estudo foram a ocorrência de grandes deformações e o surgimento de fissuras na estrutura, o que poderia provocar a inutilização do forno. O objetivo do trabalho foi a investigação das causas destas patologias, assim como o estabelecimento de propostas para reforço da estrutura. Para isso, foram realizados estudos experimentais e computacionais do comportamento da estrutura. Inicialmente, a estrutura foi monitorada utilizando-se transdutores de deslocamento e sensores de temperatura (termopares), conectados a um sistema de aquisição de dados para obtenção e armazenamento automática das amostras ao longo do tempo. Em seguida, foram desenvolvidos modelos computacionais em Elementos Finitos com auxilio do programa computacional Algor, para determinação da distribuição de temperatura e as correspondentes tensões e deformações de origem térmica na estrutura. Nestes estudos, foram realizadas análises estacionárias e transientes de condução de calor, seguidas de análises de tensões de origem térmica. Como conclusão do estudo, tem-se que a metodologia proposta para a solução do problema foi bastante satisfatória, solucionando o problema de forma precisa, porém econômica.

Modelagem de transformadores de núcleo amorfo usando método de elementos finitos

Os materiais amorfos vêm ganhando grande espaço na indústria de transformadores, devido às baixas perdas no núcleo, visto que estes possuem um ciclo de histerese mais estreito quando comparados com os núcleos tradicionais de aço silício. Entretanto, seu custo ainda tem sido um grande fator para pouca inserção deste tipo de equipamento nos sistemas elétricos de potência. Estudos sobre o custo/benefício em termos técnicos de desempenho e robustez devem ser considerados quando se deseja projetar transformadores que utilizam esse tipo de material. A grande contribuição deste trabalho encontra-se na análise do desempenho do transformador de núcleo amorfo diante de curtos-circuitos, visto que tais falhas neste equipamento ocasionam redução da receita, não apenas por gastos com manutenção, mas também porque a concessionária deixa de vender seu produto, energia elétrica, além de poder estar sujeita a penalidades por parte dos órgãos de regulação do setor elétrico. Quando em condições de curto-circuito os enrolamentos dos transformadores ficam submetidos a esforços mecânicos, produzidos por forças de Lorentz, essas forças surgem como resultado do fluxo produzido pelos próprios condutores em paralelo que transportam corrente na mesma direção. Diante disso, estudar o comportamento eletromagnético do transformador é fundamental para obtenção de tais forças. Para o desenvolvimento deste trabalho, utilizou-se o software Finite Element Method Magnetics (FEMM). Esta ferramenta se baseia no método de elementos finitos para realizar os cálculos das magnitudes eletromecânicas e, consequentemente, o cálculo das forças atuando nas espiras, as quais permitem realizar os cálculos dos esforços mecânicos. Por fim, este trabalho aborda a aplicação da ferramenta FEMM para o cálculo de esforços mecânicos e simulação do comportamento eletromagnético de um transformador de distribuição.

O emprego da geoestatística na determinação do tamanho "ótimo" de amostras aleatórias com vistas à obtenção de estimativas dos volumes dos fustes de espécies florestais em Paragominas, estado do Pará

O objetivo do trabalho foi determinar o tamanho adequado de amostra para estimar o volume de fustes de espécies florestais de uma população de árvores a serem cortadas no sistema de manejo florestal da empresa Cikel Brasil Verde Madeiras - Pará. Utilizaram-se as metodologias da amostragem sistemática e do estimador geoestatístico da krigagem ordinária com simulação sequencial, respectivamente para a escolha das amostras e estimação dos volumes dos fustes das árvores. Os resultados mostraram que os métodos podem ser utilizados no cálculo dos volumes de fustes de árvores. Entretanto, o método da krigagem apresenta um efeito de suavização, tendo como conseqüência uma subestimação dos volumes calculados. Neste caso, um fator de correção foi aplicado para minimizar o efeito da suavização. A simulação sequencial indicativa apresentou resultados mais precisos em relação à krigagem, uma vez que tal método apresentou algumas vantagens, tal como a não exigência de amostras com distribuições normais e ausência de efeito de suavização, característico dos métodos de interpolação.

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.

Estudo comparativo das técnicas de elementos finitos e equação integral na modelagem eletromagnética bidimensional

Os métodos numéricos de Elementos Finitos e Equação Integral são comumente utilizados para investigações eletromagnéticas na Geofísica, e, para essas modelagens é importante saber qual algoritmo é mais rápido num certo modelo geofísico. Neste trabalho são feitas comparações nos resultados de tempo computacional desses dois métodos em modelos bidimensionais com heterogeneidades condutivas num semiespaço resistivo energizados por uma linha infinita de corrente (com 1000Hz de freqüência) e situada na superfície paralelamente ao "strike" das heterogeneidades. Após a validação e otimização dos programas analisamos o comportamento dos tempos de processamento nos modelos de corpos retangulares variandose o tamanho, o número e a inclinação dos corpos. Além disso, investigamos nesses métodos as etapas que demandam maior custo computacional. Em nossos modelos, o método de Elementos Finitos foi mais vantajoso que o de Equação Integral, com exceção na situação de corpos com baixa condutividade ou com geometria inclinada.

Modelagem 2,5D por elementos finitos dos efeitos da topografia do terreno sobre dados obtidos com o método eletromagnético a multifrequência

Este trabalho consiste em realizar a modelagem, via elementos finitos (EF) 2,5D, do efeito da topografia do terreno sobre dados obtidos com o método eletromagnético a multi-frequência (EMMF). Este método usa como fonte uma grande espira quadrada de corrente elétrica com centenas de metros de lado, e como receptores, bobinas posicionadas na horizontal em alinhamento com o transmissor. A subsuperfície é representada por heterogeneidades bidimensionais imersas em um meio horizontalmente estratificado. A formulação, partindo das equações de Maxwell, é desenvolvida a partir da separação do campo eletromagnético em primário (campos no hospedeiro multi-estratificado) e secundário (diferença entre o campo total e o primário). O domínio discretizado é descrito por uma malha não estruturada, com elementos triangulares. Para calcular as componentes derivadas da solução de elementos finitos, em um determinado nó da malha, foi usada a média aritmética das derivadas das funções bases de EF em torno daquele nó. O código de modelagem construído permite quantificar e analisar como os gradientes topográficos influenciam as medidas dos campos eletromagnéticos gerados. A aplicação é a avaliação dessas influências sobre a componente radial do campo da espira na superfície terrestre, que é a componente empregada no método eletromagnético a multi-frequência (EMMF).

Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D

A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.

Modelagem direta bidimensional do método magnetotelúrico com o método dos elementos finitos de arestas

Fizemos a modelagem direta 2D do método magnetotelúrico (MT) com o método dos elementos finitos (MEF) de arestas em termos dos campos primários e secundários. Para usarmos modelos de maior complexidade e diminuirmos o custo computacional utilizamos malhas não estruturadas. Nas malhas utilizadas, introduzimos quatro nós em torno de cada estação MT, constituindo um quadrado alinhado nas direções dos eixos cartesianos