Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

•The gas temperature rise is non-linear in the particle thermal mass loading ratio.•Preferentially concentrated particles produce large temperature fluctuations.•Hot spots in the boundary layer reflect streamwise elongated particle structures.•Particle concentration and gas temperature decorrelate at small length scales.•The Lagrangian history of particle clusters and turbulent mixing are important.