如何解决gekko中的简单混合操作？答案

【问题标题】：how to solve a simple mixing operation in gekko?如何解决gekko中的简单混合操作？
【发布时间】：2019-07-26 11:59:48
【问题描述】：

我正在尝试解决 gekko 中的简单混合操作。混合器mx 采用两个入口流Feed1 和Feed2。预期结果是出口流mx.outlet 的质量流量应该是入口流质量流量的总和。这是我尝试过的。

from gekko import GEKKO, chemical
m = GEKKO(remote=False)    

f = chemical.Flowsheet(m)
P = chemical.Properties(m)

c1 = P.compound('Butane')
c2 = P.compound('Propane')

feed1 = f.stream()
m_feed1 = f.massflows(sn= feed1.name)
m_feed1.mdot = 200
m_feed1.mdoti = [50,150]

feed2= f.stream()
m_feed2 = f.massflows(sn= feed2.name)
m_feed2.mdot = 200
m_feed2.mdoti = [50,150]      

mx = f.mixer(ni=2)
mx.inlet = [feed1.name,feed2.name]
m.options.SOLVER = 1

mf = f.massflows(sn = mx.outlet)
m.solve()

代码运行成功。但是，在mf.mdot 上似乎输出了不正确的值 [-1.8220132454e-06]。预期值为 400。任何帮助，我的代码有什么问题？

【问题讨论】：

问题是您的质量流量mdot 和组件质量流量mdoti 变量是在需要修复时计算的。还计算了入口流上的摩尔流量变量。这是一个翻转变量规格（固定/计算）的问题。我需要让这些更容易使用，所以我要向 Gekko 添加一些新功能，这些新功能将在下一个版本中提供。我将推迟对这个问题的回答，直到准备好更简单的修复程序。化学模块是新的，需要更多的开发以方便用户使用。

标签： gekko

【解决方案1】：

以下是适用于此混合应用程序的源代码：

from gekko import GEKKO, chemical
import json

m = GEKKO(remote=False)    

f = chemical.Flowsheet(m)
P = chemical.Properties(m)

# define compounds
c1 = P.compound('Butane')
c2 = P.compound('Propane')

# create feed streams
feed1 = f.stream(fixed=False)
feed2 = f.stream(fixed=False)

# create massflows objects
m_feed1 = f.massflows(sn=feed1.name)
m_feed2 = f.massflows(sn=feed2.name)

# create mixer
mx = f.mixer(ni=2)

# connect feed streams to massflows objects
f.connect(feed1,mx.inlet[0])
f.connect(feed2,mx.inlet[1])
m.options.SOLVER = 1

mf = f.massflows(sn = mx.outlet)

# specify mass inlet flows
mi = [50,150]
for i in range(2):
    m.fix(m_feed1.mdoti[i],val=mi[i])
    m.fix(m_feed2.mdoti[i],val=mi[i])
# fix pressure and temperature
m.fix(feed1.P,val=101325)
m.fix(feed2.P,val=101325)
m.fix(feed1.T,val=300)
m.fix(feed2.T,val=305)

m.solve(disp=True)

# print results
print(f'The total massflow out is {mf.mdot.value}')

print('')

# get additional solution information
with open(m.path+'//results.json') as f:
    r = json.load(f)
for name, val in r.items():
    print(f'{name}={val[0]}')

下面是求解器的输出。这仅适用于 APM 0.9.1 和 Gekko v0.2.3（将于 2019 年 8 月发布）。 thermo 和 flowsheet 对象库随 v0.2.2 发布，还有几个功能仍在开发中。下一个版本应该会解决其中的许多问题。

 ----------------------------------------------------------------
 APMonitor, Version 0.9.1
 APMonitor Optimization Suite
 ----------------------------------------------------------------


 --------- APM Model Size ------------
 Each time step contains
   Objects      :  6
   Constants    :  0
   Variables    :  19
   Intermediates:  0
   Connections  :  44
   Equations    :  2
   Residuals    :  2

 Number of state variables:    14
 Number of total equations: -  14
 Number of slack variables: -  0
 ---------------------------------------
 Degrees of freedom       :    0

 ----------------------------------------------
 Steady State Optimization with APOPT Solver
 ----------------------------------------------

 Iter    Objective  Convergence
    0  3.86642E-16  1.99000E+02
    1  4.39087E-18  1.11937E+01
    2  8.33448E-19  6.05819E-01
    3  1.84640E-19  1.62783E-01
    4  2.91981E-20  7.21250E-02
    5  1.55439E-21  2.28110E-02
    6  5.51232E-24  1.21437E-03
    7  7.03139E-29  4.30650E-06
    8  7.03139E-29  4.30650E-06
 Successful solution

 ---------------------------------------------------
 Solver         :  APOPT (v1.0)
 Solution time  :  0.0469 sec
 Objective      :  0.
 Successful solution
 ---------------------------------------------------


v1 not found in results file
The total massflow out is [400.0]

time=0.0
feed1.h=44154989.486
feed1.x[2]=0.79815448476
feed1.vdot=104.9180373
feed1.dens=0.040621756423
feed1.c[1]=0.0081993193551
feed1.c[2]=0.032422437068
feed1.mdot=200.0
feed1.y[1]=0.25
feed1.y[2]=0.75
feed1.sfrc=0.0
feed1.lfrc=0.0
feed1.vfrc=1.0
feed2.h=44552246.421
feed2.x[2]=0.79815448476
feed2.vdot=106.66667125
feed2.dens=0.03995582599
feed2.c[1]=0.0080649042837
feed2.c[2]=0.031890921707
feed2.mdot=200.0
feed2.y[1]=0.25
feed2.y[2]=0.75
feed2.sfrc=0.0
feed2.lfrc=0.0
feed2.vfrc=1.0
mixer5.outlet.t=381.10062836
mixer5.outlet.h=44353617.96
mixer5.outlet.ndot=8.5239099109
mixer5.outlet.x[1]=0.20184551524
mixer5.outlet.x[2]=0.79815448476
mixer5.outlet.vdot=1.5797241143
mixer5.outlet.dens=5.5635215396
mixer5.outlet.c[1]=1.0891224437
mixer5.outlet.c[2]=4.3066994177
mixer5.outlet.mdot=400.0
mixer5.outlet.y[1]=0.25
mixer5.outlet.y[2]=0.75
mixer5.outlet.sfrc=0.0
mixer5.outlet.lfrc=1.0
mixer5.outlet.vfrc=0.0
v2=300.0
v3=4.2619549555
v4=0.20184551524
v5=0.79815448476
v6=101325.0
v7=305.0
v8=4.2619549555
v9=0.20184551524
v10=0.79815448476
v11=200.0
v12=50.0
v13=150.0
v14=200.0
v15=50.0
v16=150.0
v17=400.0
v18=100.0
v19=300.0

【讨论】：

谢谢@John，我会等待下一个版本。从上面的代码看来，用户必须手动固定压力和温度。或者mixer对象有什么方法可以根据热平衡解决温度问题？
是的，温度是根据热平衡来解决的。我固定了入口温度，但计算了出口温度。您可以使用：print(f"The outlet temperature is {r[mx.name+'.outlet.t']}") 显示出口温度。在求解之前，您还可以将变量连接到出口温度。这使您可以在计算入口温度以匹配固定出口的同时修复出口。需要注意的一件事是，我的两个入口流都是蒸汽，出口是液体，所以温度升高了。您需要进行 VLE 计算或指定阶段。
我将在 0.2.3 版中将所有型号的默认设置切换为液体。否则，您可以使用f.set_phase(y,phase='vapor'): 指定阶段，如文档中所示：gekko.readthedocs.io/en/latest/chemical.html
太棒了！这些都是使 gekko 成为成熟的过程模拟、控制和优化工具的真正令人兴奋的新功能。