summaryrefslogtreecommitdiff
path: root/779/CH14
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /779/CH14
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '779/CH14')
-rwxr-xr-x779/CH14/EX14.1/14_1.sce6
-rwxr-xr-x779/CH14/EX14.2/14_2.sce12
-rwxr-xr-x779/CH14/EX14.3/14_3.sce29
-rwxr-xr-x779/CH14/EX14.4/14_4.sce24
-rwxr-xr-x779/CH14/EX14.5/14_5.sce20
-rwxr-xr-x779/CH14/EX14.6/14_6.sce12
-rwxr-xr-x779/CH14/EX14.7/14_7.sce14
-rwxr-xr-x779/CH14/EX14.8/14_8.sce17
-rwxr-xr-x779/CH14/EX14.9/14_9.sce15
9 files changed, 149 insertions, 0 deletions
diff --git a/779/CH14/EX14.1/14_1.sce b/779/CH14/EX14.1/14_1.sce
new file mode 100755
index 000000000..419a01074
--- /dev/null
+++ b/779/CH14/EX14.1/14_1.sce
@@ -0,0 +1,6 @@
+T2 = 268; T1 = 308;
+COP = T2/(T1-T2);
+ACOP = COP/3; // Actual COP
+Q2 = 29; // in kW
+W = Q2/ACOP;
+disp("kW",W,"Power required to derive the plane is") \ No newline at end of file
diff --git a/779/CH14/EX14.2/14_2.sce b/779/CH14/EX14.2/14_2.sce
new file mode 100755
index 000000000..a7fc63382
--- /dev/null
+++ b/779/CH14/EX14.2/14_2.sce
@@ -0,0 +1,12 @@
+h1 = 236.04; s1 = 0.9322; s2 = s1;
+P2 = 0.8; // in MPa
+h2 = 272.05; h3 = 93.42; h4 = h3;
+m = 0.06; // mass flow rate
+Q2 = m*(h1-h4);
+Wc = m*(h2-h1);
+Q1 = m*(h2-h4);
+COP = Q2/Wc;
+disp("kW",Q2,"The rate of heat removal is")
+disp("kW",Wc,"Power input to the compressor is")
+disp("kW",Q1,"The heat rejection rate in the condenser is")
+disp(COP,"COP is")
diff --git a/779/CH14/EX14.3/14_3.sce b/779/CH14/EX14.3/14_3.sce
new file mode 100755
index 000000000..65103e935
--- /dev/null
+++ b/779/CH14/EX14.3/14_3.sce
@@ -0,0 +1,29 @@
+h1 = 183.19; h2 = 209.41; h3 = 74.59; h4 = h3;
+T1 = 313; T2 = 263;
+W = 70000/3600; // Plant capacity in kW
+w = W/(h1-h4); // Refrigerant flow rate
+v1 = 0.077;
+VFR = w*v1;
+T = 48; // in degree
+P2 = 9.6066; P1 = 2.1912;
+rp = P2/P1; // Pressure ratio
+Q1 = w*(h2-h3);
+hf = 26.87; hfg = 156.31;
+x4 = (h4-hf)/hfg;
+COP = (h1-h4)/(h2-h1);
+PI = w*(h2-h1);
+COP = T2/(T1-T2);
+COP_v = 4.14;
+r = COP_v/COP;
+disp("kg/s",w,"Refrigerant flow rate is")
+disp("m3/s",VFR,"Volume flow rate is")
+disp("degree",T,"Compressor discharge temperature is")
+disp(rp,"Pressure ratio is")
+disp("kW",Q1,"Heat rejected to the condenser is")
+disp("%",x4*100,"Flash gas percentage is")
+disp(COP,"COP is")
+disp("kW",PI,"Power required to drive the compressor is")
+disp(r,"Ratio of COP of carnot refrigerator is")
+
+
+
diff --git a/779/CH14/EX14.4/14_4.sce b/779/CH14/EX14.4/14_4.sce
new file mode 100755
index 000000000..ee7ec92d3
--- /dev/null
+++ b/779/CH14/EX14.4/14_4.sce
@@ -0,0 +1,24 @@
+h3 = 882; h2 = 1034;
+h6 = 998; h1 = 1008;
+v1 = 0.084;
+h4 = h3-h1+h6; h5 = h4;
+t4 = 25+273;
+disp("kJ/kg",h6-h5,"Refrigeration effect is")
+m = 10;
+w = (m*14000)/((h6-h5)*3600); // in kg/s
+disp("kg/s",w,"Refrigerant flow rate is")
+v1 = 0.084;
+VFR = w*3600*v1; // in kg/h
+ve = 0.8; // volumetric efficiency
+CD = VFR/(ve*60); // in m3/min
+N = 900;
+n = 2;
+D = ((CD*4)/(%pi*1.1*N*n))^(1/3); // L = 1.1D L = length D = diameter
+L = 1.1*D;
+disp("cm",D*100,"Diameter of cylinder is")
+disp("cm",L*100,"Length of cylinder is")
+COP = (h6-h5)/(h2-h1);
+PI = w*(h2-h1);
+disp("kW",PI,"Power required to drive the compresor is")
+disp(COP,"COP is")
+
diff --git a/779/CH14/EX14.5/14_5.sce b/779/CH14/EX14.5/14_5.sce
new file mode 100755
index 000000000..a20684c20
--- /dev/null
+++ b/779/CH14/EX14.5/14_5.sce
@@ -0,0 +1,20 @@
+P2 = 1554.3;
+P1 = 119.5;
+Pi = sqrt(P1*P2);
+h1 = 1404.6; h2 = 1574.3; h3 = 1443.5; h4 = 1628.1;
+h5 = 371.7; h6 = h5; h7 = 181.5;
+w = 30; // capacity of plant
+m2_dot = (3.89*30)/(h1-h7);
+m1_dot = m2_dot*((h2-h7)/(h3-h6));
+Wc_dot = m2_dot*(h2-h1)+m1_dot*(h4-h3);
+COP = w*3.89/Wc_dot;
+// single stage
+h1_ = 1404.6; h2_ = 1805.1;
+h3_ = 371.1; h4_ = h3_;
+m_dot = (3.89*30)/(h1_-h4_);
+Wc = m_dot*(h2_-h1_);
+COP_ = w*3.89/Wc;
+IW = (Wc-Wc_dot)/Wc_dot;
+ICOP = (COP-COP_)/COP_
+disp("%",IW*100,"Increase in work of compression is")
+disp("%",ICOP*100,"Increase in COP for 2 stage compression is") \ No newline at end of file
diff --git a/779/CH14/EX14.6/14_6.sce b/779/CH14/EX14.6/14_6.sce
new file mode 100755
index 000000000..14fb2b0d5
--- /dev/null
+++ b/779/CH14/EX14.6/14_6.sce
@@ -0,0 +1,12 @@
+tsat = 120.2+273; hfg = 2201.9;
+T1 = 120.2+273;
+T2 = 30+273;
+Tr = -10+273;
+COP_max = ((T1-T2)*Tr)/((T2-Tr)*T1);
+ACOP = 0.4*COP_max;
+Qe = (20*14000)/3600; // in KW
+Qg = Qe/ACOP;
+x = 0.9;
+H = x*hfg;
+SFR = Qg/H;
+disp("kg/s",SFR,"Steam flow rate required is") \ No newline at end of file
diff --git a/779/CH14/EX14.7/14_7.sce b/779/CH14/EX14.7/14_7.sce
new file mode 100755
index 000000000..c9c59465a
--- /dev/null
+++ b/779/CH14/EX14.7/14_7.sce
@@ -0,0 +1,14 @@
+T1 = 277; T3 = 273+55;
+rp = 3; // Pressure ratio
+g = 1.4; cp = 1.005;
+T2s = T1*(rp^((g-1)/g));
+T2 = T1+(T2s-T1)/0.72
+T4s = T3/(rp^((g-1)/g));
+T34 = 0.78*(T3-T4s); // T3-T4
+T4 = T3-T34;
+COP = (T1-T4)/((T2-T1)-(T3-T4));
+disp(COP,"COP of the refrigerator is")
+P = (3*14000)/(COP*3600)
+disp("kW",P,"Driving power required is")
+m = (3*14000)/(cp*(T1-T4));
+disp("kg/s",m/3600,"Mass flow rate is")
diff --git a/779/CH14/EX14.8/14_8.sce b/779/CH14/EX14.8/14_8.sce
new file mode 100755
index 000000000..8f1436c77
--- /dev/null
+++ b/779/CH14/EX14.8/14_8.sce
@@ -0,0 +1,17 @@
+P1 = 2.4; T1 = 0+273;
+h1 = 188.9; s1 = 0.7177; v1 = 0.0703;
+P2 = 9; T2 = 60+273;
+h2 = 219.37;
+h2s = 213.27;
+h3 = 71.93; h4 = h3;
+v1 = 0.0703;
+A1V1 = 0.6/60;
+m_dot = A1V1/0.0703;
+Wc_dot = m_dot*(h2-h1);
+Q1_dot = m_dot*(h2-h3);
+COP = Q1_dot/Wc_dot;
+nis = (h2s-h1)/(h2-h1);
+disp("kW",Wc_dot,"Power input is")
+disp("kW",Q1_dot,"Heating capacity is")
+disp(COP,"COP is")
+disp("%",nis*100,"The isentropic compressor efficiency is") \ No newline at end of file
diff --git a/779/CH14/EX14.9/14_9.sce b/779/CH14/EX14.9/14_9.sce
new file mode 100755
index 000000000..7a4879378
--- /dev/null
+++ b/779/CH14/EX14.9/14_9.sce
@@ -0,0 +1,15 @@
+T1 = 275; T3 = 310;
+P1 = 1 ; P2 = 4;
+T2s = T1*(P2/P1);
+nc = 0.8;
+T2 = T1 + (T2s-T1)*nc;
+pr = 0.1;
+P3 = P2-0.1;
+P4 = P1+0.08;
+PR = P3/P4;
+disp(PR,"Pressure ratio for the turbine is")
+T4s = T3*(1/PR)^(0.286);
+nt = 0.85;
+T4 = T3-(T3-T4s)*nt;
+COP = (T1-T4)/((T2-T3)-(T1-T4));
+disp(COP,"COP is")