summaryrefslogtreecommitdiff
path: root/3765/CH5/EX5.7/Ex5_7.sce
blob: a0ffe3e57c3d533fffd587aebd7eaad3c6e1db9c (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
26
clc
// Example 5.7.py
// Consider the wind tunnel described in example 5.2. Estimate the ratio of diffuser
// throat area to nozzle throat area required to allow the tunnel to start. Also, 
// assuming that the diffuser efficiency is 1.2 after the tunnel has started, calculate
// the pressure ratio across the tunnel necessary for running i.e. calculate the ratio
// of total pressure at the diffuser exit to the reservoir pressure.

// Variable declaration

M = 2.5      // mach number before the shock
eta_d = 1.2  // diffuser efficiency

// Calculations

// from table for M = 2.5
po2_by_po1 = 0.499              // po2/po1
At2_by_At1 = 1 / po2_by_po1     // At2/At1 = po1/po2
 
Pdo_by_po = eta_d * po2_by_po1  // pdo/po

// Results
printf("\n Ratio of diffuser throat area to nozzle throat area %.2f", At2_by_At1)

printf("\n Ratio of total pressure at the diffuser exit to the reservoir pressure, %.3f",(Pdo_by_po))