| 查看: 326 | 回复: 1 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 yqyq0000 的 9 个金币 | |||
[交流]
【求助】有关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 ] |
回复此楼» 猜你喜欢
278求调剂
已经有3人回复
-
化工调剂303分，过四级
已经有27人回复
-
265求调剂
已经有20人回复
-
348分环境工程·调剂
已经有6人回复
-
322求调剂
已经有3人回复
-
335求调剂
已经有7人回复
-
311（085601）求调剂
已经有15人回复
-
081700化学工程与技术 一志愿中海洋 323 求调剂学校
已经有6人回复
-
327求调剂
已经有8人回复
-
292求调剂
已经有21人回复
04nylxb
木虫 (正式写手)
- 应助: 33 (小学生)
- 金币: 2321.9
- 散金: 46
- 红花: 4
- 帖子: 824
- 在线: 262.6小时
- 虫号: 817223
- 注册: 2009-07-28
- 性别: GG
- 专业: 工程热物理相关交叉领域
2楼2009-11-18 09:15:40
