14-17 May 2018
New Orleans
US/Central timezone

Algebraic Dynamic Multilevel (ADM) method for flow in heterogeneous porous media with embedded discrete fractures

16 May 2018, 16:16
New Orleans

New Orleans

Oral 20 Minutes MS 2.13: Advances in numerical modelling of multiphase flow and transport in fractured porous media Parallel 8-C


Matteo Cusini (Delft University of Technology)


An Algebraic Dynamic Multilevel (ADM) method [1, 2] for fully implicit simulations of multiphase flow in heterogeneous fractured natural porous media is presented. The fine-scale fully-implicit (FIM) system is obtained following the Embedded Discrete Fracture Modelling (EDFM) [3] approach. A set of nested coarse grids at different resolutions (or levels) is constructed independently for each medium. At each time-step the fine-scale FIM system is mapped to a dynamically selected grid formed by grid-blocks of the previously defined grids. The grid resolution is chosen, independently for each medium, based on a front- tracking criterion that ensures that fine-scale resolution is employed only where most physical interactions take place (i.e., moving saturation front). Mapping between different grid resolutions is performed by employing sequences of restriction and prolongation (interpolation) operators. Finite-volume restriction [4] operators are employed to ensure mass conservation whereas different interpolation strategies are employed for each variable. In particular, multiscale basis functions [5] are considered as pressure interpolators to ensure an accurate interpolation of the pressure field inside coarse blocks.
The efficiency and accuracy of the proposed algorithm is shown through a set of challenging 2D and 3D test cases involving various non-linear physics. The sensitivity of the method to the choice of the coarsening and interpolation strategy is also presented.


[1] Matteo Cusini, Cor van Kruijsdijk, and Hadi Hajibeygi. Algebraic dynamic multilevel (adm) method for fully implicit simulations of multiphase flow in porous media. J. Comput. Phys., 314:60–79, 2016.
[2] Matteo Cusini, Barnaby Fryer, Cor van Kruijsdijk, and Hadi Hajibeygi. Algebraic dynamic multilevel method for compositional flow in heterogeneous porous media. J. Comput. Phys., 354:593-612, 2018.
[3] SH Lee, CL Jensen, MF Lough, et al. An efficient finite difference model for flow in a reservoir with multiple length-scale fractures. In SPE Annual Technical Conference and Exhibition. Society of Petroleum Engineers, 1999.
[4] P. Jenny, S. H. Lee, and H. A. Tchelepi. Multi-scale finite-volume method forelliptic problems in subsurface flow simulation. J. Comput. Phys., 187:47–67, 2003.
[5] T. Y. Hou and X.-H. Wu. A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys., 134:169–189, 1997.

Acceptance of Terms and Conditions Click here to agree

Primary authors

Mr Mousa HosseiniMehr (Delft University of Technology) Matteo Cusini (Delft University of Technology) Prof. Hadi Hajibeygi (TU Delft)

Presentation Materials

There are no materials yet.