summaryrefslogtreecommitdiff
path: root/551/CH14
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /551/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 '551/CH14')
-rwxr-xr-x551/CH14/EX14.1/1.sce27
-rwxr-xr-x551/CH14/EX14.10/10.sce48
-rwxr-xr-x551/CH14/EX14.11/11.sce55
-rwxr-xr-x551/CH14/EX14.12/12.sce30
-rwxr-xr-x551/CH14/EX14.13/13.sce25
-rwxr-xr-x551/CH14/EX14.14/14.sce31
-rwxr-xr-x551/CH14/EX14.15/15.sce30
-rwxr-xr-x551/CH14/EX14.16/16.sce129
-rwxr-xr-x551/CH14/EX14.17/17.sce28
-rwxr-xr-x551/CH14/EX14.18/18.sce22
-rwxr-xr-x551/CH14/EX14.19/19.sce33
-rwxr-xr-x551/CH14/EX14.2/2.sce22
-rwxr-xr-x551/CH14/EX14.20/20.sce37
-rwxr-xr-x551/CH14/EX14.21/21.sce17
-rwxr-xr-x551/CH14/EX14.22/22.sce26
-rwxr-xr-x551/CH14/EX14.23/23.sce21
-rwxr-xr-x551/CH14/EX14.24/24.sce24
-rwxr-xr-x551/CH14/EX14.25/25.sce60
-rwxr-xr-x551/CH14/EX14.26/26.sce26
-rwxr-xr-x551/CH14/EX14.27/27.sce49
-rwxr-xr-x551/CH14/EX14.3/3.sce11
-rwxr-xr-x551/CH14/EX14.4/4.sce13
-rwxr-xr-x551/CH14/EX14.5/5.sce16
-rwxr-xr-x551/CH14/EX14.6/6.sce26
-rwxr-xr-x551/CH14/EX14.7/7.sce33
-rwxr-xr-x551/CH14/EX14.8/8.sce14
-rwxr-xr-x551/CH14/EX14.9/9.sce21
27 files changed, 874 insertions, 0 deletions
diff --git a/551/CH14/EX14.1/1.sce b/551/CH14/EX14.1/1.sce
new file mode 100755
index 000000000..418a1d84f
--- /dev/null
+++ b/551/CH14/EX14.1/1.sce
@@ -0,0 +1,27 @@
+clc
+T2=235; //K
+P=1.3; //kW
+
+disp("(i) C.O.P. of Carnot refrigerator =")
+COP=14000/P/60/60;
+disp(COP)
+
+
+disp("(ii) Higher temperature of the cycle =")
+T1=T2/COP + T2;
+t1=T1-273;
+disp(t1)
+disp("0C")
+
+
+disp("(iii) Heat delivered as heat pump")
+Qabs=14000/60; //Heat absorbed
+W=P*60;
+Q=Qabs+W;
+disp("Q=")
+disp(Q)
+disp("kJ/min")
+
+COP=Q/W;
+disp("COP of heat pump =")
+disp(COP) \ No newline at end of file
diff --git a/551/CH14/EX14.10/10.sce b/551/CH14/EX14.10/10.sce
new file mode 100755
index 000000000..cd9106b66
--- /dev/null
+++ b/551/CH14/EX14.10/10.sce
@@ -0,0 +1,48 @@
+clc
+p1=1000; //kPa
+p2=100; //kPa
+p4=p1;
+p3=p2;
+E=2000; // Refrigerating effect produced in kJ/min
+T3=268; //K
+T1=303; //K
+y=1.4;
+
+disp("(i) Mass of air circulated per minute")
+T2=T1*(p2/p1)^((y-1)/y);
+e=cp*(T3-T2); //Refrigerating effect per kg; kJ/kg
+
+m=E/e;
+disp("m=")
+disp(m)
+disp("kg/min")
+
+
+disp("(ii) Compressor work (Wcomp.), expander work (Wexp.) and cycle work (Wcycle)")
+T4=T3*(p4/p3)^((y-1)/y);
+
+Wcomp=y/(y-1)*m*R*(T4-T3);
+disp("Compressor work =")
+disp(Wcomp)
+disp("kJ/min")
+
+Wexp=y/(y-1)*m*R*(T1-T2);
+disp("Expander work =")
+disp(Wexp)
+disp("kJ/min")
+
+W_cycle=Wcomp-Wexp;
+disp("Wcycle=")
+disp(W_cycle)
+disp("kJ/min")
+
+
+disp("(iii) C.O.P. and power required")
+COP=E/W_cycle;
+disp("COP =")
+disp(COP)
+
+P=W_cycle/60;
+disp("Power required =")
+disp(P)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.11/11.sce b/551/CH14/EX14.11/11.sce
new file mode 100755
index 000000000..6b75dc31d
--- /dev/null
+++ b/551/CH14/EX14.11/11.sce
@@ -0,0 +1,55 @@
+clc
+y=1.4;
+cp=1.003; //kJ/kg K
+T3=289; //K
+T1=314; //K
+p1=5.2; //bar
+p2=1; //bar
+capacity=6; //tonnes
+R=287; //J/kg K
+l=0.2; //m
+
+T4=T3*(p1/p2)^((y-1)/y);
+T2=T1*(p2/p1)^((y-1)/y);
+
+
+disp("(i) C.O.P. =")
+COP=T2/(T1-T2);
+disp(COP)
+
+disp("(ii) Mass of air in circulation")
+e=cp*(T3-T2); //Refrigerating effect per kg of air
+E=capacity*14000; //Refrigerating effect produced by the refrigerating machine in kJ/h
+
+m=E/e/60;
+disp("mass of air in circulation =")
+disp(m)
+disp("kg/min")
+
+
+disp("Piston displacement of compressor")
+V3=m*R*T3/p2/10^5;
+
+V_swept=V3/2/240;
+
+d_c=sqrt(V_swept/l/%pi*4);
+
+disp("Diameter or bore of the compressor cylinder =")
+disp(d_c*1000)
+disp("mm")
+
+disp("Piston displacement of expander")
+V2=m*R*T2/p2/10^5;
+V_swept=V2/2/240;
+
+d_c=sqrt(V_swept/l/%pi*4);
+disp("Diameter or bore of the expander cylinder =")
+disp(d_c*1000)
+disp("mm")
+
+
+disp("(v) Power required to drive the unit")
+W=capacity*14000/COP/3600;
+disp("power =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.12/12.sce b/551/CH14/EX14.12/12.sce
new file mode 100755
index 000000000..af02dde34
--- /dev/null
+++ b/551/CH14/EX14.12/12.sce
@@ -0,0 +1,30 @@
+clc
+m=6; //kg/min
+n_relative=0.50;
+cpw=4.187; //kJ/kg K
+L=335; //kJ/kg
+
+h_f2=31.4; //kJ/kg
+h_fg2=154; //kJ/kg
+h_f3=59.7; //kJ/kg
+h_fg3=138; //kJ/kg
+h_f4=59.7; //kJ/kg
+x2=0.6;
+s_f3=0.2232; //kJ/kg K
+s_f2=0.1251; //kJ/kg K
+T2=268; //K
+T3=298; //K
+
+h2=h_f2+x2*h_fg2;
+x3=((s_f2-s_f3)+x2*(h_fg2/T2))*T3/h_fg3;
+h3=h_f3+x3*h_fg3;
+h1=h_f4;
+COP_th=(h2-h1)/(h3-h2); //Theoritical COP
+COP=n_relative*COP_th;
+
+Q=cpw*(20-0) + L; //Heat extracted from 1 kg of water at 20°C for the formation of 1 kg of ice at 0°C
+
+m_ice=COP*m*(h3-h2)/Q*60*24/1000; //in 24 hours
+disp("m_ice=")
+disp(m_ice)
+disp("tonnes") \ No newline at end of file
diff --git a/551/CH14/EX14.13/13.sce b/551/CH14/EX14.13/13.sce
new file mode 100755
index 000000000..4449b3862
--- /dev/null
+++ b/551/CH14/EX14.13/13.sce
@@ -0,0 +1,25 @@
+clc
+L=335; //kJ/kg
+h3=1319.22; //kJ/kg
+h1=100.04; //kJ/kg
+h4=h1;
+s_f2=-2.1338; //kJ/kg K
+s_g2=5.0585; //kJ/kg K
+s_g3=4.4852; //kJ/kg K
+h_f2=-54.56; //kJ/kg
+h_g2=1304.99; //kJ/kg
+
+x2=(s_g3-s_f2)/(s_g2-s_f2);
+
+h2=h_f2+x2*(h_g2-h_f2);
+COP_theoritical=(h2-h1)/(h3-h2);
+COP_actual=0.62*COP_theoritical;
+RE=COP_actual*(h3-h2); //Actual refrigerating effect per kg
+Q=28*1000*L/24/3600; //Heat to be extracted per second
+
+m=Q/RE; //Mass of refrigerant circulated per second
+
+W=m*(h3-h2);
+disp("Power required =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.14/14.sce b/551/CH14/EX14.14/14.sce
new file mode 100755
index 000000000..b373abaf6
--- /dev/null
+++ b/551/CH14/EX14.14/14.sce
@@ -0,0 +1,31 @@
+clc
+h_f2=158.2; //kJ/kg
+x2=0.62;
+h_fg2=1280.8;
+h1=298.9; //kJ/kg
+h_f4=h1;
+s_f2=0.630; //kJ/kg K
+T2=268; //K
+T3=298; //K
+s_f3=1.124; //kJ/kg K
+h_fg3=1167.1; //kJ/kg
+m=6.4; //kg/min
+cp=4.187;
+L=335; //kJ/kg
+h_f3=298.9; //kJ/kg
+
+h2=h_f2+x2*h_fg2;
+x3=((s_f2-s_f3)+x2*h_fg2/T2)/h_fg3*T3;
+h3=h_f3+x3*h_fg3;
+
+COP_theoritical=(h2-h1)/(h3-h2);
+COP_actual=0.55*COP_theoritical;
+
+W1=h3-h2; //Work done per kg of refrigerant
+W=m*W1/60; //Work done per second kJ/s
+
+Q=15*cp+L;
+m_ice=W*3600*24/Q;
+disp("Amount of ice formed in 24 hours =")
+disp(m_ice)
+disp("kg") \ No newline at end of file
diff --git a/551/CH14/EX14.15/15.sce b/551/CH14/EX14.15/15.sce
new file mode 100755
index 000000000..a7fb1da78
--- /dev/null
+++ b/551/CH14/EX14.15/15.sce
@@ -0,0 +1,30 @@
+clc
+RE=5*14000/3600; //Total refrigeration produced in kg/s
+h2=183.19; //kJ/kg
+h3=209.41; //kJ/kg
+h4=74.59; //kJ/kg
+h1=h4;
+
+disp("(i) The refrigerant flow rate")
+RE_net=h2-h1; //Net refrigerating effect produced per kg
+m=RE/RE_net;
+disp("Refrigerant flow rate =")
+disp(m)
+disp("kg/s")
+
+
+disp("(ii) The C.O.P. =")
+COP=(h2-h1)/(h3-h2);
+disp(COP)
+
+
+disp("(iii) The power required to drive the compressor =")
+P=m*(h3-h2);
+disp(P)
+disp("kW")
+
+
+disp("(iv) The rate of heat rejection to the condenser =")
+rate=m*(h3-h4);
+disp(rate)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.16/16.sce b/551/CH14/EX14.16/16.sce
new file mode 100755
index 000000000..167ff6992
--- /dev/null
+++ b/551/CH14/EX14.16/16.sce
@@ -0,0 +1,129 @@
+clc
+
+disp("(i) If an expansion cylinder is used in a vapour compression system, the work recovered would be extremely small, in fact not even sufficient to overcome the mechanical friction. It will not be possible to gain any work. Further, the expansion cylinder is bulky. On the other hand the expansion valve is a very simple and handy device, much cheaper than the expansion cylinder. It does not need installation, lubrication or maintenance.")
+disp("The expansion valve also controls the refrigerant flow rate according to the requirement, in addition to serving the function of reducting the pressure of the refrigerant.")
+
+
+disp("(ii) The comparison between centrifugal and reciprocating compressors ")
+
+disp("1.Suitability")
+
+disp("Centrifugal compressor")
+disp("Suitable for handling large volumes of air at low pressures")
+
+disp("Reciprocating compressor")
+disp("Suitable for low discharges of air at high pressure.")
+
+
+disp("2.Operational speeds")
+
+disp("Centrifugal compressor")
+disp("Usually high")
+
+disp("Reciprocating compressor")
+disp("Low")
+
+
+disp("3.Air supply")
+
+disp("Centrifugal compressor")
+disp("Continuous")
+
+disp("Reciprocating compressor")
+disp("Pulsating")
+
+
+disp("4.Balancing")
+
+disp("Centrifugal compressor")
+disp("Less Vibrations")
+
+disp("Reciprocating compressor")
+disp("Cyclic vibrations occur")
+
+
+disp("5.Lubrication system")
+
+disp("Centrifugal compressor")
+disp("Generally simple lubrication systems are required.")
+
+disp("Reciprocating compressor")
+disp("Generally complicated")
+
+
+disp("6.Quality of air delivered")
+
+disp("Centrifugal compressor")
+disp("Air delivered is relatively more clean")
+
+disp("Reciprocating compressor")
+disp("Generally contaminated with oil.")
+
+
+disp("7.Air compressor size")
+
+disp("Centrifugal compressor")
+disp("Small for given discharge")
+
+disp("Reciprocating compressor")
+disp("Large for same discharge")
+
+
+disp("8.Free air handled")
+
+disp("Centrifugal compressor")
+disp("2000-3000 m3/min")
+
+disp("Reciprocating compressor")
+disp("250-300 m3/min")
+
+
+disp("9.Delivery pressure")
+
+disp("Centrifugal compressor")
+disp("Normally below 10 bar")
+
+disp("Reciprocating compressor")
+disp("500 to 800 bar")
+
+
+disp("10.Usual standard of compression")
+
+disp("Centrifugal compressor")
+disp("Isentropic compression")
+
+disp("Reciprocating compressor")
+disp("Isothermal compression")
+
+
+disp("11.Action of compressor")
+
+disp("Centrifugal compressor")
+disp("Dynamic action")
+
+disp("Reciprocating compressor")
+disp("Positive displacement")
+
+
+disp("(iii)")
+h2=344.927; //kJ/kg
+h4=228.538; //kJ/kg
+h1=h4;
+cpv=0.611; ///kJ/kg0C
+// s2=s3
+t3=39.995; //0C
+h3=363.575+cpv*(t3-30);
+Rn=h2-h1;
+W=h3-h2;
+
+COP=Rn/W;
+disp("COP =")
+disp(COP)
+
+cp=2.0935; //kJ/kg 0C
+Q=2400/24/3600*[4.187*(15-0)+335+cp*(0-(-5))];
+
+W=Q/COP;
+disp("Work required =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.17/17.sce b/551/CH14/EX14.17/17.sce
new file mode 100755
index 000000000..f69510832
--- /dev/null
+++ b/551/CH14/EX14.17/17.sce
@@ -0,0 +1,28 @@
+clc
+disp("(ii) Mass of refrigerant circulated per minute")
+h2=352; //kJ/kg
+h3=374; //kJ/kg
+h4=221; //kJ/kg
+h1=h4;
+v2=0.08; //m^3/kg
+rpm=500;
+D=0.2;
+L=0.15;
+n_vol=0.85;
+
+RE=h2-h1;
+V=%pi/4*D^2*L*rpm*2*n_vol;
+
+m=V/v2;
+disp("Mass of refrigerant circulated per minute = ")
+disp(m)
+disp("kg/min")
+
+disp("(iii) Cooling capacity in tonnes of refrigeration =")
+cc=50*(h2-h1)*60/14000;
+disp(cc)
+disp("TR")
+
+disp("(iv)COP =")
+COP=(h2-h1)/(h3-h2);
+disp(COP) \ No newline at end of file
diff --git a/551/CH14/EX14.18/18.sce b/551/CH14/EX14.18/18.sce
new file mode 100755
index 000000000..750fab604
--- /dev/null
+++ b/551/CH14/EX14.18/18.sce
@@ -0,0 +1,22 @@
+clc
+te=-10; //0C
+tc=40; //0C
+h3=220; //kJ/kg
+h2=183.1; //kJ/kg
+h1=74.53; //kJ/kg
+h_f4=26.85; //kJ/kg
+m=1; //kg
+
+disp("(i) The C.O.P. the cycle =")
+COP=(h2-h1)/(h3-h2);
+disp(COP)
+
+disp("(ii) Refrigerating capacity =")
+RC=m*(h2-h1);
+disp(RC)
+disp("kJ/min")
+
+disp("Compressor power =")
+CP=m*(h3-h2)/60;
+disp(CP)
+disp("kJ/s") \ No newline at end of file
diff --git a/551/CH14/EX14.19/19.sce b/551/CH14/EX14.19/19.sce
new file mode 100755
index 000000000..93899bc95
--- /dev/null
+++ b/551/CH14/EX14.19/19.sce
@@ -0,0 +1,33 @@
+clc
+h2=178.61; //kJ/kg
+h3a=203.05; //kJ/kg
+h_f4=74.53; //kJ/kg
+h1=h_f4;
+s3a=0.682; //kJ/kg K
+s2=0.7082; //kJ/kg K
+cp=0.747; //kJ/kg K
+T3a=313; //K
+CE=20; //Cooling effect
+C=0.03;
+v_g=0.1088;
+p_d=9.607;
+p_s=1.509;
+n=1.13;
+
+m=CE/(h2-h1);
+T3=T3a*%e^((s2-s3a)/cp)
+h3=h3a+cp*(T3-T3a);
+
+P=m*(h3-h2);
+disp("Power required by the machine =")
+disp(P)
+disp("kW")
+
+n_vol=1+C-C*(p_d/p_s)^(1/n); //Volumetric efficiency
+V1=m*v_g; //volume of refrigerant at the intake conditions
+V_swept=V1/n_vol;
+
+V=V_swept*60/300;
+disp("Piston displacement =")
+disp(V)
+disp("m^3") \ No newline at end of file
diff --git a/551/CH14/EX14.2/2.sce b/551/CH14/EX14.2/2.sce
new file mode 100755
index 000000000..a6017eb99
--- /dev/null
+++ b/551/CH14/EX14.2/2.sce
@@ -0,0 +1,22 @@
+clc
+T1=308; //K
+T2=258; //K
+capacity=12; //tonne
+
+COP=T2/(T1-T2);
+disp("(i) Co-efficient of performance =")
+disp(COP)
+
+
+disp("(ii) Heat rejected from the system per hour")
+W=capacity*14000/5.16;
+Q=capacity*14000+W;
+disp("Q=")
+disp(Q)
+disp("kJ/h")
+
+
+disp("(iii) Power required =")
+P=W/60/60;
+disp(P)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.20/20.sce b/551/CH14/EX14.20/20.sce
new file mode 100755
index 000000000..091451387
--- /dev/null
+++ b/551/CH14/EX14.20/20.sce
@@ -0,0 +1,37 @@
+clc
+h2=1450.22; //kJ/kg
+h3a=1488.57; //kJ/kg
+h_f4=366.072; //kJ/kg
+cpl2=4.556; //kJ/kg K
+cpv1=2.492; //kJ/kg K
+cpv2=2.903; //kJ/kg K
+T1=303; //K
+T2=308; //K
+s3a=5.2086; //kJ/kg K
+s2=5.755; //kJ/kg K
+T3a=308; //K
+N=1000;
+
+h_f4a=h_f4-cpl2*(T2-T1);
+h1=h_f4a;
+T3=T3a*%e^((s2-s3a)/cpv2);
+h3=h3a+cpv2*(T3-T3a);
+m=50/(h2-h1);
+
+
+disp("(i) Power required =")
+P=m*(h3-h2);
+disp(P)
+disp("kW")
+
+
+disp("(ii) Cylinder dimensions ")
+D=(m*4*60/%pi/1.2/N/0.417477)^(1/3);
+disp("Diameter of cylinder =")
+disp(D)
+disp("m")
+
+L=1.2*D;
+disp("Length of the cylinder=")
+disp(L)
+disp("m") \ No newline at end of file
diff --git a/551/CH14/EX14.21/21.sce b/551/CH14/EX14.21/21.sce
new file mode 100755
index 000000000..d38900f4c
--- /dev/null
+++ b/551/CH14/EX14.21/21.sce
@@ -0,0 +1,17 @@
+clc
+cooling_load=150; //W
+n_vol=0.8;
+N=720; //rpm
+h2=183; //kJ/kg
+h1=74.5; //kJ/kg
+v2=0.08; //m^3/kg
+
+m=cooling_load/(108.5*1000);
+disp("Mass flow rate of the refrigerant =")
+disp(m)
+disp("kJ/s")
+
+d=m*v2/n_vol;
+disp("Displacement volume of the compressor =")
+disp(d)
+disp("m^3/s") \ No newline at end of file
diff --git a/551/CH14/EX14.22/22.sce b/551/CH14/EX14.22/22.sce
new file mode 100755
index 000000000..66b87d752
--- /dev/null
+++ b/551/CH14/EX14.22/22.sce
@@ -0,0 +1,26 @@
+clc
+h2=183.2; //kJ/kg
+h3=222.6; //kJ/kg
+h4=84.9; //kJ/kg
+
+v2=0.0767; //m^3/kg
+v3=0.0164; //m^3/kg
+v4=0.00083; //m^3/kg
+
+V=1.5*1000*10^(-6); //Piston displacement volume m^3/revolution
+n_vol=0.80;
+
+disp("(i) Power rating of the compressor (kW)")
+discharge=V*1600*n_vol; //Compressor discharge
+m=discharge/v2;
+
+P=m/60*(h3-h2); //kW
+disp("Power =")
+disp(P)
+disp("kW")
+
+
+disp("(ii) Refrigerating effect =")
+RE=m/60*(h2-h4);
+disp(RE)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.23/23.sce b/551/CH14/EX14.23/23.sce
new file mode 100755
index 000000000..d139ea66c
--- /dev/null
+++ b/551/CH14/EX14.23/23.sce
@@ -0,0 +1,21 @@
+clc
+COP=6.5;
+W=50; //kW
+h3a=201.45; //kJ/kg
+h_f4=69.55; //kJ/kg
+h1=h_f4;
+h2=187.53; //kJ/kg
+cp=0.6155; //kJ/kg
+t3a=35; //0C
+
+RC=W*COP; //Refrigerating capacity
+Q1=h2-h_f4; //Heat extracted per kg of refrigerant
+rate=RC/Q1; //Refrigerant flow rate
+Q2=W/rate; //Heat input per kg
+h=h2+Q2; //Enthalpy of vapour after compression
+Q=h-h3a; //Superheat
+
+t3=Q/cp+t3a;
+disp("t3=")
+disp(t3)
+disp("°C") \ No newline at end of file
diff --git a/551/CH14/EX14.24/24.sce b/551/CH14/EX14.24/24.sce
new file mode 100755
index 000000000..e3ce463f9
--- /dev/null
+++ b/551/CH14/EX14.24/24.sce
@@ -0,0 +1,24 @@
+clc
+Q1=500; //total heating requirement of 500 kJ/min
+n_compressor=0.8;
+s1=0.7035; //kJ/kg K
+s2=0.6799; //kJ/kg K
+T2=322.31; //K
+cp=0.7; //kJ/kg K
+h_v2=206.24; //kJ/kg
+h_l2=84.21; //kJ/kg
+h_v1=182.07 //kJ/kg
+
+Q2=Q1/n_compressor; //Heat rejected by the cycle
+
+//Entropy of dry saturated vapour at 2 bar= Entropy of superheated vapour at 12 bar
+T=T2*%e^((s1-s2)/cp);
+
+H=h_v2+cp*(T-T2); //Enthalpy of superheated vapour at 12 bar
+Q3=H-h_l2; //Heat rejected per cycle
+m=Q2/Q3; //kg/min
+W=m*(H-h_v1)/60; //kW
+W_actual=W/n_compressor;
+disp("Power =")
+disp(W_actual)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.25/25.sce b/551/CH14/EX14.25/25.sce
new file mode 100755
index 000000000..3ca6233aa
--- /dev/null
+++ b/551/CH14/EX14.25/25.sce
@@ -0,0 +1,60 @@
+clc
+h2a=183.2; //kJ/kg K
+cpv=0.733; //Vapour specific heat in kJ/kg K
+cpl=1.235; //Liquid specific heat in kJ/kg K
+s2a=0.7020; //Entropy of vapour in kJ/kg K
+s3a=0.6854; //Entropy of vapour in kJ/kg K
+T2=270; //K
+T2a=263; //K
+T3a=303; //K
+h3a=199.6; //kJ/kg
+h_f4=64.6; //kJ/kg
+dT4=6; //dT4=T4-T4a
+v2a=0.0767;
+n=2; //number of cylinder
+
+h2=h2a+cpv*(T2-T2a);
+s2=s2a+cpv*log(T2/T2a);
+
+T3=T3a*%e^((s2-s3a)/cpv);
+h3=h3a+cpv*(T3-T3a);
+h_f4a=h_f4-cpl*dT4;
+h1=h_f4a;
+
+v2=v2a/T2a*T2;
+
+disp("(i) Refrigerating effect per kg =")
+RE=h2-h1;
+disp(RE)
+disp("kJ/kg")
+
+disp("(ii) Mass of refrigerant to be circulated per minute =")
+m=2400/RE;
+disp(m)
+disp("kg/min")
+
+disp("(iii) Theoretical piston displacement per minute =")
+v=m*v2;
+disp(v)
+disp("m^3/min")
+
+disp("(iv) Theoretical power required to run the compressor = ")
+P=m/60*(h3-h2);
+disp(P)
+disp("kW")
+
+disp("(v) Heat removed through the condenser per min =")
+Q=m*(h3-h_f4a);
+disp(Q)
+disp("kJ/min")
+
+disp("(vi) Theoretical bore (d) and stroke (l)")
+d=(v/n/%pi*4/1.25/1000)^(1/3)*1000;
+disp("Theroritical bore =")
+disp(d)
+disp("mm")
+
+disp("stroke =")
+l=1.25*d;
+disp(l)
+disp("mm") \ No newline at end of file
diff --git a/551/CH14/EX14.26/26.sce b/551/CH14/EX14.26/26.sce
new file mode 100755
index 000000000..3222524de
--- /dev/null
+++ b/551/CH14/EX14.26/26.sce
@@ -0,0 +1,26 @@
+clc
+h2=1597; //kJ/kg
+h3=1790; //kJ/kg
+h4=513; //kJ/kg
+h1=h4;
+t3=58; //0C
+x1=0.13;
+tc=27; //0C
+capacity=10.5; //tonnes
+
+disp("(i) Condition of the vapour at the outlet of the compressor =")
+t=t3-tc;
+disp(t)
+disp("°C")
+
+disp("(ii) Condition of vapour at entrance to evaporator =")
+disp(x1)
+
+disp("COP =")
+COP=(h2-h1)/(h3-h2);
+disp(COP)
+
+disp("(iv) Power required =")
+P=capacity*14000/COP/3600;
+disp(P)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.27/27.sce b/551/CH14/EX14.27/27.sce
new file mode 100755
index 000000000..c3db653fc
--- /dev/null
+++ b/551/CH14/EX14.27/27.sce
@@ -0,0 +1,49 @@
+clc
+h2=615; //kJ/kg
+h3=664; //kJ/kg
+h4=446; //kJ/kg
+h1=h4;
+v2=0.14; //m^3/kg
+capacity=20; //tonnes
+n=6; //number of cylinder
+
+disp("(i) Refrigerating effect per kg =")
+RE=h2-h1;
+disp(RE)
+disp("kJ/kg")
+
+disp("(ii) Mass of refrigerant to be circulated per minute =")
+m=capacity*14000/RE/60;
+disp(m)
+disp("kg/min")
+
+disp("(iii) Theoretical piston displacement =")
+v=v2*m;
+disp(v)
+disp("m^3/min")
+
+disp("(iv) Theoretical power =")
+P=m/60*(h3-h2);
+disp(P)
+disp("kW")
+
+disp("(v)COP =")
+COP=(h2-h1)/(h3-h2);
+disp(COP)
+
+disp("(vi) Heat removed through the condenser =")
+Q=m*(h3-h4);
+disp(Q)
+disp("kJ/min")
+
+disp("(vii) Theoretical displacement per minute per cylinder")
+
+d=(v/n*4/%pi/950)^(1/3)*1000;
+disp("Diameter of cylinder =")
+disp(d)
+disp("mm")
+
+l=d;
+disp("Stroke length =")
+disp(l)
+disp("mm") \ No newline at end of file
diff --git a/551/CH14/EX14.3/3.sce b/551/CH14/EX14.3/3.sce
new file mode 100755
index 000000000..243ef1939
--- /dev/null
+++ b/551/CH14/EX14.3/3.sce
@@ -0,0 +1,11 @@
+clc
+T2=268; //K
+T1=308; //K
+Q=29; //Heat leakage from the surroundings into the cold storage in kW
+COP_ideal=T2/(T1-T2);
+COP_actual=1/3*COP_ideal;
+
+W=Q/COP_actual;
+disp("Power required =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.4/4.sce b/551/CH14/EX14.4/4.sce
new file mode 100755
index 000000000..ae0563c4d
--- /dev/null
+++ b/551/CH14/EX14.4/4.sce
@@ -0,0 +1,13 @@
+clc
+T1=293; //K
+T2=265; //K
+T0=273; //K
+L=335; //Latent heat of ice in kJ/kg
+cpw=4.18;
+
+COP=T2/(T1-T2);
+Rn=cpw*(T1-T0)+L;
+m_ice=COP*3600/Rn;
+disp("ice formed per kWh =")
+disp(m_ice)
+disp("kg") \ No newline at end of file
diff --git a/551/CH14/EX14.5/5.sce b/551/CH14/EX14.5/5.sce
new file mode 100755
index 000000000..3d88afec2
--- /dev/null
+++ b/551/CH14/EX14.5/5.sce
@@ -0,0 +1,16 @@
+clc
+T1=291; //K
+T2=265; //K
+T0=273; //K
+cpw=4.18; //kJ/kg
+cpi=2.09; //kJ/kg
+L=334; //kJ/kg
+m=400; //kg
+
+COP=T2/(T1-T2);
+Rn=cpw*(T1-T0) + L + cpi*(T0-T2);
+
+W=Rn*m/COP/3600; //kJ/s
+disp("Least power =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.6/6.sce b/551/CH14/EX14.6/6.sce
new file mode 100755
index 000000000..d327a80a7
--- /dev/null
+++ b/551/CH14/EX14.6/6.sce
@@ -0,0 +1,26 @@
+clc
+cpw=4.18; //kJ/kg
+
+disp("(i) Quantity of ice produced")
+t=20; //0C
+L=335; //kJ/kg
+capacity=280; //tonnes
+
+Q1=cpw*t + L; //Heat to be extracted per kg of water (to form ice at 0°C)
+Rn=capacity*14000; //kJ/h
+
+m_ice=Rn*24/Q1/1000;
+disp("Quantity of ice produced in 24 hours =")
+disp(m_ice)
+disp("tonnes")
+
+
+disp("(ii) Minimum power required =")
+T1=298; //K
+T2=263; //K
+
+COP=T2/(T1-T2);
+W=Rn/COP/3600; //kJ/s
+disp("Power required =")
+disp(W)
+disp("kW") \ No newline at end of file
diff --git a/551/CH14/EX14.7/7.sce b/551/CH14/EX14.7/7.sce
new file mode 100755
index 000000000..cb7ab1852
--- /dev/null
+++ b/551/CH14/EX14.7/7.sce
@@ -0,0 +1,33 @@
+clc
+cp1=1.25; //kJ/kg 0C
+cp2=2.93; //kJ/kg 0C
+L=232; //kJ/kg
+T1=-3; //0C
+T2=-8; //0C
+T3=25; //0C
+
+Q1=cp2*(T3-T1) + L + cp1*(T1-T2); //Heat removed in 8 hours from each kg of fish
+
+Q=Q1*20*1000/8; //Heat removed by the plant /min
+
+disp("(i) Capacity of the refrigerating plant =")
+capacity=Q/14000; //tonnes
+disp(capacity)
+disp("tonnes")
+
+disp("(ii) Carnot cycle C.O.P. between this temperature range.")
+T1=298; //K
+T2=265; //K
+
+COP=T2/(T1-T2);
+disp("COP of reversed carnot cycle =")
+disp(COP)
+
+
+disp("(iii) Power required")
+COP_actual=1/3*COP;
+
+W=Q/COP_actual/3600; //kJ/s
+disp("Power =")
+disp(W)
+disp("kW")
diff --git a/551/CH14/EX14.8/8.sce b/551/CH14/EX14.8/8.sce
new file mode 100755
index 000000000..df63ab405
--- /dev/null
+++ b/551/CH14/EX14.8/8.sce
@@ -0,0 +1,14 @@
+clc
+T1=1273; //K
+T2=298; //K
+T3=268; //K
+T4=298; //K
+
+//Let Q2/Q1=r1, r2=Q3/Q4;
+r1=298/1273; //Q2/Q1
+r2=268/298; //Q3/Q4
+
+//Let Q4/Q1=r
+r=(1-r1)/(1-r2);
+disp("The ratio in which the heat pump and heat engine share the heating load =")
+disp(r)
diff --git a/551/CH14/EX14.9/9.sce b/551/CH14/EX14.9/9.sce
new file mode 100755
index 000000000..e76b69765
--- /dev/null
+++ b/551/CH14/EX14.9/9.sce
@@ -0,0 +1,21 @@
+clc
+y=1.4;
+n=1.35;
+cp=1.003; //kJ/kg K
+p2=1; //bar
+p1=8; //bar
+T3=282; //K
+T4=302; //K
+T1=T4;
+T4=T3*(p1/p2)^((n-1)/n);
+T2=T1*(p2/p1)^((n-1)/n);
+
+Q1=cp*(T3-T2); //Heat extracted from cold chamber per kg of air
+Q2=cp*(T4-T1); //Heat rejected in the cooling chamber per kg of air
+cv=cp/y;
+R=cp-cv;
+W=n/(n-1)*R*((T4-T3) - (T1-T2));
+
+COP=Q1/W;
+disp("COP=")
+disp(COP) \ No newline at end of file