| 查看: 1841 | 回复: 10 | |||
[求助]
关于二元混合溶剂溶解度数据利用wilson、NRTL模拟的1stOPT代码 已有2人参与
|
|
各位大神，我测得了不同比例的某物质在二元混合溶剂中的溶解度随温度的变化，欲选用wilson、NRTL方程对其回归求得参数，下面是我的数据，但是1stOp软件t程序不大会写，总是提示function有错，请求高手指教 x2=0 283.3 3.02194E-05 288.31 0.000157704 293.3 0.0004041 298.25 0.000686379 303.15 0.001026437 308.18 0.001405125 313.05 0.001810204 318.1 0.002726428 323.05 0.003769643 x2=0.05 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 x2=0.1 283.2 0.000828269 288.23 0.001597956 293.16 0.002383646 298.24 0.003548518 303.2 0.004900833 308.15 0.006296129 313.15 0.008064994 318.2 0.010216362 323.2 0.013075324 x2=0.15 283.15 0.002970424 288.15 0.003905376 293.06 0.005280804 298.15 0.006616207 303.15 0.008243561 308.15 0.010303864 313.2 0.013119118 318.15 0.016480102 323.07 0.020616259 x2=0.2 283.1 0.005204799 288.15 0.006279544 293.05 0.007921776 298.05 0.00986987 303.15 0.012566941 308.15 0.015421326 313.18 0.019498661 318.25 0.024651302 323.15 0.030872327 x2=0.25 283.22 0.008048279 288.27 0.009556376 293.3 0.011981522 298.2 0.014520227 303.2 0.01789219 308.1 0.022389299 313.05 0.027783872 318.07 0.035092411 323.1 0.044652622 x2=0.3 283.05 0.011149853 288.05 0.013651819 293.06 0.016540268 298.07 0.020332495 303.1 0.024782096 308.23 0.030544405 313.25 0.037978064 318.35 0.048406476 323.2 0.061434365 x2=0.35 283.15 0.015106036 288.22 0.018196111 293.26 0.021838514 298.2 0.026852153 303.07 0.032794157 308.15 0.040592546 313 0.050666982 318.1 0.064596043 323.1 0.083519503 x2=0.4 283.1 0.020213298 288.15 0.024960264 293.25 0.029919577 298.3 0.035835999 303.2 0.042800335 308.25 0.052897495 313.17 0.065920884 318.21 0.082267938 323.2 0.106310795 x2=0.45 283.09 0.022163167 288.1 0.028017074 293.1 0.034843674 298.25 0.0420929 303.15 0.051962363 308.15 0.065005896 313.2 0.08215546 318.25 0.104250344 323.1 0.137575799 x2=0.5 283.15 0.027972969 288.13 0.03399462 293.05 0.040605327 298.01 0.048112851 303.1 0.058477267 308.05 0.072571638 313.2 0.09047995 318.25 0.114083364 323.15 0.147697283 下面是我的程序 VarConstant x2=; Constant R=8.314,V1=140.118,V2=97.9,V3=108.1053,H=23100,Tm=349.8; Parameters g12,g21,g31,g13,g32,g23; Variable T,x1; constStr a21=V1/V2*exp(-g21/R/T),a12=V2/V1*exp(-g12/R/T),a31=V1/V3*exp(-g31/R/T),a13=V3/V1*exp(-g13/R/T),a23=V3/V2*exp(-g23/R/T),a32=V2/V3*exp(-g32/R/T); constStr lnr=1-ln(x1+a12*x2+a13*(1-x1-x2))-x1/(x1+a12*x2+a13*(1-x1-x2))-a21*x2/(x2+a21*x1+a23(1-x1-x2))-a31*(1-x1-x2)/((1-x1-x2)+a31*x1+a32*x2); constStr x1+x2+x3=1; Function x1=exp(-H/R*(1/T-1/Tm)-lnr); Data; //T,x1 283.3 3.02194E-05 288.31 0.000157704 293.3 0.0004041 298.25 0.000686379 303.15 0.001026437 308.18 0.001405125 313.05 0.001810204 318.1 0.002726428 323.05 0.003769643 data; 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 data; 283.2 0.000828269 288.23 0.001597956 293.16 0.002383646 298.24 0.003548518 303.2 0.004900833 308.15 0.006296129 313.15 0.008064994 318.2 0.010216362 323.2 0.013075324 data; 283.15 0.002970424 288.15 0.003905376 293.06 0.005280804 298.15 0.006616207 303.15 0.008243561 308.15 0.010303864 313.2 0.013119118 318.15 0.016480102 323.07 0.020616259 data; 283.1 0.005204799 288.15 0.006279544 293.05 0.007921776 298.05 0.00986987 303.15 0.012566941 308.15 0.015421326 313.18 0.019498661 318.25 0.024651302 323.15 0.030872327 data; 283.22 0.008048279 288.27 0.009556376 293.3 0.011981522 298.2 0.014520227 303.2 0.01789219 308.1 0.022389299 313.05 0.027783872 318.07 0.035092411 323.1 0.044652622 data; 283.05 0.011149853 288.05 0.013651819 293.06 0.016540268 298.07 0.020332495 303.1 0.024782096 308.23 0.030544405 313.25 0.037978064 318.35 0.048406476 323.2 0.061434365 data; 283.15 0.015106036 288.22 0.018196111 293.26 0.021838514 298.2 0.026852153 303.07 0.032794157 308.15 0.040592546 313 0.050666982 318.1 0.064596043 323.1 0.083519503 data; 283.1 0.020213298 288.15 0.024960264 293.25 0.029919577 298.3 0.035835999 303.2 0.042800335 308.25 0.052897495 313.17 0.065920884 318.21 0.082267938 323.2 0.106310795 data; 283.09 0.022163167 288.1 0.028017074 293.1 0.034843674 298.25 0.0420929 303.15 0.051962363 308.15 0.065005896 313.2 0.08215546 318.25 0.104250344 323.1 0.137575799 data; 283.15 0.027972969 288.13 0.03399462 293.05 0.040605327 298.01 0.048112851 303.1 0.058477267 308.05 0.072571638 313.2 0.09047995 318.25 0.114083364 323.15 0.147697283 不知道哪里出现问题了，请求指教！ |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有239人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 自己这一段时间过柱子的总结
已经有11人回复
2楼2014-07-14 15:08:53
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68131.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2014-07-14 15:33:47
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-07-14 19:32:07
|
大神您好，跪求指导，是这样的：我做的是固体物质A在不同比例的两种混合溶剂（B和C）中的溶解度数据回归 x1是溶液中溶质的摩尔分率，x2、x3是溶液中两种溶剂B和C的摩尔分率，即有关系x1+x2+x3=1，其实我上边列出的那一大列数据中的x2应该改为xB，是溶剂B在混合溶剂B+C中所占比例，xB+xC=1，因为随着溶质A不断加入溶液中，x1、x2、x3都在变化，不是固定值。 比如这一组数据， xB=0.05 T xA 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 是指在混合溶剂中B摩尔分率为0.05时，A在混合溶剂B+C中摩尔溶解度x1随温度变化。 单对此组数据写程序还好，但是这样的话每个组成下都会拟合出6个参数，上面的10组数据就会出现10组参数，没法统一，我想整体拟合，最终只有6个参数，请求指导 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 溶解度模型.JPG
2014-07-15 11:18:03, 46.48 K
5楼2014-07-15 11:15:29
6楼2014-07-15 11:18:07
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
7楼2014-07-15 21:16:47
8楼2014-07-16 14:14:38
9楼2014-07-18 09:52:45
萧之夏
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 796.3
- 散金: 217
- 帖子: 208
- 在线: 149.3小时
- 虫号: 2160502
- 注册: 2012-12-02
- 性别: GG
- 专业: 化工热力学和基础数据
10楼2016-03-29 20:20:30
