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 (Kg/cm2 a)                       :";INPR:P1=14.7*INPR
120 LOCATE 5,10:INPUT "Temp(øC) (If saturated steam, input 0 for temp) :";TEM:T1=(TEM*1.8)+32
130 LOCATE 6,10:INPUT "Exit pressure (Kg/cm2 a)                        :";OUTPR:P2=14.7*OUTPR
140 LOCATE 7,10:INPUT "Expansion efficiency %                          :";EF
150 LOCATE 8,10:INPUT "Enter the steam flow rate for power output(Kg/hr):";STEAM
160 IF TEM=0 THEN GOTO 190
170 P=P1:T=T1:GOSUB 540
180 H1=Z:S1=SV:GOTO 210
190 P=P1:GOSUB 650
200 H1=HV:S1=SV:T1=TSAT
210 P=P2:GOSUB 650
220 IF S1>SV THEN 270
230 X=(S1-SL)/(SV-SL):H2S=HV*X+(1-X)*HLIQ
240 H2=H1-.01*EF*(H1-H2S)
250 IF H2>HV THEN 270
260 XF=(H2-HLIQ)/(HV-HLIQ):DELH=H1-H2:T2=TSAT:GOTO 370
270 T2=TSAT:RA=1000-TSAT
280 P=P2:T=T2:GOSUB 540
290 S2C=SV:IF ABS((S1-S2C)/S1)<.002 THEN 310
300 RA=.5*RA:T2=T2+SGN(S1-S2C)*RA:GOTO 280
310 XF=1:T=T2:GOSUB 540:H2S=Z
320 H2=H1-.01*EF*(H1-H2S)
330 T2=TSAT:RA=1000-TSAT
340 P=P2:T=T2:GOSUB 540
350 H2C=Z:IF ABS((H2-H2C)/H2)<.002 THEN 370
360 RA=.5*RA:T2=T2+SGN(H2-H2C)*RA:GOTO 340
370 ASR=3413/(H1-H2):TSR=3413/(H1-H2S)
380 CLS
390 LOCATE 2,15:PRINT "Turbine Steam Rates & Properties Calculation Program"
400 LOCATE 4,10:PRINT "Inlet press    (kg/cm2 a)   :";INPR
410 LOCATE 5,10:PRINT "Inlet Temp     (øC)         :";TEM
420 LOCATE 6,10:PRINT "Inlet Enthalpy (Kcal/kg)    :";H1*.252/.454
430 LOCATE 7,10:PRINT "Inlet Entropy  (Kcal/kg/øC) :";S1*.252/.454/1.8
440 LOCATE 8,10:PRINT "Exit pressure  (Kg/cm2 a)   :";OUTPR
450 LOCATE 9,10:PRINT "Exit Temp      (øC)         :";(T2-32)*5/9
460 LOCATE 10,10:PRINT "Exit Enthalpy  (Kcal/kg)    :";H2*.252/.454
470 LOCATE 11,10:PRINT "Exit Saturated Temp (øC)    :";(TSAT-32)*5/9
480 LOCATE 12,10:PRINT "Actual steam req. (ASR, Kg/KWh)  :";ASR*.454
490 LOCATE 13,10:PRINT "Theor. steam req. (TSR, Kg/KWh)  :";TSR*.454
500 LOCATE 14,10:PRINT "Expansion efficiency        :";EF
510 LOCATE 15,10:PRINT "Exit quality                :";XF
520 LOCATE 16,10:PRINT "Power Output  (KW)          :";STEAM/(ASR*.454)
530 END
540 T=273.1+(T-32)/1.8:P=P/14.696
550 K4=80870!/T/T:L1=10^K4*(-2641.62/T):M1=1.89+L1:N1=M1*P*P/T/T
560 O1=2+(372420!/T/T):Q1=O1*L1:R1=1.89+Q1:U1=(.21828*T-126970!/T):V4=2*U1*R1-(M1/T)*126970!
570 W1=82.54-162460!/T:Y1=2*W1*R1-(M1/T)*162460!
580 F1=(((U1*M1*N1+W1)*N1/P+1)*M1+4.5504*T/P)*.0160185
590 Z=775.6+.63296*T+1.62467E-04*T*T+20.5697*LOG(T)
600 Z=Z+.043557*(R1*P+.5*N1*(Y1+M1*(W1+V4*N1)))
610 SV=(((U1*M1-2*V4)*.5*M1*N1-Y1)*.5*N1+(M1-R1)*P)/T*(-.0241983)-.355579-11.4276/T
620 SV=SV+1.8052E-04*T-.11022*LOG(P)+.35164*LOG(T)
630 P=14.696*P:T=(T-273.15)*1.8+32
640 RETURN
650 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
660 TSAT=AA(1):HLIQ=AA(5):HV=AA(6):SL=AA(7):SV=AA(8)
670 RETURN