Year: 2009
Communications in Computational Physics, Vol. 6 (2009), Iss. 1 : pp. 72–84
Abstract
We describe an operator splitting technique based on physics rather than on dimension for the numerical solution of a nonlinear system of partial differential equations which models three-phase flow through heterogeneous porous media. The model for three-phase flow considered in this work takes into account capillary forces, general relations for the relative permeability functions and variable porosity and permeability fields. In our numerical procedure a high resolution, nonoscillatory, second order, conservative central difference scheme is used for the approximation of the nonlinear system of hyperbolic conservation laws modeling the convective transport of the fluid phases. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the parabolic and elliptic problems associated with the diffusive transport of fluid phases and the pressure-velocity problem. This numerical procedure has been used to investigate the existence and stability of nonclassical shock waves (called transitional or undercompressive shock waves) in two-dimensional heterogeneous flows, thereby extending previous results for one-dimensional flow problems. Numerical experiments indicate that the operator splitting technique discussed here leads to computational efficiency and accurate numerical results.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/2009-CiCP-7672
Communications in Computational Physics, Vol. 6 (2009), Iss. 1 : pp. 72–84
Published online: 2009-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 13
Keywords: Operator splitting three-phase flow heterogeneous porous media central differencing schemes mixed finite elements.