Um estudo de percolação em barragens de terra, em regimes permanente e transiente, com a aplicação do método de elementos finitos

Este trabalho apresenta um estudo de fluxo de água em barragens de terra, em regimes permanente e transiente, com a utilização do Método de Elementos Finitos. No estudo de fluxo em regime permanente duas formas de abordar o problema são apresentadas e comparadas. A primeira considera, para a discretização da malha de elementos finitos, somente a região saturada, de maneira que a linha freática é obtida através de ajustes desta malha de elementos finitos. A segunda considera toda a região saturada-insaturada, sendo discretizado todo o domínio físico da barragem. A malha de elementos finitos não é modificada ao longo das iterações e a linha freática é obtida por interpolação dentro dos elementos, em função dos valores nodais do potencial de pressões. O desenvolvimento teórico das equações utilizadas para as duas formas de abardagem é apresentado, mostrando onde elas diferem entre si. No estudo de fluxo em regime transiente é utilizado apenas o esquema de malha fixa de elementos finitos.

Um modelo elasto-viscoplástico para análise de peças de concreto estrutural, submetidas a estados planos de tensão, através do método dos elementos finitos

O objetivo deste trabalho é desenvolver um modelo computacional, baseado no método dos elementos finitos, para o estudo de peças de concreto armado e protendido submetidas a estados planos de tensão. O estudo abrange situações de carga de curta e longa duração, onde consideram-se fluência e retração do concreto e relaxação do aço. São utilizados modelos constitutivos elasto-viscoplásticos para descrever o comportamento dos materiais. Implementou-se um modelo de camadas superpostas para melhor representar o comportamento do concreto, onde o material é composto de diversas camadas que sofrem a mesma deformação. Cada camada possui diferentes características materiais e a tensão total é obtida pela soma das diferentes contribuições de cada camada. Para a fissuração da concreto, utilizou-se um modelo de fissuras distribuídas, que leva em conta a contribuição do concreto entre fissuras. Tanto a amadura passiva como a de pratensão são introduzidas no modelo como uma linha de material mais rígido dentro do elemento de concreto. Os deslocamentos ao longo da armadura são referenciados aos deslocamentos nodais do elemento de concreto. Deste modo, obtém-se uma matriz de rigidez para a armadura com as mesmas dimensões que a matriz de rigidez do elemento de concreto, A matriz de rigidez do elemento concreto-aço é a soma das duas matrizes. Considera-se aderência perfeita entre o concreto e o aço. Os resultados obtidos com esse programa computacionai são comparados com valores experimentais disponíveis.

Simulação numérica de escoamentos hipersônicos em não-equilíbrio termo-químico através do método dos elementos finitos

o presente trabalho tem por objetivo desenvolver um código computacional destinado à simulaçáo numérica tridimensional e de escoamentos compressíveis hipersônicos em especial para analisar o processo de reentrada de veículos espaciais na atmosfera terrestre. Devido ás temperaturas elevadas encontradas neste tipo de escoamento, efeitos de altas temperaturas sáo incluídos, como excitação vibracional das moléculas e reações químícas. Utiliza-se o modelo de 2 temperaturas (translacional e vibracional);5 espécíes (N2, O2, NO, N, O) e 17 reações químicas de Park [1986 para o modelamento do fenômeno termo-químico. Introduz-se uma formulação baseada no método dos elementos finitos, através do esquema de Taylor-Calerkin, para a discretização do tempo e do espaço e para a solução do sistema de equações. Alguns exemplos e aplicações de problemas hipersônicos bi e tridimensionais, difusivos e não difusivos, são apresentados, incluindo comparações entre os resultados obtidos no presente trabalho e na literatura.

Estudo sobre modelos de fissuração de peças de concreto armado via método dos elementos finitos

o objetivo desta tese é apresentar dois modelos distintos para representar as fissuras em peças de concreto armado: um do tipo distribuído e outro do tipo incorporado. Os modelos de fissura incorporada se baseiam no conceito de descontinuidades incorporadas dentro do campo de deslocamento do elemento finito padrão. Já nos modelos de fissura distribuída a descontinuidade do campo de deslocamentos causada pela fissura é espalhada ao longo do elemento. O modelo distribuído proposto apresenta um diagrama tensão-deformação multilinear e sua calibragem é feita, segundo o Código Modelo CEB-FIP 1990, através de ensaios de tirantes de concreto armado. O modelo incorporado implementado é baseado no modelo de Dvorlcin, CuiMo e Gioia, o qual não inclui a contribuição da armadura no equilíbrio interno de forças do elemento. A inclusão da parcela da armadura é feita através de um modelo de transferência de tensão por aderência, conforme Russo, Zingone e Romano, Russo e Romano e FIE - Bulletin 10. Para representar o comportamento do concreto intacto, utiliza-se o modelo constitutivo de Ottosen. Trata-se de um modelo elástico não-linear, tridimensional, que utiliza valores secantes dos parâmetros do material. Já para simular o comportamento das barras de aço da armadura, emprega-se o modelo incorporado desenvolvido por Elwi e Hrudey. Neste modelo, permite-se uma disposição arbitrária das barras de aço no interior dos elementos de concreto. O modelo constitutivo adotado para a armadura é do tipo elasto-plástico com endurecimento. Por fim, alguns exemplos numéricos são analisados com o objetivo de comprovar a eficácia dos dois modelos propostos.

Análise de pilares esbeltos de concreto armado confinados através do método dos elementos finitos

o confinamento é uma técnica que tem sido cada vez mais aplicada no reforço e reabilitação de peças comprimidas de concreto, com o objetivo de aumentar sua capacidade resistente e recuperá-las de eventuais deteriorações. Este trabalho apresenta um modelo não-linear fisico e geométrico, baseado no Método dos Elementos Finitos, para a análise de pilares esbeltos de concreto armado, sob estado triplo de tensões. Neste estudo, são utilizados elementos finitos isoparamétricos tridimensionais para o concreto e um modelo incorporado de armadura. O comportamento do concreto é descrito por um modelo elasto-viscoplástico, que permite a análise da estrutura sob cargas de curta e longa duração. Para representação do concreto ao longo do tempo, utiliza-se um modelo de camadas superpostas. A fissuração do concreto é representada através de um modelo de fissuras distribuídas, onde se considera a contribuição do concreto entre fissuras. É admitida aderência perfeita entre os materiais. O aço é modelado como um material elastoplástico perfeito ou com endurecimento linear, e os materiais utilizados para confinamento, como os compósitos de fibras de carbono, são modelados como um caso particular do modelo elastoplástico definido para o aço. O modelo não-linear geométrico foi desenvolvido com base na formulação Lagrangeana Total. Consideram-se grandes deslocamentos e pequenas deformações. Os resultados obtidos através do modelo computacional são cOnITontados com resultados experimentais disponíveis na literatura, de modo a validar o modelo matemático e a metodologia numérica. Boa correlação é obtida entre os resultados numéricos e experimentais, que confirmam estudos prévios realizados por diversos autores, na medida em que são verificados significativos ganhos de resistência e/ou ductilidade em pilares comprimidos de concreto, devido ao confinamento.

Análise da interação solo-estrutura em fundações de concreto armado pelo método dos elementos finitos

A análise da iteração solo-estrutura em fundações é um importante campo de pesquisa que ainda tem um grande progresso a ser feito. No presente trabalho foi desenvolvido um programa computacional para a análise da interação solo-estrutura de fundações de concreto armado. Este tema abrange duas áreas da engenharia civil: estruturas e geotecnia. O método dos elementos finitos foi usado no trabalho na seqüência para resolver o problema considerando estado plano de defonnação e comportamento elastoplásti.co dos materiais estudados (solo, concreto e aço). A linguagem de programação MATLAB foi usada em toda esta pesquisa como alternativa ao FORTRAN. O MATLAB foi escolhido uma vez que é uma linguagem de programação que permite facilmente construir uma interfàce de pré e pósprocessamento amigável. Os passos para a solução completa do problema foram os seguintes: Primeiramente um programa foi desenvolvido considerando o comportamento elastoplástico com critérios de plastificação e ruptura específicos para o concreto e solo. Soluções analíticas fechadas foram usadas para checar a precisão do programa. O segundo passo foi a introdução do reforço de aço no concreto por meio de um modelo para armaduras. Logo após, um modelo de fissuras para o concreto 1racionado foi in1roduzido no programa. Na seqüência o programa de pré e pós-processamento foi desenvolvido para gerar a malha de elementos finitos (pré-processamento), distribuição tensões e deformações, mapa de fissuras, etc (pósprocessamento). Finalmente, os parâme1ros constitutivos do solo, concreto e aço foram calibrados e várias situações reais de interação do solo-concreto de fundações de concreto armado foram simuladas. Nesta dissertação são encontrados resultados para as pressões de contato sapata-solo. Diferentes diagramas de tensões de interfàce foram obtidos em função rigidez relativa do elemento estrutural de concreto armado-solo. Na análise mnnérica, rigidez relativa desempenhou uma relevante função no comportamento mecânico do elemento estrutural de concreto armado (sapata) e da base de assentamento (solo), uma vez ruptura em ambos os casos esteve diretamente relacionada a esta grandeza. São encon1rados, em função da rigidez relativa, resultados indicativos dos modos de falha da fundação, excessiva plastificação do solo em fundações com rigidez relativa alta, e a plastificação armaduras, esmagamento do concreto, formação de fissuras, bielas e confinamento concreto para fundações de rigidez relativa baixa. Na análise numérica, obteve-se resultados importantes com relação ao projeto de fundações. Estes resultados foram cOnITontadoscom normas, destacando-se as discordâncias com relação às recomendações da nonna brasileira Projeto e Execução de Fundações" NBR-6122 (1996) para os diagramas de tensões interface sapata-solo usados no dimensionamento de fundações de concreto armado.

Análise de estruturas de concreto armado sujeitas à corrosão de armadura por cloretos através do método dos elementos finitos

A corrosão da annadura é um dos principais problemas para um grande número de es1ruturas de concreto annado por causa da redução substancial da área de aço, da perda de aderência, da fonnação de fissuras e lascamento do concreto. O modelo tridimensional proposto por este trabalho é baseado no processo de corrosão iniciado por ataque de cloretos e propagado a partir de micropilhas de corrosão, utilizando o método dos elementos finitos. A influência da temperatura, da umidade relativa do concreto, da ligação dos cloretos com compostos do concreto, da fonnação de fissuras provenientes do carregamento e da resistência ôhmica do concreto são levados em consideração. O modelo é dividido em 4 módulos que interagem entre si, sendo eles: análise estrutural, transferência de calor, difusão de cloretos e corrosão eletroquímica. O modelo computacional para a análise estrutural utiliza dois procedimentos. No primeiro busca-se a resposta da estrutura para um carregamento instantâneo, considerando-se o material com comportamento elasto-viscoplástico, associado com modelos de annadura incorporada e fissuras distribuídas.No segundo, detennina-se a resposta da estrutura ao longo do tempo com a teoria da solidificação proposta por Bazant e Prasannan (1989), utilizando as cadeias Kelvin, apresentando ótimos resultados com uma grande economia computacional. O processo de transferência de calor é baseado na fonnulação tridimensional da equação de Helmholtz, com condições de contorno de fluxos provenientes da radiação solar, convecção e radiação, apresentando bons resultados em relação a fonnulações analíticas. Para análise de vigas com a variação da temperatura média diária, constatou-se que as mesmas se encontram em equilíbrio ténnico com o ambiente A análise de difusão de cloretos implementada apresentou um comportamento idêntico ao da segunda lei de Fick para casos unidimensionais. Para análises de vigas carregadas se constatou que as fissuras provenientes de cargas de serviço não influenciam o processo de difusão e que o coeficiente de difusão referencial é o principal parâmetro na análise. O modelo de corrosão eletroquímica é baseado na fonnulação proposta por Liu (1996), obtendo-se como resultado a taxa de corrosão em função do tempo.

Modelo da ligação entre concreto e armadura na análise de estruturas de concreto pelo método dos elementos finitos

Este trabalho apresenta um modelo computacional para análise tridimensional de peças de concreto armado, utilizando o método dos elementos finitos. São utilizados modelos constitutivos elasto-viscoplásticos para representar o comportamento dos materiais, quando submetidos a cargas de curta ou longa duração. Se a estrutura está submetida a cargas de curta duração, o modelo desenvolvido fornece, quando alcançado o estado estável, a solução do problema elastoplástico. Em contrapartida, se a estrutura está submetida a cargas de longa duração, o modelo pode representar o comportamento viscoelástico com envelhecimento do concreto. Foi utilizado um modelo de camadas superpostas para assegurar uma representação adequada do comportamento real do concreto ao longo do tempo. Deste modo, admite-se que o concreto é constituído por um número conveniente de camadas, que apresentam a mesma deformação. Para a fissuração do concreto, utilizou-se um modelo de fissuras distribuídas, que leva em consideração a contribuição do concreto entre fissuras O objetivo deste estudo é a modelagem da armadura e do reforço nas peças de concreto armado, através de um modelo de aderência entre o concreto e o material adjacente. A armadura é introduzida no modelo como uma linha de material mais rígido dentro do elemento de concreto. E, a consideração da degradação da aderência é realizada através de um modelo onde os efeitos da aderência imperfeita são incluídos pela introdução de graus de liberdade para os deslocamentos relativos entre os materiais. Os resultados obtidos com este programa computacional são comparados com valores experimentais disponíveis na literatura, visando validar o modelo matemático e a metodologia numérica.

Modelagem de situações de punção em lajes de concreto armado através do método dos elementos finitos

O objetivo deste trabalho é o desenvolvimento de um programa computacional, elaborado através do método dos elementos finitos, para análise tridimensional de situações de punção em lajes de concreto armado com e sem armadura de cisalhamento. São utilizados modelos constitutivos elasto-viscoplásticos para representar o comportamento do concreto. Para a simulação do concreto submetido a um carregamento de longa duração, empregou-se o modelo de camadas superpostas. Este modelo admite que o concreto é constituído por um número conveniente de camadas com a mesma deformação, porém características materiais diferentes. A tensão total é obtida pela soma das contribuições de cada camada. A armadura é introduzida no modelo como uma linha de material mais rígido no interior do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Assim, os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, a matriz de rigidez da armadura tem a mesma dimensão da matriz de rigidez do elemento de concreto, as quais são somadas para a obtenção da matriz de rigidez total Para representar o comportamento após à fissuração do concreto, é utilizado o modelo de fissuras distribuídas, onde apenas as propriedades dos materiais são modificadas e a malha de elementos finitos não necessita de atualização. Sugerem-se duas formulações para a consideração da contribuição do concreto entre fissuras e da redução na capacidade de transferência de corte devido ao aumento da abertura da fissura. Este último é incluído através de uma simplificação, que consiste em atribuir ao módulo de elasticidade transversal, correspondente ao plano fissurado, um valor reduzido. Para comprovar a validade do programa computacional desenvolvido, comparam-se os resultados numéricos com valores experimentais de ensaios realizados por Coelho (1999), Trautwein (2001) e Fusco (1988). Obtiveram-se boas aproximações para lajes com e sem armadura de cisalhamento. Desta forma, o programa computacional desenvolvido pode ser empregado para generalizar resultados experimentais em lajes cogumelo com variação na resistência do concreto, na bitola das barras de armadura, no espaçamento entre as barras e no número de camadas de armadura de cisalhamento.

Aplicação do método dos elementos finitos para a análise elasto-viscoplástica de peças de concreto armado e protendido, submetidas a estados planos de tensão

O objetivo deste trabalho é apresentar o desenvolvimento de um modelo computacional, baseado no método dos elementos finitos, para o estudo de estruturas de concreto armado e protendido, sob estados planos de tensão. O modelo computacional utiliza dois procedimentos distintos, baseados em um modelo elasto-viscoplástico. No primeiro, determina-se a resposta da estrutura ao longo do tempo. Considera-se, nesta etapa, que o material possui comportamento viscoelástico. No segundo, busca-se a resposta da estrutura para um carregamento instantâneo, considerando-se o material com um comportamento elastoplástico. Um modelo de camadas superpostas é adotado para melhor representar o comportamento viscoelástico do concreto e da armadura protendida, onde estes materiais são compostos de diversas camadas que sofrem a mesma deformação. Cada camada possui diferentes propriedades mecânicas e a tensão total é obtida pela soma das diferentes contribuições de cada camada. Para representar o concreto fissurado, utiliza-se um modelo de fissuras distribuídas, que leva em conta a contribuição do concreto entre fissuras. Tanto a armadura passiva quanto a armadura de protensão são introduzidas no modelo como uma linha de material mais rígido dentro do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, as matrizes de rigidez das armaduras têm as mesmas dimensões das matrizes de rigidez dos elementos de concreto e ambas são somadas para a obtenção das matrizes de rigidez totais. Para verificar a precisão do modelo computacional, compararam-se resultados numéricos com valores experimentais disponíveis.

Análise numérica de escoamentos turbulentos tridimensionais empregando o método de elementos finitos e simulação de grandes escalas

O presente trabalho apresenta o estudo e implementação de um algoritmo numérico para análise de escoamentos turbulentos, tridimensionais, transientes, incompressíveis e isotérmicos, através da Simulação de Grande Escalas, empregando o Método de Elementos Finitos. A modelagem matemática do problema baseia-se nas equações de conservação de massa e quantidade de movimento de um fluido quase-incompressível. Adota-se um esquema de Taylor-Galerkin, com integração reduzida e fórmulas analíticas das funções de interpolação, para o elemento hexaédrico de oito nós, com funções lineares para as componentes de velocidade e constante no elemento para a pressão. Para abordar o problema da turbulência, emprega-se a Simulação de Grandes Escalas, com modelo para escalas inferiores à resolução da malha. Foram implementados o modelo clássico de Smagorinsky e o modelo dinâmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, é proposta e empregada, para o processo de segunda filtragem do modelo dinâmico. O esquema, que utiliza elementos finitos independentes envolvendo cada nó da malha original, apresentou bons resultados com um baixo custo computacional adicional. São apresentados resultados para problemas clássicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para análise de escoamentos caracterizados por elevados números de Reynolds, é discutida no capítulo final. São apresentadas sugestões para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para análise de escoamentos tridimensionais, turbulentos e transientes .

Análise interativa de situações de reforço e recuperação de peças de concreto armado, empregando o método dos elementos finitos

O objetivo deste trabalho é desenvolver um programa computacional, baseado no método dos elementos finitos, para simular situações de reforço e recuperação de peças de concreto armado. A linguagem Matlab é a ferramenta utilizada para a elaboração do programa. É uma linguagem de alta performance para computação técnica. Integra computação, visualização e programação em um fácil ambiente para uso, onde problemas e soluções são expressos em familiar notação matemática. A característica principal deste programa é de permitir alterações na modelagem numérica durante a análise do problema, podendo-se retirar ou acrescentar elementos da estrutura, aumentar ou diminuir rigidezes, alterar materiais, etc, viabilizando-se, assim, uma avaliação das diversas etapas ou hipóteses dos procedimentos de recuperação ou reforço de estruturas. O programa permite a mudança de vinculações do elemento estrutural, durante a análise. Além disto, é permitido não só incrementos de forças como incrementos de deslocamentos impostos à estrutura. O programa computacional utiliza duas etapas distintas, baseadas em um modelo elasto-viscoplástico. Na primeira etapa, faz-se a determinação da resposta da estrutura ao longo do tempo. Considera-se, nesta etapa, que o material possui comportamento viscoelástico. Na segunda, busca-se determinar a resposta da estrutura para um carregamento instantâneo, considerando-se o material com comportamento elastoplástico Para melhor representar o comportamento do concreto, quando submetido a carregamento de longa duração, utilizou-se um modelo de camadas superpostas. A armadura é introduzida no modelo como uma linha de material mais rígido dentro do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Assim, os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Para verificar a precisão do programa computacional desenvolvido, comparam-se os resultados numéricos com resultados experimentais, cujos ensaios foram realizados no Laboratório de Ensaios e Modelos Estruturais do Programa de Pós-Graduação em Engenharia Civil da Escola de Engenharia da Universidade Federal do Rio Grande do Sul, e no Laboratório de Estruturas da Escola de Engenharia de São Carlos- USP. Nessas comparações, obtiveram-se ótimas aproximações entre os resultados numéricos e experimentais.

Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.

Análise de estruturas compostas por aduelas pré-moldadas de concreto com protenção externa, sob estado plano de tensão, através do método dos elementos finitos

Aduelas protendidas externamente utilizadas em pontes contrastam com as clássicas construções monolíticas, consistindo de “pequenos” segmentos pré-moldados protendidos por cordoalhas externas. Há muitas vantagens na utilização deste tipo de composição, como rapidez e versatilidade na construção, alto controle de qualidade e economia. No presente trabalho, discute-se uma formulação, baseada no método dos elementos finitos, para simular o comportamento de estruturas compostas por aduelas externamente protendidas. Esta formulação permite tanto a análise de situações de serviço, onde as seções permanecem completamente comprimidas, como também de situações últimas, onde se considera a abertura das juntas e a transferência de esforços através das mesmas. O modelo computacional utiliza dois procedimentos distintos, baseados em um modelo elasto-viscoplástico. No primeiro, determina-se a resposta da estrutura ao longo do tempo, considerando-se, nesta etapa, que o material possui comportamento viscoelástico. No segundo, busca-se a resposta da estrutura para um carregamento instantâneo, considerando-se o material com um comportamento elastoplástico. As juntas entre as aduelas são representadas através de elementos de contato unidimensionais. Os nós de cada elemento de contato coincidem com os nós do elemento de concreto que a ele foi associado.Um modelo de camadas superpostas é utilizado para representar o comportamento viscoelástico do concreto e da armadura protendida Para representar o concreto fissurado, utiliza-se um modelo de fissuras distribuídas, que leva em conta a contribuição do concreto entre fissuras. A armadura passiva é introduzida na análise através de um modelo incorporado e a armadura de protensão não aderente é disposta em elementos fictícios, paralelos aos elementos reais da estrutura. Para verificar a adequação do modelo computacional, compararam-se resultados numéricos com valores experimentais disponíveis na literatura.

Estudo de um modelo convectivo-difusivo-reativo em combustão no método de elementos finitos

Neste trabalho, apresenta-se um estudo numérico de um modelo convectivo-difusivo-reativo em combustão baseado no Método de Elementos Finitos. Primeiramente, apresenta-se o desenvolvimento das equações de balanço (quantidade de movimento, massa, espécie e energia) que modelam um processo de mistura molecular e reação química, irreversível, de passo único e exotérmica entre duas espécies químicas F (Combustível) e O (Oxidante). Tais espécies reagem e formam um produto P, conforme vFF +vOO ! vPP + calor, onde vF , vO e vP são os coeficientes estequiométricos molares. No modelo, considera-se que a reação é de primeira ordem com respeito a cada um dos reagentes e que a taxa de reação específica segue a cinética de Arrhenius. Em seguida, o modelo é estudado numericamente considerando-se um domínio retangular e condições de contorno do tipo Neumann. Tanto a Técnica das Diferenças Finitas como a Técnica de Elementos Finitos são utilizadas na discretização espacial das equações do modelo. Para a integração no tempo, utiliza-se a método de Runge-Kutta simplificado de três estágios. Os diferentes códigos computacionais obtidos, tanto pela Técnica de Diferenças Finitas como de Elementos Finitos, são comparados frente ao problema de interesse. Observa-se que ambas as técnicas apresentam resultados equivalentes. Além disso, os códigos desenvolvidos são robustos (capazes de lidar com vários conjuntos de parâmetros), de baixo custo e precisos. Por fim, apresenta-se uma revisão do trabalho de Zavaleta [48], no qual obtem-se uma estimativa local do erro na aproximação do problema estudado pela Técnica de Elementos Finitos.