Greetings. I hope, you are all doing good. I'm working on homogeneous simple shear simulation using the following two methods.
In first method, I'm taking a cube of side 7.5 and filling it with 483 particles each of diameter 1 and density 1000 so that packing / volume fraction is 0.6. There is only 1 atom_type and all six faces are periodic. Particle fiction is 0.1 and there is no gravity for time being. I apply fix deform xy with a known value of dimensioned shear rate, say 1/sec, so that the top plane (ymax plane) moves along positive x-direction with a velocity of 7.5. Script writes out VTK files and I observe that the simulation box / region tilts correctly and at steady state, using spatial averaging I compute x - velocity gradient which closely matches with applied shear rate. This I have tried with lower as well as higher shear rates and also for different particle friction values and it works well.
In second method, I want to mimic same phenomenon by moving top and bottom walls in opposite directions (and hence not tilting box and not using fix deform). Rest all settings are same. I first used fix wall/gran on top (ymax) and bottom (ymin) planes to create walls of same atom_type and used shear keyword along with +3.75 and -3.75 units as velocity, respectively. The particles, close to the wall move very very slowly with velocity several orders of magnitude less. This could be due to the less friction between the wall and the particles. In next attempt, I created walls of particles of same diameter, density in a staggered arrangement such that the particles just touch each other. I'm using fix move on these two walls and I observe that a bit more momentum is transferred as compared to fix wall/gran but falling short by 40%.
1. If I use only fix wall/gran and declare walls as type 2, then what should be the wall-particle friction value? Perhaps, LIGGGHTS can take any value and can this value be as high as possible?
2. In particle wall method, I'm not integrating the wall particles using nve/sphere as I'm already moving them with fixed velocity. How can I improve this setup to achieve the correct shear rate?
Please let me know if you have any additional questions. I can share the scripts if you want more details.