Dear All,
I try to simulate the Couette flow using LAMMPS. From one input script named in.flow.couette which is localed in the lammps*/example/flow sub-directory, we can see
fix 1 all nve
fix 2 flow temp/rescale 200 1.0 1.0 0.02 1.0
fix_modify 2 temp mobile
velocity lower set 0.0 0.0 0.0
velocity upper set 3.0 0.0 0.0
The linear velocity profile of the fluid can be obtained because of the different velocities of the two walls and the interaction between the walls and the fluid.
However, the velocity of the fluid is still affectted by the temperature controlling command—fix temp/rescale.
My question is:
How the fluid maintain the linear velocity profile under the influence of both the solid-fluid interraction and the fix temp/rescale command?
I try to simulate the couette flow of water. The input script is just similar to in.flow.couette. But I can’t get the linear velocity profile of the fluid unless I turn off the fix temp/rescale command. Since I deleted the temperature controlling command, I found that the fluid is heated and the temperature is rising all the time.
Any help will be highly appreciated.
Thank you all very much
2009-02-20