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

Re: [lammps-users] Reflecting surfaces momentum conservation

From: Stefan Paquay <stefanpaquay@...24...>
Date: Wed, 19 Jul 2017 04:11:34 -0400

You can calculate the momentum exchange between the wall and your atoms with compute and variable commands, since the fix calculates the forces the wall exerts. If the thermostat is momentum conserving then your setup should probably satisfy (change in momentum)  == (impulse from wall), at least,  up to numerical errors.  

On Jul 18, 2017 10:56 PM, "Vaidyanathan M.S" <msvaidyanathan1729@...24...> wrote:

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

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