清华大学《化工系统工程基础》大作业二.docxVIP

模拟上机作业二：甲苯加氢制苯工艺流程模拟稳态模型建立（1）按照题目中的工艺流程ASPEN PLUS中绘制流程图如下：图1 aspen绘制流程图（2）参数设置根据已知条件对aspen中各个设备和物流输入参数如下：图2 输入物流化学成分物性计算方法选择：根据化工热力学所学知识，对于氢气、甲烷、苯-甲苯这种非极性体系，可以采用PENG-ROB或者RK-SOLVE等状态方程模型。在模拟的时候开始采用了PENG-ROB方法，难以收敛，后采用RK-SOLVE模型，可以很好的收敛。故选用RK-SOLVE模型。图3 模型选择输入进料状态：FFH2，FFT01物流：液体甲苯（132.7kmo/h）、新鲜氢气（222.4kmol/h，氢气摩尔分数为0.97，甲烷摩尔分数为0.03）图4 FFH2物流设置图5 FFT01物流设置输入设备参数：C1塔有6块塔板，第1块为进料板，塔顶不设冷凝器，C1塔塔釜采出作为苯塔C2的进料。C2塔有27块塔板，第15块为进料板，C2塔塔釜采出作为回收塔C3的进料。C3塔有7块塔板，第5块为进料板， C3塔塔釜采出为副产品联苯。图6 C1精馏塔设置1图7C1精馏塔设置2图8C1精馏塔设置3注：精馏塔的设置中考虑到aspen不考虑再沸器的一块板，故设置的塔板数要加一，设置时蒸出比参考了附件EXCEL中的数值，计算得到；塔顶的压力采用的是附件EXCEL中塔顶出口物流D1的压力，下同。图9 C2塔设置1图10 C2塔设置2图11 C2塔设置3图12 C3塔设置1图13 C3塔设置2图14 C3塔设置3各换热器参数：（注：温度根据题目中所给条件设定，如COND出口温度应该为45 ℃；E-100出口温度是PFR的入口温度，是加热炉Fur加热到621.1℃，FEHE出口温度为125.5 ℃）图15 COND设置图16 E-100设置图17 FEHE设置分离器：（注：温度为混合物流温度，压力为CONOUT物流的压力，采用EXCEL中值）图18 分离器SEP设置分流器：图19 分流器TEE-100设置图20 分流器TEE-101设置各个泵和阀门压力设置：根据EXCEL中的出入口物流的压力为参考，给出设置的各个泵、压缩机和阀门的压力升和压力降，整理如下：（注：负号为压降）表1 泵、压缩机、阀门压力设置设备COMPP1P2P3V3V4V5V6压力差（kpa）833523344.84101.9-345-689-324-344设备V7V8V9V10VLV100VLV101VLV102压力差（kpa）-827.3-69-365.5-82.7-345-345-303PFR设置：设置为绝热，并输入PFR尺寸：管径2905mm，管长17374mm化学反应设置：首先在Reaction模块中新建反应并输入反应方程式，反应级数选择PWERLAW：图21 化学反应设置然后分别设置各个反应的平衡常数：主反应 C7H8（甲苯）＋H2—→ C6H6（苯）＋CH4（甲烷）－ 41840 kJ/kmol可以得到主反应的kE=254390 kJ/kmol更改单位units为KPa，反应相态为气相Vapor；basis为分压Partial pressure图22 主反应的平衡常数设置图23 反应级数设置副反应：2 C6H6（苯）—→ C12H10（联苯）＋H2＋8136 kJ/kmol副反应为可逆反应，设置比较特殊。根据对话框中的Powerlaw Kinetic expression的定义：r=[Kinetic factor][Drving force]将r2的形式进行转换：所以该反应的Kinetic factor=k=1142200;E=211060 kJ/kmol图24 可逆反应平衡常数设置对于Drving force中的正反应forward的设置，由r2的表达式可以看出：此时的正反应k=1，级数为2，故A和B均为0图25 Driving force 正反应设置而对于逆反应，有k=28.37；E=6150 kJ/kmol;A=lnk=3.345B=-E/R=-6510/8.314=783.017将计算所得的A和B填入，得到逆反应的反应推动力设置如下：图26 Driving force逆反应设置切割物流（1）简化信息流图将上面的流程图转化为简化的信息流图，信息流图如下：（为方便计算，将各股物流编号）图27 简化的信息流图（2）回路切割法进行切割采用回路切割法选择切割物流，对信息流图中的回路整理如下：表2回路矩阵jL123456789101112131411112111131111411111111S1S2S3S4S5S6S7S8S9S10S11S12S13S14其中；；；；；；；