blob: 9dfa644533b60a8ff569c5d5f0cf7ef074e42338 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
clc
//initialisation of variables
clear
r= 1.4
R= 53.3 //ft lbf/lbm R
g= 32.2 //ft/sec^2
T1= 410 //R
v= 2500 //ft/sec
P1= 628 //lbf/in^2
//CALCULATIONS
v1= sqrt(r*g*R*T1)
Ma1= v/v1
Ts1= T1*(1+0.5*(r-1)*Ma1^2)
Ps1= P1*(1+0.5*(r-1)*Ma1^2)^(r/(r-1))
Ps2= Ps1*((r+1)/(2*r*Ma1^2-r+1))^(1/(r-1))*(0.5*(r+1)*Ma1^2/(1+0.5*(r-1)*Ma1^2))^(r/(r-1))
//RESULTS
printf ('acoustic velocity = %.f ft/sec',v1)
printf ('\n Match number = %.2f ',Ma1)
printf ('\n Stagnition temperature = %.f R',Ts1)
printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps1)
printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps2)
|
