Com o objetivo de aumentar o lucro de plantas químicas, a Otimização em Tempo Real (RTO) é uma ferramenta que busca determinar as condições ótimas operacionais do processo em estado estacionário, respeitando as restrições operacionais estabelecidas. Neste trabalho foi realizada a implementação prática de um ciclo RTO em um processo de destilação por recompressão de vapor (VRD), propileno-propano, da Refinaria de Paulínia (Petrobras S.A.), a partir de dados históricos da planta. Foram consideradas as principais etapas de um ciclo clássico de RTO: identificação de estado estacionário, reconciliação de dados, estimação de parâmetros e otimização econômica. Essa unidade foi modelada, simulada e otimizada em EMSO (Environment for Modeling, Simulation and Optimization), um simulador de processos orientado a equações desenvolvido no Brasil. Foram analisados e comparados dois métodos de identificação de estado estacionário, um baseado no teste estatístico F e outro baseado em wavelets. Ambos os métodos tiveram resultados semelhantes e mostraram-se capazes de identificar os estados estacionários de forma satisfatória, embora seja necessário o ajuste de parâmetros na sua implementação. Foram identificados alguns pontos estacionários para serem submetidos ao ciclo RTO e foi possível verificar a importância de partir de um estado estacionário para a continuidade do ciclo, já que essa é uma premissa do método. A partir dos pontos analisados, os resultados deste estudo mostram que o RTO é capaz de aumentar o ganho econômico entre 2,5-24%, dependendo das condições iniciais consideradas, o que pode representar ganhos de até 18 milhões de dólares por ano. Além disso, para essa unidade, verificou-se que o compressor é um equipamento limitante no aumento de ganho econômico do processo. / In order to increase the profits of chemical plants, the Real-Time Optimization (RTO) is a tool that seeks to determine the steady-state optimal process operating conditions to maximize its profit under the operational restrictions. In this work, a practical implementation of a RTO cycle was implemented in a vapor recompression distillation (VRD) process, propylene-propane, from Paulínia Refinery (Petrobras S.A.), from historical plant data. The main steps of a classical RTO cycle are considered: steady-state identification, data reconciliation, parameter estimation and economical optimization. This unit was modeled, simulated and optimized in EMSO (Environment for Modeling, Simulation and Optimization), which is an equation oriented simulator conceived and developed in Brazil. Two steady state identification methods were analyzed and compared, one based on a F-like test and other based on wavelets. Both methods had similar results and showed to be able to identify the stationary states satisfactorily, although parameter tuning is necessary in their implementation. Some stationary points were identified and submitted to the RTO cycle. It was possible to verify the importance of start from a steady-state to continue the cycle, since this is a premise of the method. From the points analyzed, the results of this study show that the RTO is able to increase the economic gain between 2.5- 24%, depending on the initial conditions that is considered, which may represent gains of up to 18 million dollar per year. Furthermore, for this unit, it was found that the compressor is a limiting equipment in increasing economical gain.

Nesta tese consideramos um modelo de regressão não paramétrica, quando a variável explicativa e um processo estritamente estacionário e alpha-mixing. São estudadas as condições sobre o processo Xt e sua estrutura de dependência, assim como do domínio da função f a ser estimada. Também são feitas as adaptações necessárias aos procedimentos para obter as taxas de convergência do risco para a norma Lp, no caso de ondaletas deformadas. Em relação às ondaletas adaptativas de Haar, obtêm-se as taxas de convergência do risco do estimador proposto. Mediante estudos de simulação, e avaliado o desempenho dos procedimentos propostos quando aplicados a amostras finitas sob diferentes níveis de perturbação do sinal e diferentes tamanhos da amostra. Também são feitas aplicações a dados reais. / In this thesis we consider a nonparametric regression model, when the exploratory variables are alpha-mixing stationary processes. We obtain convergence rates for risk for Lp norm, via warped wavelets, under suitable regularity conditions. For estimation using design adapted Haar wavelets we obtain convergence rates for the risk of the proposed estimator. The performance of the estimators are assessed via simulation studies with dierent sample sizes and dierent signal-to-noise ratios. Applications to real data are also given.

Reatores de leito fluidizado circulante são caracterizados por escoamentos altamente heterogêneos e instáveis. Devido a isso, operam em regimes pseudo-estacionários ou estatisticamente permanentes, onde os parâmetros do escoamento oscilam em torno de valores médios bem definidos. Neste trabalho estuda-se o comportamento de uma solução numérica para um escoamento típico de reatores de leito fluidizado circulante na região estatisticamente permanente. A simulação numérica foi realizada com o modelo de dois fluidos. A partir dos dados de simulação são estabelecidas médias e auto-correlações para as várias variáveis de interesse, considerando-se diferentes tempos de simulação ou diferentes tempos de observação. A partir destes resultados discute-se a questão do tempo mínimo de simulação necessário para que os resultados temporais médios sejam representativos do processo. Embora os resultados sugiram que 10 segundos de fluidização sejam suficientes para a geração de resultados médios temporais representativos, uma demonstração mais rigorosa permanece por ser executada. / Circulating fluidized bed reactors are characterized by highly heterogenic and instable flow. Due to this flow, the reactors operate under either pseudo-stationary or statistically permanent patterns, where the flowing parameters vary among well defined average values. This paper aims to study the behavior of a numeric solution for a typical flow of circulating fluidized bed reactors in the statistically permanent region. The numeric simulation was performed with the two-fluid model. From the simulation data both average and self-correlations are established for several variables of interest, considering either different simulation time or different observation time. From these results, the minimum simulation time necessary is discussed, so that the average temporal results can be considered representative in the process. Although the results suggest 10 seconds of fluidization are enough to generate the average temporal results, a more rigorous demonstration is supposed to be executed.

Com o objetivo de aumentar o lucro de plantas químicas, a Otimização em Tempo Real (RTO) é uma ferramenta que busca determinar as condições ótimas operacionais do processo em estado estacionário, respeitando as restrições operacionais estabelecidas. Neste trabalho foi realizada a implementação prática de um ciclo RTO em um processo de destilação por recompressão de vapor (VRD), propileno-propano, da Refinaria de Paulínia (Petrobras S.A.), a partir de dados históricos da planta. Foram consideradas as principais etapas de um ciclo clássico de RTO: identificação de estado estacionário, reconciliação de dados, estimação de parâmetros e otimização econômica. Essa unidade foi modelada, simulada e otimizada em EMSO (Environment for Modeling, Simulation and Optimization), um simulador de processos orientado a equações desenvolvido no Brasil. Foram analisados e comparados dois métodos de identificação de estado estacionário, um baseado no teste estatístico F e outro baseado em wavelets. Ambos os métodos tiveram resultados semelhantes e mostraram-se capazes de identificar os estados estacionários de forma satisfatória, embora seja necessário o ajuste de parâmetros na sua implementação. Foram identificados alguns pontos estacionários para serem submetidos ao ciclo RTO e foi possível verificar a importância de partir de um estado estacionário para a continuidade do ciclo, já que essa é uma premissa do método. A partir dos pontos analisados, os resultados deste estudo mostram que o RTO é capaz de aumentar o ganho econômico entre 2,5-24%, dependendo das condições iniciais consideradas, o que pode representar ganhos de até 18 milhões de dólares por ano. Além disso, para essa unidade, verificou-se que o compressor é um equipamento limitante no aumento de ganho econômico do processo. / In order to increase the profits of chemical plants, the Real-Time Optimization (RTO) is a tool that seeks to determine the steady-state optimal process operating conditions to maximize its profit under the operational restrictions. In this work, a practical implementation of a RTO cycle was implemented in a vapor recompression distillation (VRD) process, propylene-propane, from Paulínia Refinery (Petrobras S.A.), from historical plant data. The main steps of a classical RTO cycle are considered: steady-state identification, data reconciliation, parameter estimation and economical optimization. This unit was modeled, simulated and optimized in EMSO (Environment for Modeling, Simulation and Optimization), which is an equation oriented simulator conceived and developed in Brazil. Two steady state identification methods were analyzed and compared, one based on a F-like test and other based on wavelets. Both methods had similar results and showed to be able to identify the stationary states satisfactorily, although parameter tuning is necessary in their implementation. Some stationary points were identified and submitted to the RTO cycle. It was possible to verify the importance of start from a steady-state to continue the cycle, since this is a premise of the method. From the points analyzed, the results of this study show that the RTO is able to increase the economic gain between 2.5- 24%, depending on the initial conditions that is considered, which may represent gains of up to 18 million dollar per year. Furthermore, for this unit, it was found that the compressor is a limiting equipment in increasing economical gain.

As estruturas aeronáuticas estão sujeitas a diversas solicitações, devido principalmente às interações com o escoamento aerodinâmico, que podem causar distúrbios e vibrações, comprometendo seu desempenho. As medidas aerodinâmicas aplicadas em uma aeronave podem ser obtidas por simulações computacionais ou testes experimentais. No entanto, podem existir imperfeições na simulação computacional, como por exemplo, se conseguir reproduzir algumas condições de vôo real. Sendo assim, diversas pesquisas vêm sendo realizadas para solucionar estes problemas. Dentre elas estão os testes experimentais feitos em túnel de vento com modelos de escala real em diversas condições de vôo. Desta forma, a construção de um modelo físico de um aerofólio em escala reduzida e a implementação de sensores a este modelo torna-se uma ferramenta bastante importante para validar resultados teóricos e experimentais. Assim, nesse trabalho realizou-se a construção de um modelo de aerofólio NACA0012, o desenvolvimento de um mecanismo de fixação do modelo ao túnel de vento e a implementação de um controlador de oscilação forçada. O modelo físico realiza oscilações harmônicas, em regime não-estacionário. O objetivo do trabalho foi mapear as medidas de pressão atuantes sobre modelo ensaiado em regime estacionário e não-estacionário e fazer a comparação entre os dois casos. / Aeronautical structures are affected by many loads, most of them given by the aerodynamic flow interactions. These flow interactions may cause vibration leading to structural failure, such as cracks and fatigue. The aerodynamic flow interactions can be measured by experiment or predicted by computational simulation. Otherwise, computational simulations on its own are not reliable and can not reproduce a real flight condition, such as the mean atmospheric turbulence dynamic. Many researches has been done to solve these problems for computational simulations. One of them are the wind tunnel experiments with a full scale models in many flight conditions for posterior comparison. For a smaller wind tunnel, a small scale physical prototype well instrumented becomes an important solution to validate theoretical and experimental results. In the present work the construction of a NACA 0012 airfoil model, the development of a constraint mechanism and the implementation of a forced oscilation control system were done. The physical model oscilates with a given frequency. The aim of present work is to map the pressure measurements actuating on the model, testing it under a steady state condition and a transient condition for posterior comparison of both conditions.

Reatores de leito fluidizado circulante são caracterizados por escoamentos altamente heterogêneos e instáveis. Devido a isso, operam em regimes pseudo-estacionários ou estatisticamente permanentes, onde os parâmetros do escoamento oscilam em torno de valores médios bem definidos. Neste trabalho estuda-se o comportamento de uma solução numérica para um escoamento típico de reatores de leito fluidizado circulante na região estatisticamente permanente. A simulação numérica foi realizada com o modelo de dois fluidos. A partir dos dados de simulação são estabelecidas médias e auto-correlações para as várias variáveis de interesse, considerando-se diferentes tempos de simulação ou diferentes tempos de observação. A partir destes resultados discute-se a questão do tempo mínimo de simulação necessário para que os resultados temporais médios sejam representativos do processo. Embora os resultados sugiram que 10 segundos de fluidização sejam suficientes para a geração de resultados médios temporais representativos, uma demonstração mais rigorosa permanece por ser executada. / Circulating fluidized bed reactors are characterized by highly heterogenic and instable flow. Due to this flow, the reactors operate under either pseudo-stationary or statistically permanent patterns, where the flowing parameters vary among well defined average values. This paper aims to study the behavior of a numeric solution for a typical flow of circulating fluidized bed reactors in the statistically permanent region. The numeric simulation was performed with the two-fluid model. From the simulation data both average and self-correlations are established for several variables of interest, considering either different simulation time or different observation time. From these results, the minimum simulation time necessary is discussed, so that the average temporal results can be considered representative in the process. Although the results suggest 10 seconds of fluidization are enough to generate the average temporal results, a more rigorous demonstration is supposed to be executed.

Neste trabalho foi desenvolvido um estudo sobre a precessão livre de onda contínua (CWFP-continuous wave free precession), uma condição especial da precessão livre no estado estacionário (SSFP), e algumas aplicações. Os resultados experimentais, e de simulações, mostraram que este sinal é atingido após a passagem por dois estágios iniciais. O primeiro é dependente da não homogeneidade do campo magnético, relacionado a razão Tp/T2*. O segundo, o qual é chamado de estado quasi-estacionáio, é dependente ambos os tempos de relaxação do sistema em estudo, T1 e T2. Este segundo estágio é responsável por levar a magnetização do estado de equilíbrio térmico a um estado verdadeiramente estacionário. Tendo como informação a amplitude da magnetização no estado de equilíbrio térmico e no estado verdadeiramente estacionário, e também do decaimento do sinal durante o estado quasi-estacionário, é possível determinar, rápida e simultaneamente, ambos os tempos de relaxação do sistema, fazendo um único experimento. Esse método oferece aplicações interessantes para o estudo de processos dinâmicos, propostas também neste trabalho. Foi possível verificar que o sinal de CWFP pode ser utilizado para acompanhar cinética de reações, e também variações de viscosidade do meio, por exemplo, em reações de polimerização, associadas à mobilidade molecular em função de mudanças estruturais. Um método para avaliar a condutividade térmica de elastômeros foi proposto, com o fornecimento de resultados quantitativos muito próximos daqueles encontrados na literatura, e que utilizam outros métodos. A dependência do sinal de CWFP com a freqüência de offset também permitiu realizar um estudo, com aplicação em RMN de alta resolução, sobre a supressão de sinal de solvente. Os resultados demonstraram que é necessário fazer alguns ajustes na largura de pulsos, juntamente com o ciclo de fases, para minimizar as anomalias de intensidade e de fase dos sinais no espectro transformado. Essa técnica foi aplicada em espectroscopia ‘in vivo’, a qual permite resolver, por exemplo, problemas como a determinação de açúcares em frutas, onde o sinal é bastante próximo do intenso sinal da água. / In this work a study of the continuous wave free precession (CWFP), a special condition of the steady-state free precession (SSFP), and some applications, was developed. The experimental results, together with simulated ones, have shown that the CWFP signal is attained after two previous stages. The first one is dependent on the non-homogeneity of the magnetic field, related to Tp/T2* ratio. The second, which is called quasi-stationary state, is dependent of both relaxation times, T1 and T2. This second stage allows leading the magnetization in the thermal equilibrium to a truly stationary state. The information of the signal amplitude in the thermal equilibrium and in the steady state, and also the signal decaying during the quasistationary state, allows the fast and simultaneous determination of the relaxation times, performed in a single experiment. This method offers interesting applications for studying dynamical processes, also proposed in this work. It was possible to verify that the CWFP signal can be used to monitor kinetics of reactions, which variations of viscosity are involved, for example, in polymerization reactions, always associated to structural changes. A method for evaluating thermal conductivity in elastomers was proposed, with results in accordance to the literature, which uses other methods. The dependence of the CWFP signal with the offset frequency has permitted to study an application to the High-Resolution NMR, about the solvent suppression. The results demonstrated that it is necessary to perform some adjustments in the pulse width, together the phase cycle, to minimize some phase and intensity anomalies of the frequency-domain signal. This technique was applied in ‘in vivo’ spectroscopy, which can solve, for example, problems with the determination of the sugar content in fresh fruits, where the sugar signal is very close to the intense signal of the water.

Este trabalho consiste em realizar um estudo detalhado das vantagens e desvantagens da utilização do FDM (Filter Diagonalization Method) para a análise de dados obtidos pela sequência de Precessão Livre no Estado Estacionário (Steady State Free Precession - SSFP) para aquisição rápida de espectros de Ressonância Magnética Nuclear (RMN). No caso de RMN de baixa resolução, o procedimento de aquisição rápida, SSFP, é uma poderosa ferramenta para melhorar a relação sinal/ruído, apresentando muitas aplicações práticas. Apesar desse sucesso em baixa resolução, a SSFP não é rotineiramente utilizada para aplicações em RMN de alta resolução, provavelmente devido ao (1) artefatos provenientes do truncamento do sinal e (2) as anomalias causadas pela mistura do FID com o eco dos sinais. Existem na literatura inúmeras possíveis técnicas para suprimir este tipo de problemas, porém, nenhuma delas é capaz de realmente eliminar as anomalias geradas devido ao procedimento de aquisição rápida da SSFP. O FDM é um método paramétrico não-linear para fitar sinais no domínio do tempo. Seu objetivo fundamental é resolver o Problema da Inversão Harmônica, HIP, tornando-se robusto e adequado para a análise espectral de sinais no domínio do tempo nos casos onde a Transformada de Fourier falha. Neste trabalho, demonstramos que o FDM pode ser implementado para análises de sinais SSFP, com mais eficiência que os obtidos pelos procedimentos padrões de TF. A temperatura ambiente, espectros de RMN 13C de amostras de brucina, obtidos com tempo entre pulsos de 100ms, podem ser reproduzidos com boa relação sinal/ruído e alta resolução por meio do FDM. A limitação da análise por FDM é mais relevante nos casos de espectros com alta densidade de picos em uma determinada região espectral. Nestes casos, o curto período de observação do sinal na janela do tempo impõe uma série de limitações na resolução obtida pelo FDM. / This work consists in a detailed study of the advantages and disadvantages of the use of the Filter Diagonalization Method, FDM, for data analysis in Steady State Free Precession, SSFP, technique, usually employed to implement fast acquisition of Nuclear Magnetic Resonance, NMR, spectra. In the case of low resolution NMR using fast acquisition procedures, SSFP is a powerful tool to improve signal-to-noise ratio, presenting several important practical applications. Despite its success in the low resolution regime, SSFP is not a routine technique for high resolution applications, so far, mainly because of (1) truncation artifacts and (2) the intrinsic anomalies caused by admixture of free-induction-decay and echo signals. The literature reports many possible techniques to solve such kind of problems, but, none of them is capable to really eliminate the generated spectra anomalies caused by the fast acquisition procedure used in SSFP. FDM is a parametric method for non-liner fitting performed in the time domain. Its main goal is to solve the Harmonic Inversion Problem, HIP, making it robust and suitable for spectral analysis of time signals in the cases where the Fourier Transform, FT, technique fail. In this work we demonstrate that FDM can be used to implement the analysis of the SSFP data, with more efficiency than that achieve by appropriated FT procedures. Room temperature 13C NMR spectra of brucine samples, obtained from pulse sequences with 100 ms repetition time, can be reproduced with good signal-to-noise ratio and high resolution by means of the FDM. The limitation of the FDM analysis is more relevant in the case of spectra with a high density of peaks in a limited spectral frequency region. In these cases, the reduced short observation time window imposes serious limitation to the resolution achieved by the FDM.

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas. / The need for solving large linear systems arising from the discretization of partial differential equations modelling physical phenomena motivates the search for scalable numerical techniques. Multigrid algorithms are instances of such techniques.In order to provide a suitable assessment of the solution obtained by such algorithms, an error estimator must be associated to the numerical solution of the discretized problem. In this context, this thesis proposes the reutilization of the hierarchical matrix structures of transfer operators and the restriction to algebraic multigrid methods to speed up the process of solving the linear systems associated with the contaminant transport equation in saturated porous media. In addition, it features the implementation of residual estimates for problems involving constant or non-constant data, the regimes of small- or large-scale advection and the proposal of employing the residual estimates associated to the source term and to the initial condition to build adaptive procedures for the problem data. The development of the computer codes of the finite element method, residual estimator and adaptive procedures were based on the FEniCS project, using the programming language PYTHONR and developed on the Eclipse platform. The implementation of the algebraic methods with reutilization relied upon the libray PyAMG. Grounding on the idea of reutilizing the hierarchical structures, fixed and automatic parameters multigrid methods were proposed and extended to non-stationary iterative methods such as GMRES and BICGSTAB. The numerical results demonstrate that the residual estimator captures the behavior of the real error of the numerical solution, and provide adaptive algorithms for the data whose output mesh yields a numerical solution alike to that obtained from a uniform mesh with more elements. Moreover, the methods with reutilization are faster than those that do not reuse the structures. Besides, the efficiency of such methods can also be observed in the solution of an auxiliary problem, which is necessary for deriving the residual estimates in the regime of large-scale advection. These results encompass both the type SA algebraic multigrid method and those pre-conditioned by them. Moreover, they involve the transport of contaminants in regime of small- and large-scale advection, structured and non-structured meshes, bi- and tridimensional problems and domains with different scales.

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas. / The need for solving large linear systems arising from the discretization of partial differential equations modelling physical phenomena motivates the search for scalable numerical techniques. Multigrid algorithms are instances of such techniques.In order to provide a suitable assessment of the solution obtained by such algorithms, an error estimator must be associated to the numerical solution of the discretized problem. In this context, this thesis proposes the reutilization of the hierarchical matrix structures of transfer operators and the restriction to algebraic multigrid methods to speed up the process of solving the linear systems associated with the contaminant transport equation in saturated porous media. In addition, it features the implementation of residual estimates for problems involving constant or non-constant data, the regimes of small- or large-scale advection and the proposal of employing the residual estimates associated to the source term and to the initial condition to build adaptive procedures for the problem data. The development of the computer codes of the finite element method, residual estimator and adaptive procedures were based on the FEniCS project, using the programming language PYTHONR and developed on the Eclipse platform. The implementation of the algebraic methods with reutilization relied upon the libray PyAMG. Grounding on the idea of reutilizing the hierarchical structures, fixed and automatic parameters multigrid methods were proposed and extended to non-stationary iterative methods such as GMRES and BICGSTAB. The numerical results demonstrate that the residual estimator captures the behavior of the real error of the numerical solution, and provide adaptive algorithms for the data whose output mesh yields a numerical solution alike to that obtained from a uniform mesh with more elements. Moreover, the methods with reutilization are faster than those that do not reuse the structures. Besides, the efficiency of such methods can also be observed in the solution of an auxiliary problem, which is necessary for deriving the residual estimates in the regime of large-scale advection. These results encompass both the type SA algebraic multigrid method and those pre-conditioned by them. Moreover, they involve the transport of contaminants in regime of small- and large-scale advection, structured and non-structured meshes, bi- and tridimensional problems and domains with different scales.