Hi,
actually I tried to figure out in the LAMMPS manual, how I could access the kinetic energy of my granular system.
I understood, that LAMMPS is calculating the kinetic energy in a MD sence (ke = 1/2 Kb T).
I wonder if the ke in my thermo-output is calculated like 1/2mv^2 for all particles if I use the LIGGGHTS package?
Thanks in advance,
Sebastian
ckloss | Mon, 06/06/2011 - 16:21
Hi Sebastian, the ke thermo
Hi Sebastian,
you are right, but nevertheless the ke thermo keyword indeed calculates the correct kinetic energy, because 3/2 kb T _is_ equal to the kinetic energy. You can check this by
compute 2 all ke/atom
compute 3 all reduce sum c_2
thermo_style custom step atoms ke c_3
(compute ke/atom uses m v^2/2 - you can check that in the source code)
The only thing that might be confusing about this is that the temperature used for this calculation does not have a meaning in a macroscopic sense, but "only" in an MD sense
Christoph
Silias | Mon, 06/06/2011 - 16:56
Hi Christoph, thanks for your
Hi Christoph,
thanks for your fast reply!
ok, wow that's surprising...
just for better understanding:
if both formulas calculate the same kinetic energy then, with
1/2 · m_particle · v_particle^2 = 3/2 · Kb · T
the temperature for 1 particel with mass=1kg and velocity=1m/s the Temperature of that particle would be
2.41 · 10^22 Kelvin (with Kb = 1,3806504 · 10−23 J/K).
So, if I would do a MD simulation this hypothetical temperature would be the temperature of my atoms, but as I am doing simulations with granular materials I can forget about that value!?
Thanks
Sebastian
ckloss | Mon, 06/06/2011 - 17:58
>>I can forget about that
>>I can forget about that value!?
Well, depending on the circumstances, this value will make sense (in the sense of statistical mechanics) or not.
For one particles it won't - for other situations it will. You may want to google "kinetic theory of granular gases" for details
Christoph