Hello,
I would like to collect all the dissipated collision energies from my simulation case (and possibly split into different modes normal and tangential). May I know how to do this? I mean which fix or compute and fix combinations to use?
I see that there are following lines defined in pair_gran.h, which I assume are for this purpose:
// stuff for tracking energy
int energytrack_enable;
class FixPropertyAtom* fppaCPEn; //collision potential energy normal
class FixPropertyAtom* fppaCDEn; //collision dissipation energy normal
class FixPropertyAtom* fppaCPEt; //collision potential energy tang
class FixPropertyAtom* fppaCDEVt; //collision dissipation energy viscous tang
class FixPropertyAtom* fppaCDEFt; //collision dissipation energy friction tang
class FixPropertyAtom* fppaCTFW; //collision tangential force work
class FixPropertyAtom* fppaDEH; //dissipation energy of history term (viscous and friction, accumulated over time)
double *CPEn, *CDEn, *CPEt, *CDEVt, *CDEFt, *CTFW, *DEH;
Regards,