A Parallel Second-Order Adaptive Mesh Algorithm for Incompressible Flow in Porous Media
In this paper we present a second-order accurate adaptive algorithm for solving multiphase, incompressible flows in porous media. We assume a multiphase form of Darcy's law with relative permeabilities given as a function of the phase saturation. The remaining equations express conservation of mass for the fluid constituents. In this setting the total velocity, defined to be the sum of the phase velocities, is divergence-free. The basic integration method is based on a total-velocity splitting approach in which we solve a second-order elliptic pressure equation to obtain a total velocity. This total velocity is then used to recast component conservation equations as nonlinear hyperbolic equations. Our approach to adaptive refinement uses a nested hierarchy of logically rectangular grids with simultaneous refinement of the grids in both space and time. The integration algorithm on the grid hierarchy is a recursive procedure in which coarse grids are advanced in time, fine grids areadvanced multiple steps to reach the same time as the coarse grids and the data atdifferent levels are then synchronized. The single grid algorithm is described briefly,but the emphasis here is on the time-stepping procedure for the adaptive hierarchy. Numerical examples are presented to demonstrate the algorithm's accuracy and convergence properties and to illustrate the behavior of the method.
- Research Organization:
- Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, CA (US)
- Sponsoring Organization:
- Computing Sciences Directorate
- DOE Contract Number:
- AC02-05CH11231
- OSTI ID:
- 957033
- Report Number(s):
- LBNL-176E
- Journal Information:
- Proceedings of the Royal Society, Journal Name: Proceedings of the Royal Society Journal Issue: 1907 Vol. 367; ISSN 1364-503X
- Country of Publication:
- United States
- Language:
- English
Similar Records
An adaptive projection method for the incompressible Euler equations
An adaptive projection method for the incompressible Navier-Stokes equations