Gekko и CoolProp - PullRequest
Новая
       107

Gekko и CoolProp

2 голосов
/ 10 июля 2020

Я использую GEKKO и CoolProp для моделирования тепловых систем. При попытке использовать функции CoolProp в уравнениях модели (как показано ниже для расширения isentropi c) я получаю сообщение об ошибке относительно типа переменной: « должно быть действительным числом, а не GKVariable ». Может ли кто-нибудь помочь мне с этой проблемой? 

from gekko import GEKKO
import CoolProp.CoolProp as CP 
#
p1 = 2e5
T1 = 300.0 + 273.15
p2 = 1e5
eta = 0.80
fluid = 'H2O'
#
h1 = CP.PropsSI('H','T',T1,'P',p1,fluid)
s1 = CP.PropsSI('S','T',T1,'P',p1,fluid)
#
m = GEKKO()
h2 = m.Var()
h2s = m.Var()
T2 = m.Var()
#
m.Equation(eta * (h1 - h2) - (h1 - h2s) == 0)
m.Equation(h2s - CP.PropsSI('H','S',s1,'P',p2,fluid) == 0)
m.Equation(h2 - CP.PropsSI('H','T',T2,'P',p2,fluid) == 0)
#
m.options.IMODE = 1  #Steady state
m.options.SOLVER = 3 # solver (IPOPT) 
m.solve(disp=False)
#
print(T2.value[0])

Заранее спасибо.

1 Ответ

2 голосов
/ 14 июля 2020

Gekko необходимо иметь уравнения CoolProp для выполнения c автоматического дифференцирования. Если у него нет уравнения, вы можете использовать cspline ( cubi c spline documentation ). 

from gekko import GEKKO
import CoolProp.CoolProp as CP 
import numpy as np
#
p1 = 2e5
T1 = 300.0 + 273.15
p2 = 1e5
eta = 0.80
fluid = 'H2O'
# constants
h1 = CP.PropsSI('H','T',T1,'P',p1,fluid)
s1 = CP.PropsSI('S','T',T1,'P',p1,fluid)
c2 = CP.PropsSI('H','S',s1,'P',p2,fluid)
# gekko model
m = GEKKO(remote=False)
h2 = m.Var()
h2s = m.Var()
T2 = m.Var()
# build cubic spline
n = 100
T = np.linspace(373,1000,100)
h = [CP.PropsSI('H','T',Ti,'P',p1,fluid) for Ti in T]
m.cspline(T2,h2,T,h)
#
m.Equation(eta * (h1 - h2) - (h1 - h2s) == 0)
m.Equation(h2s - c2 == 0)
#
m.options.IMODE = 1  #Steady state
m.options.SOLVER = 3 # solver (IPOPT) 
m.solve(disp=False)
#
print(T2.value[0])

Это дает решение: 

T2 = 468.52459939
...