Geometrically realistic particles are getting more and more available in discrete (DEM) simulations thanks to a variety of dedicated DEM implementations, enabling one to properly address the known influence of particle shape in geomechanics [1]. In addition to the multi-sphere approach, whereby an appropriate number of overlapping spheres are clumped together and act as one rigid body, one can mention the possibility to define a particle's surface as a polyhedra with its set of vertices, edges and facets or through a Level Set description [2]. The latter relies on a discrete field for the signed distance function to the surface at hand, whose zero level set implicitly describes that surface. For the purpose of contact treatment, a surface discretization in terms of boundary nodes is deduced from the distance field and also enters the method as a second ingredient.