Temperature exponential

Riccardo Maione's picture
Submitted by Riccardo Maione on Tue, 04/01/2014 - 10:32

Hi all,

I would like to write a fix who store the function exp(1/T), any advice about what fix or compute I could take as a start?

Thanks all

Forums:

richti83's picture

richti83 | Thu, 04/03/2014 - 22:00

Here is a very simple fix which updates the velocity of all particles in a given region (ignores forces, but before position integration, nice to geneate a given output velocity of a hopper):

fix_set.cpp
http://pastebin.com/Y5778GiL

fix_set.h:
http://pastebin.com/R73Zu07V

the timestep is update->dt;
the temperature seems to be a property/atom introduced by fix_heat_gran.cpp, see code there how to access it (I never used heat transfer).

Best,
Christian.

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

ckloss's picture

ckloss | Fri, 04/04/2014 - 10:01

Hi Riccardo,

you could also do this with a variable command!

Christoph