10 SCREEN 0,0,0:KEY OFF
20 CLS
30 DIM A(8),B(8),C(8),D(8),E(8),F(8),G(8),HW(8,28),VW(8,28)
40 FOR I=1 TO 8:READ A(I),B(I),C(I),D(I),E(I),F(I),G(I):NEXT
42 LOCATE 1,20:PRINT "Turbine Efficiency Calculation program"
50 LOCATE 2,10:INPUT "Enter the Inlet Pressure (Kg/cm 2 a) :";PIN:PIN1=(PIN*14.7)-14.7  'psig
60 LOCATE 3,10:INPUT "Enter the Inlet Temperature (øC)     :";TIN
70 LOCATE 4,10:INPUT "Enter the Exhaust Pressure(Kg/cm2 a) :";POUT:POUT1=POUT*14.7 'psia
80 LOCATE 5,20:PRINT "Type of Turbine"
90 LOCATE 7,15:PRINT "1. Single Stage, and Single valve multistage, back pressure"
100 LOCATE 8,15:PRINT "2. Single Stage, and Single valve multistage, condensing"
110 LOCATE 9,15:PRINT "3. Multivalve, Multistage, back pressure"
120 LOCATE 10,15:PRINT "4. Multivalve, Multistage, Condensing"
130 LOCATE 12,10:INPUT "Enter your type of turbine (1-4) :";TURTYPE
140 IF TURTYPE <1 OR TURTYPE >4 THEN 130
150 IF TURTYPE =>3 THEN 170
160 PRRATIO=.92:GOTO 180
170 PRRATIO=.95
180 GOVERNPR=PRRATIO*(PIN1+15) 'psia
190 GOVERNPR1=GOVERNPR/14.7 'Kg/cm2 a
200 IF TURTYPE =3 OR TURTYPE=1 THEN EFFCORR=.88:GOTO 220
210 EFFCORR=.94
220 'h1,h2s,h2sdash
230 P1=PIN1+14.7:P2=POUT1:T1=(TIN*1.8)+32 'øF
240 GOSUB 880:H1=H1:H2S=H2S:ENTHIN=H1:ISENENTHOUT=H2S  'Btu/lb
250 P2=GOVERNPR
260 GOSUB 880:H2SDASH=H2S:ISENENTHGOV=H2S  'Btu/lb
270 'LOCATE 9,10:INPUT "Enter the Enthalpy at inlet(Kcal/Kg)                :";H1:H1=H1*.454/.252
280 'LOCATE 10,10:INPUT "Enter the Isentropic Enthalpy at exhaust(Kcal/Kg)   :";H2S:H2S=H2S*.454/.252
290 'LOCATE 20,10:COLOR 10:PRINT "The Governor down stream pressure is :";GOVERNPR1;"Kg/cm2 a":COLOR 15
300 'LOCATE 11,10:PRINT "Enter the Isentropic Enthalpy based on  (Kcal/Kg)   :"
310 'LOCATE 12,10:PRINT "     governor d/s pressure."
320 'LOCATE 11,10:INPUT "Enter the Isentropic Enthalpy based on  (Kcal/Kg)   :";H2SDASH:H2SDASH=H2SDASH*.454/.252
330 'LOCATE 20,10:PRINT "                                                           "
340 LOCATE 13,10:INPUT "Enter the rotational speed     :";RPM
350 LOCATE 14,10:INPUT "Enter the wheel dia.(Outside) (m)  :";DIA:DIA1=DIA*100/2.54   'in.
360 LOCATE 15,10:INPUT "Enter the Power output required(MW):";POWER:POWER1=(1*1000000!/746)*POWER  'hp
362 'IF TURTYPE=1 OR TURTYPE=2 THEN 410 ELSE 560
370 LOCATE 16,10:INPUT "Whether it is a multistage/single stage turbine (M/S) :";TYPE$
380 IF TYPE$="m" OR TYPE$="M" THEN 560
390 IF TYPE$="s" OR TYPE$="S" THEN 410
400 GOTO 370
410 EFFCORR=.94:STAGES=1
420 BLADEVEL=RPM*(22/7)*DIA1/720 'ft/s
430 DELHSDASH=ENTHIN-H2SDASH
440 STEAMVEL=224*(DELHSDASH^.5)
450 VELRATIO=BLADEVEL/STEAMVEL
460 IF VELRATIO>=.2 THEN BASEEFF=.6:GOTO 490
470 IF VELRATIO>=9.000001E-02 THEN BASEEFF=1.359*(VELRATIO^.5080001):GOTO 490
480 BASEEFF=3.68*(VELRATIO^.922)
490 CORRBASEEFF=EFFCORR*BASEEFF
500 H2=ENTHIN-(CORRBASEEFF*(ENTHIN-ISENENTHOUT)):DUMMY=H2
510 TUREFF=.98*(ENTHIN-H2)/(H1-ISENENTHOUT)
520 STEAMRATE=2545/(TUREFF*(ENTHIN-ISENENTHOUT)) 'lbm/hp/h
530 STEAMFLOW=STEAMRATE*POWER1
540 GOTO 730
550 '
560 'Multistage
570 BLADEVEL=RPM*(22/7)*DIA1/720 'ft/s
580 DELHSDASHPERSTAGE=(BLADEVEL/100.8)^2
590 DELHSDASH=ENTHIN-H2SDASH
600 STAGES=CINT((ENTHIN-ISENENTHOUT)/DELHSDASHPERSTAGE)
610 IF STAGES<1 THEN STAGES=1
620 STEAMVEL=224*((DELHSDASH/STAGES)^.5)
630 VELRATIO=BLADEVEL/STEAMVEL
640 IF VELRATIO>=.4 THEN BASEEFF=.85:GOTO 680
650 IF VELRATIO>=.25 THEN BASEEFF=1.277*(VELRATIO^.444):GOTO 680
660 IF VELRATIO>=.125 THEN BASEEFF=2.055*(VELRATIO^.787):GOTO 680
670 BASEEFF=3.011*(VELRATIO^.971)
680 CORRBASEEFF=EFFCORR*BASEEFF
690 H2=ENTHIN-(CORRBASEEFF*(ENTHIN-ISENENTHOUT)):DUMMY=H2
700 TUREFF=.98*(ENTHIN-H2)/(ENTHIN-ISENENTHOUT)
710 STEAMRATE=2545/(TUREFF*(ENTHIN-ISENENTHOUT)) 'lbm/hp/h
720 STEAMFLOW=STEAMRATE*POWER1
730 'results
732 P1=PIN1+14.7:P2=POUT1:T1=(TIN*1.8)+32:EF=100*TUREFF
734 GOSUB 2000:TOUT=T2  'øF
740 CLS
750 LOCATE 2,25:PRINT "Turbine Efficiency Calculation Program"
760 LOCATE 4,15:PRINT "Number of stages                             :";STAGES
770 LOCATE 5,15:PRINT "Inlet Enthalpy (Kcal/Kg)                     :";ENTHIN*.252/.454
780 LOCATE 6,15:PRINT "Exhaust Isentropic Enthalpy(Kcal/Kg)         :";ISENENTHOUT*.252/.454
790 LOCATE 7,15:PRINT "Isentropic Enthalpy at governor d/s(Kcal/Kg) :";ISENENTHGOV*.252/.454
800 LOCATE 8,15:PRINT "Turbine efficiency (%)                       :";TUREFF*100
810 LOCATE 9,15:PRINT "Steam rate (Kg/(Kw)(h))                      :";STEAMRATE*.454/.746
820 LOCATE 10,15:PRINT "Steam flow rate (Kg/h)                     :";STEAMFLOW*.454
830 LOCATE 11,15:PRINT "Outlet Enthalpy (Kcal/Kg)                  :";DUMMY*.252/.454
840 LOCATE 12,15:PRINT "Governor d/s pr. (Kg/cm2 a)                :";GOVERNPR1
850 LOCATE 13,15:PRINT "Theor. Steam Velocity (ft/s)               :";STEAMVEL;"(";STEAMVEL*.3048;"m/s)"
860 LOCATE 14,15:PRINT "Blade velocity (ft/s)                      :";BLADEVEL;"(";BLADEVEL*.3048;"m/s)"
862 LOCATE 15,15:PRINT "Exhaust Temperature (øC)                   :";(TOUT-32)/1.8
863 LOCATE 16,15:PRINT "Exhaust Saturtated Temperature (øC)        :";(TSAT-32)/1.8
870 END
880 'Enthalpy/Properties calculations routine
890 'IF T1=0 THEN GOTO 980
900 P=P1:T=T1:GOSUB 1260
910 H1=Z:S1=SV:GOTO 940
920 'P=P1:GOSUB 1430
930 'H1=HV:S1=SV:T1=TSAT
940 P=P2:GOSUB 1370
950 IF S1>SV THEN 1000
960 X=(S1-SL)/(SV-SL):H2S=HV*X+(1-X)*HLIQ:RETURN
970 'H2=H1-.01*EF*(H1-H2S)
980 'IF H2>HV THEN 1060
990 'XF=(H2-HLIQ)/(HV-HLIQ):DELH=H1-H2:T2=TSAT:GOTO 1160
1000 T2=TSAT:RA=1000-TSAT
1010 P=P2:T=T2:GOSUB 1260
1020 S2C=SV:IF ABS((S1-S2C)/S1)<.002 THEN 1040
1030 RA=.5*RA:T2=T2+SGN(S1-S2C)*RA:GOTO 1010
1040 XF=1:T=T2:GOSUB 1260:H2S=Z:RETURN
1050 'H2=H1-.01*EF*(H1-H2S)
1060 'T2=TSAT:RA=1000-TSAT
1070 'P=P2:T=T2:GOSUB 1320
1080 'H2C=Z:IF ABS((H2-H2C)/H2)<.002 THEN 1160
1090 'RA=.5*RA:T2=T2+SGN(H2-H2C)*RA:GOTO 1130
1100 'ASR=3413/(H1-H2):TSR=3413/(H1-H2S)
1110 'CLS
1120 'LOCATE 2,15:PRINT "Turbine Steam Rates & Properties Calculation Program"
1130 'LOCATE 4,10:PRINT "Inlet press (psia)          :";P1
1140 'LOCATE 5,10:PRINT "Inlet Temp (F)              :";T1
1150 'LOCATE 6,10:PRINT "Inlet Enthalpy              :";H1
1160 'LOCATE 7,10:PRINT "Inlet Entropy               :";S1
1170 'LOCATE 8,10:PRINT "Exit pressure               :";P2
1180 'LOCATE 9,10:PRINT "Exit Temp                   :";T2
1190 'LOCATE 10,10:PRINT "Exit Enthalpy               :";H2
1200 'LOCATE 11,10:PRINT "Exit Saturated Temp         :";TSAT
1210 'LOCATE 12,10:PRINT "Actual steam req. (ASR)     :";ASR
1220 'LOCATE 13,10:PRINT "Theor. steam req. (TSR)     :";TSR
1230 'LOCATE 14,10:PRINT "Expansion efficiency        :";EF
1240 'LOCATE 15,10:PRINT "Exit quality                :";XF
1250 'END
1260 T=273.1+(T-32)/1.8:P=P/14.696
1270 K4=80870!/T/T:L1=10^K4*(-2641.62/T):M1=1.89+L1:N1=M1*P*P/T/T
1280 O1=2+(372420!/T/T):Q1=O1*L1:R1=1.89+Q1:U1=(.21828*T-126970!/T):V4=2*U1*R1-(M1/T)*126970!
1290 W1=82.54-162460!/T:Y1=2*W1*R1-(M1/T)*162460!
1300 F1=(((U1*M1*N1+W1)*N1/P+1)*M1+4.5504*T/P)*.0160185
1310 Z=775.6+.63296*T+1.62467E-04*T*T+20.5697*LOG(T)
1320 Z=Z+.043557*(R1*P+.5*N1*(Y1+M1*(W1+V4*N1)))
1330 SV=(((U1*M1-2*V4)*.5*M1*N1-Y1)*.5*N1+(M1-R1)*P)/T*(-.0241983)-.355579-11.4276/T
1340 SV=SV+1.8052E-04*T-.11022*LOG(P)+.35164*LOG(T)
1350 P=14.696*P:T=(T-273.15)*1.8+32
1360 RETURN
1370 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
1380 TSAT=AA(1):HLIQ=AA(5):HV=AA(6):SL=AA(7):SV=AA(8)
1390 RETURN
1400 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
1410 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
1420 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
1430 DATA -1.476933e-4,1.2617946e-3,3.44201e-3,-.08494128,6.89138e-8,-2.4941e-11,1.97364
2000 '
2150 P=P1:T=T1:GOSUB 2510
2160 H1=Z:S1=SV:GOTO 2190
2190 P=P2:GOSUB 2620
2200 IF S1>SV THEN 2250
2210 X=(S1-SL)/(SV-SL):H2S=HV*X+(1-X)*HLIQ
2220 H2=H1-.01*EF*(H1-H2S)
2230 IF H2>HV THEN 2250
2240 XF=(H2-HLIQ)/(HV-HLIQ):DELH=H1-H2:T2=TSAT:RETURN
2250 T2=TSAT:RA=1000-TSAT
2260 P=P2:T=T2:GOSUB 2510
2270 S2C=SV:IF ABS((S1-S2C)/S1)<.002 THEN 2290
2280 RA=.5*RA:T2=T2+SGN(S1-S2C)*RA:GOTO 2260
2290 XF=1:T=T2:GOSUB 2510:H2S=Z
2300 H2=H1-.01*EF*(H1-H2S)
2310 T2=TSAT:RA=1000-TSAT
2320 P=P2:T=T2:GOSUB 2510
2330 H2C=Z:IF ABS((H2-H2C)/H2)<.002 THEN 2350
2340 RA=.5*RA:T2=T2+SGN(H2-H2C)*RA:GOTO 2320
2350 RETURN
2360 CLS
2510 T=273.1+(T-32)/1.8:P=P/14.696
2520 K4=80870!/T/T:L1=10^K4*(-2641.62/T):M1=1.89+L1:N1=M1*P*P/T/T
2530 O1=2+(372420!/T/T):Q1=O1*L1:R1=1.89+Q1:U1=(.21828*T-126970!/T):V4=2*U1*R1-(M1/T)*126970!
2540 W1=82.54-162460!/T:Y1=2*W1*R1-(M1/T)*162460!
2550 F1=(((U1*M1*N1+W1)*N1/P+1)*M1+4.5504*T/P)*.0160185
2560 Z=775.6+.63296*T+1.62467E-04*T*T+20.5697*LOG(T)
2570 Z=Z+.043557*(R1*P+.5*N1*(Y1+M1*(W1+V4*N1)))
2580 SV=(((U1*M1-2*V4)*.5*M1*N1-Y1)*.5*N1+(M1-R1)*P)/T*(-.0241983)-.355579-11.4276/T
2590 SV=SV+1.8052E-04*T-.11022*LOG(P)+.35164*LOG(T)
2600 P=14.696*P:T=(T-273.15)*1.8+32
2610 RETURN
2620 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
2630 TSAT=AA(1):HLIQ=AA(5):HV=AA(6):SL=AA(7):SV=AA(8)
2640 RETURN