summaryrefslogtreecommitdiff
path: root/3685/CH17/EX17.4/Ex17_4.sce
blob: 0774f53336a5fca100b008a08c72390b4fcf3160 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
clc
Px = 16 // pressure in kPa
Poy = 70 //pressure in kPa 
Mx = 1.735 // Mach number
Pyx = 3.34 // Pressure ratio
rho_yx = 2.25 // Density ratio
Tyx = 1.483  // Temperature ratio
Poyox = 0.84 // pressure ratio
My = 0.631 // Mach number
g = 1.4 // Ratio of heat capacities
Tox = 573 // stagnation temperature in K 
Toy = Tox // temperature equivalence
Tx = Tox/(1+((g-1)/2)*Mx^2) // temperature at x
Ty = Tyx*Tx // temperature at y
Pox = Poy/Poyox  // total pressure 
// From table
Mx = 1.735

printf("\n Example 17.4\n")
printf("\n Mach number of the tunnel is %f ",Mx)