关于使用AspenPlus回归NRTL二元交互作用参数的实例

1、 1 / 7 关于使用关于使用 aspen plus 回归回归 nrtl 二元交互作用参数的实例二元交互作用参数的实例 天津大学化工学院 08 级化学工程与工艺 2 班 马晨皓 2011-12-12 于天津大学鹏翔公寓 4 斋 本实例采用文献本实例采用文献1所提供的苯所提供的苯-1-庚烯庚烯-n-甲基吡咯烷酮常压下的汽液平衡数据使用甲基吡咯烷酮常压下的汽液平衡数据使用 aspen plus 进行进行 nrtl（非随机两液相）方程二元交互作用参数的回归。该法同样适用于对于其他（非随机两液相）方程二元交互作用参数的回归。该法同样适用于对于其他 例如例如 wilson 方程等交互作用参数的回归计算。

2、方程等交互作用参数的回归计算。 表格 1 苯-1-庚烯-nmp 体系常压下的汽液平衡数据 p=101.3kpa 编号编号 温度温度/ c xb xh yb yh 1 82.91 0.56278 0.06114 0.91421 0.07862 2 83.05 0.85528 0.10407 0.92901 0.05046 3 84.09 0.78580 0.09188 0.88271 0.10812 4 84.32 0.09282 0.31493 0.17092 0.81576 5 85.18 0.77799 0.18362 0.86716 0.11315 6 85.27 0.77434 0.0

3、8763 0.89007 0.10694 7 87.19 0.76951 0.07935 0.94171 0.03957 8 87.46 0.75724 0.10041 0.93659 0.04779 9 87.81 0.54620 0.30181 0.70249 0.28632 10 89.76 0.13965 0.15684 0.35165 0.63094 11 90.48 0.29694 0.43669 0.41571 0.57227 12 90.61 0.30425 0.10845 0.62972 0.34806 13 90.92 0.36833 0.07552 0.75608 0.2

4、1899 14 91.22 0.63803 0.24418 0.82459 0.15863 15 91.61 0.12762 0.39204 0.21003 0.77601 16 92.27 0.03668 0.15026 0.12553 0.84945 17 93.01 0.18450 0.09615 0.53547 0.44722 18 96.38 0.28909 0.06265 0.86106 0.09886 19 96.52 0.21820 0.56594 0.17621 0.78649 20 97.09 0.12362 0.08623 0.62788 0.33372 21 97.98

5、 0.16018 0.31524 0.25134 0.71198 2 / 7 22 98.03 0.09332 0.12175 0.43046 0.53981 1. 选择一个你想用的选择一个你想用的 aspen plus 模板，我通常使用的是通用米制单位模板，选择运行类型模板，我通常使用的是通用米制单位模板，选择运行类型 （run type）为数据回归（）为数据回归（data regression） ；） ； 2. 在组分（在组分（components）中输入各组分。本例使用的是苯）中输入各组分。本例使用的是苯-1-庚烯庚烯-nmp 的的 vle 数据，其中数据，其中 1 表示苯（表示苯（be

6、nzene） ，） ，2 表示表示 1-庚烯（庚烯（1-heptene） ，） ，3 表示表示 n-甲基吡咯烷酮（甲基吡咯烷酮（nmp） ；） ； 3 / 7 3. 在物性 （在物性 （properties） 中全局变量 （） 中全局变量 （global） 的物性方法 （） 的物性方法 （property method） 中选择） 中选择 nrtl； 4. 在数据（在数据（data）中新建一个数据）中新建一个数据 d-1，类型为混合物（，类型为混合物（mixture） ；） ； 4 / 7 5. 在在 data-setup 中导入所有组分，选择数据类型，本例使用的常压下的中导入所有组分，选择数

7、据类型，本例使用的常压下的 txy，故数据类型为，故数据类型为 txy，规定压力为，规定压力为 101.3kpa，根据文献，根据文献1以质量分数为计算基准；以质量分数为计算基准； 6. 在在 data-data 中黏贴数据表格（可以直接将中黏贴数据表格（可以直接将 excel 表黏贴至此） ，注意，表黏贴至此） ，注意，std-dev 一行是标一行是标 准方差；准方差； 5 / 7 7. 在在 regression 中新建一个中新建一个回归计算回归计算 r-1； 8. 在在 regression-setup 中选择模型为中选择模型为 nrtl，将数据，将数据 d-1 填入，选择计算类型为回归计

8、算填入，选择计算类型为回归计算 （regression） ；） ； 6 / 7 9. 在在regression-parameters中选择参数类型为二元参数 （中选择参数类型为二元参数 （binary parameters） ， 名字为） ， 名字为nrtl， 元素为元素为 1（表示回归的参数为标量，（表示回归的参数为标量，nrtl 二元交互作用参数为标量） ，定义组分二元交互作用参数为标量） ，定义组分 i 和和 j，最后，最后 一行的一行的 aij=aji 可以选择可以选择 yes 也可以不选（默认为也可以不选（默认为 no） ；） ； 10. 运行，计算成功后，在运行，计算成功后，在 parameters-binary interaction 中查看回归结果。中查看回归结果。 本例的回归计算结果如下表所示 7 / 7 表格 2 苯-1-庚烯-nmp 体系的 nrtl 二元交互作用参数回归结果一览表 组分组分 i 苯苯 苯苯 1-庚烯庚烯 组分组分 j 1-庚烯庚烯 nmp nmp 温度单位温度单位 c c c aij -0.8133673010 5.756207120 4.010662490 aji 0.0 -3.298700000 0.0 bij 267.4873000 -12