Problem with fix deform

LAMMPS LAMMPS Mailing List Mirror
1

Dear lammps users,

Now I’m using fix deform and Sllod to simulate a shear flow.
Using the code shown below, I can get a linear velocity profile and the velocity of the fluid ranges from 0 to 10, for instance.
Fig.1 is the velocity profile I can obtain now.
My question is how can I get a velocity profile in which the velocity ranges from -5 to 5?
Fig.2 is the velocity profile I want to obtain.

My code is like:

fix 2 all nvt/sllod temp 300.0 300.0 0.1
fix 3 ar deform 1 xy erate 0.1 remap v

I’m looking forward to having your help.

Best wishes
Ruo-Yu Dong

Fig.1.jpg

Fig.2.jpg

2

My question is how can I get a velocity profile in which the velocity ranges
from -5 to 5?

You can't. LAMMPS does the deformation asymmetrically
by moving the top box surface in one direction.

But what does it matter? 0 to 10 versus -5 to 5 are equivalent.

Steve