Spurious horizontal movement with mesh walls

Submitted by tapsab on Wed, 09/12/2012 - 04:27

Hi,
Sorry to flood the bug report forum. Here's a synopsis of a bug that i recently encountered that causes spurious motion of particles that interact with mesh walls. The simulation (which can be recreated using the attached files) involves just one particle falling under gravity. The particles falls down vertically on a bottom wall. I tried three different ways to model the bottom and here is a summary of what I observe in each case:

1. Bottom wall modeled using the "primitive" option with "fix wall/gran/...": This case works as it should. The particle bounces a bit and settles down.

2. Bottom wall discretized using a coarse two-element mesh (use file mesh2.stl): As soon as the particle hits the bottom mesh it gains a spurious velocity in the x and z direction and starts to move diagonally. The forces on the particle are negligible (~1E-15), however the velocity in the x and z direction is around 1 mm/s and stays constant after the first few thousand iterations. This particular behavior is processor independent.

3. Bottom wall discretized using a fine mesh (use file mesh1.stl): The particle hits the bottom mesh it gains a spurious velocity in the x and z direction and starts to move diagonally. However unlike the previous case it oscillates back and forth as if bounded by some invisible wall. This behavior is processor dependent. I see it if I run it on 4 processors but not if I run it on a single processor.

Thanks,
Tapan

AttachmentSize
Binary Data Horizontal_Movement_Bug.tar_.gz4.02 KB

cstoltz | Wed, 09/12/2012 - 22:00

Very odd. I downloaded your deck and ran it and see similar problems to what you describe, though nothing leaps out at me. I would recommend increasing your z-direction bounds as you effectively have the simulation tight to the particle diameter, but even after doing that, I still see problems.

I've run many of the same sort of cases using a similar deck, but w/ SI units, and they run fine. I don't know if it is something related to the LJ units system, or possibly a roundoff error that is becoming a problem, or what.

Regards,
Chris

tapsab | Wed, 09/12/2012 - 22:30

Chris,
Thanks for taking the time to run the simulation. I think rather that it being a problem of units or roundoff the bug is somehow connected to how the mesh walls are handled because if you run the same simulation without the mesh walls using instead a "primitive" wall with the exact same properties none of the problems come into picture.

cstoltz | Thu, 09/13/2012 - 04:46

Could be. The examples I ran at my end also compared primitive walls, mesh walls, and dropping balls on the edges/corners of boxes or pyramid shaped structures. I haven't run into any issues as yet with version 2.1.1 on any of these simple examples.

That said, I'm still having some issues with some other problems in my personal test harness where particles are leaking out of some different types of batch mixers. It is a rare occurrence (which makes debugging that much harder), but is somehow related to the mesh as well I assume. I haven't tried really pushing the mesh discretization on my simple problems to see if I can prompt a failure that way. Will try it out and repost if anything comes of it.

Regards,
Chris

cstoltz | Thu, 09/13/2012 - 15:03

I repeated my runs with more finely discretized meshes and was able to reproduce a similar problem with the particles suddenly bouncing laterally. This happened if I dropped onto a node or slightly offset from the node.

Regards,
Chris

ckloss's picture

ckloss | Sun, 09/16/2012 - 23:51

Hi Tapan and Chris,

Thanks for reporting this! I can reproduce the issue. . We'll try our best to make this the last iteration cycle on meshes in 2.X

Thanks, Christoph