constr log(aveopConsumer(exp(pExp*(1+P/consumerPower)^2)))/pExp Constraint pipePrice 2*intopHot(pipeCost(pipeD)) "Pipe price [EUR]" power1 intopInlet1((p-1[atm])*nipfl.Qv) "Pumping power" power2 intopInlet2((p-1[atm])*nipfl.Qv) "Pumping power" pumpPrice pumpWeight*(power1+power2) "Pump price [EUR]" obj pipePrice+pumpPrice "Objective function" flow1in intopInlet1(nipfl.Qm) "Mass flow rate, Inlet 1" flow2in intopInlet2(nipfl.Qm) "Mass flow rate, Inlet 2" Pin flow1in*intopOutlet1((T1-T)*mat1.def.Cp(T))+flow2in*intopOutlet2((T2-T)*mat1.def.Cp(T)) "Energy input" Pconsumer intopConsumers(-P/Lz) "Total consumer power" Efficiency Pconsumer/Pin "Network efficiency" Pout intopAll(-qHeat) "Total heat loss"