Dear All
Who can help me on this "fix viscosity"?
1) in the lammps example, the unit is lj. we also can use another unit in the system, right?
2) variable, dVx equal f_5[11][3] - f_5[1][3] Can you explain clearly on this equation?
3) for units reduced, does this relate to unit in the system or something else?
4) variable visc equal -(f_4/(2*(step*0.005-125)*lx+1.0e-10))/(v_dVx/(ly/2))
from lammps help website, 1.0e-10 is a fudge factor to avoid a divide-by-0.0 on the 1st time step, right?
ly/2 is the distance over which that delta occurs, so dividing by it gives a slope, right?
if 3 dimension here, what are the parameters?
Thanks for your time.
Look forward to your reply.
Example from lammps:
# sample LAMMPS input script for viscosity of 2d LJ liquid
# Muller-Plathe via fix viscosity
# settings
variable x equal 20
variable y equal 20
variable rho equal 0.6
variable t equal 1.0
variable rc equal 2.5
# problem setup
units lj
dimension 2
atom_style atomic
neigh_modify delay 0 every 1
lattice sq2 ${rho}
region simbox block 0 $x 0 $y -0.1 0.1
create_box 1 simbox
create_atoms 1 box
pair_style lj/cut ${rc}
pair_coeff * * 1 1
mass * 1.0
velocity all create $t 97287
fix 1 all nve
fix 2 all langevin $t $t 0.1 498094
fix 3 all enforce2d
# equilibration run
thermo 1000
run 5000
unfix 2
# turn on Muller-Plathe driving force and equilibrate some more
velocity all scale $t
fix 4 all viscosity 100 x y 20
fix 5 all ave/spatial 20 50 1000 y center 0.05 vx &
units reduced file profile.mp.2d
# equilibration run
variable dVx equal f_5[11][3]-f_5[1][3]
thermo 1000
thermo_style custom step temp epair etotal press f_4 v_dVx
run 20000
# data gathering run
# reset fix viscosity to zero flux accumulation
unfix 4
fix 4 all viscosity 100 x y 20
variable visc equal -(f_4/(2*(step*0.005-125)*lx+1.0e-10))/(v_dVx/(ly/2))
fix vave all ave/time 1000 1 1000 v_visc ave running start 26000
thermo_style custom step temp f_4 v_dVx v_visc f_vave
# only need to run for 10000 steps to make a good 100-frame movie
#dump 1 all custom 100 dump.mp.2d id type x y z vx
#dump 1 all image 100 image.*.jpg vx type zoom 1.6 adiam 1.2
#dump_modify 1 pad 5 amap -0.7 0.7 ca 0.0 2 min blue max red
run 50000