This paper presents an Euler-Euler Large-Eddy Simulation (LES) approach for the numerical modeling of non isothermal dispersed turbulent two-phase flows. The proposed approach is presented and validated by a priori tests from an Euler-Lagrange database, provided using discrete particle simulation (DPS) of the particle phase coupled with direct numerical simulation (DNS) of the turbulent carrier flow, in a non isothermal particle-laden temporal jet configuration. A statistical approach, the Mesoscopic Eulerian Formalism (MEF) [Fe´vrier et al., J. Fluid Mech., 2005, vol. 533, pp. 1–46], is used to write local and instantaneous Eulerian equations for the dispersed phase and then, by spatial averaging, to derive the LES equations governing the filtered variables. In this work, the MEF approach is extended to scalar variables transported by the particles in order to develop LES for reactive turbulent dispersed two-phase flows with mass and heat turbulent transport. This approach leads to separate the instantaneous particle temperature distribution in a Mesoscopic Eulerian field, shared by all the particles, and a Random Uncorrelated distribution which may be characterized in terms of Eulerian fields of particle moments such as the uncorrelated temperature variance. In this paper, the DPS-DNS numerical database is presented, LES Eulerian equations for the dispersed phase are derived in the frame of the Mesoscopic approach and models for the unresolved subgrid and random uncorrelated terms are proposed and a priori tested using the DPS-DNS database.