couette flow simulation result differ from the analytical result

Dear Sir,

    When i run couette flow examples with Lammps, i got the the average velocities of atoms along the direction normal to the flow through c code. When i plot the velocities, they show a relationship of curve lline, instead of straight line. Actually, the relation should be straight line according to the Couette Fomulla:

[u (y) = u_0\frac{y}{h}]

<http://upload.wikimedia.org/wikipedia/commons/9/93/Laminar_shear.svg> I do not know where is the problem althouth i try changing almost every index. So i really appreciate your reply!

    Best,

   Eric

Laminar_shear.bmp (257 KB)

Dear Sir,

When i run couette flow examples with Lammps, i got the the average velocities of atoms along the direction normal to the flow through c code. When i plot the velocities, they show a relationship of curve lline, instead of straight line. Actually, the relation should be straight line according to the Couette Fomulla:

[u (y) = u_0\frac{y}{h}]

<http://upload.wikimedia.org/wikipedia/commons/9/93/Laminar_shear.svg&gt; I do not know where is the problem althouth i try changing almost every index. So i really appreciate your reply!

most of the examples in the LAMMPS distribution
are supposed to show, how something works in
principle (i.e. what commands can be combined how),
serve as a simple test case and run fast.

they are *not* meant to give an example for how
to do something properly. you are expected to
read up on how to do such simulations and then
create meaningful inputs with more reasonable
choices of parameters instead.

cheers,
    axel.