Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Análise da trajetória de poços horizontais em reservatório de óleo pesado

The application of thermal methods, to increase the recovery of heavy oil in mature fields through drainage with multilateral and horizontal wells, has been thoroughly studied, theorically, experimentally, testing new tools and methods. The continuous injection of steam, through a steam injector well and a horizontal producer well in order to improve horizontal sweep of the fluid reservoir, it is an efficient method. Starting from an heterogeneous model, geologically characterized, modeling geostatistics, set history and identification of the best path of permeability, with seismic 3D, has been dubbed a studying model. It was studied horizontal wells in various directions in relation to the steam and the channel of higher permeability, in eight different depths. Into in the same area were studied, the sensitivity of the trajectories of horizontal wells, according to the depth of navigation. With the purpose of obtaining the highest output of oil to a particular flow, quality, temperature and time for the injection of steam. The wells studied showed a significant improvement in the cumulative oil recovery in one of the paths by promoting an alternative to application in mature fields or under development fields with heavy oil

Otimização da injeção cíclica de vapor em reservatório de óleo pesado

Thermal methods made heavy oil production possible in fields where primary recovery failed. Throughout the years steam injection became one of the most important alternatives to increase heavy oil recovery. There are many types of steam injection, and one of them is the cyclic steam injection, which has been used with success in several countries, including Brazil. The process involves three phases: firstly, steam is injected, inside of the producing well; secondly, the well is closed (soak period); and finally, the well is put back into production. These steps constitute one cycle. The cycle is repeated several times until economical production limit is reached. Usually, independent of reservoir type, as the number of cycles increases the cyclic injection turns less efficient. This work aims to analyze rock and reservoir property influence in the cyclic steam injection. The objective was to study the ideal number of cycles and, consequently, process optimization. Simulations were realized using the STARS simulator from the CMG group based in a proposed reservoir model. It was observed that the reservoir thickness was the most important parameter in the process performance, whilst soaking time influence was not significant

Análise da trajetória de poços horizontais em reservatório de óleo pesado

The application of thermal methods, to increase the recovery of heavy oil in mature fields through drainage with multilateral and horizontal wells, has been thoroughly studied, theorically, experimentally, testing new tools and methods. The continuous injection of steam, through a steam injector well and a horizontal producer well in order to improve horizontal sweep of the fluid reservoir, it is an efficient method. Starting from an heterogeneous model, geologically characterized, modeling geostatistics, set history and identification of the best path of permeability, with seismic 3D, has been dubbed a studying model. It was studied horizontal wells in various directions in relation to the steam and the channel of higher permeability, in eight different depths. Into in the same area were studied, the sensitivity of the trajectories of horizontal wells, according to the depth of navigation. With the purpose of obtaining the highest output of oil to a particular flow, quality, temperature and time for the injection of steam. The wells studied showed a significant improvement in the cumulative oil recovery in one of the paths by promoting an alternative to application in mature fields or under development fields with heavy oil

Métodos de simulação para ciclos de Rankine

O presente trabalho é dedicado ao estudo de métodos de simulação para ciclos de Rankine. O trabalho é iniciado com a modelagem de um ciclo de Rankine simples e segue evoluindo para configurações mais complexas tal como o ciclo de Rankine com reaquecimento e regeneração. São adotadas as considerações mais convencionais da prática de projeto de centrais termelétricas cujos sistema térmicos baseiam-se no ciclo de Rankine, incluindo-se queda de pressão em tubulações do circuito além de outras perdas. Em seguida, são estabelecidas as expressões matemáticas que possibilitam a determinação das propriedades termodinâmicas da água em seus mais diversos estados ao longo do ciclo. Por último, são desenvolvidos métodos de simulação, chamados neste trabalho de Substituição Sucessiva e Bloco Único, que caracterizam-se pela resolução simultânea do conjunto de equações algébricas dos ciclos elaborados. As simulações são efetuadas através de programas escritos na linguagem Fortran. Os métodos de simulação são aplicados para a obtenção dos resultados considerados mais importantes na análise de sistemas térmicos de potência, tais como rendimento térmico do ciclo, título na saída da turbina, vazões mássicas pelo sistema, potência nas bombas e calor trocado no gerador de vapor e no condensador Na maioria das simulações, estes resultados apresentam-se como funções da: (1) potência elétrica requerida, eficiência isentrópica e pressões na turbina; (2) eficiência térmica, pressão e temperatura no gerador de vapor; (3) pressão e grau de subresfriamento do líquido saturado no condensador e (4) eficiência isentrópica das bombas. São obtidos os mesmos resultados para os métodos de simulação utilizados. O método da Substituição Sucessiva apresentou menor tempo computacional, principalmente para configurações de ciclo mais complexas. Uma aplicação alternativa do método de Bloco Único demonstrou ser inconveniente para ciclos de configurações mais complexas devido ao elevado tempo computacional, quando todas as equações de cálculo das propriedades termodinâmicas são incluídas no sistema de equações a ser resolvido. Melhores rendimentos térmicos e título na saída da turbina foram obtidos para configurações de ciclo de Rankine com reaquecimento e regeneração.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Análise comparativa de diferentes configurações de poços no processo de combustão in-situ

With an increasing number of mature fields, heavy oil recovery has performed one of the great challenges of the oil industry. The Brazilian Northeast, for example, has numerous heavy oil reservoirs are explored with the use of thermal methods. Among the types of methods used for heavy oil, there is the method of in-situ combustion, a technique in which heat is produced within the container, unlike the injection of heated fluid when the heat is generated at the surface and transported to the reservoir. In this type of process, it is common to use vertical wells as injectors and producers. However, methods which use horizontal wells like oil producers are increasingly studied because of greater contact area between the formation and combustion front. Thus, the main objective of this work was to study the different configurations of wells (CIS THAITM and CAGD) in the process of in-situ combustion in oil recovery using a semi-synthetic tank with Brazilian Northeast features. The method "toe-to-heel air injection" (THAITM) is a process of enhanced oil recovery, which is the integration of in-situ combustion with technological advances in drilling horizontal wells. This method uses horizontal wells such as oil producers, keeping vertical injection wells for injecting air. The oil drain process by differential gravitational assisted with combustion (CAGD) is an integrated, in this configuration the horizontal injector well is drilled at the top formation with a horizontal production well in the lower section. The simulations were performed in a commercial program of thermal processes, called "STARS" (Steam, Thermal, and Advanced Processes Reservoir Simulator), the company CMG (Computer Modelling Group). An analysis of the air flow injection was performed and it was found that each method had a maximum injection to the base model, a show that through this air injection limit was reduced cumulative production of oil. Analyses of operating parameters were used: injection flow, configuration and completion of wells. In the sensitivity analysis we found that the air injection flow showed greater influence on THAI method, since the CIS method the completion of the wells was the most influential parameter and CAGD configuration wells showed the greatest influence in the recovered fraction. The economic results have shown that the best case obtained in CAGD method because, despite having higher initial cost showed the best financial return compared to the best cases the CIS and THAI.

Análise comparativa de diferentes configurações de poços no processo de combustão in-situ

With an increasing number of mature fields, heavy oil recovery has performed one of the great challenges of the oil industry. The Brazilian Northeast, for example, has numerous heavy oil reservoirs are explored with the use of thermal methods. Among the types of methods used for heavy oil, there is the method of in-situ combustion, a technique in which heat is produced within the container, unlike the injection of heated fluid when the heat is generated at the surface and transported to the reservoir. In this type of process, it is common to use vertical wells as injectors and producers. However, methods which use horizontal wells like oil producers are increasingly studied because of greater contact area between the formation and combustion front. Thus, the main objective of this work was to study the different configurations of wells (CIS THAITM and CAGD) in the process of in-situ combustion in oil recovery using a semi-synthetic tank with Brazilian Northeast features. The method "toe-to-heel air injection" (THAITM) is a process of enhanced oil recovery, which is the integration of in-situ combustion with technological advances in drilling horizontal wells. This method uses horizontal wells such as oil producers, keeping vertical injection wells for injecting air. The oil drain process by differential gravitational assisted with combustion (CAGD) is an integrated, in this configuration the horizontal injector well is drilled at the top formation with a horizontal production well in the lower section. The simulations were performed in a commercial program of thermal processes, called "STARS" (Steam, Thermal, and Advanced Processes Reservoir Simulator), the company CMG (Computer Modelling Group). An analysis of the air flow injection was performed and it was found that each method had a maximum injection to the base model, a show that through this air injection limit was reduced cumulative production of oil. Analyses of operating parameters were used: injection flow, configuration and completion of wells. In the sensitivity analysis we found that the air injection flow showed greater influence on THAI method, since the CIS method the completion of the wells was the most influential parameter and CAGD configuration wells showed the greatest influence in the recovered fraction. The economic results have shown that the best case obtained in CAGD method because, despite having higher initial cost showed the best financial return compared to the best cases the CIS and THAI.

Simulações de problemas inversos com aplicações em engenharia nuclear usando técnicas de transporte de partículas neutras monoenergéticas na formulação unidimensional de ordenadas discretas

Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese.

Análise matemática de soluções descontínuas de leis de conservação hiperbólicas e resoluções numéricas para a captura de ondas de choque em escoamentos multifásicos em meios porosos

O processo de recuperação secundária de petróleo é comumente realizado com a injeção de água ou gás no reservatório a fim de manter a pressão necessária para sua extração. Para que o investimento seja viável, os gastos com a extração precisam ser menores do que o retorno financeiro obtido com a produção de petróleo. Objetivando-se estudar possíveis cenários para o processo de exploração, costuma-se utilizar simulações dos processos de extração. As equações que modelam esse processo de recuperação são de caráter hiperbólico e não lineares, as quais podem ser interpretadas como Leis de Conservação, cujas soluções são complexas por suas naturezas descontínuas. Essas descontinuidades ou saltos são conhecidas como ondas de choque. Neste trabalho foi abordada uma análise matemática para os fenômenos oriundos de leis de conservação, para em seguida utilizá-la no entendimento do referido problema. Foram estudadas soluções fracas que, fisicamente, podem ser interpretadas como ondas de choque ou rarefação, então, para que fossem distinguidas as fisicamente admissíveis, foi considerado o princípio de entropia, nas suas diversas formas. As simulações foram realizadas nos âmbitos dos escoamentos bifásicos e trifásicos, em que os fluidos são imiscíveis e os efeitos gravitacionais e difusivos, devido à pressão capilar, foram desprezados. Inicialmente, foi feito um estudo comparativo de resoluções numéricas na captura de ondas de choque em escoamento bifásico água-óleo. Neste estudo destacam-se o método Composto LWLF-k, o esquema NonStandard e a introdução da nova função de renormalização para o esquema NonStandard, onde obteve resultados satisfatórios, principalmente em regiões onde a viscosidade do óleo é muito maior do que a viscosidade da água. No escoamento bidimensional, um novo método é proposto, partindo de uma generalização do esquema NonStandard unidimensional. Por fim, é feita uma adaptação dos métodos LWLF-4 e NonStandard para a simulação em escoamentos trifásicos em domínios unidimensional. O esquema NonStandard foi considerado mais eficiente nos problemas abordados, uma vez que sua versão bidimensional mostrou-se satisfatória na captura de ondas de choque em escoamentos bifásicos em meios porosos.

Simulação computacional do comportamento de estruturas de aço sob incêndio

A análise do comportamento estrutural sob incêndio constitui uma parte importante da engenharia de proteção contra incêndio, especialmente no caso de estruturas de aço, por sua alta condutividade térmica e relativa esbeltez das seções. As dificuldades econômicas e práticas, associadas à avaliação do comportamento estrutural, por meio de ensaios em escala real, têm estimulado o desenvolvimento e uso de métodos de simulação numérica. Esta tese trata da simulação numérica do comportamento de estruturas de aço sob condições de incêndio e se divide em três partes. As duas primeiras partes foram desenvolvidas na Universidade de Liége, na Bélgica, usando-se o programa SAFIR como ferramenta numérica. A terceira parte foi desenvolvida de forma independente, na Universidade Federal do Rio Grande do Sul. Na primeira parte, é feito um estudo comparativo entre o uso de elementos finitos de viga e de casca, na modelagem de vigas simplesmente apoiadas, sujeitas a flambagem lateral por torção. Os esforços de torção, presentes no caso de flambagem lateral, podem levar a uma plastificação da seção transversal e à conseqüente redução da rigidez à torção da seção. Da mesma forma, a degradação das propriedades do material, com o aumento da temperatura, leva à redução da rigidez à torção Havia dúvidas se o modelo com elementos de viga, com uma rigidez à torção constante, poderia fornecer uma resposta aceitável. O estudo mostrou que uma resposta com boa precisão pode ser conseguida, usando-se elementos de viga, desde que o módulo de elasticidade transversal seja ajustado para refletir seu valor correspondente à temperatura de falha. Isso implica um processo iterativo, uma vez que a temperatura de falha não é previamente conhecida. Por outro lado, a degradação da rigidez à torção, por efeitos de plastificação, pode ser ignorada. Na segunda parte, é feita a comparação entre as modelagens bidimensional e tridimensional, de galpões industriais de um andar, sob incêndio. Comumente, a estrutura de galpões industriais é composta por pórticos-tipo, dispostos em paralelo. A análise desses galpões é comumente feita pela simulação no plano do pórtico de aço ou, simplesmente, da treliça da cobertura Na análise bidimensional, importantes efeitos fora do plano são ignorados, como a redistribuição de esforços devido à degradação do material ou à expansão térmica, ou instabilidade lateral dos elementos. A importância desses efeitos e a adequabilidade do modelo 2D para representar o comportamento real são discutidas. Na terceira parte, um modelo numérico para a simulação tridimensional do comportamento de estruturas de aço sob incêndio é apresentado. O modelo é baseado no conceito de rótulas plásticas generalizadas, com modificações para melhor representar a formação e expansão da plastificação no elemento. A descrição cinemática adotada permite obter bons resultados, mesmo com o uso de poucos elementos. A determinação do vetor de esforços internos no elemento, incluindo os efeitos da temperatura, é detalhada. O procedimento foi validado por comparação com ensaios e com modelos numéricos mais sofisticados, como o programa SAFIR. Os resultados demonstram que o procedimento proposto pode ser usado como uma forma alternativa de análise 3D de estruturas sob incêndio, com precisão razoável e baixo esforço computacional.

Proposição de um modelo para análise de fluxo de carga, curto-circuito e simulação de relés de produção

The objective of the present work is develop a model to simulate electrical energy networks in transient and stead states, using the software ATP (Alternative Transient Program), able to be a way to join two distinct themes, present in classical methodology planning networks: short circuit analysis and load flow theory. Beyond that, using a tool for relay simulation, this paper intend to use the new developed model to investigate the influence of transient phenomenon in operation of protection relays, and calibrate the enterprise's protections relays. For testing the model, some relays, actually, installed at COSERN were used

Remoção de óleo da água de produção por flotação em coluna utilizando tensoativos de origem vegetal

In the petroleum industry, water is always present in the reservoir formation together with petroleum and natural gas and this fact provokes the production of water with petroleum, resulting in a great environmental impact. Several methods can be applied for treatment of oily waters, such as: gravitational vases, granulated media filtration systems, flotation process, centrifugation process and the use of hydrocyclones, which can also be used in a combined way. However, the flotation process has showed a great efficiency as compared with other methods, because these methods do not remove great part of the emulsified oil. In this work was investigated the use of surfactants derived from vegetable oils, OSS and OGS, as collectors, using the flotation process in a glass column with a porous plate filter in its base for the input of the gaseous steam. For this purpose, oil/water emulsions were prepared using mechanical stirring, with concentrations around 300 ppm. The air flow rate was set at 700 cm3/min and the porous plate filter used for the generation of the air bubbles has pore size varying from 16 to 40 Pm. The column operated at constant volume (1500mL). A new methodology has been developed to collect the samples, where, instead of collecting the water phase, it was collected the oil phase removed by the process in the top of the flotation column. It has been observed that it is necessary to find an optimum surfactant concentration to achieve enhanced removal efficiency. Being for OSS 1.275 mmol/L and for OGS 0.840 mmol/L, with removal efficiencies of 93% and 99%, respectively, using synthetic solutions. For the produced water, the removal in these concentrations was 75% for OSS and 65% for OGS. It is possible to remove oil from water in a flotation process using surfactants of high HLB, fact that is against the own definition of HLB (Hydrophile-Lipophile Balance). The interfacial tension is an important factor in the oil removal process using a flotation process, because it has direct interference in the coalescence of the oil drops. The spreading of the oil of the air bubble should be considered in the process, and for the optimum surfactant concentrations it reached a maximum value. The removal kinetics for the flotation process using surfactants in the optimum concentration has been adjusted according to a first order model, for synthetic water as for the produced water.

Tipologia e freqüência de métodos de ensino em disciplinas de formação profissionalizante em engenharia de produção

A globalização dos mercados trouxe, para o cenário brasileiro, impactos políticos, socioecônomicos e culturais enfatizando a figura de um homem, denominado cidadão do mundo. Devido às próprias limitações humanas, em um mundo em que as mudanças são constantes, esse cidadão tem como característica ser um eterno aprendiz. A universidade, locus de produção de identidade do país e cuja função social é formar indivíduos neste novo contexto global, se encontra em um dilema existencial: pesquisa e ensino; a superação dessa dicotomia poderá derivar de novas relações humanas na formação do sujeito-coletivo que constrói sua identidade através das possibilidades de diversidades do mundo que o circunda. Estando na aprendizagem, para a aquisição de conhecimentos ou o desenvolvimento de habilidades e atitudes pelas experiências educativas, o papel do docente como mediador em organizar estratégias (métodos e técnicas de ensino) para que o discente conheça e crie a cultura. O objetivo desta pesquisa é mapear a tipologia e freqüência de métodos de ensino em disciplinas de formação profissionalizante em engenharia de produção. Identificar, através de um estudo teórico e revisão bibliográfica, os tipos e métodos mais utilizados no ensino superior do curso de engenharia de produção, bem como as vantagens, limitações e habilidades desenvolvidas segundo a utilização adequada dessas técnicas. No estudo de caso, se analisam os métodos de ensino e as competências necessárias nas disciplinas profissionalizantes do curso de engenharia de produção mecânica da escola de engenharia de São Carlos cujos métodos de coleta foram através de questionário para os discentes e entrevistas junto aos docentes. Embora haja uma variedade de técnicas socializantes que deveriam atender aos objetivos específicos de cada disciplina, os avanços tornam-se despercebidos pela descrição da técnica. Ao menosprezar a conduta e observação docente (dimensão técnica/ética/política), durante o processo, o docente abdica de sua autoridade no processo de ensino/aprendizagem; as limitações, como a aula expositiva tradicional (centrando o conhecimento no professor), têm se tornado fonte de estagnação, desfavorecendo a amplitude do desenvolvimento do potencial humano e fazendo com que o discente não se torne co-autor de seu processo de aprender a aprender, pois limita sua fonte de investigação, criatividade e anseio de desafio.

Estudo das forças capilar, viscosa e gravitacional através de grupos adimensionais na simulação da drenagem gravitacional assistida por gás (GAGD)

The distribution and mobilization of fluid in a porous medium depend on the capillary, gravity, and viscous forces. In oil field, the processes of enhanced oil recovery involve change and importance of these forces to increase the oil recovery factor. In the case of gas assisted gravity drainage (GAGD) process is important to understand the physical mechanisms to mobilize oil through the interaction of these forces. For this reason, several authors have developed physical models in laboratory and core floods of GAGD to study the performance of these forces through dimensionless groups. These models showed conclusive results. However, numerical simulation models have not been used for this type of study. Therefore, the objective of this work is to study the performance of capillary, viscous and gravity forces on GAGD process and its influence on the oil recovery factor through a 2D numerical simulation model. To analyze the interplay of these forces, dimensionless groups reported in the literature have been used such as Capillary Number (Nc), Bond number (Nb) and Gravity Number (Ng). This was done to determine the effectiveness of each force related to the other one. A comparison of the results obtained from the numerical simulation was also carried out with the results reported in the literature. The results showed that before breakthrough time, the lower is the injection flow rate, oil recovery is increased by capillary force, and after breakthrough time, the higher is the injection flow rate, oil recovery is increased by gravity force. A good relationship was found between the results obtained in this research with those published in the literature. The simulation results indicated that before the gas breakthrough, higher oil recoveries were obtained at lower Nc and Nb and, after the gas breakthrough, higher oil recoveries were obtained at lower Ng. The numerical models are consistent with the reported results in the literature