Hi all,
I am trying to model a 3D couette flow with lj fluids (not the 2D couette flow as in one of LAMMPS examples). I am wondering if there is a way to rotate a cylindrical wall made of atoms. Thank you very much for your help.
Shihai
No - but either the fix wiggle or fix rigid commands could
be extended to do this with relatively little effort I think.
Steve
Actually, I think if you put all the cylinder wall atoms
in one rigid body with fix rigid, and assign velocities
to the atoms so that the body has a
net angular momentum with the rotation you want,
and net translation of 0.0, then if you use the force
and torque keywords of fix rigid to turn off external force/torque
on those atoms, then the cylinder will rotate as you desire.
Steve