Skip to content

Commit 6c3290f

Browse files
committed
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';
1 parent 9bc6333 commit 6c3290f

File tree

1 file changed

+1
-2
lines changed

1 file changed

+1
-2
lines changed

tutorials/cfdemSolverPiso/settlingTestMPI/DEM/in.liggghts_run

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -24,7 +24,6 @@ fix m1 all property/global youngsModulus peratomtype 5.e6
2424
fix m2 all property/global poissonsRatio peratomtype 0.45
2525
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
2626
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
27-
#fix m5 all property/global characteristicVelocity scalar 2.0
2827

2928
# pair style
3029
pair_style gran model hertz tangential history
@@ -64,7 +63,7 @@ compute rke all erotate/sphere
6463
thermo_style custom step atoms ke c_rke vol
6564
thermo 1000
6665
thermo_modify lost ignore norm no
67-
compute_modify thermo_temp dynamic yes
66+
compute_modify thermo_temp dynamic no extra 0
6867

6968
# insert the first particles so that dump is not empty
7069
run 1

0 commit comments

Comments
 (0)