Have a look at DEFINE_DPM_HEAT_MASS too. DPM's interaction with the fluid flow is complicated, in part because of the Lagragian-Eulerian converversion and in part because of the sheer number of models that interact with DPM and the fluid phase. You may find you need to be creative with how you define things.