C*************************************************************** C TEST2.for :- test program #2 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 real ZL, ZV, VL, VV, HL, HV, SL, SV C -------------------------------------------------------------- C C -------------------------------------------------------------- call KPAMES(3,1) 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 -------------------------------------------------------------- C CRPM C -------------------------------------------------------------- TC1 = CRPM(1, 'T') PC1 = CRPM(1, 'P') VC1 = CRPM(1, 'V') TC2 = CRPM(2, 'T') PC2 = CRPM(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], ', f7.3, '[bar]', + e10.3, '[m3]') C ************************************************************ C C MAIN LOOP C C ************************************************************ ZZZMAX = 10000.0 do 10 ZZZ = 0.0, ZZZMAX, 500.0 Z0 = ZZZ / ZZZMAX TB = 0.0 do 11 P0 = 1, 16, 3 IP = 1 do 12 T0 = -20, 70, 3 V = 1 VV = 1 VL = 1 ZL = 0 ZV = 0 C -------------------------------------------------------------- C IPHASE, QMIX C -------------------------------------------------------------- P = P0 T = T0 Z = Z0 if (Z . le. 0.0) goto 12 if (IP .eq. 3 ) goto 12 C already vapor if (T .lt. TB ) goto 12 IP = IPHASE(T, P, Z) C write(*,*) Z, P, ' T=', T, IP, ', TB=', TB if (IP .eq. 3) then TB = T - 3.0 goto 12 endif if (IP .ne. 2) goto 12 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) 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 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 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 write(*,*) write(*,*) 12 continue 11 continue 10 continue END