甲基二氯硅烷-甲基三氯硅烷—二元体系汽液平衡的模拟计算解析.docx

计算机在化学化工中的应用结课作业专业与班级化工302 姓名 zcf 序号 成绩甲基二氯硅烷-甲基三氯硅烷 —二元体系汽液平衡的模拟计算一．问题描述利用新型泵式沸点仪，在101.3KPa下测定了甲基二氯硅烷-甲基三氯硅烷苯体系以及3个二元体系在不同液组成时的沸点；对所测的二元体系数据，采用间接法由TPx推算了与之平衡的气相组成，再用最小二乘法分别求出了二元体系的Wilson、NRTL、Margules和van Laar4种模型的最佳配偶液相活度系数模型参数，并由二元最佳配偶参数推算甲基二氯硅烷-甲基三氯硅烷二元体系的气相组成和跑点温度，与实验测得的结果进行比较。由乙炔通过催化硅氢加成反应生产甲基乙烯二氯硅烷的反应中，甲基二氯硅烷是反应物，由于存在副产物甲基三氯硅烷和未反应的甲基二氯硅烷，故在分离过程设计中，有关甲基二氯硅烷和甲基三氯硅烷的汽液平衡数据是必不可少的。表一甲基二氯硅烷-甲基三氯硅烷汽液平衡实测数数表二 用于拟合安托因常数的温度-饱和蒸汽压数据T/K P1/Pa P2/Pa316 106682.7 45845.5318 113871.9 49262.8320 121436.7 52883.2322 129390.5 56715.2324 137747.1 60767.9426 146520.6 65050.2328 155725.5 69571.7330 165376.4 74341.7332 175488.2 79370.1334 186076.1 84666.7 336 197155.5 90241.6 338 208742.2 96105.2二．解题思路1.一般中低压下（1500~2000KPa），气相作为理想气体处理时，汽液相平衡方程可以写为：(且ΣXi=1,ΣYi=1,i=1,2,3..)其中：P为体系压力，Xi‘Yi为气相摩尔分数，γi为活度系数，Ps为饱和蒸汽压根据安托因方程：其中A,B,C均为安托因常数，T为温度由“纯物质化学性质查询软件”可以查得拟合ABC所需要的数据即表2并运用Matlab软件通过运行包含三（1）和（2）中的语句即可获得环己醇与环己酮的安托因常数,便可以求得饱和蒸汽压Ps.此处会用到的拟合函数lsqnonlin的语法。2.求解活度系数，即算出Wilson活度系数模型的配偶参数，Wilson方程如下：由变成计算三（3）拟合出Wilson活度系数模型的配偶参数，代入方程即可求得活度系数，再将求得的Ps,活度系数，以及体系压力P，物质的液相实验值即可求得汽相组成计算值Yic与相对偏差。根据由上即已求出所有量，作图即可。在Matlab中主要函数的说明：A.function ：创建函数命令。定义的一般格式：function [y1,y2,…,yn] = FuncName(x1,x2,…,xn) % 函数声明语句y1 = … % （表达式 1 ）y2 = … % （表达式 2 ）…yn = … % （表达式 n ）其中， 输入参数为 x1，x2，…，xn, 输出参数为y1，y2，…，yn 。 各参数可以是标量、 向量或矩阵。若函数只有一个返回结果，声明语句可以写为：function y = FuncName(x1,x2,…,xn)若没有返回结果，声明语句写为：function FuncNam