sphere particle rotation

Submitted by tala on Tue, 11/13/2018 - 16:50

Dear all
I am wondering if there is any way to calculate sphere particle's cumulative rotation. is it possible to dump particle's rotation through run? if the answer is no, how it can be done by postprocessing?
any help is appreciated.

arnom's picture

arnom | Thu, 11/15/2018 - 15:30

It's a bit tricky to define one rotation as a particle can tumble and rotate back and forth around changing axes. So you need to have a good definition of what you actually want to measure. Can you clarify?

DCS team member & LIGGGHTS(R) core developer

tala | Tue, 11/20/2018 - 12:19

I mean for one particle all back and forth rotations be summed. and the resultant rotation be dumped in specified step.
Thank you for attention.
Talat

zumack | Sun, 11/25/2018 - 19:15

I think you could use the compute/erotate/sphere comand to calculate the kinetic energy for each atom. After that you could dump it with dump/custom. And after that you will perform a post processing with dump archive. Because the result that you have will be ke=1/2 I w^2, since that result it is posible to have the angular velocity and whit the numerical integration, the angular displacement at each time step. I wondering it is possible to posprocess in Matlab.

If this way is possible, please write me, because i will need to do something like that in the future with mi project.

Best regards

tala | Fri, 11/30/2018 - 11:49

Hi
To dump particles angular velocity there is a simple way. By dump custom omgx omgy omgz you could have particles angular velocity. But to integrate them over time you need to dump velocity in every timestep. By now I am trying to dump mean angular velocities, if it work for me, i will inform you.
Regards,