diff options
Diffstat (limited to '2144/CH6')
34 files changed, 639 insertions, 0 deletions
diff --git a/2144/CH6/EX6.1/ex6_1.sce b/2144/CH6/EX6.1/ex6_1.sce new file mode 100755 index 000000000..f2aebdf01 --- /dev/null +++ b/2144/CH6/EX6.1/ex6_1.sce @@ -0,0 +1,11 @@ +// Exa 6.1
+clc;
+clear;
+close;
+// Given data
+P = 10;// in bar
+P = P *10^5;// in N/m^2
+V = 2;//volume of water in m^3
+W = P * V;// in J
+W = W * 10^-6;// in MJ
+disp(W,"Work done in MJ is");
diff --git a/2144/CH6/EX6.10/ex6_10.sce b/2144/CH6/EX6.10/ex6_10.sce new file mode 100755 index 000000000..0031c2e68 --- /dev/null +++ b/2144/CH6/EX6.10/ex6_10.sce @@ -0,0 +1,20 @@ +// Exa 6.10
+clc;
+clear;
+close;
+// Given data
+T = 0;// in degree C
+T = T + 273;// in K
+T_sat = 179.9;// in degree C
+T_sat = T_sat + 273;// in K
+x = 0.8;
+h_fg = 2013.8;// in kJ/kg
+c_f = 4.188;
+Phi_wet = c_f*log(T_sat/T)+x*h_fg/T_sat;// in kJ/kg-K
+disp(Phi_wet,"The entropy of wet steam in kJ/kg-K is");
+Phi_g = (c_f*(log(T_sat/T))) + (h_fg/T_sat);// in kJ/kg-K
+disp(Phi_g,"The entropy of dry saturated steam in kJ/kg-K is");
+C_P = 2.3;
+T_sup = 200+273;// in K
+Phi = c_f *log(T_sat/T) + h_fg/T_sat+ C_P*log(T_sup/T_sat);// in kJ/kg-K
+disp(Phi,"The entropy of superheated steam in kJ/kg-K is");
diff --git a/2144/CH6/EX6.11/ex6_11.sce b/2144/CH6/EX6.11/ex6_11.sce new file mode 100755 index 000000000..64e39c143 --- /dev/null +++ b/2144/CH6/EX6.11/ex6_11.sce @@ -0,0 +1,15 @@ +// Exa 6.11
+clc;
+clear;
+close;
+// Given data
+T_sat = 198.3;// in degree C
+T_sat = T_sat + 273;// in K
+T_sup = 300;// in degree C
+T_sup = T_sup + 273;// in K
+c_f = 4.188;
+h_fg = 1945;// in kJ/kg-K
+T = 273;// in K
+C_P = 2.093;// in kJ/kg-K
+Phi_sup =c_f *log(T_sat/T)+h_fg/T_sat+C_P*log(T_sup/T_sat);// in kJ/kg-K
+disp(Phi_sup,"The value of specific entropy in kJ/kg-K is");
diff --git a/2144/CH6/EX6.12/ex6_12.sce b/2144/CH6/EX6.12/ex6_12.sce new file mode 100755 index 000000000..b2f45a282 --- /dev/null +++ b/2144/CH6/EX6.12/ex6_12.sce @@ -0,0 +1,10 @@ +// Exa 6.12
+clc;
+clear;
+close;
+//Given data
+P = 16;// in bar
+m_w = 73;// in gm
+m_s = 980;// in gm
+x = m_s/(m_s+m_w);
+disp(x,"Dryness fraction of steam is");
diff --git a/2144/CH6/EX6.13/ex6_13.sce b/2144/CH6/EX6.13/ex6_13.sce new file mode 100755 index 000000000..54e0db878 --- /dev/null +++ b/2144/CH6/EX6.13/ex6_13.sce @@ -0,0 +1,15 @@ +// Exa 6.13
+clc;
+clear;
+close;
+// Given data
+P1 = 7;// in bar
+P2 = 1.2;// in bar
+h_f1 = 696.9;// in kJ/kg
+h_fg1 = 2065;// in kJ/kg
+h_g2 = 2684.9;// in kJ/kg
+T_sup = 112;// in degree C
+T_sat = 104.77;// in degree C
+C_P = 2.1;// in kJ/kg
+x1 = (h_g2+(C_P*(T_sup-T_sat))-h_f1)/h_fg1;
+disp(x1,"Dryness fraction of steam is");
diff --git a/2144/CH6/EX6.14/ex6_14.sce b/2144/CH6/EX6.14/ex6_14.sce new file mode 100755 index 000000000..1f717a650 --- /dev/null +++ b/2144/CH6/EX6.14/ex6_14.sce @@ -0,0 +1,12 @@ +// Exa 6.14
+clc;
+clear;
+close;
+// Given data
+P1 = 14;// in bar
+P2 = 1.2;// in bar
+h_f1 = 830;
+h_fg1 = 1958;
+h_g2 = 2684.9;
+x = (h_g2-h_f1)/h_fg1;
+disp(x,"Dryness fraction of steam is");
diff --git a/2144/CH6/EX6.15/ex6_15.sce b/2144/CH6/EX6.15/ex6_15.sce new file mode 100755 index 000000000..29128955c --- /dev/null +++ b/2144/CH6/EX6.15/ex6_15.sce @@ -0,0 +1,17 @@ +// Exa 6.15
+clc;
+clear;
+close;
+// Given data
+m_s = 2.2;// in kg
+m_w = 0.18;// in kg
+x1 = m_s/(m_s+m_w);
+h_f1 = 743;
+h_fg1 = 2031;
+h_g2 = 2685;
+C_P = 2;
+T_sup = 115;// in degree C
+T_sat = 104.8;// in degree C
+x2 = (h_g2 + (C_P*(T_sup-T_sat)) - h_f1)/h_fg1;
+x = x1 * x2;
+disp(x,"The dryness fraction of steam is");
diff --git a/2144/CH6/EX6.16/ex6_16.sce b/2144/CH6/EX6.16/ex6_16.sce new file mode 100755 index 000000000..95e4f3644 --- /dev/null +++ b/2144/CH6/EX6.16/ex6_16.sce @@ -0,0 +1,17 @@ +// Exa 6.16
+clc;
+clear;
+close;
+// Given data
+h_f1 = 232;// in kJ/kg
+h_fg = 2369;// in kJ/kg
+x = 0.8;
+h_f2 = 167.5;// in kJ/kg
+H_wet1 = h_f1 + (x*h_fg);// in kJ/kg
+H_wet = H_wet1 - h_f2;// in kJ/kg
+T1 = 38;// in degree C
+T2 = 25;// in degree C
+T = T1-T2;// in degree C
+SpeHeat = 4.188;// in kJ/kg-K
+m = H_wet/(T*SpeHeat);// in kJ/kg
+disp(round(m),"The quantity of circulating water required of condensed steam in kJ/kg is");
diff --git a/2144/CH6/EX6.17/ex6_17.sce b/2144/CH6/EX6.17/ex6_17.sce new file mode 100755 index 000000000..2ce1a922b --- /dev/null +++ b/2144/CH6/EX6.17/ex6_17.sce @@ -0,0 +1,28 @@ +// Exa 6.17
+clc;
+clear;
+close;
+// Given data
+V1 = 0.4;// volume of dry saturated steam
+P1 = 1.5;// in MN/m^2
+disp("Part (i) : For Isothermal operation :")
+Vs = 0.1318;// specific volume of dry steam
+m = V1/Vs;// quantity of steam present in the vessel in kg
+h_f1= 844.6;// in kJ/kg
+x1= 0.5;// dryness fraction
+h_fg1= 1945.2;// in kJ/kg
+Specific_Enth= h_f1+x1*h_fg1;// in kJ/kg
+En= Specific_Enth*m;// kJ
+disp(En,"Enthalpy of the fluid in kJ is : ")
+HeatLost= m*(1-x1)*h_fg1;// in kJ
+disp(HeatLost,"The loss of heat during the constant temperature process in kJ is : ");
+disp("Part (ii) : For Hyperbolic operation :")
+h_f2= 1008.3;// in kJ/kg
+h_fg2= 1794;// in kJ/kg
+Vs= 0.0659;// Specific volume after compression in m^3/kg
+Vs1= 0.0666;// Specific volume of dry saturated steam in m^3/kg
+x2=Vs/Vs1;
+H_wet= h_f2+x2*h_fg2;// in kJ/kg
+H= m*H_wet;// in kJ
+disp(H,"Enthalpy of the fluid in kJ is :")
+
diff --git a/2144/CH6/EX6.18/ex6_18.sce b/2144/CH6/EX6.18/ex6_18.sce new file mode 100755 index 000000000..60f70350f --- /dev/null +++ b/2144/CH6/EX6.18/ex6_18.sce @@ -0,0 +1,17 @@ +// Exa 6.18
+clc;
+clear;
+close;
+// Given data
+P = 13.5;// power developed by engine in kW
+P1 = 12;// Steam consumption of the engine in kg/kWh
+S_C = P*P1;//steam consumed per hour in kg
+S_C = S_C/60;// in kg/min
+x = 0.85;
+V_g = 1.430;
+Volume = x * V_g;// in m^3/kg
+Volume = S_C * Volume;// in m^3
+d = 15*10^-2;//diameter of exhaust pipe in meter
+A = (%pi/4) * (d)^2;// in m^2
+C = Volume/A;// in meter/minute
+disp(C,"The velocity of steam in metre/minute is");
diff --git a/2144/CH6/EX6.19/ex6_19.sce b/2144/CH6/EX6.19/ex6_19.sce new file mode 100755 index 000000000..a1ac76495 --- /dev/null +++ b/2144/CH6/EX6.19/ex6_19.sce @@ -0,0 +1,11 @@ +// Exa 6.19
+clc;
+clear;
+close;
+// Given data
+P = 2;//pressure of steam in bar
+m = 0.1;//mass of steam in kg
+V = 0.080;//volume of steam in m^3
+V1 = 0.8872;//volume of 1kg dry saturated steam in m^3
+x = V/(m*V1);
+disp(x,"Dryness fraction of steam is");
diff --git a/2144/CH6/EX6.2/ex6_2.sce b/2144/CH6/EX6.2/ex6_2.sce new file mode 100755 index 000000000..1bd6b982b --- /dev/null +++ b/2144/CH6/EX6.2/ex6_2.sce @@ -0,0 +1,25 @@ +// Exa 6.2
+clc;
+clear;
+close;
+// Given data
+P = 1.013;// atm pressure in bar
+P = P * 10^5;// in N/m^2
+area= 1000*10^-4;// in m^2
+L_w = 1000;// in N
+P_L = L_w/area;// Pressure due to load in N/m^2
+PressOnPiston = P_L+P;// absolute pressure to piston in N/m^2
+a = 10^-3;// in m^2
+disp("Energy required to pump 1 kg of water at 0°C into the cylinder in joules is:")
+U = PressOnPiston*a;// in Joules
+disp(U);
+// Part (b)
+absPressure= 111.3*10^3;// in N/m^2
+increaseInVol= (1.02-1)*10^-3;// in m^3
+disp("Energy required to effect the change in volume in joules is : ")
+u_f= increaseInVol*absPressure;// in joules
+disp(u_f);
+// Part (c)
+increaseInVol= (1.52-0.001);// in m^3
+ExternalWorkDone= absPressure*increaseInVol;// in joules
+disp(ExternalWorkDone*10^-3,"External work done in kJ is :")
diff --git a/2144/CH6/EX6.20/ex6_20.sce b/2144/CH6/EX6.20/ex6_20.sce new file mode 100755 index 000000000..cd7af8bca --- /dev/null +++ b/2144/CH6/EX6.20/ex6_20.sce @@ -0,0 +1,15 @@ +// Exa 6.20
+clc;
+clear;
+close;
+// Given data
+P1 = 7+ 1;// in bar
+H = 2767;// Enthalpy in kJ/kg
+P2 = 1.5+1;// in bar
+H1 = 2717;// enthalpy of 1kg of dry steam in kJ/kg
+H_sup = H - H1;// Superheated of 1kg of steam in kJ
+S1 = 2.17;// super heated steam in kJ/kg-K
+theta = H_sup/S1;// in degree C
+T_sat = 127.4;// in degree C
+T_sup = T_sat + theta;// in degree C
+disp(T_sup,"The super heated temperature in degree C is");
diff --git a/2144/CH6/EX6.21/ex6_21.sce b/2144/CH6/EX6.21/ex6_21.sce new file mode 100755 index 000000000..d74e1fe74 --- /dev/null +++ b/2144/CH6/EX6.21/ex6_21.sce @@ -0,0 +1,15 @@ +// Exa 6.21
+clc;
+clear;
+close;
+// Given data
+T_sat = 99.6;// in degree C
+h_fg = 2258;// in kJ/kg
+m = 1;// steam output of the boiler in (assumed)
+m1 = 0.03;// exhaust steam
+x = 0.9;
+T1 = 21;// in degree C
+Cp = 4.187;// kJ/kg-K
+// Formula m1*(Cp*(T_sat-t)+x*h_fg)= m*Cp*(t-T1)
+t= (m1*(Cp*T_sat+x*h_fg)+m*Cp*T1)/(Cp*(m+m1))
+disp(t,"Temperature of the feed water leaving the heater in degree C is");
diff --git a/2144/CH6/EX6.22/ex6_22.sce b/2144/CH6/EX6.22/ex6_22.sce new file mode 100755 index 000000000..44b08cba5 --- /dev/null +++ b/2144/CH6/EX6.22/ex6_22.sce @@ -0,0 +1,12 @@ +// Exa 6.22
+clc;
+clear;
+close;
+// Given data
+T = 20;// in degree C
+H1 = 3039;// Enthalpy in kJ/kg
+H2 = 2725;// Enthalpy of 1kg dry saturated steam
+H_sup = H1-H2;// superheat of 1kg of steam in kJ/kg
+H= 2621.4;// heat required for 1kg or water in kJ
+m = H_sup/H;// in kg
+disp(m,"Quantity of water in kg is");
diff --git a/2144/CH6/EX6.23/ex6_23.sce b/2144/CH6/EX6.23/ex6_23.sce new file mode 100755 index 000000000..68b4baab7 --- /dev/null +++ b/2144/CH6/EX6.23/ex6_23.sce @@ -0,0 +1,12 @@ +// Exa 6.23
+clc;
+clear;
+close;
+// Given data
+x = 0.9;
+h_f1 = 1087.4;// in kJ/kg
+h_fg1 = 1712.9;// in kJ/kg
+H_wet1 = h_f1 + (x*h_fg1);// in kJ/kg
+H_sup2 = 3095;// in kJ/kg
+H = H_sup2 - H_wet1;// in kJ/kg
+disp(H,"Heat recieved in kJ/kg is");
diff --git a/2144/CH6/EX6.24/ex6_24.sce b/2144/CH6/EX6.24/ex6_24.sce new file mode 100755 index 000000000..25dcccfd1 --- /dev/null +++ b/2144/CH6/EX6.24/ex6_24.sce @@ -0,0 +1,16 @@ +// Exa 6.24
+clc;
+clear;
+close;
+// Given data
+V_fg =0.1632;// in m^3
+T_sup = 200;// in degree C
+T_sup = T_sup + 273;// in K
+T_sat = 188;// in degree C
+T_sat = T_sat + 273;// in K
+V_sup = (V_fg*T_sup)/T_sat;// in m^3/kg
+V = 0.24;// Capacity of the vessel in m^3
+Q = V/V_sup;// in kg
+V1 = 0.9774;//volume of 1kg dry saturated steam in m^3
+x = V_sup/V1;
+disp(x,"Dryness fraction of steam is");
diff --git a/2144/CH6/EX6.25/ex6_25.sce b/2144/CH6/EX6.25/ex6_25.sce new file mode 100755 index 000000000..0b13ce74a --- /dev/null +++ b/2144/CH6/EX6.25/ex6_25.sce @@ -0,0 +1,20 @@ +// Exa 6.25
+clc;
+clear;
+close;
+// Given data
+V = 0.6;// in m^3
+P2 = 2*10^2;// in kN/m^2
+P1 = 10*10^2;// in kN/m^2
+m = V/0.1946;// in kg
+V_s = 0.8872;// Specific volume of dry saturated steam in m^3
+x = 0.1946/V_s;
+h_f1 = 505;// in kJ/kg
+h_fg1 = 2202;// in kJ/kg
+H2 = m*(h_f1 + (x*h_fg1));// in kJ
+H1 = m*2776;// in kJ
+Q = (H2-H1) - (V*(P2-P1));// in kJ
+disp(m,"The mass of steam in the vessel in kg is : ")
+disp(x,"The dryness fraction of steam in the vessel is : ")
+disp(Q,"The amount o heat transferred in kJ is");
+disp("Thus during cooling process there is loss of heat")
diff --git a/2144/CH6/EX6.26/ex6_26.sce b/2144/CH6/EX6.26/ex6_26.sce new file mode 100755 index 000000000..eeceddc78 --- /dev/null +++ b/2144/CH6/EX6.26/ex6_26.sce @@ -0,0 +1,28 @@ +// Exa 6.26
+clc;
+clear;
+close;
+// Given data
+x1 = 0.95;
+P1 = 9;// in bar
+P1= P1*10^2;// in kN/m^2
+h_f1 = 743;// in kJ/kg
+h_fg1 = 2030;// in kJ/kg
+V = 0.204;// in m^3
+x2 = 0.544
+P2 = 5;// in bar
+P2= P2*10^2;// in kN/m^2
+h_f2 = 640;// in kJ/kg
+h_fg2 = 2108;// in kJ/kg
+H_wet1 = h_f1 + (x1*h_fg1);// in kJ/kg
+disp(H_wet1,"Total energy in kJ/kg is");
+U1 = H_wet1 - P1*V;// in kJ/kg
+disp(U1,"The internal energy in kJ/kg is : ")
+V_g1 = 0.204;// in m^3
+V1 = 0.3753;//volume of 1kg of dry stream in m^3
+x2 = V_g1/V1;
+H_wet2 = h_f2 + (x2*h_fg2);// in kJ
+U2 = H_wet2 - P2*V;// in kJ
+del_U = U1-U2;// in kJ
+H = del_U/V;// in kJ
+disp(H,"Heat removed from 1 m^3 of steam in kJ is");
diff --git a/2144/CH6/EX6.27/ex6_27.sce b/2144/CH6/EX6.27/ex6_27.sce new file mode 100755 index 000000000..56cd4baaa --- /dev/null +++ b/2144/CH6/EX6.27/ex6_27.sce @@ -0,0 +1,33 @@ +// Exa 6.27
+clc;
+clear;
+close;
+// Given data
+P1 = 2.1;// in MN/m^2
+P1= P1*10^3;//in kN/m^2
+P2 = 0.7;// in MN/m^2
+P2= P2*10^3;//in kN/m^2
+V1 = 0.1281;// in m^3
+x = 0.9;
+n = 1.25;
+h_f1= 920;// in kJ/kg
+h_fg1= 1878.6;// in kJ/kg
+h_f2= 697.0;// in kJ/kg
+h_fg2= 2065.0;// in kJ/kg
+V_wet1 = x * 0.0949;// in m^3/kg
+m = V1/V_wet1;// in kg
+disp(m,"Mass of steam in kg is");
+V2 = V1*((P1/P2)^(1/n));//in m^3
+del_W = (P1*V1-P2*V2)/(n-1);// in kJ
+disp(del_W,"External work done in kJ is");
+V_2 = V2/m;// in m^3/kg
+x2 = V_2/0.273;
+H1= h_f1+x*h_fg1;// in kJ/kg
+U1= H1-P1*V_wet1;// in kJ/kg
+H2= h_f2+x2*h_fg2;// in kJ/kg
+U2= H2-P2*V_2;// in kJ/kg
+del_E = m*(U2-U1);//in kJ
+disp(del_E,"Change in internal energy in kJ is");
+Q = del_W +del_E;// in kJ
+disp(Q,"Heat exchange in kJ is");
+disp("Heat is lost to the surroundings.")
diff --git a/2144/CH6/EX6.28/ex6_28.sce b/2144/CH6/EX6.28/ex6_28.sce new file mode 100755 index 000000000..f8b094aaf --- /dev/null +++ b/2144/CH6/EX6.28/ex6_28.sce @@ -0,0 +1,24 @@ +// Exa 6.28
+clc;
+clear;
+close;
+// Given data
+h_f1 = 670;// in kJ/kg
+h_fg1 = 2085;// in kJ/kg
+h_f2 = 475;// in kJ/kg
+h_fg2 = 2221;// in kJ/kg
+P2 = 6*10^2;// in kJ/kg
+P1 = 1.6*10^2;// in kJ/kg
+n = 1.1;
+x1 = 0.9;
+V1 = 0.3159;// in m^3
+V2 = 1.092;// in m^3
+H_wet = h_f1 + (x1*h_fg1);// in kJ/kg
+V_wet1 = x1*V1;// in m^3
+V_wet2 = V_wet1*(P2/P1)^(1/n);// in m^3
+x2 = V_wet2/V2;
+H_wet2 = h_f2 + (x2*h_fg2);// in kJ/kg
+U2= H_wet2-H_wet+P2*V_wet1-P1*V_wet2;// in kJ/kg
+W= (P2*V_wet1-P1*V_wet2)/(n-1);// in kJ/kg
+Q= U2+W;// in kJ/kg
+disp(Q,"Heat recieved by steam in kJ/kg is");
diff --git a/2144/CH6/EX6.29/ex6_29.sce b/2144/CH6/EX6.29/ex6_29.sce new file mode 100755 index 000000000..06d8215ab --- /dev/null +++ b/2144/CH6/EX6.29/ex6_29.sce @@ -0,0 +1,15 @@ +// Exa 6.29
+clc;
+clear;
+close;
+// Given data
+P1 = 0.85*10^3;// in kN/m^2
+P2 = 0.17*10^3;// in kN/m^2
+n = 1.13;
+x1 = 0.95;
+V1 = x1*0.227;// in m^3/kg
+V2 = V1 * ((P1/P2)^(1/n));// in m^3/kg
+x2 = V2/1.032;
+disp(x2,"Final dryness fraction of steam is");
+W = (P1*V1-P2*V2)/(n-1);// in kJ/kg
+disp(W,"Change in internal energy in kJ/kg is");
diff --git a/2144/CH6/EX6.3/ex6_3.sce b/2144/CH6/EX6.3/ex6_3.sce new file mode 100755 index 000000000..2fd3a9667 --- /dev/null +++ b/2144/CH6/EX6.3/ex6_3.sce @@ -0,0 +1,11 @@ +// Exa 6.3
+clc;
+clear;
+close;
+// Given data
+m_s = 92.3;// mass of steam in kg
+m_w = 0.78;// mass of water in kg
+m = m_s + m_w;// total mass in kg
+D_s = 92.3;// Dry steam in kg
+D_F = D_s/m;// Dryness fraction
+disp(D_F,"Dryness fraction is");
diff --git a/2144/CH6/EX6.30/ex6_30.sce b/2144/CH6/EX6.30/ex6_30.sce new file mode 100755 index 000000000..1dd856512 --- /dev/null +++ b/2144/CH6/EX6.30/ex6_30.sce @@ -0,0 +1,17 @@ +// Exa 6.30
+clc;
+clear;
+close;
+// Given data
+Cp= 2.3;// in kJ/kg-K
+T_sat= 179.9;// in °C
+T_sat= T_sat+273;// in K
+H= 3052;// enthalpy in kJ/kg
+P= 10*10^2;// in kN/m^2
+h_f= 763;// in kJ/kg
+h_fg= 2015;// in kJ/kg
+V= 0.1944;// in m^3
+// Formula H= h_f+h_fg*Cp*(t_sup-T_sat)-P*V*(t_sup/T_sat)
+t_sup= (h_f+h_fg-Cp*T_sat-H)/(P*V/T_sat-Cp);// in K
+t_sup= t_sup-273;// in °C
+disp(t_sup,"The final temperature of the steam in °C is : ")
diff --git a/2144/CH6/EX6.31/ex6_31.sce b/2144/CH6/EX6.31/ex6_31.sce new file mode 100755 index 000000000..c7039aa2d --- /dev/null +++ b/2144/CH6/EX6.31/ex6_31.sce @@ -0,0 +1,18 @@ +// Exa 6.31
+clc;
+clear;
+close;
+// Given data
+m1 = 3;// in kg
+m2 = 2;// in kg
+T1 = 10;// in degree C
+T2 = 80;// In Degree C
+T = ((m1*T1)+(m2*T2))/(m1+m2);// in degree C
+T = T + 273;// in K
+T1 = T1 + 273;// in K
+T2 = T2 + 273;// in K
+c_f = 4.188;
+del_phi1 = m1 * c_f*log(T/T1);// in kJ/K
+del_phi2 = m2 * c_f*log(T/T2);// in kJ/K
+Phi = del_phi1 + del_phi2;// in kJ/K
+disp(Phi,"Total change in entropy in kJ/K is")
diff --git a/2144/CH6/EX6.32/ex6_32.sce b/2144/CH6/EX6.32/ex6_32.sce new file mode 100755 index 000000000..a6db419d6 --- /dev/null +++ b/2144/CH6/EX6.32/ex6_32.sce @@ -0,0 +1,50 @@ +// Exa 6.32
+clc;
+clear;
+close;
+// Given data
+P1 = 15;// in bar
+P2 = 0.15;// in bar
+T_sat = 198.3;// in degree C
+T_sat = T_sat + 273;// in K
+h_fg1 = 1947;// in kJ/kg
+h_fg2= 2369;// in kJ/kg
+h_g1 = 845;// in kJ/kg
+h_f2 = 232;// in kJ/kg
+f_g2 = 7.985;// in kJ/kg-K
+x1 = 0.8;
+Phi_f1 = 2.315;// in kJ/kg-K
+Phi_f2 = 0.772;// in kJ/kg-K
+Phi1 = Phi_f1 + ((x1*h_fg1)/T_sat);// in kJ/kg-K
+H1 = h_g1 + (x1*h_fg1);// in kJ/kg-K
+Phi2 = Phi1;// in kJ/kg-K
+x2 = (Phi1 - Phi_f2)/(f_g2 - Phi_f2);
+H2 = h_f2 + (x2*h_fg2);// in kJ/kg
+Eta = ((H1-H2)/(H1-h_f2))*100;// in %
+disp("Part (i) When the steam supply is wet and dryness fraction is 0.8")
+disp(Eta,"Rankine efficiency in % is");
+delH = H1-H2;//theoretical work of steam in kJ/kg
+W = delH*60/100;// in kJ/kg
+Energy_Equivalent= 3600;// in kJ
+Steam_C = Energy_Equivalent/W;// Steam consumption in kg
+disp(Steam_C,"Steam consumption per kW-hr in kg is :")
+disp("Part (ii) When the steam supply is dry and saturated")
+H_1 = 2792;// in kJ/kg
+Phi_g1 = 6.445;// in kJ/kg-K
+x_2 = (Phi_g1-Phi_f2)/(f_g2-Phi_f2);
+H_2 = h_f2 + (x_2*h_fg2);// in kJ/kg
+Eta1 = (H_1-H_2)/(H_1-h_f2);
+disp("Rankine efficiency is "+string(Eta1)+" or "+string(Eta1*100)+" %");
+W1 = (H_1-H_2)*60/100;// in kJ/kg
+Steam_C= Energy_Equivalent/W1;// in kg
+disp(Steam_C,"Steam consumption per kW-hr in kg is :")
+disp("Part (iii) When steam is superheated and temperature is 300°C")
+H_1 = 3039;// in kJ/kg
+Phi_1 = 6.919;// in kJ/kg-K
+x_2 = (Phi_1 - Phi_f2)/(f_g2-Phi_f2);
+H_2 = h_f2 + (x_2 * h_fg2);// in kJ/kg
+Eta = (H_1 - H_2)/(H_1-h_f2);
+disp("Rankine efficiency is "+string(Eta)+" or "+string(Eta*100)+" %");
+W2 = (H_1-H_2)*60/100;// in kJ/kg
+Steam_C= Energy_Equivalent/W2;// in kg
+disp(Steam_C,"Steam consumption per kW-hr in kg is :")
diff --git a/2144/CH6/EX6.33/ex6_33.sce b/2144/CH6/EX6.33/ex6_33.sce new file mode 100755 index 000000000..2d147146b --- /dev/null +++ b/2144/CH6/EX6.33/ex6_33.sce @@ -0,0 +1,18 @@ +// Exa 6.33
+clc;
+clear;
+close;
+// Given data
+T1 = 400;// in degree C
+T1 = T1 + 273;// in K
+T2 = 72.7;// in degree C
+T2 = T2 + 273;// in K
+Eta = ((T1-T2)/T1)*100;// in %
+disp("For carnot cycle : ")
+disp(Eta,"Rankine efficiency in % is : ")
+H1 = 3248;// in kJ/kg
+h_f2 = 304.5;// in kJ/kg
+del_H = 809.2;// in kJ/kg
+Eta = (del_H/(H1-h_f2))*100;// in %
+disp("For Rankine cycle : ")
+disp(Eta,"Rankine efficiency in % is");
diff --git a/2144/CH6/EX6.34/ex6_34.sce b/2144/CH6/EX6.34/ex6_34.sce new file mode 100755 index 000000000..6ded34259 --- /dev/null +++ b/2144/CH6/EX6.34/ex6_34.sce @@ -0,0 +1,27 @@ +// Exa 6.34
+clc;
+clear;
+close;
+// Given data
+P1 = 15;// in bar
+H1 = 3039;// in kJ/kg
+V_g1 = 0.1697;// in m^3/kg
+Phi1 = 6.919;// in kJ/kg-K
+P2_desh = 3.5*10^2;// in kN/m^2
+Phi_g2 = 6.941;// in kJ/kg-K
+Phi_f2 = 1.727;// in kJ/kg-K
+P2 = 0.15*10^2;// in kN/m^2
+h_f2 = 232;// in kJ/kg
+x = (Phi1-Phi_f2)/(Phi_g2 - Phi_f2);
+h_f = 584;// in kJ/kg
+h_fg = 2148;// in kJ/kg
+H2 = h_f + (x*h_fg);// in kJ/kg
+V = 0.5241;// in m^3
+V2=x*V;// in m^3/kg
+W = (H1-H2) + (P2_desh-P2)*V2;//work output of the cycle in kJ/kg
+Eta = W/(H1-h_f2)*100;// in %
+disp(Eta,"The efficiency of the cycle in % is");
+Energy_equivalent= 3600;// in kJ
+S_consumption = Energy_equivalent/W;// in kg
+V = S_consumption* V_g1;// in m^3
+disp(V,"Total volume of steam in m^3 is");
diff --git a/2144/CH6/EX6.4/ex6_4.sce b/2144/CH6/EX6.4/ex6_4.sce new file mode 100755 index 000000000..0bcfe60af --- /dev/null +++ b/2144/CH6/EX6.4/ex6_4.sce @@ -0,0 +1,11 @@ +// Exa 6.4
+clc;
+clear;
+close;
+// Given data
+L = 693.3;// Liquid heat in kJ/kg
+L1 = 125.7;// Liquid heat of feed water in kJ/kg
+m = 2;// mass of water in kg
+Q = m * (L-L1);// in kJ
+disp(Q,"Heat required to raise temperature in kJ is");
+disp("The water is still liquid at the end of the heat supply")
diff --git a/2144/CH6/EX6.5/ex6_5.sce b/2144/CH6/EX6.5/ex6_5.sce new file mode 100755 index 000000000..6f2cd350b --- /dev/null +++ b/2144/CH6/EX6.5/ex6_5.sce @@ -0,0 +1,22 @@ +// Exa 6.5
+clc;
+clear;
+close;
+// Given data
+x = 0.9;
+h_f = 762.2;// in kJ/kg
+h_fg = 2013.8;// in kJ/kg
+H_wet = h_f + (x*h_fg);// in kJ/kg
+En = 125.7;// Enthapy of liquid in kJ/kg
+H_wet = H_wet - En;// in kJ
+disp(H_wet,"When dry fraction is 0.9, Heat required, to convert in kJ is");
+// Part (b) when dry fraction is saturated
+H_sat = h_f + h_fg;// in kJ/kg
+H_sat = H_sat-En;// in kJ
+disp(H_sat,"Heat required when steam is dry and saturated in kJ is");
+C_P = 2.093;// in kJ/kg-K
+t_sup = 300;// in degree C
+t_sat = 179.9;// in degree C
+H_sup = h_f + h_fg + C_P*(t_sup-t_sat);// in kJ
+H_sup1 = H_sup - En;// in kJ
+disp(H_sup1,"Heat required when the steam is super heated in kJ is");
diff --git a/2144/CH6/EX6.6/ex6_6.sce b/2144/CH6/EX6.6/ex6_6.sce new file mode 100755 index 000000000..84bf842a9 --- /dev/null +++ b/2144/CH6/EX6.6/ex6_6.sce @@ -0,0 +1,19 @@ +// Exa 6.6
+clc;
+clear;
+close;
+// Given data
+x = 0.95;
+v_f = 0.001;
+v_g = 0.1238;// in m^3/kg
+V_wet = ((1-x)*v_f)+(x*v_g);// in m^3 correction little diff in ans
+disp(V_wet,"Specific volume of wet steam in m^3 is");
+disp(v_g,"When the steam is dry saturated, the specific volume in m^3/kg is");
+T_sat = 201.3;// in degree C
+T_sat = T_sat + 273;// in K
+T_sup = 300;// in degree C
+T_sup = T_sup + 273;// in K
+V_sup = v_g * (T_sup/T_sat);// in m^3
+disp(V_sup,"When the steam is superheated, the specific volume in m^3 is");
+
+
diff --git a/2144/CH6/EX6.7/ex6_7.sce b/2144/CH6/EX6.7/ex6_7.sce new file mode 100755 index 000000000..1308f0af1 --- /dev/null +++ b/2144/CH6/EX6.7/ex6_7.sce @@ -0,0 +1,22 @@ +// Exa 6.7
+clc;
+clear;
+close;
+// given data
+h_f = 720.7;// in kJ
+h_fg = 2046.6;// in kJ
+v_g = 0.2405;// in m^3
+x = 0.9;
+P = 8*10^2;// in kN/m^2
+U_sat = h_f+x*h_fg-x*v_g*P;// in kJ
+disp(U_sat,"When the steam is wet, the internal energy in kJ is");
+En = 2767.3;// Enthalpy of dry saturated stream
+U_sat = En-(v_g*P);//in kJ/kg
+disp(U_sat,"When the steam is dry and saturated, the internal energy in kJ/kg is");
+C_P = 2.093;
+del_s = 100;// in degree C
+H_sup = h_f + h_fg + (C_P*del_s);// in kJ/kg
+t_sat = 170.4+273;// in K
+V_sup = (v_g*(t_sat+del_s))/t_sat;// in m^3
+U_sup = H_sup - P*V_sup;// in kJ/kg
+disp(U_sup,"When the steam is super heated, the internal energy in kJ/kg is");
diff --git a/2144/CH6/EX6.8/ex6_8.sce b/2144/CH6/EX6.8/ex6_8.sce new file mode 100755 index 000000000..821236c78 --- /dev/null +++ b/2144/CH6/EX6.8/ex6_8.sce @@ -0,0 +1,14 @@ +// Exa 6.8
+clc;
+clear;
+close;
+// Given data
+x = 0.88;// dryness fraction
+h_fg = 2392.7;// in kJ/kg
+H_wet = x * h_fg;// in kJ/kg
+Vs = 14.67;// Specific volume in m^3/kg
+V_wet = x * Vs;// in m^3/kg
+Q = H_wet/V_wet;// in kJ/m^3
+disp(Q,"Heat to be extracted in kJ/m^3 is");
+
+
diff --git a/2144/CH6/EX6.9/ex6_9.sce b/2144/CH6/EX6.9/ex6_9.sce new file mode 100755 index 000000000..f7ba11d21 --- /dev/null +++ b/2144/CH6/EX6.9/ex6_9.sce @@ -0,0 +1,22 @@ +// Exa 6.9
+clc;
+clear;
+close;
+// Given data
+P = 12*10^2;// in kN/m^2
+h_f = 798.1;// in kJ/kg
+h_fg = 1984.5;// in kJ/kg
+x = 0.8;
+H_wet = h_f + (x*h_fg);// in kJ/kg
+v_f = 0.001;// in m^3
+v_g = 0.1684;// in m^3
+V_wet = ((1-x)*v_f) + (x*v_g);// in m^3
+En = H_wet/V_wet;// kJ/m^3
+disp(En,"The enthalpy in kJ/m^3 is");
+U_wet = H_wet - ( V_wet * P );// in kJ
+U_wet1 = (U_wet/V_wet);// in kJ/m^3
+disp(U_wet1,"Internal energy in kJ/m^3 is");
+
+// Note: There is calculation error to find the value of V_wet.( the correct value of V_wet is 0.13492 not 0.1308), so there is some difference between the output of coding and the answer of the book
+
+
|