Implementação de Visualização de Dados Tridimensionais de Malhas Irregulares no Processador Cell Broadband Engine.

A renderização de volume direta tornou-se uma técnica popular para visualização volumétrica de dados extraídos de fontes como simulações científicas, funções analíticas, scanners médicos, entre outras. Algoritmos de renderização de volume, como o raycasting, produzem imagens de alta qualidade. O seu uso, contudo, é limitado devido à alta demanda de processamento computacional e o alto uso de memória. Nesse trabalho, propomos uma nova implementação do algoritmo de raycasting que aproveita a arquitetura altamente paralela do processador Cell Broadband Engine, com seus 9 núcleos heterogêneos, que permitem renderização eficiente em malhas irregulares de dados. O poder computacional do processador Cell BE demanda um modelo de programação diferente. Aplicações precisam ser reescritas para explorar o potencial completo do processador Cell, que requer o uso de multithreading e código vetorizado. Em nossa abordagem, enfrentamos esse problema distribuindo a computação de cada raio incidente nas faces visíveis do volume entre os núcleos do processador, e vetorizando as operações da integral de iluminação em cada um. Os resultados experimentais mostram que podemos obter bons speedups reduzindo o tempo total de renderização de forma significativa.

Cenários de mudanças climáticas usando modelagem dinâmica na Bacia do Alto Taquari.

A região Centro-Oeste do Brasil tornou-se nos últimos 40 anos grande produtora de grãos e carne bovina. As condições edafoclimáticas, o sistema de manejo do solo e o descumprimento de leis ambientais trouxeram conseqüências drásticas à região como o agravamento do processo hídrico erosivo, principalmente na Bacia do Alto Taquari (BAT). Cerca de 90% da BAT localiza-se na porção norte do estado de Mato Grosso do Sul (MS), porém os efeitos do transporte de sedimentos e volume de água são refletidos a jusante dos rios, na Bacia do Pantanal. Utilizando-se pressupostos do Painel Intergovernamental de Mudanças Climáticas (IPCC) foram estabelecidos cenários de mudanças climáticas na Bacia do Alto Taquari, visando identificar áreas com maior vulnerabilidade ao processo erosivo em função de pressões de uso da terra. Usando a modelagem dinâmica no TerraME (Environment Modeling) foram gerados cenários topopluviais até 2100, considerando-se para a temperatura do ar média anual um aumento de 1C, em cenário otimista e, em pessimista, elevações térmicas de 3C. Para a precipitação pluvial média anual um cenário foi com aumento de 15% e outro com reduções de 15%. Os dados foram espacializados no ArcGis 9.2 e exportados para o TerraView 3.2, criando-se espaços celulares e integrando-se com as informações do modelo digital do terreno do Shuttle Radar Topography Mission (SRTM) para geração dos mapas topoclimáticos e simulações de cenários no TerraMe. Os resultados apontam que 85% da área da BAT nas condições atuais as temperaturas médias variam entre 23,6 a 25,7C. As simulações térmicas no cenário otimista indicam que em 40 anos as temperaturas tendem a superar o maior limite térmico médio nas áreas ao longo do rio Taquari, no sentido Oeste-Leste. Esses valores evidenciam elevações nas taxas evapotranspiratórias de matas ciliares, indicando reduções na vazão do Taquari. Em cenário pessimista essas temperaturas antecipam sua ocorrência, em um prazo de 20 anos. Os cenários com acréscimo de 15% na precipitação pluvial mostram aumentos no volume de água precipitada na parte norte da Bacia, região mais vulnerável aos problemas de erosão hídrica. Cenários do regime térmico-hídrico apontam áreas mais sensíveis às mudanças climáticas na parte oeste da BAT e impactos ambientais também na Bacia do Pantanal. Conclui-se que o TerraME é indicado para gerar cenários de mudanças climáticas em bacias hidrográficas.

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.

Estabilidade de talude com efeito sísmico a partir dos métodos de equilíbrio limite e de elementos finitos.

A Costa Peruana apresenta alta atividade sísmica, tornando imprescindível a execução de análises de estabilidade que considerem os eventos sísmicos. Com o desenvolvimento de novas ferramentas numéricas, a análise dinâmica está se tornando cada vez mais importante e usual na fase de projeto, não se justificando mais a execução de análises estáticas em locais tão vulneráveis a sismos. A presente dissertação tem como objetivo realizar a análise da estabilidade dos Penhascos de Lima, Peru, considerando quatro taludes distintos, em uma região bastante afetada por abalos sísmicos. As análises foram executadas pelos métodos de equilíbrio limite, elementos finitos e pseudo-estático, buscando-se uma comparação entre os diferentes métodos. O trabalho apresenta uma descrição completa dos taludes em estudo, uma avaliação da condição sismológica da região e finalmente a estabilidade dos Penhascos de Lima, fazendo uso dos programas computacionais Slide (método do equilíbrio limite) e Plaxis (método dos elementos finitos) em 2D. Os resultados mostraram que os três métodos adotados forneceram fatores de segurança compatíveis, principalmente quando se considera perfis menos estratificados. Para perfis homogêneos, as diferenças obtidas foram da ordem de 0,5 a 1,0 %. As análises ressaltaram a importância de considerar a condição dinâmica, e mostraram-se bastante sensíveis aos valores de carga sísmica adotado.

Enfoque da estatística espacial em modelos dinâmicos de mudança do uso do solo.

O atual nível das mudanças uso do solo causa impactos nas mudanças ambientais globais. Os processos de mudanças do uso e cobertura do solo são processos complexos e não acontecem ao acaso sobre uma região. Geralmente estas mudanças são determinadas localmente, regionalmente ou globalmente por fatores geográficos, ambientais, sociais, econômicos e políticos interagindo em diversas escalas temporais e espaciais. Parte desta complexidade é capturada por modelos de simulação de mudanças do uso e cobertura do solo. Uma etapa do processo de simulação do modelo CLUE-S é a quantificação da influência local dos impulsores de mudança sobre a probabilidade de ocorrência de uma classe de uso do solo. Esta influência local é obtida ajustando um modelo de regressão logística. Um modelo de regressão espacial é proposto como alternativa para selecionar os impulsores de mudanças. Este modelo incorpora a informação da vizinhança espacial existente nos dados que não é considerada na regressão logística. Baseado em um cenário de tendência linear para a demanda agregada do uso do solo, simulações da mudança do uso do solo para a microbacia do Coxim, Mato Grosso do Sul, foram geradas, comparadas e analisadas usando o modelo CLUE-S sob os enfoques da regressão logística e espacial para o período de 2001 a 2011. Ambos os enfoques apresentaram simulações com muito boa concordância, medidas de acurácia global e Kappa altos, com o uso do solo para o ano de referência de 2004. A diferença entre os enfoques foi observada na distribuição espacial da simulação do uso do solo para o ano 2011, sendo o enfoque da regressão espacial que teve a simulação com menor discrepância com a demanda do uso do solo para esse ano.

Localização colaborativa em robótica de enxame.

Diversas das possíveis aplicações da robótica de enxame demandam que cada robô seja capaz de estimar a sua posição. A informação de localização dos robôs é necessária, por exemplo, para que cada elemento do enxame possa se posicionar dentro de uma formatura de robôs pré-definida. Da mesma forma, quando os robôs atuam como sensores móveis, a informação de posição é necessária para que seja possível identificar o local dos eventos medidos. Em virtude do tamanho, custo e energia dos dispositivos, bem como limitações impostas pelo ambiente de operação, a solução mais evidente, i.e. utilizar um Sistema de Posicionamento Global (GPS), torna-se muitas vezes inviável. O método proposto neste trabalho permite que as posições absolutas de um conjunto de nós desconhecidos sejam estimadas, com base nas coordenadas de um conjunto de nós de referência e nas medidas de distância tomadas entre os nós da rede. A solução é obtida por meio de uma estratégia de processamento distribuído, onde cada nó desconhecido estima sua própria posição e ajuda os seus vizinhos a calcular as suas respectivas coordenadas. A solução conta com um novo método denominado Multi-hop Collaborative Min-Max Localization (MCMM), ora proposto com o objetivo de melhorar a qualidade da posição inicial dos nós desconhecidos em caso de falhas durante o reconhecimento dos nós de referência. O refinamento das posições é feito com base nos algoritmos de busca por retrocesso (BSA) e de otimização por enxame de partículas (PSO), cujos desempenhos são comparados. Para compor a função objetivo, é introduzido um novo método para o cálculo do fator de confiança dos nós da rede, o Fator de Confiança pela Área Min-Max (MMA-CF), o qual é comparado com o Fator de Confiança por Saltos às Referências (HTA-CF), previamente existente. Com base no método de localização proposto, foram desenvolvidos quatro algoritmos, os quais são avaliados por meio de simulações realizadas no MATLABr e experimentos conduzidos em enxames de robôs do tipo Kilobot. O desempenho dos algoritmos é avaliado em problemas com diferentes topologias, quantidades de nós e proporção de nós de referência. O desempenho dos algoritmos é também comparado com o de outros algoritmos de localização, tendo apresentado resultados 40% a 51% melhores. Os resultados das simulações e dos experimentos demonstram a eficácia do método proposto.

Modelo de Alocação de Mamógrafos Móveis: Um estudo para a Região Serrana do Estado do Rio de Janeiro.

O Estado do Rio de Janeiro possui indicadores de produção muito baixos na realização de exames de câncer de mama. Na tentativa de melhorar o acesso aos exames, principalmente em regiões com baixa densidade populacional onde a aquisição de mamógrafos não é custo-efetiva, o uso da mamografia móvel é uma alternativa para aumentar a execução de exames de rastreamento de câncer de mama. O objetivo desta pesquisa é a construção de um modelo computacional para definir a alocação de mamógrafos móveis. O Modelo considera as variáveis associadas com os custos e prazos, indicando quando, onde e por quanto tempo, as unidades móveis de mamografia devem permanecer em cada cidade. O modelo foi construído no software de modelagem e simulação Anylogic, usando técnicas de modelagem baseada em agentes. O principal resultado é determinar o percurso de cada veículo disponível, para oferecer a cobertura desejada em cada cidade. Todas as entradas são parametrizadas, permitindo simular diferentes cenários e fornecer informações importantes para o processo de tomada de decisão. O horizonte de tempo, número de mamógrafos (fixos e móveis), a cobertura desejada da população, a capacidade de produção de cada dispositivo, a adesão da população urbana e rural, entre outras variáveis, foram consideradas no modelo. Os dados da Região Serrana do Rio de Janeiro foram usados nas simulações, onde menos de metade das cidades possuem mamógrafos fixos. Com o modelo proposto foi possível determinar a distribuição de cada dispositivo físico e o número ótimo de unidades móveis de mamografia para oferecer cobertura à totalidade da população no ciclo de dois anos. O número de mamógrafos para oferecer cobertura de toda a população da região poderia ser reduzido pela metade com o modelo de alocação proposto neste trabalho. A utilização de mamografia móvel, em conjunto com a rede existente de mamógrafos fixos, procura maximizar a disponibilização de exames de testes de diagnóstico de câncer de mama no estado do Rio de Janeiro. O desenvolvimento de um modelo de roteamento que aperfeiçoa a cobertura de rastreio do câncer de mama é apresentado como um complemento importante na tentativa de melhorar o acesso à população residente em áreas urbanas e rurais dos municípios.

Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.

Análise de desempenho de um protocolo BitTorrent ciente de localização em redes corporativas.

Hoje em dia, distribuições de grandes volumes de dados por redes TCP/IP corporativas trazem problemas como a alta utilização da rede e de servidores, longos períodos para conclusão e maior sensibilidade a falhas na infraestrutura de rede. Estes problemas podem ser reduzidos com utilização de redes par-a-par (P2P). O objetivo desta dissertação é analisar o desempenho do protocolo BitTorrent padrão em redes corporativas e também realizar a análise após uma modificação no comportamento padrão do protocolo BitTorrent. Nesta modificação, o rastreador identifica o endereço IP do par que está solicitando a lista de endereços IP do enxame e envia somente aqueles pertencentes à mesma rede local e ao semeador original, com o objetivo de reduzir o tráfego em redes de longa distância. Em cenários corporativos típicos, as simulações mostraram que a alteração é capaz de reduzir o consumo médio de banda e o tempo médio dos downloads, quando comparados ao BitTorrent padrão, além de conferir maior robustez à distribuição em casos de falhas em enlaces de longa distância. As simulações mostraram também que em ambientes mais complexos, com muitos clientes, e onde a restrição de banda em enlaces de longa distância provoca congestionamento e descartes, o desempenho do protocolo BitTorrent padrão pode ser semelhante a uma distribuição em arquitetura cliente-servidor. Neste último caso, a modificação proposta mostrou resultados consistentes de melhoria do desempenho da distribuição.

Morfologia e topografia de superfícies de aço inoxidável duplex UNS S32205, submetidas à erosão por impactos de partículas de alumina em fluxo de ar.

Os aços inoxidáveis duplex vêm sendo cada vez mais utilizados em aplicações que demandam maior resistência à corrosão que os aços inoxidáveis tradicionais. Suas aplicações em estruturas ao ar livre, ou submersas, motivam estudos sobre seu desempenho em variadas situações, como no caso de exposição ao desgaste erosivo.O presente trabalho estuda os efeitos da erosão por impacto de partículas sólidas em superfícies de amostras de um aço inoxidável duplex UNS S32205, tanto do ponto de vista morfológico quanto do topográfico. Superfícies de amostras polidas e posteriormente submetidas a ataques erosivos com partículas de alumina em fluxo de ar, foram examinadas por microscopia eletrônica de varredura, para verificação das formas de desgaste obtidas, e submetidas a ensaios de rugosidade, para determinação do perfil e avaliação da evolução do processo erosivo. Foi feito um exame mais detalhado de formatos de crateras, para comparação com mecanismos propostos em modelos descritos na literatura e com resultados de simulações pelo método dos elementos finitos. Nesse aspecto, o método utilizado mostrou-se tão eficaz quanto aos métodos tradicionais de estudos balísticos com uma só partícula, além de mais abrangente e econômico. As medidas de rugosidade superficial e as informações estatísticas derivadas sinalizaram a possibilidade de utilização dessa ferramenta para avaliar as condições de efetiva perda de massa, que caracterizam a erosão.

Avaliação de modelos de permeabilidade em meios porosos não consolidados

As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais.

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.

Controle a estrutura variável para atenuar modos ressonantes em servomecanismos.

Nesta dissertação é mostrada experimentalmente a existência de modos ressonantes em servomotores de corrente contínua sem núcleo de ferro. Até onde se sabe, não já publicações ou documentos técnicos de fabricantes que relatam esse efeito neste tipo de motor. Foi observado também que a frequência da ressonância é variável e que possivelmente depende da temperatura interna do motor. Conforme e mostrado em simulações técnicas tradicionais para o controle de motores podem não ser eficazes para lidar com esses modos ressonantes. Com o objetivo de amortecer a ressonância e obter o controle de velocidade e posição desses motores, são desenvolvidos controladores por modo deslizante baseados em funções de chaveamento obtidas por observadores de alto ganho. Esse controle e robusto a incertezas paramétricas e perturbações. Com o auxílio de simulações são mostrados os pontos fortes e limitações de cada um dos controladores.

Modelagem Numérica da Construção de Aterro Instrumentado na Baixada Fluminense, Rio de Janeiro.

A complexidade na previsão do comportamento de aterros sobre solos compressíveis representa um desafio para os engenheiros geotécnicos, mesmo os mais experientes. As limitações das teorias de análise e a dificuldade na seleção de parâmetros muitas vezes acarretam em estimativas de recalque incoerentes com a observação de campo. O presente trabalho tem como objetivo reproduzir a seqüência construtiva de extenso aterro, instrumentado, executado na Baixada Fluminense, Rio de Janeiro, para implantação da indústria Rio Polímeros, utilizando o programa PLAXIS. Os parâmetros geotécnicos foram definidos a partir de ensaios oedométricos, tendo sido realizado um estudo prévio sobre a qualidade das amostras. Estes parâmetros foram confrontados com a experiência local, que resume mais de 30 anos de investigações na argila mole da Baixada Fluminense, no Rio de Janeiro. Apesar da distância entre os locais de investigações, não houve diferenças significativas entre os parâmetros geotécnicos. As análises numéricas foram realizadas considerando-se diferentes hipóteses para o coeficiente de permeabilidade, de forma a reproduzir adequadamente as leituras de placa de recalque e poropressão registradas no campo. Este parâmetro teve grande influência na resposta da modelagem numérica, tendo sido obtidos bons resultados com valores de permeabilidade variando em função da profundidade.

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.