Diversos procesos naturales e industriales de interés medioambiental se modelan a través de la ecuación de convección-difusión-reacción transitoria. Dos aplicaciones tecnológicas que han motivado esta tesis son el funcionamiento de filtros de carbón activo y la dispersión de contaminantes en la atmósfera. Para que la modelización numérica de estos problemas sea eficaz es indispensable contar con un solver lineal eficiente para resolver los sistemas de ecuaciones obtenidos al discretizar la ecuación en derivadas parciales, mediante elementos finitos.Por ello, el objetivo de esta tesis es resolver de forma eficiente los grandes sistemas de ecuaciones, simétricos definidos positivos (SDP), tipo sparse asociados a los problemas de convección-difusión transitorios. Con este fin se estudian los precondicionadores tanto explícitos como implícitos, así como los métodos de descomposición de dominios (DD). La tesis se estructura en tres partes. En la primera se elabora un análisis computacional detallado del comportamiento de dos familias de factorizaciones incompletas de Cholesky (FIC): de memoria prescrita y de umbral. Estas técnicas se utilizan para precondicionar el método iterativo de gradientes conjugados (PCG). En la segunda parte se construye una inversa aproximada sparse simétrica (SSPAI) basada en la minimización en la norma de Frobenius. El precondicionador explícito se diseña para resolver en paralelo grandes sistemas de ecuaciones sparse, SDP, tridiagonales por bloques con múltiples lados derechos. Finalmente, se desarrolla el método multiplicativo de Schwarz (MSM) en dominios activos, es decir, DD solapados con la innovación de activar y desactivar dominios. Se estudia el comportamiento de esta estrategia al resolver los subproblemas mediante: (1) el método directo de Cholesky y (2) PCG + FIC de umbral.De los resultados numéricos presentados se concluye que es preferible utilizar el método directo de Cholesky para sistemas con menos de 30,000 variables. Para sistemas mayores y hasta 80,000 incógnitas se sugiere emplear una FIC de umbral. Y para sistemas aún más grandes, el MSM en dominios activos + PCG + FIC de umbral propuesto es el más eficiente usando un solo procesador. Por su parte, la SSPAI presentada podría superar a las FIC de umbral si se trabaja en paralelo. / Many natural and industrial processes of environmental interest are modeled through the transient convection-diffusion-reaction equation. Two technological applications that have motivated this thesis are the operation of activated-carbon filters and the dispersion of pollutants in the atmosphere. In order to ensure the effectiveness of numerical modeling of these problems it is necessary to have an efficient linear solver to solve the systems obtained when discretizing the partial differential equation by means of the finite element method.For that reason, the goal of this thesis is to solve in an efficient way the large sparse symmetric positive definite (SPD) systems of linear equations associated to the transient convection-diffusion problems. With this purpose, we have studied the explicit and implicit preconditioners, as well as domain decomposition (DD) methods.The thesis is structured in three parts. In the first one we have elaborated a detailed analysis of the numerical performance of two families of incomplete Cholesky factorizations (ICF): drop tolerance and prescribed-memory strategies. These techniques are used to precondition conjugate gradient iterations (PCG). In the second part a symmetric sparse approximate inverse (SSPAI) based on the minimization of the Frobenius norm is built. The explicit preconditioner is designed to solve in parallel large block tridiagonal SPD systems with multiple right hand sides.Finally, the multiplicative Schwarz method (MSM) in active domains is developed, which consists of overlapped domain decomposition with the innovation to activate and to deactivate domains. The behavior of this strategy is studied when solving the subproblems by means of: (1) the Cholesky direct solver and (2) PCG + drop-tolerance ICF. According to our numerical experiments we conclude that it is preferable to use the Cholesky direct solver for systems with less than 30,000 variables. For larger systems and up to 80,000 equations we suggest to use a drop tolerance ICF. And for even lager systems, the proposed MSM in active domains is the most efficient when using a single processor. On the other hand, the presented SSPAI could overcome the drop-tolerance ICF if one works in parallel.

[pt] Companhias de petróleo investindo no desenvolvimento de campos de hidrocarboneto dependem de estudos de reservatórios para realizarem previsões de produção e quantificarem os riscos associados à economicidade dos projetos. Neste sentido, a área de modelagem de reservatórios é de suma importância, sendo responsável por prever o desempenho futuro do reservatório sob diversas condições operacionais. Considerando que a solução dos sistemas de equações construídos a cada passo de tempo de uma simulação, durante o ciclo de linearização, é a parte que apresenta a maior demanda computacional, esta dissertação foca na análise de diferentes técnicas de solvers numéricos que podem ser aplicadas a simuladores, para mensurar seus desempenhos. Os solvers numéricos mais adequados para a solução de grandes sistemas de equações, tais como os encontrados em simulações de reservatórios, são os denominados solvers iterativos, que gradativamente aproximam a solução de um dado problema por meio da combinação de um método iterativo e um precondicionador. Os métodos iterativos avaliados nesta pesquisa foram o Gradiente Biconjugado Estabilizado (BiCGSTAB), Mínimos Resíduos Generalizado (GMRES) e Minimização Ortogonal (ORTHOMIN). Além disso, três técnicas de precondicionamento foram implementadas para auxiliar os métodos iterativos, sendo estas a Decomposição LU Incompleta (ILU), Fatoração Aninhada (NF) e Pressão Residual Restrita (CPR). A combinação destes diferentes métodos iterativos e precondicionadores permite a avaliação de diversas configurações distintas de solvers, em termos de seus desempenhos em um simulador. Os testes numéricos conduzidos neste trabalho utilizaram um novo simulador de reservatórios que está sendo desenvolvido pela Pontifícia Universidade Católica (PUC-Rio) em conjunto com a Petrobras. O objetivo dos testes foi analisar a robustez e eficiência de cada um dos solvers quanto à sua capacidade de resolver as equações de escoamento multifásico no meio poroso, visando assim auxiliar na seleção do solver mais adequado para o simulador. / [en] Petroleum companies investing in the development of hydrocarbon fields rely upon a variety of reservoir studies to perform production forecasts and quantify the risks associated with the economics of their projects. Integral to these studies is the discipline of reservoir modeling, responsible for predicting future reservoir performance under various operational conditions. Considering that the most time-demanding aspect of reservoir simulations is the solution of the systems of equations that arise within the linearization cycles at each time-step, this research focuses on analyzing different numerical solver techniques to be applied to a simulator, in order to assess their performance. The numerical solvers most suited for the solution of very large systems of equations, such as those encountered in reservoir simulations, are the so-called iterative solvers, which gradually approach the solution to a problem by combining an iterative strategy with a preconditioning method. The iterative methods examined in this research were the Stabilized Biconjugate Gradient (BiCGSTAB), the Generalized Minimum Residual (GMRES), and the Orthogonal Minimization (ORTHOMIN) methods. Furthermore, three preconditioning techniques were implemented to aid the iterative methods, namely the Incomplete LU Factorization (ILU), the Nested Factorization (NF), and the Constrained Pressure Residual (CPR) methods. The combination of these different iterative methods and preconditioners enables the appraisal of several distinct solver configurations, in terms of their performance in a simulator. The numerical tests conducted in this work made use of a new reservoir simulator currently under development at Pontifical Catholic University of Rio de Janeiro (PUC-Rio), as part of a joint project with Petrobras. The objective of these tests was to assess the robustness and efficiency of each solver in the solution of the multiphase flow equations in porous media, and support the selection of the solver most suited for the simulator.