An efficient hybrid local nonmatching method for multiphase flow simulations in heterogeneous fractured media

This paper presents simulation methodology that combines a local nonmatching grid with a discrete fracture model. Designed for 2D and 3D multiphase flow simulations in standard simulators, the method handles matrix–matrix, fracture–fracture, and matrix–fracture connections in the context of an unstr...

Full description

Saved in:
Bibliographic Details
Published in:Engineering with computers 2015-04, Vol.31 (2), p.347-360
Main Author: Mustapha, Hussein
Format: Article
Language:English
Subjects:
Accuracy
CAE) and Design
Calculus of Variations and Optimal Control
Optimization
Classical Mechanics
Computer Science
Computer simulation
Computer-Aided Engineering (CAD
Control
Fracture mechanics
Joints
Math. Applications in Chemistry
Mathematical and Computational Engineering
Mathematical models
Multiphase flow
Original Article
Simulators
Systems Theory
Three dimensional
Two dimensional
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This paper presents simulation methodology that combines a local nonmatching grid with a discrete fracture model. Designed for 2D and 3D multiphase flow simulations in standard simulators, the method handles matrix–matrix, fracture–fracture, and matrix–fracture connections in the context of an unstructured, local nonmatching grid. The grid is generated at the fracture intersections, enabling accurate modeling of small control volumes between connecting fractures. Grids are obtained simply by redistributing the volume of small control volumes surrounding the small control volumes, making the method computationally efficient. A unified method to calculate the interblock transmissibility is used for both matching and nonmatching mesh. An unstructured finite-volume graph-based reservoir simulator with a two-point flux approximation reads the new grid by making a simple modification to the graph of connections between the control volumes. The method requires no special treatment of fracture–fracture or matrix–fracture transmissibility calculations and has the flexibility to simulate any flow problem efficiently. Several 2D and 3D numerical examples demonstrate the method’s performance and accuracy. Both simple and complex fracture configurations are presented with various levels of geologic and fluid complexity. The numerical results are in good agreement with those of a reference solution obtained on a finely structured grid.
ISSN:0177-0667
1435-5663
DOI:10.1007/s00366-014-0355-0