LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Reflecting surfaces momentum conservation
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Reflecting surfaces momentum conservation

From: "Vaidyanathan M.S" <msvaidyanathan1729@...24...>
Date: Tue, 18 Jul 2017 15:55:10 -0500


I am trying to simulate polymers between two surfaces. For a test case, I simulated a few polymer particles between two reflecting surfaces (fix wall/lj93). Polymers are homopolymers with LJ-FENE interactions.

Q: Even when I start with zero momentum initial condition, when I check the result after some timesteps (ran using Nose Hoover thermostat - which is supposed to have the momentum conserving property) the momentum in the direction perpendicular to the surface is not zero. I could guess that the momentum in the direction perpendicular is not zero since there is an unbalanced force from the surface. Question is whether fix momentum will be a good fix (or it will just cure the symptoms and not the problem?) for this or there is something else that needs to be taken care of LAMMPS which is equivalent to the force holding the surface in place?

PS: If required I can send a small datafile and an input file to generate the scenario. The datafile I am having is big to load here. Sorry about that


Vaidyanathan M S
PhD Candidate
University of Texas at Austin