Avaliação de um biorreator piloto de lodo ativado no controle de odores.

Apesar da importância das estações de tratamento de efluentes industriais (ETEIs) na conservação dos ecossistemas, estas podem ser geradoras de gases com maus odores contendo compostos orgânicos voláteis - COVs. Os maus odores têm sido motivos de protestos e reclamações por parte da população circunvizinha às fontes emissoras. Em virtude da conscientização ambiental, e dos impactos sobre a saúde do homem, o objetivo geral deste trabalho é avaliar a eficiência de um biorreator aeróbio piloto no controle de gases odoríferos emitidos em estação de tratamento de efluentes de indústrias de alimentos. Foi desenvolvido um sistema de difusão de ar odorífero em um reator aeróbio de lodo ativado, hermeticamente fechado, operado no regime de batelada sequencial, durante os dias 14, 21, 23, e 30 do mês de julho.Foram realizadas análises dos parâmetros físico-químicos do lodo ativado utilizado no reator aeróbio piloto, como determinação dos sólidos, DBO5, DQO, OD, pH, temperatura e IVL. A atividade da biomassa do lodo ativado foi avaliada por meio do teste de Respirometria. A eficiência do reator quanto à redução da DQO dos gases foi analisada por meio da absorção dos gases em solução de dicromato de potássio. Para avaliação da mensuração dos odores utilizou-se o método de cromatografia gasosa e espectrometria de massa, quantificando amostras de gases odorantes adsorvidas em tubos de carvão ativado, na entrada e na saída dobiorreator. Os resultados obtidos confirmaram o potencial do sistema de difusão em lodos ativados para o tratamento de gases odoríferos em ETEI, com eficiência de remoção dos COV`s variando de 97,3% a 98,9%.

Estudo da efetividade e segurança de um novo esquema posológico da solução saturada de iodeto de potássio no tratamento da esporotricose linfocutânea e cutânea fixa

Fundamentos: A primeira escolha para o tratamento da esporotricose cutânea é o itraconazol oral, contudo o aumento de casos ocorridos numa epidemia de transmissão zoonótica fez com que se buscassem alternativas efetivas e seguras de tratamento. Objetivo: Avaliar uma nova posologia do iodeto de potássio (KI) como alternativa para o tratamento das formas cutâneas limitadas da esporotricose. Métodos: Foram incluídos 102 pacientes com esporotricose, diagnosticados através do isolamento do Sporothrix sp. e divididos em dois grupos que receberam doses diferentes de KI: grupo A recebeu a dose anterior (média 4,83g/dia) e grupo B, a dose reduzida (média 2,52g/dia). Os critérios de cura se basearam em dados clínicos e sorológicos. Resultados: Setenta e nove pacientes (77,4%) obtiveram cura clínica, 70,6 e 84,3% respectivamente nos grupos A e B. Dezesseis pacientes (15,6%) perderam o seguimento e sete trocaram de medicamento: cinco no grupo A e dois no grupo B. A incidência de eventos adversos foi a mesma nos dois grupos (64,7%) com predomínio do gosto metálico (44%), seguido por intolerância gastrointestinal leve e erupção acneiforme (10,7% cada). Nenhum evento adverso grave ocorreu e não houve recidivas. A análise dos desfechos não demonstrou diferença estatisticamente significativa entre os grupos (p=0,9255). A melhora dos títulos sorológicos foi significativa em ambos grupos. Conclusões: O iodeto de potássio, sob a forma de solução saturada, administrado em dose e frequência reduzidas, pode ser usado como alternativa efetiva e segura para o tratamento da esporotricose cutânea. Através de análises estatísticas, a posologia habitual não demonstrou ser superior à proposta por este estudo. A sorologia para esporotricose pode ser utilizada como ferramenta valiosa para acompanhamento clínico destes pacientes.

Materiais compósitos biocidas baseados em resinas reticuladas contendo nanopartículas de prata

Nesta dissertação, duas resinas reticuladas comerciais (denominadas, AmberliteGT73 e AmberliteIRC748) foram empregadas para suportar partículas de prata. Os grupos tiol da GT73 e ácido iminodiacético da IRC748 presentes nessas resinas foram empregados para a fixação de íons Ag+ a partir de solução aquosa. Posteriormente, os íons Ag+ foram reduzidos pelo emprego de três redutores diferentes em pH alcalino, denominados hidrazina, hidroxilamina e formaldeído (pH  12). A morfologia e a impregnação de prata dos materiais binários assim obtidos foram avaliadas por meio de microscópio eletrônico de varredura equipado com detector de elétrons retro-espalhados (SEM-BSE). O detector de espectrometria de energia dispersiva de raios-X (EDAX) acoplado ao SEM permitiu a observação de partículas de prata. Os espectros de raios-X revelaram a presença do metal nas superfícies interna e externa das microesferas dos compósitos. A quantidade de prata incorporada foi determinada pelo método titulométrico, empregando solução padrão de tiocianato de potássio. As características antibactericidas dos compósitos foram avaliadas em colunas contendo pérolas de resina por onde foram percoladas suspensões da bactéria Escherichia coli auxotrópica AB1157 (tipo selvagem) nas concentrações de 103 a 107 células/mL. A avaliação biocida mostrou que estes materiais foram completamente bactericidas, sendo efetivos na eliminação da bactéria em poucos minutos. Esta ação biocida foi atribuída à combinação da atuação da prata e dos grupos funcionais das resinas

Avaliação da remoção de níquel em solução aquosa pelo processo de osmose inversa

O objetivo deste trabalho foi investigar o desempenho de uma membrana comercial na remoção de um metal pesado (níquel) de efluentes sintéticos, por osmose inversa. Na primeira etapa foi realizada uma comparação com os resultados obtidos com soluções de alimentação contendo sais, como NaCl, NaNO3 e Ni(NO3)2.6 H2O, nas concentrações de 50, 100 e 200 ppm, e nas pressões de 10, 20 e 26 bar. Os resultados mostraram que a influência da concentração e da pressão aplicada ao sistema não afetaram as rejeições de forma significativa. Na segunda etapa, como os parâmetros não influíram significativamente na rejeição dos sais, optou-se, pela aplicação de uma pressão de operação de 10 bar, para avaliar a eficiência de remoção de níquel. A membrana utilizada, constituída de poliamida, modelo HR98PP e fornecida pela DOW/Filmtec, apresentou uma boa permeabilidade hidráulica. Os resultados mostraram que para todas as concentrações testadas, as rejeições de níquel ultrapassaram 96%, comprovando a boa seletividade deste tipo de membrana na rejeição do referido metal, com fluxos de permeado variando entre 4,78 e 5,55 L/h.m2 , sob pressão de operação de 10 bar. Para estudar o efeito do tamanho iônico na rejeição da membrana, o níquel foi complexado pela adição de um agente quelante na solução de alimentação. O agente escolhido foi o Na2EDTA, devido à formação de um complexo estável com o níquel e por ser um agente não prejudicial à saúde humana. Os resultados com adição de EDTA indicaram um aumento na rejeição de níquel, atingindo o índice máximo de 98,22 %, partindo-se de uma solução com 40,39 ppm de Ni2+, e confirmam que o processo de osmose inversa com a membrana HR98PP é altamente adequado para o tratamento de efluentes contendo níquel

Fixação de complexo de ferro II em matriz de poli(acrilato de sódio)

Neste trabalho estudou-se a complexação do poli(acrilato de sódio), PAS, com aminpentacianoferrato de sódio,APCF. As reações de complexação ocorreram imediatamente com a formação de uma coloração amarela com total solubilização do polímero na solução 0,032 M molar do complexo. As soluções mostraram-se estáveis por 48 h e após esse período observou-se alterações na natureza do complexo polimérico com a formação de precipitado. O UV foi usado como ferramenta de caracterização do complexo. O máximo de absorção obtido após dissolução imediata foi de 405 nm com desvio para o azul (398 nm) e um pequeno efeito hipercrômico. As amostras mantidas à temperatura ambiente por mais de 48 h deram origem a precipitados, que como a solução, absorveram com máximo de 364, 389 e 398 nm. O Complexo PAS-APCF foi também caracterizado através de FTIR por ATR e apresentou pequenas variações no espectro do material de partida (PAS). Um incremento na intensidade da deformação axial assimétrica do grupo carboxilato (1651 cm-1) e a presença do estiramento em 2055 cm-1 do grupo cianeto, diferentemente do APCF (2048 cm-1), confirmaram a formação do complexo PAS-APCF. As freqüências de absorção observadas para o complexo foram compatíveis com a presença de estruturas mono e bidentadas de complexação. As análises de TGA e DSC também foram utilizadas para a caracterização das estruturas. O estudo modelo envolvendo a complexação de sais sódicos de diácidos orgânicos de diferentes tamanhos de cadeia (oxalato, malonato, succinato, glutarato e adipato), diferentemente do PAS, promoveu um desvio para o vermelho na frequência máxima de absorção junto de um pequeno efeito hipercrômico (421 nm). Esta variação pode também ser observada quando do emprego de acetato de sódio, indicando, provavelmente, apenas a formação de estruturas monodentadas

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

As competências legislativas na Constituição de 1988: parâmetros para sua interpretação e para a solução de seus conflitos

O presente estudo procura examinar criticamente a forma como as competências legislativas são interpretadas no Brasil. Em especial, pretende-se demonstrar que o tema pode e deve se beneficiar das modernas técnicas e instrumentos desenvolvidos pela dogmática do Direito Constitucional contemporâneo. O trabalho se estrutura em três partes. Na Primeira Parte, serão expostas algumas premissas teóricas sobre a interpretação constitucional, o federalismo e a sindicabilidade judicial das competências, que nortearão o desenvolvimento do estudo. Na Segunda Parte, examinam-se os processos de qualificação das leis e de interpretação das competências legislativas. A partir do esboço de uma teoria das competências legislativas, será defendida a aplicação de parâmetros segundo os quais, em princípio, todos os dispositivos de competência devem ser interpretados da forma mais ampliativa possível, sendo as eventuais restrições, impostas por outras regras de competência, consideradas e justificadas argumentativamente. Em sua Terceira Parte, e última, o estudo identificará o fenômeno dos conflitos de competências legislativas em geral, esquecidos pela doutrina brasileira , examinando, na sequência, alguns critérios para sua solução. Afastada a possibilidade de recurso à supremacia do direito federal e ao princípio da subsidiariedade, bem como a preferências de mérito, serão desenvolvidos dois parâmetros formais e um material para a solução das inconsistências insolúveis entre competências.

Implementação de um algoritmo numérico para solução da equação de Christoffel generalizada em acustoelasticidade.

Extensos estudos realizados nas últimas décadas sobre a propagação de ondas ultrassônicas em sólidos levaram ao desenvolvimento de técnicas não destrutivas para a avaliação da segurança e integridade de estruturas e componentes industriais. O interesse na aplicação de técnicas ultrassônicas para medição de tensões aplicadas e residuais decorre da mudança mensurável da velocidade das ondas ultrassônicas na presença de um campo de tensões, fenômeno conhecido como efeito acustoelástico. Uma teoria de acustoelasticidade fornece um meio atrativo e não destrutivo de medir a tensão média ao longo do caminho percorrido pela onda. O estudo da propagação das ondas ultrassônicas em meios homogêneos anisotrópicos sob tensão conduz a um problema não linear de autovalores dado pela equação de Christoffel generalizada. A característica não linear deste problema decorre da interdependência entre as constantes elásticas efetivas do material e as tensões atuantes. A medição experimental de tensões por técnicas ultrassônicas é um problema inverso da acustoelasticidade. Esta dissertação apresenta a implementação de um algoritmo numérico, baseado no método proposto por Degtyar e Rokhlin, para solução do problema inverso da acustoelasticidade em sólidos ortotrópicos sujeitos a um estado plano de tensões. A solução da equação de Christoffel generalizada apresenta dificuldades de natureza numérica e prática. A estabilidade e a precisão do algoritmo desenvolvido, bem como a influência das incertezas na medição experimental das velocidades das ondas ultrassônicas, foram então investigadas. Dados sintéticos para as velocidades das ondas ultrassônicas de incidência oblíqua em uma placa sujeita a um estado plano de tensões foram gerados pela solução direta da equação de Christoffel generalizada para ilustrar a aplicação do algoritmo desenvolvido. O objetivo maior desta dissertação é a disponibilização de uma nova ferramenta de cálculo para suporte às atividades experimentais de medição de tensões por ultrassom no país.

Os conflitos tributários internacionais e sua possível solução pela via arbitral

O fenômeno da globalização teve o condão de aproximar os diversos povos, cada um com seus interesses e culturas próprios. A existência de um consenso internacional na definição de princípios a serem seguidos quando das relações externas não consegue impedir, contudo, o surgimento de possíveis conflitos e divergências, tendo em vista a pluralidade cultural das diversas nações mundiais, fato que induziu a sociedade internacional a desenvolver meios que pudessem dirimir pacificamente as controvérsias, porventura, surgidas entre elas. A adoção dos meios para solução pacífica dos conflitos internacionais encontra-se incentivada pela Constituição da República Federativa do Brasil (art. 4., incisos VI e VII) e sua utilização não importa qualquer renúncia ao exercício da soberania, nem tampouco à imunidade de jurisdição. Para que se tenha uma eficácia maior da submissão dos conflitos surgidos no âmbito externo aos meios admitidos para resolvê-los, é importante que os países envolvidos no litígio possuam orientação interna no sentido de privilegiar o Direito Internacional frente à sua legislação infraconstitucional doméstica. A eventual primazia do direito interno pode resultar na inocuidade da adoção dos meios pacíficos de solução de controvérsias internacionais, uma vez que as autoridades dos países litigantes poderão se esquivar do cumprimento do acordo ou decisão alegando uma possível contrariedade com os ditames legais domésticos. Nesse contexto, a seara tributária tem despertado constantes divergências internacionais. As diferentes interpretações conferidas pelas diversas nações, dentre elas o Brasil, quando da aplicação dos tratados por elas firmados e que tenham vertente fiscal, em especial aqueles que visam evitar a dupla imposição fiscal da renda, ou garantir o livre trânsito de bens, pessoas e serviços, acaba trazendo grande insegurança àqueles investidores que possuem operações conectadas a dois ou mais sistemas tributários diferentes. Assim, ganham cada vez mais corpo, os debates em torno da extensão dos mecanismos pacíficos para resolução de divergências, também ao âmbito de aplicação de todo e qualquer tratado que verse sobre a matéria tributária. Tal fato propicia a busca de uma possível uniformização dos métodos hermenêuticos aplicáveis àqueles fatos geradores tributáveis que se encontrem vinculados a dois ou mais entes soberanos. É nesse contexto que se apresenta o presente estudo, o qual aborda a possibilidade de a República Federativa do Brasil submeter ao procedimento arbitral aquelas controvérsias de cunho tributário que eventualmente decorram da interpretação divergente das convenções internacionais das quais seja parte e que tratem de matéria fiscal.

Avaliação da remoção do íon cádmio em solução aquosa utilizando processos de osmose inversa e nanofiltração

Os efluentes contendo cádmio podem ser gerados dos por várias indústrias. A aplicação de ambas as tecnologias osmose inversa (OI) e nanofiltração (NF) para o tratamento de efluentes contendo o íon cádmio para reduzir a degradação ambiental foram estudadas no presente trabalho. Amostras sintéticas contendo o íon cádmio em várias concentrações foram preparadas e submetidas aos tratamentos por OI e NF em escala laboratorial. Assim, dividiu-se o processo de remoção em três etapas: 1- analisou-se os processos de osmose inversa e nanofiltração em termos da rejeição do íon cádmio, onde se correlacionou e avaliou os principais parâmetros de processos (concentração da solução alimentação, fluxo de permeado e pressão) e comparou-se os sistemas entre si. Nas etapas 2 e 3 avaliou-se os efeitos da salinidade utilizando os contra-íons zinco e sódio respectivamente, sobre a rejeição do íon cádmio no sistema de osmose inversa em soluções aquosas contendo misturas de sulfato de cádmio e sulfato de zinco, e sulfato de cádmio e sulfato de sódio. Mediu-se as concentrações de cádmio e zinco por espectroscopia de absorção atômica e as concentrações de sulfato e sódio por cromatografia de íons. Os resultados obtidos na primeira etapa dos experimentos mostram que a osmose inversa e nanofiltração são eficientes na remoção de cádmio, onde se obteve as taxas médias de rejeição de 92,4% e 96,6% respectivamente, sendo esses valores iguais em magnitude, tendo em vista os valores da incerteza padrão relativa associada a esses resultados. Em ambos os sitemas, os fluxos de permeado e as rejeições não sofreram influência do aumento da concentração de cádmio na solução alimentação.Na etapa 2, os resultados mostraram a diminuição da rejeição de cádmio com o aumento da concentração de sulfato de zinco na alimentação, sendo esse fenômeno atribuído a semelhança existente entre os íons Cd2+ e Zn2+ associados aos seus raios iônicos de hidratação e as energias de hidratação. Foram observadas as quedas do fluxo de permeado e das rejeições com o aumento da concentração de alimentação provavelmente devido ao efeito da polarização. Na etapa 3, a rejeição do íon cádmio não foi afetada pelo aumento da concentração de sulfato de sódio na alimentação, assim como os fluxos de permeado. Outro aspecto relacionado à rejeição é que a difusão do sódio através da membrana é favorecida frente à difusão do cádmio, sendo atribuído ao menor raio de hidratação do sódio

Preparação e avaliação de resinas biocidas impregnadas com iodo a partir de resinas comerciais de estireno e divinilbenzeno

Neste trabalho, foram empregadas duas resinas comerciais de caráter fortemente básico. Ambas tendo como base copolímeros de estireno e divinilbenzeno (DVB), sendo que a VPOC 1950 contém em sua composição grupos quaternários de amônio do tipo 1, que apresentam três grupos metila e a VPOC 1960 grupos quaternários de amônio do tipo 2, que apresentam um grupo etanoíla e dois grupos metila. As resinas comerciais citadas foram escolhidas por apresentarem grande capacidade de troca iônica, estabilidade e grupos funcionais de interesse para a impregnação com iodo. As resinas foram impregnadas com iodo por meio de três metodologias distintas, uma solução metanólica de iodo 0,08 mol/L com e sem iodeto de potássio 0,14 mol/L, e iodo 0,08 mol/L em heptano. As resinas foram caracterizadas por área específica, volume do poro, grau de inchamento, microscopia ótica, espectrometria de infravermelho por transformada de Fourier (FTIR), análise elementar, termogravimetria, microscopia eletrônica de varredura e determinação do iodo fixado por iodometria. A avaliação da atividade biocida foi realizada através do método da contagem em placas, utilizando-se uma cepa de Escherichia coli ATCC11775 em concentrações de 103 a 107 células/mL. Todas as resinas impregnadas mostraram atividade bactericida significante devido à presença de iodo correlacionada às características da resina, tais como: grupos funcionais, tamanho e formato dos poros. Para efeito de comparação, foram realizados ensaios bactericidas com as resinas de partida para comprovação ou não da ação bactericida ser atribuída somente ao iodo

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Percepção sobre o papel do Código de Ética Profissional na solução de dilemas éticos: um estudo empírico na Controladoria Geral do Município do Rio de Janeiro (CGM)

O grande impacto dos escândalos que afetam organizações das mais diversas áreas, inclusive as públicas, tem difundido um interesse generalizado em comportamentos éticos e antiéticos. A administração ideal baseia-se em um modelo de gestão, cujo processo utiliza mecanismos que assegurem sua eficiência e eficácia e, conseqüentemente, propiciem a redução de conflitos dentro da organização. A observância de um código de ética profissional favorece essa redução de conflitos, pois seu objetivo central é a formação da consciência sobre padrões de conduta. Nesse contexto, os servidores públicos devem constituir e observar determinadas regras de conduta que norteiam seu exercício profissional, pois a consciência dos princípios morais, seja no exercício do cargo ou função, ou fora dele, refletirá no exercício vocacional do poder do estado. Seus atos, comportamentos e atitudes deverão preservar a honra e a tradição do Funcionalismo Público diante da sociedade. Entretanto, a efetividade do cumprimento de um Código de Ética não se baseia nas leis administrativas e nem com estas se confunde, mas se apóia no sentimento de adesão moral e de convicção íntima de cada servidor público. Diante dos desafios encontrados na gestão de organizações públicas e de sua importância para a sociedade em geral, o objetivo desta pesquisa é verificar os fatores que influenciam a percepção dos servidores públicos civis municipais quanto ao valor do Código de Ética Profissional da Prefeitura do Rio de Janeiro como auxiliar na solução de dilemas éticos. A amostra foi constituída por 90 (noventa) servidores lotados na Controladoria Geral do Município (CGM). A metodologia utilizada para a análise exploratória e confirmatória dos dados, necessitou de técnicas de análise fatorial na redução de dados e regressão linear múltipla para testar as cinco hipóteses oriundas do problema de pesquisa. Os resultados obtidos evidenciaram que a percepção dos servidores sobre os valores morais professados no Código de Ética e o conhecimento deles sobre o Código influenciam a percepção sobre o seu valor como auxiliar na resolução de dilemas éticos. Os resultados também sugerem que o servidor municipal percebe, em um nível baixo, o compromisso da Administração com o código.

Um algoritmo em paralelo para solução de equações diferenciais evolutivas

Este trabalho que envolve matemática aplicada e processamento paralelo: seu objetivo é avaliar uma estratégia de implementação em paralelo para algoritmos de diferenças finitas que aproximam a solução de equações diferenciais de evolução. A alternativa proposta é a substituição dos produtos matriz-vetor efetuados sequencialmente por multiplicações matriz-matriz aceleradas pelo método de Strassen em paralelo. O trabalho desenvolve testes visando verificar o ganho computacional relacionado a essa estratégia de paralelização, pois as aplicacações computacionais, que empregam a estratégia sequencial, possuem como característica o longo período de computação causado pelo grande volume de cálculo. Inclusive como alternativa, nós usamos o algoritmo em paralelo convencional para solução de algoritmos explícitos para solução de equações diferenciais parciais evolutivas no tempo. Portanto, de acordo com os resultados obtidos, nós observamos as características de cada estratégia em paralelo, tendo como principal objetivo diminuir o esforço computacional despendido.

Avaliação dos algoritmos de Picard-Krylov e Newton-Krylov na solução da equação de Richards

A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards.