blob: 2c1a4a17874d8941bb2cf509f024be718baec51a (
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
27
28
29
|
// Exa 3.24
clc;
clear;
close;
// Given data
C_P = 1.005;// in kJ/kg-K
C_V = 0.718;// in kJ/kg-K
R = C_P-C_V;// in kJ/kg-K
P1 = 20;//in bar
P2 = 12;// in bar
T1 = 200;//in degree C
T1 = T1 + 273;// in K
T2 = 125;//in degree c
T2 = T2 + 273;// in K
V1 = (R*10^3*T1)/(P1*10^5);// in m^3
V2 = (R*10^3*T2)/(P2*10^5);// in m^3
W = 10^5 * integrate('-293*V + 40','V',0.0679,0.0952);//in Joules
W = round(W * 10^-3);// in kJ
disp(W,"Work done in kJ is");
m = 1;// in kg
del_U = m*C_V*(T2-T1);//change in internal energy in kJ
disp(del_U,"Change in internal energy in kJ is");
disp("Negative sign indicates that there is decrease in internal energy of the gas. ")
C_Enthalpy = m*C_P*(T2-T1);//change in enthalpy in kJ
disp(C_Enthalpy,"The change in enthalpy in kJ is :")
disp("Negative sign indicates that there is decrease in enthalpy of the gas")
Q = W+ del_U;// in kJ
disp(Q,"Heat transfer in kJ is");
disp("Negative sign indicates that the heat is rejected by the air")
|
