diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /779/CH14 | |
download | Scilab-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-x | 779/CH14/EX14.1/14_1.sce | 6 | ||||
-rwxr-xr-x | 779/CH14/EX14.2/14_2.sce | 12 | ||||
-rwxr-xr-x | 779/CH14/EX14.3/14_3.sce | 29 | ||||
-rwxr-xr-x | 779/CH14/EX14.4/14_4.sce | 24 | ||||
-rwxr-xr-x | 779/CH14/EX14.5/14_5.sce | 20 | ||||
-rwxr-xr-x | 779/CH14/EX14.6/14_6.sce | 12 | ||||
-rwxr-xr-x | 779/CH14/EX14.7/14_7.sce | 14 | ||||
-rwxr-xr-x | 779/CH14/EX14.8/14_8.sce | 17 | ||||
-rwxr-xr-x | 779/CH14/EX14.9/14_9.sce | 15 |
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")
|