In this dissertation, the evolution of a pressure wave driven by a harmonic signal on the boundary during gas combustion is studied. The problem is modeled by a nonlinear, hyperbolic partial differential equation. Steady-state behavior is investigated using the perturbation method to ensure that enough time has passed for any transient effects to have dissipated. The zeroth, first and second-order perturbation solutions are obtained and their moduli are plotted against frequency. It is seen that the first and second-order corrections have unique maxima that shift to the right as the frequency decreases and to the left as the frequency increases. Dispersion relations are determined and their limiting behavior investigated in the low and high frequency regimes. It is seen that for low frequencies, the medium assumes a diffusive-like nature. However, for high frequencies the medium behaves similarly to one exhibiting relaxation. The phase speed is determined and its limiting behavior examined. For low frequencies, the phase speed is approximately equal to sqrt[ω/(n+1)] and for high frequencies, it behaves as 1/(n+1), where n is the mode number. Additionally, a maximum allowable value of the perturbation parameter, ε = 0.8, is determined that ensures boundedness of the solution. The location of the peak of the first-order correction, xmax, as a function of frequency is determined and is seen to approach the limiting value of 0.828/sqrt(ω) as the frequency tends to zero and the constant value of 2 ln 2 as the frequency tends to infinity. Analytic expressions are obtained for the approximate general perturbation solution in the low and high-frequency regimes and are plotted together with the perturbation solution in the corresponding frequency regimes, where the agreement is seen to be excellent. Finally, the solution obtained from the perturbation method is compared with the long-time solution obtained by the finite-difference scheme; again, ensuring that the transient effects have dissipated. Since the finite-difference scheme requires a right boundary, its location is chosen so that the wave dissipates in amplitude enough so that any reflections from the boundary will be negligible. The perturbation solution and the finite-difference solution are found to be in excellent agreement. Thus, the validity of the perturbation method is established.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / 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. / The process of secondary oil recovery is commonly accomplished by injecting water or gas into the reservoir to maintain the necessary pressure for their extraction. So that the investment is viable spending extraction must be smaller than the financial return to oil production. Aiming to study possible scenarios for the exploration process, it is customary to use simulations of extraction processes. The equations that model this process of recovery are hyperbolic and nonlinear, which can be interpreted as Conservation Laws , whose solutions are complex by their discontinuous nature . These discontinuities or jumps are known as shock waves. Due to this importance, this work will be discussed a mathematical analysis of the phenomena arising from conservation laws, to then use it in the understanding of this problem. Weak solutions that physically can be interpreted as shock waves or rarefaction, so that they might be distinguished physically admissible were studied, was considered the principle of entropy, in its various forms. The simulations were performed in the fields of two-phase and three-phase flow, in which the fluids are immiscible and gravitational and diffusive effects due to capillary pressure, were discarded. Initially a comparative study of numerical resolutions in the capture of shock waves in water-oil two-phase flow was made. This study highlights LWLF k Composite method and Nonstandard. Was also presented a new renormalization function for nonstandard scheme with satisfactory results, especially in regions where the oil viscosity is much higher than the viscosity of the water. In twodimensional flow, a new method will be presented. The same is a generalization of onedimensional nonstandard schema. Finally, the adaptation of nonstandard and LWLF-4 methods for simulating in three-phase one-dimensional flows. In general, the nonstandard scheme was considered the most efficient method in problems addressed, since its twodimensional version was satisfactory in capturing shock waves in two-phase flow in porous media.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / 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. / The process of secondary oil recovery is commonly accomplished by injecting water or gas into the reservoir to maintain the necessary pressure for their extraction. So that the investment is viable spending extraction must be smaller than the financial return to oil production. Aiming to study possible scenarios for the exploration process, it is customary to use simulations of extraction processes. The equations that model this process of recovery are hyperbolic and nonlinear, which can be interpreted as Conservation Laws , whose solutions are complex by their discontinuous nature . These discontinuities or jumps are known as shock waves. Due to this importance, this work will be discussed a mathematical analysis of the phenomena arising from conservation laws, to then use it in the understanding of this problem. Weak solutions that physically can be interpreted as shock waves or rarefaction, so that they might be distinguished physically admissible were studied, was considered the principle of entropy, in its various forms. The simulations were performed in the fields of two-phase and three-phase flow, in which the fluids are immiscible and gravitational and diffusive effects due to capillary pressure, were discarded. Initially a comparative study of numerical resolutions in the capture of shock waves in water-oil two-phase flow was made. This study highlights LWLF k Composite method and Nonstandard. Was also presented a new renormalization function for nonstandard scheme with satisfactory results, especially in regions where the oil viscosity is much higher than the viscosity of the water. In twodimensional flow, a new method will be presented. The same is a generalization of onedimensional nonstandard schema. Finally, the adaptation of nonstandard and LWLF-4 methods for simulating in three-phase one-dimensional flows. In general, the nonstandard scheme was considered the most efficient method in problems addressed, since its twodimensional version was satisfactory in capturing shock waves in two-phase flow in porous media.