Numerical approximation of non-isothermal multi-component, multi-phase field systems