科研之友
微信
新浪微博
Facebook
分享链接
摘要
We present a computationally fast Invasion Percolation (IP) algorithm. IP is a numerical approach for generating realistic fluid distributions for quasi-static (i.e., slow) immiscible fluid invasion in porous media. The algorithm proposed here uses a binary-tree data structure to identify the site (pore) connected to the invasion cluster that is the next to be invaded. Gravity is included. Trapping is not explicitly treated in the numerical examples but can be added, for example, using a Hoshen-Kopelman algorithm. Computation time to percolation for a 3D system having total sites and invaded sites at percolation goes as for the proposed binary-tree algorithm and as for a standard implementation of IP that searches through all of the uninvaded sites at each step. The relation between and is , where is the fractal dimension of an infinite cluster and is Euclidean space dimension. In numerical practice, on finite-sized cubic lattices with invasion structures influenced by the injection boundary and boundary conditions lateral to the flow direction, we observe the scaling in 3D (valid through the second decimal place) instead of based on the infinite cluster fractal dimension D = 2.53.
