Hi All,
I has a liggghts compression code and want to visulaize force chain in paraview. To visulaize force chain, I used these two command lines ("compute_pair_gran_local.cpp" and "dump_local_gran_vtk.cpp").
But I have a issue that force chain cross the simulation domain along x and y direction (see attached pictures). Boundary condition is periodical along x an y direction. There are two walls along z direction. Top wall move down along z direction to create compression.
Could anyone know the reason of this issue and help me solve this issue?
I am wondering if these source codes "compute_pair_gran_local.cpp" and "dump_local_gran_vtk.cpp" are able to avoid force lines that cross the simulation domain along the periodic boundary?
Thank you in advance!
Best wishes,
chrisqian
| Attachment | Size |
|---|---|
 particle_under_compression.png | 841.68 KB |
| force_chain_line.png | 362.44 KB |
| force_chain_width.png | 550.15 KB |
Salieri | Mon, 01/22/2024 - 09:48
Have you solve this problem？
Hi chrisqian,
More than half a year has passed since your question, do you have any solutions to this problem?
In my research, I also need to use periodic boundary conditions, and I also encountered the same problem as you. As for my guess, this phenomenon is caused by the fact that some aggregates are still at the end when a aggregate returns from the end to the initial end, but I do not know how to solve it. If you have a good solution now, would you please give me some suggestions.
Looking forward to your reply！
Sincerely,
Salieri