LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Fwd: Thermal/conductivity
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Fwd: Thermal/conductivity

From: Neda Rafiee <ne.rafiee@...444...>
Date: Fri, 18 Aug 2017 18:51:08 +0430

Actually yes, without swapping atoms, energy is conserved in my system. When I use thermal conductivity, kinetic energy of the system increases and total energy does not remain constant. I cannot find the problem!

On Thu, Aug 17, 2017 at 8:16 AM, Neda Rafiee via lammps-users <> wrote:
Dear users,
I am using fix Thermal/conductivity to measure thermal conductivity of a box of water using Muller-Plathe algorithm. Actually I am swapping one atom every 10 timesteps. when I draw temperature of the first and the middle layers, I see an interesting thing happening... Both layers are increasing in temperature with time.  As you know, I am working with an NVE ensemble and I expect it to conserve energy, but total kinetic energy is increasing in time. Anyone knows why this happens?

​does your simulation conserve energy well without swapping atoms?



Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list

Dr. Axel Kohlmeyer  akohlmey@...24...
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.