ISLEC-Web集成溶液平衡计算系统用户手册
1.0版
李东东 博士 编制
中国科学院青海盐湖研究所
2018.08.12
版权所有 盗版必究
1. ISLEC-Web集成溶液平衡计算系统简介
ISLEC-Web在线计算服务授权用户可以通过Web浏览器(IE、Firefox、Google chrome等)访问网址www.islec.net 使用ISLEC-Web进行计算。在访问到ISLEC-Web主服务器后(图1.1),用鼠标点击LOGIN(用户登录)菜单,进入登录页面(图1.2)。之后输入授权协议中获得的用户名和密码,进入ISLEC-Web在线计算页面(图1.3)。
图1.1 ISLEC-Web主服务器访问页面
图1.2 ISLEC-Web用户登录页面
图1.3 ISLEC-Web在线计算服务页面
2. 在线计算服务模型及使用方法
ISLEC-Web在线计算服务授权用户可以获得基于ISLEC主程序和SLMB-G.db热力学数据库所提供的全部服务。SLMB-G.db热力学数据库包含9种化学元素(H, O, Li, Na, Mg, S, Cl, K和Ca),8种水溶液物种(Li+(aq), Na+(aq), K+(aq), Mg2+(aq), Ca2+(aq), Cl−(aq), SO42−(aq)和H2O),以及54种固相(见表2.1)。
表2.1 SLMB-G.db热力学数据库中包含的全部固相物种
H2O(cr,Ih) |
LiCl·H2O(cr) |
LiCl·2H2O(cr) |
LiCl·3H2O(cr) |
LiCl·5H2O(cr) |
NaCl(cr) |
NaCl·2H2O(cr) |
KCl(cr) |
MgCl2·12H2O(cr) |
MgCl2·8H2O(cr) |
MgCl2·6H2O(cr) |
MgCl2·4H2O(cr) |
MgCl2·2H2O(cr) |
CaCl2·6H2O(cr) |
CaCl2·4H2O(cr,a) |
CaCl2·4H2O(cr,b) |
CaCl2·4H2O(cr,g) |
CaCl2·2H2O(cr) |
CaCl2·H2O(cr) |
CaCl2·1/3H2O(cr) |
Li2SO4·H2O(cr) |
Na2SO4·10H2O(cr) |
Na2SO4(cr,V) |
K2SO4(cr) |
MgSO4·11H2O(cr) |
MgSO4·7H2O(cr) |
MgSO4·6H2O(cr) |
MgSO4·5H2O(cr) |
MgSO4·4H2O(cr) |
MgSO4·H2O(cr) |
CaSO4·2H2O(cr) |
CaSO4·1/2H2O(cr) |
CaSO4(cr) |
LiCl·MgCl2·7H2O(cr) |
LiCl·CaCl2·5H2O(cr) |
KCl·MgCl2·6H2O(cr) |
KCl·CaCl2(cr) |
2MgCl2·CaCl2·12H2O(cr) |
LiKSO4(cr) |
Li2SO4·Na2SO4(cr) |
NaK3(SO4)2(cr) |
K2SO4·MgSO4·6H2O(cr) |
K2SO4·MgSO4·4H2O(cr) |
K2SO4·CaSO4·H2O(cr) |
K2SO4·5CaSO4·H2O(cr) |
Li2SO4·3Na2SO4·12H2O(cr) |
Na2SO4·MgSO4·4H2O(cr) |
Na2SO4·CaSO4(cr) |
2Na2SO4·CaSO4·2H2O(cr) |
Na2SO4·5CaSO4·3H2O(cr) |
K2SO4·2MgSO4(cr) |
3Na2SO4·MgSO4(cr) |
Na2SO4·MgSO4·2.5H2O(cr) |
K2SO4·MgSO4·2CaSO4·2H2O(cr) |
该热力学数据库的适用温度范围为−50 ○C~100 ○C,个别情况的温度适用范围略有差异。该数据库的检验的结果详见[Li-Na-K-Mg-Ca-Cl-SO4-H2O体系多温热力学相平衡模型开发及其应用, 2015, 博士学位论文, 中国科学院大学]。尽管如此,受限于热力学数据库构建时实验数据的缺乏、热力学模型外推能力、平衡求解算法的局限性、以及不全面的检验,我们无法确保该模型在任何所声明的条件下都能够给出完全无误的结果。此外,由于该热力学数据库的构建利用了大量不同来源的实验数据及数据评估方法,可能与用于所能找到的实验数据存在差异。
基于ISLEC桌面和服务器应用软件4.1版和SLMB-G.db热力学数据库,ISLEC-Web为用户提供四种独立在线计算模块(图2.1)。
图2.1 ISLEC-Web在线计算服务四大功能模块
2.1 单点平衡计算(Single Point Equilibration Calculation)
用户可以通过鼠标点击“Single Point Equilibration Calculation”选项进入单点平衡在线计算页面(图2.2)。该部分计算功能对应于ISLEC桌面版和服务器版软件中的“EQUILIBRIUM”指令功能,以实现对给定化学体系在给定温度和压力条件下的化学和相平衡进行计算。
当进行单点平衡计算时,ISLEC-Web自动略去了经常出现过饱和状态的固相(包括MgSO4·H2O、CaSO4·1/2H2O、CaSO4、Na2SO4·MgSO4·4H2O、Na2SO4·CaSO4、2Na2SO4·CaSO4·2H2O、Na2SO4·5CaSO4·3H2O、K2SO4·2MgSO4、3Na2SO4·MgSO4、Na2SO4·MgSO4·2.5H2O、K2SO4·MgSO4·2CaSO4·2H2O共11个),即执行该种假定条件下的亚稳平衡计算。如需计算体系的完全热力学平衡态,请参见2.4节。
图2.2 ISLEC-Web在线计算服务单点平衡页面
ISLEC-Web用户界面简洁,便于理解,使用过程简单。用户只需要输入温度、压力、体系总电荷(默认值为0,推荐接受默认值;几乎所有水样分析实验结果的总电荷均不为0,需人为设定差减组分浓度)、以及体系个组分的摩尔数(如图2.3),之后点击”计算(Calculate)“按钮,ISLEC-Web自动以标准输出文本格式给出平衡计算结果(如图2.4)。用户可以根据自己需要对文本输出进行数据抽取、后续加工及图形化。图形化工具可选excel、origin、gnuplot等。
尽管压力”Pressure (Unit is Kilopascal, kPa) ”作为可变参数开放给用户,但是由于SLMB-G.db热力学数据库仅仅考虑体系的常压(101.325 kPa)变温性质,因此无论用户如何改变压力值,ISLEC-Web均视为在101.325 kPa下进行平衡计算。
图2.3 ISLEC-Web在线计算服务单点平衡输入示例
图2.4 ISLEC-Web在线计算服务单点平衡输出示例
2.2 卤水等温蒸发计算(Brine Isothermal Evaporation Calculation)
卤水等温蒸发计算模块的ISLEC软件和热力学数据库设置与单点平衡计算完全相同。当进行等温蒸发计算时,ISLEC-Web也自动略去了经常出现过饱和状态的固相(包括MgSO4·H2O、CaSO4·1/2H2O、CaSO4、Na2SO4·MgSO4·4H2O、Na2SO4·CaSO4、2Na2SO4·CaSO4·2H2O、Na2SO4·5CaSO4·3H2O、K2SO4·2MgSO4、3Na2SO4·MgSO4、Na2SO4·MgSO4·2.5H2O、K2SO4·MgSO4·2CaSO4·2H2O共11个),即执行该种假定条件下的亚稳平衡计算。
在利用等温蒸发在线计算模块时,用户操作界面如图2.5。用户不仅需要提供蒸发温度、压力和卤水(或卤水与盐的混合物)的初始体系组成,还需要给定蒸发结束时体系中的水量(该值需大于0)。ISLEC-Web将用户设置的始末水量按照相同间隔等分为十份,并对相应的十个阶段进行独立平衡计算。在每步的计算过程中,考虑平衡结晶,而非分步结晶,即后一步只考虑前一步蒸发掉水的情况,而不将之前可能析出的固体分离出体系。这时之前析出的固相仍然允许和卤水重新平衡。如需模拟分步结晶的情况,建议用户采用单点平衡计算模块,手动将析出固相从体系中移除。
图2.5 ISLEC-Web在线计算服务等温蒸发输入示例
由ISLEC-Web产生的连续等温蒸发平衡结果仍然以ISLEC标准文本输出形式给出。用户可以将结果区输(图2.6)出保存为文本文件(鼠标邮件点击结果区->查看框架源代码->Ctrl+A/Ctrl+C复制所有内容至文本文件),并采用文本处理程序,如awk,进行后续加工和格式化。之后,再利用绘图软件对结果进行图形表达。
图2.6 ISLEC-Web在线计算服务等温蒸发输出示例
2.3 卤水连续冷冻计算(Brine Freezing Calculation)
卤水连续冷冻计算模块的ISLEC软件和热力学数据库设置与单点平衡计算和等温蒸发计算完全相同。当进行连续冷冻计算时,ISLEC-Web也自动略去了经常出现过饱和状态的固相(包括MgSO4·H2O、CaSO4·1/2H2O、CaSO4、Na2SO4·MgSO4·4H2O、Na2SO4·CaSO4、2Na2SO4·CaSO4·2H2O、Na2SO4·5CaSO4·3H2O、K2SO4·2MgSO4、3Na2SO4·MgSO4、Na2SO4·MgSO4·2.5H2O、K2SO4·MgSO4·2CaSO4·2H2O共11个),即执行该种假定条件下的亚稳平衡计算。
在利用连续冷冻在线计算模块时,用户操作界面如图2.7。用户不仅需要提供初始体系温度、压力和组成,还需要指定冷冻结束温度(系统不对该值进行检测,其有效性须由用户自己确定,不建议小于223.15 K)。ISLEC-Web将用户设置的始末温度按照相同间隔等分为十份,并对相应的十个阶段进行独立平衡计算。在每步的计算过程中,考虑平衡结晶,而非分步结晶,即后一步只考虑系统改变温度的情况,而不将之前可能已经析出的固体分离出体系。这时之前析出的固相仍然允许和卤水重新平衡。如需模拟分步结晶的情况,建议用户采用单点平衡计算模块,手动将析出固相从体系中移除。
图2.7 ISLEC-Web在线计算服务连续冷冻输入示例
图2.8 ISLEC-Web在线计算服务连续冷冻输出示例
由ISLEC-Web产生的连续冷冻平衡结果仍然以ISLEC标准文本输出形式给出。用户可以将结果区输(图2.8)出保存为文本文件(鼠标邮件点击结果区->查看框架源代码->Ctrl+A/Ctrl+C复制所有内容至文本文件),并采用文本处理程序,如awk,进行后续加工和格式化。之后,再利用绘图软件对结果进行图形表达。
2.4 变质作用单点计算(Brine/Mineral Metamorphism Calculation)
变质作用单点计算模块原理和使用方法与单点平衡计算(Single Point Equilibration Calculation)模块完全相同,只是在该模块中体系包含所有的固相全部存在于热力学数据库中,即执行热力学稳定平衡计算。
3. ISLEC-Web应用程序接口
ISLEC-Web可以为用户提供基于PHP的应用程序接口(API),接口协议和加密方式需协商确定后定制开发和开放授权。
4. ISLEC-Web用户培训案例
- 案例一:单点平衡计算
1、访问ISLEC-Web在线计算服务网站:https://www.islec.net
2、点击LOGIN(用户登录),输入ISLEC-Web授权账号和密码进行登录;选择LOGIN(用户登录)菜单下MY ISLEC(我的ISLEC)进入ISLEC在线计算功能选择。
3、选择单点计算(Single Point Equilibration Calculation)模块,进入计算页面。
4、输入计算的温度和压强参数,默认温度为298.15K,压强为101.325KPa。这里以默认温度和压强为例,温度和压强的适用范围参考用户手册相应内容。
5、体系默认的总电荷为零,这里若无特殊计算需要不建议更改。每次参数输入完成后应检查体系电荷是否平衡,若电荷不平衡,会有提示字样出现或无法正常计算。
6、输入体系中各组分的含量,包含水以及各种离子的物质的量(摩尔数)。这里输入水的物质的量为55.51 mol,Na+的物质的量为10 mol,K+的物质的量为1.0 mol,Cl−的物质的量为11 mol。输入完成后,点击Calculate进行计算。
7、输出的计算结果包括平衡液相中各组分的含量、活度系数和活度,平衡固相的组成成分和含量。如上例,输出结果中平衡液相的组成信息如下:
S | Z | n(mol) | m(mol/kgw) | x | a | r |
H2O(aq) | +0 | 5.551000e+01 | 5.551000e+01 | 8.068715e-01 | 7.381805e-01 | 9.148676e-01 |
Na+(aq) | +1 | 5.643291e+00 | 5.643291e+00 | 8.202865e-02 | 5.798319e+00 | 1.027484e+00 |
K+(aq) | +1 | 1.000000e+00 | 1.000000e+00 | 1.453561e-02 | 5.063719e-01 | 5.063782e-01 |
Cl–(aq) | -1 | 6.643291e+00 | 6.643291e+00 | 9.656426e-02 | 6.608013e+00 | 9.947023e-01 |
表中S表示溶液相物种,Z表示电荷,n表示物质的量,m表示质量摩尔浓度,x表示摩尔分数,a表示活度,r表示活度系数。
输出结果中平衡固相的组成信息如下:
S | Z | n(mol) |
NaCl(cr) | +0 | 4.356709e+00 |
其中S表示物相,Z表示电荷,n表示物质的量。
如上表所示,析出的固相为无水氯化钠,物质的量4.356709 mol。若无平衡固相的输出结果,说明无固相能够达到饱和析出的条件,体系平衡时处于溶液状态。
- 案例二:卤水等温蒸发计算
1、访问ISLEC-Web在线计算服务网站:https://www.islec.net
2、点击LOGIN(用户登录),输入ISLEC-Web授权账号和密码进行登录;选择LOGIN(用户登录)菜单下MY ISLEC(我的ISLEC)进入ISLEC在线计算功能选择。
3、选择卤水等温蒸发计算(Brine Isothermal Evaporation Calculation)模块,进入计算页面。
4、输入计算的卤水蒸发温度和压强参数,默认温度为298.15K,压强为101.325KPa。这里以默认温度和压强为例,温度和压强的适用范围参考用户手册相应内容。
5、体系默认的总电荷为零,这里若无特殊计算需要不能更改。每次输入完成后应检查电荷是否平衡,若电荷不平衡,会有提示字样出现或无法正常计算。
6、输入卤水蒸发前的体系组成,包含水以及各种离子的物质的量,以及剩余水的物质的量。输入蒸发前的体系物分组成,水的物质的量为55.51 mol,Na+的物质的量为3 mol,K+的物质的量为3 mol,Cl−的物质的量为6 mol。输入剩余水的物质的量为20 mol。输入完成后,点击Calculate进行计算。
7、该计算过程,将水的蒸发量等分为十份,对每一步进行优化计算并输出结果。每一步等温蒸发的输出结果,包括是否优化成功,蒸发前平衡液相的组分含量、活度系数和活度,蒸发前平衡固相的组成成分和含量。
第一步的计算结果展示如下:
优化结果 | SUCCESSED IPOPT OPTIMIZATION
SUCCESSED IPOPT OPTIMIZATION |
如果优化失败会输出”SUCCESSED IPOPT OPTIMIZATION”和” ERROR OCCURRED DURING IPOPT OPTIMIZATION WITH WARM START “字样,同时不输出其它结果,继续进行下一步优化计算。
初始体系平衡液相的组成信息如下:
S | Z | n(mol) | m(mol/kgw) | x | a | r |
H2O(aq) | +0 | 5.551000e+01 | 5.551000e+01 | 8.222485e-01 | 7.837040e-01 | 9.531230e-01 |
Na+(aq) | +1 | 3.000000e+00 | 3.000000e+00 | 4.443786e-02 | 2.599406e+00 | 8.664797e-01 |
K+aq) | +1 | 3.000000e+00 | 3.000000e+00 | 4.443786e-02 | 1.650184e+00 | 5.500684e-01 |
Cl–(aq) | -1 | 6.000000e+00 | 6.000000e+00 | 8.887573e-02 | 4.732940e+00 | 7.888332e-01 |
其中S表示物相,Z表示电荷,n表示物质的量,m表示浓度,x表示摩尔分数,a表示活度,r表示活度系数。
第一步蒸发前平衡固相的组成信息如下:
S | Z | n(mol) |
SOLID 1,H2O(cr,Ih) | +0 | 0.000000e+00 |
SOLID 6,NaCl(cr) | +0 | 0.000000e+00 |
SOLID 7,NaCl*2H2O(cr) | +0 | 0.000000e+00 |
SOLID 8,KCl(cr) | +0 | 0.000000e+00 |
其中S表示物相,Z表示电荷,n表示物质的量。在平衡固相结果输出行列,共有43个固相及相应的含量。你可以依据输入的物质元素查找相应的固相及含量。
如上表所示,四种相应的固相均为零,说明第一步蒸发前无固相析出。如果想要直接得到卤水等温蒸发的最终结果,可直接浏览最后一步的输出结果即可。
- 案例三:卤水冷冻计算
1、访问ISLEC-Web在线计算服务网站:https://www.islec.net
2、点击LOGIN(用户登录),输入ISLEC-Web授权账号和密码进行登录;选择LOGIN(用户登录)菜单下MY ISLEC(我的ISLEC)进入ISLEC在线计算功能选择。
3、选择卤水冷冻计算(Brine Freezing Calculation)模块,进入计算页面。
4、输入卤水冷冻的初始温度、终点温度以及压强,默认初始温度和终点温度均为298.15K,压强为101.32kpa。这里以默认初始温度、终点温度为248.15K,默认压强为例。
5、体系默认的总电荷为零,这里若无特殊计算需要不能更改。每次输入完成后应检查电荷是否平衡,若电荷不平衡,会有提示字样出现或无法正常计算。
6、输入卤水冷冻的体系物分组成,包含水以及各种离子的物质的量。输入水的物质的量为55.51 mol,Na+的物质的量为3 mol,K+的物质的量为3 mol,Cl−的物质的量为6 mol。输入完成后,点击Calculate进行计算。
7、该计算过程,将降温的温度范围等分为十份,对每一步进行优化计算并输出结果。每一步冷冻的输出结果,包括是否优化成功,降温前平衡液相的组分含量、活度系数和活度,降温前平衡固相的组成成分和含量。
第一步的计算结果展示如下:
优化结果 | SUCCESSED IPOPT OPTIMIZATION
SUCCESSED IPOPT OPTIMIZATION |
如果优化失败会输出”SUCCESSED IPOPT OPTIMIZATION”和” ERROR OCCURRED DURING IPOPT OPTIMIZATION WITH WARM START “字样,同时不输出其它结果,继续进行下一步优化计算。
第一步降温前平衡液相的组成信息如下:
S | Z | n(mol) | m(mol/kgw) | x | a | r |
H2O(aq) | +0 | 5.551000e+01 | 5.551000e+01 | 8.222485e-01 | 7.837040e-01 | 9.531230e-01 |
Na+(aq) | +1 | 3.000000e+00 | 3.000000e+00 | 4.443786e-02 | 2.599406e+00 | 8.664797e-01 |
K+aq) | +1 | 3.000000e+00 | 3.000000e+00 | 4.443786e-02 | 1.650184e+00 | 5.500684e-01 |
Cl–(aq) | -1 | 6.000000e+00 | 6.000000e+00 | 8.887573e-02 | 4.732940e+00 | 7.888332e-01 |
其中S表示物相,Z表示电荷,n表示物质的量,m表示浓度,x表示摩尔分数,a表示活度,r表示活度系数。
第一步降温前平衡固相的组成信息如下:
S | Z | n(mol) |
SOLID 1,H2O(cr,Ih) | +0 | 0.000000e+00 |
SOLID 6,NaCl(cr) | +0 | 0.000000e+00 |
SOLID 7,NaCl*2H2O(cr) | +0 | 0.000000e+00 |
SOLID 8,KCl(cr) | +0 | 0.000000e+00 |
其中S表示物相,Z表示电荷,n表示物质的量。在平衡固相结果输出行列,共有43个固相及相应的含量。
用户可以依据输入的物质元素查找相应的固相及含量。如上表所示,四种相应的固相均为零,说明第一步降温前无固相析出。如果想要直接得到卤水冷冻的最终结果,可直接浏览最后一步的输出结果。
- 案例四:变质作用单点计算
1、访问ISLEC-Web在线计算服务网站:https://www.islec.net
2、点击LOGIN(用户登录),输入ISLEC-Web授权账号和密码进行登录;选择LOGIN(用户登录)菜单下MY ISLEC(我的ISLEC)进入ISLEC在线计算功能选择。
3、选择卤水变质作用计算(Brine/Mineral Metamorphism Calculation)模块,进入计算页面。
4、输入计算的温度和压力参数,默认温度为298.15K,压强为101.325KPa。这里以默认温度和压强为例,温度和压强的适用范围参考用户手册相应内容。
5、体系默认的总电荷为零,这里若无特殊计算需要不能更改。每次输入完成后应检查电荷是否平衡,若电荷不平衡,会有提示字样出现或无法正常计算。
6、输入体系物分组成,包含水以及各种离子的物质的量。输入水的物质的量为55.51 mol,Na+的物质的量为10 mol,K+的物质的量为1 mol,Cl−的物质的量为11 mol。输入完成后,点击Calculate进行计算。
7、计算的输出结果包括是否优化成功,平衡液相中各组分的含量、活度系数和活度,平衡固相的组成成分和含量。如上例,输出的结果如下:
优化结果 | SUCCESSED IPOPT OPTIMIZATION
SUCCESSED IPOPT OPTIMIZATION |
如果优化失败会输出” ERROR OCCURRED DURING IPOPT OPTIMIZATION”字样,同时不输出其它结果。
平衡液相的组成信息如下:
S | Z | n(mol) | m(mol/kgw) | x | a | r |
H2O(aq) | +0 | 5.551000e+01 | 5.551000e+01 | 8.068715e-01 | 7.381805e-01 | 9.148676e-01 |
Na+(aq) | +1 | 5.643291e+00 | 5.643291e+00 | 8.202865e-02 | 5.798319e+00 | 1.027484e+00 |
K+(aq) | +1 | 1.000000e+00 | 1.000000e+00 | 1.453561e-02 | 5.063719e-01 | 5.063782e-01 |
Cl–(aq) | -1 | 6.643291e+00 | 6.643291e+00 | 9.656426e-02 | 6.608014e+00 | 9.947025e-01 |
其中S表示物相,Z表示电荷,n表示物质的量,m表示浓度,x表示摩尔分数,a表示活度,r表示活度系数。
平衡固相的组成信息如下:
S | Z | n(mol) |
SOLID 1,H2O(cr,Ih) | +0 | 0.000000e+00 |
SOLID 6,NaCl(cr) | +0 | 4.356709e+00 |
SOLID 7,NaCl*2H2O(cr) | +0 | 0.000000e+00 |
SOLID 8,KCl(cr) | +0 | 0.000000e+00 |
其中S表示物相,Z表示电荷,n表示物质的量。在平衡固相结果输出行列,共有54个固相及相应的含量。
用户可以依据输入的物质元素查找相应的固相及含量。如上表所示,析出的固相为无水氯化钠,物质的量为4.356709 mol。
5. 常见问题与反馈
5.1 初始水量不合理
当进入各计算模块界面时,体系温度、压力、总电荷以及各组分摩尔数均有默认值。默认温度为298.15 K,默认压力为101.325 kPa,默认体系总电荷为0。默认各物种的摩尔数为0。当用户忘记输入体系中实际水的摩尔数时,默认值0将会被发送至服务器端调用ISLEC软件进行计算,而只有盐没有水的情况被ISLEC认为是不合理的。所以会提示错误。因此,只要给定的水量大于0,这个问题就会被避免。
5.2 体系电荷不平衡
ISLEC-Web默认用户给定体系处于电中性状态,体系总电荷为0。因此,当用户输入的体系所有物种的总电荷不为0时,则提示错误信息。这时,用户需要重新核验体系是否严格满足电中性。有时数据的有效数字和进位方式需要特别注意。例如,当输入的体系各物种摩尔数来源于excel计算结果时,由于有效数字保留形式的不同,可能看上去计算出的总电荷为0,但实际存在计算误差。
当遇到类似问题时,可以将由excel计算出的结果复制到文本文件中进行电中性检验,以消除excel有效数字保留导致的微小电荷误差。
另外,由于实验误差的存在,用户从实验获得的样品分析结果通常也难以严格的满足电中性。这时,需要用户调整某一组分的量,以使得给定体系严格满足电中性。通常,这个被选定的调整组分应该是体系中存在的量最多的组分。
5.3 优化过程错误
优化过程错误是指ISLEC所采用的Gibbs自由能最小化优化算法无法获得问题的解,而给出的错误提示信息“ERROR OCCURRED DURING IPOPT OPTIMIZATION.”。该错误的产生主要来源于两方面:一是最优化算法本身;二是体系中水溶液相的量相对于固相过少。如果问题来源于前者,则可以尝试略微改变体系中水的量。如果问题来源于后者,则表明ISLEC软件在该高浓度条件下热力学函数产生扭曲,无法求解。这时需要考虑增加体系中水的含量。
ISLEC软件无法处理纯固相平衡,而是始终假定水溶液相的存在。因此,当水溶液相的量极少时(如共结点或共晶点附近),优化过程错误将频繁出现。
该错误信息在单点平衡计算(Single Point Equilibration Calculation)模块不会出现,而只是出现在等温蒸发计算、连续冷冻计算以及变质作用计算三个模块中。
5.4 计算结果返回空白页面
该错误见于单点平衡计算(Single Point Equilibration Calculation)模块。其产生原因与优化过程错误相同。该模块添加了后处理脚本以选出所有析出固相。当由于优化过程错误导致返回出错信息时,后处理脚本无法获得有效的计算结果,因此导致空白页面的出现。这时需要考虑体系是否过浓,或在确保电中性条件的前提下略微调整某些组分的量(通常可以改变水的量)。
5.5 问题反馈
如果用户在使用ISLEC-Web在线计算服务时遇到任何问题请与ddong_li@hotmail.com或+86 13007766878联系,我们会在最短时间内就问题进行回复、做出解释、并尽力对ISLEC-Web在线计算服务质量进行改进。
6. 权利与免责声明
本文档中的全部内容由ISLEC开发团队提供,任何人不允许以商业性目的复制、分发和传播该文档。在同意使用条款的前提下,限于教育和科学研究的复制、分发和传播不受限制。我们保留对侵犯版权行为诉诸法律的权利。
由于热力学模拟结果解读的经验性和背景相关性,对于ISLEC-Web所产生结果的可靠性和合理解读由用户自身负责,ISLEC-Web开发团队不会对其及其产生的结果,包括其中的错误、准确性、完整性和对于特定问题的可用性承担任何负责。ISLEC-Web开发团队对于任何直接或间接由ISLEC程序导致的损失不附有法律责任,也不对ISLEC提供的信息的可靠性和完整性提供保障。
7. 致谢
感谢中国科学院(“百人计划”和“西部之光”项目)、国家自然科学基金项目(U1707602, U1407131, U1607103 and 41703064)给予的资助。感谢雷锦顺同志参与整理和编写了用户培训案例。