blob: 9635931a783835f4d8fbad501aba09e047bba946 (
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
37
38
|
clc
clear
//INPUT DATA
hf1=100.4;//enthalpy in kJ/kg
hf2=-54.56;//enthalpy in kJ/kg
hg1=1319.22;//enthalpy in kJ/kg
hg2=1304.99;//enthalpy in kJ/kg
sf1=0.3473;//entropy in kJ/kgK
sf2=-0.2134;//entropy in kJ/kgK
sg1=4.4852;//entropy in kJ/kgK
sg2=5.0585;//entropy in kJ/kgK
t1=20;//temperature in Degree C
t2=-15;//temperature in Degree C
mi=30;//mass of ice
hfgi=335;//enthalpy in kJ/kg
cpw=4.1868;//specific pressure of water in kJ/kg K
//CALCULATIONS
x1=((sg1-sf2)/(sg2-sf2));//fraction
h1=hf2+x1*(hg2-hf2);//enthalpy in kJ/kg
copt=((h1-hf1)/(hg1-h1));//COP of the system
copa=copt*0.6;//actual cop
Qa=mi*10^3*(cpw*t1+hfgi)/(24*3600);//heat removed from water in kJ/s
w=Qa/copa;//Power required to drive compressor in kW
//OUTPUT
printf('power required to drive the compresor i %3.2f kW',w)
|
