aspen二元体系共沸组分_aspen模拟时有共沸物。

换了一个思路,请lz,专家们看一下是否可行.方法为 (1)用aspen检查一下系统的vlle状态 (2)利用系统特点,先进行ll分离去水相以求节能和简化下一步精馏处理 (3)精馏变为常规stripping塔 塔进料,塔板数塔压可以进一步优化,方法为常规方法.在aspen相图生成时用了dcpd-h2o近似参数,由溶解度导出.下面为proii输入文件.供参考 $ generated by pro/ii keyword generation system $ generated on: fri oct 31 21:36:16 2008 title dimension metric, pres=kpag, stdtemp=0, stdpres=0 sequence simsci calculation rvpbasis=apin, tvp=37.778 component data libid 1,toluene/2,h2o/3,dcpd, bank=simsci,process assay conversion=api94, curvefit=improved, kvreconcile=tails thermodynamic data method system(vlle)=nrtl, set=nrtl01, default method system(vlle)=uniq, set=uniq01 method system(vlle)=unif, set=unif01 method system=nrtl, set=nrtl02 stream data property stream=feed, temperature=26.667, pressure=0, phase=m, & rate(wt)=1000, composition(wt)=1,50/2,38/3,12, normalize $ original input data for stream reflux $ property stream=reflux, temperature=26.667, pressure=-55.158, & $ phase=m, composition(m,kgm/h)=1,2.70979/2,0.0072102/ & $ 3,0.0646412 prop strm=reflux ,temp= 26.6666870 , pres=-55.1580811 ,& comp= & 1, 10.2083763 / 2,2.635874175e-02 / 3,5.953210990e-03 unit operations flash uid=f1 feed feed product v=s2, l=hc_phase, w=s9 adiabatic method set=unif01 flash uid=f2 feed hc_phase product l=twr_feed bubble column uid=t1 parameter tray=10,io=50 feed twr_feed,3/reflux,1 product btms(m)=dicp_prod,0.905126, ovhd(m)=s5,15.7149, & supersede=on duty 1,10,0.151328,reboiler $pspec ptop=-55.158, dpcolumn=13.79 print proptable=part estimate model=conventional, rratio=3 $estimate model=conventional, rratio=3 $temperature 1,26.667 pressure 1,-55.1581/2,-53.6259/3,-52.0937/4,-50.5616/ & 5,-49.0294/6,-47.4972/7,-45.9651/8,-44.4329/ & 9,-42.9007/10,-41.3686 spec id=col1spec2, stream=dicp_prod, rate(wt,kg/h), comp=3,wet, & divide, reffeed, rate(wt,kg/h), comp=3,wet, & value=0.99 vary dname=reboiler reboiler type=kettle method set=nrtl02 temp & 1, 85.153969 / 2, 87.026347 / 3, 92.225530 / & 4, 105.354766 / 5, 125.135469 / 6, 138.452611 / & 7, 144.381337 / 8, 146.993789 / 9, 148.465600 / & 10, 149.566463 vapo & 1, 15.7149517 / 2, 18.4262159 / 3, 18.1425448 / & 4, 16.6170905 / 5, 15.4714604 / 6, 15.3177297 / & 7, 15.7366616 / 8, 16.0143252 / 9, 16.1395949 / & 10, 16.2024316 liqu & 1, 12.9853293 / 2, 12.7016582 / 3, 17.5222160 / & 4, 16.3765860 / 5, 16.2228552 / 6, 16.6417871 / & 7, 16.9194507 / 8, 17.0447205 / 9, 17.1075572 / & 10, 0.9051256 flash uid=f3 feed s5 product w=s7, l=s8 iso temperature=26.667 splitter uid=sp1 feed s8 product m=tln_prod, m=reflux operation option=fill spec stream=reflux, rate(lbm/h),total,wet, divide, reffeed, & rate(lbm/h),wet, value=0.5 controller uid=rflx_cntl spec stream=reflux, rate(wt,kg/h),total,wet, divide, & stream=tln_prod, rate(wt,kg/h),total,wet, & value=1.87 vary splitter=sp1, spec(1) cparameter iprint, nostop mixer uid=m1 feed s9,s7 product m=waste_h2o calculator uid=ca1 define p(1) as stream=dicp_prod, rate(wt,kg/h), comp=3,wet define p(2) as stream=feed, rate(wt,kg/h), comp=3,wet procedure r(1)=p(1)/p(2)*100 return end [ ]