Hi All,
I'm a new user for liggghts. I've outputted the torque using from my simulation of a shear cell using the following command:
compute tor all wall/gran/local torque
dump 8 all local 100 tor.txt index c_tor[1] c_tor[2] c_tor[3]
My question is how would I process the dumped file for all the time steps?
Part of the dumped file:
ITEM: TIMESTEP
4200
ITEM: NUMBER OF ENTRIES
18
ITEM: BOX BOUNDS pp pp pp
-0.1 0.1
-0.1 0.1
-0.1 0.1
ITEM: ENTRIES index c_tor[1] c_tor[2] c_tor[3]
1 3.81974e-08 -6.52956e-08 0
2 6.51987e-08 8.83472e-08 0
3 1.068e-07 -4.40003e-08 0
4 1.46978e-09 -1.79701e-08 -1.06896e-07
5 -1.24187e-07 6.10127e-08 0
6 -8.67085e-09 4.32198e-09 0
7 -2.45119e-09 -2.44386e-09 0
8 -4.93194e-08 -2.75871e-08 1.20003e-07
9 -2.09081e-06 -5.67872e-07 8.84947e-07
10 3.20943e-06 1.46568e-06 -8.07744e-08
11 -1.78395e-10 -2.77594e-10 -1.37877e-09
12 3.19757e-08 -9.59372e-08 0
13 1.50163e-09 2.19332e-09 0
14 -8.69585e-08 -2.9615e-07 1.45874e-07
15 -2.52398e-11 4.61804e-10 0
16 5.5543e-08 -1.8916e-07 9.6018e-08
17 4.4508e-07 -2.73369e-07 1.18427e-08
18 -3.55666e-10 8.71359e-09 3.75247e-10
| Attachment | Size |
|---|---|
 tmp.jpg | 162.05 KB |
medvedeg | Wed, 05/17/2017 - 10:51
What do you mean by "post-process" torques?
Hello, sorry for late reply.
Torque is defined for every particle at every time step. You want to visualize particles in Paraview in color them with respect to the torue components/magnitude? Or plot a graph as a function of time? In the second case you can write an octave script that will iterate all the dump files and read them. Please, explain more what you want. Then I can help you.