10 SCREEN 0,0,0:KEY OFF
20 CLS
30 '
40 DIM A(8),B(8),C(8),D(8),E(8),F(8),G(8),HW(8,28),VW(8,28)
50 FOR I=1 TO 8:READ A(I),B(I),C(I),D(I),E(I),F(I),G(I):NEXT
60 DATA -.17724,3.83836,11.48345,31.1311,8.762969e-5,-2.78794e-8,86.594,-5.280126e-7,2.99461e-5,1.521874e-4,6.62512e-5,8.408856e-10,1.86401e-14,.01596,-.48799,304.717614,9.8299035
70 DATA -16.455274,9.474745e-4,-1.363366e-6,19.53953,2.662e-3,457.5802,-.176959,.826862,-4.601876e-7,6.3181e-11,-2.3928,-.15115567,3.671404,11.622558,30.832667,8.74117e-5
80 DATA -2.62306e-8,54.55,-.14129,2.258225,3.4014802,14.438078,4.222624e-5,-1.569916e-8,1100.5,-1.67772e-4,4.272688e-3,.01048048,.05801509,9.101291e-8,-2.7592e-11,.11801
90 DATA -1.476933e-4,1.2617946e-3,3.44201e-3,-.08494128,6.89138e-8,-2.4941e-11,1.97364
100 LOCATE 2,15:PRINT "Turbine Steam Rates & Properties Calculation Program"
110 LOCATE 4,10:INPUT "Inlet Pressure (psia)                           :";P1
120 LOCATE 5,10:INPUT "Temp(øF) (If saturated steam, input 0 for temp) :";T1
130 LOCATE 6,10:INPUT "Exit pressure (psia)                            :";P2
140 LOCATE 7,10:INPUT "Expansion efficiency %                          :";EF
150 IF T1=0 THEN GOTO 180
160 P=P1:T=T1:GOSUB 520
170 H1=Z:S1=SV:GOTO 200
180 P=P1:GOSUB 630
190 H1=HV:S1=SV:T1=TSAT
200 P=P2:GOSUB 630
210 IF S1>SV THEN 260
220 X=(S1-SL)/(SV-SL):H2S=HV*X+(1-X)*HLIQ
230 H2=H1-.01*EF*(H1-H2S)
240 IF H2>HV THEN 260
250 XF=(H2-HLIQ)/(HV-HLIQ):DELH=H1-H2:T2=TSAT:GOTO 360
260 T2=TSAT:RA=1000-TSAT
270 P=P2:T=T2:GOSUB 520
280 S2C=SV:IF ABS((S1-S2C)/S1)<.002 THEN 300
290 RA=.5*RA:T2=T2+SGN(S1-S2C)*RA:GOTO 270
300 XF=1:T=T2:GOSUB 520:H2S=Z
310 H2=H1-.01*EF*(H1-H2S)
320 T2=TSAT:RA=1000-TSAT
330 P=P2:T=T2:GOSUB 520
340 H2C=Z:IF ABS((H2-H2C)/H2)<.002 THEN 360
350 RA=.5*RA:T2=T2+SGN(H2-H2C)*RA:GOTO 330
360 ASR=3413/(H1-H2):TSR=3413/(H1-H2S)
370 CLS
380 LOCATE 2,15:PRINT "Turbine Steam Rates & Properties Calculation Program"
390 LOCATE 4,10:PRINT "Inlet press (psia)          :";P1
400 LOCATE 5,10:PRINT "Inlet Temp (F)              :";T1
410 LOCATE 6,10:PRINT "Inlet Enthalpy              :";H1
420 LOCATE 7,10:PRINT "Inlet Entropy               :";S1
430 LOCATE 8,10:PRINT "Exit pressure               :";P2
440 LOCATE 9,10:PRINT "Exit Temp                   :";T2
450 LOCATE 10,10:PRINT "Exit Enthalpy               :";H2
460 LOCATE 11,10:PRINT "Exit Saturated Temp         :";TSAT
470 LOCATE 12,10:PRINT "Actual steam req. (ASR)     :";ASR
480 LOCATE 13,10:PRINT "Theor. steam req. (TSR)     :";TSR
490 LOCATE 14,10:PRINT "Expansion efficiency        :";EF
500 LOCATE 15,10:PRINT "Exit quality                :";XF
510 END
520 T=273.1+(T-32)/1.8:P=P/14.696
530 K4=80870!/T/T:L1=10^K4*(-2641.62/T):M1=1.89+L1:N1=M1*P*P/T/T
540 O1=2+(372420!/T/T):Q1=O1*L1:R1=1.89+Q1:U1=(.21828*T-126970!/T):V4=2*U1*R1-(M1/T)*126970!
550 W1=82.54-162460!/T:Y1=2*W1*R1-(M1/T)*162460!
560 F1=(((U1*M1*N1+W1)*N1/P+1)*M1+4.5504*T/P)*.0160185
570 Z=775.6+.63296*T+1.62467E-04*T*T+20.5697*LOG(T)
580 Z=Z+.043557*(R1*P+.5*N1*(Y1+M1*(W1+V4*N1)))
590 SV=(((U1*M1-2*V4)*.5*M1*N1-Y1)*.5*N1+(M1-R1)*P)/T*(-.0241983)-.355579-11.4276/T
600 SV=SV+1.8052E-04*T-.11022*LOG(P)+.35164*LOG(T)
610 P=14.696*P:T=(T-273.15)*1.8+32
620 RETURN
630 FOR I=1 TO 8:AA(I)=A(I)*P+B(I)/P+C(I)*P^.5+D(I)*LOG(P)+E(I)*P*P+F(I)*P^3+G(I):NEXT I
640 TSAT=AA(1):HLIQ=AA(5):HV=AA(6):SL=AA(7):SV=AA(8)
650 RETURN