Computing Profiles of Porosity and Surface/Volume Ratio in a Packed Bed Using Monte Carlo Method
Porosity and surface/volume ratio are two important parameters for a packed bed. In cylindrical packed beds at low tube-to-particle diameter ratio, they vary greatly in the radial direction. In the existing literature, radial profiles of porosity and surface/volume ratio have been computed using some analytical equations which involve elliptic integrals. In this paper, a Monte Carlo method is used to compute those profiles. To the authors’ knowledge, the method has never been employed in this context. The procedure of using this novel method is explained in detail. Through a rigorous error analysis based on statistics, the accuracy of the simulation result can be controlled. Before any simulation, the number of sampling points needed in the Monte Carlo simulation can be determined given an accuracy requirement. The method is completely general and can be used to compute profiles of porosity and surface/volume ratio in any packed bed with any shape of packing elements.