10:40
Modelling Porous Media (Part 2)
Chair: Iuliu Sorin Pop
10:40
25 mins
|
Dynamic Multilevel Simulation of Fractured Reservoirs with Generic Conductivity Contrasts
Mousa HosseiniMehr, Cornelis Vuik, Hadi Hajibeygi
Abstract: A dynamic multilevel method for simulation of flow in fractured reservoirs with Projection-based Embedded Discrete Fracture Model (pEDFM) is presented. The proposed method develops a scalable and efficient simulation framework for large-scale heterogeneous fractured reservoirs. At the same time, pEDFM provides a consistent formulation approach to capture complex flow physics in presence of explicit fractures which can act from highly conductive manifolds to flow barriers.
pEDFM formulation of [Tene et al., 2017] is extended to account for generic 3D fracture geometries, within a robust fully-coupled fully-implicit scheme. This fine-scale system is mapped into a dynamic multilevel system (Cusini et al., 2016, HosseiniMehr et al., 2018), with automatically determined resolution to minimize the error. Sequences of restriction and prolongation operators are developed to generate an accurate map between finescale and hierarchical multilevel multiscale systems. These multiscale operators allow for avoiding up-scaling when constructing coarser-resolution systems, and as such lead to significantly more accurate results. Local mass is conserved across all resolutions by using finite-volume restriction operators.
The performance of our method is studied for a range of 2D and 3D test cases including fractures with higher and lower conductivities compared with that of the matrix. It is shown that our method provides accurate and scalable simulation framework for highly fractured heterogeneous reservoirs. Importantly, we study the sensitivity of the production estimations with respect to the fracture density, orientation and properties. In addition, the scalability of the simulation approach for large-scale geo-models is discussed. Our results demonstrate that the proposed method is capable of simulating large-scale fractured reservoir models.
References:
[Cusini et al., 2016] Cusini, M., van Kruijsdijk, C., and Hajibeygi, H. (2016). Algebraic dynamic multilevel (ADM) method for fully implicit simulations of multiphase flow in porous media. Journal of Computational Physics, 314:60-79.
[HosseiniMehr et al., 2018] HosseiniMehr, M., Cusini, M., Vuik, C., and Hajibeygi, H. (2018). Algebraic dynamic multilevel method for embedded discrete fracture model (F-ADM). Journal of Computational Physics, 373:324-345.
[Tene et al., 2017] Tene, M., Bosma, S. B., Kobaisi, M. S. A., and Hajibeygi, H. (2017). Projection-based embedded discrete fracture model (pEDFM). Advances in Water Resources, 105:205-216.
|
11:05
25 mins
|
Global random walk solution for flow and transport in porous media
Nicolae Suciu
Abstract: Global random walk (GRW) algorithms, unconditionally stable and free of numerical
diffusion, solve parabolic partial differential equations by moving computational particles on regular lattices according to random walk rules [1, 2]. Solutions of the 1D elliptic equation for the pressure head in saturated porous media can be obtained as steady-state limit of the staggered finite difference (FD) scheme
|
11:30
25 mins
|
Modelling of the influence of vegetative barrier on particulate matter concentration using OpenFOAM
Hynek Reznicek
Abstract: High concentration of atmospheric dust is a well known risk factor to
human health. Vegetative barriers are one of the most popular ways how to substan-
tially reduce the high pollution concentration. Correct modelling of air flow inside
the Atmospheric Boundary Layer (ABL) is essential to accurately predict concentra-
tion of the passive scalar (dust). The question whether the CFD toolbox OpenFOAM
is capable of modelling of this type of problems is tested in the contribution. The
results obtained from OpenFOAM were compared simultaneously with the experi-
mental data and the CFD results of the program Atifes, developed at CTU for ABL
simulations. It is shown that the recommended setting of OpenFOAM’s atmospheric
library has several limitations. Special attention is paid to the different wall functions
used in both solvers and the differences are discussed.
|
|
