Simple shear flow with lammps

LAMMPS LAMMPS Mailing List Mirror
1

Dear all,

I want to have a simple shear flow in my simulation box which means zero velocity in the middle of the simulation box.
I am using :

fix 1 all nvt/sllod temp 1.35 1.35 0.5 tchain 1
fix 2 all deform 1 xz erate 0.03 remap v

but this is producing a linear flow with the zero and the maximum velocity in the bottom and top of the simulation box respectively.

Is there any idea?

Best regards,
Atefeh

2

That's the way fix deform and fix nvt/sllod work. See their doc
pages. It is equivalent to a flow with 0 in the middle and +/- v/2
at the top and bottom.

Steve