This article presents a new numerical scheme designed to solve for any scalar equation coupled with a Lattice-Boltzmann solver (in so-called hybrid methods). Its most direct application is to solve an energy equation, in parallel with a Lattice-Boltzmann solver dealing with mass and momentum conservation.The numerical scheme is specifically designed to compute the energy flux consistently with the mass and momentum flux (as is done, for instance, using Riemann solvers).This scheme effectively lifts a major limitation of current compressible hybrid Lattice-Boltzmann, in which the energy conservation is tackled under nonconservative form, leading to discretization errors on jump conditions across shocks.Combined with our recently presented pressure-based solver [G. Farag et al, Physics of Fluids, vol. 32, no. 6, p. 066106, (2020)], the resulting hybrid Lattice-Boltzmann scheme is, to the authors' knowledge, the first to numerically conserve simultaneously mass, momentum and total energy.