Adding Temperature to icoUncoupledKinematicParcelFoam
 

Hi!

I am trying to create a model that has a TEqn.H with the following equation:

ddt(T)+U. div(T) = laplacian(alpha.T) +heat transfer(particle-particle) +heat transfer(wall-particle)

But, the last two terms would require particle temperatures at each time step (depends on T(1st particle)-T(2nd particle)). How do I get these values? and what would the initialization look like?

Thank you!!