Multi-Energy Systems (MES) represent an increasingly important element in the undergoing energy transition, as they can effectively supply different energy vectors to their users through the coordinated local management of a set of different energy sources. This work presents an innovative Mixed Integer Linear Programming (MILP) formulation for the operational planning of MES featuring thermal generators arranged in series and/or parallel with different and/or variable heat delivery temperatures and filled-tank thermal energy storages. The proposed model can approximate the non-isothermal mixing processes occurring in water collectors and storage systems through the definition of a set of virtual headers at constant temperature levels. To a...