Changes of ¨Cohesive Energy Density¨ with Temperature during the simulation

Dear everyone,

I am trying to model the litter bed drying process in a poultry house. So I would like to dynamically link the temperature of litter bed into the cohesive mode and update the Cohesive Energy Density based on the average temperature of all particles during the CFDEM simulation. Hence, I have added the following lines into my code and I am able to determine the average particle temperatures during the simulation.

But I am not able to replace the updated ¨cohesionEnergyDensity¨ during the simulation. Also, some additional descriptions are written inside the code!

1) Does anyone have a recommendation or a practical suggestion to solve this issue?
2) Is it possible to dynamically update the values of property/global during the cfdem simulation after a certain number of iterations?

Thanks for your help.
Sayed

fix integr all nve/sphere

compute 2 all reduce ave f_Temp[0]
fix 2 all ave/time 1 1 1 c_2 mode scalar ave running overwrite file ../DEM/averagetemp.txt

¨in the averagetemp.txt file the right average temp is written¨

# output settings, include total thermal energy
compute rke all erotate/sphere
thermo_style custom step atoms ke c_rke vol f_heattransfer cpu c_2 f_2
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

variable ced file ../DEM/averagetemp.txt
¨The value on the first column is read by liggghts while the value in the second column is required !!!¨

print "Average particles tempreture is equal to :${ced} centigrade."
¨liggghts prints 0¨

fix m6 all property/global cohesionEnergyDensity peratomtypepair 1 ${ced} # restartvalues yes

dump dmp all custom/vtk 1000 ../DEM/post/multi_*.vtk id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius density f_Temp[0] f_heatFlux[0]

restart 100000 ../DEM/restart/liggghts*.restart

run 0

print "Average particles tempreture is equal to :${ced} centigrade."
¨liggghts prints a right value¨