-
papaya2: 2D Irreducible Minkowski Tensor computation
Authors:
Fabian M. Schaller,
Jenny Wagner,
Sebastian C. Kapfer
Abstract:
A common challenge in scientific and technical domains is the quantitative description of geometries and shapes, e.g. in the analysis of microscope imagery or astronomical observation data. Frequently, it is desirable to go beyond scalar shape metrics such as porosity and surface to volume ratios because the samples are anisotropic or because direction-dependent quantities such as conductances or…
▽ More
A common challenge in scientific and technical domains is the quantitative description of geometries and shapes, e.g. in the analysis of microscope imagery or astronomical observation data. Frequently, it is desirable to go beyond scalar shape metrics such as porosity and surface to volume ratios because the samples are anisotropic or because direction-dependent quantities such as conductances or elasticity are of interest. Minkowski Tensors are a systematic family of versatile and robust higher-order shape descriptors that allow for shape characterization of arbitrary order and promise a path to systematic structure-function relationships for direction-dependent properties. Papaya2 is a software to calculate 2D higher-order shape metrics with a library interface, support for Irreducible Minkowski Tensors and interpolated marching squares. Extensions to Matlab, JavaScript and Python are provided as well. While the tensor of inertia is computed by many tools, we are not aware of other open-source software which provides higher-rank shape characterization in 2D.
△ Less
Submitted 28 October, 2020;
originally announced October 2020.
-
Pomelo, a tool for computing Generic Set Voronoi Diagrams of Aspherical Particles of Arbitrary Shape
Authors:
Simon Weis,
Philipp W. A. Schönhöfer,
Fabian M. Schaller,
Matthias Schröter,
Gerd E. Schröder-Turk
Abstract:
We describe the development of a new software tool, called "Pomelo", for the calculation of Set Voronoi diagrams. Voronoi diagrams are a spatial partition of the space around the particles into separate Voronoi cells, e.g. applicable to granular materials. A generalization of the conventional Voronoi diagram for points or monodisperse spheres is the Set Voronoi diagram, also known as navigational…
▽ More
We describe the development of a new software tool, called "Pomelo", for the calculation of Set Voronoi diagrams. Voronoi diagrams are a spatial partition of the space around the particles into separate Voronoi cells, e.g. applicable to granular materials. A generalization of the conventional Voronoi diagram for points or monodisperse spheres is the Set Voronoi diagram, also known as navigational map or tessellation by zone of influence. In this construction, a Set Voronoi cell contains the volume that is closer to the surface of one particle than to the surface of any other particle. This is required for aspherical or polydisperse systems.
Pomelo is designed to be easy to use and as generic as possible. It directly supports common particle shapes and offers a generic mode, which allows to deal with any type of particles that can be described mathematically. Pomelo can create output in different standard formats, which allows direct visualization and further processing. Finally, we describe three applications of the Set Voronoi code in granular and soft matter physics, namely the problem of packings of ellipsoidal particles with varying degrees of particle-particle friction, mechanical stable packings of tetrahedra and a model for liquid crystal systems of particles with shapes reminiscent of pears
△ Less
Submitted 10 March, 2017;
originally announced March 2017.
-
Cuddling Ellipsoids: Densest local structures of uniaxial ellipsoids
Authors:
Fabian M. Schaller,
Robert F. B. Weigel,
Sebastian C. Kapfer
Abstract:
Connecting the collective behavior of disordered systems with local structure on the particle scale is an important challenge, for example in granular and glassy systems. Compounding complexity, in many scientific and industrial applications, particles are polydisperse, aspherical or even of varying shape. Here, we investigate a generalization of the classical kissing problem in order to understan…
▽ More
Connecting the collective behavior of disordered systems with local structure on the particle scale is an important challenge, for example in granular and glassy systems. Compounding complexity, in many scientific and industrial applications, particles are polydisperse, aspherical or even of varying shape. Here, we investigate a generalization of the classical kissing problem in order to understand the local building blocks of packings of aspherical grains. We numerically determine the densest local structures of uniaxial ellipsoids by minimizing the Set Voronoi cell volume around a given particle. Depending on the particle aspect ratio, different local structures are observed and classified by symmetry and Voronoi coordination number. In extended disordered packings of frictionless particles, knowledge of the densest structures allows to rescale the Voronoi volume distributions onto the single-parameter family of $k$-Gamma distributions. Moreover, we find that approximate icosahedral clusters are found in random packings, while the optimal local structures for more aspherical particles are not formed.
△ Less
Submitted 13 October, 2016;
originally announced October 2016.
-
Cell shape analysis of random tessellations based on Minkowski tensors
Authors:
Michael A. Klatt,
Günter Last,
Klaus Mecke,
Claudia Redenbach,
Fabian M. Schaller,
Gerd E. Schröder-Turk
Abstract:
To which degree are shape indices of individual cells of a tessellation characteristic for the stochastic process that generates them? Within the context of stochastic geometry and the physics of disordered materials, this corresponds to the question of relationships between different stochastic models. In the context of image analysis of synthetic and biological materials, this question is centra…
▽ More
To which degree are shape indices of individual cells of a tessellation characteristic for the stochastic process that generates them? Within the context of stochastic geometry and the physics of disordered materials, this corresponds to the question of relationships between different stochastic models. In the context of image analysis of synthetic and biological materials, this question is central to the problem of inferring information about formation processes from spatial measurements of resulting random structures. We address this question by a theory-based simulation study of shape indices derived from Minkowski tensors for a variety of tessellation models. We focus on the relationship between two indices: an isoperimetric ratio of the empirical averages of cell volume and area and the cell elongation quantified by eigenvalue ratios of interfacial Minkowski tensors. Simulation data for these quantities, as well as for distributions thereof and for correlations of cell shape and volume, are presented for Voronoi mosaics of the Poisson point process, determinantal and permanental point processes, and Gibbs hard-core and random sequential absorption processes as well as for Laguerre tessellations of polydisperse spheres and STIT- and Poisson hyperplane tessellations. These data are complemented by mechanically stable crystalline sphere and disordered ellipsoid packings and area-minimising foam models. We find that shape indices of individual cells are not sufficient to unambiguously identify the generating process even amongst this limited set of processes. However, we identify significant differences of the shape indices between many of these tessellation models. Given a realization of a tessellation, these shape indices can narrow the choice of possible generating processes, providing a powerful tool which can be further strengthened by density-resolved volume-shape correlations.
△ Less
Submitted 24 June, 2016;
originally announced June 2016.
-
Non-universal Voronoi cell shapes in amorphous ellipsoid packings
Authors:
Fabian M. Schaller,
Sebastian C. Kapfer,
James E. Hilton,
Paul W. Cleary,
Klaus Mecke,
Cristiano De Michele,
Tanja Schilling,
Mohammad Saadatfar,
Matthias Schröter,
Gary W. Delaney,
Gerd E. Schröder-Turk
Abstract:
In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure plays a pivotal role in determining the macroscopic physical properties. Here, we analyse local structure metrics derived from the Voronoi diagram of configurations of oblate ellipsoids, for various aspect ratios $α$ and global volume fractions $φ_g$. We focus on jammed static configurat…
▽ More
In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure plays a pivotal role in determining the macroscopic physical properties. Here, we analyse local structure metrics derived from the Voronoi diagram of configurations of oblate ellipsoids, for various aspect ratios $α$ and global volume fractions $φ_g$. We focus on jammed static configurations of frictional ellipsoids, obtained by tomographic imaging and by discrete element method simulations. In particular, we consider the local packing fraction $φ_l$, defined as the particle's volume divided by its Voronoi cell volume. We find that the probability $P(φ_l)$ for a Voronoi cell to have a given local packing fraction shows the same scaling behaviour as function of $φ_g$ as observed for random sphere packs. Surprisingly, this scaling behaviour is further found to be independent of the particle aspect ratio. By contrast, the typical Voronoi cell shape, quantified by the Minkowski tensor anisotropy index $β=β_0^{2,0}$, points towards a significant difference between random packings of spheres and those of oblate ellipsoids. While the average cell shape $β$ of all cells with a given value of $φ_l$ is very similar in dense and loose jammed sphere packings, the structure of dense and loose ellipsoid packings differs substantially such that this does not hold true. This non-universality has implications for our understanding of jamming of aspherical particles.
△ Less
Submitted 10 April, 2015;
originally announced April 2015.
-
Local analysis of the history dependence in tetrahedra packings
Authors:
N. Nirmal Thyagu,
Max Neudecker,
Simon Weis,
Fabian M. Schaller,
Matthias Schröter
Abstract:
The mechanical properties of a granular sample depend frequently on the way the packing was prepared. However, is not well understood which properties of the packing store this information. Here we present an X-ray tomography study of three pairs of tetrahedra packings prepared with three different tapping protocols. The packings in each pair differs in the number of mechanical constraints $C$ imp…
▽ More
The mechanical properties of a granular sample depend frequently on the way the packing was prepared. However, is not well understood which properties of the packing store this information. Here we present an X-ray tomography study of three pairs of tetrahedra packings prepared with three different tapping protocols. The packings in each pair differs in the number of mechanical constraints $C$ imposed on the particles by their contacts, while their bulk volume fraction $φ_{global}$ is approximately the same. We decompose $C$ into the contributions of the three different contact types possible between tetrahedra -- face-to-face (F2F), edge-to-face (E2F), and point contacts -- which each fix a different amount of constraints. We then perform a local analysis of the contact distribution by grouping the particles together according to their individual volume fraction $φ_{local}$ computed from a Voronoi tessellation. We find that in samples which have been tapped sufficiently long the number of F2F contacts becomes an universal function of $φ_{local}$. In contrast the number of E2F and point contacts varies with the applied tapping protocol. Moreover, we find that the anisotropy of the shape of the Voronoi cells depends on the tapping protocol. This behavior differs from spheres and ellipsoids and posses a significant constraint for any mean-field approach to tetrahedra packings.
△ Less
Submitted 1 March, 2018; v1 submitted 19 January, 2015;
originally announced January 2015.
-
Local origin of global contact numbers in frictional ellipsoid packings
Authors:
Fabian M. Schaller,
Max Neudecker,
Mohammad Saadatfar,
Gary Delaney,
Gerd E. Schröder-Turk,
Matthias Schröter
Abstract:
In particulate soft matter systems the average number of contacts $Z$ of a particle is an important predictor of the mechanical properties of the system. Using X-ray tomography, we analyze packings of frictional, oblate ellipsoids of various aspect ratios $α$, prepared at different global volume fractions $φ_g$. We find that $Z$ is a monotonously increasing function of $φ_g$ for all $α$. We demons…
▽ More
In particulate soft matter systems the average number of contacts $Z$ of a particle is an important predictor of the mechanical properties of the system. Using X-ray tomography, we analyze packings of frictional, oblate ellipsoids of various aspect ratios $α$, prepared at different global volume fractions $φ_g$. We find that $Z$ is a monotonously increasing function of $φ_g$ for all $α$. We demonstrate that this functional dependence can be explained by a local analysis where each particle is described by its local volume fraction $φ_l$ computed from a Voronoi tessellation. $Z$ can be expressed as an integral over all values of $φ_l$: $Z(φ_g, α, X) = \int Z_l (φ_l, α, X) \; P(φ_l | φ_g) \; dφ_l$. The local contact number function $ Z_l (φ_l, α, X)$ describes the relevant physics in term of locally defined variables only, including possible higher order terms $X$. The conditional probability $P(φ_l | φ_g)$ to find a specific value of $φ_l$ given a global packing fraction $φ_g$ is found to be independent of $α$ and $X$. Our results demonstrate that for frictional particles a local approach is not only a theoretical requirement but also feasible.
△ Less
Submitted 13 March, 2015; v1 submitted 4 December, 2013;
originally announced December 2013.
-
Minkowski Tensors of Anisotropic Spatial Structure
Authors:
Gerd E. Schröder-Turk,
Walter Mickel,
Sebastian C. Kapfer,
Fabian M. Schaller,
Boris Breidenbach,
Daniel Hug,
Klaus Mecke
Abstract:
This article describes the theoretical foundation of and explicit algorithms for a novel approach to morphology and anisotropy analysis of complex spatial structure using tensor-valued Minkowski functionals, the so-called Minkowski tensors. Minkowski tensors are generalisations of the well-known scalar Minkowski functionals and are explicitly sensitive to anisotropic aspects of morphology, relevan…
▽ More
This article describes the theoretical foundation of and explicit algorithms for a novel approach to morphology and anisotropy analysis of complex spatial structure using tensor-valued Minkowski functionals, the so-called Minkowski tensors. Minkowski tensors are generalisations of the well-known scalar Minkowski functionals and are explicitly sensitive to anisotropic aspects of morphology, relevant for example for elastic moduli or permeability of microstructured materials. Here we derive explicit linear-time algorithms to compute these tensorial measures for three-dimensional shapes. These apply to representations of any object that can be represented by a triangulation of its bounding surface; their application is illustrated for the polyhedral Voronoi cellular complexes of jammed sphere configurations, and for triangulations of a biopolymer fibre network obtained by confocal microscopy. The article further bridges the substantial notational and conceptual gap between the different but equivalent approaches to scalar or tensorial Minkowski functionals in mathematics and in physics, hence making the mathematical measure theoretic method more readily accessible for future application in the physical sciences.
△ Less
Submitted 6 August, 2013; v1 submitted 13 September, 2010;
originally announced September 2010.
-
Local Anisotropy of Fluids using Minkowski Tensors
Authors:
Sebastian C. Kapfer,
Walter Mickel,
Fabian M. Schaller,
Markus Spanner,
Christian Goll,
Tomoaki Nogawa,
Nobuyasu Ito,
Klaus Mecke,
Gerd E. Schröder-Turk
Abstract:
Statistics of the free volume available to individual particles have previously been studied for simple and complex fluids, granular matter, amorphous solids, and structural glasses. Minkowski tensors provide a set of shape measures that are based on strong mathematical theorems and easily computed for polygonal and polyhedral bodies such as free volume cells (Voronoi cells). They characterize the…
▽ More
Statistics of the free volume available to individual particles have previously been studied for simple and complex fluids, granular matter, amorphous solids, and structural glasses. Minkowski tensors provide a set of shape measures that are based on strong mathematical theorems and easily computed for polygonal and polyhedral bodies such as free volume cells (Voronoi cells). They characterize the local structure beyond the two-point correlation function and are suitable to define indices $0\leq β_ν^{a,b}\leq 1$ of local anisotropy. Here, we analyze the statistics of Minkowski tensors for configurations of simple liquid models, including the ideal gas (Poisson point process), the hard disks and hard spheres ensemble, and the Lennard-Jones fluid. We show that Minkowski tensors provide a robust characterization of local anisotropy, which ranges from $β_ν^{a,b}\approx 0.3$ for vapor phases to $β_ν^{a,b}\to 1$ for ordered solids. We find that for fluids, local anisotropy decreases monotonously with increasing free volume and randomness of particle positions. Furthermore, the local anisotropy indices $β_ν^{a,b}$ are sensitive to structural transitions in these simple fluids, as has been previously shown in granular systems for the transition from loose to jammed bead packs.
△ Less
Submitted 15 October, 2010; v1 submitted 3 September, 2010;
originally announced September 2010.
