Pore-to-core simulations of flow with large velocities using continuum models and imaging data

Public Deposited

Downloadable Content

Download PDF


Attribute NameValues
  • We consider computational modeling of flow with small and large velocities at porescale and at corescale, and we address various challenges in simulation, upscaling, and modeling. While our focus is on voxel-based data sets from real porous media imaging, our methodology is verified first on synthetic geometries, and we analyze various scaling and convergence properties. We show that the choice of a voxel-based grid and REV size can lead up to 10-20% difference in calculated conductivities. On the other hand, the conductivities decrease significantly with flow rates, starting in a regime usually associated with the onset of inertia effects. This is accompanied by deteriorating porescale solver performance, and we continue our experiments up until about 50% reduction in conductivities, i.e., to Reynolds number just under 1. To account for this decrease, we propose a practical power-based fully anisotropic non-Darcy model at corescale for which we calculate the parameters by upscaling.
  • Keywords: Upscaling, Inertia effects, Anisotropy, Forchheimer model, Flow in porous media, 76S05, 76M45, Navier–Stokes equations, Convergence, Porescale simulations, 76M50, Non-Darcy flow
Resource Type
Date Available
Date Issued
  • Peszynska, M., & Trykozko, A. (2013). Pore-to-core simulations of flow with large velocities using continuum models and imaging data. Computational Geosciences, 17(4), 623-645. doi:10.1007/s10596-013-9344-4
Journal Title
Journal Volume
  • 17
Journal Issue/Number
  • 4
Academic Affiliation
Rights Statement
Funding Statement (additional comments about funding)
  • M. Peszynska was partially supported by the grant NSF DMS-1115827, and A. Trykozko was in part supported by PL-Grid infrastructure. This research was carried out with the support of the ”HPC Infrastructure for Grand Challenges of Science and Engineering” Project, co-financed by the European Regional Development Fund under the Innovative Economy Operational Programme. Computations on cluster halo2 were performed with grant nr G35-12.
Peer Reviewed



This work has no parents.