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 = 1500
pin = 25
pout = 1
deltap = 1
kappa = 1.4
cp = 1004
PRINT ""
PRINT "                             Expansion Calculation"
PRINT ""
PRINT "                             Enter the Input Data!"
PRINT "                          Followed by  the enter key!"
PRINT ""
INPUT " Expansion inlet temperature (default value is 1,500 [K])! ", tin
PRINT ""
INPUT " Expansion inlet pressure (default value is 25 [atm])! ", pin
PRINT ""
INPUT " Expansion outlet pressure (default value is 1 [atm])! ", pout
PRINT ""
INPUT " Expansion 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 = 1500
IF pin = 0 THEN pin = 25
IF pout = 0 THEN pout = 1
IF deltap = 0 THEN deltap = 1
IF kappa = 0 THEN kappa = 1.4
IF cp = 0 THEN cp = 1004
PRINT #1, " Expansion inlet temperature [K]: ", TAB(45); USING "##,###"; tin;
PRINT #1, ""
PRINT #1, " Expansion inlet pressure [atm]: ", TAB(45); pin;
PRINT #1, ""
PRINT #1, " Expansion outlet pressure [atm]: ", TAB(45); pout;
PRINT #1, ""
PRINT #1, " Expansion 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(10); USING "#,###"; tin;
PRINT TAB(19); pin; TAB(29); USING "##,###.##"; t2(i); TAB(41);
PRINT p2(i); TAB(48); 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 * (tin - t2(i + 1)) / 1000
PRINT #1, TAB(1); i; TAB(9); tin; TAB(18); pin; TAB(29); USING "##,###.##"; t2(i); TAB(41);
PRINT #1, p2(i); TAB(50); 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