packing2porenet 0.1.2

Description:

packing2porenet 0.1.2

packing2porenet
Python package to compute pore network based on given particle packing.
Overview
Computes pores and throats and their properties based on given packing.
spherical pores from spherical packing
Firstly, the weighted (also denoted as regular) Voronoi tesselation is done based on the input spherical packing. Weights are spheres' radii.
From the weighted Voronoi tesselation, Delaunay tetrahedralization is determined, resulting in the set of tetrahedrons and their connectivity.
Pores are assigned for each Delaunay tetrahedron.
Size and position of the pore is assigned as an inscribed spheres between the four "vertex" spheres.
Throats between pores are assigned radius equal to circle in face plane and inscribed between three "vertex" spheres.
It is possible to merge overlapping pores.
Then the new pore is such that:

volume is equal to the volume of spheres union
position is weighted average of positions, with weights being input volumes
Throats related to merged pores has radius such that the area of new throat is equal to the area of both merged throats.

Usage
Example
from packing2porenet import Sphere, Voronoi, PoreNetwork, export, ymport

# testing sphere pack
spheres = (
Sphere((0.00, 0.00, 0.00), 1.00),
Sphere((2.10, 0.10, 0.05), 1.20),
Sphere((4.05, 0.20, 0.10), 1.10),
Sphere((1.15, 2.15, 0.15), 1.15),
Sphere((3.05, 2.05, 0.05), 1.05),
Sphere((2.10, 4.20, 0.10), 1.25),
Sphere((1.05, 1.05, 2.10), 1.00),
Sphere((3.15, 1.10, 2.20), 1.10),
Sphere((2.05, 3.10, 2.05), 1.20),
Sphere((2.10, 2.15, 4.10), 1.15),
)

# pore network
export.spheres2vtk(spheres, "vtk-packing")
voro = Voronoi(spheres)
export.network2vtk(voro, "vtk-voro")
network = PoreNetwork().fromVoronoi(voro)
export.spheres2vtk([pore.asSphere() for pore in network.cells], "vtk-pores")
export.network2vtk(network, "vtk-throats")
network.mergeIntersectingPores()
export.spheres2vtk([pore.asSphere() for pore in network.cells], "vtk-pores-merged")
export.network2vtk(network, "vtk-throats-merged")

Compatibility
Tested on Ubuntu 18.04 LTS and Python 3.6
Installation


Using pip
python3 -m pip install packing2porenet


Using make (calls python3 -m pip install ., i.e. option 3, internally):
make install


Using pip locally
python3 -m pip install .


Contribution
Merge Requests
Are welcome
Bug reporting
In case of any question or problem, please leave an issue at the GitLab page of the project.
Contributors

Jan Stránský

License
This project is licensed under the MIT License - see the license file for details.
TODO

documentation
tests
periodicity
other shapes
other properties and functionality