Total Kinetic Energy

LAMMPS LAMMPS Mailing List Mirror
1

Hi,

I need to compute the total kinetic energy for a polydisperse system in a 2D Couette flow as it is being sheared over time in order to determine whether the system has reached steady state conditions. I plan to plot the total kinetic energy versus number of timesteps. Do I simply do this:

sumOver_i [ 0.5(mass_i x velMagnitude_i ^2) ]

Where:

velMagnitude_i = sqrt(vx_i^2 + vy_i^2)

For all particles i?

Or is there a simpler way?

Kind regards,

Andrew

2

Hi,

I need to compute the total kinetic energy for a polydisperse system in a
2D Couette flow as it is being sheared over time in order to determine
whether the system has reached steady state conditions. I plan to plot the
total kinetic energy versus number of timesteps. Do I simply do this:

sumOver_i [ 0.5(mass_i x velMagnitude_i ^2) ]

Where:

velMagnitude_i = sqrt(vx_i^2 + vy_i^2)

For all particles i?

Or is there a simpler way?

​or use the "ke" thermodynamic variable.

axel.​