CLS
DIM i AS INTEGER
DIM j AS INTEGER
DIM p(1000)
DIM t(1000)
DIM v(1000)
DIM vs(1000)
DIM m(1000)
OPEN "Output.txt" FOR OUTPUT AS #1
tin = 298
pin = 1
deltav = 1
kappa = 1.4
cp = 1004
r = 286.7
PRINT ""
PRINT "                             Diffuser  Calculation"
PRINT ""
PRINT "                             Enter the Input Data!"
PRINT "                          Followed by  the enter key!"
PRINT ""
INPUT " Static inlet temperature (default value is 298 [K])! ", tin
PRINT ""
INPUT " Static inlet pressure (default value is 1 [atm])! ", pin
PRINT ""
INPUT " Inlet Mach number (default value is 0.85 [/])! ", mout
PRINT ""
INPUT " Velocity delta step (default value is 10 [m/s])! ", deltav
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 ""
INPUT " Working fluid gas constant (default value is 286.7 [J/kg*K])! ", r
PRINT ""
PRINT " Step    T        p           v          Mach No."
PRINT " [/]     [K]      [atm]       [m/s]      [/]   "
PRINT ""
IF tin = 0 THEN tin = 298
IF pin = 0 THEN pin = 1
IF mout = 0 THEN mout = .85
IF deltav = 0 THEN deltav = 10
IF kappa = 0 THEN kappa = 1.4
IF cp = 0 THEN cp = 1004
IF r = 0 THEN r = 286.7
PRINT #1, " Static inlet temperature [K]: ", TAB(44); USING "#,###"; tin;
PRINT #1, ""
PRINT #1, " Static inlet pressure [atm]: ", TAB(45); pin;
PRINT #1, ""
PRINT #1, " Inlet Mach Number [/]: ", TAB(45); USING "##.##"; mout;
PRINT #1, ""
PRINT #1, " Velocity delta step [m/s]: ", TAB(45); deltav;
PRINT #1, ""
PRINT #1, " Woorking fluid kappa [/]: ", TAB(45); kappa;
PRINT #1, ""
PRINT #1, " Working fluid cp [J/kg*K]: "; TAB(45); USING "##,###"; cp;
PRINT #1, ""
PRINT #1, " Working fluid gas constant [J/kg*K]: "; TAB(45); r
PRINT #1, ""
PRINT #1, ""
PRINT #1, " Step    Tin     pin      Tout    pout        v          Mach No."
PRINT #1, " [/]     [K]     [atm]    [K]     [atm]       [m/s]      [/]     "
PRINT #1, ""
p(1) = pin
t(1) = tin
v(1) = 0
vs(1) = 0
m(1) = 0
FOR i = 1 TO 1000
IF i > 1000 THEN GOTO finish:
IF m(i) > mout THEN GOTO test:
PRINT TAB(1); i; TAB(8); USING "#,###"; t(i);
PRINT TAB(18); USING "##.#####"; p(i);
PRINT TAB(28); USING "##,###.##"; v(i); TAB(42);
PRINT USING "#.#####"; m(i);
test:
IF (i + 1) > 1000 THEN GOTO finish:
v(i + 1) = v(i) + deltav
IF m(i) > mout THEN GOTO finish:
t(i + 1) = tin + (v(i + 1) ^ 2) / (2 * cp)
p(i + 1) = p(i) * ((t(i + 1) / t(i)) ^ (kappa / (kappa - 1)))
vs(i + 1) = (kappa * r * tin) ^ .5
m(i + 1) = v(i + 1) / vs(i + 1)
PRINT #1, TAB(1); i; TAB(8); USING "#,###"; tin;
PRINT #1, TAB(17); USING "##.##"; pin;
PRINT #1, TAB(24); USING "##,###"; t(i);
PRINT #1, TAB(34); USING "##.#####"; p(i);
PRINT #1, TAB(44); USING "##,###.##"; v(i); TAB(58);
PRINT #1, USING "#.#####"; m(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