“…Models for this coupling are commonly based on continuum mechanics combined with a discrete representation of cracks [1,2,3,4,5]. Alternatively, discrete approaches, such as discrete element method, lattice and network models, have been proposed to model these processes [6,7,8,9,10,11,12,13,14,15].One network approach, based on the Delaunay tessellation of a random set of points, has been shown to be suitable for modelling fracture [16,17,18,19] and mass transport [20], providing mesh insensitive results. In this approach, the physical processes are modelled by a multi-dimensional network of one-dimensional elements, which are placed on the Delaunay edges ( Figure 1a); the element properties are determined by the corresponding Voronoi tessellation.…”