cohesionEnergyDensity

Submitted by paulbrumby on Tue, 09/25/2012 - 03:20

I just had a quick question. I hope it's okay. I am trying to do a comparison of the various forces present within my granular simulations and want to plot these in an excel spreadsheet.

I am using the gran/hertz/history pair_style but I am having problems seeing how the cohesionEnergyDensity is used. I cannot find mention of it in the manual so if anyone could point me in the right direction, or better yet let me know the equation that it appears in that would be really helpful for me!

Thanks!

Paul

Forums:

moritzhoefert | Wed, 09/26/2012 - 08:26

Hi Paul,

check doc/pair_gran.html for information about the cohesion energy density model equation.

Cheers,
Moritz

paulbrumby | Mon, 10/01/2012 - 09:56

Dear Moritz,

Many thanks for your reply. I have taken a look through the manual and the page you suggested, as well as the paper of Brilliantov et al. (1996), but I have to confess I am still at a loss to see where and how the cohesionEnergyDensity variable is used within the code. If you could point me in the right direction I'd really appreciate it!

Best regards,

Paul

ckloss's picture

ckloss | Mon, 10/01/2012 - 12:11

Have a look at pair_gran_hooke_history.cpp!

Cheers, Christoph

paulbrumby | Tue, 10/02/2012 - 03:40

Thanks Christoph,

I think I see what is going on now! I expect that my choice of contact model would use pair_gran_hertz_history.cpp but I see that from within there PairGranHookeHistory is somehow linked. So that's how cohEnergyDens is used with the hertz model?

I'm not so good at reading cpp but I think I can get what I need from here. Thanks again!

Paul