From cf26690feb5cd0d4e35ad3ab2679d33f5b84d1cc Mon Sep 17 00:00:00 2001 From: athomps Date: Tue, 12 Apr 2016 14:56:22 +0000 Subject: [PATCH] Corrected formula for case lx != ly git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@14820 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- examples/VISCOSITY/in.nemd.2d | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/examples/VISCOSITY/in.nemd.2d b/examples/VISCOSITY/in.nemd.2d index e27fe83327..0ff9c1e321 100644 --- a/examples/VISCOSITY/in.nemd.2d +++ b/examples/VISCOSITY/in.nemd.2d @@ -10,7 +10,7 @@ variable rho equal 0.6 variable t equal 1.0 variable rc equal 2.5 -variable srate equal 2.7 +variable srate equal 2.7 # velocity of top edge # problem setup @@ -22,7 +22,7 @@ neigh_modify delay 0 every 1 # problem setup lattice sq2 ${rho} -region simbox prism 0 $y 0 $y -0.1 0.1 0.0 0.0 0.0 +region simbox prism 0 $x 0 $y -0.1 0.1 0.0 0.0 0.0 create_box 1 simbox create_atoms 1 box @@ -48,7 +48,9 @@ unfix 2 velocity all scale $t -variable xyrate equal ${srate}/lx +# shear rate defined relative to perpendicular dimension + +variable xyrate equal ${srate}/ly fix 1 all nvt/sllod temp $t $t 0.1 fix 2 all deform 1 xy erate ${xyrate} remap v