| 查看: 1047 | 回复: 0 | ||||
| 当前主题已经存档。 | ||||
[交流]
aspen plus反应精馏模拟子程序编译
|
||||
|
用aspen-plus radfrac模块做反应精馏模拟其动力学子程序的编译(我用fortran编译了一些哪位高手帮我完成一下) 反应式:乙酸+乙醇——水+乙酸乙酯 反应速率:R=1.38kx12 x1醋酸摩尔分率 k=3.7-2710/T 汽液平衡常数:醋酸:lgK1=0.0225T-1.666 乙醇:lgK2=6.5-2300/T 水:lgK3=6.4-2300/T 乙酸乙酯:lgK4=6.8-2300/T 如下: SUBROUTINE USRKNT (N, NCOMP, NR, NRL, NRV, 2 T, TLIQ, TVAP, P, PHFRAC, 3 F, X, Y, IDX, NBOPST, 4 KDIAG, STOIC, IHLBAS, HLDLIQ, TIMLIQ, 5 IHVBAS, HLDVAP, TIMVAP, NINT, INT, 6 NREAL, REAL, RATES, RATEL, RATEV, 7 NINTB, INTB, NREALB, REALB, NIWORK, 8 IWORK, NWORK, WORK) C IMPLICIT NONE C C DECLARE VARIABLES USED IN DIMENSIONING C INTEGER NCOMP, NR, NRV, NINT, NINTB, + NREALB,NIWORK,NWORK C #include "ppexec_user.cmn" EQUIVALENCE (RMISS, USER_RUMISS) EQUIVALENCE (IMISS, USER_IUMISS) C C C*********************************************************************** C C C C DECLARE ARGUMENTS C INTEGER NRL(3),IDX(NCOMP),NBOPST(6), + INT(NINT),INTB(NINTB), + IWORK(NIWORK),N,KDIAG, IHLBAS, + IHVBAS,NREAL REAL*8 PHFRAC(3),X(NCOMP,3),Y(NCOMP), + STOIC(NCOMP,NR),RATES(NCOMP), + RATEL(1),RATEV(NRV), + REALB(NREALB),WORK(NWORK), T,TLIQ, + TVAP,P,F,HLDLIQ,TIMLIQ REAL*8 HLDVAP,TIMVAP C C DECLARE LOCAL VARIABLES C INTEGER IMISS REAL*8 REAL(NREAL), RMISS C C BEGIN EXECUTABLE CODE C INTEGER I REAL*8 K,RATE K=3.7-2710/T RATES=K*X(1,1)*X(1,1) RETURN END |
» 收录本帖的淘帖专辑推荐
Aspen Plus 模拟化学平衡和吸收反应 |
» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有5人回复
看来今天不会放榜了?
已经有10人回复
时间戳又变了
已经有19人回复
filecode,4个jtjc了
已经有14人回复
让我中一个面上吧!
已经有13人回复
人气不行了
已经有7人回复
今日不放榜?网传国自然预计 8 月 27 日可查结果
已经有18人回复
感觉是下周放榜了
已经有12人回复
我面上完蛋了
已经有8人回复
只有每年这种时候来逛逛小木虫
已经有25人回复









回复此楼