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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
////
//Variable declaration
n = 2.5 //moles of ideal gas
R = 0.08314 //Ideal gas constant, bar.L/(mol.K)
cvm = 20.79 //Heat Capacity at constant volume, J/(mol.K)
p1 = 16.6 //Pressure at point 1, bar
v1 = 1.00 //Volume at point 1, L
p2 = 16.6 //Pressure at point 2, bar
v2 = 25.0 //Volume at point 2, L
v3 = 25.0 //Volume at point 3, L
//Calculations
T1 = p1*v1/(n*R)
T2 = p2*v2/(n*R)
T3 = T1 //from problem statement
//for path 1-2
DU12 = n*cvm*(T2-T1)
w12 = -p1*1e5*(v2-v1)*1e-3
q12 = DU12 - w12
DH12 = DU12 + n*R*(T2-T1)*1e2
//for path 2-3
w23 = 0.0
DU23=n*cvm*(T3-T2)
;q23=n*cvm*(T3-T2)
;
DH23 = -DH12
//for path 3-1
DU31 = 0.0 //Isothemal process
DH31 = 0.0
w31 = -n*R*1e2*T1*log(v1/v3)
q31 = -w31
DU = DU12+DU23+DU31
w = w12+w23+w31
q = q12+q23+q31
DH = DH12+DH23+DH31
//Results
printf("\n For Path q w DU DH ')
printf("\n 1-2 %7.2f %7.2f %7.2f %7.2f",q12,w12,DU12,DH12)
printf("\n 2-3 %7.2f %7.2f %7.2f %7.2f",q23,w23,DU23,DH23)
printf("\n 3-1 %7.2f %7.2f %7.2f %7.2f",q31,w31,DU31,DH31)
printf("\n Overall %7.2f %7.2f %7.2f %7.2f",q,w,DU,DH)
printf("\n all values are in J')
|
