Two-way coupled Cloud-In-Cell modeling of non-isothermal particle-laden flows: a Subgrid Particle-Averaged Reynolds Stress-Equivalent (SPARSE) formulation