I am moduling a case of reactive flow which involves a specie created by myself. When the case converged, I found that all variables' distribution were normal except the temperature distribution. It seems the energy created by the chemical reaction couldn't be added to the flowfield. I have checked each specie's material property(especially the Cp) and they are all correct. Can anyone give some advice where I should check for my case? Thanks!

Run your simulation with and without the reaction and see if there is any change in the energy contours.