blob: 5f8edaff2d77a8d8b30b1b389ddbdc7bd8a8e82c (
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
30
31
32
33
34
35
36
|
clc
r=10; //kg/min
p1=1.5*10^5; //N/m^2
rho1=26; //kg/m^3
C1=110; //m/s
u1=910; //kJ/kg
p2=5.5*10^5; //N/m^2
rho2=5.5; //kg/m^3
C2=190; //m/s
u2=710; //kJ/kg
Q=55; //kJ/s
h=55; //m
g=9.81; //m/s^2
v2=1/rho2;
v1=1/rho1;
disp("(i) Change in enthalpy")
dh=u2-u1+ (p2*v2-p1*v1)/10^3;
disp(dh)
disp("kJ/kg")
disp("(ii) Work done during the process (W).")
Q=330; //kJ/kg
KE=(C2^2-C1^2)/2/10^3; //kJ
PE=g*h/10^3; //kJ
W=-Q-KE-PE-dh;
disp("Work done = ")
disp(W)
disp("kJ")
disp("Work done per second = ")
P=W*10/60;
disp(P)
disp("kW")
|
