Estudo, proposta e previsão de desempenho de um filtro de energia para experimentos com alta resolução de comprimento de onda em fontes pulsantes de nêutrons

Em fontes de nêutrons por spallation para pesquisa de materiais, o comprimento de onda dos nêutrons é geralmente determinado pelos tempos de percurso (TOF) dos nêutrons desde a fonte até o detector. A precisão atingível é limitada pelo fato de o tempo de emissão do pulso característico do sistema alvo/moderador ser diferente de zero (a situação ideal fictícia seria a emissão se ocorresse na forma de um impulso). “Moderadores acoplados” (elementos usados para produzir feixes de alta intensidade com nêutrons de baixa energia) apresentam um decaimento de intensidade em função do tempo muito longo, ao longo de todo o espectro usado nos experimentos. Por este motivo, “moderadores desacoplados”, os quais produzem feixes com intensidade mais reduzida, são freqüentemente usados para instrumentos que requerem alta resolução. Neste trabalho, propusemos e analisamos uma nova técnica de filtragem dinâmica de feixes de nêutrons polarizados de baixa energia para experimentos que utilizam TOF na determinação do comprimento de onda. O dispositivo consiste de um sistema ótico polarizador/analisador e um inversor de spin seletivo em energia, o qual funciona por ressonância espacial do spin. Variando a condição de ressonância em sincronia com a estrutura temporal do pulso de nêutrons (através do controle de campos magnéticos), o filtro pode ajustar a resolução de energia (ou de comprimento de onda) de pulsos de banda larga em tais experimentos, separando os nêutrons com a correta relação “TOF/comprimento de onda” dos demais Um método para o cálculo de desempenho do sistema foi apresentado em um Trabalho Individual (TI) (PARIZZI et al., 2002 - i), do qual se fará aqui uma breve revisão além de amplo uso na otimização dos parâmetros a serem ajustados para o projeto do filtro. Os resultados finais mostram que ganhos consideráveis em resolução podem ser obtidos com a aplicação desta técnica em experimentos de reflectometria por tempo de percurso, sem que para tal seja necessário comprometer a intensidade do feixe usado pelo mesmo instrumento quando operado em um modo de baixa resolução, dando ao usuário do instrumento a opção de escolher a relação ótima entre intensidade e resolução para seu experimento. Como parte da conclusão desta dissertação, é apresentada uma proposta de parâmetros para a construção deste tipo de filtro e previsão de desempenho da configuração proposta, baseada no software de modelamento desenvolvido.

Interpretação integrada de sísmica de alta-resolução e da morfologia submarina da costa de fiordes na patagônia central - Chile

Este estudo apresenta o uso integrado da sísmica de alta resolução e da morfologia submarina para interpretar a evolução do ambiente glacimarinho da costa de fiordes da Patagônia Central, Chile. Foram analisados registros de perfilador de fundo e subfundo 3,5 kHz e modelos submarinos 3D de sete fiordes adjacentes ao campo de gelo Patagônico Sul (Eyre, Falcon, Penguin, Europa Peel, Calvo e Amalia) e parte do canal Icy. Os registros, com cerca de 300 km de levantamento acústico, foram obtidos pelo Servicio Hidrográfico y Oceanográfico de la Armada de Chile (SHOA), durante o cruzeiro de investigação científica Campo de Hielo Sur, realizado em 1995. Foram identificadas as principais fácies acústicas e geoformas sedimentares. A morfologia submarina e das bacias subaéreas adjacentes foi analisada com a elaboração de modelos tridimensionais subaéreos e submarinos e de perfis batimétricos longitudinais a partir de cartas náuticas, também do SHOA. Foram utilizadas imagens Landsat ETM+ na interpretação da geomorfologia glacial da área de entorno subaérea dos fiordes. O eixo longitudinal dos fiordes exibe morfologia irregular com bacias profundas e mostra fácies acústicas associadas ao sistema de depósitos de zonas de linha de encalhe (grounding line), às línguas de gelos flutuantes, e aos icebergs e ao gelo marinho. Refletores acústicos distinguem duas fácies principais, segundo sua configuração interna e geometria externa: caóticas e estratificadas. A geometria dos depocentros e as características dos refletores acústicos indicam a importante influência da batimetria e da topografia pré-existentes na dinâmica das geleiras e nos conseqüentes processos de sedimentação. Devido às grandes profundidades das bacias, as frentes das geleiras poderiam ser flutuantes ou aterradas ao fundo marinho. Em todo o caso, predomina um regime glacial onde o gelo está perto do ponto de fusão. É sugerido que os fiordes estudados resultam de um continuum de formas e de processos, que vão dos lineamentos e sistemas de falhas pré-existentes e controlados tectonicamente desde o Mioceno Inferior (ca. 25 milhões de anos AP), processos de denudação fluvial e de vertentes, vales fluviais modificados glacialmente e, finalmente, canais e fiordes erodidos pela ação do gelo. A tectônica originou uma topografia favorável para o desenvolvimento das geleiras, ocorreram processos de retroalimentação positiva entre aumento o da precipitação e expansão do campo de gelo no lado oeste da cordilheira. Provavelmente, as geleiras durante o Último Máximo Glacial não foram suficientemente espessas para aprofundar os fiordes, tampouco para dragar seus depósitos sedimentares para fora deles. Os processos de deposição de sedimentos caóticos acusticamente visíveis ocorreram durante o recuo das geleiras, já no Holoceno, onde elas alcançaram pontos de estabilidade, ajudadas pela morfologia dos fiordes, mesmo em águas profundas. Os depósitos estratificados localizados nas bacias intra-sills e em maiores profundidades, por sua vez, tiveram sua origem, provavelmente, das partes flutuantes do gelo, onde predominam os processos de desprendimento de icebergs (calving) e também do gelo marinho. Nessas mesmas bacias, depósitos mais profundos, não visíveis ao sistema de alta-resolução, poderiam estar preservados no fundo, pois não foram erodidos pelos sucessivos ciclos de avanço e recuo das geleiras, contendo assim informações sobre a evolução do campo de gelo ao longo do Quaternário.

Metodologia de análise de dioxinas em cinzas obtidas da incineração de couro Wet-Blue

O presente trabalho trata da validação de uma metodologia de análise para a determinação de dibenzo-p-dioxinas policloradas (PCDDs) em cinzas obtidas através de um processo de incineração de aparas de couro wet-blue. As dioxinas, compostos de extrema toxicidade, podem ser formadas na superfície das partículas de cinzas em reações de combustão. O mecanismo de formação de PCDDs requer uma atmosfera rica em oxigênio, uma fonte de cloro e a presença na cinza de um metal capaz de catalizar uma reação Deacon. O couro wet-blue apresenta em sua composição os precursores para a formação de dioxinas tais como o cloro e o metal catalisador sódio, cuja presença destes elementos foi comprovada por análises de superfície realizadas no couro e nas cinzas usando as técnicas de Rutherford Back-scattering Spectrometry (RBS) e Microscopia Eletrônica de Varredura (MEV). Desta forma, torna-se necessária a análise de dioxinas em cinzas de couro para especificar seu tratamento e destino final. As cinzas de couro foram obtidas numa unidade de bancada de incineração com reator de leito fixo, processo conhecido como GCC (Gaseificação e Combustão Combinadas). A metodologia utilizada para analisar dioxinas nas cinzas de couro foi baseada no Método 8280B da Agência de Proteção Ambiental dos Estados Unidos da América (EPA-USA). Este método descreve os procedimentos de extração de matrizes específicas, o clean-up para os analitos específicos e a determinação de dioxinas através das técnicas de cromatografia gasosa de alta resolução e espectrometria de massa de baixa resolução. O sucesso da análise depende da preparação dos extratos das amostras, etapa conhecida como clean-up, a qual é descrita minuciosamente neste trabalho, incluindo a montagem das colunas cromatográficas. A eficiência de extração das dioxinas foi obtida com os resultados comparativos das análises entre matrizes contaminadas com solução padrão de dioxinas e o padrão analítico. A validação deste método para análise de dioxinas policloradas em cinzas de couro wet-blue foi considerada satisfatória para os objetivos do trabalho, uma vez que a eficiência de extração de dioxinas obtida ficou dentro dos critérios estabelecidos pelo método.

Estudo da difusão de Ag e Al em [alfa]-Ti utilizando as técnicas de espectrometria de retroespalhamento Rutherford e reação nuclear

O objetivo da presente dissertação é o estudo da difusão de Ag e Al em uma matriz de α-Ti. Sua motivação principal se origina do fato de haver na literatura resultados contraditórios sobre o comportamento difusional desses elementos. Além disso, este estudo é necessário para dar continuidade à investigação sistemática da difusão de impurezas substitucionais em α-Ti, que vem sendo realizada pelo grupo de implantação iônica, a fim de estabelecer uma relação entre tamanho, valência e solubilidade dos elementos e seus coeficientes de difusão. A dependência do coeficiente de difusão com a temperatura para ambos elementos foi estudada nos intervalos de temperatura de 823 a 1073 K para a Ag e 948 a 1073 K para o Al. No caso da Ag foi usada a técnica de RBS para determinar os perfis de concentração, a qual tem uma alta resolução em profundidade (tipicamente 10nm), o que permite a determinação das baixas difusividades esperadas no presente experimento ( D ≤ 10-7 m2/s). No caso do Al foi usada a técnica de NRA, que preenche os mesmos requisitos citados anteriormente, com uma resolução em profundidade de aproximadamente 10 Å. As medidas realizadas mostraram que, para ambos os casos, os coeficientes de difusão seguem uma representação linear de Arrhenius. Foram encontrados os seguintes parâmetros característicos de difusão: Do = (1 ± 0,75) x10-4 m2s-1 e Q = 279 ± 6 kJ/mol para a Ag e Do = (1,4 ± 1,2) x10-2 m2s-1 e Q = 326 ± 10 kJ/mol para o Al, os quais são típicos de um comportamento substitucional normal. A comparação desse trabalho com trabalhos prévios não mostra evidência de relacionamento entre tamanho, valência e solubilidade dos elementos e seus coeficientes de difusão.

Estudo de difusão de In e Pd em [alfa]-Ti utilizando a técnica de espectrometria de retroespalhamento de Rutherford e canalização

O objetivo da presente dissertação é o estudo de difusão de In e Pd na matriz de a-Ti. O interesse deste estudo baseia-se na investigação sistemática, realizada pelo grupo de Implantação Iônica, de difusão de impurezas substitucionais em a-Ti com o fim de determinar uma relação entre tamanho, valência e solubilidade e os respectivos coeficientes de difusão. Considerações de tamanho, diferença de valência e solubilidade indicam que o In deverá difundir substitucionalmente. Por outro lado, não está claro qual será o mecanismo de difusão do Pd, pois a sua solubilidade não é muito alta e seu raio atômico é menor que o da matriz. A dependência do coeficiente de difusão desses dois elementos com a temperatura foi estudada entre 823 e 1073 K no caso do In, e entre 673 e 1073 K no caso do Pd. Em ambos os casos, a técnica de RBS foi utilizada para determinar os perfis de concentração. Esta técnica tem uma alta resolução em profundidade (tipicamente 10 nm), o que permite a determinação das baixas difusividades esperadas no presente experimento (D ≤ 10-17 m² s-¹). Adicionalmente, utilizamos para o caso do Pd a técnica de canalização, com a finalidade de determinar o grau de substitucionalidade dos átomos de Pd, tanto imediatamente após a implantação, quanto depois de efetuados os recozimentos e, portanto, determinar o mecanismo de difusão do mesmo em α - Ti. No caso do In, os coeficientes de difusão seguem uma lei de Arrhenius com parâmetros de difusão Q = (260 ± 10) kJ/mol e Do = (2,0 ± 1,3) x 10-6 m² s-¹ , valores estes, típicos de um comportamento difusivo substitucional. No entanto, para o Pd, os coeficientes de difusão também seguem um comportamento tipo Arrhenius, com Q = (264 ± 9) kJ/mol e Do = (2,8 ± 1,3) x 10-3 m² s-¹. Quando se comparam os resultados do Pd com outros substitucionais, tais como Pb, Au, Sn e In, os resultados indicam que o Pd é um difusor rápido, porém mais lento que os intersticiais do tipo Fe, Co ou Cu. Por outro lado, os experimentos de canalização indicam que no intervalo de temperatura estudado, 30% dos átomos de Pd ficam em posição intersticial na rede. Isto é um forte indício de que o mecanismo de difusão é do tipo misto.

Densitovolumetria pulmonar em crianças com bronquiolite obliterante pós-infecciosa : correlação do volume de zonas hipoatenuadas com função pulmonar e teste da caminhada de seis minutos

Objetivo: Avaliar a correlação do volume de pulmão com densidade anormal através da Densitovolumetria (DV) com características clínicas e testes de função pulmonar em crianças com Bronquiolite Obliterante (BO). Métodos: Realizou-se um estudo transversal em 19 crianças, com idade entre 7 e 15 anos, com diagnóstico clínico-evolutivo e tomografia de tórax característicos de BO. Foram excluídas ou-tras doenças que cursam com obstrução ao fluxo aéreo. Todas as crianças fizeram o Teste da cami-nhada de seis minutos com monitorização da saturação de oxigênio da hemoglobina. A espirometria foi feita em repouso, após o Teste da caminhada e após a administração de broncodilatador. A Den-sitovolumetria foi realizada em um tomógrafo computadorizado helicoidal de pista simples, marca Toshiba, modelo Xvision EX, com pós-processamento de imagem em estação de trabalho O2 da Si-licon Graphics, com programa de computação Alatoview®. Cada paciente foi submetido a 9 aquisi-ções tomográficas eqüidistantes de alta resolução e a duas varreduras helicoidais, cobrindo toda a extensão do tórax em pausa da respiração no final da inspiração e expiração profundas. Para separar parênquima normal de parênquima com diminuição de atenuação utilizou-se o limiar -950 UH e um outro limiar escolhido de forma subjetiva pelo radiologista Resultados: O volume de parênquima pulmonar com densidade anormalmente baixa na inspiração variou de 0,03 a 8,67 % e, na expiração, de zero a 7,27% do volume pulmonar total. O volume de pulmão hipoatenuado teve boa correlação com os testes de função pulmonar; na inspiração com VEF1% (r= -0,56) e com VEF1/CVF%(r= -0,75). O percentual de zonas hipoatenuadas na inspiração apresentou correlação com VEF1%(r= -0,64) e com VEF1/CVF%(r= -0,71). Na expiração, houve correlação com VEF1% (r= -0,50) e não houve com o VEF1/CVF%. Na comparação com o previsto em adultos, a correlação foi melhor na inspiração. A saturação de oxigênio em repouso não apresentou associação com o volume de áreas hipoatenuadas, enquanto que a saturação mínima durante o exercício apresentou correlação negativa forte com o volume de zonas com hipoatenuação na inspiração (r= -0,60) e na expiração (r= -0,61). Os pacientes com volumes pulmonares maiores percorreram uma distância maior (r=0,53), e a distância percorrida não foi afetada significativamente pelo volume de áreas hipoatenuadas. Conclusão: Em crianças com BO, o volume de zonas hipoatenuadas correlaciona-se com o VEF1% e VEF1/CVF% e com a queda na saturação durante o exercício, mas não se correlaciona com a saturação em repouso e a distância percorrida.

Estimativa da área impermeável dentro da bacia hidrográfica do Arroio Diluvio (porto Alegre/RS) através de técnicas de sensoriamento remoto e geoprocessamento

Neste trabalho foram realizadas classificações utilizando-se as bandas 1 a 5 e 7 dos sensores Landsat 5 TM (1987) e Landsat 7 ETM+ (2000). A caracterização espectral dos materiais foi realizada em laboratório utilizando um espectrorradiômetro, e através das bandas 1 a 5 e 7 dos sensores Landsat 5 TM (1987) e Landsat 7 ETM+ (2000). A transformação dos dados multiespectrais de imagens de sensoriamento remoto é uma maneira de reduzir o volume de dados através da identificação de classes de interesse numa imagem digital. No intuito de verificar condições de melhoramento na classificação de alvos urbanos em imagens digitais, identificados por procedimentos já conhecidos, como a classificação pela Máxima Verossimilhança, escolheu-se um classificador baseado na lógica fuzzy. O classificador utilizado foi o Fuzzy Set Membership classification - Fuzclass, que faz parte de um conjunto de classificadores não-rígidos disponíveis no programa Idrisi 32. Uma vez que informações sobre o desempenho de produtos deste classificador em áreas urbanas são escassas, foram conduzidos ensaios de comparação de resultados obtidos por este classificador com a verdade terrestre, representada por uma imagem de alta resolução espacial do satélite QuickBird. As áreas teste selecionadas desta imagem atendem ao critério de inalterância das condições de ocupação para o intervalo temporal considerado A comparação feita, permite concluir que o classificador apresenta limitações na classificação de áreas urbanas devido ao comportamento espectral semelhante dos materiais que fazem parte dessa cobertura. A utilização de uma classe única para identificar áreas impermeáveis foi a solução adotada para contornar este óbice. O emprego de áreas teste possibilitou acertar a escolha do grau de possibilidade de presença da classe no pixel (PPCP). Uma comparação entre os resultados apresentados na classificação de áreas impermeáveis, com base nos classificadores Máxima Verossimilhança e Fuzclass, demonstrou um desempenho melhor do classificador fuzzy, em função do nível de PPCP ajustado durante a análise comparativa Landsat e Quickbird nas áreas teste. Um procedimento alternativo de estimativa de áreas impermeáveis em bacias urbanas é apresentado no final.

Desenvolvimento de métodos analíticos para determinação de níquel e vanádio em petróleo por espectrometria de absorção atômica em forno de grafite

Neste trabalho foram investigados diferentes métodos de preparação de amostra para a determinação de níquel e vanádio em óleo cru usando a espectrometria de absorção atômica com forno de grafite (GF AAS). Investigou-se o uso de xileno como diluente, o qual foi inadequado devido à baixa estabilidade de níquel nas soluções e resultando na contaminação do equipamento. As soluções diluídas em isobutilmetilcetona apresentaram melhor estabilidade, mas não houve reprodutibilidade das medidas com os padrões orgânicos de níquel. Melhores resultados foram obtidos com emulsões óleo-em-água, usando Triton X-100 como surfactante. Também foi estudado o uso de padrões aquosos (inorgânicos) para o preparo das emulsões. Após o estabelecimento de curvas de pirólise para padrão e amostra de óleo cru, observou-se diferença no comportamento térmico do níquel na matriz de óleo cru. Foi então realizada uma investigação usando um espectrômetro de absorção atômica com fonte contínua de alta resolução (HR-CS AAS), equipamento com elevada capacidade de correção de fundo, e constatou-se que cerca de 50 % do analito é perdido em temperaturas abaixo de 400 ºC, enquanto que o restante é termicamente estável até pelo menos 1200 ºC. A fim de resolver o problema de perda de níquel volátil, foi investigado o uso de paládio como modificador químico. Através de curvas de pirólise com e sem paládio, observou-se que a mesma perda também ocorreu para o vanádio. Após a otimização da melhor quantidade de paládio, verificou-se que foi necessária uma quantidade de 20 µg de paládio para estabilização de níquel e vanádio. A determinação de níquel e vanádio total em amostras de óleo cru foi realizada usando-se paládio como modificador químico. A especiação destes elementos foi possível através da determinação das espécies estáveis, provavelmente compostos salinos não porfirínicos de níquel e vanádio, sem a adição de paládio, e a fração de compostos voláteis (provavelmente porfirinas de níquel e vanádio) foi calculada por diferença. Estabeleceram-se as figuras de mérito obtendo-se valores de massa característica de 19 e 33 pg e limite de detecção de 43 e 113 pg para níquel e vanádio, respectivamente. Esses resultados são compatíveis com dados reportados na literatura. A exatidão deste procedimento foi verificada pela análise dos materiais de referência certificados de metais traço em óleo residual (SRM 1634c) para níquel e vanádio e o petróleo (RM 8505) somente para vanádio e também por comparação com os resultados obtidos por HR-CS AAS para várias amostras de óleo cru. A aplicação do teste estatístico t-student aos resultados de níquel e vanádio total obtidos por GF AAS e HR-CS AAS mostrou que os mesmos não são significativamente diferentes a um nível de 95% de confiança.

Ambiente de alto desempenho com alta exatidão para a resolução de problemas

Este trabalho visa a disponibilização de um ambiente de alto desempenho, do tipo cluster de computadores, com alta exatidão, obtida através da utilização da biblioteca C–XSC. A alta exatidão na solução de um problema é obtida através da realização de cálculos intermediários sem arredondamentos como se fossem em precisão infinita. Ao final do cálculo, o resultado deve ser representado na máquina. O resultado exato real e o resultado representado diferem apenas por um único arredondamento. Esses cálculos em alta exatidão devem estar disponíveis para algumas operações aritméticas básicas, em especial as que possibilitam a realização de somatório e de produto escalar. Com isso, deseja-se utilizar o alto desempenho através de um ambiente de cluster onde se tem vários nodos executando tarefas ou cálculos. A comunicação será realizada por troca de mensagens usando a biblioteca de comunicação MPI. Para se obter a alta exatidão neste tipo de ambiente, extensões ou adaptações nos programas paralelos tiveram que ser disponibilizadas para garantir que a qualidade do resultado final realizado em um cluster, onde vários nodos colaboram para o resultado final do cálculo, mantivesse a mesma qualidade do resultado que é obtido em uma única máquina (ou nodo) de um ambiente de alta exatidão. Para validar o ambiente proposto foram realizados testes básicos abordando o cálculo do produto escalar, a multiplicação entre matrizes, a implementação de solvers intervalares para matrizes densas e bandas e a implementação de alguns métodos numéricos para a resolução de sistemas de equações lineares com a característica da alta exatidão. Destes testes foram realizadas análises e comparações a respeito do desempenho e da exatidão obtidos com e sem o uso da biblioteca C–XSC, tanto em programas seqüenciais como em programas paralelos. Com a conseqüente implementação dessas rotinas e métodos será aberto um vasto campo de pesquisa no que se refere ao estudo de aplicações reais de grande porte que necessitem durante a sua resolução (ou em parte dela) da realização de operações aritméticas com uma exatidão melhor do que a obtida usualmente pelas ferramentas computacionais tradicionais.

Processamento de materiais carbonáceos por pulsos intensos de laser em alta pressão

Neste trabalho foi desenvolvida uma câmara de alta pressão com janela de safira para processamento de filmes finos com pulsos de laser de alta potência, num regime de resfriamento ultra-rápido e geometria confinada. As amostras estudadas consistiram de filmes finos de carbono amorfo depositados sobre substratos de cobre. Os processamentos foram realizados com um laser pulsado Nd:YAG com energia de até 500 mJ por pulso, com duração de 8 ns, focalizada numa região de cerca de 1,5mm2, gerando uma região de elevada temperatura na superfície da amostra durante um intervalo de tempo bastante curto, da ordem do tempo de duração do pulso do laser. Para evitar a evaporação do filme de carbono, aplicava-se através da câmara, uma pressão de 0,5 a 1,0 GPa, confinando a amostra e eliminando o efeito da ablação. Este sistema tornou possível produzir taxas de resfriamento extremamente elevadas, com supressão da formação de uma pluma durante a incidência do laser, sendo o calor dissipado rapidamente pelo contato com os substratos de cobre e safira, ambos com elevada condutividade térmica. As amostras processadas foram analisadas por micro-espectroscopia Raman e os resultados revelaram a formação de estruturas com cadeias lineares de carbono, “carbynes”, caracterizadas pela presença de um pico Raman intenso na região de 2150 cm-1. Outro conjunto de picos Raman foi observado em 996 cm-1, 1116 cm-1 e 1498 cm-1 quando o filme fino de carbono amorfo foi processado dentro da câmara, com uma seqüência de mais de três pulsos consecutivos de laser. Várias tentativas foram feitas para investigar a natureza da fase que origina estes picos. Apesar da similaridade com o espectro Raman correspondente ao poliacetileno (CnHn), não foi possível constatar evidências experimentais sobre a presença de hidrogênio nos filmes de carbono processados. Estes picos Raman não foram observados quando o filme de carbono era depositado sobre outros substratos metálicos, a não ser em cobre. O conjunto de resultados experimentais obtidos indica que estes picos estariam relacionados a pequenos aglomerados lineares de átomos de carbono, diluídos numa matriz de átomos de cobre, formados durante os pulsos subseqüentes de laser e retidos durante o resfriamento ultra-rápido da amostra. A comparação dos resultados experimentais com a simulação do espectro Raman para diferentes configurações, permite propor que estes aglomerados seriam pequenas cadeias lineares, com poucos átomos, estabilizadas frente à formação de grafenos pela presença de átomos de cobre em abundância. Foram também realizados processamentos de materiais carbonáceos por pulsos de laser em meios líquidos, através de câmaras especialmente construídas para este fim. Os resultados, em diversos materiais e configurações, mostraram apenas a formação de estruturas grafíticas, sem evidência de outras fases.

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.