交换塔设计.DOC

郎五可 化工学院 1011207111 分离B10同位素化学交换精馏塔理论板数计算 一、背景 天然硼其中重要 10B对中子具有非常强的吸收能力塔顶回流（络合产生液相） 交换塔 在B10同位素更倾向于形成络合物，在反复多次接触后，B10在液相富集。从而得到高丰度B10 塔底回流（裂解产生气相） 2.模型假设 过程处于稳态，恒摩尔流； 交换仅发生在液相； 理论板上气-液相处于化学平衡； 恒温操作。 液相全部是1:1络合物 苯甲醚没有损耗 三、建立理论级稳态模型 1.模型图 图1 总流程简图 2. 建立模型 全塔模型如图1；由于B10在塔底富集所以进料以下称富集段；以上称提取或剥淡段。富集段中第j块理论板i组分的模型如图2。 剥淡段第j块理论板的模型如图3。 全塔 （1） （2） 顶回流: （3） （4） 剥淡段第j板： （5） 富集段第j板： （6） 对底回流： （7） ， （8） 分离系数方程: （9） 化学交换反应平衡方程：（10） 归一方程：（11）， （12） 3.模型化简 对于二组分体系，归一化方程可以省略，i可以省略；交换塔主体是恒温操作，K是常数，K=α，则每个板的分离系数α恒定，j可以省略。 令，塔底回流比； 塔顶回流比； 液气比 这些参数有助于指导实际操作。 已知 yf，B, yb, yt；TT（温度）, P 求 F, T, N, Rb, Rt, cb, R ;NT, NB; 方程： ； ； ； K=α 平衡方程 分离系数 上段（剥淡段） 下段（富集段） 说明：上段和下端的y有两种表达方式，前面是本模型的逐板浓度表示；另一种是对整个上段的组分做物料衡算得到，相当于精馏塔的操作线方程。后者可以用来确定最小回流比。在进行程序实现时，这两种方程都能达到结果。 四、程序实现及实例计算 1.计算框图 2.程序实现 函数M 文件 1)ban.m function x=ban(a,y) x=a*y/((a-1)*y+1); end 2)rminb.m function R=rminb(yb,yf,a) R=(yb-ban(a,yf))/(ban(a,yf)-yf); end 3)upp.m function y=upp(Rt,yj,xj,xjj) y=yj+Rt/(Rt+1)*(xj-xjj); end 4)downer.m function y=downer(Rb,yj,xj,xjj) y=yj+(Rb+1)/Rb*(xj-xjj); end 主程序M文件，inicial1.m clear all b=input(b) yf=input(yf) yb=input(yb) yt=input(yt) tt=input(tt) TT=273+tt a=exp(24.1867/TT-0.0504) B=b*1000/67.81 T=B*(yb-yf)/(yf-yt) F=B+T cb=B*yb/(F*yf) RBMIN=rminb(yb,yf,a) RTMIN=(RBMIN+1)*B/T Rtmin=rmint(yf,yt,a); Rb=1.5*RBMIN L=(Rb+1)*B V=B*Rb Rt=L/T R=Rt*(Rb+1)/(Rt+Rb+1) R1=L/F; y(1)=yt y(2)=yt x(1)=yt j=1 while y(j)yf j=j+1; x(j)=ban(a,y(j)); y(j+1)=upp(Rt,y(j),x(j),x(j-1)); end NT=j j=j while y(j)yb j=j+1; x(j)=feval(ban,a,y(j)); y(j+1)=downer(Rb,y(j),x(j),x(j-1)); end N=j+1 NB=N-NT disp( x1 y2 x1 y1), CC=[x(:,1:N-1), y(:,2:N),x(:,1:N-1),y(:,1:N-1)]; plot(x(:,1:N-1),