A sequential explicit-implicit algorithm for computing discontinuous flows in porous media
- Univ. of California, Berkeley, CA (United States)
A novel numerical algorithm for computing incompressible, discontinuous, two-phase flows in two-dimensional, inhomogeneous, and isotropic porous media is presented. The algorithm uses Colella et al`s. hybrid sequential explicit-implicit approach for both accuracy and speed of the calculations. The explicit part uses a high-order Godunov scheme with a modified Van Leer geometrical slope limiter, similar to those used in shock dynamics. The implicit part is a two-step solver: the first step is a Crank-Nicolson saturation solver and the second one is a Poisson solver for the phase pressure. Both use fast multilevel multigrid solvers with the number of operations of the order of {var_theta}[Nlog(N)], where N is the number of grid points. For an implicit simulator, the number of operations is {var_theta}(N{sup 3}) per time step. Two numerically stiff reservoir engineering problems are presented to demonstrate the low numerical dispersion and second-order accuracy of our method.
- OSTI ID:
- 468280
- Report Number(s):
- CONF-961003-; TRN: 96:006577-0085
- Resource Relation:
- Conference: 71. annual technical conference and exhibition of the Society of Petroleum Engineers, Denver, CO (United States), 6-9 Oct 1996; Other Information: PBD: 1996; Related Information: Is Part Of 1996 SPE annual technical conference and exhibition: Production operations and engineering/general; PB: 892 p.
- Country of Publication:
- United States
- Language:
- English
Similar Records
Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. II. Implicit Methods
Implicit–explicit multirate infinitesimal stage-restart methods