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 .

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.

Variação da estrutura e da biomassa fitoplanctônica na Lagoa Itapeva (Litoral Norte do Rio Grande do Sul) em função da hidrodinâmica

As lagoas costeiras em geral exibem reduzidas profundidades máximas, sendo portanto, ambientes altamente vulneráveis às oscilações dos fatores climáticos, em particular o vento. Movimentos da água em zonas pelágicas são fortemente influenciados por mudanças no tempo, as quais promovem mudanças no plâncton. A biomassa é uma importante variável em ecossistemas aquáticos uma vez que representa a quantidade de energia estocada em certo nível trófico da comunidade biológica. O trabalho teve como hipótese principal do trabalho que o vento (principal fator determinante na hidrodinâmica da Lagoa Itapeva), influencia diretamente a biomassa fitoplanctônica. Períodos de ausência de ventos ou ventos fracos provocam uma homogeneidade da comunidade, devido á estabilidade da coluna d’água, propiciando eventos de florações com alta biomassa e baixa diversidade. Outra hipótese formulada é que há variação espaço-temporal da biomassa, da diversidade de espécies e das diferentes frações de tamanho na Lagoa Itapeva, no período estudado. O objetivo do trabalho foi determinar o efeito do vento, principal fator hidrodinâmico da Lagoa Itapeva, sobre a comunidade fitoplanctônica em termos de biomassa, diversidade de espécies e categorias de tamanho. A Lagoa Itapeva é a primeira lagoa do Sistema Lagunar de Tramandaí no sentido N→S. Foram realizadas análises espaciais, formando-se três pontos de amostragem (Norte, Centro e Sul) ao longo do maior eixo longitudinal da lagoa, aliada a uma escala temporal diária (turnos de amostragem: 6h, 10h, 14h e 18h) e sazonal (quatro estações do ano), procurando estabelecer padrões e relações com os fatores hidrodinâmicos (velocidade e direção de vento, e nível d’água). Durante o estudo a Divisão Cyanophyta foi a responsável pelos elevados valores de biomassa observados nos pontos de amostragem. Este fato deve-se principalmente às florações de Anabaena circinalis durante as campanhas sazonais. A biomassa foi a variável que melhor se relacionou com os fatores físicos (nível d’água e velocidade média do vento) pela ação dos ventos de quadrantes predominantes que resultavam num fetch efetivo longo. Neste estudo foram apresentados índices de diversidade de espécies calculados através de indivíduos (H’ind) e de biomassa (H’ biom). Estes índices apresentaram estreitas relações positivas entre si. Os índices também apresentaram relações diretas e significativas com os fatores físicos (nível d’água e velocidade média do vento). Na Lagoa Itapeva, a contribuição das diferentes categorias de tamanho em função da biomassa (estimada através do biovolume), ressaltou três frações que melhor representaram a comunidade fitoplanctônica: nanoplâncton, microplâncton e netplâncton. A fração netplanctônica predominou sobre as demais frações de tamanho. A predominância desta fração sobre as outras estudadas deu-se principalmente às espécies de colônias filamentosas Aulacoseira granulata, Anabaena circinalis e A. spiroides. O comportamento da comunidade fitoplanctônica na Lagoa Itapeva sugere estar de acordo com a Hipótese de Distúrbio Intermediário de Connell (1978). Os distúrbios intensos (provocados pelos forte ventos e fetch longo) ou a falta deste (criando um ambiente estável), nos pontos de amostragem, resultaram em florações da cianobactéria Anabaena circinalis. Essas florações proveram uma densa biomassa em todo o período da pesquisa, além de promoverem a baixa diversidade através do domínio desta cianobactéria, e provocar sucessões das frações de tamanho do fitoplâncton.

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.

Uma ferramenta para automação da geração do leiaute de circuitos analógicos sobre uma matriz de transistores MOS pré-difundidos

Este trabalho apresenta o LIT, uma ferramenta de auxílio ao projeto de circuitos integrados analógicos que utiliza a técnica da associação trapezoidal de transistores (TAT) sobre uma matriz digital pré-difundida. A principal característica é a conversão de cada transistor simples de um circuito analógico em uma associação TAT equivalente, seguido da síntese automática do leiaute da associação séria-paralela de transistores. A ferramenta é baseada na matriz SOT (sea-of-transistors), cuja arquitetura é voltada para o projeto de circuitos digitais. A matriz é formada somente por transistores unitários de canal curto de dimensões fixas. Através da técnica TAT, entretanto, é possível criar associações série-paralelas cujo comportamento DC aproxima-se dos transistores de dimensões diferentes dos unitários. O LIT é capaz de gerar automaticamente o leiaute da matriz SOT e dos TATs, além de células analógicas básicas, como par diferencial e espelho de corrente, respeitando as regras de casamento de transistores. O cálculo dos TATs equivalentes também é realizado pela ferramenta. Ela permite a interação com o usuário no momento da escolha da melhor associação. Uma lista de possíveis associações é fornecida, cabendo ao projetista escolher a melhor. Além disso, foi incluído na ferramenta um ambiente gráfico para posicionamento das células sobre a matriz e um roteador global automático. Com isso, é possível realizar todo o fluxo de projeto de um circuito analógico com TATs dentro do mesmo ambiente, sem a necessidade de migração para outras ferramentas. Foi realizado também um estudo sobre o cálculo do TAT equivalente, sendo que dois métodos foram implementados: aproximação por resistores lineares (válida para transistores unitários de canal longo) e aproximação pelo modelo analítico da corrente de dreno através do modelo BSIM3. Três diferentes critérios para a escolha da melhor associação foram abordados e discutidos: menor diferença de corrente entre o TAT e o transistor simples, menor número de transistores unitários e menor condutância de saída. Como circuito de teste, foi realizado o projeto com TATs de um amplificador operacional de dois estágios (amplificador Miller) e a sua comparação com o mesmo projeto utilizando transistores full-custom. Os resultados demonstram que se pode obter bons resultados usando esta técnica, principalmente em termos de desempenho em freqüência. A contribuição da ferramenta LIT ao projeto de circuitos analógicos reside na redução do tempo de projeto, sendo que as tarefas mais suscetíveis a erro são automatizadas, como a geração do leiaute da matriz e das células e o roteamento global. O ambiente de projeto, totalmente gráfico, permite que mesmo projetistas analógicos menos experientes realizem projetos com rapidez e qualidade. Além disso, a ferramenta também pode ser usada para fins educacionais, já que as facilidades proporcionadas ajudam na compreensão da metodologia de projeto.

Paralelizações de métodos numéricos em clusters empregando as bibliotecas MPICH, DECK e Pthread

Este trabalho tem como objetivo desenvolver e empregar técnicas e estruturas de dados agrupadas visando paralelizar os métodos do subespaço de Krylov, fazendo-se uso de diversas ferramentas e abordagens. A partir dos resultados é feita uma análise comparativa de desemvpenho destas ferramentas e abordagens. As paralelizações aqui desenvolvidas foram projetadas para serem executadas em um arquitetura formada por um agregado de máquinas indepentes e multiprocessadas (Cluster), ou seja , são considerados o paralelismo e intra-nodos. Para auxiliar a programação paralela em clusters foram, e estão sendo, desenvolvidas diferentes ferramentas (bibliotecas) que visam a exploração dos dois níveis de paralelismo existentes neste tipo de arquitetura. Neste trabalho emprega-se diferentes bibliotecas de troca de mensagens e de criação de threads para a exploração do paralelismo inter-nodos e intra-nodos. As bibliotecas adotadas são o DECK e o MPICH e a Pthread. Um dos itens a serem analisados nestes trabalho é acomparação do desempenho obtido com essas bibliotecas.O outro item é a análise da influência no desemepnho quando quando tulizadas múltiplas threads no paralelismo em clusters multiprocessados. Os métodos paralelizados nesse trabalho são o Gradiente Conjugação (GC) e o Resíduo Mínmo Generalizado (GMRES), quepodem ser adotados, respectivamente, para solução de sistemas de equações lineares sintéticos positivos e definidos e não simétricas. Tais sistemas surgem da discretização, por exemplo, dos modelos da hidrodinâmica e do transporte de massa que estão sendo desenvolvidos no GMCPAD. A utilização desses métodos é justificada pelo fato de serem métodos iterativos, o que os torna adequados à solução de sistemas de equações esparsas e de grande porte. Na solução desses sistemas através desses métodos iterativos paralelizados faz-se necessário o particionamento do domínio do problema, o qual deve ser feito visando um bom balanceamento de carga e minimização das fronteiras entre os sub-domínios. A estrutura de dados desenvolvida para os métodos paralelizados nesse trabalho permite que eles sejam adotados para solução de sistemas de equações gerados a partir de qualquer tipo de particionamento, pois o formato de armazenamento de dados adotado supre qualquer tipo de dependência de dados. Além disso, nesse trabalho são adotadas duas estratégias de ordenação para as comunicações, estratégias essas que podem ser importantes quando se considera a portabilidade das paralelizações para máquinas interligadas por redes de interconexão com buffer de tamanho insuficiente para evitar a ocorrência de dealock. Os resultados obtidos nessa dissertação contribuem nos trabalhos do GMCPAD, pois as paralelizações são adotadas em aplicações que estão sendo desenvolvidas no grupo.

Alcoxisilanos como doadores de elétrons externo na polimerização de propeno

Este estudo teve como objetivo avaliar o efeito do grupo hidrocarboneto do doador de elétrons externo (DE), do tipo alcoxisilano, no desempenho do catalisador MgCl2.TiCl4.DI-AlEt3/DE na polimerização de propeno. Os DE testados apresentam estrutura genérica do tipo R2Si(OMe)2 e os mesmos foram avaliados com um catalisador comercial do tipo MgCl2.TiCl4.DI na polimerização de propeno em escala de laboratório. Foram empregados quatro doadores de elétrons externos, sendo um usualmente utilizado em processos industriais, o Ciclohexilmetildimetoxisilano (CHMDS). O Diciclopentildimetoxisilano (DCPDS) foi escolhido pela sua alta atividade catalítica e controle da estereoespecificidade, o 3,3,3-trifluorpropilmetildimetoxisilano (FPMDS) pela elevada sensibilidade ao hidrogênio. O quarto DE utilizado foi o Di-n-propildimetoxisilano (DPDS), sendo o mesmo escolhido pela sua elevada atividade catalítica. Os grupos hidrocarbonetos diciclopentil conferiram ao DCPDS o envenenamento seletivo dos sítios catalíticos bem como boa afinidade com o catalisador em estudo. O DCPDS propiciou elevada atividade catalítica, sendo esta em torno de 43 kg PP/g cat.2 h, bem como boa estereoespecificidade, tornando possível a produção de polipropileno com teor de solúveis em xileno da ordem de 1,7 %. O grupo hidrocarboneto 3,3,3-trifluorpropil, presente no FPMDS demonstrou favorecer a ação do hidrogênio, possibilitando a produção de polipropileno com taxa de fluidez superior em duas vezes a dos polímeros obtidos com o uso do CHMDS sob as mesmas condições de polimerização. O DPDS apresentou atividade catalítica similar a do CHMDS, ou seja de 32 kg PP/g cat.2 h, com a desvantagem de apresentar baixo controle da estereoregularidade do catalisador, indicando que os grupos hidrocarbonetos lineares n-propil presentes no DPDS não conferem boa performance a este DE Desta forma, os testes realizados ratificam as informações constantes na literatura para o DCPDS e o FPMDS, ou seja, apresentam determinadas peculiaridades e, quando estas forem requeridas, estes DE devem ser empregados em substituição ao CHMDS. Cabe salientar que o mesmo não ocorreu para o DPDS, porém, neste caso, verificou-se recentemente junto ao fornecedor deste produto que o mesmo foi produzido apenas em escala experimental devido ao seu mau desempenho. No pré-contato entre o catalisador e o cocatalisador, verificou-se que é preferível o emprego de baixas temperaturas (5°C) e pequenos tempos de pré-contato (5 min) entre ambos. Testes realizados indicaram que ao passar de 5 para 10 °C ocorre uma queda da atividade e da estereoespecificidade do catalisador. Para temperaturas de pré-contato da ordem de 10 a 35°C praticamente não houve variação da atividade catalítica, porém, observa-se uma diminuição da estereoregularidade do polipropileno. Intervalos de tempo na faixa de 15 a 120 min ocasionam uma pequena queda na atividade catalítica e os valores de taxa de fluidez deslocam-se para valores superiores. A adição de pequenas quantidades de hidrogênio possibilitou um aumento significativo da atividade catalítica. Neste caso, quando esta foi elevada de 0,004 para 0,014 % em massa (C3H6/H2) houve um acréscimo de 250 % na atividade do sistema catalítico empregado (MgCl2.TiCl4.DIBP- TEA/FPMDS).

Racionalidade imperfeita e inércia inflacionária : uma nova estimação da Curva de Phillips para o Brasil

O presente trabalho estima uma nova relação entre emprego e inflação para o Brasil, tendo como pano de fundo hipóteses propostas pelo arcabouço neo-keynesiano. Quatro hipóteses são testadas e sustentadas ao longo de todo o trabalho: i) os agentes não possuem racionalidade perfeita; ii) a imperfeição no processo de formação das expectativas dos agentes pode ser um fator determinante no componente inercial da inflação brasileira; iii) a inflação possui um componente inercial autônomo, que não depende de choques verificados em mercados isolados; e, iv) relações não-lineares entre inflação e desemprego são capazes de fornecer melhores explicações para o comportamento da economia nos últimos 12 anos. Enquanto as duas primeiras hipóteses são verificadas através de um modelo com mudanças markovianas, as duas últimas são testadas a partir da estimação de uma Curva de Phillips convexa, estimadas pelo Filtro de Kalman. Entretanto, mesmo fazendo uso destas estimativas, as funções de resposta da política monetária apresentam as mesmas propriedades de estimativas tradicionais para o Brasil.

Modelagem estocástica de sistemas hierárquicos de estoques

Esta dissertação apresenta a Modelagem Estocástica de Sistemas Hierárquicos de Estoques (SHEs) como um instrumento de uso gerencial, orientado para a maximização do nível de serviço oferecido aos consumidores finais e a minimização do nível global de estoques em uma cadeia de abastecimento. Tal modelagem envolve otimização de sistemas não-lineares com integrais indefinidas, normalmente tratados através de simulação, dada a impossibilidade de estabelecer uma solução analítica para o problema. Neste trabalho, diferentemente, a solução é obtida através da otimização do sistema, a partir de uma modelagem que estabelece os níveis de estoque-alvo, os pontos de reposição de estoques, os níveis de serviço oferecidos aos consumidores finais e o nível ótimo global de estoques. A sua operacionalização utiliza um aplicativo computacional, projetado e desenvolvido especificamente para este fim. Estruturas hierárquicas, quando otimizadas localmente, em cada um de seus níveis, podem apresentar resultados subótimos em nível global. Assim, não é suficiente identificar os ótimos locais para compor uma solução otimizadora para o sistema. Para resolver esse problema, a maior parte dos pesquisadores utiliza simulação. A qualidade das soluções assim obtidas geralmente é inferior à qualidade da solução obtida mediante otimização global do sistema A presente modelagem tem como objetivo gerar soluções mais qualificadas para o problema da otimização de SHEs. Nesta dissertação, as principais abordagens para otimização de SHEs são apresentadas e analisadas, com destaque para o trabalho de Rosenbaum (1981), cujos resultados, obtidos através de simulação, são comparados com os resultados deste trabalho.

Previsão em tempo atual de cheias com uso de sistema especialista difuso

Sistemas de previsão de cheias podem ser adequadamente utilizados quando o alcance é suficiente, em comparação com o tempo necessário para ações preventivas ou corretivas. Além disso, são fundamentalmente importantes a confiabilidade e a precisão das previsões. Previsões de níveis de inundação são sempre aproximações, e intervalos de confiança não são sempre aplicáveis, especialmente com graus de incerteza altos, o que produz intervalos de confiança muito grandes. Estes intervalos são problemáticos, em presença de níveis fluviais muito altos ou muito baixos. Neste estudo, previsões de níveis de cheia são efetuadas, tanto na forma numérica tradicional quanto na forma de categorias, para as quais utiliza-se um sistema especialista baseado em regras e inferências difusas. Metodologias e procedimentos computacionais para aprendizado, simulação e consulta são idealizados, e então desenvolvidos sob forma de um aplicativo (SELF – Sistema Especialista com uso de Lógica “Fuzzy”), com objetivo de pesquisa e operação. As comparações, com base nos aspectos de utilização para a previsão, de sistemas especialistas difusos e modelos empíricos lineares, revelam forte analogia, apesar das diferenças teóricas fundamentais existentes. As metodologias são aplicadas para previsão na bacia do rio Camaquã (15543 km2), para alcances entre 10 e 48 horas. Dificuldades práticas à aplicação são identificadas, resultando em soluções as quais constituem-se em avanços do conhecimento e da técnica. Previsões, tanto na forma numérica quanto categorizada são executadas com sucesso, com uso dos novos recursos. As avaliações e comparações das previsões são feitas utilizandose um novo grupo de estatísticas, derivadas das freqüências simultâneas de ocorrência de valores observados e preditos na mesma categoria, durante a simulação. Os efeitos da variação da densidade da rede são analisados, verificando-se que sistemas de previsão pluvio-hidrométrica em tempo atual são possíveis, mesmo com pequeno número de postos de aquisição de dados de chuva, para previsões sob forma de categorias difusas.

Marcação a mercado dos fundos de investimento financeiros

O tema central deste trabalho é a marcação a mercado das cotas dos fundos de investimento financeiros. Seu objetivo principal é analisar o procedimento da marcação e como ocorrem as oscilações nas cotas dos fundos. Aborda os conceitos teóricos relacionados aos fundos, bem como estudo dos fatores utilizados para a formação das taxas de juros, do risco, que afeta o retorno dos ativos e da rentabilidade das aplicações financeiras. Analisa a performance, o retorno e as variações, através de instrumentos estatísticos, de um fundo de investimento administrado por instituição financeira no período de fevereiro de 1999 a fevereiro de 2003. A conclusão do trabalho mostra que a marcação a mercado veio trazer maior transparência, credibilidade e segurança à indústria dos fundos, gerando maior estabilidade no mercado financeiro.

Modelagem dos processos envolvidos nos sistemas de secagem e armazenamento de grãos

A produção de soja é uma das principais atividades econômicas na Região Noroeste do Estado do Rio Grande do Sul. As perdas de produto em condições de comercialização ocasionadas nas atividades de secagem e armazenamento são significativas, justificando a pesquisa e aprimoramento destes processos. Nesta tese foram pesquisados dois problemas: 1. Modelamento matemático dos processos de secagem, utilizando parâmetros conhecidos de soja e 2. Modelamento matemático do problema de aeração para o cálculo da distribuição da pressão e da velocidade do ar na massa de grãos em unidades de armazenamento de soja. No problema de secagem foi desenvolvido um sistema composto de quatro equações diferenciais parciais hiperbólicas acopladas não-lineares, que descreve o comportamento da temperatura e do teor de umidade do ar e dos grãos em função do tempo. Para resolver o sistema foram utilizados os métodos das diferenças finitas (p. ex., métodos de MacCormack e Crank- Nicolson.) e o método dos volumes finitos. A análise dos resultados permitiu recomendar o método mais adequado para cada tipo do problema. Para determinação da intensidade do fluxo de massa e de calor foram utilizados os dados experimentais de camada fina obtidos da literatura e complementados com dados experimentais desta tese. Foi desenvolvido um equipamento para obtenção das curvas de secagem de grãos em secador de leito fixo, a fim de identificar o modelo para secagem em camada espessa. A comparação entre os resultados experimentais e das simulações numéricas mostrou que o modelo descreve razoavelmente a dinâmica de secagem No problema de aeração foi desenvolvido um modelo matemático que descreve o escoamento do ar em sistemas de armazenamento de grãos, baseado em relações experimentais entre velocidade e gradiente de pressão. Para resolver o problema de aeração foi utilizado o método dos elementos finitos e desenvolvido um programa computacional. Um teste realizado com o programa mostrou que os resultados da solução numérica convergem para uma solução analítica conhecida. As simulações realizadas mostraram que o programa computacional pode ser usado como instrumento auxiliar para o projeto de silos, possibilitando o cálculo e a visualização gráfica da distribuição das pressões e das linhas de corrente em diferentes seções do armazém.

Ondas e instabilidades em plasmas com inomogeneidades suaves na densidade e no campo magnético de equilíbrio

No presente trabalho, estudamos a absorção e amplificação de ondas eletromagnéticas que se propagam em plasmas com densidade e temperatura fracamente inomogêneas, imersos em um campo magnético também inomogêneo, tendo como base a teoria cinética, dentro do contexto da aproximação local. Esse estudo se dá efetivamente a partir da obtenção do tensor dielétrico do plasma, que deve ser empregado na relação de dispersão. Iniciamos com uma revisão dos conceitos básicos sobre plasmas homogêneos e inomogêneos. Os fundamentos da teoria cinética também foram abordados. Apresentamos uma revisão de trabalhos anteriores que enfocam o mesmo tema, embora descrevendo separadamente os dois tipos de inomogeneidades. A partir desses trabalhos, obtivemos um tensor dielétrico geral, que descreve de forma simultânea as inomogeneidades do campo magnético de equilíbrio e da função distribuição de equilíbrio.Tal tensor foi obtido a partir de um sólido desenvolvimento teórico, que garante a correta descrição da troca de energia entre as ondas e as partículas do plasma. Abordamos os aspectos gerais das instabilidades de deriva, direcionando o estudo à faixa de frequência das ondas híbridas inferiores, e às instabilidades LHDI e MTSI (IWI). Utilizamos perfis lineares de inomogeneidades de campo magnêtico ambiente e densidade para modelar a região da magnetosfera conhecida como neutral sheet. Particularizamos o tensor dielétrico para o estudo específico das instabilidades LHDI e MTSI (IWI), para o tipo de perfil citado acima. Apresentamos uma nova rela»c~ao de dispers~ao para plasmas inomogêneos, que incorpora explicitamente as derivadas espaciais do tensor dielétrico do plasma. Usamos o tensor que unifica os tratamentos das inomogeneidades do campo e densidade nessa relação de dispersão, e obtivemos uma descrição unificada das instabilidades LHDI e MTSI (IWI).

Espalhamento de luz estático e dinâmico em polímeros do tipo polimetacrilato, fluorescentes por transferência protônica intramolecular no estado eletrônico excitado (TPIEE)

As técnicas de espalhamento de luz estático e dinâmico foram utilizadas para a caracterização dos copolímeros de PMMA-benzazolas e do PMMA, em clorofórmio e THF nos regimes diluído e semi-diluído. Os copolímeros foram obtidos pela polimerização do metil-metacrilato na presença dos corantes orgânicos do tipo benzazolas, que caracterizam-se por apresentar uma intensa emissão de fluorescência através de um mecanismo de transferência protônica intramolecular no estado eletrônico excitado. Através da técnica de espalhamento de luz estático foram obtidos parâmetros macromoleculares como a massa molar ponderal média, o raio de giro e o segundo coeficiente virial, bem como o módulo osmótico reduzido e o parâmetro g, relacionado com a arquitetura do polímero em solução. Estes parâmetros indicam que os copolímeros e o PMMA em solução diluída comportam-se como cadeias lineares flexíveis e apresentam-se como esferas homogêneas em solução. Através da espectroscopia de correlação de fótons foram obtidas funções normalizadas de correlação temporal de intensidade correspondentes a um único processo dinâmico tanto para o PMMA como para os copolímeros em regime diluído, independentemente do solvente utilizado. Esta dinâmica corresponde à difusão das cadeias poliméricas em solução. Para as soluções em regime semi-diluído (5£c£30 g·L-1), as funções de correlação temporal de intensidade apresentaram um único decaimento exponencial correspondente a difusão cooperativa dos entrelaçamento das cadeias poliméricas em solução pode ser observado em ambos os solventes. O tempo de relaxação ( t) obtido sugere que o copolímero apresenta, neste regime de diluição, uma dinâmica diferente daquela observada para o PMMA. Portanto, a incorporação da benzazola na cadeia polimérica afeta a dinâmica do polímero tanto em THF como em clorofórmio. Para soluções contendo PMMA e o Copolímero 6 na concentração de 60 g·L-1 em THF observou-se o aparecimento de um segundo movimento, mais lento, nas funções normalizadas de correlação temporal de intensidade, sugerindo a formação de uma estrutura contínua em solução. O mesmo não ocorre utilizando-se clorofórmio como solvente.