清华大学《化工系统工程基础》第二题-孙悦.docxVIP

模拟上机作业二：甲苯加氢制苯工艺流程模拟建立稳态模型流程图仿照题中所给工艺流程图，在Aspen Plus中绘制流程图如上所示。模拟参数设定根据题目中的条件，以及作业excel中的数据，设定模拟参数。(1) 化合物定义根据题目设置如下：(2) 物性方法选择根据化工热力学所学知识，对于氢气、甲烷、苯-甲苯这种非极性体系，可以采用PENG-ROB或者RK-SOVAE等状态方程模型。其中计算一些相对较为简单地物质更为适合，这里便使用RK-SOVAE状态方程模型。(3) 进料物流设置液体甲苯（132.7kmo/h）、新鲜氢气（222.4kmol/h，氢气摩尔分数为0.97，甲烷摩尔分数为0.03）FFH2：FFT01：(4) 设备参数设定设备参数设定过程中较为复杂。由于设备较多，单位错误等问题，极易出现。因而要格外小心。精馏塔C1塔有6块塔板，第1块为进料板，塔顶不设冷凝器，C1塔塔釜采出作为苯塔C2的进料。C2塔有27块塔板，第15块为进料板，C2塔塔釜采出作为回收塔C3的进料。C3塔有7块塔板，第5块为进料板， C3塔塔釜采出为副产品联苯。C1塔：这里蒸出物流与进料物流的比，是C1的出料D1与进料V5OUT的比，由excel中数据计算得到。C2、C3同理。这里压力设置，参考excel中D1数据的压力项为1135kPa。C2、C3同理。C2塔：C3塔：压缩机、泵和阀门由题目所给excel的数据中进出阀门、泵或压缩机的压力及压力差来进行设定。COMP、P1、P2、P3：COMP：类似于COMP，P1、P2、P3的压力变化值均是由excel数据进出泵物流的压差来确定。表1 泵的Pressure increase ( kPa )P1P2P3523344.84101.9V3~V10、VLV-02、VLV-100、VLV-102：阀门设置过程中，出现了一个问题，就是将压力设置选择为Pressure drop时，进行后期计算的时候经常会现warnings，然而若将其改为Outlet pressure，这个问题就就解决了。分析其原因，可能是由于在后来的收敛过程中由于回流比蒸出比等不同，物料在经过回路一遍一遍的计算过程中，可能会与题目中给出的excel中数据偏差越来越大，导致后来压力过大或者过小，进程无法进行，出现错误。而改为Outlet pressure之后，压力变化后相会对稳定一些。V3：同理，类似于V3，V4~V10、VLV-02、VLV-100、VLV-102的压力设定均是由excel中阀门的出口物流的压力设定的。如下表所示：表2 阀门的Outlet pressure ( kPa )V4V5V6V7V833763431790.8335.7239.2V9V10VLV-02VLV-100VLV-102315253345240654410换热器换热器参数设置由题目中条件得出：经换热器FEHE与反应器出口物流换热升温，再进入加热炉Fur加热到621.1℃，进入一绝热管式反应器PFR。反应器出口温度667.7℃，与70.76kmol/h冷流体混合经急冷（以防结焦）到606.2℃进入换热器FEHE冷却到125.5℃。反应器出口流经水冷器进一步冷却到45℃，经闪蒸罐Sep分离。压力参数，则由出口物流的压力设定。COND：E-100：FEHE：混合器由excel数据可知道混合器压降为0。反应器PFR参数设定：题目条件为绝热管式反应器PFR，反应器内无催化剂，管径2905mm，管长17374mm。化学反应设定：在反应器中发生如下反应：主反应 C7H8（甲苯）＋H2 —→ C6H6（苯）＋CH4（甲烷） -41840 kJ/kmol副反应2 C6H6（苯）—→ C12H10（联苯）＋H2＋8136 kJ/kmol两反应速率方程如下：式中，反应速率单位为kmol/(s·m3)；，，，分别为甲苯、氢气、苯、联苯的压力，单位为kpa；活化能单位为kJ/kmol；温度单位为K；R=8.314kj/(kmol·K)。反应动力学设定主反应为：Driving Force：副反应为：可逆反应，r=[ Kinetic factor][Driving force]。r2可化为故设置如下：Driving Force：由下式推出由coefficients for driving force constant (A、B、C、D)定义可推出A、B、C、D的值。闪蒸罐SEP由题目：反应器出口流经水冷器进一步冷却到45℃，经闪蒸罐Sep分离。及excel中数据设定如下：分流器利用excel中给出的出口物流流量比可得出分流率。TEE-100：TEE-101：至此，稳态模型建立完毕。物流切割及流程收敛切割物流(1) 信息流图省略压缩机交换器等。可画出信息流图如下：节点分别表示：1