Output-Temperature of particles in a certain region (LIGGGHTS)

Submitted by Volkan on Mon, 07/04/2022 - 17:34

Hello everyone,

i wanted to write the temperatures of the particles in a certain region into a text or csv file. Getting the mass, the number or other many properties of the particles in a certain region is possible easily with region functions, however i could not find any region function with temperature. For example for mass i wrote this and it works.

variable m1 equal mass(all,outlet)
variable t equal step*dt
fix data all print 100000 "$t,${m1}" file ../DEM/Messdaten/mass.csv title "Zeit(s),Masse_in(kg)"

I need the same thing but for the temperature.

Normally i can get the dump files with temperature of particles with f_Temp[0] syntex but it's from the all particles everywhere. i want to get the temperature only of those particles, which are in the region that is defined as my outlet.

(It's about an extruder screw. I want to see the temperatures of the particles at the end of the screw, so that i can see the difference clearly)

Is that somehow possible? I would really appreciate if you guys can help me out :)

Thanks in advance.

Best regards,
Volkan Atar

mschramm | Thu, 07/14/2022 - 00:52

Hello,
Please look into compute/reduce/region command. As an example using heatTransfer_1 and first moving the command {region halfbed block 0 INF INF INF INF INF units box} higher into the script.

# Create the needed compute
compute aveTemp all reduce/region halfbed ave f_Temp[0]

# print the average temperature in the halfbed region to screen
thermo_style custom step atoms ke c_rke f_heattransfer vol c_aveTemp