blob: 89200fe6a6d5b3541ee2c04d752bdb4b0cee0df6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
clc
P1 = 0.5 // Initial pressure in MPa
V1 = 0.2 // Initial volume in m^3
V2 = 0.05 // Final volume in m^3
n = 1.3 // Polytropic index
printf("\n Example 7.7")
P2 = P1*(V1/V2)^n
function y = f(p)
y = ((P1*V1^n)/p)^(1/n)
endfunction
H = integrate('f','p',P1,P2) // H = H2-H1
U = H-(P2*V2-P1*V1)
W12 = -U
printf("\n Change in enthalpy is %f kJ",H*1e3)
printf("\n Change in internal energy is %f kJ",U*1000)
printf("\n The change in entropy and heat transfer are is %d kJ",0)
printf("\n The work transfer during the process is %f kJ",W12*1000)
//The answers vary due to round off error
|
