blob: aecc14046a5609364473f1275c4d49ccd1a07ae1 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
clc;
clear;
p1=14.7;//psi(abs)
V1=1;//ft^3
V2=0.5;//ft^3
//for isentropic compression, (p1(d1^k))=(p2/(d2^k))
//volume*density=constant(mass)
ratd=V1/V2;
p2=((ratd)^1.66)*p1;//psi(abs)
disp("psi(abs)",p2,"final pressure p2=")
i=1;
ratV=0.01:0.01:1.0;
for j=0.01:0.01:1.0
pres(i)=p1/((j)^1.66);
i=i+1;
end
plot2d(ratV,pres,rect=[0,0,1,1000])
xtitle('p2 vs V2/V1','V2/V1','p2 psi')
|
