blob: 4572106523f36c4bc7ce9bfe59c54ee8c86725cc (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
clc
clear
//input data
M1=0.25 //Initial Mach number
M2=0.75 //Final mach number
P1=1.5 //Inlet pressure in bar
T1=300 //Inlet temperature in K
k=1.4 //Adiabatic constant
R=0.287 //Gas constant in kJ/kg-K
//calculation
a1=sqrt(k*R*10^3*T1) //Sound velocity in m/s, R in J/kg
C1=a1*M1 //air velocity in m/s
p1=4.3615 //Pressure ratio at entry from gas tables @M1,k
Pt=P1/p1 //Static critical pressure in kPa
c1=0.272 //Velocity ratio from gas tables (fanno flow tables,k=1.4,M1)
Ct=C1/c1 //Critical velocity in m/s
p2=1.385 //Pressure ratio at exit from gas tables @M2,k
P2=Pt*p2 //Exit pressure in bar
c2=0.779 //Velocity ratio at exit from gas tables @M2,k
C2=Ct*c2 //Exit velocity in m/s
//output
printf('Final pressure and velocity are %3.4f bar and %3.2f m/s',P2,C2)
|
