An Improved Porosity Calculation Algorithm for Particle Flow Code

The widely used discrete-element particle flow software PFC’s (PFC 7.0 and previous versions) algorithm for calculating porosity is not sufficiently accurate. Because of this, when the particles are densely packed, the solution to the equation produces an algorithm exception for odd calculations of porosity, which results in the inability to calculate the results. This paper, based on a Darcy seepage model of fluid flow through a granular bed, analyzed the shortcomings of the two porosity calculation methods of PFC and the function analysis method. Combining this analysis with the theory of computer graphics, a new and efficient porosity calculation algorithm was proposed. The result showed that the new proposed porosity calculation algorithm calculated a more accurate and reasonable porosity field and made the iterative solution of the CFD equation more stable. This method makes porosity-related models of PFC more accurate. The algorithm can be not only used to calculate porosity, but also applied to other fields.