O interesse na medição do perfil de velocidade na área de fluidodinâmica tem crescido nos últimos anos devido a evolução das técnicas de medição. Nesse contexto, a técnica ultrassônica tem se destacado por ser não intrusiva, não invasiva e funcionar mesmo em fluidos opacos. Neste trabalho foi investigado a técnica de ultrassom Doppler para medição de perfis de velocidade, utilizando equipamentos ultrassônicos para ensaios-não-destrutivos (END). Tais equipamentos são mais acessíveis do que equipamentos convencionais de medição de velocidade por ultrassom. Também foi proposto o uso da técnica de autocorrelação estendida (EAM), para medição de velocidades além do limite de Nyquist. Essa restrição existe na grande maioria dos medidores, que utilizam a técnica convencional de autocorrelação (ACM). O EAM combina o ACM com o método da correlação cruzada (CCM), outro estimador amplamente conhecido, mas que não é muito utilizado devido a seu alto custo computacional. Desta forma, o EAM consegue medir velocidades maiores com um custo computacional intermediário, que não é tão baixo quanto o ACM e nem tão alto quanto o CCM. Para adquirir e processar os dados obtidos, foi desenvolvido um sistema para aquisição e processamento dos sinais baseado na linguagem LabView. O pulsador ultrassônico END e os estimadores de velocidade por EAM, CCM e ACM foram validados medindo o perfil de velocidade em um cilindro girante, capaz de fornecer velocidades controladas de fácil solução analítica. Os resultados mostram erros médios quadráticos abaixo de 2%, validando o equipamento e a técnica. O EAM também mediu velocidades acima no limite de Nyquist com um desempenho computacional de 9 vezes maior do que o CCM. Na segunda parte deste trabalho, a técnica ultrassônica de medição de velocidade foi aplicada para medição de escoamentos multifásicos em tubulações, tema de grande interesse da indústria de petróleo e gás. Um escoamento vertical líquido-gás-sólido foi analisado com o mesmo equipamento END. Primeiramente, o perfil de velocidade do escoamento liquido-sólido em regime laminar, foi medido e validado utilizando um equipamento Coriolis como referência. Em seguida, foram realizados testes adicionando gás ao escoamento. As velocidades superficiais de líquido e gás foram variadas para gerar os padrões de escoamento tipo bolhas dispersas, intermitente e intermitente aerado. Os resultados foram comparados a imagens de uma filmadora de alta velocidade. Foram obtidos parâmetros como perfil de velocidade das bolhas dispersas, velocidade do filme de líquido e velocidade da mistura dependendo do padrão de escoamento analisado. Assim, a medição de velocidade por ultrassom Doppler, utilizando um equipamento de END, foi aplicado com sucesso em dois problemas de fluidodinâmica. / Interest in knowing the instantaneous velocity profile in fluid dynamics has grown in recent years as new flow visualization techniques are improving. In this context, the ultrasonic Doppler velocity profiler (UVP) has desirable characteristics, as it is non-invasive, works with opaque liquids, and it is portable and easy to install if compared with other velocity profiler methods. In this work, the use of nondestructive ultrasonic devices in the UVP field is investigated. NDT systems are widely available and have lower cost if compared to traditional ultrasonic velocity profiler systems. The use of an extended autocorrelation method (EAM) for ultrasonic velocity estimation beyond Nyquist limit are also evaluated. The Nyquist limit causes a restriction on the maximum measurable velocity of the traditional autocorrelation method (ACM), present in most of ultrasonic velocity profiler systems. EAM combines the ACM technique with cross-correlation method (CCM) which is a well-established velocity estimator that does not suffer with Nyquist limit. However, the technique has a high computational cost that limits real time applications. EAM has the advantage of measure velocities beyond the Nyquist limit but with a lower computational cost than CCM. To evaluate the NDT device and the velocity estimation techniques ACM, CCM and EAM, a data acquisition system and a signal-processing unit based on LabView language were developed. The velocity profile of a rotating cylinder was used to validate all measurements. The techniques ACM, CCM and EAM successfully measured velocities within Nyquist limit with less than 2% deviation, validating the NDT system. EAM also measured velocities beyond Nyquist limit with a computational performance 9 times faster than CCM. The ultrasonic technique was also applied to measure the velocity profile of a multiphase flow in a pipeline, which are of great interest in oil and gas industry. Tests within a multiphase flow composed by different combinations of oil/sand/nylon-particles and gas were conducted in a vertical rig. A high-speed camera was used to validate the measurements. First the ultrasonic velocity profile measured was validated in a liquid-solid flow by using a Coriolis flowmeter as a reference. Next, superficial liquid and gas velocity were controlled to obtain different flow patterns such as bubbly flow and slug flow. The technique measured the bubbles velocity, the mixture velocity and the liquid film velocity depending on the flow pattern. Finally, the ultrasonic NDT system was successfully applied to investigate two different fluid engineering problems.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.

Técnicas de injeção de traçadores são bastante utilizadas nos estudos de escoamentos em meios porosos heterogêneos, principalmente em problemas relacionados à simulação numérica de escoamentos miscíveis em reservatórios de petróleo e à dispersão de contaminantes em aqüíferos. Neste trabalho apresentamos novos algoritmos para a aproximação numérica do problema de injeção de traçadores. Apresentaremos desenvolvimentos recentes do método Forward Integral-Tube Tracking (FIT) que foi originalmente apresentado em Aquino et al. (2007a). O FIT é um método lagrangeano localmente conservativo utilizado na resolução de problemas de transporte linear. Este método não faz o uso de soluções de problemas de Riemann e baseia-se na construção dos tubos integrais introduzidas em Douglas Jr. et al. (2000b). Além disso, ele possui excelente eficiência computacional e é virtualmente livre de difusão numérica. Resultados numéricos são apresentados com o objetivo de comparar a precisão das soluções fornecidas por novas implementações do método FIT na resolução do problema do traçador em reservatórios de petróleo. / The injection of tracers are used in the investigation of flows in heterogeneous porous media, in studies related to the simulation of miscible dispacements in petroleum reservoirs and the dispersion of contaminants in aquifers. In this work we present new algorithms for the numerical approximation of tracer injection problems. We discuss recent developments of the Forward Integral-Tube Tracking (FIT) scheme which was introduced in Aquino et al. (2007a). The FIT is a locally conservative lagrangian scheme for the approximation of the linear transport problems. This scheme does not use analytic solutions of Riemann problems and is based on the construction of the integral tubes introduced in Douglas Jr. et al. (2000b). The FIT scheme is computationally very eficient and is virtually free of numerical diffusion. Numerical results are presented to compare the accuracy of the solutions provided by new implementation of the FIT scheme for the injection of tracers in petroleum reservoirs.

Petróleo Brasileiro S.A. / O presente trabalho aborda o comportamento da planta de processamento primário com alimentação por fluxo em padrão de golfadas. O fluxo no sistema de tubulações é descrito por um modelo de parâmetros concentrados, fornecendo as características principais necessárias para o controle da planta, e a resposta dinâmica desta pode então ser analisada. Usando a estratégia de controle tradicional verifica-se que as oscilações de fluxo são transmitidas para as vazões de saída de líquido e gás, para obter uma vazão de saída mais estável é permitida a flutuação de carga no separador dentro de uma tolerância, isto é possível reduzindo a atuação do controlador e estabelecendo um controle adicional diretamente na válvula de entrada. / The present work addresses the behavior of a primary processing plant subjected to slug flow pattern at its entrance. The flow in a pipeline system is described by a simplified concentrated parameter model, which preserves the main physical features that are important to control the plant. The dynamic response of the plant is then analyzed. Using a standard control strategy for the gas liquid separator, it is seen that the flow oscillations are transmitted to the liquid and gas outlets. In order to obtain a more stable outlet flow, the liquid level in the separator is then allowed to fluctuate within a given range, by reducing the effect of the controller constants, and establishing an additional control directly on the inlet entrance valve.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho é apresentado um novo procedimento numérico para o upscaling de permeabilidade utilizando condições de contorno periódicas. Este procedimento combina decomposição de domínio com elementos finitos mistos na discretização do problema local de pressão-velocidade necessário para se encontrar as permeabilidades equivalentes. / A new numerical method is proposed for the permeabilities upscaling take into consideration periodic boundary conditions. This method combines domain decomposition with mixed finite elements in discretization of the local problem of pressure-velocity necessary to meet the equivalent permeabilities.

CAPES / Vários estudos numéricos investigaram cavidades fechadas sob o efeito da convecção natural preenchidas com fluidos newtonianos generalizados (FNG) nos últimos anos pelas aplicações diretas em trocadores de calor compactos, no resfriamento de sistemas eletrônicos e na engenharia de polímeros. Neste trabalho é realizada a investigação numérica do processo de convecção natural de fluidos de lei de Potência e de Bingham em cavidades fechadas, aquecidas lateralmente e preenchidas com meios heterogêneos e bloco centrado. O meio heterogêneo é constituído de blocos sólidos, quadrados, desconectados e condutores de calor. Como parâmetros são utilizados a faixa de Rayleigh de 104 à 107, índice de potência n de 0, 6 à 1, 6, número de Bingham de 0, 5 até Bimax , sendo investigado da influência do número de Prandtl para cada modelo de fluido. Nas cavidades com meio heterogêneo são utilizadas as quantidades de blocos de 9, 16, 36 e 64, mantendo-se a razão entre a condutividade térmica do sólido e do fluido κ = 1. Para as cavidades com bloco centrado, são utilizados os tamanhos adimensionais de 0, 1 à 0, 9 com κ = 0, 1; 1 e 10. A modelagem matemática é realizada pelas equações de balanço de massa, de quantidade de movimento e de energia. As simulações são conduzidas no programa comercial ANSYS FLUENT R . Inicialmente são resolvidos problemas com fluidos newtonianos em cavidade limpa, seguida de cavidade preenchida com meio heterogêneo e posteriormente bloco centrado para validação da metodologia de solução. Na segunda etapa é realizada o estudo com os modelos de fluidos de lei de Potência e de Bingham seguindo a mesma sequência. Os resultados são apresentados na forma de linhas de corrente, isotermas e pelo número de Nusselt médio na parede quente. De maneira geral, a transferência de calor na cavidade é regida pelo número de Rayleigh, tamanho e condutividade térmica dos blocos, pelo índice de potência para o modelo de lei de Potência e do número de Bingham para o modelo de Bingham. O número de Prandtl tem grande influência nos dois modelos de fluidos. O meio heterogêneo reduz a transferência de calor na cavidade quando interfere na camada limite térmica para ambos os fluidos, sendo feita uma previsão analítica para o fluido de lei de Potência. Para bloco centrado, a interferência na camada limite com fluido de lei de Potência também foi prevista analiticamente. A transferência de calor aumentou com bloco de baixa condutividade térmica e pouca interferência e com bloco de alta condutividade térmica e grande interferência, para ambos os fluidos. / Many studies have been carried out in square enclosures with generalized Newtonian fluids with natural convection in past few years for directly applications in compact heat exchangers, cooling of electronics systems and polymeric engineering. The natural convection in square enclosures with differently heated sidewalls, filled with power-law and Bingham fluids in addition with heterogeneous medium and centered block are analyzed in this study. The heterogeneous medium are solid, square, disconnected and conducting blocks. The parameters used are the Rayleigh number in the range 104 - 107 , power index n range of 0, 6 - 1, 6, Bingham number range of 0, 5 - Bimax , being the influence of Prandtl number investigated for each fluid model. The number of blocks for heterogeneous medium are 9, 16, 36 and 64, keeping constant solid to fluid conductive ratio, κ = 1. For enclosures with centered block are used the nondimensional block size from 0, 1 to 0, 9, with solid to fluid conductive ratio in range κ = 0, 1; 1 and 10. Mathematical modeling is done by mass, momentum and energy balance equations. The solution of equations have been numerically solved in ANSYS FLUENT R software. Firstly, numerical solutions for validation with Newtonian fluids in clean enclosures are conducted, followed by enclosures with heterogeneous medium and centered block. Subsequently, numerical solutions of power-law and Bingham fluids with same enclosures configurations are conducted. The results are reported in the form of streamlines, isotherms and average Nusselt number at hot wall. In general, the heat transfer process in enclosure is governed by Rayleigh number, size and thermal conductivity of the blocks, power index n for power-law fluid and Bingham number for Bingham fluid. Both fluid models are very sensitive with Prandtl number changes. Heterogeneous medium decrease heat transfer in enclosure when affects thermal boundary layer for both fluid models. One analytical prediction was made for power-law fluid. An increase in heat transfer occurs with low thermal conductivity block and few interference and with high thermal conductivity block and great interference, for both fluids.

Durante a perfuração de poços de petróleo é comum que a pressão do fluido de perfuração no poço seja maior que a pressão no interior da formação, o que pode fazer com que o fluido escoe no sentido da formação, num fenômeno denominado invasão. O escoamento no sentido da formação porosa carrega partículas que são retidas pelo substrato poroso, através do mecanismo de filtração, formando um leito de partículas na parede do poço. Neste trabalho, o modelo matemático e numérico do escoamento particulado através de um canal poroso posicionado verticalmente é investigado. O aumento da perda de carga provocado pela eventual obstrução do escoamento sobre o meio poroso devido à deposição de material particulado é investigado. O substrato poroso é representado através do modelo heterogêneo, cujo domínio sólido é descrito através de cilindros desconectados. A porosidade ao longo do domínio poroso é variada na direção axial do escoamento de acordo com a dimensão dos obstáculos, os quais são alocados num arranjo alternado. A análise é realizada em duas etapas. Na primeira o escoamento monofásico do fluido permite a determinação da permeabilidade e da perda de carga do meio poroso. A segunda etapa trata do processo de deposição de partículas presentes no escoamento sólido-fluido sobre o substrato poroso. A formulação matemática e a modelagem numérica para o escoamento particulado são representadas por uma abordagem Euler-Lagrange. A solução acoplada das fases discreta (partículas) e contínua (fluido) é realizada através da combinação dos modelos Dense Discrete Phase Model (DDPM) e Discrete Element Method (DEM). Resultados são obtidos para um fluido de massa específica e viscosidade semelhante à de um fluido de perfuração, com propriedades de uma mistura de 37,3% de água e 73,7% de glicerina, as partículas têm diâmetro variado de 0,6 a 0,8 mm e a porosidade do meio poroso varia entre 0,4 e 0,7. Resultados mostram que o aumento do diâmetro das partículas promove a redução da permeabilidade do meio resultante formado pelo meio poroso e o leito de partículas. O acréscimo da concentração de partículas injetadas implica em um aumento na queda de pressão, consequentemente na redução da permeabilidade através do canal. É possível observar também que existe um crescimento na espessura do leito com o aumento da concentração de partículas injetadas. / During oil well drilling commonly the drilling fluid pressure in the well is greater than the pressure within the formation, which may cause the flow towards the porous substrate, a phenomenon referred to as invasion. The flow to porous formation carries particles which are retained by the substrate, through the filtration mechanism, originating a packed-bed of particles in the wall of the wellbore. In this paper, the mathematical and numerical model of particle flow through a porous channel positioned vertically is proposed. Increase in pressure loss caused by any obstruction of the flow through porous media due to deposition of particulate material is investigated. The porous substrate is represented by a heterogeneous model, which solid domain is described by disconnected cylinders. Porosity throughout the porous domain varies in the axial direction of the flow according to obstacle sizes, which are allocated in a staggered array. The analysis is performed in two steps. In the first, single-phase fluid flow allows the determination of permeability and pressure drop of the porous media. The second step is the process of particles deposition observed in two-phase flow into the porous substrate. Mathematical formulation and numerical modeling for particulate flow are represented by a Euler-Lagrange approach. The coupled solution of discrete phases (particles) and continuous (fluid) is performed by combining the models Dense Discrete Phase Model (DDPM) and Discrete Element Method (DEM). Results are obtained for a fluid with density and viscosity similar to a drilling fluid, with properties of a mixture with 37.3% of water and 73.7% of glycerin, the particles have diameters varying from 0.6 to 0.8 mm and the porosity of the porous medium is between 0.4 and 0.7. Results show that the diameter of the particles increase promotes the permeability reduction of the resulting medium formed by the porous medium and the particles bed. The increase of the injected particles concentration implies an increase in pressure drop, thus reducing the permeability through the channel. It is also observed that there is growth in the thickness of the bed with the increase of the concentration of injected particles.

CAPES / Vários estudos numéricos investigaram cavidades fechadas sob o efeito da convecção natural preenchidas com fluidos newtonianos generalizados (FNG) nos últimos anos pelas aplicações diretas em trocadores de calor compactos, no resfriamento de sistemas eletrônicos e na engenharia de polímeros. Neste trabalho é realizada a investigação numérica do processo de convecção natural de fluidos de lei de Potência e de Bingham em cavidades fechadas, aquecidas lateralmente e preenchidas com meios heterogêneos e bloco centrado. O meio heterogêneo é constituído de blocos sólidos, quadrados, desconectados e condutores de calor. Como parâmetros são utilizados a faixa de Rayleigh de 104 à 107, índice de potência n de 0, 6 à 1, 6, número de Bingham de 0, 5 até Bimax , sendo investigado da influência do número de Prandtl para cada modelo de fluido. Nas cavidades com meio heterogêneo são utilizadas as quantidades de blocos de 9, 16, 36 e 64, mantendo-se a razão entre a condutividade térmica do sólido e do fluido κ = 1. Para as cavidades com bloco centrado, são utilizados os tamanhos adimensionais de 0, 1 à 0, 9 com κ = 0, 1; 1 e 10. A modelagem matemática é realizada pelas equações de balanço de massa, de quantidade de movimento e de energia. As simulações são conduzidas no programa comercial ANSYS FLUENT R . Inicialmente são resolvidos problemas com fluidos newtonianos em cavidade limpa, seguida de cavidade preenchida com meio heterogêneo e posteriormente bloco centrado para validação da metodologia de solução. Na segunda etapa é realizada o estudo com os modelos de fluidos de lei de Potência e de Bingham seguindo a mesma sequência. Os resultados são apresentados na forma de linhas de corrente, isotermas e pelo número de Nusselt médio na parede quente. De maneira geral, a transferência de calor na cavidade é regida pelo número de Rayleigh, tamanho e condutividade térmica dos blocos, pelo índice de potência para o modelo de lei de Potência e do número de Bingham para o modelo de Bingham. O número de Prandtl tem grande influência nos dois modelos de fluidos. O meio heterogêneo reduz a transferência de calor na cavidade quando interfere na camada limite térmica para ambos os fluidos, sendo feita uma previsão analítica para o fluido de lei de Potência. Para bloco centrado, a interferência na camada limite com fluido de lei de Potência também foi prevista analiticamente. A transferência de calor aumentou com bloco de baixa condutividade térmica e pouca interferência e com bloco de alta condutividade térmica e grande interferência, para ambos os fluidos. / Many studies have been carried out in square enclosures with generalized Newtonian fluids with natural convection in past few years for directly applications in compact heat exchangers, cooling of electronics systems and polymeric engineering. The natural convection in square enclosures with differently heated sidewalls, filled with power-law and Bingham fluids in addition with heterogeneous medium and centered block are analyzed in this study. The heterogeneous medium are solid, square, disconnected and conducting blocks. The parameters used are the Rayleigh number in the range 104 - 107 , power index n range of 0, 6 - 1, 6, Bingham number range of 0, 5 - Bimax , being the influence of Prandtl number investigated for each fluid model. The number of blocks for heterogeneous medium are 9, 16, 36 and 64, keeping constant solid to fluid conductive ratio, κ = 1. For enclosures with centered block are used the nondimensional block size from 0, 1 to 0, 9, with solid to fluid conductive ratio in range κ = 0, 1; 1 and 10. Mathematical modeling is done by mass, momentum and energy balance equations. The solution of equations have been numerically solved in ANSYS FLUENT R software. Firstly, numerical solutions for validation with Newtonian fluids in clean enclosures are conducted, followed by enclosures with heterogeneous medium and centered block. Subsequently, numerical solutions of power-law and Bingham fluids with same enclosures configurations are conducted. The results are reported in the form of streamlines, isotherms and average Nusselt number at hot wall. In general, the heat transfer process in enclosure is governed by Rayleigh number, size and thermal conductivity of the blocks, power index n for power-law fluid and Bingham number for Bingham fluid. Both fluid models are very sensitive with Prandtl number changes. Heterogeneous medium decrease heat transfer in enclosure when affects thermal boundary layer for both fluid models. One analytical prediction was made for power-law fluid. An increase in heat transfer occurs with low thermal conductivity block and few interference and with high thermal conductivity block and great interference, for both fluids.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.

Técnicas de injeção de traçadores são bastante utilizadas nos estudos de escoamentos em meios porosos heterogêneos, principalmente em problemas relacionados à simulação numérica de escoamentos miscíveis em reservatórios de petróleo e à dispersão de contaminantes em aqüíferos. Neste trabalho apresentamos novos algoritmos para a aproximação numérica do problema de injeção de traçadores. Apresentaremos desenvolvimentos recentes do método Forward Integral-Tube Tracking (FIT) que foi originalmente apresentado em Aquino et al. (2007a). O FIT é um método lagrangeano localmente conservativo utilizado na resolução de problemas de transporte linear. Este método não faz o uso de soluções de problemas de Riemann e baseia-se na construção dos tubos integrais introduzidas em Douglas Jr. et al. (2000b). Além disso, ele possui excelente eficiência computacional e é virtualmente livre de difusão numérica. Resultados numéricos são apresentados com o objetivo de comparar a precisão das soluções fornecidas por novas implementações do método FIT na resolução do problema do traçador em reservatórios de petróleo. / The injection of tracers are used in the investigation of flows in heterogeneous porous media, in studies related to the simulation of miscible dispacements in petroleum reservoirs and the dispersion of contaminants in aquifers. In this work we present new algorithms for the numerical approximation of tracer injection problems. We discuss recent developments of the Forward Integral-Tube Tracking (FIT) scheme which was introduced in Aquino et al. (2007a). The FIT is a locally conservative lagrangian scheme for the approximation of the linear transport problems. This scheme does not use analytic solutions of Riemann problems and is based on the construction of the integral tubes introduced in Douglas Jr. et al. (2000b). The FIT scheme is computationally very eficient and is virtually free of numerical diffusion. Numerical results are presented to compare the accuracy of the solutions provided by new implementation of the FIT scheme for the injection of tracers in petroleum reservoirs.