Microcomputer Gas Reservoir Simulation Using Finite Element Methods

Author(s):  
R.W. Logan ◽  
R.L. Lee ◽  
M.R. Tek
1981 ◽  
Vol 21 (01) ◽  
pp. 115-128 ◽  
Author(s):  
Larry C. Young

Abstract Several previous studies have applied finite-element methods to reservoir simulation problems. Accurate solutions have been demonstrated with these methods; however, competitiveness with finite difference has not been established for most nonlinear reservoir simulation problems. In this study a more efficient finite-element procedures is presented and tested. The method is Galerkin-based, and improved efficiency is obtained by combining Lagrange trial functions with Lobatto quadrature in a particular way. The simulation of tracer performance, ion exchange preflush performance, and adverse mobility ratio miscible displacements is considered. For the problems considered, the method is shown to yield accurate solutions with less computing expense than finite differences or previously proposed finite-element techniques. For the special case of linear trial functions, the method reduces to a five-point central difference approximation. In contrast to previously reported results, this approximation is found to simulate adverse mobility ratio displacements without grid orientation sensitivity, provided a sufficiently fine grid is used. Introduction In the past few years several studies have investigated the use of finite-element methods in reservoir simulation. These include single-phase two-component simulations in one1–3 and two4,5 spatial dimensions and two-phase immiscible calculations in both one6–8 and two9,10 dimensions. These studies have demonstrated that the method is capable of giving accurate solutions, particularly for small slug problems and adverse mobility ratio displacements. All these studies used what we term conventional Galerkin finite-element techniques,1 and, unfortunately, these methods have not proved to be cost competitive with finite differences for most nonlinear reservoir simulation problems. A reduction in computing requirements is, therefore, necessary to make finite-element methods truly useful for reservoir simulation. Relative to finite differences, the increased computing requirements of conventional Galerkin-based methods are due to the following.The approximation of time-derivative terms involves the same number of surrounding grid points as the approximation of flow terms; thus, implicit-pressure/explicit-saturation (IMPES) techniques are not possible (see Ref. 12, Chap. 7).The matrices which result from the approximation of flow terms are not nearly so sparse as in finite differences; thus, the solution of matrix problems requires more computation.The computational work required to generate matrix coefficients is considerably greater than with finite differences due to the number of numerical integrations which must be performed.


2013 ◽  
Vol 41 (2) ◽  
pp. 127-151
Author(s):  
Rudolf F. Bauer

ABSTRACT The benefits of a tire's equilibrium profile have been suggested by several authors in the published literature, and mathematical procedures were developed that represented well the behavior of bias ply tires. However, for modern belted radial ply tires, and particularly those with a lower aspect ratio, the tire constructions are much more complicated and pose new problems for a mathematical analysis. Solutions to these problems are presented in this paper, and for a modern radial touring tire the equilibrium profile was calculated together with the mold profile to produce such tires. Some construction modifications were then applied to these tires to render their profiles “nonequilibrium.” Finite element methods were used to analyze for stress concentrations and deformations within all tires that did or did not conform to equilibrium profiles. Finally, tires were built and tested to verify the predictions of these analyses. From the analysis of internal stresses and deformations on inflation and loading and from the actual tire tests, the superior durability of tires with an equilibrium profile was established, and hence it is concluded that an equilibrium profile is a beneficial property of modern belted radial ply tires.


1983 ◽  
Author(s):  
W. HABASHI ◽  
M. HAFEZ ◽  
P. KOTIUGA

Sign in / Sign up

Export Citation Format

Share Document