Time dependent Temperature

Submitted by Mateus Kostka on Fri, 03/26/2021 - 09:23

Dear Ladies and Gentlemen,

is it possible to assign in Liggghts time dependent temperature profiles to the particles?

My plan ist to calculate the thermal expansion of my particles in dependence ofhe thermal expansion coefficent and the timedependent temperature.
So far I assign the temperature in that way:

fix 1 all heat/g
ran initial_temperature 373.15 contact_area overlap

Is it possible to print the temeprature of each particle in a text file?
Maybe in that way?

dump dmp all custom 5000 dumpSinternE*.liggghts id type x y z vx vy vz f_Temp[0] radius

I tried to implement a heat source, but it had no influence on the printed Temperature f_Temp[0]:

fix 1 all heat/gran initial_temperature 273.15 heatSource 500 contact_area overlap

Best regards and thanks in advance

Mateus Kostka