Touching neighbor mesh elements and unexpected behaviour

msbentley's picture
Submitted by msbentley on Thu, 08/02/2012 - 15:55

Hi all,

I have a simulation comprising 250k particles of ~100 µm diameter packed loosely between two concentric cylinders. The packing is read from a file with read_data, and I simply want to let the packing settle under gravity for now (gravity is along the axis of the cylinder). I have two mesh walls defining the inner and outer cylinders, and as in the attached images no particles are in contact with the walls.

If I don't import the inner mesh, the particles collapse fine and residual energy and motion damps out quickly. However, with the inner mesh in place I get warnings such as:

INFO: more than 3 touching neighbor mesh elements found, growing contact history

and shortly after this the typical "time step too large or skin too small" warnings from check/timestep/grain appear.

The simulation seems to keep getting "kicks" of energy, and the near-surface particles don't quickly damp down.

I have tried with different mesh element sizes, but to no avail (see attachment). Can anyone explain what the mesh warning above means? Is it related to the mesh itself, or particle-mesh interactions?

Thanks! Mark

AttachmentSize
Image icon cylinder_setup.png705.56 KB
Image icon cylinder_setup2.png727.95 KB
Image icon cylinder_setup3.png53.83 KB
Image icon all_meshes.png33.32 KB
ckloss's picture

ckloss | Sun, 08/05/2012 - 14:15

Hi Mark,

the info message

INFO: more than 3 touching neighbor mesh elements found, growing contact history

means that some particle is touching more than 3 mesh elements at a time. Per default, memory is allocated for 3 element contacts per particle. If more contacts are detected, additional memory for particle-mesh element contact history is allocated (if contact history is used).

So this just informs you about the memory re-allocation. Nothing to worry about. I just put this message so the user can evaluate if the number of particle-element contacts are reasonable (e.g. if the wrong mesh scaling factor was chosen or alike)

Cheers, Christoph