I see that several people in the forum has experienced particles penetrating through the wall, but has anyone experienced a case in which particles do/don't penetrate the wall depending on the computational domain size?
In my case (I cannot share my files unfortunately), I have an L-shaped pipe, and in one end of this pipe I have another geometry inside the pipe where I start inserting particles to settle inside. When I make the computational domain to surround just that inner structure, then there is no issue. Once I make the computational domain to include all of the structures, then I start seeing about half of the particles that are well within the computational domain start penetrating through the walls, while the other half stay settled on the surface of the inner structure. Below is an example of what I mean in terms of my input file:
atom_style granular
atom_modify map array sort 0 0
comm_style tiled
comm_modify mode single vel yes
stiffnessScaling 3.177672
boundary f f f
newton off
units cgs # was si
processors * * *
region reg block -1.6 1.6 -0.1 16.4 -11.8 1.6 units box
The above region is the computational domain that includes the L-shaped pipe, in which half of my particles penetrate through the walls.
If I change the computational domain to be tighter and to surround just the inner structure inside one end of the L-shaped pipe:
region reg block -1.6 1.6 10.0 16.4 -11.8 1.6 units box
then I no longer experience any penetrations. The diameter of the largest particle in my simulation is 70 microns, (i.e. 0.007 in the input file, which is in cgs). Any suggestions would be very helpful. Thank you.