24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 yqyq0000 的 9 个金币

yqyq0000

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 231.8
帖子: 16
在线: 48小时
虫号: 736275
注册: 2009-03-31

[交流] 【求助】有关ASPEN的FORTRAN问题

请哪位高手帮忙看看
我把方程编辑好了通过编译  但是在ASPEN中运行出现问题
很大程度上我认为是单位不一致造成的  在ASPEN中反应速率单位默认为是Kg-mole/sec
但是我们通常的动力学方程单位不是这个  因此会造成一些错误
附件上是我的方程还有动力学方程及相关的数据，希望同行高手们帮忙编辑运行试试，能给我解决下问题。
请哪位高手帮忙看看
C User Kinetics Subroutine for gasification
C nov  10 2009
C
   SUBROUTINE GASFC (SOUT, NSUBS,  IDXSUB, ITYPE,  NINT,
   2                INT, NREAL,  REAL,    IDS, NPO,
   3                NBOPST, NIWORK, IWORK, NWORK,  WORK,
   4                NC,    NR,    STOIC, RATES,  FLUXM,
   5                FLUXS,  XCURR,  NTCAT, RATCAT, NTSSAT,
   6                RATSSA, KCALL,  KFAIL, KFLASH, NCOMP,
   7                IDX, Y,    X,       X1,    X2,
   8                NRALL,  RATALL, NUSERV, USERV,  NINTR,
   9                INTR, NREALR, REALR, NIWR, IWR,
   1                NWR, WR)
C
   IMPLICIT NONE
C
C    DECLARE VARIABLES USED IN DIMENSIONING
C
   INTEGER NSUBS, NINT,  NPO, NIWORK,NWORK,
   +       NC, NR, NTCAT, NTSSAT,NCOMP,
   +       NRALL, NUSERV,NINTR, NREALR,NIWR,
   +       NWR
C
#include "ppexec_user.cmn"
   EQUIVALENCE (RMISS, USER_RUMISS)
   EQUIVALENCE (IMISS, USER_IUMISS)
C
C
C
C
C.....RCSTR...
#include "rcst_rcstri.cmn"
#include "rxn_rcstrr.cmn"
C
C.....RPLUG...
#include "rplg_rplugi.cmn"
#include "rplg_rplugr.cmn"
   EQUIVALENCE (XLEN, RPLUGR_UXLONG)
   EQUIVALENCE (DIAM, RPLUGR_UDIAM)
C
C.....RBATCH...
#include "rbtc_rbati.cmn"
#include "rbtc_rbatr.cmn"
C
C.....PRES-RELIEF...
#include "rbtc_presrr.cmn"
C
C.....REACTOR (OR PRES-RELIEF VESSEL) PROPERTIES...
#include "rxn_rprops.cmn"
   EQUIVALENCE (TEMP, RPROPS_UTEMP)
   EQUIVALENCE (PRES, RPROPS_UPRES)
   EQUIVALENCE (VFRAC, RPROPS_UVFRAC)
   EQUIVALENCE (BETA, RPROPS_UBETA)
   EQUIVALENCE (VVAP, RPROPS_UVVAP)
   EQUIVALENCE (VLIQ, RPROPS_UVLIQ)
   EQUIVALENCE (VLIQS, RPROPS_UVLIQS)
C
#include "shs_stwork.cmn"
   EQUIVALENCE (MKBAS, STWORK_NDUM)
   EQUIVALENCE (MKPHAS, STWORK_NBLM)
   EQUIVALENCE (MTAPP, STWORK_NCOVAR)
   EQUIVALENCE (MKBASS, STWORK_NWR)
   EQUIVALENCE (MTAPPS, STWORK_NIWR)
   EQUIVALENCE (SSALT, STWORK_RDUM1)
   EQUIVALENCE (VSALT, STWORK_RDUM2)
   EQUIVALENCE (FSALT, STWORK_FFSALT)
#include "pputl_ppglob.cmn"
#include "dms_ncomp.cmn"
#include "dms_plex.cmn"
   EQUIVALENCE (IB(1), B(1))
C
C    DECLARE ARGUMENTS
C
   INTEGER IDXSUB(NSUBS),ITYPE(NSUBS), INT(NINT),
   +       IDS(2),NBOPST(6,NPO),IWORK(NIWORK),
   +       IDX(NCOMP), INTR(NINTR),  IWR(NIWR),
   +       NREAL, KCALL, KFAIL, KFLASH,I
   REAL*8 SOUT(1),    WORK(NWORK),
   +    STOIC(NC,NSUBS,NR),  RATES(1),
   +    FLUXM(1),    FLUXS(1),    RATCAT(NTCAT),
   +    RATSSA(NTSSAT),    Y(NCOMP),
   +    X(NCOMP),    X1(NCOMP), X2(NCOMP)
   REAL*8 RATALL(NRALL), USERV(NUSERV),
   +    REALR(NREALR), WR(NWR),    XCURR
C
C
C    DECLARE LOCAL VARIABLES
C
   INTEGER IMISS, MKBAS, MKPHAS, MTAPP, MKBASS,
   +       MTAPPS, LMW, KDIAG, KPHI, KER
   REAL*8 REAL(NREAL), B(1),  RMISS, XLEN,
   +    DIAM,  TEMP,  PRES,  VFRAC, BETA,
   +    VVAP,  VLIQ,  VLIQS, SSALT, VSALT,
   +    FSALT,  RTEMP,  RPRES, FUGA(NC), PHI(NCOMP),
   +    DPHI(NCOMP), RATE, K1, K2, K3, K4, K5,WC,
   +    RPAR, YL(NC), PHIL(NC), XC1WT, XC2WT,
   +    RATE1, RATE2, RATE3
#include "dms_ipoff1.cmn"
C
C    INITIALIZE RATES
C
C
C  STATEMENT FUNCTIONS FOLLOW
C
C
   DO 100 I = 1, NC
      RATES(I) = 0D0
YL(I) = 0D0
PHIL(I) = 0D0
100  CONTINUE
C
C
C    The structure in the array SOUT is as follows:
C
C    SOUT(1) - SOUT(NCC) : Component flowrates(kg-moles/sec)
C    SOUT(NCC+1)       : Total flowrates(kg-moles/sec)
C    SOUT(NCC+2)       : Temperature(K)
C    SOUT(NCC+3)       : Pressure(N/SQM)
C    SOUT(NCC+4)       : Mass enthalpy(J/KG)
C    SOUT(NCC+5)       : Molar vapor fraction
C    SOUT(NCC+6)       : Molar liquid fraction
C    SOUT(NCC+7)       : Mass entropy(J/KG-K)
C    SOUT(NCC+8)       : Mass density(KG/CUM)
C    SOUT(NCC+9)       : Molecular Weight
C
C
C    Set Reactor Temperature & Pressure
C
   RTEMP = SOUT(NCOMP_NCC+2)
RPRES = SOUT(NCOMP_NCC+3)

C
C    Compute Vapor fugacity
C
   KDIAG=4
KPHI=1
CALL PPMON_FUGV(RTEMP, RPRES, Y, NCOMP, IDX,
   +             NBOPST, KDIAG, KPHI, PHI, DPHI,
   +             KER)
   DO 101 I = 1, NCOMP
      YL(IDX(I)) = Y(I)
      PHIL(IDX(I)) = PHI(I)
101  CONTINUE

DO 102 I=1, NC
FUGA(I)=RPRES * YL(I) * PHIL(I)
102 CONTINUE

C    Compute RATES OF REACTIONS:
C
C
   realr1=8710
realr2=14938.3567
realr3=247
realr4=17509.9789

K1 = REALR1 * EXP( REALR2*10000 /PPGLOB_RGAS / RTEMP )
K2 = REALR3 * EXP( REALR4*10000 / PPGLOB_RGAS / RTEMP )
XC1WT = 1-(1-(1/3)*2.5*6000000*K1*FUGA(1)/RPRES)**3
XC2WT = 1-(1-(1/3)*2.5*6000000*K2*FUGA(2)/RPRES)**3
C
C
C    CALCULATE RATES OF ALL COMPONENTS
C
C    RATES - The rates for all components
C
C
C
RATES(1) =K1 * FUGA(1) *(1-XC1WT)**(2/3)
   RATES(2) =K2 *FUGA(2) * (1-XC2WT)**(2/3)

   END

[ Last edited by yqyq0000 on 2009-11-19 at 13:10 ]

回复此楼

» 猜你喜欢

AI论文写作工具：是科研加速器还是学术作弊器？已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2009-11-10 20:19:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

04nylxb

木虫 (正式写手)

应助: 33 (小学生)
金币: 2321.9
散金: 46
红花: 4
帖子: 824
在线: 262.6小时
虫号: 817223
注册: 2009-07-28
性别: GG
专业: 工程热物理相关交叉领域

★ ★
nono2009(金币+1,VIP+0):谢谢应助！ 11-18 14:04
nono2009(金币+0,VIP+0):欢迎去投票，寻找更多同行，同时领取金币。http://emuch.net/bbs/viewthread.php?tid=1649328 11-18 14:06
yqyq0000(金币+1,VIP+0): 11-19 13:07

用固体solid模型试试看，自定义煤的CHONS

赞一下(22人)

回复此楼

集中精力发文章

2楼2009-11-18 09:15:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 yqyq0000 的主题更新

返回列表