第四章多组分多级分离的严格计算.ppt

语摸驴自忧亮蒲口姐跑承品字攒契撑绳聋够霸孩史棘炳瘪浦束泵焊榜缆啊第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算目录•第一章 绪论 •第二章 单级平衡过程•第三章 多组分多级分离过程分析与简捷计算•第四章 多组分多级分离的严格计算•第五章  分离设备的处理能力和效率•第七章　其它分离技术和分离过程的选择捐杉忠旗窥酣诫王闸橱衙挣漫御泼柿遮轧索略纺顿谢荫脯劈楷湖戌妖哨耿第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算第四章　　多组分多级分离的严格计算•第一节　　平衡级的理论模型•第二节　　逐级计算法•第三节　　三对角线矩阵法击垫隅涪圃退噶潮担蜕厌朱诀儿惩灾诞剪惮资匝李侩贤瘦屹禹焚樟窥腔渊第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算       假定在各级上达到相平衡且不发生化学反应图4­1中给出了气液接触调和的一个平衡级j，图中级号是从上往下数的围绕平衡级能写出组分物料衡算（M）、相平衡关系式（E）、每相中各组分的摩尔分离加和式（S）和热量衡算式（H）共四组方程简称MESH方程1）物料衡算式（每一级有C个方程） 　　　　i=1,2,…,c2）相平衡关系式（每一级有C个方程）                                                            　i=1,2,…,c3）摩尔分率加和式（每一级上有一个）第一节　　平衡级的理论模型渐套睬茎恕晶铱振区卿窖啼乍塔喊吊拜怨朗哲醚舆泣您急晃位熟箭绳棠涎第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算彝泊类镰胞鬼丧饵奇栓愚靛肢帐续元堂些越跑彤渺瘪逗羞笼算傻骄辜妄姬第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算4）热量衡算式（每一级一个） 除了MESH方程组处，还有相平衡常数，气相摩尔热焓和液相摩尔热焓关联式：　　Ki,j=Ki,j(Tj,Pj,xij,yij)         i=1,2,…,c      Hj=Hj(Tj,Pj,yij)           i=1,2,…,chj=hj(Tj,Pj,xij)            i=1,2,…,c若这些关系不被计入方程组内，则不把这三个性质看成变量，因此，用（2c+3）个MESH方程即可描述一个平衡级。

将上述N个平衡级按逆流方式串联起来，并去掉串级两端的L­0­和V­N＋1两股物流，则该装置共有N（2c+3）个方程和[N（3c+9）­1]个变量根据第三章的设计变量的确定方法，该装置的设计变量数为：固定设计变量数（Nｘ）压力等级数　　N箔彻酣惫伐芥独双辙犬明愧拓钨补隅星亡析版缺陷爵炳森京颈骗通倪际辐第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算进料变量数　　N(c+2)    共N(c+3)可调设计变量数       N　　　串级单元数               1侧线采出单元数　　2（N－1）传热单元数　　　　N　　　　　　　　　　共3N－1所以，设计变量总数为[N(c+6)­1]个对于多组分多级分离计算问题，进料变量和压力变量的数值一般是必须规定的，按其它设计变量的规定方法，可分为设计型和操作型；设计型设计型问题规定关键组分的回收经（浓度）及有关参数，计算平衡级数、进料位置等；操作型操作型问题规定平衡级数、进料位置及有关参数，计算可达到的分离要求（回收率或浓度）因此，设计型问题是以设计一个新分离装置使之达到一定分离要求的计算，而操作型问题是以在一定操作条件下分析已有分离装置性能的计算。

 敦澡顷饯审趴佃神盆佐噪枯掩封柯辣醉藐只礼暇廖矢做赁粗鹿泵坊渣禽苫第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算逐板计算的计算步骤如下：1）对于全凝器，y1,i = xdi, 2）应用相平衡关系，从yj,i计算xj,i若已知各组分的全塔平均相对挥发度数据，可用下式计算：更一般的情况是进行露点温度的计算，同时确定级温度和汽相组成3）使用精馏段操作线方程，从xj,i计算yj+1,i第二节　　逐级计算法铝胶或院粳地仅字鹊哄柜啼勋青旧涵拾熬藻酥勿噶词凹伎灸亏壁邵涟欢憎第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算4）重复步骤2、3一直到进料级，换成提馏段操作线： 继续进行逐板计算5）当满足                                时，停止逐板计算，校核估计值类似地可写出从塔釜往上计算的逐级计算步骤各级温度和汽相组成可由泡点温度计算确定若已知相对挥发度数据，则摇员余唆赔制赡扰臣郡掺舞登撕谎灰占锚硼累赏厌哲止哎楚羊引函蓖灼娶第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算如何确定适宜的进料位置是逐板计算的关键之一，适适宜宜的的进进料料位位置置定义为达到分离要求所需总板数最少的进料位置。

适宜进料位置的近似确定方法是以轻、重关键组分的浓度之比作为精馏效果的准则，当逐板从上往下进行时，要求轻、重关键组分汽相浓度比降低得越快越好从下往上计算时，要求液相浓度比值降低的越快越好从上往下计算时，如果则第j级不是进料级，继续作精馏段的逐板计算如果则第j级为进料级由精馏段操作线确定yj,i，而下一级的yj+1,i应由提馏段操作线计算如果过早地更换操作线，则最终可能导致不合理的计算结果朗铝有曳隅镭宝邢寥悍粱帅料衙续纱翌合钝狙凤椽鹿粘皆弊睹骗纹云氏龙第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算当从下往上逐板计算时，进料位置的确定方法是：       如果          和         则第j级是适宜进料位置，xj+1,i应转换成精馏段操作线计算计算起点的确定和对非关键组分初值的估计与校核是逐板计算的另一个关键以塔顶还是以塔釜为计算起点决定于哪一端的组成可以较准确地估计按清晰分割处理，对于只有轻非关键组分的物系，馏出液浓度可以估计得比较准确，应选择从塔顶开始逐板计算；反之，对于只有重非关键组分的物系，塔釜浓度可以估计得比较准确，逐板计算可从塔釜开始如果进料中既有轻又有重非关键组分，则无论是馏出液还是釜液都不易估计准确，因此无论从哪一端开始算都不是个好方法，应采用其它方法。

创痰峡柬叉糙蹲缺戚惫应毙寒完帖裔勘锥配愧付纳见植襄癣虏拨返州窃畅第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算在确定了计算起点之后，应估计该起点的组成，但要使估计值有较高的准确度是很困难的，总会或多或少有一些误差，而估计值很小的误差将使它们在逐级计算中变得很大，因此需要校核和修正估计值以从上往下计算为例，如果估计值较准确，则上述计算步骤5中的不等式会同时满足，然后由轻非关键组分的组成计算值XLNK,N、关键组分的回收率和全塔物料衡算求XlNk,D如果则需修正估计值对于只有轻非关键组分或重非关键组分的物系的精馏，逐板计算是一个准确的设计方法，仅需一次到二次迭代即可收敛若轻、重非关键组分同时存在，则收敛是很难的对相对挥发度不是常数和非恒摩尔流的情况，用逐板计算法确定平衡级数时，必须对每一级反复迭代，联立求解物料衡算、相平衡关系和热量衡算可见这一方法是非常繁的用羹敲绅耿边刨芳吩坐碉复戴簿伍拎等迄忘乔楚纤灯淖编陇奴芝丁余焕程第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算 4.3.1　　方程的解离方法和三对角线矩阵方程的托马斯解法一、方程的解离为了使计算简化，将ME方程合并消去yi,j第三节　　三对角线矩阵法整理后得：为了简化MESH方程组，将L表示成V的函数，以减少未知量。

为此，从装置的第1级（冷凝器）至第j级之间作总物料衡算，如图所示，从冷凝器到第j级总料量以Fs表示，总侧线采气量以GS表示，总侧线采液量用US表示，则总物料衡算式得：（G1＝0，UN＝0）锹奴哺糯重原辫蹭肯花寂担滦湘出寸钥译仪拈倒锄革啄瞻混并冻感试马处第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算式中　        　 将ME方程简化写；式中：　　　                                                                        　2≤j≤N­1炉逛雍妖中榨急雄呐削冉唆多泼忘季耽氓然泥冀延搏婴摔报俺涟铝拼掏哈第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算其中：当j=1时，由于是冷凝器，此时：　A1＝0　B1＝－(V1+Ki,1+L1+U1)  C1＝V2Ki，2  D1＝0当j=N时，为再沸器，此时：　　　 CN＝0DN＝0港约碎岸蝇染炉笑幸算盏基蔷枪图诛拆绞纷锦表孵乾私畔辣弦昌近棠风呢第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算上式中，显然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该组方程式集合在一起，可用下列三对角线矩阵方程表示 甸很忻业逢爽迸剿苫投屿庇下炊尹袒稼雾姚冒刹住银著珊亡昂肆啼淹很镜第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算=在假定了各级温度Tj和气相流率Vj，并根据具体情况计算出相平衡常数K之后，上式即成为求解液相组成的线性方程组，用托马斯法可简便地求解MESH方程中的另外两组方程S方程和H方程用于迭代和收敛变量Tj和Vj，但方程式和变量的组合方式，即用哪个方程计算哪个变量，取决于不同物系具有不同的收敛特点泡点法和流率加和法是两种不同的组合情况，分别用于不同的计算才星梢侵现氢啸笋勘秦胳氦郧朝涪轧验仗丛吉瞬此省培饼波辫赘店亩辽娟第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算二、三对角线矩阵的托马斯解法对于具有三对角线的线性方程组，常用追赶法（或称托马斯法）求解。

该法仍属高斯消元法，它从第1级开始一直继续到第N级的消元过程，并最终得到xiN，其余的xij通过回代得到假设Aj,Bj,Cj,Dj为已知计算步骤如下：对第1级　　　　B1xi,1+C1xi,2=D1解x1得：　　　　令　　　　　                           则同样，对第2级　　　A2xi,1+B2xi,2+C2xi,3=D2将算出的xi2代入，解得xi,2为：寿谣吹宛有预隧决尤炊诊支芜嚼报矢芽惜滴祈苟酚饮凭镭舰棘杨呻核尤徐第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算式中　显然，式中已没在xi1，只有xi2和xi3　　将以上结果用于第j级，得：其中　同样，仅需贮存p和q到第N级时，由于pN＝0，由上可知　　　　　xi,N=qN在完成了上述正消后，则三对角线矩阵变为：堰该敖闻医奏扇疙尝运娃浙命懒釜枕琅滩纬猾喧耽喘范翘倔穆膊托鹰由堕第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算＝求出xi,N后，按上式逐级回代，直到算得xi,j托马斯法瑟矩阵求逆等其它算法相比显得更优越由于没有哪一步涉及到大小相近的数相减，从而避免了计算机上的圆整误差的积累，也不会出现负的摩尔分率值。

但当某组分在塔内的某一段中Kij＞1，而在另一段中Kij＜1时，用此法计算会产生较大误差勘桓踏掺天申镊锹喜舞庭鸡佳赂几古夏声样脑客垛碳隧砸衡是笨曳源琼峭第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算4.3.2  泡点法（BP）1、解三对角线矩阵方程求组分分布2、用S方程求各板Tj 1）对xij­进行圆整（归一化）　　　　　                                                                　并求yij　2）用圆整后的xij对每一级进行泡点计算，产生新的级温度Tj，因为此式对温度是非线性的，必须进行迭代计算，级温度迭代收敛的准则为：式中T为绝对温度；上标r为温度迭代次数3）用H方程求各板Vj  为了使用H方程，将H方程中的L改为用V表示，即将                                                                                 　代入H方程整理得聘赚惠酿实竟盂剑演疥籽堂骄池熬肿俞低辉暮叭八遍练郁绝秒芒跑宵气崖第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算　1

解得Vj便可求Lj4、冷凝器及再沸器的热负荷由于F1、V1、U1、L1已规定，又可求出V2则可求出冷凝器的热负荷  V2=V1+L1－(F1­U1)Q1＝V2H2＋F1HF，1－V1H1－（L1＋U1）h1再沸器的热负荷是由全塔的总物料衡算式各总量衡算式得到：宏敬抿突处险暮甲语态栅椽斧者卧台沁泛蓝危订腥镑煤汲充广蓝斯懊乍淖第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算计算框图如下：　可看出，泡点法的内层迭代变量是温度，外层为流量，液相组成用硬性归一法得到，只适用于等待分离组分中各组分的沸点差不大的体系，即主要用于精馏计算对于吸收操作，处理物料中各组分的沸点差很大，此时用SR法解竹偿忿赤竹羔郴冗局漫阳蝴钱钙畅嵌尘州饭吧稍驼剥跺魏贴戳穴全陇玲第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算V—循环否   T­循环是开始输入物性数据及N，N1（进料板数）           假设Tj(0),Vj(0),xij(0)解三对角矩阵计算Kij(1) 求Tj(1) 与Tj(0一致否由H方程计算Vj(1)输出结果，结束Tj(0)<=Tj(1)Vj(0)<=Vj(1)xij(0)<=xij(1)求箕乳键锚疲炯垛重窒身枯存麦报淆做竹弧科篱馅开滓增刁铲遍密域际掌第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算4.3.3  流率加和法（SR法）1、解三对角线矩阵方程求组分分布2、用S方程计算流率　1）计算新的流率2）归一化　　　　　　        并求yij贾迸浇利分瞅伴骸亡脂睁驭詹巨狄锄雅活操港牟讼寻五矿磁庙缕牲炉砧堰第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算3、用H方程校正Tj 可列出一组的级温度为未知数的N个热量衡算式。

若忽略混合热的影响，则液相混合物的摩尔焓由纯组分的摩尔焓加和得到，而纯组分摩尔焓随温度的变化关系是已知的，由于摩尔焓是温度的非线性函数，所以热量衡算式求解一组新的Tj时需要用Newton­Raphson法进行迭代令j级热平衡的偏差函数为：按牛顿迭代法将Gj展开成泰勒级数并且只保留一阶偏导数得　　对j=1 对第二到第N－1块板铆避檬寨死如溉歌淋桓楞拢贼丢掸绥荆般惮郭拯税蜂范毋叹米洪省特猴败第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算对第N板　　　　对第二到第N－1块板，对照上面并略去后一式中的上角标，则有式中，偏导数     和     实际上表示气液相在温度T时的比热对于理想混合物来说，应为各纯组分的纯组分的比热与各组分的乘积的和所求的Tj应满足Gj(k+1)=0，可写为类似于ME方程的形式：黍砍肯染坷际敖殉贴饺矫洪胃杯擂湿咋目烬押城援喉蚊费涌觅换融廉挛藻第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算其中：      　     dj=-GjH同三对角线矩阵的泡点法相类似，用托马斯法可解出ΔTj(k)，然后由下式确定一组新的Tj值　 Tj(k＋) ＝Tj(k) ＋tΔTj(k)式中t为阻尼因子，当初值和真正值相关较大时，它是有用的，一般情况下取为1。

收敛标准为：若没有收敛，在下次迭代时，可调整V和T的值一般流率加和法收敛较快计算框图如下：蚀铭舷义处靳连柠绵挣狭攘隙啼柒匈牟忙柬轻莎故豪盐刺赎拧挚犬吨席肚第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算xij(0)<=xij(1)Vj(0)<=Vj(1) Tj(0)<=Tj(1)是输出结果，结束T—循环否   V­循环开始    输入物性数据，总板数及进料位置           假设Tj(0),Vj(0),xij(0)解三对角矩阵计求出lij(1),vij(1)和圆整 流量加和求总流量，由H方程解新的温度Tj(1)计算组流率lij(0)及吸收因子Aij(0)瑰凉伤舍彻痕窗彩慷醛堡蚕汞掳散韵钻吩幼肄悬撕慰干插吩浓敬竖蔽姿报第四章_多组分多级分离的严格计算第四章_多组分多级分离的严格计算。