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
License
For personal and professional use. You cannot resell or redistribute these repositories in their original state.
Customer Reviews
There are no reviews.
