24小时热门版块排行榜    

查看: 930  |  回复: 0
当前主题已经存档。

zoeeden

新虫 (初入文坛)

[交流] 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 模拟化学平衡和吸收反应

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zoeeden 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见