A. Lemmer, R. Hilfer
Journal of Computational Physics 281, 970-981 (2015)
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