C*************************************************************** C TEST.for :- test program #1 for BWR Equation of PROPATH 10 C by S.MOMOKI (May.1997) C*************************************************************** real PR1(5), PR2(5), CP1(10), CP2(10), T, P, Z, ZZZ real ZL, ZV, VL, VV, HL, HV, SL, SV character*40, IDENTM, ID C -------------------------------------------------------------- C C -------------------------------------------------------------- call KPAMES(1,1) call STNKAS(1,1) C -------------------------------------------------------------- C START2() C -------------------------------------------------------------- C PR(Molecular weight and Critical Constants) for R22 C JAR, Thermophysical Properties of Refrigerants(R22) PR1(1) = 86.469 PR1(2) = 369.3 PR1(3) = 0.0 PR1(4) = 1/(513.0/PR1(1)) C PR(5) : acentric factor PR1(5) = 0.205 C CP(cp at ideal gas condition) for R22 C Tokita,K., Refrigeration No.52 Vol.600 (in Japanese) CP1(1) = 3 CP1(2) = 250 CP1(3) = 350 ZZZ = PR1(1) * 1000 CP1(4) = -917.3 * ZZZ CP1(5) = 28.755 * ZZZ CP1(6) = -2.030E-02 * ZZZ CP1(7) = 2.5351E-03 * ZZZ CP1(8) = -2.185E-06 * ZZZ CP1(9) = 7.1E-10 * ZZZ C PR(Molecular weight and Critical Constants) for R114 C JAR, Thermophysical Properties of Refrigerants(R114) PR2(1) = 170.922 PR2(2) = 418.78 PR2(3) = 0.0 PR2(4) = 1.0/(576.0/PR2(1)) C acentric factor PR2(5) = 0.234 C CP(cp at ideal gas condition) for R114 C Tokita,K., Refrigeration No.52 Vol.600(in Japanese) CP2(1) = 2 CP2(2) = 250 CP2(3) = 400 ZZZ = PR2(1) * 1000 CP2(4) = -0.12916 * ZZZ CP2(5) = 4.7009E-03 * ZZZ CP2(6) = -8.612E-06 * ZZZ CP2(7) = 6.05E-09 * ZZZ CP2(8) = 0.0 * ZZZ CP2(9) = 0.0 * ZZZ C call START2(ECODE, PR1, PR2, CP1, CP2, 0.976) C C -------------------------------------------------------------- C IDENTM C -------------------------------------------------------------- ID = 'aaaaaaaaa' ID = IDENTM(1,'V') write(*,*) 'Version ', ID C -------------------------------------------------------------- C CRPM, FCM C -------------------------------------------------------------- TC1 = CRPM(1, 'T') PC1 = CRPM(1, 'P') VC1 = CRPM(1, 'V') TC2 = CRPM(2, 'T') PC2 = FCM(2, 'P') VC2 = CRPM(2, 'V') write(*,900) 900 format (' Critical constants of each component') write(*,901) 1, TC1, PC1, VC1 write(*,901) 2, TC2, PC2, VC2 901 format(' ', I2, ':', f7.2, '[C], ', f8.3, '[bar]', + e10.3, '[m3]') C ************************************************************ C C MAIN LOOP C C ************************************************************ do 10 IZZZ = -500, 10500, 500 ZZZ = IZZZ Z0 = ZZZ / 10000.0 do 11 P0 = 1, 16, 3 do 12 T0 = -20, 70, 20 V = 1 VV = 1 VL = 1 ZL = 0 ZV = 0 write(*,3000) 3000 format(' T[C] P[bar] PHASE Qlty', + ' Z[kg/kg] X Y RHOb,l,v[kg/m3]', + ' Hb,l,v[kJ/kg] Sl,v[kJ/kgK] ') 3001 format(a12, ' ', f6.1, f6.2, i5, f8.3, + f6.3, f6.3, f6.3, f8.2, f6.0, f6.2, + f7.1, f7.1, f7.1, f7.3, f7.3, f7.3) C -------------------------------------------------------------- C IPHASE, QMIX C -------------------------------------------------------------- P = P0 T = T0 Z = Z0 IP = IPHASE(T, P, Z) XX = QMIX(T,P,Z) call SUBXY(ECODE, T, P, ZL, ZV, VL, VV, HL, HV, SL, SV) write(*,3001) 'SUBXY', + T, P, IP, XX, Z, ZL, ZV, -1.0, 1/VL, 1/VV, + -1.0, HL/1e3, HV/1e3, -1.0, SL/1e3, SV/1e3 call SUBMIX(1,ECODE, T, P, Z, V, H, S) write(*,3001) 'SUBMIX_1', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 T = 0 call SUBMIX(2,ECODE, T, P, Z, V, H, S) write(*,3001) 'SUBMIX_2', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 C -------------------------------------------------------------- C TBP, TDP, SUBTD, SUBTB C -------------------------------------------------------------- P = P0 T = 0 Z = Z0 call SUBTB(ECODE, T, P, Z, V, H, S) T2 = TBP(P, Z) write(*,3001) 'SUBTB', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 T = 0 T2 = TDP(P, Z) call SUBTD(ECODE, T, P, Z, V, H, S) write(*,3001) 'SUBTD', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 C -------------------------------------------------------------- C PBT, PDT, SUBPD, SUBPB C -------------------------------------------------------------- P = 0 T = T0 Z = Z0 P2 = PBT(T, Z) call SUBPB(ECODE, T, P, Z, V, H, S) write(*,3001) 'SUBPB', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 P = 0 P2 = PDT(T, Z) call SUBPD(ECODE, T, P, Z, V, H, S) write(*,3001) 'SUBPD', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 C -------------------------------------------------------------- C SUBPUR C -------------------------------------------------------------- P = P0 T = T0 call SUBPUR(1,ECODE, T, P, V, H, S) write(*,3001) 'SUBPUR1', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 call SUBPUR(2,ECODE, T, P, V, H, S) write(*,3001) 'SUBPUR2', + T, P, IP, XX, Z, -1.0, -1.0, 1/V, -1.0, -1.0, + H/1e3, -1.0, -1.0, S/1e3, -1.0, -1.0 C -------------------------------------------------------------- C SUBTSP, TSPM C -------------------------------------------------------------- P = P0 I = 1 T = 0 T2 = TSPM(I,P) call SUBTSP(I, ECODE, T, P, VL, VV, HL, HV, SL, SV) write(*,3001) 'SUBTSP1', + T, P, IP, XX, Z, ZL, ZV, -1.0, 1/VL, 1/VV, + -1.0, HL/1e3, HV/1e3, -1.0, SL/1e3, SV/1e3 I = 2 T = 0 T2 = TSPM(I,P) call SUBTSP(I, ECODE, T, P, VL, VV, HL, HV, SL, SV) write(*,3001) 'SUBTSP2', + T, P, IP, XX, Z, ZL, ZV, -1.0, 1/VL, 1/VV, + -1.0, HL/1e3, HV/1e3, -1.0, SL/1e3, SV/1e3 C C -------------------------------------------------------------- C SUBPST, PSTM C -------------------------------------------------------------- T = T0 I = 1 P = 0 P2 = PSTM(I,T) call SUBPST(I, ECODE, T, P, VL, VV, HL, HV, SL, SV) write(*,3001) 'SUBPST1', + T, P, IP, XX, Z, ZL, ZV, -1.0, 1/VL, 1/VV, + -1.0, HL/1e3, HV/1e3, -1.0, SL/1e3, SV/1e3 I = 2 P = 0 P2 = PSTM(I,T) call SUBPST(I, ECODE, T, P, VL, VV, HL, HV, SL, SV) write(*,3001) 'SUBPST2', + T, P, IP, XX, Z, ZL, ZV, -1.0, 1/VL, 1/VV, + -1.0, HL/1e3, HV/1e3, -1.0, SL/1e3, SV/1e3 write(*,*) write(*,*) 12 continue 11 continue 10 continue END