Espécies de moscas frugívoras (Diptera: Tephritidae e Lonchaeidae), quantificação de danos e avaliação de medidas para o seu manejo em pomares orgânicos de citros

Aspectos relativos às moscas frugívoras (Diptera: Tephritidae e Lonchaeidae), em sistemas de produção familiares e com princípios ecológicos, são pouco conhecidos, assim como alternativas para o seu manejo. Este estudo foi realizado com o objetivo de reconhecer as espécies de moscas frugívoras e sua variação numérica, quantificar os danos e avaliar as alternativas de manejo utilizadas pelos citricultores no sistema orgânico de produção de citros, na região do vale do rio Caí, Rio Grande do Sul. O experimento foi instalado em pomares de laranjeira ‘Céu’ (Citrus sinensis), de janeiro a maio de 2003 e 2004, e do tangoreiro ‘Murcott’ (Citrus reticulata x Citrus sinensis), de junho a setembro de 2003 e junho a agosto de 2004, na fase de maturação. Para o monitoramento populacional das moscas foram utilizadas armadilhas tipo McPhail, contendo suco de uva a 25%, verificadas semanalmente. Os danos foram avaliados por meio da porcentagem de frutos danificados. As temperaturas máxima, mínima e média, a precipitação, a umidade relativa do ar e a coloração dos frutos, foram registradas. Os métodos de controle (calda sulfocálcica, soro de leite e ensacamento) foram analisados pela redução dos danos e os impactos na composição dos táxons dos pomares. A viabilidade econômica destes métodos também foi calculada. Constatou-se que Anastrepha fraterculus é a espécie de maior importância para a região. A precipitação é o principal fator climático que afeta o número de moscas frugívoras Em laranjeira ‘Céu’ e em tangoreiro ‘Murcott’, no período que antecede a mudança de coloração dos frutos, os danos ocorrem de forma atrasada em relação ao pico populacional, após este período a resposta é imediata. O ensacamento é o método que promove a maior redução dos danos e com menores impactos na composição de táxons, no entanto, a sua viabilidade econômica é dependente do número de moscas-das-frutas e do valor a ser recebido pela produção.

Influência de parâmetros físicos no dimensionamento de chaminés de equilíbrio simples de usinas hidrelétricas

Tendo em vista que as chaminés de equilíbrio são estruturas importantes na proteção dos circuitos hidráulicos das usinas hidrelétricas contra os fenômenos transientes que ocorrem rotineiramente durante sua operação, realizou-se uma investigação sobre a influência dos parâmetros físicos no desempenho e conseqüentemente no dimensionamento de chaminés de equilíbrio simples localizadas a jusante das turbinas, caso aplicado a Usina Hidrelétrica de Serra da Mesa. Os parâmetros estudados foram: área da chaminé, comprimento e área do túnel de fuga, perda de carga, vazão, nível no canal de fuga e tempo de manobra. Foram coletados dados na Usina Hidrelétrica de Serra da Mesa e em uma estrutura laboratorial existente no LAHE, ambas pertencentes a FURNAS Centrais Elétricas S.A, referentes a situações de regimes permanentes (operação normal) e de regimes transientes. Estes dados foram utilizados na determinação da perda de carga linear nos túneis e das características das oscilações em suas chaminés de equilíbrio. Para o estudo da perda de carga linear em túneis, também foram feitas medições, em regime permanente, na Usina Hidrelétrica de Salto Forqueta – CERTEL - Cooperativa de Eletrificação Teutônia, e utilizados os dados da Usina Hidrelétrica Bugres fornecidos pela CEEE – Companhia Estadual de Energia Elétrica –RS. Com estes dados foram estabelecidas relações entre o coeficiente de atrito e o número de Reynolds para cada estrutura e realizada a sua comparação com valores teóricos Verificou-se, através da comparação com dados experimentais que os quatro modelos numéricos disponíveis estavam aptos à aplicação na seqüência da pesquisa, porém, elegeu-se o modelo CHAMINE-IPH devido a sua maior versatilidade e simplicidade de uso. Com a análise dos resultados das simulações numéricas realizadas foi possível a identificação da influência dos parâmetros físicos no comportamento do nível d’água em chaminés de equilíbrio simples durante regimes transientes. Constatou-se que para a reprodução numérica de fenômenos transitórios em chaminés de equilíbrio a jusante das turbinas, é fundamental a consideração da oscilação do nível no canal de fuga e dos coeficientes de perda de carga adequados.

Algoritmos adaptativos para o método GMRES(m)

Nesse trabalho apresentamos algoritmos adaptativos do M´etodo do Res´ıduo M´ınimo Generalizado (GMRES) [Saad e Schultz, 1986], um m´etodo iterativo para resolver sistemas de equa¸c˜oes lineares com matrizes n˜ao sim´etricas e esparsas, o qual baseia-se nos m´etodos de proje¸c˜ao ortogonal sobre um subespa¸co de Krylov. O GMRES apresenta uma vers˜ao reinicializada, denotada por GMRES(m), tamb´em proposta por [Saad e Schultz, 1986], com o intuito de permitir a utiliza¸c˜ao do m´etodo para resolver grandes sistemas de n equa¸c˜oes, sendo n a dimens˜ao da matriz dos coeficientes do sistema, j´a que a vers˜ao n˜ao-reinicializada (“Full-GMRES”) apresenta um gasto de mem´oria proporcional a n2 e de n´umero de opera¸c˜oes de ponto-flutuante proporcional a n3, no pior caso. No entanto, escolher um valor apropriado para m ´e dif´ıcil, sendo m a dimens˜ao da base do subespa¸co de Krylov, visto que dependendo do valor do m podemos obter a estagna¸c˜ao ou uma r´apida convergˆencia. Dessa forma, nesse trabalho, acrescentamos ao GMRES(m) e algumas de suas variantes um crit´erio que tem por objetivo escolher, adequadamente, a dimens˜ao, m da base do subespa¸co de Krylov para o problema o qual deseja-se resolver, visando assim uma mais r´apida, e poss´ıvel, convergˆencia. Aproximadamente duas centenas de experimentos foram realizados utilizando as matrizes da Cole¸c˜ao Harwell-Boeing [MCSD/ITL/NIST, 2003], que foram utilizados para mostrar o comportamento dos algoritmos adaptativos. Foram obtidos resultados muito bons; isso poder´a ser constatado atrav´es da an´alise das tabelas e tamb´em da observa ¸c˜ao dos gr´aficos expostos ao longo desse trabalho.

Um modelo analítico e computacional para consideração de efeitos de envelhecimento em estruturas de material compósito

Os materiais compósitos reforçados por fibras apresentam vantagens quando comparados aos materiais de construção mais tradicionais como concreto e aço. Por outro lado, devido ao fato destes materiais serem relativamente recentes no mercado, questões a respeito de sua durabilidade são ainda objeto de discussão e faz-se necessária intensa pesquisa sobre o envelhecimento dos compósitos. Como conseqüência, recentemente têm surgido inúmeros trabalhos à respeito da degradação dos compósitos considerando efeitos como temperatura, oxidação, radiação UV, condições de carregamento, etc. A maioria destas pesquisas, no entanto, são realizadas a nível de material e não são diretamente aplicáveis à situações de projeto. Desta forma, existe grande demanda por novos estudos e dados compatíveis com aplicações estruturais. Neste trabalho apresenta-se um modelo analítico-numérico adequado para, interpretação e aplicação destes dados experimentais em análise e projeto de estrutural. A formulação proposta inclui relações constitutivas elásticas anisotrópicas com envelhecimento, relações constitutivas viscoelásticas anisotrópicas com envelhecimento em termos de variáveis de estado, análise de falhas com critério de degradação ajustado à idade do material e considera-se grandes deslocamentos e pequenas deformações. As diferenças essenciais entre os processos de envelhecimento em endurecimento e amolecimento são descritos juntamente com as relações constitutivas para cada caso. Estas equações são deduzidas na forma adequada para análise numérica via método dos elementos finitos usando uma solução incremental-iterativa com consideração de efeitos pos-críticos. Vários exemplos são apresentados, incluindo análises elásticas, viscoelásticas e de falha com envelhecimento.

Dinâmica monetária eficiente sob encontros aleatórios: uma classe de métodos numéricos que exploram concavidade

A dificuldade em se caracterizar alocações ou equilíbrios não estacionários é uma das principais explicações para a utilização de conceitos e hipóteses que trivializam a dinâmica da economia. Tal dificuldade é especialmente crítica em Teoria Monetária, em que a dimensionalidade do problema é alta mesmo para modelos muito simples. Neste contexto, o presente trabalho relata a estratégia computacional de implementação do método recursivo proposto por Monteiro e Cavalcanti (2006), o qual permite calcular a sequência ótima (possivelmente não estacionária) de distribuições de moeda em uma extensão do modelo proposto por Kiyotaki e Wright (1989). Três aspectos deste cálculo são enfatizados: (i) a implementação computacional do problema do planejador envolve a escolha de variáveis contínuas e discretas que maximizem uma função não linear e satisfaçam restrições não lineares; (ii) a função objetivo deste problema não é côncava e as restrições não são convexas; e (iii) o conjunto de escolhas admissíveis não é conhecido a priori. O objetivo é documentar as dificuldades envolvidas, as soluções propostas e os métodos e recursos disponíveis para a implementação numérica da caracterização da dinâmica monetária eficiente sob a hipótese de encontros aleatórios.

Modelagem e simulação do campo de velocidades do vento em tormentas elétricas

Esta tese propõe um modelo para descrever o campo de velocidade do vento originado numa tormenta elétrica. O escoamento descendente que ocorre durante uma tormenta elétrica madura é modelado como um campo de fluxo axial simétrico. A evolução da velocidade do vento com o tempo é considerada por meio de uma função de modulação empírica. As flutuações da velocidade ao redor da velocidade média em evolução são representadas por um processo aleatório do tipo ruído branco em faixa. Todos os parâmetros utilizados no modelo estão relacionados a variáveis meteorológicas, tais como a altura e diâmetro da nuvem cumulonimbus causadora do evento, a diferença de pressão e a distância entre o solo e a base da nuvem. Adicionalmente considera-se que a nuvem é transportada numa direção pré escolhida pelo vento de fundo, que pode ser devido à circulação geral na região ou a um ciclone extra-tropical. O campo de velocidades do vento resulta da superposição de ambos campos. O modelo é validado comparando as velocidades do vento determinadas em eventos simulados, com registros de estações meteorológicas. Adicionalmente, perfis verticais previstos pelo modelo são comparados com vários perfis para tormentas elétricas propostos na literatura. Finalmente, séries de velocidade máximas anuais a altura padrão de 10m foram simuladas para três estações meteorológicas no sul do Brasil, para as quais dados experimentais são conhecidos, em termos do número médio de tormentas elétricas por ano e das funções distribuição de probabilidade que caracterizam o vento de fundo. A comparação com séries de registros máximos existentes é satisfatória, dando credibilidade adicional ao modelo proposto.

Análise de métodos numéricos para precificação de opções

Analisa e compara, sob as óticas de acurácia e custo, os principais métodos numéricos (simulação de Monte carlo, métodos de diferenças finitas, e modelos baseados em lattice para precificação de opções financeiras, exóticas e reais. Aponta as situações onde cada método deve ser empregado, e apresenta os algoritmos necessários para implementação dos métodos numéricos na prática.

Aplicação de sensoriamento remoto orbital no mapeamento de unidades vulcano-sedimentar no Platô da Ramada, Vila Nova do Sul, RS

Os produtos de sensoriamento remoto gerados pelos novos sensores orbitais, aliados ao desenvolvimento de sistemas computacionais e de técnicas de processamento de imagens possibilitam a manipulação numérica das imagens, visando o realce de informações de interesse para a interpretação visual ou automática. O uso dessas técnicas traz boas perspectivas no mapeamento geológico básico, pois possibilita uma visão sinóptica de amplas áreas e a integração de dados litológicos e estruturais de várias fontes. Este trabalho selecionou as imagens do sensor ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) devido as características de suas resoluções espectral e espacial, como as mais promissoras para identificação dos diferentes litotipos da área do Platô da Ramada, em especial as unidades vulcânicas e hipabissais relacionadas e sedimentares subjacentes. O Platô da Ramada, região de Vila Nova do Sul (RS) constitui uma feição geomorfológica de destaque na área, facilmente identificável nas imagens orbitais, sendo formado por uma seqüência vulcânica caracterizada por depósitos efusivos e piroclásticos, de composição dominantemente ácida. Representa uma fração significativa do magmatismo alcalino sódico neoproterozóico no Escudo Sul-rio-grandense e marca um dos ciclos vulcânicos do período pós-colisional do Ciclo Brasiliano Este estudo testou e avaliou diversos processamentos das imagens multiespectrais para diferenciação litológica, a identificação de alvos e a definição de morfoestruturas da área do Platô da Ramada. A integração de dados geológicos existentes com os produtos obtidos neste estudo possibilitou novas informações anexadas ao mapa geológico utilizado como verdade terrestre. O processamento utilizando a técnica de Transformação por Componentes Principais Seletivas proporcionou os melhores resultados realçando diferenças espectrais existentes entre as rochas vulcânicas e hipabissais e as rochas sedimentares. A partir dessa técnica, selecionou-se as imagens CP2 dos pares das bandas 4-5 (R), 3-4 (G) e 8-9 (B) na geração de uma composição colorida. Esta imagem permitiu a diferenciação espectral entre as rochas vulcânicas do Platô da Ramada e as rochas sedimentares do Grupo Maricá, bem como a individualização, no Grupo Maricá, de duas subunidades, levando-se em conta a densidade de intrusões de diques riolíticos ao norte da área. Identificou-se, ainda, um corpo diorítico com forma elíptica na borda SW do Platô da Ramada, não registrado anteriormente.

Solução semi-analítica da equação de Langevin assintótica para o deslocamento aleatório pelo método Picard

Neste trabalho é desenvolvida uma solução semi-analítica para a Equação de Langevin assintótica (Equação de Deslocamento Aleatório) aplicada à dispersão de poluentes na Camada Limite Convectiva (CLC). A solução tem como ponto de partida uma equação diferencial de primeira ordem para o deslocamento aleatório, sobre a qual é aplicado o Método Iterativo de Picard. O novo modelo é parametrizado por um coeficiente de difusão obtido a partir da Teoria de Difusão Estatística de Taylor e de um modelo para o espectro de turbulência, assumindo a supersposição linear dos efeitos de turbulência térmica e mecânica. A avaliação do modelo é realizada através da comparação com dados de concentração medidos durante o experimento de dispersão de Copenhagen e com resultados obtidos por outros quatro modelos: modelo de partículas estocástico para velocidade aleatória (Modelo de Langevin), solução analítica da equação difusão-advecção, solução numérica da equação difusão-advecção e modelo Gaussiano. Uma análise estatística revela que o modelo proposto simula satisfatoriamente os valores de concentração observados e apresenta boa concordância com os resultados dos outros modelos de dispersão. Além disso, a solução através do Método Iterativo de Picard pode apresentar algumas vantagem em relação ao método clássico de solução.

Uso do método das interfaces coesivas na análise da simulação do comportamento de fratura em materiais frágeis

Este trabalho apresenta uma abordagem numérica da fratura baseada no Método dos Elementos Finitos na qual são considerados materiais frágeis como objeto de estudo. As simulações realizadas são possíveis através de um algoritmo computacional que consiste na implementação de molas entre as faces dos elementos finitos. Esta composição constitui o método denominado Interfaces Coesivas. As molas fornecem o efeito de coesão e é representado por leis constitutivas do tipo bi-linear. Sua descrição é apresentada, bem como sua implementação na análise. Tanto o modo I como o modo II de ruptura são considerados. Com o objetivo de verificar o presente modelo com outros resultados e soluções, são analisados exemplos numéricos e experimentais extraídos da literatura. Tais exemplos tratam da nucleação, propagação e eventuais paradas de trincas em peças estruturais submetidas aos mais diversos carregamentos e condições de contorno. É ainda discutido como os parâmetros de fratura do material interferem no processo de ruptura. Uma tentativa é feita de determinar estes parâmetros, em especial no caso de concreto.

Estudo de recapeamentos asfálticos de pavimentos através de ensaios acelerados

O objetivo deste estudo é avaliar o desempenho de recapeamentos asfálticos, aplicados sobre pavimentos severamente trincados, através do acompanhamento da degradação ocorrida em função do tráfego. Foram testados dois tipos de recapeamentos: uma camada de 5,6 cm de espessura, em concreto asfáltico convencional e um tratamento superficial duplo, modificado por polímero SBS. Os recapeamentos asfálticos foram aplicados sobre uma pista experimental já trafegada e realizou-se, artificialmente, um trincamento padronizado buscando reproduzir efeitos similares àqueles provocados pelo fenômeno de fadiga. Utilizou-se a técnica de ensaios acelerados de pavimentos através da utilização de um simulador linear de tráfego, com rodado duplo. A carga de semi-eixo de 50 kN adotada proporcionou a degradação do pavimento três meses de ensaio para cada recapeamento. O monitoramento do pavimento, em termos estruturais e funcionais, fez-se pelo levantamento periódico de deflexões, bacias deflectométricas, afundamentos de trilha de roda e o controle do surgimento de trincas na superfície. Para melhor quantificar e compreender os fenômenos de degradação estrutural, instrumentou-se a pista experimental com sensores (células de tensão total e strain gages) na interface entre pavimento trincado e recapeamentos, obtendo-se tensões e deformações nestes locais, considerados críticos no desempenho de pavimentos recapeados. Durante os períodos de ensaio, foram monitoradas as condições pluviométricas e temperaturas do ar e pavimento. A determinação dos módulos elásticos das camadas das estruturas testadas foi efetuada pela aproximação das bacias obtidas em campo com os dados gerados em uma simulação numérica, através do programa computacional FLAPS. Percebeu-se que, além do tráfego, o trincamento no recapeamento em tratamento superficial foi devido à condição estrutural inicial do pavimento existente, severidade do trincamento artificial incompatível com a espessura desse tipo de restauração e baixas temperaturas. Quanto ao recapeamento em concreto asfáltico, percebeu-se que o maior trincamento surgiu em áreas sem trincas artificiais subjacentes, partindo do fundo da camada asfáltica do pavimento existente, em função da pouca rigidez da camada granular.

O custo de capital em operações de securitização de recebíveis de empresas não-financeiras por meio da emissão de quotas de Fundos de Investimento em Direitos Creditórios - FIDCs

Esta dissertação investiga, em caráter exploratório, se há evidência de que o custo de capital proporcionado a empresas brasileiras não-financeiras em operações de securitização por meio da estruturação de Fundos de Investimento em Direitos Creditórios – FIDCs (espécie de fundo de investimento regulamentada em 2001 pela Comissão de Valores Mobiliários – CVM) pode ser inferior ao custo de capital associado a operações de crédito tradicionais. A investigação baseia-se em comparação das notas de rating atribuídas a FIDCs estruturados para adquirir direitos creditórios de empresas não-financeiras e ao endividamento de longo prazo dessas mesmas empresas. A comparação entre notas de rating é estruturada a partir da premissa de que há forte correlação entre rating e spread (principal componente do custo de capital de terceiros), conforme documentado em diversos estudos empíricos precedentes. Foram selecionadas 23 operações de emissão de cotas de FIDCs, desde sua regulamentação até 2007, que atendem aos seguintes critérios: terem sido objeto de distribuição pública registrada na Comissão de Valores Mobiliários; terem mais de 5 cotistas seniores; e terem sido constituídos para adquirir recebíveis de uma única empresa não-financeira. A análise das 23 ocorrências fornece evidência de que a estruturação de FIDCs pode resultar na emissão de títulos securitizados com risco de crédito inferior àquele associado ao endividamento de longo prazo da empresa cedente. A redução do risco de crédito é observada por meio do diferencial de notas de rating em 8 das 23 operações analisadas. A ocorrência de 5 operações em que não há diferencial sugere, porém, que a redução do custo de capital não é a única razão pela qual os FIDCs vêm sendo estruturados. As limitações do rating como índice do custo de capital e as limitações numérica e temporal da amostra inviabilizam afirmações conclusivas acerca do impacto potencial do FIDC sobre o custo de capital. Não obstante, a análise sugere novas linhas de pesquisa que podem lançar mais luz sobre o assunto.

Dissipação de energia em estruturas com utilização de polímeros termoplásticos

Esta dissertação tem por objetivo contribuir à compreensão dos mecanismos e efeitos de dissipação de energia em estruturas mediante a utilização de materiais viscoelásticos. São obtidas experimentalmente as propriedades destes polímeros termoplásticos enfatizando seu uso como meio de redução das amplitudes de vibrações em estruturas submetidas a ações dinâmicas. Inicialmente, apresenta-se uma breve resenha histórica sobre o uso de materiais viscoelásticos com finalidade de amortecimento estrutural. São descritos modelos reológicos referentes ao comportamento mecânico, químico e térmico e também modelos matemáticos relacionados ao fenômeno do amortecimento necessários para o desenvolvimento do projeto dos amortecedores. A seguir, descreve-se detalhadamente o projeto, construção e análise de dissipadores de energia confeccionados com polímeros termoplásticos, mediante aplicações de carregamentos harmônicos cíclicos, é análisado seu comportamento. Finalmente, faz-se uma análise numérica da influência da incorporação de dissipadores numa estrutura real com problemas de vibrações excessivas. Conclusões sobre a aplicação deste tipo de amortecedor em estruturas em geral são apresentadas.

Aplicação de um modelo de dispersão atmosférico na região de Candiota-RS e sua relação com as condições meteorológicas

Este trabalho tem por objetivo, estudar a dispersão de poluentes na região da Usina Termelétrica Presidente Médici (UTPM), localizada no município de Candiota-RS, utilizando a modelagem numérica (mais especificamente o modelo ISCST - Industrial Surce Complex Term Short Term) e relacionar o seu comportamento com as condições meteorológicas. Esse modelo de dispersão utilizado é recomendado pela EPA (Environmental Protection Agency) para tratamento de dispersão de poluentes emitidos por fontes industriais como refinarias, termelétricas, siderúrgicas, etc. O período de estudo foi de 19/01/03 à 28/01/03, quando se realizou a primeira campanha experimental, pertencente ao projeto "Estudo Ambiental Aplicando a Modelagem Numérica na Região de Candiota-RS". Os resultados simulados pelo modelo foram comparados com os dados das concentrações de Material Particulado (MP) e Dióxido de Enxofre (SO 2 ), medidos através dos pontos receptores em duas estações de qualidade do ar próximas à Usina Termoelétrica. Uma correlação entre as concentrações de poluentes modelados e os dados coletados foi realizada, assim como, uma comparação com os padrões de qualidade do ar estipulados CONAMA (Conselho Nacional do Meio Ambiente). Os resultados denotam que o modelo ISCST subestima as concentrações de MP e SO 2 observadas nas estações de qualidade do ar, o que provavelmente ocorre, devido à presença de outras fontes de emissão (como indústrias de cimento e calcário) nas proximidades das estações. A situação com a melhor correlação nos dez dias, entre modelo e amostrador, apresentou-se com 90%. Em apenas um dos dez dias, foi excedido o padrão secundário (100 m g/m 3 ) com o poluente SO 2 , caracterizado a qualidade do ar da região como boa para esse período.

Métodos multigrid paralelos em malhas não estruturadas aplicados à simulação de problemas de dinâmica de fluidos computacional e transferência de calor

Fenômenos naturais, tecnológicos e industriais podem, em geral, ser modelados de modo acurado através de equações diferenciais parciais, definidas sobre domínios contínuos que necessitam ser discretizados para serem resolvidos. Dependendo do esquema de discretização utilizado, pode-se gerar sistemas de equações lineares. Esses sistemas são, de modo geral, esparsos e de grande porte, onde as incógnitas podem ser da ordem de milhares, ou até mesmo de milhões. Levando em consideração essas características, o emprego de métodos iterativos é o mais apropriado para a resolução dos sistemas gerados, devido principalmente a sua potencialidade quanto à otimização de armazenamento e eficiência computacional. Uma forma de incrementar o desempenho dos métodos iterativos é empregar uma técnica multigrid. Multigrid são uma classe de métodos que resolvem eficientemente um grande conjunto de equações algébricas através da aceleração da convergência de métodos iterativos. Considerando que a resolução de sistemas de equações de problemas realísticos pode requerer grande capacidade de processamento e de armazenamento, torna-se imprescindível o uso de ambientes computacionais de alto desempenho. Uma das abordagens encontradas na literatura técnica para a resolução de sistemas de equações em paralelo é aquela que emprega métodos de decomposição de domínio (MDDs). Os MDDs são baseados no particionamento do domínio computacional em subdomínios, de modo que a solução global do problema é obtida pela combinação apropriada das soluções obtidas em cada um dos subdomínios Assim, neste trabalho são disponibilizados diferentes métodos de resolução paralela baseado em decomposição de domínio, utilizando técnicas multigrid para a aceleração da solução de sistemas de equações lineares. Para cada método, são apresentados dois estudos de caso visando a validação das implementações. Os estudos de caso abordados são o problema da difusão de calor e o modelo de hidrodinâmica do modelo UnHIDRA. Os métodos implementados mostraram-se altamente paralelizáveis, apresentando bons ganhos de desempenho. Os métodos multigrid mostraram-se eficiente na aceleração dos métodos iterativos, já que métodos que utilizaram esta técnica apresentaram desempenho superior aos métodos que não utilizaram nenhum método de aceleração.