fix output of kinetic energy in settling test
in case of a single particle, we need to specify the option 'extra 0'; since we just create one particle at the start of the simulation, we can also set 'dynamic' to 'no';
This commit is contained in:
@ -24,7 +24,6 @@ fix m1 all property/global youngsModulus peratomtype 5.e6
|
||||
fix m2 all property/global poissonsRatio peratomtype 0.45
|
||||
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
|
||||
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
|
||||
#fix m5 all property/global characteristicVelocity scalar 2.0
|
||||
|
||||
# pair style
|
||||
pair_style gran model hertz tangential history
|
||||
@ -64,7 +63,7 @@ compute rke all erotate/sphere
|
||||
thermo_style custom step atoms ke c_rke vol
|
||||
thermo 1000
|
||||
thermo_modify lost ignore norm no
|
||||
compute_modify thermo_temp dynamic yes
|
||||
compute_modify thermo_temp dynamic no extra 0
|
||||
|
||||
# insert the first particles so that dump is not empty
|
||||
run 1
|
||||
|
||||
Reference in New Issue
Block a user