A Python package for geometry construction in particle-based simulation.
I mainly use this package for geometry construction in LAMMPS, with some examples provided in the repository; of course, it can also be used for other software.
pip install geoparticleDownload and enter the source code directory, then
pip install .
The documentation is hosted on Read the Docs.
Particles of specified geometries are typically created by the lattice command in LAMMPS, which can lead to rough surfaces when the particle spacing is not small enough. However, too small spacing can result in too many particles and thus increase the computational cost.
The case is the same when one creates atoms based on an external STL file (an example STL file exported by COMSOL is shown below):
To resolve this problem, I developed this package for easy construction of geometries where smooth surfaces are required. Miscellaneous geometries are provided, including 1D geometries (lines and curves):
2D geometries (rectangles and circles):
3D geometries (blocks, cylinders, tori, and spheres):
all of which can be surface, thick shells, or filled bodies.
Diverse operations are also provided, including translation, mirror, rotation, stack, clipping, union, intersection, and subtraction.
Some utility functions are also provided.
The resultant geometries can be exported to various formats, including LAMMPS data &restart & dump file, and SPHinXsys restart file.
examples/gallery.py provided detailed scripts to yield the geometries above.
Two more examples are given to shown how to couple geoparticle with LAMMPS. The first example is the 2D gas-liquid dam break, while the second is a 3D intestine connected with an inlet and an outlet.
See the documentation for details.
The particle spacings may be not exactly as specified in order to create a smooth surface.
Resultant geometries of boolean operations can have more particles than expected in some cases, because
- For intersection and subtraction, only particles with distances smaller than
rmaxwill be identified the same. Users should align particles of different geometries to get the expected results. - For union, particles of all the given geometries will be collected to yield the union. Users should ensure no particles are overlapped.