Modelagem matemática da transferência de calor numa placa plana sob o efeito de uma fonte pontual externa de radiação térmica.

Este trabalho apresenta uma modelagem matemática para o processo de aquecimento de um corpo exposto a uma fonte pontual de radiação térmica. O resultado original que permite a solução exata de uma equação diferencial parcial não linear a partir de uma seqüência de problemas lineares também é apresentado. Gráficos gerados com resultados obtidos pelo método de diferenças finitas ilustram a solução do problema proposto.

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

Identificação de areias siliciclásticas para a recuperação de praias em erosão através de métodos geofísicos acústicos

Devido a sua grande extensão latitudinal, a costa brasileira é influenciada por diferentes regimes climáticos e oceanográficos. Adicionalmente, a distribuição da população brasileira é caracterizada por uma alta concentração nas capitais litorâneas. Todos esses fatores levam à construção de inúmeras estruturas de engenharia que podem de alguma forma impactar o transporte de sedimento e consequentemente o balanço sedimentar de algumas praias. Uma das formas mais eficientes de recuperar esse balanço sedimentar é a alimentação artificial de praias, com sedimentos provenientes da plataforma continental com características semelhantes. Os métodos geofísicos acústicos permitem mapear de forma eficiente o fundo e o subfundo marinho para a busca de áreas fontes para a lavra de forma eficiente. O objetivo deste estudo é a identificação de padrões geoacústicos a partir de analises quantitativas e qualitativas, para a caracterização de áreas fontes de areias siliciclásticas compatíveis com sedimentos de praias em erosão na cidade do Rio de Janeiro. O mapeamento foi realizado na plataforma continental interna do Rio de Janeiro, em área adjacente à praia de Itaipuaçu (Maricá, RJ). As análises evidenciaram depósitos, localizados entre 19 e 30 metros, de areias com granulometria, textura e selecionamento em condições ideias para a recomposição de praias do Rio de Janeiro.

Identificação de indicadores globais para o monitoramento da segurança de processos de plataformas de produção de petróleo e gás natural: estudo de caso da indústria brasileira

Este trabalho verificou o tema de segurança de processos no segmento de plataformas de produção de petróleo e gás no mar de forma a propor um modelo de indicadores globais a fim de facilitar a identificação de fragilidades no gerenciamento de riscos de processos da indústria de produção de petróleo e gás natural no mar. O resultado esperado é o fornecimento de dados para a otimização de recursos regulatórios na busca da melhoria da segurança operacional. Para tal, buscando conhecer a postura das empresas no cumprimento da necessidade regulamentar da implementação de um Sistema de Gestão de Segurança de Processos (SGSP), foram identificadas as principais práticas preventivas para a o gerenciamento de riscos. Por conseguinte, mensurou-se o esforço na implementação deste SGSP como correlação direta para a prevenção de acidentes (foco preventivo). Já a análise dos resultados deste gerenciamento de riscos (foco reativo) foi obtido através incidentes operacionais classificados de acordo com sua respectiva importância/impacto, dentro de linhas de corte pré-estabelecidas. Assim, o indicador Global preventivo (IGPRV) e o Indicador Global reativo (IGRTV) foram elaborados agregando diferentes indicadores específicos de cada segmento para simplificar a análise e incluir a articulação de conceitos de segurança de processos à atuação regulatória. Os indicadores globais propostos quantificam conceitos subjetivos e complexos para o estabelecimento de metas e permitem aferir o desempenho das ações e projetos da regulação da indústria. Também permitem a identificação e a disseminação das melhores práticas de gerenciamento de riscos, comparando-as em diferente níveis de gerenciamento. Na análise dos resultados, pôde-se perceber que a utilização dos indicadores globais propostos pode responder com sucesso às necessidades regulatórias identificadas

Método numérico de Matriz Resposta acoplado a um esquema de reconstrução espacial analítica para cálculos unidimensionais de transporte de nêutrons na formulação de ordenadas discretas multigrupo de energia com fonte fixa

Um método de Matriz Resposta (MR) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas multigrupo de transporte de nêutrons com fonte fixa e em geometria unidimensional na formulação de ordenadas discretas (SN). Portanto, o método multigrupo MR com esquema iterativo de inversão nodal parcial (NBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações multigrupo SN, afora os erros de arredondamento da aritmética finita computacional. É também desenvolvido um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em cada grupo de energia em um intervalo qualquer do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Gestão de Custos através do Custeio Baseado em Atividades (ABC): um estudo de caso em um projeto de desenvolvimento de software

O contexto da era da informação exige novas habilidades para competir com sucesso, tanto para empresas industriais como para as de serviços. Desde o seu surgimento, a Engenharia de Software tem enfrentado o desafio em estimar, controlar e gerenciar custos no processo de desenvolvimento de software. Os custos são parte de um projeto de software e nem sempre são claramente mensurados. Algumas desenvolvedoras de software são caracterizadas como prestadoras de serviços e as organizações que prestam serviços possuem as mesmas dificuldades em termos gerenciais que as indústrias, elas precisam de sistemas de gestão adequados para gerir seus custos, logo seus recursos. Nesse contexto é que os sistemas de informações contábeis estão inseridos, fornecendo respostas que os gestores necessitam para tomar decisões. Os custos são exemplos desse tipo de informação, pois o seu conhecimento é indispensável para uma boa gerência dos recursos organizacionais. O Custeio Baseado em Atividades (ABC) é uma ferramenta contábil útil para o gerenciamento estratégico dos custos uma vez que os sistemas de custeio tradicionais já não atendem essa necessidade das organizações. Diante desse cenário, o objetivo geral dessa dissertação é aplicar a metodologia de custeio baseado em atividades em um projeto de desenvolvimento de software com a finalidade de gestão de custos. Este pesquisa caracteriza-se como qualitativa descritiva e faz uso do método do estudo de caso. Através desse estudo foi possível comprovar a utilidade do ABC como instrumento de auxílio na gestão de custos nas organizações que prestam serviços em desenvolvimento de software, pois possibilita a visualização da maneira que os recursos são consumidos pelas atividades que agregam e não agregam valor ao negócio e ao cliente.

A estandardização dos contratos internacionais e o modelo FIDIC - Silver Book para contratos EPC

No presente trabalho acadêmico se apresentarão as principais características do processo de globalização, passando-se à exposição da relevância do Direito Comparado, até se chegar à análise da adoção de modelos internacionais de contratos como esforço uniformizador, de modo a reduzir riscos e a minorar custos de transação. Na sequência, investigar-se-á o modelo contratual EPC (acrônimo das palavras inglesas Engineering, Procurement and Construction - Engenharia, Gestão de Compra e Construção), originário da prática anglo-saxã, focando-se na análise e na qualificação tipológica do Contrato EPC, comparando-o com institutos existentes na legislação brasileira. Delinear-se-á o contexto de disseminação do Contrato EPC no exterior e sua consolidação no Brasil, amparadas, em larga medida, na necessidade de que - em vultosos projetos de infraestrutura, sobretudo em áreas de investimento em relação às quais os empreendedores desconheçam o ambiente regulatório e a realidade socioeconômica se tenha estatuto privado a transladar ao construtor a maior parte dos riscos atinentes a serviços complexos de engenharia. Enfrentar-se-ão as características essenciais deste modelo contratual, tomando-se como padrão o Conditions of Contract for EPC Turnkey Projects - general conditions, guidance for the preparation of particular conditions, forms of letter of tender, contract agreement and dispute adjudication agreement, recomendado pela Fédération Internationale des Ingénieurs Conseils - FIDIC. Apresentar-se-á como o direito estrangeiro equilibra a assunção dos riscos pelo construtor (Contractor ou Builder), inclusive aqueles, referentes a eventos não antecipáveis (unforeseen risks), inobstante a preservação pelo contratante (Owner) de lato poder de monitoramento e fiscalização (overseeing attributions, key personnel and contract manager approval, step-in rights).

Efeito da temperatura de secagem de dois tipos de silano e da supressão ou redução do tempo de condicionamento com ácido hidrofluorídrico na resistência de união entre uma cerâmica de dissilicato de lítio e um cimento resinoso

Este estudo avaliou o efeito de diferentes métodos de silanização e aplicação do ácido hidrofluorídrico (HF) sobre a resistência à microtração de uma cerâmica de dissilicato de lítio a um cimento resinoso. Quarenta blocos de IPS e.max Press /Ivoclar Vivadent (5x5x6mm) foram cimentados a blocos de resina Z250/3M ESPE (5x5x6mm) usando o cimento resinoso RelyX ARC/3M ESPE de acordo com os seguintes métodos de tratamento superficial: G1: 20s de ácido fluorídrico (HF) + silano não hidrolisado Primer-Activactor/Dentsply (SNH) seco à temperatura ambiente; G2: 20s HF + silano pré-hidrolisado RelyX Ceramic-Primer/3M ESPE (SPH) seco à temperatura ambiente; G3: 10s HF + SNH seco com ar quente (50oC-2min); G4: 10sHF + SPH seco com ar quente (50oC-2min); G5: sem ácido, SNH seco com ar quente (50oC-2min); G6: sem ácido, SPH seco com ar quente (50oC-2min); G7: sem ácido, SNH seco à temperatura ambiente; G8: sem ácido, SPH seco à temperatura ambiente. Antes de cada método de silanização, os blocos cerâmicos receberam acabamento com lixas de carbeto de silício (220-600) e limpeza com ácido fosfórico 37% (1min). A cimentação foi realizada com carga vertical de 1kg por 10min. Os conjuntos de cerâmica/cimento/resina foram armazenados em água destilada (37C) por 24 horas e depois seccionados em máquina de corte Isomet 1000 a fim de obter palitos (n = 40) de 1mm2 de área da seção transversal, que foram submetidos ao teste de microtração em máquina de ensaio universal Emic (v = 0,5mm/min). O modo de fratura foi avaliado em microscópio eletrônico de varredura. A análise estatística foi realizada utilizando ANOVA / Dunnett (p-valor = 0,000). As médias MPa e desvio padrão foram: G1-21,5 (8,9) BC; G2-30,5 (7,2) A; G3-19.4 (9.1) BC; G4-24,0 (9,0) B; G5-8.1 (3.2) D; G6 -18,0 (6,2) C; G7-7.8 (2,6) D; G8-6.3 (2,5) D. Grupos 2, 3, 4 e 6 não tiveram falhas prematuras dos palitos contra os grupos 1, 5, 7 e 8, que apresentaram 2,2; 44,4; 75,6 e 33,3% de perdas prematuras, respectivamente. O teste de correlação foi realizado apresentando significância estatística, com valor de -0,736 (p-valor = 0,000), mostrando que, a medida que o percentual de perda prematura aumenta, a média da MPa diminui. Quanto ao modo de fratura, observou-se 44,97% de falhas do tipo mista, 51,70% de falhas do tipo adesiva, 3,33% de falhas do tipo coesiva do cimento. Quando é realizada a supressão do condicionamento com HF como pré-tratamento da superfície cerâmica IPS e.max Press, a aplicação de silano SPH, associada ao seu tratamento térmico, deve ser o método de silanização recomendado, embora os valores mais elevados de resistência de união tenham sido os obtidos quando utilizado o condicionamento com HF por 20s. Quando é realizada a redução do tempo de condicionamento com HF para 10s, a aplicação do silano (SPH ou SNH) deve ser sempre associada ao seu tratamento térmico. O SNH só deve ser usado se as superfícies das cerâmicas IPS e.max Press forem tratadas com HF.

Otimização do ciclo de liofilização de uma vacina viral

A liofilização - ou secagem a frio (freeze drying em inglês) - é um complexo processo multiestágios, onde o produto é primeiramente congelado e sua secagem feita através de sublimação. Devido a esta forma de secagem a liofilização se torna um processo atrativo particularmente importante para a estabilização de componentes orgânicos hidratados e lábeis, de origem biológica. O processo de liofilização é amplamente empregado na indústria farmacêutica, porém em termos de gestão de processo, deve-se evitar a liofilização a todo custo, pois o processo possui diversas desvantagens como: equipamentos de alto investimento, alta demanda energética, processo que demanda tempos longos e produtos com facilidade de hidratar e frágeis, necessitando ser cuidadosamente embalados e armazenados. Este trabalho tem como objetivo a diminuição do ciclo de liofilização de uma vacina viral e analisar a possibilidade de carregamento desse produto no liofilizador a temperaturas negativas, de forma a possibilitar o aumento de produtividade. Para tal, foram realizados três experimentos com ciclos de liofilização com 17 e 20h a menos que o ciclo de liofilização da vacina comercial. Os experimentos foram realizados com a mesma formulação do lote comercial e utilizando um liofilizador piloto. As modificações foram realizadas nas propriedades físicas do ciclo de liofilização atual (temperatura, pressão e tempo) e na temperatura de carga do produto, sem alteração na formulação da vacina ou embalagem primária. Amostras do produto liofilizado experimental foram analisadas quanto ao seu aspecto, desempenho, umidade residual, potência e termoestabilidade acelerada segundo os Mínimos Requerimentos da Organização Mundial da Saúde. Todos os resultados analisados estiveram dentro das especificações e próximos ou melhores quando comparados aos lotes comerciais de origem

Estudo de pirólise catalítica de biomassa em escala piloto para melhoramento da qualidade do bio-óleo

A pirólise rápida é um processo para conversão térmica de uma biomassa sólida em altos rendimentos de um produto líquido chamado de bio-óleo. Uma das alternativas para geração de um bio-óleo com menor teor de oxigênio é uso de catalisadores nos reatores de pirólise, ao invés de um inerte, num processo chamado de pirólise catalítica. O objetivo deste trabalho foi testar catalisadores comerciais, um ácido e outro básico, em uma unidade piloto de leito fluidizado circulante. O catalisador ácido utilizado foi o Ecat, proveniente de uma unidade industrial de craqueamento catalítico fluido (FCC), e como catalisador básico foi utilizado uma hidrotalcita. Os resultados foram comparados com testes utilizando um material inerte, no caso uma sílica. Uma unidade piloto de FCC do CENPES foi adaptada para realizar os testes de pirólise catalítica. Após fase de modificação e testes de condicionamento, foi comprovada a viabilidade na utilização da unidade piloto adaptada. Contudo, devido a limitações operacionais, maiores tempos de residência tiveram que ser aplicados no reator, configurando o processo como pirólise intermediária. Foram então realizados testes com os três materiais nas temperaturas de 450C e 550C. Os resultados mostraram que o aumento do tempo de residência dos vapores de pirólise teve um impacto significativo nos rendimentos dos produtos quando comparada com o perfil encontrado na literatura para pirólise rápida, pois devido ao incremento das reações secundárias, produziu maiores rendimentos de coque e água, e menores rendimentos de bio-óleo. O Ecat e a hidrotalcita se apresentaram mais efetivos em termos de desoxigenação. O primeiro apresentou maiores taxas de desoxigenação via desidratação e a hidrotalcita apresentou maior capacidade para descarboxilação. Contudo, o uso de Ecat e hidrotalcita não se mostrou adequado para uso em reatores de pirólise intermediária, pois acentuou ainda mais as reações secundárias, gerando um produto com alto teor de água e baixo teor de compostos orgânicos no bio-óleo, além de produzirem mais coque. À temperatura de 450C estes efeitos foram mais pronunciados. Em termos de caracterização química, a condição de pirólise intermediária apontou para a produção de bio-óleos com perfil fenólico, sendo a sílica o que proporcionou os melhores rendimentos, principalmente a temperatura de 550C, sendo superiores aos encontrados na literatura. Analisando as composições dos bio-óleos sob a ótica da produção de biocombustíveis, nenhum dos materiais testados apresentou rendimentos consideráveis em hidrocarbonetos. De maneira geral, a sílica foi o que proporcionou os melhores resultados em termos de rendimento e qualidade do bio-óleo. Sua menor área superficial e sua característica de inerte se mostraram mais adequados para o processo de pirólise intermediária, onde a contribuição das reações secundárias em fase gasosa é elevada em função do tempo de residência no reator

Desenvolvimento de técnicas de otimização da programação das limpezas de redes de trocadores de calor

Trocadores de calor são equipamentos muito utilizados na indústria de processos com o objetivo de modificar a temperatura e/ou o estado físico de correntes materiais. Uma rede de trocadores de calor pode ser definida como um grupo de trocadores de calor interligados, a fim de reduzir as necessidades de energia de um sistema. No entanto, durante a operação de uma rede, a eficiência térmica dos trocadores de calor diminui devido à deposição. Esse efeito promove o aumento dos custos de combustível e das emissões de carbono. Uma alternativa para mitigar este problema baseia-se no estabelecimento de uma programação das limpezas dos trocadores de calor durante a operação de rede. Este tipo de abordagem ocasiona uma situação na qual ocorre um conflito de escolha: a limpeza de um trocador de calor pode recuperar a sua eficiência térmica, mas implica custos adicionais, tais como, mão-de-obra, produtos químicos, etc. Além disso, durante a limpeza, o trocador de calor tem de ser contornado por uma corrente de by-pass, o que aumenta temporariamente o consumo de energia. Neste contexto, o presente trabalho tem como objetivo explorar diferentes técnicas de otimização envolvendo métodos estocásticos e heurísticos. Com este objetivo foi desenvolvido um conjunto de códigos computacionais integrados que envolvem a simulação pseudo-estacionária do comportamento da rede relacionado com incrustações e a otimização da programação das limpezas deste tipo de sistema. A solução do problema indica os períodos de tempo para a limpeza de cada trocador de calor. Na abordagem estocástica empregada, os parâmetros do algoritmo genético, como probabilidade de crossover e probabilidade de mutação, foram calibrados para o presente problema. A abordagem heurística desenvolvida se deu através da sequência do conjunto de movimentos zero, um e dois. De forma alternativa, desenvolveu-se a metodologia heurística recursiva na qual os conjuntos de movimentos um e dois foram empregados recursivamente. Também foi desenvolvida a abordagem híbrida que consistiu em diferentes combinações da metodologia estocástica e heurística. A análise comparativa entre as metodologias empregadas teve como objetivo avaliar a abordagem mais adequada para o presente problema da programação das limpezas em termos de função objetivo e esforço computacional. O desempenho da abordagem proposta foi explorado através de uma série de exemplos, incluindo uma refinaria real brasileira. Os resultados foram promissores, indicando que as técnicas de otimização analisadas neste trabalho podem ser abordagens interessantes para operações que envolvam redes de trocadores de calor. Dentre as abordagens de otimização analisadas, a metodologia heurística desenvolvida no presente trabalho apresentou os melhores resultados se mostrando competitiva frente às abordagens comparadas da literatura

Efeito de surfactantes na estimativa da densidade bacteriana em amostras de petróleo

Foi examinado o efeito dos surfactantes polisorbato 60 (Tween 60), polisorbato 80 (Tween 80), brometo de cetil trimetil amônio (CTAB) e lauril sulfato de sódio (SDS) na estimativa da densidade de Bactérias Redutoras de Sulfato (BRS) e Bactérias Anaeróbias Heterotróficas Totais (BANHT) em amostras de petróleo. Para a realização dos experimentos, foram selecionadas três amostras com diferentes proporções de óleo e água de forma a representar amostras reais. A primeira amostra contém uma alta proporção de óleo, a segunda uma proporção média e a última amostra uma baixa proporção de óleo. A densidade bacteriana foi estimada através do método do Número Mais Provável (NMP). As concentrações dos surfactantes empregadas neste estudo foram estabelecidas através de estudo anterior. Os resultados demonstram que nas amostras com alta e média proporção de óleo, a adição dos surfactantes não foi favorável a um aumento na quantificação de BRS. Por outro lado, o Tween 60 e o Tween 80 mostraram um aumento significativo na quantificação de BANHT quando aplicados na concentração de 0,01% e 0,03% m/v, respectivamente. O CTAB favoreceu o crescimento de BRS e BANHT na amostra com baixa proporção de óleo quando sua concentração foi de 0,001% m/v e 0,0001% m/v, respectivamente

Medição de espessura de filmes eletrodepositados em suporte metálico por EDXRF

Esse trabalho apresenta a técnica de fluorescência de raios X por dispersão de energia (EDXRF) para calcular a espessura de filmes metálicos finos por ser de custo relativamente baixo, não destrutiva e que permite também uma análise qualitativa e até quantitativa dos elementos do revestimento. As espessuras do Níquel, Zinco e Cromo sobre aço, foram calculadas utilizando cinco procedimentos diferentes e todas as relações matemáticas obtidas foram baseadas na lei de atenuação de Beer-Lambert, que relaciona as intensidades Kα e/ou Kβ e as relações específicas do material (coeficiente de atenuação, μm). Os resultados experimentais para as espessuras do Zinco e Níquel, utilizando somente a intensidade Kα, apresentaram valores próximos dos valores estimados. Quando foi considerada a razão Kα/Kβ, os resultados experimentais divergiram bastante do estimado. Para calcular a espessura do Cromo foi necessário utilizar um método em que as camadas são calculadas separadamente, pois o revestimento de Cromo é composto por múltiplas camadas (Cromo, Cobre e Níquel), Os resultados para a espessura do Cromo foram satisfatórios.

Monitoramento microbiológico e físico-químico de tanques de armazenamento de óleo e água

A injeção da água do mar nos campos marítimos (offshore), processo este conhecido como recuperação secundária de petróleo, gera muitos resíduos e efluentes. Dentre estes, pode-se destacar a água produzida, que consiste de água de formação, água naturalmente presente na formação geológica do reservatório de petróleo, e água de injeção, aquela normalmente injetada no reservatório para aumento de produção. Sete tanques de armazenamento de água/óleo de um terminal foram monitorados quanto à presença de micro-organismos e teores de sulfato, sulfeto, pH e condutividade. Particularmente, as bactérias redutoras de sulfato (BRS), que agem às expensas da atividade de outras espécies, reduzindo sulfato à sulfeto, constituindo-se num problema-chave. Os tanques de óleo codificados como Verde, Ciano, Roxo, Cinza, Vermelho, Amarelo e Azul, apresentaram comportamentos distintos quanto aos parâmetros microbiológicos e físico-químicos. Após este monitoramento, de acordo com valores referência adotados, e levando-se em conta como principais parâmetros classificatórios concentrações de BRS, bactérias anaeróbias totais e sulfeto, os dois tanques considerados mais limpos do monitoramento foram os tanques roxo e ciano. Analogamente, por apresentarem os piores desempenhos frente aos três principais parâmetros, os tanques amarelo e cinza foram considerados os mais sujos de todo o monitoramento. Após esta segregação, esses três principais parâmetros, mais a concentração de sulfato, foram inter-relacionados a fim de se corroborar esta classificação. Foi possível observar que o sulfeto instantâneo não foi o parâmetro mais adequado para se avaliar o potencial metabólico de uma amostra. Por este motivo, foram verificados os perfis metabólicos das BRS presentes nas amostras, confirmando a segregação dos tanques, baseada em parâmetros em batelada

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.