InterPore2023

Direct Numerical Simulation of weak-inertia single-phase flow in porous materials using SPH

25 May 2023, 15:30

1h 30m

Poster Presentation (MS09) Pore-scale modelling Poster

Speaker

Mr David Krach (University of Stuttgart)

Description

Instrinsic permeability of a porous material is a crucial material parameter in various application fields like e.g. geosciences, materials science and mechanical engineering. In these disciplines intrinsic permeability is a widley used input parameter for numerical simulations in the framework of continuum-scale models that employ Darcy’s law. However, the assumption of linearity in constitutive relations for the momentum exchange is subjected to restrictions and only valid in the so-called creeping flow regime (Reynolds number Re < 1). We present an approach focusing on 3-D pore-scale-resolved computations of single-phase fluid flow through porous media from moderate to higher Re-numbers (1 < Re < 1000) aiming to show a smooth transition of effective properties from creeping flow to the weak inertia regime. While this has already been shown for 2-D artificial domains [4], we seek to generalize the phenomena by investigating spherepackings and natural porous materials with small to moderate porosities (φ ≤ 0.2). For representative 3-D simulations based on XRCT-scans, with voxels in the order of 10003 massively parallel direct numerical simulation methods are required. Therefore, we choose fully-Lagrangian Smoothed Particle Hydrodynamics (SPH) as a simulation method to model pore-scale-resolved flow by means of the weakly compressible Navier-Stokes equations. The solver is implemented on top of the software framework HOOMD-Blue [1, 2] since this allows for massively parallel CPU and GPU computations. A sufficient scaling behavior as well as the numerical accuracy in the Darcy regime is demonstrated [3].

References

[1] J. A. Anderson, C. D. Lorenz, and A. Travesset. General purpose molecular dynamics simulations fully implemented on graphics processing units. Journal of Computational Physics, 227(10):5342–5359, 2008.
[2] J. Glaser, T. D. Nguyen, J. A. Anderson, P. Lui, F. Spiga, J. A. Millan, D. C. Morse, and S. C. Glotzer. Strong scaling of general-purpose molecular dynamics simulations on GPUs. Computer Physics Communications, 192:97–107, 2015.
[3] M. Osorno, M. Schirwon, N. Kijanski, R. Sivanesapillai, H. Steeb, and D. Göddeke. A cross-platform, high-performance sph toolkit for image-based flow simulations on the pore scale of porous media. Computer Physics Communications, 267:108059, 2021.
[4] R. Sivanesapillai, H. Steeb, and A. Hartmaier. Transition of effective hydraulic properties from low to high reynolds number flow in porous media. Geophysical Research Letters, 41(14):4920–4928, 2014.