4.7 Article

Higher order FE-FV method on unstructured grids for transport and two-phase flow with variable viscosity in heterogeneous porous media

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 241, Issue -, Pages 416-444

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2012.12.017

Keywords

Finite-element finite volume method; Flow in porous media; Transport; Slope limiter; Monotone upwind scheme for conservation laws; Unstructured grid

Funding

  1. ExxonMobil Research Alliance 'Fundamental Controls of Flow in Carbonates' (FC)2
  2. CMG Reservoir Simulation Foundation
  3. Edinburgh Collaborative of Subsurface Science and Engineering, a joint research institute of the Edinburgh Research Partnership in Engineering and Mathematics (ECOSSE)

Ask authors/readers for more resources

This paper presents higher order methods for the numerical modeling of two-phase flow with simultaneous transport and adsorption of viscosifying species within the individual phases in permeable porous media. The numerical scheme presented addresses the three major challenges in simulating this process. Firstly, the component transport is strongly coupled with the viscous and capillary forces that act on the movement of the carrier phase. The discretization of the capillary parts is especially difficult since its effect on flow yields non-linear parabolic conservation equations. These are amenable to non-linear finite elements (FEs), while the capillary contribution on the component transport is first-order hyperbolic, where classical FEs are unsuitable. We solve this efficiently by a Strang splitting that uses finite volumes (FVs) with explicit time-stepping for the viscous parts and a combined finite element-finite volume (FEFV) scheme with implicit time-stepping for the capillary parts. Secondly, the components undergo hydrodynamic dispersion and discerning between numerical and physical dispersion is essential. We develop higher-order formulations for the phase and component fluxes that keep numerical dispersion low and combine them with implicit FEs such that the non-linearities of the dispersion tensor are fully incorporated. Thirdly, subsurface permeable media show strong spatial heterogeneity, with coefficients varying over many orders of magnitude and geometric complexity that make the use of unstructured grids essential. In this work, we employ node-centered FVs that combine their ability to resolve flow with the flexibility of FEs. Numerical examples of increasing complexity are presented that demonstrate the convergence and robustness of our approach and prove its versatility for highly heterogeneous, and geometrically complex fractured porous media. (C) 2013 Elsevier Inc. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available