This study presents and calibrates a Discrete Element Method (DEM) contact model for wet granular materials in the pendular regime. The model extends a previously calibrated dry contact formulation by incorporating liquid bridges that generate capillary adhesion between particles, while liquid migration is represented through evolving bridge volumes. Two reactor-grade polypropylene powders with different particle size distributions, bulk densities, and surface morphologies are investigated, resulting in distinct wetting behavior. A schematic framework is introduced to relate increasing liquid content to the transition from dry to wet contacts using two key parameters: the minimum liquid film volume and the maximum liquid bridge volume. These parameters are calibrated using dynamic angle of repose measurements from rotating drum experiments. The calibrated model reproduces the experimental flow behavior of both powders: full agreement is obtained for the coarser, more porous powder across all liquid contents, while for the finer, denser powder, agreement is achieved at low to moderate liquid contents. At higher liquid contents, discrepancies arise due to agglomeration effects amplified by particle scaling. These results demonstrate the effectiveness of the dynamic angle of repose as a calibration target and highlight the limitations of particle scaling for strongly cohesive wet granular systems. The proposed framework provides a practical basis for DEM-based modeling of wet powder flow in rotating and continuously sheared industrial systems operating in the pendular regime.

In this paper we introduce the open-source code MercuryDPM: a code for simulating discrete particles. The paper discusses software and management issues that may be interesting for the developers of other open-source codes. Then we review the new features that have been added since the last publication: an improved Hertz-Mindlin model; a new liquid bridge model of Lian and Seville; a droplet-spray model; better support for re-creating complex, measured particle size distributions; a new implementation of rigid clumps; an implementation of elastic membranes; a wear model for walls; a soft-kill feature and a cloud-deployment interface for AWS.

We introduce the open-source package MercuryDPM, which we have been developing over the last few years. MercuryDPM is a code for discrete particle simulations. It simulates the motion of particles by applying forces and torques that stem either from external body forces, (gravity, magnetic fields, etc.) or particle interactions. The code has been developed extensively for granular applications, and in this case these are typically (elastic, plastic, viscous, frictional) contact forces or (adhesive) short-range forces. However, it could be adapted to include long-range (molecular, self-gravity) interactions as well. MercuryDPM is an object-oriented algorithm with an easy-to-use user interface and a flexible core, allowing developers to quickly add new features. It is parallelised using MPI and released under the BSD 3-clause licence. Its open-source developers’ community has developed many features, including moving and curved walls; state-of-the-art granular contact models; specialised classes for common geometries; non-spherical particles; general interfaces; restarting; visualisation; a large self-test suite; extensive documentation; and numerous tutorials and demos. In addition, MercuryDPM has three major components that were originally invented and developed by its team: an advanced contact detection method, which allows for the first time large simulations with wide size distributions; curved (non-triangulated) walls; and multicomponent, spatial and temporal coarse-graining, a novel way to extract continuum fields from discrete particle systems. We illustrate these tools and a selection of other MercuryDPM features via various applications, including size-driven segregation down inclined planes, rotating drums, and dosing silos. Program summary: Program Title: MercuryDPM Program Files doi: http://dx.doi.org/10.17632/n7jmdrdc52.1 Licensing provisions: BSD 3-Clause Programming language: C++, Fortran Supplementary material: http://mercurydpm.org Nature of problem: Simulation of granular materials, i.e. conglomerations of discrete, macroscopic particles. The interaction between individual grains is characterised by a loss of energy, making the behaviour of granular materials distinct from atomistic materials, i.e. solids, liquids and gases. Solution method: MercuryDPM (Thornton et al., 2013, 2019; Weinhart et al., 2016, 2017, 2019) is an implementation of the Discrete Particle Method (DPM), also known as the Discrete Element Method (DEM) (Cundall and Strack, 1979). It simulates the motion of individual particles by applying forces and torques that stem either from external forces (gravity, magnetic fields, etc.) or from particle-pair and particle–wall interactions (typically elastic, plastic, dissipative, frictional, and adhesive contact forces). DPM simulations have been successfully used to understand the many unique granular phenomena – sudden phase transitions, jamming, force localisation, etc. – that cannot be explained without considering the granular microstructure. Unusual features: MercuryDPM was designed ab initio with the aim of allowing the simulation of realistic geometries and materials found in industrial and geotechnical applications. It thus contains several bespoke features invented by the MercuryDPM team: (i) a neighbourhood detection algorithm (Krijgsman et al., 2014) that can efficiently simulate highly polydisperse packings, which are common in industry; (ii) curved walls (Weinhart et al., 2016) making it possible to model real industrial geometries exactly, without triangulation errors; and (iii) MercuryCG (Weinhart et al., 2012, 2013, 2016; Tunuguntla et al., 2016), a state-of-the-art analysis tool that extracts local continuum fields, providing accurate analytical/rheological information often not available from experiments or pilot plants. It further contains a large range of contact models to simulate complex interactions such as elasto-plastic deformation (Luding, 2008), sintering (Fuchs et al., 2017), melting (Weinhart et al., 2019), breaking, wet and dry cohesion (Roy et al., 2016, 2017), and liquid migration (Roy et al., 2018), all of which have important industrial applications.