19 files changed, 307 insertions, 0 deletions
diff --git a/2006/CH5/EX5.1/ex5_1.sce b/2006/CH5/EX5.1/ex5_1.sce
new file mode 100755
index 000000000..fd8daef5d
--- /dev/null
+++ b/2006/CH5/EX5.1/ex5_1.sce
@@ -0,0 +1,5 @@
+clc;
+Q12=-250; // Heat transfer during Discharging of battery in kcal
+W21=-0.53; // Consumption of electricity dring Charging process in kWh
+Q21=(W21*3600)-(Q12*4.1868); // First law of thermodynamics
+disp ("kJ",Q21,"Heat loss from battery during charging process");
diff --git a/2006/CH5/EX5.10/ex5_10.sce b/2006/CH5/EX5.10/ex5_10.sce
new file mode 100755
index 000000000..d9e32fa75
--- /dev/null
+++ b/2006/CH5/EX5.10/ex5_10.sce
@@ -0,0 +1,10 @@
+clc;
+p1=1; // initial pressure of air in piston cylinder arrangement in bar
+T=300; // Temperature of air in piston cylinder arrangement in kelvin
+p2=10; // Final pressure of air in piston cylinder arrangement in bar
+R=0.287; // Characteristic gas constant of air in kJ/kg K
+disp("The change in internal energy during the isothermal process is zero");
+w=R*T*log (p1/p2); // Work done
+disp ("kJ/kg",w,"Work done = ");
+q=w; // From first law of thermodynamics
+disp ("kJ/kg",q,"Heat transfer = ");
diff --git a/2006/CH5/EX5.12/ex5_12.sce b/2006/CH5/EX5.12/ex5_12.sce
new file mode 100755
index 000000000..0c2cd44f5
--- /dev/null
+++ b/2006/CH5/EX5.12/ex5_12.sce
@@ -0,0 +1,32 @@
+clc;
+p1=65;// (Error in textbook) // Pressure of air at state 1 in bar
+v1=0.0135; // Volume of air at state 1 in m^3
+v2=0.1; // Volume of air at state 2 in m^3
+R=0.287; // Characteristic gas constant of air in kJ/kg K
+Cvo=0.7165; // Specific heat at constant volume in kJ/kg K
+// (a).Adiabatic process
+k=1.4; // Index of adiabatic process
+p2=p1*(v1/v2)^k; // Pressure of air at state 2 
+T1=p1*10^2*v1/R;// Tempewrature of air at state 1
+T2=p2*10^2*v2/R;// Tempewrature of air at state 2
+w=R*(T2-T1)/(1-k); // work done
+q=0; // Adiabatic expansion process
+delta_u=Cvo*(T2-T1); // Change in internal energy of air
+disp ("kJ",q,"Heat Ineraction = ","kJ",delta_u,"Change in internal energy of air = ","kJ",w,"Work done = ","K",T2,"Final Temperature = ","(a).Adiabatic Process");
+// (b).Polytropic  process
+n=1.3; // Index of adiabatic process
+p2=p1*(v1/v2)^n; // Pressure of air at state 2 
+T1=p1*10^2*v1/R;// Tempewrature of air at state 1
+T2=p2*10^2*v2/R;// Tempewrature of air at state 2
+w=R*(T2-T1)/(1-n); // work done
+delta_u=Cvo*(T2-T1); // Change in internal energy of air
+q=delta_u+w; // Adiabatic expansion process
+disp ("kJ",q,"Heat Ineraction = ","kJ",delta_u,"Change in internal energy of air = ","kJ",w,"Work done = ","K",T2,"Final Temperature = ","(b).Polytropic Process");
+// (c).Isothermal process
+T1=p1*10^2*v1/R;// Tempewrature of air at state 1
+T2=T1; // Tempewrature of air at state 2
+p2=p1*(v1/v2); // Pressure of air at state 2
+w=R*T1*log (v2/v1); // work done
+delta_u=Cvo*(T2-T1); // Change in internal energy of air
+q=delta_u+w; // Adiabatic expansion process
+disp ("kJ",q,"Heat Ineraction = ","kJ",delta_u,"Change in internal energy of air = ","kJ",w,"Work done = ","K",T2,"Final Temperature = ","(c).Isothermal Process");
diff --git a/2006/CH5/EX5.13/ex5_13.sce b/2006/CH5/EX5.13/ex5_13.sce
new file mode 100755
index 000000000..c6d5b847c
--- /dev/null
+++ b/2006/CH5/EX5.13/ex5_13.sce
@@ -0,0 +1,12 @@
+clc;
+m=10; // mass flow rate of chilled water in kg/s
+V1=50; //velocity of chilled water at section 1 in m/s
+z1=30; // Elevation of section 1 in m
+V2=10; //velocity of chilled water at section 2 in m/s
+z2=60; // Elevation of section 2 in m
+h1=21; // Enthalpy of chilled water at section 1 in kJ/kg
+h2=43; // Enthalpy of chilled water at section 2 in kJ/kg
+W=35; // Rate of work done by pump in kW
+g=9.80665; // Acceleration due to gravity in m/s^2
+Q=m*(h2-h1)+(m*(V2^2-V1^2)/2000)+(m*g*(z2-z1)/1000)-W; // Steady flow energy equation
+disp ("kW",Q,"The rate of Heat Transfer From Building (Error in textbook)= ");
diff --git a/2006/CH5/EX5.14/ex5_14.sce b/2006/CH5/EX5.14/ex5_14.sce
new file mode 100755
index 000000000..335fe328f
--- /dev/null
+++ b/2006/CH5/EX5.14/ex5_14.sce
@@ -0,0 +1,19 @@
+clc;
+p1=1; // Suction pressure of air in bar
+p2=5; // Delivery pressure of air in bar
+T1=310; // Suction Temperature of air in kelvin
+Cpo=1.0035; // Specific heat at constant pressure in kJ/kg K
+// (a) Polytropic compression
+T2=475; // Delivery Temperature of air in kelvin
+Q=-0.15; // Heat loss to the cooling water in kW
+Wpoly=-5.3; // Power consumption of compressor in kW
+m=(-Wpoly+Q)/(Cpo*(T2-T1)); // mass flow rate of air from SSSF energy equation
+n=1/((1-((log (T2/T1))/(log (p2/p1))))); // Index of polytropic process
+disp (n,"Index of polytropic process = ","kg/s",m,"mass flow rate of air = ","(a).Polytropic compression");
+// (b) Adiabatic compression
+k=1.4; // Index of adiabatic process
+Wad=-m*Cpo*T1*((p2/p1)^((k-1)/k)-1); // power consumption of compressor
+disp ("kW",Wad,"Power consumption of compressor = ","(b).Adiabatic compression");
+// (c).Difference between specific work
+wdiff=(abs(Wad)-abs(Wpoly))/m; // Difference between specific work
+disp ("kJ/kg",wdiff,"Difference between specific work = ","(c).Difference between specific work");
diff --git a/2006/CH5/EX5.15/ex5_15.sce b/2006/CH5/EX5.15/ex5_15.sce
new file mode 100755
index 000000000..33775fa72
--- /dev/null
+++ b/2006/CH5/EX5.15/ex5_15.sce
@@ -0,0 +1,13 @@
+clc;
+W=500; // Power output from steam turbine in MW
+Q=10; // Heat loss to surroundings in MW
+p1=12.5; // Pressure of staem at steam turbine inlet in MPa
+p2=10; // Pressure of staem at steam turbine outlet in kPa
+V1=50; // Velocity of steam at steam turbine inlet in m/s
+V2=100; // Velocity of steam at steam turbine outlet in m/s
+x2=0.85; // Quality of steam at steam turbine outlet
+h1=3341.8; // Specific enthalpy of staem at inlet from steam table in kJ/kg
+hf2=191.83; hg2=2584.7;// Specific enthalpies of fluid and steam at outlet from steam table in kJ/kg
+h2=(1-x2)*hf2+x2*hg2;// Specific enthalpy of staem at outlet in kJ/kg
+m=(W-Q)*10^3/((h1-h2)+(V1^2-V2^2)/2000); // Mass flow rate of steam
+disp ("kg/s",m,"Mass flow rate of steam = ");
diff --git a/2006/CH5/EX5.16/ex5_16.sce b/2006/CH5/EX5.16/ex5_16.sce
new file mode 100755
index 000000000..55d20839e
--- /dev/null
+++ b/2006/CH5/EX5.16/ex5_16.sce
@@ -0,0 +1,23 @@
+clc;
+p1=3; // Pressre of air at state 1 in bar
+p2=p1; // constant pressure process
+T1=450; // Temperature of air at state 1 in kelvin
+T2=1250; // Temperature of air at state 2 in kelvin
+T3=1000; // Temperature of air at state 3 in kelvin
+V3=50; // Velocity of air at state 3 in m/s
+T4=800; // Temperature of air at state 4 in kelvin
+Cpo=1.0035; // Specific heat at constant pressure in kJ/kg K
+// (a).Combustion chamber
+q=Cpo*(T2-T1); // Heat added to air
+disp ("kJ/kg   (round off error)",q,"Heat added to air = ","(a).Combustion chamber");
+// (b).Turbine 
+k=1.4; // Index of adiabatic process
+w=Cpo*(T2-T3)-V3^2/2000; // Work done 
+disp ("kJ/kg   (round off error)",w,"Work done = ",("(b).Turbine)"));
+// (c).Nozzle
+V4=sqrt (2*Cpo*10^3*(T3-T4)+V3^2); // Velocity of air leaving the nozzle
+disp ("m/s   (round off error)",V4,"Velocity of air leaving the nozzle = ","(c).Nozzle");
+// (d).Pressure drop
+p3=p2*(T3/T2)^(k/(k-1)); // Pressure of air leaving turbine
+p4=p3*(T4/T3)^(k/(k-1)); // Pressure of air leaving nozzle
+disp ("bar   ",p4,"Pressure of air leaving nozzle = ","bar",p3,"Pressure of air leaving turbine = ","(d).Pressure drop");
diff --git a/2006/CH5/EX5.17/ex5_17.sce b/2006/CH5/EX5.17/ex5_17.sce
new file mode 100755
index 000000000..a1271553e
--- /dev/null
+++ b/2006/CH5/EX5.17/ex5_17.sce
@@ -0,0 +1,13 @@
+clc;
+V=1000; // Speed of aircraft in kmph
+p1=0.35; // Ambient pressure in bar
+T1=258; // Ambient temperature in kelvin
+V1=V*1000/3600; // unit conversion kmph into m/s
+p=1.01325; // Atmospheric pressure in bar
+Cpo=1.0035; // Specific heat at constant pressure in kJ/kg K
+k=1.4; // Index of compression process
+T2=T1+(V1^2)/(2*Cpo*10^3); // The temperature after leaving inlet diffuser
+p2=p1*(T2/T1)^(k/(k-1)); // Pressure after leaving inlet diffuser
+r=p/p2; // Pressre ratio of compressor required for pressurization
+disp ("K",T2,"The temperature after leaving inlet diffuser = ");
+disp (r,"Pressre ratio of compressor required for pressurization = ");
diff --git a/2006/CH5/EX5.18/ex5_18.sce b/2006/CH5/EX5.18/ex5_18.sce
new file mode 100755
index 000000000..26f15cf5a
--- /dev/null
+++ b/2006/CH5/EX5.18/ex5_18.sce
@@ -0,0 +1,16 @@
+clc;
+p1=50; // Pressre of steam at diffuser inlet in kPa
+T1=150; // Temperature of steam at diffuser inlet in degree celcius
+V1=180; // Velocity of steam at diffuser inlet in m/s
+A1=1000; // area of diffuser inlet in cm^2
+V2=90; // Velocity of steam at diffuser outlet in m/s
+p2=1; // Pressre of steam at diffuser outlet in bar
+Q=120; // Heat loss to the surroundings in kW
+v1=3.24; // Specific volume of steam from superheated steam table in m^3/kg at inlet
+h1=2645.9; //  // Specific enthalpy of steam from superheated steam table in m^3/kg at inlet
+m=V1*A1*10^-4/v1; // Mass flow rate of steam
+q=Q/m; // Heat transfer per unit mass of steam
+h2=q+h1+(V1^2-V2^2)/2000; //  Specific enthalpy of steam from SSSF energy equationat outlet
+v2=1.704; // Specific volume of steam from superheated steam table in m^3/kg at outlet
+A2=m*v2/V2; // Area of diffuser exit 
+disp ("cm^2",A2*10^4,"Area of diffuser exit (Error in textbook)= ");
diff --git a/2006/CH5/EX5.19/ex5_19.sce b/2006/CH5/EX5.19/ex5_19.sce
new file mode 100755
index 000000000..93e1a688c
--- /dev/null
+++ b/2006/CH5/EX5.19/ex5_19.sce
@@ -0,0 +1,12 @@
+clc;
+p=0.15; // Pressure of Freon-12 gas in MPa
+T1=10; // Temperature of Freon-12 gas in degree celcius
+T2=35; // Temperature of Freon-12 gas in degree celcius
+h1=196.762; // Specific enthalpy of Freon-12 from table at 1 in kJ/kg 
+h2=69.49; // Specific enthalpy of Freon-12 from table at 2 in kJ/kg 
+h3=178.54; // Specific enthalpy of Freon-12 from table at 3 in kJ/kg 
+// For solving Conservation of mass and SSSF energy equations
+A=[1 -1; h3 -h2];
+B=[1 ; h1];
+M=A\B;
+disp ("kg/s",M(1),"Mass of saturated vapour at 0.15 MPa = ","kg/s",M(2),"Mass of saturated liquid at 35 oC = ");
diff --git a/2006/CH5/EX5.2/ex5_2.sce b/2006/CH5/EX5.2/ex5_2.sce
new file mode 100755
index 000000000..5d9aa959c
--- /dev/null
+++ b/2006/CH5/EX5.2/ex5_2.sce
@@ -0,0 +1,12 @@
+clc;
+m=5; // Mass of water in a tank in kg
+T1=30; // Temperature of water at initial state (1) in degree celcius
+T2=95; // Temperature of water at final state (2) in degree celcius
+Qout=70; // Heat transfer from the water tank to the surrounding air in kJ
+W=75; //Electric energy input to a stirrer inside water in kJ
+mf=32.3; // Mass of fel in bomb in grams
+u1=125.78; // Internal energy of water from steam table (uf at T1) in kJ/kg
+u2=397.88; // Internal energy of water from steam table (uf at T2) in kJ/kg
+Qf=m*(u2-u1)-W+Qout; // From First law of thermodynamics
+qf=Qf/(mf*10^-3); // Heat consumption per unit mass of fuel 
+disp ("kJ/kg",qf,"Heat consumption per unit mass of fuel =");
diff --git a/2006/CH5/EX5.20/ex5_20.sce b/2006/CH5/EX5.20/ex5_20.sce
new file mode 100755
index 000000000..176ff3958
--- /dev/null
+++ b/2006/CH5/EX5.20/ex5_20.sce
@@ -0,0 +1,12 @@
+clc;
+p1=15;// Inlet pressure of condenser in kPa
+x=0.92; // Quality of steam
+ms=25; // mass flow rate of steam in kg/s
+d_t=12; // temperature rise of water 
+Cpw=4.1868; // Specific heat of water in kJ/kg K
+hfg=225.94; h4_3=2599.1; // specific enthalpy in kJ/kg
+h1=(1-x)*hfg+x*h4_3; // specific enthalpy 
+h2=225.94; // specific enthalpy in kJ/kg
+Q=ms*(h1-h2); // Heat transfer rate in condenser
+mw=Q/(Cpw*d_t); // Mass flow rate of water 
+disp ("kg/s",mw,"Mass flow rate of water = ","kW",Q,"Heat transfer rate in condenser = ");
diff --git a/2006/CH5/EX5.3/ex5_3.sce b/2006/CH5/EX5.3/ex5_3.sce
new file mode 100755
index 000000000..84212704b
--- /dev/null
+++ b/2006/CH5/EX5.3/ex5_3.sce
@@ -0,0 +1,21 @@
+clc;
+V=50; // Volume of water in a tank in litres
+T1=120; // Temperature of water at initial state (1) in degree celcius
+x1=0.6; // Dryness fraction at initial state (1)
+T2=-10; // Temperature of water at final state (2) in degree celcius
+vf1=0.00106; // specific volume of water from steam tables at T1 in m^3/kg
+vg1=0.8919; // specific volume of water from steam tables at T1 in m^3/kg
+v1=(1-x1)*vf1+x1*vg1; // Specific volume of misture of liquid and water at state (1)
+m=(V*10^-3)/v1; // Mass of water in the tank
+vs2=0.0010891; // Specific volume of saturated ice at T2 in m^3/kg
+vg2=466.7; // Specific volume of water vapour at T2 in m^3/kg
+v2=v1; // constant specific volume during cooling process
+x2=(v2-vs2)/(vg2-vs2); // Dryness fraction at state (2)
+uf1=503.5; // Specific internal energy at state (1) in kJ/kg
+ug1=2529.3; // Specific internal energy at state (1) in kJ/kg
+us2=-354.09; // Specific internal energy at state (2) in kJ/kg
+ug2=2361.4; // Specific internal energy at state (2) in kJ/kg
+u1=(1-x1)*uf1+x1*ug1; // Total Specific internal energy at state (1) in kJ/kg
+u2=(1-x2)*us2+x2*ug2; // Total Specific internal energy at state (2) in kJ/kg
+Q12=m*(u2-u1); // Heat transfer during cooling pocess
+disp ("kJ",Q12,"Heat transfer during cooling pocess = ");
diff --git a/2006/CH5/EX5.4/ex5_4.sce b/2006/CH5/EX5.4/ex5_4.sce
new file mode 100755
index 000000000..9f2dc7ed3
--- /dev/null
+++ b/2006/CH5/EX5.4/ex5_4.sce
@@ -0,0 +1,24 @@
+clc;
+V1=0.3; // Initial volume of water upto stop 1 in m^3
+p1=1; // Initial pressure of water in bar
+x1=0.2; // Dryness fraction at initial state (1)
+p2=3; // Pressur required to lift the piston in bar
+V4=0.45; // Volume of water upto stop 2 in m^3
+vf1=0.001043; // Specific volume at state (1) from steam table in m^3/kg
+vg1=1.694; // Specific volume at state (1) from steam table in m^3/kg
+v1=vf1+x1*(vg1-vf1); // Total Specific volume at state (1) from steam table in m^3/kg
+m=V1/v1; // Mass of water
+v3=V4/m; // Specific volume at stop 2
+v2=v1; p3=p2; v4=v3; V3=V4; V2=V1; // From process diagram
+//  (a)
+p4=0.361; // Final Pressure at v4 from steam table in Mpa
+disp ("MPa",p4,"Fianl pressure = ","(a)");
+//  (b)
+W14=p2*10^2*(V3-V2); // Work done in process 
+uf1=417.36; // Specific internal energ at initial state in kJ/kg
+ufg1=2088.7; // Specific internal energ at initial state in kJ/kg
+u1=uf1+x1*ufg1;  // Total Specific internal energr at initial state in kJ.kg
+u4=2550.2; // Specific internal energ at final state in kJ/kg
+Q14=m*(u4-u1)+W14; // From first law of thermodynamics
+disp ("kJ",W14,"Work done during the process = ","(b)");
+disp ("kJ",Q14,"Heat transfer during the process = ");
diff --git a/2006/CH5/EX5.5/ex5_5.sce b/2006/CH5/EX5.5/ex5_5.sce
new file mode 100755
index 000000000..e69c58fa0
--- /dev/null
+++ b/2006/CH5/EX5.5/ex5_5.sce
@@ -0,0 +1,16 @@
+clc;
+V1=0.01; // Initial Volume of Freon 12 vapour in cylinder in m^3
+T1=15; // Initial Temperature of Freon 12 vapour in degree celcius
+p1=4.914; // Initial pressure (Psat at T1) in bar
+p2=9; // Final pressure of Freon 12 vapour after compression in bar
+T2=65; // Final temperature of Freon 12 vapour after compression in degree celcius
+Q=-0.5; // Heat lost to surroundings during compresson process in kJ
+v1=0.035413; // Initial specific volume of Freon 12 vapour from table in m^3/kg
+m=V1/v1; // Mass of vapour
+hg1=193.644; // specific enthalpy of Freon 12 vopour at state 1 in kJ/kg
+u1=hg1-(p1*10^2*v1); // Total Specific internal energy at state 1
+h2=222.9; // specific enthalpy of Freon 12 vapour at state 2 in kJ/kg
+v2=0.022537; // specific volume of Freon 12 vapour at state 2 in m^3/kg
+u2=h2-(p2*10^2*v2); // Total Specific internal energy at state 2
+W=-m*(u2-u1)+Q; // From first law of thermodynamics
+disp ("kJ",W,"Work of compression = ");
diff --git a/2006/CH5/EX5.6/ex5_6.sce b/2006/CH5/EX5.6/ex5_6.sce
new file mode 100755
index 000000000..69fbdc2bb
--- /dev/null
+++ b/2006/CH5/EX5.6/ex5_6.sce
@@ -0,0 +1,12 @@
+clc;
+tS=-10; // initial temperature of ice in degree celcius
+ts=0; // Melting temperature of ice in degree celcius
+Cs=1.94; // Specific heat of ice in kJ/kg k
+m=1000; // Mass of ice in kg
+tF=10; // Temperature of water after 24 hours in degree celcius
+tf=0; // Freezing temperature of water in degree celcius
+Cf=4.1868; // Specific heat of water in kJ/kg K
+hsg=335; // Latent heat of fusion of ice in kJ/kg
+Q=m*(Cs*(ts-tS)+hsg+Cf*(tF-tf)); // Heat gain of water
+Q_dot=Q/(24*3600); // Rate of cooling
+disp ("kW",Q_dot,"Rate of cooling = ");
diff --git a/2006/CH5/EX5.7/ex5_7.sce b/2006/CH5/EX5.7/ex5_7.sce
new file mode 100755
index 000000000..df242d37b
--- /dev/null
+++ b/2006/CH5/EX5.7/ex5_7.sce
@@ -0,0 +1,12 @@
+clc;
+T1=300; // Temperature of air at state 1 in Kelvin
+T2=500; // Temperature of air at state 2 in Kelvin
+m_=28.966; // Molecular weight oh air in kg
+Cpoav=1.017; // Average value of specific heat of air in kJ/kg K
+// (a). change in enthalpy
+h_=27.43*(T2-T1)+3.09*10^-3*(T2^2-T1^2)-0.2296*10^-6*(T2^3-T1^3); //change in enthalpy during process in kJ/kmol
+h=h_/m_; // change in enthalpy during process in kJ/kg
+disp ("kJ/kg",h,"(a).change in enthalpy during process = ");
+// (b).change in enthalpy
+h=Cpoav*(T2-T1); // change in enthalpy in kJ/kg
+disp ("kJ/kg",h,"(b).change in enthalpy during process with average specfic heat = ");
diff --git a/2006/CH5/EX5.8/ex5_8.sce b/2006/CH5/EX5.8/ex5_8.sce
new file mode 100755
index 000000000..8387e5fcd
--- /dev/null
+++ b/2006/CH5/EX5.8/ex5_8.sce
@@ -0,0 +1,14 @@
+clc;
+m=0.1; // mass of nitrogen gas in kg
+V1=0.1; // Initial volme of nitrogen gas in m^3
+p1=1.2; // Initial pressure of nitrogen gas in bar
+V2=0.075; // Final volume of nitrogen gas in m^3
+Cpo=1.041; // Specific heat at constant pressure of nitrogen in kJ/kg K
+R=0.2969393; // Characteristic gas constant of nitrogen in Kj/kg K
+T1=(p1*10^2*V1)/(m*R); // Initial temperature of nitrogen gas 
+T2=T1*(V2/V1); // Final temperature of nitrogen gas (constant pressure process)
+Q=m*Cpo*(T2-T1); // Heat transfer to surroundings
+W=p1*10^2*(V2-V1); // Work done
+disp ("K",T2,"Final Temperature of nitrogen gas = ");
+disp ("kJ",Q,"Heat transfer to surroundings = ");
+disp ("The Work is done on the gas","kJ",W,"Work done = ");
diff --git a/2006/CH5/EX5.9/ex5_9.sce b/2006/CH5/EX5.9/ex5_9.sce
new file mode 100755
index 000000000..7b464b297
--- /dev/null
+++ b/2006/CH5/EX5.9/ex5_9.sce
@@ -0,0 +1,29 @@
+clc;
+p=1; // pressure inside piston cylinder arrangement in MPa
+// stae 1 = saturated liquid
+// state 2 = saturated vapour
+// state 3 = superheated vapour
+v1=0.001127; // specific volume at state 1 in m^3/kg
+v2=0.19444; // specific volume at state 2 in m^3/kg
+v3=0.4011; // specific volume at state 3 in m^3/kg
+u1=761.68; // specific internal energy at state 1 in kK/kg
+u2=2583.6; // specific internal energy at state 2 in kK/kg
+u3=3296.8; // specific internal energy at state 3 in kK/kg
+h1=762.81; // specific enthalpy at state 1 in kJ/kg
+h2=2778.1; // specific enthalpy at state 2 in kJ/kg
+h3=3697.9; // specific enthalpy at state 3 in kJ/kg
+w12=p*10^3*(v2-v1); // Work done during process 1-2
+w23=p*10^3*(v3-v2); // Work done during process 2-3
+wtotal=w12+w23; // Total work done
+disp ("kJ/kg",wtotal,"Work done = ");
+// Calculation of heat transfer
+// Method I
+q12=(u2-u1)+w12; // Heat transfer during process 1-2
+q23=(u3-u2)+w23; // Heat transfer during process 2-3
+qtotal=q12+q23; // Total Heat transfer
+disp ("kJ/kg",qtotal,"Heat Transfer = ","Method I","Calculation of Heat Transfer");
+// Method II
+q12=h2-h1; // Heat transfer during process 1-2
+q23=h3-h2; // Heat transfer during process 2-3
+qtotal=q12+q23; // Total Heat transfer
+disp ("kJ/kg",qtotal,"Heat Transfer = ","Method II");