2011多组分精馏-简洁计算.ppt

②提馏段： 与精馏段的比流量公式的推导相同，可以得到： 相平衡关系： 物料平衡关系： （蒸出因子） 为提馏段用比流量值关联的物料平衡关系方程式和相平衡相结合的方程式。 可一直向上算，一直计算到进料板上面一块板lM+2 /d为止。 算塔釜：上式变为： 已知： L1 , V1 , T1是开始假设 的初值，可算出s1，进而算出l2 /w的值。 对于第二块板： L2 , V2 , T2同样是假设的初值， 由算出了l2 /w ，进而算出l3 /w的值。 (3) 计算各组分的wi /di值 3、比流量法（Thiese-Geddes法） 由塔顶向下，由塔釜向上均可算至第M+2板： F F=FV+FL FV FL VM+1 VM+1 LM+2 LM+2 M+1 M+2 ①若饱和蒸汽进料:FL=0 则:LM+2=LM+2 即: lM+2=lM+2 则: ②若饱和液体进料:FV=0 则:VM+1=VM+1 即: vM+2=vM+2 则: 由式 可知： 由式 可知： 代入①式： ① F F=FV+FL FV FL VM+1 +1 VM+1 LM+2 LM+2 M+1 M+2 ③对于汽液混合进料: 则有: F F=FV+FL FV FL VM+1 VM+1 LM+2 LM+2 M+1 M+2 推导过程： 对任意组分： f=d+w 根据物料平衡： vM+1=lM+2+d 由图可知: ① 代入②式： ② ③ 将①代入式③ ： 将上式解出w/d的显函数形式： 将下式代入上式 因为： 当饱和液体进料时： 当饱和蒸汽进料时： 4、θ法加速契合 通过比流量法可以得到上面推导出的三个方程，求出w/d,令其为c 。根据物料平衡： w/d = c f = w+d 可解出： 似乎我们的计算基本可以结束了！！！？？？ 计算到这里还远没有结束！！！！ 吸收因子：Ani，蒸出因子：Smi Ln , Vn , tn Lm , Vm , tm 它们分别是 的函数 如何判断它们的假设是正确的？？？ 满足此二等式，则假设是正确的 喂！！！修正的工作量相当的大啊！ 1959年，lyster等人提出了θ法加速契合的方法，该法不是直接去修正初设值，而是修正计算出的还不正确的w/d值。 原理:将上一次计算出的还不正确的w/d值乘一系数θ 进行修正。 d’为修正后的值。 同时 θ的函数形式 为θ的隐函数形式，难解析，需要试差法。 试分析其图形： 令： θ g(θ) W -D 0 θ所求解 θ =0为初值 当θ =0时， g(θ)= F -D=W 当θ =∞时， g(θ)= -D g(θ)随θ单调递减,可用牛顿切线法来试差。试差一般从θ=0开始。 5、计算新的温度分布和流量分布（t,V,L断面） 求出θ 求出比较上次正确的d’,w’ 求出新的温度及流量断面 (1) 由di’ 和wi’求出各板的液相组成（物料平衡） Lnxni= lni 精馏段： 提馏段： (3) 计算新的流量断面（热平衡） 求出各板的xni和xmi，可用试差法求各板的泡点温度tn和tm，即为各板新的断面温度。 (2) 计算新的温度断面（相平衡） 精馏段：进入=流出 n VN HN qc Vn Hn Ln+1 hn+1 n+1 VnHn=Ln+1hn+1+DhD+qc D ① 将Ln+1=Vn-D代入①式中 提馏段：与精馏段的处理方法相同，可得： qr=DhD+WhW+qc-FHF 其中： 用上式从塔顶逐板一直计算到进料板，即得精馏段的新的流量断面Vn , Ln。 用上式从塔釜逐板向上一直计算到进料板，即得精 馏段的新的流量断面Vn,Ln 。 计算过程中常会发生震荡现象怎么办啊？？？ Lyster提出很好的建议啦！！ (1)取新旧两次tn的平均值作为下次试差的温度分布值。 (2)最初两轮的计算采用相同的Vn , Ln分布。 (3)第3至第7轮试差，Vn的调整幅度不超过20%，第8至第12 轮小于10%，以后各轮则小于2.5%。 (4)采取上述约束措施后，一般经过7-8次试差便能达到收敛。 收敛判据为： 允许误差规定：D指定和D计算在有效位数第四位上相差不超过5。 6、全塔收敛判据 当全塔完全契合时，即收敛时，应该有di = di’, θ=1,当然，因为前一次的di 与后一次的di 相同，那么前一次的t,V,L断面与后一次的相同，一般7-8次即可收敛。 7、比流量法计算程序 ? ? 简捷法求D,W xni xmi R N Nm JF 原始数据 求Ani Sni 否 是 输出 算wi