Comparative analysis of porosity coarse-graining techniques for discrete element simulations of dense particulate systems

Moris Kalderon,Edward Smith,Catherine O’Sullivan
DOI: https://doi.org/10.1007/s40571-021-00402-4
2021-06-23
Computational Particle Mechanics
Abstract:The discrete element method (DEM) is a well-established approach to study granular materials in numerous fields of application; each granular particle is modelled individually to predict the overall behaviour. This behaviour can be then extracted by averaging, or coarse graining, the sample using a suitable method. The choice of appropriate coarse-graining method entails a compromise between accuracy and computational cost, especially in the large-scale simulations typically required by industry. A number of coarse-graining methods have been proposed in the literature, and these are reviewed and categorized in this work. Within this contribution, two novel porosity coarse-graining strategies are proposed including a voxel method where a secondary dense grid of "pixel cells" is implemented adopting a binary logic for the coarse graining and a hybrid method where both analytical formulas and pixels are utilized. The proposed methods are compared with four coarse-graining schemes that have been documented in the literature, including the particle centroid method, an analytical method, a method which solves the diffusion equation and an approach which employs averaging using kernels. The novel methods are validated for problems in both two and three dimensions through comparison with the "accurate" analytical method. It is shown that, once validated, both the proposed schemes can approximate the exact solutions quite accurately; however, there is a high computational cost associated with the voxel method. The accuracy of both methods can be adjusted allowing the user to decide between accuracy and computational time. A detailed comparison is then presented for all six schemes considering "accuracy", "smoothness" and "computational cost". Optimal parameters are obtained for all six methods, and recommendations for coarse-graining DEM samples are discussed.
mechanics,mathematics, interdisciplinary applications
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to effectively coarsen the porosity of particle systems in Discrete Element Method (DEM) simulations for use in Computational Fluid Dynamics (CFD) - DEM coupled simulations. Specifically, the paper focuses on how to reduce computational costs while ensuring accuracy, maintain the conservation of the total mass of the solid phase, and reduce large - scale fluctuations (noise) in numerical solutions. The paper addresses these issues by comparing existing porosity coarsening techniques and proposing two new porosity coarsening strategies, namely the voxel - based method and the hybrid method. ### Background of the Paper The Discrete Element Method (DEM) is a widely used method for studying granular materials. It predicts the overall behavior by modeling each particle individually. In order to extract useful information from these simulation data, an appropriate coarsening method, i.e., sample averaging, is usually required. Selecting an appropriate coarsening method requires a trade - off between accuracy and computational cost, especially in large - scale simulations common in industrial applications. ### Existing Methods The paper reviews several porosity coarsening methods proposed in the literature, which are mainly divided into two categories: 1. **Binning Methods**: This method allocates the particle volume to Eulerian grid cells (bins) and calculates the porosity independently within each cell. 2. **Smoothing Methods**: This method generates an Eulerian field and uses Gaussian or Lucy functions to weight the particle positions, thereby distributing the particle density over the support domain of the entire weight function. ### New Methods The paper proposes two new porosity coarsening strategies: 1. **Voxel Method**: A dense secondary grid is introduced, and binary logic is used to determine whether each voxel is covered by particles. By calculating the volume of voxels covered by particles, it is allocated to the corresponding Eulerian grid cell. 2. **Hybrid Method**: When circular / spherical particles overlap only one cell face or are completely within an Eulerian cell, an analytical method is used to calculate the area / volume; otherwise, when spherical particles overlap multiple cell faces, the voxel - summing method is used for calculation. ### Method Evaluation The paper makes a detailed comparison of these six methods (two existing methods and four new methods), considering accuracy, smoothness, and computational cost. Through the verification of 2 - D and 3 - D problems, the paper shows the advantages of the new methods in terms of accuracy, but also points out the high computational cost of the voxel - based method. ### Conclusion Finally, the paper provides the optimal parameters for each method and discusses the recommended methods in DEM sample coarsening. Through these analyses, the paper provides valuable references for researchers when choosing suitable porosity coarsening methods.