expr += (
                    list(n.heat_output.values())[0].Q_CW_min[t] * self.Y[n, t]