Stability of Piston-Like Displacements of Water by Steam and Nitrogen in Porous Media

1982 ◽  
Vol 22 (05) ◽  
pp. 625-634 ◽  
Author(s):  
David A. Krueger

Abstract Downhole steam generation leads to consideration of reservoir fluid displacement by a mixture of steam and nitrogen. The linear stability analysis of the steam condensation front has been generalized to include a noncondensing gas. Roughly speaking, the addition of nitrogen increases the likelihood of having fingers, but, compared with the no-nitrogen case, the fingers will grow more slowly. Introduction The theory of the stability of flows through porous media has been a subject of interest for more than 25 years, dating back to the pioneering work of Dietz, Chuoke et al., and Saffman and Taylor. They considered injecting one fluid (e.g., water) to force a second fluid (e.g., oil) out of a porous medium. The primary result was that instabilities (fingering) occurred when the driving fluid was more mobile than the driven fluid. Hagoort included multiple fluid phases. Miller generalized the original work to include steam driving water (liquid). He showed that the thermodynamic phase transition (steam to water) introduces two stabilizing effects. The first effect introduces a water/steam velocity ratio as a multiplier of the mobility ratio. This factor is less than one because of the volume change upon condensation. The second effect is the cooling of incipient steam fingers by the surrounding water, which retards their growth. Baker anticipated these effects in a qualitative way to explain his experiments, which showed a more stable displacement by steam than was expected on the basis of mobility ratios alone. Armento and Miller also have considered the stability of the in-situ combustion front in porous media. Their work deals with a region where steam is generated. This paper reformulates Miller's results for a condensation front in a more useful form including general numerical results and extends the theory to include injection of a noncondensing gas (e.g., nitrogen) together with the steam. Depending on the particular situation, the presence of nitrogen can be either stabilizing or destabilizing. The motivation for the generalization comes from enhanced oil recovery projects where the exhaust gases from the steam generator are injected into the reservoir along with the steam. This paper considers perturbations on a flat condensation front that is perpendicular to its velocity. The gravitational force along this velocity is included, but the component of the gravitational force perpendicular to the velocity is not. Thus we include the effect of gravity on fingering, but we do not discuss the gravity override problem. In Stability Analysis we present two steps:determination of the motion of a flat condensation front (details are in the Appendix) andevaluation of the characteristic time for growth or decay of a perturbation of that front. In Results wegive the results for a specific reservoir;discuss the sensitivity of these results to the important reservoir parameters (flow velocities and absolute permeabilities),show that, if surface tension and gravitation are unimportant, the stability condition is independent of the absolute permeability and absolute flow rates, anddiscuss the longest wavelength for a stable perturbation. In the final section we discuss the main conclusions. Stability Analysis We consider a homogeneous porous medium with fluids in two regions as illustrated in Fig. 1. A steam/nitrogen mixture is injected at the left, and water (liquid) and nitrogen are produced at the fight. The linear stability analysis proceeds in two main stages and follows the general methods as discussed by Chandrasekhar and the specific application of Miller. First, we assume that the condensation front is flat, moves with constant velocity, v, and has properties that vary with z alone. SPEJ P. 625^

Author(s):  
K.B. Tsiberkin

The stability of incompressible fluid plane-parallel flow over a layer of a saturated porous medium is studied. The results of a linear stability analysis are described at different porosity values. The considered system is bounded by solid wall from the porous layer bottom. Top fluid surface is free and rigid. A linear stability analysis of plane-parallel stationary flow is presented. It is realized for parameter area where the neutral stability curves are bimodal. The porosity variation effect on flow stability is considered. It is shown that there is a transition between two main instability modes: long-wave and short-wave. The long-wave instability mechanism is determined by inflection points within the velocity profile. The short-wave instability is due to the large transverse gradient of flow velocity near the interface between liquid and porous medium. Porosity decrease stabilizes the long wave perturbations without significant shift of the critical wavenumber. Simultaneously, the short-wave perturbations destabilize, and their critical wavenumber changes in wide range. When the porosity is less than 0.7, the inertial terms in filtration equation and magnitude of the viscous stress near the interface increase to such an extent that the Kelvin-Helmholtz analogue of instability becomes the dominant mechanism for instability development. The stability band realizes in narrow porosity area. It separates the two branches of the neutral curve.


Author(s):  
Cody S. Dowd ◽  
Joseph W. Meadows

Lean premixed (LPM) combustion systems are susceptible to thermoacoustic instability, which occurs when acoustic pressure oscillations are in phase with the unsteady heat release rates. Porous media has inherent acoustic damping properties, and has been shown to mitigate thermoacoustic instability; however, theoretical models for predicting thermoacoustic instability with porous media do not exist. In the present study, a 1-D model has been developed for the linear stability analysis of the longitudinal modes for a series of constant cross-sectional area ducts with porous media using a n-Tau flame transfer function. By studying the linear regime, the prediction of acoustic growth rates and subsequently the stability of the system is possible. A transfer matrix approach is used to solve for acoustic perturbations of pressure and velocity, stability growth rate, and frequency shift without and with porous media. The Galerkin approximation is used to approximate the stability growth rate and frequency shift, and it is compared to the numerical solution of the governing equations. Porous media is modeled using the following properties: porosity, flow resistivity, effective bulk modulus, and structure factor. The properties of porous media are systematically varied to determine the impact on the eigenfrequencies and stability growth rates. Porous media is shown to increase the stability domain for a range of time delays (Tau) compared to similar cases without porous media.


2018 ◽  
Vol 141 (4) ◽  
Author(s):  
Cody S. Dowd ◽  
Joseph W. Meadows

Lean premixed (LPM) combustion systems are susceptible to thermoacoustic instability, which occurs when acoustic pressure oscillations are in phase with the unsteady heat release rates. Porous media has inherent acoustic damping properties and has been shown to mitigate thermoacoustic instability; however, theoretical models for predicting thermoacoustic instability with porous media do not exist. In the present study, a one-dimensional (1D) model has been developed for the linear stability analysis of the longitudinal modes for a series of constant cross-sectional area ducts with porous media using a n-Tau flame transfer function (FTF). By studying the linear regime, the prediction of acoustic growth rates and subsequently the stability of the system is possible. A transfer matrix approach is used to solve for acoustic perturbations of pressure and velocity, stability growth rate, and frequency shift without and with porous media. The Galerkin approximation is used to approximate the stability growth rate and frequency shift, and it is compared to the numerical solution of the governing equations. Porous media is modeled using the following properties: porosity, flow resistivity, effective bulk modulus, and structure factor. The properties of porous media are systematically varied to determine the impact on the eigenfrequencies and stability growth rates. Porous media is shown to increase the stability domain for a range of time delays (Tau) compared to similar cases without porous media.


2008 ◽  
Vol 612 ◽  
pp. 1-19 ◽  
Author(s):  
P. A. LAKSHMI NARAYANA ◽  
P. V. S. N. MURTHY ◽  
RAMA SUBBA REDDY GORLA

The stability of Soret-driven thermosolutal convection in a shallow horizontal layer of a porous medium subjected to inclined thermal and solutal gradients of finite magnitude is investigated theoretically by means of a linear stability analysis. The horizontal components of these gradients induce a Hadley circulation, which becomes unstable when vertical components are sufficiently large. We employed a two-term Galerkin approximation for various modes of instability. The effect of the Soret parameter on the mechanism of instability of the thermosolutal convection is investigated. Results are presented for various values of the governing parameters of the flow. It is observed that the Soret parameter has a significant effect on convective instability and this is discussed.


2001 ◽  
Author(s):  
Davide Valtorta ◽  
Khaled E. Zaazaa ◽  
Ahmed A. Shabana ◽  
Jalil R. Sany

Abstract The lateral stability of railroad vehicles travelling on tangent tracks is one of the important problems that has been the subject of extensive research since the nineteenth century. Early detailed studies of this problem in the twentieth century are the work of Carter and Rocard on the stability of locomotives. The linear theory for the lateral stability analysis has been extensively used in the past and can give good results under certain operating conditions. In this paper, the results obtained using a linear stability analysis are compared with the results obtained using a general nonlinear multibody methodology. In the linear stability analysis, the sources of the instability are investigated using Liapunov’s linear theory and the eigenvalue analysis for a simple wheelset model on a tangent track. The effects of the stiffness of the primary and secondary suspensions on the stability results are investigated. The results obtained for the simple model using the linear approach are compared with the results obtained using a new nonlinear multibody based constrained wheel/rail contact formulation. This comparative numerical study can be used to validate the use of the constrained wheel/rail contact formulation in the study of lateral stability. Similar studies can be used in the future to define the limitations of the linear theory under general operating conditions.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1980 ◽  
Vol 96 (2) ◽  
pp. 375-393 ◽  
Author(s):  
R. Mckibbin ◽  
M. J. O'Sullivan

The formalism required to determine the criterion for the onset of convection in a multi-layered porous medium heated from below is developed using a straightforward linear stability analysis. Detailed results for two- and three-layer configurations are presented. These results show that large permeability differences between the layers are required to force the system into an onset mode different from a homogeneous system.


1997 ◽  
Vol 352 ◽  
pp. 265-281 ◽  
Author(s):  
A. M. H. BROOKER ◽  
J. C. PATTERSON ◽  
S. W. ARMFIELD

A non-parallel linear stability analysis which utilizes the assumptions made in the parabolized stability equations is applied to the buoyancy-driven flow in a differentially heated cavity. Numerical integration of the complete Navier–Stokes and energy equations is used to validate the non-parallel theory by introducing an oscillatory heat input at the upstream end of the boundary layer. In this way the stability properties are obtained by analysing the evolution of the resulting disturbances. The solutions show that the spatial growth rate and wavenumber are highly dependent on the transverse location and the disturbance flow quantity under consideration. The local solution to the parabolized stability equations accurately predicts the wave properties observed in the direct simulation whereas conventional parallel stability analysis overpredicts the spatial amplification and the wavenumber.


2013 ◽  
Vol 736 ◽  
pp. 464-494 ◽  
Author(s):  
P. Pearce ◽  
J. Daou

AbstractWe investigate the Rayleigh–Bénard convection problem within the context of a diffusion flame formed in a horizontal channel where the fuel and oxidizer concentrations are prescribed at the porous walls. This problem seems to have received no attention in the literature. When formulated in the low-Mach-number approximation the model depends on two main non-dimensional parameters, the Rayleigh number and the Damköhler number, which govern gravitational strength and reaction speed respectively. In the steady state the system admits a planar diffusion flame solution; the aim is to find the critical Rayleigh number at which this solution becomes unstable to infinitesimal perturbations. In the Boussinesq approximation, a linear stability analysis reduces the system to a matrix equation with a solution comparable to that of the well-studied non-reactive case of Rayleigh–Bénard convection with a hot lower boundary. The planar Burke–Schumann diffusion flame, which has been previously considered unconditionally stable in studies disregarding gravity, is shown to become unstable when the Rayleigh number exceeds a critical value. A numerical treatment is performed to test the effects of compressibility and finite chemistry on the stability of the system. For weak values of the thermal expansion coefficient $\alpha $, the numerical results show strong agreement with those of the linear stability analysis. It is found that as $\alpha $ increases to a more realistic value the system becomes considerably more stable, and also exhibits hysteresis at the onset of instability. Finally, a reduction in the Damköhler number is found to decrease the stability of the system.


Sign in / Sign up

Export Citation Format

Share Document