Projetos de controladores para sistemas de potência utilizando LMI'S

A presente dissertação tem como objetivo estudar e aprimorar métodos de projetos de controladores para sistemas de potência, sendo que esse trabalho trata da estabilidade dinâmica de sistemas de potência e, portanto, do projeto de controladores amortecedores de oscilações eletromecânicas para esses sistemas. A escolha dos métodos aqui estudados foi orientada pelos requisitos que um estabilizador de sistemas de potência (ESP) deve ter, que são robustez, descentralização e coordenação. Sendo que alguns deles tiveram suas características aprimoradas para atender a esses requisitos. A abordagem dos métodos estudados foi restringida à análise no domínio tempo, pois a abordagem temporal facilita a modelagem das incertezas paramétricas, para atender ao requisito da robustez, e também permite a formulação do controle descentralizado de maneira simples. Além disso, a abordagem temporal permite a formulação do problema de projeto utilizando desigualdades matriciais lineares (LMI’s), as quais possuem como vantagem o fato do conjunto solução ser sempre convexo e a existência de algoritmos eficientes para o cálculo de sua solução. De fato, existem diversos pacotes computacionais desenvolvidos no mercado para o cálculo da solução de um problema de inequações matriciais lineares. Por esse motivo, os métodos de projeto para controladores de saída buscam sempre colocar o problema na forma de LMI’s, tendo em vista que ela garante a obtenção de solução, caso essa solução exista.

Identificação de sistemas contínuos utilizando o método dos momentos de Poisson

A maioria dos métodos de síntese e sintonia de controladores, bem como métodos de otimização e análise de processos necessitam de um modelo do processo em estudo. A identificação de processos é portanto uma área de grande importância para a engenharia em geral pois permite a obtenção de modelos empíricos dos processos com que nos deparamos de uma forma simples e rápida. Mesmo não utilizando leis da natureza, os modelos empíricos são úteis pois descrevem o comportamento específico de determinado processo. Com o rápido desenvolvimento dos computadores digitais e sua larga aplicação nos sistemas de controle em geral, a identificação de modelos discretos foi amplamente desenvolvida e empregada, entretanto, modelos discretos não são de fácil interpretação como os modelos contínuos pois a maioria dos sistema com que lidamos são de representação contínua. A identificação de modelos contínuos é portanto útil na medida que gera modelos de compreensão mais simples. A presente dissertação estuda a identificação de modelos lineares contínuos a partir de dados amostrados discretamente. O método estudado é o chamado método dos momentos de Poisson. Este método se baseia em uma transformação linear que quando aplicada a uma equação diferencial ordinária linear a transforma em uma equação algébrica evitando com isso a necessidade do cálculo das derivadas do sinais de entrada e saída Além da análise detalhada desse método, onde demonstramos o efeito de cada parâmetro do método de Poisson sobre o desempenho desse, foi realizado também um estudo dos problemas decorrentes da discretização de sinais contínuos, como por exemplo o efeito aliasing decorrente da utilização de tempos de amostragem muito grandes e de problemas numéricos da identificação de modelos discretos utilizando dados com tempos de amostragem muito pequenos de forma a destacar as vantagens da identificação contínua sobre a identificação discreta Também foi estudado um método para compensar a presença de offsets nos sinais de entrada e saída, método esse inédito quando se trata do método dos momentos de Poisson. Esse trabalho também comprova a equivalência entre o método dos momentos de Poisson e uma metodologia apresentada por Rolf Johansson em um artigo de 1994. Na parte final desse trabalho são apresentados métodos para a compensação de erros de modelagem devido à presença de ruído e distúrbios não medidos nos dados utilizados na identificação. Esses métodos permitem que o método dos momentos de Poisson concorra com os métodos de identificação discretos normalmente empregados como por exemplo ARMAX e Box-Jenkins.

Obtenção de alfa-olefinas em reações de oligomerização do eteno utilizando complexos de níquel em meio homogêneo e bifásico

Neste trabalho são descritos estudos comparativos de reações de oligomerização do etileno utilizando complexos do tipo: Ni-[P^O] e Ni-[Diimina] (na presença de MAO como cocatalisador). A pesquisa foi realizada com o objetivo de verificar e comparar o comportamento dos dois complexos em estudo, quando estes passam do meio homogêneo para o meio bifásico (na presença de um líquido iônico). O líquido iônico utilizado é o tetrafluorborato de 1-butil-3-metilimidazol. Foram avaliados parâmetros reacionais como pressão, temperatura, tempo de reação, velocidade de agitação e volume de líquido iônico para o catalisador Ni-[P^O] nos meios homogêneo e bifásico. Para o precursor Ni-[Diimina], os parâmetros avaliados foram pressão e temperatura nos dois meios. Os produtos obtidos são olefinas lineares e ramificadas, principalmente α-olefinas; com faixa de distribuição de C4 a C20 (Distribuição Shulz-Flowry). As atividades obtidas nos ciclos de reações variam entre 21 h-1 a 32.772 h-1, e os resultados são interpretados envolvendo os mecanismos de reação. Devido as diferenças existentes entre nos complexos em estudo, verifica-se diferentes atividades e diferentes faixas de produtos oligoméricos. O comportamento obtido para cada complexo depende diretamente dos parâmetros reacionais utilizados, além do meio ao qual foram submetidos durante as reações.

Número de desempenho robusto não-linear : uma nova medida do grau de não-linearidade

nRPN, nRPNSTAT e nRPNDYN são novos índices, propostos nesta dissertação, que quantificam o grau de não-linearidade total, estática e dinâmica, respectivamente, de sistemas dinâmicos não- lineares. Estes novos índices permitem responder se é necessária a utilização de controladores não-lineares ou se simples controladores lineares apresentarão desempenho satisfatório. Além disso, se um controlador não- linear se fizer necessário, quais características este deve apresentar. Nesta dissertação é realizado um estudo acerca das metodologias para determinação do grau de não- linearidade de plant as industriais. Entre as encontradas, a que apresenta as melhores potencialidades é estudada em detalhes, tendo seus atributos e deficiências explorados. As deficiências foram sanadas e um novo conjunto de índices é proposto (nRPN, nRPNSTAT e nRPNDYN), o qual, além de apresentar resultados conclusivos em termos quantitativos, também traz novas informações valiosas para a escolha do controlador adequado, como a parcela estática e dinâmica da não- linearidade. A metodologia proposta é aplicada e testada em um conjunto bastante amplo e complementar de possíveis sistemas dinâmicos não-lineares (tais como: pH, reatores químicos e coluna de destilação de alta pureza) onde o comportamento nãolinear já foi explorado na literatura O desenvolvimento e o estudo do comportamento dinâmico de uma planta laboratorial de cinco tanques acoplados são apresentados. Este sistema possui dinâmica e ganhos variáveis, descontinuidades, inversão no sinal do ganho multivariável e diferentes graus de acoplamento, ilustrando de forma simples vários possíveis comportamentos dinâmicos encontrados em plantas industriais. Devido a esta grande versatilidade a planta foi construída e testada, confirmando os comportamentos previstos nos estudos teóricos.

Planejamento de perturbações para a identificação de modelos dinâmicos

A identificação de modelos é determinante no sucesso das modernas técnicas de controle avançado de processos. Um modelo para o sistema pode ser obtido através de modelagem rigorosa, baseada em equações governantes do sistema ou através da modelagem empírica e estimação de parâmetros. Embora mais rápida e fácil, a modelagem empírica necessita de alguns testes de identificação nos quais as variáveis manipuladas são variadas de modo que resultem em variações nas variáveis controladas. Os testes de identificação podem apresentar custos muito elevados tendo em vista que o sistema pode sair de seu ponto normal de operação, gerando produtos com folga de especificação. Este fato ocorre porque usualmente as perturbações aplicadas nas variáveis manipuladas nas indústrias de processos são independentes umas das outras, aumentando a duração do teste de identificação. Desta forma, neste trabalho foi desenvolvida uma nova metodologia de projeto de perturbações simultâneas para a identificação de modelos dinâmicos baseada na direcionalidade do sistema, com o objetivo de fornecer dados mais ricos para se capturar corretamente o comportamento multivariável do sistema e manter o processo no ponto de operação normal. As perturbações são projetadas conforme as características de um modelo simplificado do processo, ou pré-modelo. Este modelo inicial é obtido essencialmente de dados históricos de planta, selecionados através de uma sistemática análise de correlação desenvolvida neste trabalho A metodologia proposta é composta de duas partes: a primeira parte diz respeito à análise dos dados históricos de planta para obtenção de informações prelimirares as quais são utilizadas no planejamento de perturbações, tais como amplitude do ruído de medida, correlação entre as variáveis de processo, constante de tempo do sistema e matriz de ganhos. E a segunda parte consiste no cálculo da amplitude das perturbações baseado nos resultados da primeira etapa do planejamento. Para sistemas mal-condicionados verificou-se que as perturbações planejadas pela metodologia removem menos a planta de seu ponto de operação gerando resultados mais consistentes em relação às perturbações tradicionais. Já para sistemas bem-condicionados, os resultados são semelhantes. A metodologia foi aplicada em uma unidade piloto experimental e numa unidade de destilação da PETROBRAS, cujos resultados apontam pouca remoção dos sistemas do ponto de operação e modelos consistentes. A validação dos modelos também foi contemplada na dissertação, uma vez que foi proposto um novo critério de validação que considera a derivada dos dados de planta e a do modelo e não apenas os dados de planta e os dados da simulação das saídas do modelo.

Low cost BIST techniques for linear and non-linear analog circuits

With the ever increasing demands for high complexity consumer electronic products, market pressures demand faster product development and lower cost. SoCbased design can provide the required design flexibility and speed by allowing the use of IP cores. However, testing costs in the SoC environment can reach a substantial percent of the total production cost. Analog testing costs may dominate the total test cost, as testing of analog circuits usually require functional verification of the circuit and special testing procedures. For RF analog circuits commonly used in wireless applications, testing is further complicated because of the high frequencies involved. In summary, reducing analog test cost is of major importance in the electronic industry today. BIST techniques for analog circuits, though potentially able to solve the analog test cost problem, have some limitations. Some techniques are circuit dependent, requiring reconfiguration of the circuit being tested, and are generally not usable in RF circuits. In the SoC environment, as processing and memory resources are available, they could be used in the test. However, the overhead for adding additional AD and DA converters may be too costly for most systems, and analog routing of signals may not be feasible and may introduce signal distortion. In this work a simple and low cost digitizer is used instead of an ADC in order to enable analog testing strategies to be implemented in a SoC environment. Thanks to the low analog area overhead of the converter, multiple analog test points can be observed and specific analog test strategies can be enabled. As the digitizer is always connected to the analog test point, it is not necessary to include muxes and switches that would degrade the signal path. For RF analog circuits, this is specially useful, as the circuit impedance is fixed and the influence of the digitizer can be accounted for in the design phase. Thanks to the simplicity of the converter, it is able to reach higher frequencies, and enables the implementation of low cost RF test strategies. The digitizer has been applied successfully in the testing of both low frequency and RF analog circuits. Also, as testing is based on frequency-domain characteristics, nonlinear characteristics like intermodulation products can also be evaluated. Specifically, practical results were obtained for prototyped base band filters and a 100MHz mixer. The application of the converter for noise figure evaluation was also addressed, and experimental results for low frequency amplifiers using conventional opamps were obtained. The proposed method is able to enhance the testability of current mixed-signal designs, being suitable for the SoC environment used in many industrial products nowadays.

Alternativas de manejo para os campos de cima da serra

Os experimentos, foram conduzidos em uma área de campo natural, localizada na Fazenda Boa Vista, no Município de São José dos Ausentes, região dos Campos de Cima da Serra, Rio Grande do Sul. Os objetivos foram descrever o ambiente de pastagem natural da região, determinar a capacidade de suporte destas pastagens e avaliar alternativas de manejo para o inverno. As diferentes ofertas de forragem proporcionaram aumentos lineares no massa de forragem da pastagem, na taxa de acúmulo de matéria seca e no percentual de material morto. Nos parâmetros de resposta animal, as ofertas de forragem geraram modelos lineares para carga animal (kg/ha), GMD (kg/animal/dia), para a avaliação de produção por área foi obtida resposta quadrática, sendo o maior valor 69,7 kg/ha/ano. Para o levantamento florístico constatou-se que o campo é composto basicamente por gramíneas cespitosas estivais, sendo Andropogon lateralis a espécie mais abundante. Após a aplicação de quatro níveis de oferta de forragem (5, 9, 13 e 17% do PV), observou-se um aumento na participação de cyperaceas, principalmente Bulbostylis sphaerocephala. Para a avaliação dos suplementos, observou-se diferença entre os tratamentos, sendo o pastejo temporário superior em GMD à ração e às misturas múltiplas. Verificou-se que a aplicação de calcário em superfície é recomendável e soluciona os problemas decorrentes dos altos teores de Al trocável. A introdução de espécies sobre o campo natural da região é indicada.

Por uma pop'escrita acadêmica educacional

Trata da escrita acadêmica educacional como problema. Desde então, de pensar e praticar tal escrita a partir de algumas aberturas propostas por conceitos desenvolvidos pela Filosofia da Diferença de Gilles Deleuze e Félix Guattari, bem como por procedimentos literários contemporâneos. Aqui, filosofia e literatura em trama instigam a invenção de um espaço para experimentação da escrita. Trata de uma tentativa de lidar com a produção e a criação de uma escrita acadêmica enquanto modo de problematizar e indagar com uma educação que se ocupe do como referir-se aos estados intensivos que agem sobre o linguageiro, incitadores de diferenças e de variações na própria escrita. Trata, também, de um exercício acolhedor de múltiplas linguagens, movendo-se naquilo que aí se agita desde agenciamentos entre formas de conteúdo e formas de expressão potencializadores de irrupturas nos efeitos de um curso regular, previsível, homogêneo e identitário do regime dominante da língua E trata de uma ação estético-política comprometida com relações entre as formas que conjuga ao se abrirem naquilo que vigora o heterogêneo, tornando visíveis múltiplos arranjos de forças que a compõem. Menos um princípio moral e um juízo da razão quando trata de um empreendimento de saúde enquanto crítica e clínica de interposição no traçado de linhas de fuga como fluxos liberadores do desejo. Busca, isto sim, fabricar com a literatura e a filosofia um território de escritura com brechas, hospedeiro de atitudes contingentes que coloquem a fugir o soberano no modo institucional. Uma pop’escrita acadêmica educacional como uma aposta, um jogo ao acaso potente para produzir novos códigos por meio de um movimento de territorialização e desterritorialização com o vigor da experimentação que contém a pura concepção afirmativa da vida em seus desejos.

Solução de uma classe de equações diferenciais parciais em espaços de Niko'skii através do método de Galerkin

Neste trabalho apresenta-se uma solu c~ao para um problema abstrato de Cauchy. Basicamente, d a-se uma formula c~ao abstrata para certos tipos de equa c~oes diferenciais parciais n~ao lineares de evolu c~ao em espa cos de Nikol'skii, tais espa cos possuem boas propriedades de regularidade e resultados de imers~ao compacta, num certo sentido s~ao intermedi arios entre os espa cos de Holder e os espa cos de Sobolev. Aplicando o m etodo de Galerkin, prova-se resultados de exist^encia global de solu c~oes fracas, como tamb em a exist^encia de solu c~oes fracas com a propriedade de reprodu c~ao. E impondo mais hip oteses sobre os operadores envolvidos demonstra-se unicidade de solu c~oes fracas.

Análise experimental e numérica do comportamento térmico de um coletor solar acumulador

Esta tese apresenta um estudo do comportamento térmico de um coletor solar acumulador e desenvolve uma metodologia para medir a sua eficiência diária. O coletor solar acumulador está instalado na face norte do prédio de Energia Solar da UFRGS e possui cerca de 26 m2. É constituído de uma massa espessa de concreto com uma superfície absorvente feita de tijolos, possuindo uma cobertura dupla de vidros colocada de modo a deixar um espaço para a circulação de ar. Os raios solares atravessam a cobertura de vidro e aquecem a massa absorvente de tijolo, a qual aquece o ar que é introduzido no interior da construção por efeito de termossifão. Uma das principais características do coletor solar acumulador consiste no fato de que a resposta do coletor é defasada no tempo. Este fenômeno permite que o coletor entregue calor ao ambiente mesmo após o término da radiação solar. Essa defasagem dos picos de energia térmica ocorre devido ao baixo valor da difusividade térmica do concreto. Este trabalho foi dividido em duas etapas. A primeira etapa consistiu na montagem de um calorímetro para controle e monitoração das variáveis envolvidas. No interior do calorímetro foram instaladas 36 garrafas com água. As temperaturas dos conteúdos das garrafas, do coletor solar e as radiações envolvidas foram monitoradas através de 26 sensores de temperatura de CI, 8 sensores resistivos PT100 e dois sensores de radiação fotovoltaicos. Para obter as medidas dos sensores instalados foi feita a montagem de um sistema de aquisição de dados interfaceado a um microcomputador A segunda etapa consistiu na produção de um programa computacional, escrito em linguagem Fortran 90, para simular o comportamento térmico dos diversos elementos constituintes do coletor, determinar a potência térmica do coletor solar e sua eficiência diária. Para a simulação numérica do coletor solar acumulador, adotou-se um modelo simplificado bidimensional do mesmo. Foi integrada, através do Método dos Volumes Finitos, a equação de difusão de calor transiente em 2 dimensões. Na formulação das equações lineares optou-se pelo emprego das diferenças centrais no espaço e formulação explícita no tempo. Ao todo foram produzidas 4 malhas computacionais, com distintos refinamentos e foi realizado o estudo da estabilidade numérica das diversas malhas. Através da montagem experimental obtiveram-se várias características térmicas do comportamento do sistema, entre as quais, a transmitância da cobertura, curvas de temperatura do ar fornecido ao calorímetro e curva da eficiência diária do coletor solar . Através da simulação numérica foi possível determinar a potência térmica que o coletor entrega para o laboratório, a eficiência do coletor, os campos de temperatura e a vazão mássica nos diversos canais interiores do coletor solar.

Centros urbanos vitais : configuração, dinâmica funcional e caráter das ruas comerciais de Porto Alegre

Este trabalho trata basicamente de investigar a estrutura de polarização da cidade de Porto Alegre. É proposto que a cidade constitui um sistema de diferenciação espacial primário, no qual emergem centros de comércio e serviços em espaços-chave da malha urbana - trechos lineares pertencentes às vias estruturadoras - à medida que esta cresce e se consolida, equilibrando permanentemente a equação entre núcleos de concentração e tecidos recorrentes. Os centros são definidos conceitualmente e seu processo formativo é brevemente analisado. Critérios são elencados para sua seleção e inclusão no corpo de estudo. Níveis intermediários de diferenciação interno à categoria das centralidades são então trabalhados, com base na quantificação das características da forma e das propriedades do arranjo espacial e na leitura qualitativa da dinâmica funcional. Finalmente é construído um nível final de diferenciação entre os centros, combinando os atributos espaciais quantificados, as variáveis funcionais e a própria mecânica de relacionamento entre essas duas esferas. O resultado desta combinação é denominado personalidade urbana, a qual é confrontada com um parâmetro definidor do seu caráter e papel no sistema evolutivo da cidade: o Espaço Urbano Espesso. As referências conceituais e metodológicas situam-se no campo dos estudos configuracionais urbanos, especificamente a Sintaxe Espacial, que trata de investigar a realidade urbana a partir do entendimento de que a forma física e os padrões espaciais interferem ou influenciam os demais fenômenos que manifestam-se na cidade.

Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.

Simulação física de processos gravitacionais subaquasos: uma aproximação para o entendimento da sedimentação marinha profunda

O tema turbiditos tem causado muita controvérsia nos últimos anos. A nosso ver, isto ocorre principalmente devido à diminuição das pesquisas sobre os mecanismos que envolvem a iniciação, o transporte e a deposição deste tipo de rocha. A proposta deste trabalho é avaliar o potencial da simulação física de correntes de turbidez em prever e explicar feições sedimentares em seus depósitos. Foi escolhido como protótipo um sistema turbidítico antigo situado na margem oriental brasileira. A geometria complexa do protótipo foi simplificada para construção do modelo nas instalações do Pavilhão Fluvial do Instituto de Pesquisas Hidráulicas. Foram desenvolvidos doze ensaios onde as observações realizadas sofisticaram-se a partir dos conhecimentos adquiridos nas etapas anteriores. Como resultado dos experimentos identificaram-se novos aspectos geométricos e dinâmicos das correntes de densidade não-conservativas e suas conseqüências na sedimentação. Constatou-se um caráter ondulatório no fluxo, que teve sua origem associada à geração de ondas internas às correntes associadas ao desprendimento de vórtices a partir da cabeça da corrente e sua propagação ao longo da porção superior da corrente. Esta dinâmica implica mudanças na taxa de sedimentação ou mesmo erosão pela corrente, associadas a variações da amplitude e freqüência daquelas ondas. Na cabeça da corrente, verificou-se uma distribuição homogênea de sedimento em suspensão desde a base até o topo da corrente. No corpo ela se divide em duas camadas, uma basal com maior concentração de sedimentos e outra, superior marcada pela expansão do fluxo. Nas quebras de declive do modelo, que ocorrem no meio do canal e no ponto em que a corrente perde o confinamento, foram observadas acelerações localizadas no fluxo. Dentre os parâmetros analisados nos experimentos, constatou-se que a vazão de alimentação tem grande influência nas características dos depósitos. De um modo geral, um aumento da vazão implica um deslocamento do pico deposicional no sentido da corrente e um aumento no conteúdo de frações mais grossas. Constatou-se, nos sedimentos depositados no canal, uma distribuição seqüenciada de formas de leito, que varia entre ripples de crista reta e ripples lingüóides. Reconheceu-se uma correlação entre a amplitude e o comprimento nestas formas de leito. Identificou-se em todos os casos em que houve o extravasamento da corrente a formação de ripples na lateral do canal com cristas lineares que indicam uma direção do fluxo próxima à que ocorre no canal. Foram desenvolvidos depósitos alongados no sentido do fluxo na área onde a corrente perde o confinamento. Observou-se uma grande similaridade entre os depósitos gerados nos experimentos e aqueles identificados em sistemas turbidíticos atuais e do registro geológico, tanto em afloramentos como em dados de subsuperfície.

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.

Integração semi-Lagrangeana num modelo oceano-atmosfera e ajuste geostrófico no domínio tempo

Este trabalho visa o uso da função de Green de valor inicial no ajuste geostrófico e do método Semi-Lagrangeano na integração de um modelo acoplado oceano-atmosfera descrito pelas equações de águas rasas. O ajuste geostrófico é considerado atravées de perturbações na pressão e do vento. No caso de sistemas sem rotação, é discutida a relação da equação hidrostática com ondas longas não-dispersivas. Com rotação, a conservação da vorticidade potencial permite escolher a elevação correspondente a um estado de equilíbrio geostrófico. O sistema de equações de águas rasas é desacoplado em equações de Klein-Gordon com valores iniciais e termos não-homogêneos acoplados. A resposta dinâmica formada pela resposta transiente e a resposta forçada é obtida para uma perturbação inicial da elevação. A ação do vento como forçante nas equações de momento 2D, através do transporte de Eckman, conduz a uma equação de águas rasas forçada. Uma decomposição da resposta forçada é realizada com uma resposta permanente, que satisfaz a equação de Helmholtz , e com o uso da base dinâmica gerada pela resposta impulso. Um modelo hidrodinâmico 3D introduzido por Casulli e governado por equações não-lineares de águas rasas é integrado na vertical para a obtenção de um modelo 2D. Com isto, as condições de contorno devido a tensão do vento e a fricção devido a topografia do fundo, transformam-se em forçantes do modelo. O modelo foi integrado com um método semi-implícito em diferenças finitas, utilizando-se o método Semi-Lagrangeano para a parte advectiva. Simulações simbólicas foram realizadas para o ajuste geostrófico devido a perturbações de duração infinita e finita para a elevação e para o efeito da tensão do vento. Foram realizadas simulações numéricas para variadas geometrias, em particular a Baia de Guanabara e a Lagoa do Patos.