CLS
DIM i AS INTEGER
DIM j AS INTEGER
DIM p2(1000)
DIM t2(1000)
DIM w(1000)
OPEN "Output.txt" FOR OUTPUT AS #1
tin = 298
pin = 1
pout = 25
deltap = 1
kappa = 1.4
cp = 1004
PRINT ""
PRINT "                            Compression Calculation"
PRINT ""
PRINT "                             Enter the Input Data!"
PRINT "                          Followed by  the enter key!"
PRINT ""
INPUT " Compression inlet temperature (default value is 298 [K])! ", tin
PRINT ""
INPUT " Compression inlet pressure (default value is 1 [atm])! ", pin
PRINT ""
INPUT " Compression outlet pressure (default value is 25 [atm])! ", pout
PRINT ""
INPUT " Compression delta pressure step (default value is 1 [atm])! ", deltap
PRINT ""
INPUT " Working fluid kappa (default value is 1.4 [/])! ", kappa
PRINT ""
INPUT " Working fluid cp (default value is 1,004 [J/kg*K])! ", cp
PRINT ""
PRINT " Step    T1      p1       T2        p2       w        "
PRINT " [/]     [K]     [atm]    [K]       [atm]    [kW/kg/s]"
PRINT ""
IF tin = 0 THEN tin = 298
IF pin = 0 THEN pin = 1
IF pout = 0 THEN pout = 25
IF deltap = 0 THEN deltap = 1
IF kappa = 0 THEN kappa = 1.4
IF cp = 0 THEN cp = 1004
PRINT #1, " Compression inlet temperature [K]: ", TAB(45); tin;
PRINT #1, ""
PRINT #1, " Compression inlet pressure [atm]: ", TAB(45); pin;
PRINT #1, ""
PRINT #1, " Compression outlet pressure [atm]: ", TAB(45); pout;
PRINT #1, ""
PRINT #1, " Compression delta pressure step [atm]: ", TAB(45); deltap;
PRINT #1, ""
PRINT #1, " Woorking fluid kappa [/]: ", TAB(45); kappa;
PRINT #1, ""
PRINT #1, " Working fluid cp [J/kg*K]: "; TAB(45); cp;
PRINT #1, ""
PRINT #1, ""
PRINT #1, ""
PRINT #1, " Step    Tin      pin      Tout      pout      w        "
PRINT #1, " [/]     [K]      [atm]    [K]       [atm]     [kW/kg/s]"
PRINT #1, ""
p2(1) = pin
t2(1) = tin * (pin / pin) ^ ((kappa - 1) / kappa)
w(1) = cp * (tin - tin) / 1000
FOR i = 1 TO 1000
IF p2(i) > pout THEN GOTO test:
PRINT TAB(1); i; TAB(9); tin; TAB(17); pin; TAB(24); USING "##,###.##"; t2(i); TAB(36);
PRINT p2(i); TAB(42); USING "###,###.##"; w(i);
test:
IF (i + 1) > 1000 THEN GOTO finish:
p2(i + 1) = p2(i) + deltap
IF p2(i) > pout THEN GOTO finish:
t2(i + 1) = tin * (p2(i + 1) / pin) ^ ((kappa - 1) / kappa)
w(i + 1) = cp * (t2(i + 1) - tin) / 1000
PRINT #1, TAB(1); i; TAB(9); tin; TAB(18); pin; TAB(25); USING "##,###.##"; t2(i); TAB(37);
PRINT #1, p2(i); TAB(44); USING "###,###.##"; w(i);
NEXT
finish:
PRINT ""
PRINT ""
PRINT "                    Output.txt file contains the output data!"
PRINT "                                The Run is Over!"
PRINT "                         Hit anykey to close the program!"
PRINT ""
PRINT "                               Engineeering Software"
PRINT "                          http://wwww.engineering-4e.com"
CLOSE #1
DO
LOOP UNTIL INKEY$ <> ""
END