Dear LAMMPS users,
I am using the latest version of LAMMPS. I’d like to calculate viscosity using NEMD method. The problem is when I am trying to calculate the viscosity of a polymeric liquid under different shear rates, I see shear thinning while I know that this polymer should has a Newtonian behavior(polymer = PEG400).
(shear thinning: viscosity decreases by increasing shear rate)
(Newtonian behavior: viscosity is constant under different shear rates)
And about my velocity profile; despite I use the real units, Is it ok to haven’t desired smooth velocity profile because of very low number of shear rates?
I have sent the input script and data file. I would appreciate if any one can give me some advice about my mistakes.
thanks in advance,
Eliar
mycell30.data (580 KB)
run.in (2.19 KB)