blob: 3f01fd7d52fe20c4280a58e42a2f66e5bff2269f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
clc
T1=300 //temperature in kelvin
P1=100 //pressure in kPa
P2=2 //pressure in MPa
gama=1.4 //Cp/Cv ratio
T2=T1*(((P2*10^6)/(P1*10^3))^((gama-1)/gama))
mprintf("T2=%fK\n",T2)//ans vary due to roundoff error
Cr=(T1*P2*10^6)/(P1*10^3*T2)
mprintf("Compression ratio=%f\n",Cr)//ans vary due to roundoff error
R=8.314
W=R*(T1-T2)/(gama-1)
mprintf("W=%fkJ/mol",W/1000)//ans vary due to roundoff error
|