pxx, pyy, pzz, pxz, pyz, pxy in thermo_style

Submitted by Weijing on Wed, 12/07/2016 - 23:44

Hi everyone,

I am a new user to liggghts and applying the packing problem in the examples to practise. In a case, I want to know the equivalent, or saying apparent, pressure of the packing system in x, y, z direction. I find that in the thermo_style command the pij may give the quantity. however, in the source code, it doesn't have those option actually, resulting in a error when I add them as the keyword.

Can anyone suggest me how to get the equivalent pressure or force for a packing in xyz directions.

Thanks
Weijing

ckloss's picture

ckloss | Mon, 12/12/2016 - 21:59

Hi Weijing,

you can use the fix ave/euler command - this will give you the stress tensor of the packing!

Christoph