第4章多组分多级分离的严格计算.pptVIP

化工分离工程;目录;第四章　　多组分多级分离的严格计算; 假定在各级上达到相平衡且不发生化学反应。图4-1中给出了气液接触调和的一个平衡级j，图中级号是从上往下数的。围绕平衡级能写出组分物料衡算（M）、相平衡关系式（E）、每相中各组分的摩尔分离加和式（S）和热量衡算式（H）共四组方程。简称MESH方程。 1）物料衡算式（每一级有C个方程） ;;4）热量衡算式（每一级一个） ;;逐板计算的计算步骤如下： 1）对于全凝器，y1,i = xdi, 2）应用相平衡关系，从yj,i计算xj,i。若已知各组分的全塔平均相对挥发度数据，可用下式计算： ;4）重复步骤2、3一直到进料级，换成提馏段操作线： ;如何确定适宜的进料位置是逐板计算的关键之一，适宜的进料位置定义为达到分离要求所需总板数最少的进料位置。适宜进料位置的近似确定方法是以轻、重关键组分的浓度之比作为精馏效果的准则，当逐板从上往下进行时，要求轻、重关键组分汽相浓度比降低得越快越好。从下往上计算时，要求液相浓度比值降低的越快越好。 从上往下计算时，如果 ;当从下往上逐板计算时，进料位置的确定方法是： 如果 ;在确定了计算起点之后，应估计该起点的组成，但要使估计值有较高的准确度是很困难的，总会或多或少有一些误差，而估计值很小的误差将使它们在逐级计算中变得很大，因此需要校核和修正估计值。 以从上往下计算为例，如果估计值较准确，则上述计算步骤5中的不等式会同时满足，然后由轻非关键组分的组成计算值XLNK,N、关键组分的回收率和全塔物料衡算求XlNk,D。如果 ;? 4.3.1　　方程的解离方法和三对角线矩阵方程的托马斯解法 一、方程的解离 为了使计算简化，将ME方程合并消去yi,j ;式中　 　 ;其中：当j=1时，由于是冷凝器，此时： 　A1＝0 　B1＝－(V1+Ki,1+L1+U1) C1＝V2Ki，2 D1＝0 当j=N时，为再沸器，此时： 　 　 　 ;上式中，显然Bj、Cj、Dj应与组分i有关是，只是为了方便起见将下标i略去。 于是，ME方程可表示为： 第1级　　B1xi,1+C1xi,2=D1 第2级 A2xi,1+B2xi,2+C2xi,3=D2 ? 第j级 Ajxi,j-1+Bjxi,j+Cjxi,j+1=Dj ? 第N-1级 AN-1xi,N-2+BN-1xi,N-1+CN-1xi,N=DN-1 第N级 ANxi,N-1+BNxi,N=DN 该组方程式集合在一起，可用下列三对角线矩阵方程表示 ;=;二、三对角线矩阵的托马斯解法 对于具有三对角线的线性方程组，常用追赶法（或称托马斯法）求解。该法仍属高斯消元法，它从第1级开始一直继续到第N级的消元过程，并最终得到xiN，其余的xij通过回代得到。 假设Aj,Bj,Cj,Dj为已知。计算步骤如下： 对第1级　　　　B1xi,1+C1xi,2=D1 解x1得：　　　　 ;式中　 ;＝ ;4.3.2 泡点法（BP） 1、解三对角线矩阵方程求组分分布 2、用S方程求各板Tj 1）对xij-进行圆整（归一化） 　　　　　 　并求yij ;　1jN ;计算框图如下： 　可看出，泡点法的内层迭代变量是温度，外层为流量，液相组成用硬性归一法得到，只适用于等待分离组分中各组分的沸点差不大的体系，即主要用于精馏计算。对于吸收操作，处理物料中各组分的沸点差很大，此时用SR法。 ;V—循环;4.3.3 流率加和法（SR法） 1、解三对角线矩阵方程求组分分布 2、用S方程计算流率 　1）计算新的流率 ;3、用H方程校正Tj 可列出一组的级温度为未知数的N个热量衡算式。若忽略混合热的影响，则液相混合物的摩尔焓由纯组分的摩尔焓加和得到，而纯组分摩尔焓随温度的变化关系是已知的，由于摩尔焓是温度的非线性函数，所以热量衡算式求解一组新的Tj时需要用Newton-Raphson法进行迭代。 令j级热平衡的偏差函数为： ;对第N板 　　　　 ;其中： ;xij(0)=xij(1) Vj(0)=Vj(1) Tj(0)=Tj(1)