第五章+传质分离过程的严格模拟计算2.ppt

三、求解xij 根据Tj、Vj求出Kij 列出三对角矩阵 用Thomas 法求解xij 四、泡点法计算新的Tj 五、H方程计算Q1、Q5、新的Vj 5.2.2 流率加和法(SR) SR法的基本要点是利用组分流率概念，将ME方程化成关于组分流率的方程形式。 解ME方程求出组分流率lij或vij，组分流率加和得到Lj和Vj， 由ME方程校正总流量Lj和Vj， 然后用热平衡方程(H方程)校正温度Tj。 若流量主要取决于板上组成，而温度主要 受热量的影响， 如宽沸点范围混合物的精馏和吸收计算适宜用SR法。 * * * * * * 第五章 传质分离过程的严格模拟计算2 上次课内容 1.物料平衡方程——M方程 2.相平衡方程——E方程 3.归一化方程——S方程 4.能量衡算方程——H方程 设计变量数N(C+6)-1 ①级数（N） 1 ②进料变量数 N(C+2) ③各级压力 N ④侧线采出单元数 2(N-1) ⑤各级传热量 N 待求变量 各平衡级温度： Tj N 塔内气液相： Lj,xij,Vj,yij 2N(C+1) 合计 N(2C+3) 总变量数、设计变量、待求变量 三对角矩阵与Thomas 方法 本次课内容 5.2.1泡点法 练习 5.2.1 泡点法 泡点法(Bubble Point Method)是三对角矩阵多组分多级计算法中的一种计算方案，由于它是在迭代过程中利用三对角矩阵解出各级组成后，通过求泡点的方法确定新一轮各级温度，因而得名。 通常，为使迭代计算有稳定的收敛，泡点温度对组成不宜过分敏感，平衡常数K值范围窄的分离物系，由于其组成变化不会引起泡点温度大幅度变化，泡点法对它的计算特别有效，精馏过程多涉及此种物系。 5.2.1 泡点法（BP法）Bubble-point Method 在此情况下，各平衡级上的传质过程主要依赖于两相流体的部分汽化和部分冷凝，平衡级温度就是泡点温度或露点温度，它们主要取决于两相组成，所以用组分的摩尔分数加和式，即S-方程来检验平衡级温度Tj是否正确。 精馏系统内的热量传递主要由潜热的变化引起，由此也引起两相流率的变化，所以用热量衡算方程，即H-方程来检验流率Vj是否正确。 规定： 进料：Fj,zi,j,TFj,PFj 压力：pj 侧采：Uj,Wj 热负荷：Qj（除Q1和QN） 级数：N 回流量：L1 气相馏出量：V1 开 始 设定Tj、Vj初值 解三对角线矩阵方程，求xi,j 归一化xi,j 泡点计算，求新的Tj，Vj 计算冷凝器和再沸器的热负荷（Q1和QN） H-eq.计算新的Vj ；计算Lj） 调整Tj和Vj 结 束 yes no 规定设计变量 是否满足迭代收敛准则 ＊BP法计算框图＊ 待求变量 各平衡级温度： Tj Q1、QN和 Lj、Vj（j=2~N)、xij、yij 一、迭代变量Tj、Vj初值的给出 1、Vj: 用指定回流比（L1）、馏出量（U1）、进料量、侧线采出量，按恒摩尔流假设给出一组Vj的初值。 一、迭代变量Tj、Vj初值的给出 2、Tj: （1）塔顶： 气相采出： 液相采出： 气、液相混合： 露点温度 泡点温度 泡、露点之间的温度 （2）塔釜： 釜液泡点温度 线性内插，得到中间各级温度初值。 二、归一化 Normalization Thomas 法计算xij 由于求三对角矩阵方程时没有考虑S-eq.的约束，必须对得到的xi,j归一化。 Σ = C j i j i x x = i 1 j i x , , , 三、计算新的Tj 一种重新设定Tj的计算方法 三 用H方程求各板Vj 1.计算新一组流率之前，先要用已求出的Tj及xil,yij与求焓关系式hi(Tj)、Hj(Tj) 求出每股物流的摩尔焓hj，Hj 2.计算V2，并用H方程计算冷凝器负荷Q1 3.计算QN（再沸器的热负荷）：全塔总物料衡算式和总热量衡算式得到 4. 修正的H方程求解Vj（j=2~N) 修正的H方程如下： Lj Lj-1 方程： 修正的 1 = + \ + j j j j j V V H g b a ） （ 起） （从 36 6 3 1 1 1 1 - - = \ - - - - V j j j j j V V b a g 6-29 （6-33） 五、迭代收敛的标准 或更简单的： 小结 Vj初值：等摩尔流假设 Tj初值：顶抽出设定T1，釜液设定TN，其他线性插值 Thomas 法