精选固定床反应器的工艺设计
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
CA = PFyo(1-x)/(RG*T(6.5555Fyo+ 急冷量)) CB = P(Fyo(5-x)+0.9急冷量)/(RG*T (6.5555Fyo+急冷量))若把这些CA和CB项代入式(6-32),则: rv = - dCA/dθ = kCACB0.5 = Aexp(-E/RGT) (P/(RG*T(6.5555Fyo+ 急冷量)))1.5 Fyo(1-x) (Fyo(5-x)+0.9急冷量)0.5将上式与(6-8)合并:dx/dz = AcAexp(-E/RGT)(P/RG*T(6.5555Fyo+急冷量))1.5(1-x) (Fyo(5-x)+0.9急冷量)0.5
压降方程
开始可计算出∆P/z值,而后在每一步进行适当的压力校正。
实例,恒温反应器设计
H2C=CH2十H2 → CH3一CH3 (6-9)这是一个放热反应,但它可在很小的恒温反应器中进行。铜—氧化镁为催化剂,把这些催化剂装在一个内径为9.35mm、长为280mm带水夹套的黄铜制的反应器管中。由于把37.85 l/min的水循环通过夹套,因而可得到大于5.67kJ/(m2 S K)的水膜传热系数和恒定的管壁温度。在此反应器内,9—79℃范围内的各种恒定温度下进行的许多试验证明,对氢气来说,式(6-9)表示的反应是一级的。通过给定下列条件可写出反应速率方程:rv = kCH (6-10)该系统中, η = 1
有时反应速率方程用分压而不是用浓度表示:rvp = A exp (-E/RT)PaαPbβ 式中,P是反应物或产物的分压,bar rv和rvp之间的关系是:rvp = rv (RT)α+β 式中R = 0.08314 m3 bar/(kmol K)
对于乙烯加氢:1克分子乙烯 十 1克分子氢气 → 1克分子乙烷 (6-18) a = 1, b = 1, p = 1, q = 0
NT = (1-yoxA)NTo
又因为,
对于乙烯加氢这个反应:
设计方程变为:
将上式积分得:V = (RTFyo/(kP))[xA + ln((yo-1+yoxA)/(yo-1))] 现在用上式计算一个试验室反应器。yo = 0.49 F(总进料)=2.0×10-8 kmo1/sR(气体常数)= 0.08314 m3 bar/(kmo1.K)R(气体常数)=8.314kJ/(kmol.K)T=67℃十273.2=340.2Kk = 5.96×106exP(一55731/RT) kmol/(s. m3)希望:xA=1.0 假定P=1.0 bar,
设计方程: Fyodx = kCHdV 由于反应涉及到体积变化,对此,必须加以考虑 。例如: aA + bB → pP +qQ (6-12)设A为数量有限的反应物。开始时, NTo = NAo +NBo (6-13) NAo = yoNTo NBo = (1-yo)NTo = (1-yo)NAo/yo (6-14)在达到转化率分数xA后:NT = NA + NB + NP + NQ = NAo (1-xA) + NBo – (b/a)NAoxA + (p/a)NAoxA + (q/a)NAoxA = NTo – NToyoxA·((a+b-p-q)/(a+b))·(a+b)/a (6-15)令 d’ = (a+b-p-q)/(a+b) (6-16) 一维 Nhomakorabea流模型
绝热气体反应器最简单的模型是一维塞流模型, 为了适应此模型,反应器必须符合下列条件:·反应器和周围环境不应有热交换(绝热条件)。·通过固定床催化剂床的流动应该是在床层内任何一点反应物流的线速度都是相同的。·不应有轴向扩散。不应有径向扩散。
一维模型绝热固定床气体反应器的设计
Fyodx = rvdV = rvAcdz Fyodx(-∆H) = ∑miCpidT
发表的活化能数值为35.4kcal/g. mol,但得不到发表的频率固子A的数值,该数值取决于具体的催化剂, 现假定A=5.73 x l06,反应器入口的反应温度假定为To=600℃十273.2=873.2K,反应器压力可从35到70绝对大气压,假定约50绝对大气压或50 bar。根据芳烃加氢的一般经验,氢与芳烃之比可假定为5.0. 氢: Cp=20.786 甲烷:Cp=0.04414T十27.87 苯: Cp=0.1067T十103.18 甲苯:Cp=0.03535T十124.85 进入加氢脱烷基装置的混合氢气物流的纯度为90%,其余的为甲烷。
每一反应物或产物的克分子数是:甲苯:Fyo(1-x)苯: Fyox氢: Fyo(5-x)+0.9急冷量甲烷:Fyo(5/9+x) + 0.1急冷量 总克分子数:Fyo(1+5+5/9)+急冷量,或6.5555Fyo+ 急冷量每一组分的克分子分数为:甲苯: y1=Fyo(1-x)/(6.5555Fyo+ 急冷量)=yA苯: y2=Fyox/(6.5555Fy。十急冷量)氢: y3=(Fyo(5一x)+0.9急冷量)/(6.5555Fyo+急冷量)=yB甲烷: y4=(Fyo(5/9+x)+0.1急冷量)/(6.5555Fyo+急冷量)浓度项可写成:C = N/V = Np/(NRG*T) = p/(RG*T)因为: p = Py 所以: C = yP/(RG*T)
设计方程
dV = Acdz
进入的克分子数=出来的克分子数 十 积聚的克分子数 十 通过反应损耗的克分子数
在稳态操作时,反应器中没有物料积聚,或没有物料从反应器中放出,Fyo = Fyo(1-dx) + rvdV (6-6)整理得: Fyodx = rvdV (6-7)将(6-4)代入(6-7)得: Fyodx = rvAcdz (6-8)
式中, (-∆H)——整个反应系统的反应热; mi—每一种反应物(包括反应物和产物)的克分子流量, kmo1/s; Cpi——每一种反应物的克分子热容量,kmol/(kmo1.K)。
设计用于放热反应的固定床气体反应器
对绝热反应器的限制是由限定的出口温度t最大决定的,达到此温度标志着某些不理想的过程如副反应、选择性不良和催化剂严重结垢开始出现。如末达到所要求的转化率x,温度就已上升到限制温度,此时,可采用几个方法中的一个在温度不超过t最大的情况下来提高x使其达到所要求的值。方法1 (1)把催化剂体积分成两个或两个以上的床层,以便使每一床层都达到t最大; (2)用间接换热器冷却两相邻床层之间的气体反应物使其回到第一床层入口温度to
以空塔为基础的空塔流速um的方程为:
若催化剂的形状是直径和高度相等的圆柱体,令rp=颗粒半径,则,
若给定2rp=3mm,则dp=0.003 m。
反应混合物的粘度μ可用以下的关系式从各个反应物的粘度μ1,μ2,μ3,μ4计算出来.此方程是:
ij项可以是12、13、14、21、23、24、31、32、34.41、42、43,Mi和Mj是各组分的分子量。
现在计算每一床层(除最后床层外)末端处的急冷量:
急冷量= ∑miCpi(T-To)/(Cpi急冷(T-298.2)) Cpi 急冷 = 0.9x20.786 + 0.1x66.413 = 25.349 kJ/Kmol K 利用Ergun方程计算压降:在Ergun方程中,G为质量流率.其单位为kg/(s m2)
如果没有一个方法能妥善地解决放出的热量问题,则必须考虑气体反应在下列反应器中的一种内进行: 流化床反应器; 带有一种惰性液体的气-液相反应器 非恒温非绝热固定床气体反应器
绝热固定床气体反应器设计实例:
甲苯加氢脱烷基反应器的设计
反应速率方程为:rv = - dCA/dθ = kCACB0.5 = Aexp(-E/RGT)CACB0.5 式中, CA——甲苯浓度,kmo1/m3; CB——氢浓度,kmo1/m3; k—甲苯加氢脱烷基的反应速率常数; A—频率因子; E—活化能 148114kJ/kmo1; ∆H— - 49974kJ/kmo1转化的甲苯; θ—停留时间,s。
方法4 第四个方法是在反应器进料中加入一种惰性气或一种过量的反应物,其作用好象一个受热器一样,吸收大部分放出的热量,这样,对给定的转化率而言,就减小了温升。和其他方法相比,此法有一严重缺点,即为了容纳惰性气或过量的反应物,需要增大反应器的尺寸。 第四个方法派生出来的一个方法是加入一种惰性液体,它在反应器的通常压力下,其沸点接近反应温度,反应放出的热使惰性液体蒸发,这样,其所吸收的大量热量就满足了潜热的需要。在某些情况下,和采用纯气相方没相比,这种方法反应速率比较低。此方法的缺点是需要的反应器大而贵。
方法2 如果所用间接换热器反应温度太高的话,则可采用第二个方法,把一种冷的反应物喷入两相邻床层之间的反应物流中,把反应物冷却到To.由于补充了一种被消耗的反应物,这种直接冷却或急冷方法又有使反应朝着所要求的方向进行的优点。
方法3 第三个方法是一种老的德国分段急冷的技术,虽然目前仍在使用,但其缺点是波动很大,且难于设计。
如反应器内径为1cm,则反应器高度将为 89.1 cm。
6.2 绝热气体反应器
定义绝热反应器是一种和周围没有热交换的反应器。因此,如果放热反应是在绝热固定床反应器中进行,则气流的温度将沿着入口到出口的方向增加。另一方面,如果吸热反应是在同样的反应器中进行,则气流的温度将会沿着反应器的长度降低。
假定内径为2m
dT/dz = (dx/dz) Fyo(-∆H)/(∑miCpi) dx/dz = AcAexp(-/RGT)(P/RG*T(6.5555Fyo+急冷量))1.5(1-x) (Fyo(5-x)+0.9急冷量)0.5联立解微分方程可求出完成一定生产任务的反应器高度z。电算结果如表6-1所示: