I am really appreciate for your reply. However, there are still some unclear places.
1.My task is to verify the term of collisional dissipation of energy (the source term of granular temperature) deduced by ourselves. In the start, I use the PDE of granular temperature Model, the interface of granular temperature source can be found in Cell Zone Conditions（As shown in the picture）, but at this time, the calculation is always divergent. So I want to use the Phase property of granular temperature Model to have a try, however, the interface of granular temperature source in Cell Zone Conditions disappear ,and then I don't know where should I add the collisional dissipation of energy (the source term of granular temperature) deduced by ourselves.
2.I need to set initial stacking area in adapt, maybe the mesh of here cause errors. It means if i set initial stacking area in adapt and then i won't get the gradient of granular temperature by the above method? Is there any other way to obtain the gradient of granular temperature？