Porous Media Simulations

Parallel domain decomposition method with non-blocking communication for flow through porous media

A. Lemmer, R. Hilfer

Journal of Computational Physics 281, 970-981 (2015)

submitted on
Thursday, December 19, 2013

This paper introduces a domain decomposition method for numerically solving the Stokes equation for very large, complex geometries. Examples arise from realistic porous media. The computational method is based on the SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm which uses a finite-differences approach for discretizing the underlying equations. It achieves comparable speed and efficiency as lattice Boltzmann methods. The domain decomposition method splits a large three-dimensional region into slices that can be processed in parallel on multi-processor computation environments with only minimal communication between the computation nodes. With this method, the flow through a porous medium with grid sizes up to 2048 x 2048 x 2048 voxel has been calculated.

For more information see