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