Vibrações naturais em um sistema de interação viga-água incluindo o efeito de onda de superfície

Este trabalho tem como objetivo o estudo modal de um sistema estruturafluído, modelado pela equação de Euler-Bernoulli, para uma viga elástica, sujeita a pressão da água, e por ondas de superfície livre. O sistema acoplado possui condições para o domínio sólido (viga fixa-livre), para o domínio fluído (impermeabilidade e rigidez inferior), com ondas de superfície, e de interfacec fluido-estrutura (condições de continunidade na deflexão, ângulo de rotação, força de corte interno e momento de curvatura). Para deteerminar as vibrações livres e deslocamento no seco e no molhado, utiliza-se o método espectral para eliminar a dependência oscilatória temporal e concentrar-se na determinação dos modos através do estudo de problemas de contornos espaciais. Os modos podem ser calculados com o uso da base clássica de Euler ou da base dinâmica gerada pela resposta impulso. Foram feitas simulações para um material específico, e apresentados os resultados obtidos.

Características dinâmicas do escoamento turbulento em bancos de tubos inclinados

Bancos de tubos ou de barras são encontrados nas indústrias nucleares e de processos, constituindo-se nos dispositivos mais comumente utilizados no desenvolvimento de trocadores de calor. O procedimento usual na análise de bancos de tubos considera o escoamento incidindo uniformemente no banco ao longo de toda a extensão dos tubos. Em aplicações reais de trocadores de calor este fato não é observado, pois o escoamento não é uniforme devido a presença de defletores e condições de entrada e saída. A proposta deste trabalho é apresentar um estudo experimental das flutuações de pressão e velocidade e sua interdependência do escoamento turbulento incidindo em arranjos de cilindros inclinados, com diferentes ângulos de incidência (30º, 45º, 60º e 90º) simulando bancos de tubos inclinados com arranjo quadrangular e razão de espaçamento P/D=1,26. A seção de testes é um canal retangular com 146 mm de altura e largura de 193 mm. O fluido de trabalho é o ar, impulsionado por um ventilador centrífugo, passando por uma câmara de tranquilização e um conjunto de grades e telas, que tem a finalidade de uniformizar o escoamento e reduzir a intensidade de turbulência para a ordem de 2%, antes de alcançar a seção onde são realizadas as medições. Medições da distribuição de pressão na parede e flutuações de velocidade e pressão em torno de um dos tubos no interior dos bancos foram realizadas para diferentes ângulos de incidência, utilizando a técnica de anemometria de fio quente e um transdutor de pressão. Os resultados experimentais são apresentados na forma adimensional A distribuição média de pressão em uma das paredes laterais do canal, com a presença do banco de tubos, mostrou que as linhas isobáricas não são paralelas ao eixo dos tubos, indicando que o gradiente de pressão não é perpendicular aos tubos, nem paralelo a direção principal do escoamento. As médias quadráticas das flutuações de pressão indicam uma distribuição mais uniforme das cargas dinâmicas em torno dos tubos do bancos com inclinações de 30º e 90º. A análise dos resultados dos espectros das flutuações de pressão mostra que conforme o ângulo de incidência do escoamento aumenta, picos aparecem para freqüência adimensionais cada vez menores. O comportamento dos espectros de flutuação de velocidade mostram um decaimento suave das curvas quando o ângulo de incidência do escoamento aumenta. A correlação cruzada entre as flutuações de velocidade e pressão indicam que, quando o ângulo de incidência do escoamento diminui, importantes efeitos tridimensionais aparecem. Esses efeitos, caracterizados pelas oscilações na correlação cruzada para o ângulo de incidência de 45º, parecem desaparecer se o ângulo de incidência continua-se a ser reduzido.

Análise teórico experimental dos evaporadores de uma bomba de calor

Faz-se um estudo da construção de uma bomba de calor residencial de baixo custo utilizando um condicionador de ar do tipo janela. Com tal equipamento operando fora das condições de projeto é preciso fazer uma avaliação das novas condições de operação, onde os evaporadores passam a ser objetivo de estudo do presente trabalho. A análise foi conduzida de duas formas. Uma tomada de dados experimental, feita em duas épocas distintas do ano, nas condições mais severas de consumo, no inverno, e nas condições mais severas de temperatura, no verão. A outra, uma análise teórica, gerando dados que podem ser comparados com a literatura. Inicialmente, é feita uma avaliação do coeficiente de performance visando a competitividade do equipamento no que tange o custo operacional. Como objetivo principal é avaliada a troca de calor nos evaporadores e os parâmetros que a influenciam, bem como o comportamento do fluido refrigerante na entra e saída. Foi preciso desenvolver uma metodologia para trocadores aletados que apresentam mudança de fase já que esta não é vista com facilidade na literatura. Faz-se uso do método NUT e- , que é adaptado de forma a ser função somente dos dados psicrométricos do ar. Para a efetividade média externa e o coeficiente de convecção externo é utilizada a metodologia desenvolvida por Schmidt em 1949 Como resultados, observou-se que o coeficiente de performance fica na faixa dos equipamentos comercias (2,6 na média anual) e o custo operacional é competitivo com o uso do gás liquefeito de petróleo, em aquecedores de passagem. A efetividade do trocador e o coeficiente de convecção externo apresentam um valor baixo, atribuído a grande área de troca de calor externa. O coeficiente de contato, apresenta um valor alto se comparado ao coeficiente de convecção externo. Quanto ao fluido refrigerante, este apresenta condições normais nas saídas dos evaporadores e título baixo na entrada de um dos evaporadores. Também é observado que um dos evaporadores apresenta baixo uso efetivo, sendo o problema atribuído a baixa vazão de refrigerante.

Análise numérica de escoamentos turbulentos tridimensionais empregando o método de elementos finitos e simulação de grandes escalas

O presente trabalho apresenta o estudo e implementação de um algoritmo numérico para análise de escoamentos turbulentos, tridimensionais, transientes, incompressíveis e isotérmicos, através da Simulação de Grande Escalas, empregando o Método de Elementos Finitos. A modelagem matemática do problema baseia-se nas equações de conservação de massa e quantidade de movimento de um fluido quase-incompressível. Adota-se um esquema de Taylor-Galerkin, com integração reduzida e fórmulas analíticas das funções de interpolação, para o elemento hexaédrico de oito nós, com funções lineares para as componentes de velocidade e constante no elemento para a pressão. Para abordar o problema da turbulência, emprega-se a Simulação de Grandes Escalas, com modelo para escalas inferiores à resolução da malha. Foram implementados o modelo clássico de Smagorinsky e o modelo dinâmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, é proposta e empregada, para o processo de segunda filtragem do modelo dinâmico. O esquema, que utiliza elementos finitos independentes envolvendo cada nó da malha original, apresentou bons resultados com um baixo custo computacional adicional. São apresentados resultados para problemas clássicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para análise de escoamentos caracterizados por elevados números de Reynolds, é discutida no capítulo final. São apresentadas sugestões para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para análise de escoamentos tridimensionais, turbulentos e transientes .

Simulação numérica de esteiras em transição utilizando o método dos contornos virtuais

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Modelo matemático de um agitador eletromagnético

O objetivo deste trabalho é estudar os efeitos eletromagnéticos e fluido-dinâmicos induzidos no aço, decorrentes do uso de um agitador eletromagnético. Para tal, foi proposta a construção de um modelo numérico que resolva, de forma acoplada, os problemas de eletromagnetismo e fluido-dinâmica. O modelo numérico do problema eletromagnético, em elementos finitos, foi construído utilizando-se o software Opera-3d/Elektra da Vector Fields. O mesmo foi validado com medidas experimentais de densidade de fluxo magnético feitas na usina. O escoamento decorrente da agitação eletromagnética foi resolvido fazendo-se o acoplamento das forças de Lorentz com as equações de Navier-Stokes. Essas últimas foram resolvidas pelo método de volumes finitos, usando-se o software CFX-4 da AEA Technology. O modelo eletromagnético mostrou que existe um torque máximo dependente da freqüência do campo magnético. Também foi observado que a força magnética aumenta em quatro vezes seu valor, quando a corrente é duplicada. O perfil de escoamento produzido no molde, sob agitação eletromagnética, indica, que as situações de lingotamento testadas, não propiciam o arraste da escória. A velocidade crítica de arraste, determinada via modelo físico, não foi atingida para nenhum caso testado. O modelo fluido-dinâmico e térmico apresentou um aumento do fluxo de calor cedido pelo fluido para a casca solidificada com o uso do agitador eletromagnético. Como conseqüência, observou-se uma queda na temperatura do banho. Também foi observado, que o uso do agitador propicia a remoção de inclusões das camadas mais externas do tarugo. Ao mesmo tempo, notou-se que o uso do agitador aumenta o índice de remoção de inclusões para as duas seções de molde analisadas.

Simulação de escoamentos incompressíveis não newtonianos em dutos com expansão brusca

O objetivo deste trabalho é a simulação numérica de escoamentos incompressíveis bidimensionais em dutos com expansão brusca, considerando o raio de expansão de 3 : 1. As equações governantes implementadas são as de Navier, que junto com relações constitutivas para a tensão visam representar comportamentos não newtonianos. A integração temporal é feita usando o esquema explícito de Runge-Kutta com três estágios e de segunda ordem; as derivadas espaciais são aproximadas pelo método de diferenças finitas centrais. Escoamentos em expansões bruscas para fluidos newtonianos apresentam um número de Reynolds crítico, dependente do raio de expansão, na qual três soluções passam a ser encontradas: uma solução sim étrica instável e duas soluções assimétricas rebatidas estáveis. Aumentando o número de Reynolds, a solução passa a ser tridimensional e dependente do tempo. Dessa forma, o objetivo é encontrar as diferenças que ocorrem no comportamento do fluxo quando o fluido utilizado possui características não newtonianas. As relações constitutivas empregadas pertencem à classe de fluidos newtonianos generalizados: power-law, Bingham e Herschel-Bulkley. Esses modelos prevêem comportamentos pseudoplásticos e dilatantes, plásticos e viscoplásticos, respectivamente. Os resultados numéricos mostram diferenças entre as soluções newtonianas e não newtonianas para Reynolds variando de 30 a 300. Os valores de Reynolds críticos para o modelo power-law não apresentaram grandes diferenças em comparação com os da solução newtoniana. Algumas variações foram percebidas nos perfis de velocidade. Entretanto, os resultados obtidos com os modelos de Bingham e Herschel-Bulkley apresentaram diferenças significativas quando comparados com os newtonianos com o aumento do parâmetro adimensional Bingham; à medida que Bingham é aumentado, o tamanho dos vórtices diminui. Além disso, os perfis de velocidade apresentam diferenças relevantes, uma vez que o fluxo possui regiões onde o fluido se comporta como sólido.

Estudo dos compostos voláteis do alecrim utilizando as técnicas de microextração em fase sólida (SPME), hidrodestilação e extração com fluído supercrítico (SFE)

Resumo não disponível.

Aspectos epidemiológicos da paratuberculose bovina no Rio Grande do Sul

Descrevem-se os achados clínicos e patológicos da paratuberculose em uma criação intensiva de bovinos de leite no municí pio de Capela de Santana, RS. Os sinais clínicos foram observados em oito de um total de 345 bovinos, consistindo em diarréia crônica refratária aos tratamentos, emagrecimento progressivo e queda na produção de leite. As principais lesões macroscópicas, observadas nos oito animais necropsiados, incluíam intestino delgado com acentuado espessamento da parede e superfície mucosa de aspecto reticulado, semelhante às circunvoluções cerebrais, lesão essa perceptível, através da serosa. A luz intestinal estava preenchida com conteúdo fluido e de aspecto leitoso. Os vasos linfáticos do mesentério mostravam-se mais evidentes sendo que alguns tinham aspecto varicoso. Os linfonodos mesentéricos estavam aumentados de volume e, ao corte, fluía grande quantidade de líquido leitoso. Focos de mineralização foram observados na íntima das artérias, nas válvulas cardíacas e na serosa do rúmen. As principais lesões macroscópicas incluíam enterite, linfadenite e linfangite granulomatosa caracterizada por infiltrado inflamatório com macrófagos, células gigantes de Langhans que continham grande quantidade de bacilos álcool-ácido-resistentes. As lesões vasculares consistiam em degeneração e mineralização das túnicas í ntimas e média das artérias de grande calibre associada a proliferação de colágeno. Havia calcificação da serosa do rúmen atrofia hepatocelular difusa e hepatite granulomatosa multifocal. O M. avium subsp. paratuberculosis (Map) foi isolado em amostras de intestino e linfonodos de 8 vacas Holandesas (3,5%) com doença de Johne, dentre 229 amostras cultivadas provenientes de um rebanho leiteiro. Amostras inoculadas em HEYM com micobactina produziram colônias identificadas como Map, segundo as caracterí sticas fenotí picas próprias como: crescimento lento, coloração álcool-ácido-resistente (A.A.R.) e dependência a micobactina. O laboratório de Referência da OIE confirmou a amostra isolada. O M. avium subsp. paratuberculosis (Map) foi isolado em amostras de intestino e linfonodos de 8 vacas Holandesas (3,5%) com doença de Johne, dentre 229 amostras cultivadas provenientes de um rebanho leiteiro Amostras inoculadas em HEYM com micobactina produziram colônias identificadas como Map, segundo as caracterí sticas fenotí picas próprias como: crescimento lento, coloração álcool-ácido-resistente (A.A.R.) e dependência a micobactina. O laboratório de Referência da OIE confirmou a amostra isolada. Não houve isolamento do agente em 221 amostras intestinais quando processadas, após 2 anos de sua colheita. O teste de IDGA aplicado como “screening”, detectou 26 vacas (11,4%) positivas, dentre 228 animais testados e sacrificados em matadouro. O teste de ELISA adsorvido, utilizando o antí geno PPA-3 detectou 125 (39,8%) amostras positivas. O ELISA não adsorvido detectou mais 32 (10,1%) reagentes positivos, dentre os 314 bovinos testados. A prevalência da infecção causada pelo Map em 36 rebanhos leiteiros procedentes de 25 municí pios do Rio Grande do Sul foi estimada em 44,6% das 1316 amostras testadas. A infecção foi identificada em 35 (97,2%) dos 36 rebanhos testados e presentes em todos os municí pios incluí dos. A ocorrência da doença de Johne foi enfatizada, tanto a forma clí nica quanto a infecção subclí nica no Rio Grande do Sul, sugerindo a adoção de medidas de controle sejam aplicadas na proteção dos rebanhos leiteiros nacionais.

Análise energética e exergética de um ciclo rankine com aquecimento distrital : estudo de uma planta termelétrica

Este trabalho tem por objetivo simular e analisar uma usina termelétrica a carvão em várias condições de funcionamento. A usina simulada neste trabalho é a AVV 1 localizada em Copenhague, Dinamarca. A AVV 1 é uma usina de geração de potência e aquecimento distrital que pode funcionar em várias condições de carga. A simulação da usina supracitada foi tema de um concurso de simuladores proposto no congresso ECOS 2003 realizado em Copenhague, Dinamarca. Para a realização deste trabalho foi construído um programa na linguagem FORTRAN 90. Cada componente da usina é modelado através de equações de balanço de massa e energia, e o sistema completo tem sua solução obtida pelo método de substituição sucessiva. Para viabilizar essa solução é necessário também implementar uma rotina de cálculo de propriedades do fluido de trabalho. No caso estudado, o fluido de trabalho da usina é a água e a formulação utilizada para o cálculo de suas propriedades nos diversos estados é a IAPWS IF-97. A usina é simulada em dois modos de operação: modo de condensação, onde ocorre apenas geração de eletricidade, e em modo de contrapressão, onde há geração de eletricidade e aquecimento distrital, conforme nomenclatura sugerida pela organização do concurso No modo de condensação, são feitas quatro séries de simulações variando a carga de 100% a 40%. Cada série contém um conjunto de hipóteses quanto à variação das eficiências isoentrópicas e pressões das turbinas em função da vazão mássica. No modo de contrapressão, a usina é simulada funcionando com 100% da carga. O programa desenvolvido calcula as propriedades para qualquer ponto de trabalho ao longo da planta, assim como a eficiência da mesma, a potência gerada, e todas as vazões mássicas pertinentes. Além disso, é feita também uma análise exergética da planta. A simulação demonstrou que a planta possui uma eficiência global de 42,02% com uma geração de 250,2 MW em 100% de carga no modo de condensação. Nessas mesmas condições, do ponto de vista exergético, a eficiência encontrada é de 37,21%. No modo de contrapressão, a usina apresenta uma eficiência exergética de 40,19% com um aproveitamento energético de 90,55%. Por fim, é possível também verificar a comportamento da eficiência da planta e a variação de água de resfriamento do condensador com a carga. Os resultados gerados são próximos àqueles encontrados pelos diversos pesquisadores que abordaram o problema.

Desenvolvimento de um programa de simulação computacional de sistemas de aquecimento solar para água

Esta Tese apresenta uma análise do comportamento térmico de um sistema de aquecimento solar operando por termossifão. Neste tipo de sistema o fluido no coletor solar é circulado por convecção natural, que acontece devido à diferença de massa específica da água ao longo circuito. Nestes sistemas a vazão mássica varia ao longo do dia e do ano, dependendo, dentre outros fatores, da irradiância solar absorvida, do perfil de temperaturas da água no sistema, da geometria, do volume e do perfil de demanda de água quente. Para uma avaliação detalhada do comportamento térmico de aquecedores solares operando por termossifão foram realizados ensaios experimentais e cálculos teóricos. Os resultados dos experimentos concordaram com aqueles apresentados na literatura e sua análise fundamentou o desenvolvimento do aplicativo TermoSim, um programa de simulação computacional do comportamento térmico de sistemas de aquecimento de água com energia solar. O tratamento matemático adotado no TermoSim compreende a modelagem dos coletores solares de acordo com a teoria de Hottel-Bliss-Whillier. O reservatório térmico é modelado com estratificação térmica, convecção e condução entre as camadas. A vazão mássica é obtida a partir do balanço da quantidade de movimento no circuito. Os modelos matemáticos empregados na construção do aplicativo TermoSim foram validados através do confronto dos resultados simulados com medidas experimentais. Foi demonstrado que a utilização destes modelos é adequada e permite reproduzir com precisão o comportamento térmico dos coletores solares e do reservatório térmico. Além do programa TermoSim, foi também desenvolvido o programa TermoDim, que é uma ferramenta para o dimensionamento de sistemas de aquecimento solar, que requer apenas o conhecimento dos parâmetros geométricos do sistema, dados meteorológicos em média mensal e informação a respeito do volume de demanda. O TermoDim é apropriado para estimar o desempenho de aquecedores solares operando por termossifão com tanques verticais e horizontais. O método de dimensionamento do TermoDim é baseado na correlação para a eficiência média mensal obtida neste trabalho a partir de um grande número de simulações.

Ferramenta computacional para a medição de campos de velocidade utilizando processamento digital de imagens

A proposta deste trabalho, consiste na elaboração de uma ferramenta computacional para a medição de campos de velocidades em escoamentos com baixas velocidades (< 0,5 m/s) utilizando o processamento digital de imagens. Ao longo dos anos, inúmeras técnicas foram desenvolvidas com este objetivo. Para cada tipo de aplicação, uma técnica se aplica com maior ou menor eficiência do que outras. Para o caso de estudos em fluídos transparentes, onde o escoamento pode ser visualizado, técnicas que utilizam processamento digital de imagens vêm ganhando um grande impulso tecnológico nos últimos anos. Este impulso, é devido a fatores como: câmaras vídeo filmadoras de última geração, dispositivos de aquisição de imagens e componentes de processamento e armazenamento de dados cada vez mais poderosos. Neste contexto, está a velocimetria por processamento de imagens de partículas cuja sigla é PIV (particle image velocimetry). Existem várias formas de se implementar um sistema do tipo PIV. As variantes dependem, basicamente, do equipamento utilizado. Para sua implementação é necessário, inicialmente, um sistema de iluminação que incide em partículas traçadoras adicionadas ao fluido em estudo. Após, as partículas em movimento são filmadas ou fotografadas e suas imagens adquiridas por um computador através de dispositivos de captura de imagens. As imagens das partículas são então processadas, para a obtenção dos vetores velocidade. Existem diferentes formas de processamento para a obtenção das velocidades. Para o trabalho em questão, devido às características dos equipamentos disponíveis, optou-se por uma metodologia de determinação da trajetória de partículas individuais, que, apesar de limitada em termos de módulo de velocidade, pode ser aplicada a muitos escoamentos reais sob condições controladas Para validar a ferramenta computacional desenvolvida, imagens ideais de partículas foram simuladas como se estivessem em escoamento, através do deslocamento conhecido de vários pixels. Seguindo o objetivo de validação, foi utilizada ainda uma imagem real de partículas, obtida com o auxílio de um plano de iluminação de luz coerente (LASER) e câmaras de vídeo tipo CCD. O programa desenvolvido foi aplicado em situações de escoamento real e os resultados obtidos foram satisfatórios dentro da escala de velocidades inicialmente presumida.

Processos e componentes mantélicos no norte da Patagônia (Argentina) e relações com a subducção Andina : Evidências petrográficas, geoquímicas e isotópicas em xenólitos ultramáficos mantélicos

Xenólitos ultramáficos, carregados até a superfície da Terra por magmatismo básico alcalino intraplaca, fornecem evidências diretas da natureza e processos envolvidos em modificações do manto litosférico subcontinental, como fusão parcial e metassomatismo. Estes xenólitos têm sido utilizados para identificar processos relacionados a evolução da litosfera continental, estimar a composição original do manto e a escala das heterogeneidades mantélicas. Raramente xenólitos ultramáficos mantélicos são encontrados em ambientes convergentes, no entanto na Patagônia (sul da América do Sul), diversas ocorrências são identificadas em basaltos alcalinos na região de arco e back arc da Cordilheira do Andes. Estes xenólitos oportunizam o estudo dos processos de interação entre a cunha mantélica, a placa oceânica subductada e a astenosfera. Nesta dissertação são apresentados dados petrográficos, mineralógicos, geoquímicos e de isótopos de Sr e O em 22 xenólitos ultramáficos de dois centros vulcânicos Mioceno-Holoceno distintos: Cerro del Mojon (41°06’S-70°13’W) e Estancia Alvarez (40°46’S-68°46’W), localizados na borda NW do Platô de Somuncura, norte de Patagônia (Argentina). A suíte de xenólitos ultramáficos do Cerro del Mojon consiste de espinélio dunitos e harzburgitos mantélicos anidros (Grupo1), espinélio lherzolitos mantélicos anidros (Grupo 2a) e hidratados (Grupo 2b) e espinélio clinopiroxenitos crustais (Grupo 3). Os xenólitos mantélicos do Grupo 1 são depletados (empobrecidos em ETR pesados e HFSE, com baixas razões 87Sr/86Sr em Cpx – 0,7028-0.7037), de alta PT (16-19 kbar, 950-1078 ºC), e têm evidências de metassomatismo críptico (enriquecimento em K, Na ETR leves) atribuído a componentes derivados de sedimentos da placa oceânica subductada, EM 2 (87Sr/86Sr em RT até 0,7126, 87Rb/86Sr até 1,66 e δ18O até +6.78‰). Estes valores anomalamente altos foram obtidos em amostras com bolsões de reação ao redor do espinélio, induzidos pela percolação de fluidos metassomáticos sob altas pressões seguida por descompressão (Sp+fluido→ Cpx+Ol+Sp+melt-andesítico-traquítico). Os xenólitos do Grupo 2 são moderadamente depletados (empobrecidos em HFSE, com baixas a altas razões 87Sr/86Sr em RT e Cpx – 0,7031-0,7045 e δ18O +5-6.2‰), de baixa PT (14-15 kbar, 936-942 ºC), e têm evidências de metassomatismo modal e críptico (enriquecimento em ETR leves, Na, K, Ti, Sr, Hf e anfibólio modal). O metassomatismo críptico parece estar relacionado ao mesmo agente metassomatizante do Grupo 1, mas ocorre em menor intensidade. O metassomatismo modal, no entanto, tem outra origem, podendo ser derivado de fontes mantélicas profundas. A quebra do anfibólio metassomático durante a descompressão e ascensão dos xenólitos até a superfície formou bolsões de reação ao redor do anfibólio (Anf→Cpx+Ol+Sp+melts-basálticos & andesíticos). A suíte de xenólitos ultramáficos de Estancia Alvarez consiste de espinélio harzburgitos anidros mantélicos (Grupo 1a) e espinélio dunitos anidros crustais (Grupo 1b), ambos com veios de serpentina. Os xenólitos mantélicos do Grupo 1a são depletados (empobrecidos em ETR leves e HFSE), de profundidades variadas (P entre 11-18 kbar) e de baixa T (877-961 ºC). Estes xenólitos têm enriquecimento em ETR leves, B, Rb e K, que pode estar relacionado a percolação de fluidos ricos em H2O e LILE gerados pela desidratação de filossilicatos dos sedimentos da placa oceânica (EM 2), como evidenciado por veios de serpentina e altas razões 87Sr/86Sr (0,7046-0,7298), 87Rb/86Sr (0,07-5,63) e valores de δ18O (até +7,88‰).

Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.