LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] how to calculate stress tensor of a bin

# Re: [lammps-users] how to calculate stress tensor of a bin

 From: Eric Murphy Date: Mon, 2 Oct 2017 13:54:44 -0500

#2. You don't have to care about voronoi volumes at all if you don't want the stress per atom - it's volume is irrelevant in your case.  If you use #1, you will be adding samples and volumes in which they live/occur haphazardly - and it won't really represent anything at all.

Nonetheless, binning and using the bin volume can present the exact same problems on your averages as using the voronoi volumes - scale dependence from subdividing your domain. Should be less severe, however.

Regards,
Eric Murphy
Iowa State University

On Fri, Sep 29, 2017 at 11:18 AM, Yan, Rui wrote:

Dear lammps users,

I divide my simulation box into small layers along z direction using compute chunk/atom command.

I want to calculate the stress tensor of each bin. I don’t know which of the following ways I should take:

1.     I know that I can use compute stress/atom command to get the stress*volume of each atom. Here volume is the volume of an atom. Just as said in lammps manual, an individual atom's volume is not well defined. One choice is to use compute voronoi/atom command for one possible way to estimate a per-atom volume. With results from compute stress/atom command and compute voronoi/atom command, I can get the stress tensor of each atom and then sum them up within a bin to get the stress tensor of a bin.

2.     The second way is to use compute stress/atom command to get the stress*volume of each atom, then sum up the results of compute stress/atom command within a bin. Finally the sum  is divided by the volume of a bin to get the stress of a bin.

Because the sum of each atom’s voronoi volume in a bin may not equal the volume of the bin, so these two ways may give different results. Could you tell me which way I should use? Or they are both wrong and there is another way to do it?

Rui

I know that I can easily get the potential energy of each atom and sum them together within a bin to get the potential energy of a bin, but I don't know how to calculate the temperature and stress of a bin. Can you give me some suggestions?

Thanks.

Rui

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most