Miracle of PBC and Meshes

richti83's picture
Submitted by richti83 on Fri, 04/10/2015 - 14:21

I needed to simulate a simple troughed chain conveyor, so I made a box from -1 to 1 and set PBC in conveying direction.
After the first plate goes through the PBC it starts to flicker, this means it is one time at the end and the other time at the beginning of the simulation box.

Furthermore simulation time increases by a factor of ten after that step (in the example it happens between step 44.000 and 45.000).
BUT (and now the magic happens) when I start a fresh run after let's say 65.000 iterations for some reasons, the simulation time becomes normal until the next mesh reaches the boarder.
3rd wonder: It takes more time when I devide the mesh in 8 plates á 2 triangles instead of having all plates in one mesh (1x16 triangles)

The timing tells me, that it is the modify time which significantly increases. I tried to use a profiler application (zoom in this case) but it only pointed me to Neighbour::bin_atoms() line 1980.

simple testcase attached (first 65.000 iterations bad, 2nd 65.000 iterations perfect, last 65.000 iterations horrible performance)

Tested with LIGGGHTS 3.1.0 git commit b3030a8fd1f722c77f69150c7b2d24b386eab848

Have a Nice Weekend,
Christian.

AttachmentSize
Binary Data pbc_test.tar_.gz1.86 KB
ckloss's picture

ckloss | Fri, 04/17/2015 - 09:19

Hi Christian,

thanks for posting this. Actually, PBC and movng meshes is actually "not meant to be". I wonder why I did not add an error message to prohibit that combination in the first place. I'll need to wrap my head around it before the next release, see what the restrictions are, and will then comment on it here what would be needed to get it going

Cheers
Christoph

ckloss's picture

ckloss | Tue, 04/28/2015 - 20:11

Hi Christian,

since this release was quite busy, I did not have a chance to follow up with this. I'll put it on the agenda for the next release

best wishes
Christoph

ckloss's picture

ckloss | Fri, 04/22/2016 - 15:56

Hi Christian,

the next release will include an update for this! Feel free to check if this now improves the situation for your case!

best wishes
Christoph