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 /2339/CH4 | |
| 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 '2339/CH4')
41 files changed, 1196 insertions, 0 deletions
diff --git a/2339/CH4/EX4.1.1/Ex4_1.sce b/2339/CH4/EX4.1.1/Ex4_1.sce new file mode 100755 index 000000000..de6e050b8 --- /dev/null +++ b/2339/CH4/EX4.1.1/Ex4_1.sce @@ -0,0 +1,60 @@ +clc +clear + +//Case 1 +Vg=0.132; +SV=0.12; //Specific Volume +//As SV is less than Vg, steam is wet + +x=SV/Vg; + +printf('\n For Case 1 \n'); +printf('Part of wet steam: %2.2f',x); +printf('\n'); + +//Case 2 +T=200; +Tsat=179.9; //Satuaration Temperature +//Steam is superheated as T > Tsat +D_sh=T-Tsat; + +printf('\n For Case 2 \n'); +printf('Degree of Superheat: %2.1f C',D_sh); +printf('\n'); + +//Case 3 +P=20; //Pressure in bars +Hf=908.8; //kJ/kg +Hfg=1890.7; //kJ/kg +Hg=2799.5; //kJ/kg +H=2650; + +//Steam is wet as Specific enthalpy is less than Hg + +x=(H-Hf)/Hfg; + +printf('\n For Case 3 \n'); +printf('Part of wet steam: %2.2f',x); +printf('\n'); + +//Case 4 +T=150; //in Celcius +SV=0.3928; //Specific Volume in m^3/kg +Vg=0.3928; //in m^3/kg + +printf('\n For Case 4 \n'); +printf('As SV=Vg, steam is dry saturated'); +printf('\n'); + +//Case 5 +P=10; //in bars +S=5.697; +Sf=2.319; +Sfg=4.448; +Sg=6.623; +//As Sample specific entropy is less than Sg and more than Sf, steam is wet + +x=(S-Sf)/Sfg; +printf('\n For Case 5 \n'); +printf('Part of wet steam: %2.1f',x); +printf('\n'); diff --git a/2339/CH4/EX4.10.1/Ex4_10.sce b/2339/CH4/EX4.10.1/Ex4_10.sce new file mode 100755 index 000000000..078237bc9 --- /dev/null +++ b/2339/CH4/EX4.10.1/Ex4_10.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+//Conditions at 4 bar
+P1=4; //in bars
+Hf=604.7; //in kJ/kg
+Hfg=2133.8; //in kJ/kg
+Vg=0.463; //in m^3/kg
+x1=0.9;
+H1=Hf+(x1*Hfg);
+V1=x1*Vg;
+
+//Now at 12 bar pressure
+P2=12; //in bars
+V2=(P1*V1)/P2;
+Vg=0.163; //in m^3/kg
+printf('At 12 bar, V2: %3.3f kJ/kg',V2);
+printf('\n');
+printf('As Vg>V2, steam is wet');
+printf('\n');
+x2=V2/Vg;
+printf('The dryness fraction at 12 bars: %3.2f ',x2);
+printf('\n');
+
+Hf=798.6; //in kJ/kg
+Hfg=1986.2; //in kJ/kg
+H2=Hf+(x2*Hfg);
+printf('The Final enthalpy of steam: %3.1f kJ/kg',H2);
+printf('\n');
diff --git a/2339/CH4/EX4.11.1/Ex4_11.sce b/2339/CH4/EX4.11.1/Ex4_11.sce new file mode 100755 index 000000000..a91524858 --- /dev/null +++ b/2339/CH4/EX4.11.1/Ex4_11.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+//At 20 degree Celcius
+Cpw=4.187; //in kJ/kg
+Tw=20;
+H1=Cpw*Tw;
+
+//At 8 bar condition
+m=4; //mass in kg
+Cps=2.1; //in kJ/kg
+Tsat=170.4+273; //in K
+Hg=2769.1; //in kJ/kg
+Tsup=200+273; //in K
+H2=Hg+(Cps*(Tsup-Tsat));
+Q=m*(H2-H1);
+printf('Heat to be added: %3.1f kJ',Q);
+printf('\n');
diff --git a/2339/CH4/EX4.12.1/Ex4_12.sce b/2339/CH4/EX4.12.1/Ex4_12.sce new file mode 100755 index 000000000..6fd9acf29 --- /dev/null +++ b/2339/CH4/EX4.12.1/Ex4_12.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+//Combined Seperating and Throttling Calorimeter
+m1=2; //mass of water seperated in kg
+m=20.5; //Steam discharged from calorimeter in kg
+mt=m1+m; //Steam inlet in kg
+
+x1=m/(mt); //Dryness fraction
+
+//At 12 bar pressure
+Hf=798.6; //in kJ/kg
+Hfg=1986.2; //in kJ/kg
+
+P_bar=760; //Pressure in mm
+P_fin=5; //Pressure in mm
+P=(P_bar+P_fin)*1.01325/P_bar; //Absolute Pressure
+
+//Now at 1.02 bar
+Cp=2.2; //in kJ/kg K
+Hg=2676.34; //in kJ/kg
+Tsat=99.66+273; //in K
+Tsup=110+273; //in K
+H2=Hg+(Cp*(Tsup-Tsat));
+x2=(H2-Hf)/Hfg;
+x=x1*x2;
+printf('The Dryness Fraction: %3.3f',x);
+printf('\n');
diff --git a/2339/CH4/EX4.13.1/Ex4_13.sce b/2339/CH4/EX4.13.1/Ex4_13.sce new file mode 100755 index 000000000..c5647d7d1 --- /dev/null +++ b/2339/CH4/EX4.13.1/Ex4_13.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+//At 7 bar and 300 Celcius
+P=7; //in bars
+Cps=2.1;
+Tsup=300+273; //in K
+Tsat=165+273; //in K
+Hg=2763.5; //in kJ/kg
+H1=Hg+(Cps*(Tsup-Tsat));
+
+x2=0.9; //Dryness Fraction
+Hf=697.2; //in kJ/kg
+Hfg=2066.3; //in kJ/kg
+H2=Hf+(x2*Hfg);
+m=(H1-Hg)/(Hg-H2);
+printf('The mass flow rate of wet steam: %3.3f kg/kg',m);
+printf('\n');
diff --git a/2339/CH4/EX4.14.1/Ex4_14.sce b/2339/CH4/EX4.14.1/Ex4_14.sce new file mode 100755 index 000000000..ecbe00ad5 --- /dev/null +++ b/2339/CH4/EX4.14.1/Ex4_14.sce @@ -0,0 +1,65 @@ +clc
+clear
+
+//Conditions at 10 bar
+P=10; //in bar
+Tsat=179.9+273; //in K
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Hg=2778.1; //in kJ/kg
+Vg=0.194; //in m^3/kg
+Sf=2.139; //in kJ/kg K
+Sg=6.623; //in kJ/kg K
+Sfg=4.448; //in kJ/kg K
+x=0.91; //Dryness Fraction
+m=3; //in kg
+
+//Now for wet steam
+H=Hf+(x*Hfg);
+H_final=m*H;
+printf('The total Enthalpy: %3.1f kJ',H_final);
+printf('\n');
+V=x*Vg;
+U=H-(P*100*V);
+U_final=m*U;
+printf('The Internal Energy: %3.1f kJ',U_final);
+printf('\n');
+S=Sf+(x*Sfg);
+S_final=m*S;
+printf('The Entropy: %3.3f kJ/K',S_final);
+printf('\n \n');
+
+//Now Case 2
+printf('Now for Case 2 \n');
+Tsat=179.9+273; //in K
+Tsup=200+273; //in K
+Cp=2.1; //in kJ/kg K
+H=Hg+(Cp*(Tsup-Tsat));
+H_final=m*H;
+printf('The Enthalpy: %3.1f kJ',H_final);
+printf('\n');
+Vsup=(Tsup*Vg)/Tsat;
+U=H-(P*100*Vsup);
+U_final=m*U;
+printf('The change in internal energy: %3.1f kJ',U_final);
+printf('\n');
+S=Sg+(Cp*log(Tsup/Tsat));
+S_final=m*S;
+printf('The Entropy: %3.1f kJ/K',S_final);
+printf('\n');
+
+//Now Case 3
+printf('\n Now for case 3 \n');
+H=Hg;
+H_final=m*H; //in kJ
+printf('The total enthalpy: %3.1f kJ',H_final);
+printf('\n');
+V=Vg;
+U=H-(P*100*V);
+U_final=m*U;
+printf('The change in internal energy: %3.1f kJ',U_final);
+printf('\n');
+S=Sg;
+S_final=m*S;
+printf('The total entropy: %3.3f kJ/kg',S_final);
+printf('\n');
diff --git a/2339/CH4/EX4.15.1/Ex4_15.sce b/2339/CH4/EX4.15.1/Ex4_15.sce new file mode 100755 index 000000000..84c9194c5 --- /dev/null +++ b/2339/CH4/EX4.15.1/Ex4_15.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+//At 15 bar condition
+Tsat=198.3+273; //in K
+m=7; //in kg
+Hg=2792.2; //in kJ/kg
+Tsup=300+273; //in K
+Cps=2.1; //in kJ/kg K
+H1=Hg+(Cps*(Tsup-Tsat));
+Cpw=4.187; //in kJ/kg K
+H2=Cpw*50;
+Q=m*(H1-H2);
+printf('The total amount of heat required: %3.1f kJ',Q);
+printf('\n');
+Sg=6.445; //in kJ/kg K
+S2=Sg+(Cps*log(Tsup/Tsat));
+Sf=0.704; //in kJ/kg K
+S1=Sf;
+S=m*(S2-S1);
+printf('The change in Entropy: %3.2f kJ/K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.16.1/Ex4_16.sce b/2339/CH4/EX4.16.1/Ex4_16.sce new file mode 100755 index 000000000..5c9eff104 --- /dev/null +++ b/2339/CH4/EX4.16.1/Ex4_16.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+//Conditions at 10 bar
+P=10; //in bar
+Tsat=179.9+273; //in K
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Hg=2778.1; //in kJ/kg
+Vg=0.194; //in m^3/kg
+x=0.7; //Dryness Fraction
+V=x*Vg;
+m=0.2/V; //mass in kg
+mf=2/V; //mass in kg
+H=Hf+(x*Hfg);
+H_tot=H*mf;
+printf('The total enthalpy: %3.1f kJ',H_tot);
+printf('\n');
+U=H-(P*100*V);
+U_tot=U*mf;
+printf('The internal energy: %3.1f kJ',U_tot);
+printf('\n');
+W=P*100*V;
+W_tot=W*mf;
+printf('The external work of evaporation: %3.1f kJ',W_tot);
+printf('\n');
diff --git a/2339/CH4/EX4.17.1/Ex4_17.sce b/2339/CH4/EX4.17.1/Ex4_17.sce new file mode 100755 index 000000000..fd33bc6df --- /dev/null +++ b/2339/CH4/EX4.17.1/Ex4_17.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+//Conditions at 10 bar pressure
+P=10; //in bar
+Tsat=179.9+273; //in K
+Tsup=350+273;
+x=0.9; //Dryness Fraction
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Hg=2778.1; //in kJ/kg
+Vg=0.194; //in m^3/kg
+Cps=2.1; //in kJ/kg K
+Ha=Hg+(Cps*(Tsup-Tsat));
+Hb=Hf+(x*Hfg);
+H_mix=(Ha+Hb)/2;
+Tsupe=((H_mix-Hg)/Cps)+Tsat;
+Tsuper=Tsupe-273;
+printf('Temperature of superheated steam: %3.0f Celcius',Tsuper);
+printf('\n');
diff --git a/2339/CH4/EX4.18.1/Ex4_18.sce b/2339/CH4/EX4.18.1/Ex4_18.sce new file mode 100755 index 000000000..05ab1560d --- /dev/null +++ b/2339/CH4/EX4.18.1/Ex4_18.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+//Now at 10 bar pressure
+V=1.5; //Volume in m^3
+P=10; //Pressure in bar
+x=0.91; //Dryness fraction
+Vg=0.194; //in m^3/kg
+m=V/Vg;
+
+Vf=x*Vg;
+m_f=V/Vf;
+printf('Amount of water to be placed in container: %2.2f kg',m);
+printf('\n');
+printf('Mass of water required: %2.2f kg',m_f);
+printf('\n');
diff --git a/2339/CH4/EX4.19.1/Ex4_19.sce b/2339/CH4/EX4.19.1/Ex4_19.sce new file mode 100755 index 000000000..ae4ad1e86 --- /dev/null +++ b/2339/CH4/EX4.19.1/Ex4_19.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+//Conditions at 7 bat
+P=7; //in bar
+Tsat=165+273; //in K
+Hf=697.2; //in kJ/kg
+Hfg=2066.3; //in kJ/kg
+Hg=2763.5; //in kJ/kg
+Vg=0.273; //in m^3/kg
+D=0.02; //in m
+vel=17; //in m/s
+Cps=4.187; //in kJ/kg K
+Tw1=25; //in Celcius
+Tw2=100; //in Celcius
+Vfr=(22/7)*D*D*vel*(1/4)*60; //Volume flow rate in m^3/min
+x=0.9; //Dryness Fraction
+V=x*Vg;
+
+Mfr=Vfr/V; //Mass flow rate
+printf('The mass flow rate of steam: %2.2f kg/min',Mfr);
+printf('\n');
+H1=Hf+(x*Hfg);
+H2=Cps*100;
+Mw=(Mfr*(H1-H2))/(Cps*(Tw2-Tw1));
+printf('The mass flow rate of water: %2.2f kg/min',Mw);
+printf('\n');
diff --git a/2339/CH4/EX4.2.1/Ex4_2.sce b/2339/CH4/EX4.2.1/Ex4_2.sce new file mode 100755 index 000000000..dc5ae20a5 --- /dev/null +++ b/2339/CH4/EX4.2.1/Ex4_2.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+//At 10 bar pressure
+P=10; //in bars
+x=0.8;
+Vg=0.194; //in kJ/kg
+W=P*100*x*Vg;
+printf('External Work Done: %3.2f kJ/kg',W);
+printf('\n');
+
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+H=Hf+(x*Hfg);
+U=H-W;
+printf('Internal energy: %3.2f kJ/kg',U);
+printf('\n');
+
+Vf=0.001127; //in m^3/kg
+Uf=Hf-(P*100*Vf);
+Ux=U-Uf;
+printf('Internal Heat of Evaporation: %3.2f kJ/kg',Ux);
+printf('\n');
+
+Sf=2.139; //in kJ/kg K
+Sfg=4.448; //in kJ/kg K
+S=Sf+(x*Sfg);
+printf('Entropy of steam: %3.3f kJ/kg',S);
+printf('\n');
diff --git a/2339/CH4/EX4.20.1/Ex4_20.sce b/2339/CH4/EX4.20.1/Ex4_20.sce new file mode 100755 index 000000000..77fbd7ddc --- /dev/null +++ b/2339/CH4/EX4.20.1/Ex4_20.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+//Conditions at 9 bar
+P=9; //in bar
+Tsat=175.4+273; //in K
+Vg=0.215; //in m^3/kg
+Hf=742.8; //in kJ/kg
+Hfg=2031.1; //in kJ/kg
+Hg=2773.9; //in kJ/kg
+T2=250+273; //in K
+x=0.91; //Dryness Fraction
+V1=x*Vg;
+V2=0.2696; //From Steam Table
+W=P*100*(V2-V1);
+printf('The Work Output: %2.2f kJ/kg',W);
+printf('\n');
+H1=Hf+(x*Hfg);
+H2=2946.3; //From steam table in kJ/kg
+Q=H2-H1;
+printf('The heat supplied to steam: %2.2f kJ/kg',Q);
+printf('\n');
+U=Q-W;
+printf('The internal energy of steam increases by: %2.2f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH4/EX4.21.1/Ex4_21.sce b/2339/CH4/EX4.21.1/Ex4_21.sce new file mode 100755 index 000000000..c56cca3dc --- /dev/null +++ b/2339/CH4/EX4.21.1/Ex4_21.sce @@ -0,0 +1,40 @@ +clc
+clear
+
+//Conditions at 16 bar
+P=16; //in bar
+Vov=0.015; //Volume of Vessel
+Mos=0.1; //Mass of steam
+SV=Vov/Mos; //Specific Volume
+Vg=0.124; //in m^3/kg
+Tsat=201.4+273; //in K
+Tsup=(SV/Vg)*Tsat;
+printf('The temperature of steam: %2.2f K',Tsup);
+printf('\n');
+
+//Now cooling takes place
+Tsat=191.16; //From steam table
+printf('After cooling, temperature of steam: %2.2f K',Tsat);
+printf('\n');
+
+//Now cooled to 10 bar pressure
+P1=16; //in bar
+Vg=0.194; //in m^3/kg
+v=0.15; //in m^3/kg
+x=v/Vg; //Dryness Fraction
+
+//For constant Volume process W=0
+Hg=2794.0; //in kJ/kg
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Cps=2.1; //in kJ/kg K
+Tsup=300.84; //in C
+Tsat=201.4; //in C
+H1=Hg+(Cps*(Tsup-Tsat));
+U1=H1-(P1*100*v);
+P2=10; //in bar
+H2=Hf+(x*Hfg);
+U2=H2-(P2*100*v);
+Q=U2-U1;
+printf('Heat rejected by system: %2.2f kJ/kg',Q);
+printf('\n');
diff --git a/2339/CH4/EX4.22.1/Ex4_22.sce b/2339/CH4/EX4.22.1/Ex4_22.sce new file mode 100755 index 000000000..b7ac9e36f --- /dev/null +++ b/2339/CH4/EX4.22.1/Ex4_22.sce @@ -0,0 +1,73 @@ +clc
+clear
+
+//Isothermal process
+P=10; //in bar
+Tsat=179.9+273; //in K
+Vg=0.194; //in m^3/kg
+Hf=762.6; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Hg=2778.1; //in kJ/kg
+x1=1; //Dryness Fraction
+Sf=2.139; //in kJ/kg K
+Sfg=4.448; //in kJ/kg K
+Sg=6.623; //in kJ/kg K
+V=0.3; //in m^3
+m=V/Vg; //in kg
+V2=Vg/2;
+x2=V2/Vg; //Dryness Fraction
+W=P*100*(V2-Vg)*m;
+printf('Work Done: %2.2f kJ',W);
+printf('\n');
+H1=Hg;
+H2=Hf+(x2*Hfg);
+Q=m*(H2-H1);
+printf('Change in Enthalpy: %2.2f kJ',Q);
+printf('\n');
+U=(Q-W);
+printf('Change in total Internal Energy: %2.2f kJ',U);
+printf('\n');
+S1=Sg;
+S2=Sf+(x2*Sfg);
+S=m*(S2-S1);
+printf('Change in Entropy: %2.2f kJ/K',S);
+printf('\n');
+
+//Now for case 2 where PV=C
+printf('Now for case 2');
+printf('\n');
+
+V01=0.097;
+V02=0.5*V01;
+P1=10; //in bars
+P2=(P1*V01)/V02; //in bars
+
+//Now at 20 bars
+Vg1=0.0996; //in m^3/kg
+V2=0.097;
+x2=V2/Vg1; //Dryness Fraction
+Hf=908.8; //in kJ/kg
+Hfg=1890.7; //in kJ/kg
+H2=Hf+(x2*Hfg);
+H=m*(H2-Hg);
+printf('Change in Enthalpy: %2.2f kJ',H);
+printf('\n');
+
+W=m*P1*100*Vg*(log(V02/V01));
+printf('Total work done: %2.2f kJ',W);
+printf('\n');
+
+U=H; //as P1 V1= P2 V2
+Q=U+W;
+printf('Change in Enthalpy: %2.2f kJ',Q);
+printf('\n');
+
+//Now at 20 bar pressure
+Sf=2.447; //in kJ/kg K
+Sfg=3.894; //in kJ/kg K
+Sg1=6.341; //in kJ/kg K
+S2=Sf+(x2*Sfg);
+S1=Sg;
+S=m*(S2-S1)
+printf('Change in Entropy: %2.3f kJ/K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.23.1/Ex4_23.sce b/2339/CH4/EX4.23.1/Ex4_23.sce new file mode 100755 index 000000000..449342b19 --- /dev/null +++ b/2339/CH4/EX4.23.1/Ex4_23.sce @@ -0,0 +1,49 @@ +clc
+clear
+
+//Initial conditions at 7 bar pressure
+P1=7; //in bars
+Vg1=0.273; //in m^3/kg
+V1=Vg1; //in m^3/kg
+Hg1=2763.5; //in kJ/kg
+H1=Hg1;
+Tsat=165+273; //in K
+Sf=1.992; //in kJ/kg K
+Sfg=4.716; //in kJ/kg K
+Sg=6.708; //in kJ/kg K
+n=1.1;
+
+//Final conditions at 0.5 bar
+P2=0.5; //in bars
+V2=((P1*(V1^1.1))/P2)^(1/1.1); //using P(V)^1.1=Constant
+
+W=((P1*100*V1)-(P2*100*V2))/(n-1);
+printf('Work Done: %3.2f kJ',W);
+printf('\n');
+
+Hf2=340.6; //in kJ/kg
+Hfg2=2305.4; //in kJ/kg
+Vg2=3.24; //in m^3/kg
+x2=V2/Vg2; //Dryness Fraction
+
+H2=Hf2+(x2*Hfg2);
+
+U1=H1-(P1*100*V1);
+U2=H2-(P2*100*V2);
+U=U2-U1;
+printf('Change in Internal Energy: %3.2f kJ/kg',U);
+printf('\n');
+
+Q=U+W; //From First law of Thermodynamics
+printf('Heat Transferred: %3.2f kJ/kg',Q);
+printf('\n');
+
+S1=Sg;
+//At 0.5 bar
+Sf2=1.091; //in kJ/kg K
+Sfg2=6.503; //in kJ/kg K
+Sg2=7.594; //in kJ/kg K
+S2=Sf2+(x2*Sfg2);
+S=S2-S1;
+printf('Change in Entropy: %3.2f kJ/kg K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.24.1/Ex4_24.sce b/2339/CH4/EX4.24.1/Ex4_24.sce new file mode 100755 index 000000000..beec1031f --- /dev/null +++ b/2339/CH4/EX4.24.1/Ex4_24.sce @@ -0,0 +1,38 @@ +clc
+clear
+
+//At state 1
+P1=20; //in bar
+V=2;
+Vg1=0.0996; //in m^3/kg
+Tsat1=212.4+273; //in K
+Tsup1=573; //in K
+V1=Vg1*(Tsup1/Tsat1);
+m=V/V1;
+
+//At state 2
+V2=V1;
+Vg2=V2;
+P2=16.9; //From Steam Table
+
+//Calculations
+Hg1=2799.5; //in kJ/kg
+Cps=2.1; //in kJ/kg K
+H1=m*(Hg1+(Cps*(Tsup1-Tsat1)));
+U1=H1-(P1*100*V);
+
+Hg2=2795.5; //in kJ/kg from Steam table
+H2=m*Hg2;
+U2=H2-(P2*100*V);
+
+Q=U2-U1;
+printf('Heat Transferred: %3.1f kJ',Q);
+printf('\n');
+
+Sg1=6.341; //in kJ/kg K
+S1=Sg1+(Cps*log(Tsup1/Tsat1));
+
+S2=6.4022; //From Steam Table
+S=m*(S2-S1);
+printf('Change in Entropy: %3.3f kJ/K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.25.1/Ex4_25.sce b/2339/CH4/EX4.25.1/Ex4_25.sce new file mode 100755 index 000000000..d5192987f --- /dev/null +++ b/2339/CH4/EX4.25.1/Ex4_25.sce @@ -0,0 +1,46 @@ +clc
+clear
+
+//For Throttling process, H1=H2
+//At 15 bar pressure
+P1=15; //in bar
+Hf1=844.9; //in kJ/kg
+Hfg1=1947.3; //in kJ/kg
+x1=0.73; //Dryness Fraction
+
+//At 1 bar pressure
+P2=1; //in bar
+Hf2=417.5; //in kJ/kg
+Hfg2=2258.0; //in kJ/kg
+Hg2=2675.5; //in kJ/kg
+H2=2266.4; //in kJ/kg
+
+H1=Hf1+(x1*Hfg1);
+x2=(H2-Hf2)/Hfg2;
+
+//Now if x1=0.95
+H1=Hf1+(0.95*Hfg1);
+H2=H1;
+
+//At 1 bar
+Hg=2675.5;
+Cps=2.1;
+x=0.93; //New dryness fraction
+T=(H2-Hg)/Cps; //Temperature difference
+Tsat=99; //in Celcius
+Tsup=Tsat+T;
+printf('Temperature of superheated steam: %3.1f Celcius',Tsup);
+printf('\n');
+
+//Now at 15 bar
+Sf=2.315; //in kJ/kg K
+Sfg=4.130; //in kJ/kg K
+Sg=6.445; //in kJ/kg K
+S1=Sf+(x*Sfg);
+
+//Now at 1 bar
+Sg1=7.360; //in kJ/kg K
+S2=Sg1+(Cps*log((Tsup+273)/(Tsat+273)));
+S=S2-S1;
+printf('Change in Entropy: %3.2f kJ/kg K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.26.1/Ex4_26.sce b/2339/CH4/EX4.26.1/Ex4_26.sce new file mode 100755 index 000000000..43a406b79 --- /dev/null +++ b/2339/CH4/EX4.26.1/Ex4_26.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+//Heat lost by Steam=Heat gained by water and calorimeter
+
+ms=2; //in kg
+Hf1=697.2; //in kJ/kg
+Hfg1=2066.3; //in kJ/kg
+Hf2=146.7; //in kJ/kg
+T2=35; //in Celcius
+T1=15; //in Celcius
+mg=56; //in kg
+Cpw=4.187; //in kJ/kg K
+H_gained=mg*Cpw*(T2-T1);
+x=(((H_gained)/2)+(Hf2-Hf1))/Hfg1;
+printf('The dryness fraction is %2.2f ',x);
+printf('\n');
diff --git a/2339/CH4/EX4.27.1/Ex4_27.sce b/2339/CH4/EX4.27.1/Ex4_27.sce new file mode 100755 index 000000000..fe070887c --- /dev/null +++ b/2339/CH4/EX4.27.1/Ex4_27.sce @@ -0,0 +1,9 @@ +clc
+clear
+
+//Calculating dryness fraction
+Ms=10;
+Mw=1;
+x=(100*Ms)/(Ms+Mw);
+printf('The Dryness Fraction of steam is %2.1f percent',x);
+printf('\n');
diff --git a/2339/CH4/EX4.28.1/Ex4_28.sce b/2339/CH4/EX4.28.1/Ex4_28.sce new file mode 100755 index 000000000..f821d3b39 --- /dev/null +++ b/2339/CH4/EX4.28.1/Ex4_28.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+P1=11; //in bar
+P2=1.1; //in bar
+T2=130+273; //in K
+Cps=2.1; //in kJ/kg K
+
+//At 11 bar
+Hf1=781.3; //in kJ/kg
+Hfg1=2000.4; //in kJ/kg
+
+//At 1.1 bar
+Hg2=2679.7; //in kJ/kg
+Tsat=102.3+273; //in K
+Tsup=130+273;
+
+//Now for throttling process, H1=H2
+H2=Hg2+(Cps*(Tsup-Tsat));
+x=((H2-Hf1)*100)/Hfg1;
+printf('The dryness fraction of steam: %2.1f',x);
+printf('\n');
diff --git a/2339/CH4/EX4.29.1/Ex4_29.sce b/2339/CH4/EX4.29.1/Ex4_29.sce new file mode 100755 index 000000000..c1c5ae46c --- /dev/null +++ b/2339/CH4/EX4.29.1/Ex4_29.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+//Combined seperating and throttling calorimeter
+Ms=5; //in kg
+Mw=0.5; //in kg
+Cps=2.1; //in kJ/kg K
+Man=166.8; //in mm of Hg
+Bar=733.6; //in mm of Hg
+
+x1=Ms/(Ms+Mw);
+P=Man+Bar;
+P_bar=(1.01325*P)/760; //Pressure in bar
+
+//From steam table
+Hf1=742.8; //in kJ/kg
+Hfg1=2031.1; //in kJ/kg
+Tsat=104.8+273; //in K
+Tsup=110.3+273; //in K
+Hg=2683.5; //in kJ/kg
+
+H2=Hg+(Cps*(Tsup-Tsat));
+x2=(H2-Hf1)/Hfg1;
+x=x1*x2;
+printf('The dryness fraction of steam: %2.3f',x);
+printf('\n');
diff --git a/2339/CH4/EX4.3.1/Ex4_3.sce b/2339/CH4/EX4.3.1/Ex4_3.sce new file mode 100755 index 000000000..b1fe4f34a --- /dev/null +++ b/2339/CH4/EX4.3.1/Ex4_3.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+//Condition at 10 bar pressure
+//Steam is wet
+
+x=0.95;
+P=10; //in bars
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+H=Hf+(x*Hfg);
+printf('Enthalpy : %3.2f kJ/kg',H);
+printf('\n');
+
+//Now we calculate Work Done
+Vg=0.194; //in m^3/kg
+W=P*100*x*Vg;
+U=H-W;
+printf('Internal energy: %3.0f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH4/EX4.30.1/Ex4_30.sce b/2339/CH4/EX4.30.1/Ex4_30.sce new file mode 100755 index 000000000..851fb9dca --- /dev/null +++ b/2339/CH4/EX4.30.1/Ex4_30.sce @@ -0,0 +1,33 @@ +clc
+clear
+
+//Combined seperating and throttling calorimeter
+
+Mw=8; //in kg
+M=63; //in kg
+Ms=M-Mw; //in kg
+P1=81.5; //Pressure after throttling in mm
+P2=754; //Barometer reading in mm
+SD=13.6; //Specific Density of Hg
+
+x1=Ms/(Ms+Mw); //Dryness Fraction
+P=(P1/SD)+P2; //Pressure in mm
+P=1.01325; //Pressure in bar
+
+//Now at 7.5 bar pressure
+Hf1=709.2; //in kJ/kg
+Hfg1=2057.0; //in kJ/kg
+
+//Now at 1.01325 bar
+Hg2=2676.0; //in kJ/kg
+Tsat=100+273; //in K
+Cps=2.1; //in kJ/kg K
+Tsup=110+273; //in K
+
+//For throttling H1=H2
+H2=Hg2+(Cps*(Tsup-Tsat));
+x2=(H2-Hf1)/Hfg1;
+
+x=x1*x2;
+printf('The dryness fraction of steam: %2.3f',x);
+printf('\n');
diff --git a/2339/CH4/EX4.31.1/Ex4_31.sce b/2339/CH4/EX4.31.1/Ex4_31.sce new file mode 100755 index 000000000..9fbe8de7a --- /dev/null +++ b/2339/CH4/EX4.31.1/Ex4_31.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+//At 9.2 bar pressure
+x1=0.96; //Dryness Fraction
+Sf1=2.1038; //in kJ/kg K
+Sg1=6.6151; //in kJ/kg K
+
+//At 3.55 bar pressure
+Sf2=1.7327; //in kJ/kg K
+Sg2=6.9358; //in kJ/kg K
+Vg2=0.5173; //in m^3/kg
+
+//Now at 0.36 bar pressure
+Vg3=4.408; //in m^3/kg
+
+S1=Sf1+(x1*(Sg1-Sf1));
+
+//As process is adiabatic
+S2=S1;
+
+//From steam table, Sg=6.9358 > S2
+
+x2=(S2-Sf2)/(Sg2-Sf2);
+V2=x2*Vg2;
+
+//As volume remains constant
+V3=V2;
+x3=V3/Vg3;
+printf('The dryness fraction of steam: %2.3f',x3);
+printf('\n');
diff --git a/2339/CH4/EX4.32.1/Ex4_32.sce b/2339/CH4/EX4.32.1/Ex4_32.sce new file mode 100755 index 000000000..4939a1475 --- /dev/null +++ b/2339/CH4/EX4.32.1/Ex4_32.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+//At 10 bar pressure
+m=1/(0.9*0.194);
+
+Hf1=762.6; //in kJ/kg
+x1=0.9; //Dryness Fraction
+Hfg1=2013.6; //in kJ/kg
+H1=Hf1+(x1*Hfg1);
+
+Hf2=640.1; //in kJ/kg
+Hfg2=2107.4; //in kJ/kg
+x2=(H1-Hf2)/Hfg2;
+Vg2=0.375;
+
+Ms=(1/(x2*Vg2));
+Vg3=0.462;
+//Now mass of steam blown off
+M=m-Ms;
+
+printf('Mass of steam blown off: %2.3f kg',M);
+printf('\n');
+
+V=1; //Volume in m^3
+x3=V/(Ms*Vg3);
+printf('Dryness fraction of steam: %2.3f ',x3);
+printf('\n');
diff --git a/2339/CH4/EX4.33.1/Ex4_33.sce b/2339/CH4/EX4.33.1/Ex4_33.sce new file mode 100755 index 000000000..0ad5367dc --- /dev/null +++ b/2339/CH4/EX4.33.1/Ex4_33.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+//At 25 bar pressure
+P=25; //Pressure in bar
+x=0.8; //Dryness fraction
+Hf=962.1; //in kJ/kg
+Hfg=1841; //in kJ/kg
+Vg=0.0801; //in m^3/kg
+H=Hf+(x*Hfg);
+printf('Enthalpy: %2.1f kJ/kg',H);
+printf('\n');
+
+U=H-(P*100*x*Vg);
+printf('Internal Energy: %2.1f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH4/EX4.34.1/Ex4_34.sce b/2339/CH4/EX4.34.1/Ex4_34.sce new file mode 100755 index 000000000..a848a827c --- /dev/null +++ b/2339/CH4/EX4.34.1/Ex4_34.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+Ms=20; //in kg
+Mw=2; //in kg
+Cps=2.1; //in kJ/kg K
+x1=Ms/(Ms+Mw); //Dryness fraction
+
+//At 12 bar pressure
+Hf1=798.6; //in kJ/kg
+Hfg1=1986.2; //in kJ/kg
+
+//At 1 bar pressure
+Hg2=2675.5; //in kJ/kg
+Tsup=110+273; //in K
+Tsat=99+273; //in K
+
+//For throttling, H1=H2
+H2=Hg2+(Cps*(Tsup-Tsat));
+x2=(H2-Hf1)/Hfg1;
+
+x=x1*x2;
+printf('Dryness fraction of steam: %2.4f kJ',x);
+printf('\n');
diff --git a/2339/CH4/EX4.35.1/Ex4_35.sce b/2339/CH4/EX4.35.1/Ex4_35.sce new file mode 100755 index 000000000..2f3ca2b41 --- /dev/null +++ b/2339/CH4/EX4.35.1/Ex4_35.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+V=0.15; //in m^3
+P=4; //in bar
+x=0.8; //Dryness fraction
+
+//Now at 4 bar pressure
+P=4; //in bar
+Vg=0.463; //in m^3/kg
+
+SV=x*Vg;
+Mos=V/SV; //Mass of Steam
+
+//Now if Volume is 1 m^3
+
+Ms=1/SV; //in kg
+//At 4 bar pressure
+Hf=604.7; //in kJ/kg
+Hfg=2133.8; //in kJ/kg
+H=Ms*(Hf+(x*Hfg));
+printf('Enthalpy of 1 m^3 steam: %2.2f kJ',H);
+printf('\n');
diff --git a/2339/CH4/EX4.36.1/Ex4_36.sce b/2339/CH4/EX4.36.1/Ex4_36.sce new file mode 100755 index 000000000..6a0d8a9b9 --- /dev/null +++ b/2339/CH4/EX4.36.1/Ex4_36.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+P1=9; //in bar
+P2=1; //in bar
+T2=115+273; //in K
+m=1.8; //in kg
+m1=0.2; //in kg
+x1=m/(m+m1); //Dryness fraction
+
+//Now from steam table
+Hf=742.8; //in kJ/kg
+Hfg=2031.1; //in kJ/kg
+Hg=2675.5; //in kJ/kg
+Tsat=99+273; //in K
+Tsup=115+273; //in K
+Cps=2.1; //in kJ/kg K
+H2=Hg+(Cps*(Tsup-Tsat));
+x2=(H2-Hf)/Hfg;
+x=x1*x2;
+printf('The dryness fraction: %2.4f kJ',x);
+printf('\n');
diff --git a/2339/CH4/EX4.37.1/Ex4_37.sce b/2339/CH4/EX4.37.1/Ex4_37.sce new file mode 100755 index 000000000..3b0beac29 --- /dev/null +++ b/2339/CH4/EX4.37.1/Ex4_37.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+m1=0.45; //in kg
+m=7; //in kg
+P1=12; //in bar
+Bar=760; //mm of Hg Barometer reading
+Man=180; //mm of Hg Manometer Reading
+Cps=2.1; //in kJ/kg K
+P=Bar+Man;
+P2=(P*1.01325)/760; //Pressure in bar
+Tsup=140+273; //in K
+x1=m/(m+m1);
+
+//Now at 12 bar pressure
+Hf=798.6; //in kJ/kg
+Hfg=1986.2; //in kJ/kg
+
+//At 1.25 bar pressure
+Hg=2685.3; //in kJ/kg
+Tsat=106+273; //in K
+//For throttling H1=H2
+H2=Hg+(Cps*(Tsup-Tsat));
+x2=(H2-Hf)/Hfg;
+
+x=x1*x2;
+printf('The dryness fraction: %2.3f ',x);
+printf('\n');
diff --git a/2339/CH4/EX4.38.1/Ex4_38.sce b/2339/CH4/EX4.38.1/Ex4_38.sce new file mode 100755 index 000000000..f43c09c32 --- /dev/null +++ b/2339/CH4/EX4.38.1/Ex4_38.sce @@ -0,0 +1,43 @@ +clc
+clear
+
+//Case 1
+P=10; //in bar
+Cps=2.1; //in kJ/kg K
+x=0.85; //Dryness fraction
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+Vg=0.194; //in m^3/kg
+Hg=2778.1; //in kJ/kg
+
+H=Hf+(x*Hfg);
+printf('Case 1: When x=0.85 \n \n');
+printf('Enthalpy of steam: %2.2f kJ',H);
+printf('\n');
+
+U=H-(P*100*x*Vg);
+printf('Internal Energy of steam: %2.2f kJ',U);
+printf('\n');
+
+//Case 2
+H=Hg; //in kJ/kg
+printf('\n \nCase 2: When steam is dry and saturated \n \n');
+printf('Enthalpy of steam: %2.2f kJ',H);
+printf('\n');
+
+U=H-(P*100*Vg);
+printf('Internal Energy of steam: %2.2f kJ',U);
+printf('\n');
+
+//Case 3
+Tsup=300+273; //in K
+Tsat=179.9+273; //in K
+H=Hg+(Cps*(Tsup-Tsat));
+printf('\n \nCase 3: When steam is superheated to 300 C \n \n');
+printf('Enthalpy of steam: %2.2f kJ',H);
+printf('\n');
+
+Vsup=(Tsup/Tsat)*Vg;
+U=H-(P*100*Vsup);
+printf('Internal Energy of steam: %2.2f kJ',U);
+printf('\n');
diff --git a/2339/CH4/EX4.39.1/Ex4_39.sce b/2339/CH4/EX4.39.1/Ex4_39.sce new file mode 100755 index 000000000..6db9d44d3 --- /dev/null +++ b/2339/CH4/EX4.39.1/Ex4_39.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+Ms=5; //in kg
+P=5; //in bar
+Tsup=250+273; //in K
+Cps=2.1; //in kJ/kg K
+Tf=30; //in C
+Cpw=4.187; //in kJ/kg K
+H1=Cpw*Tf;
+
+//At 5 bar pressure
+Tsat=151.9+273; //in K
+Hg=2748.7; //in kJ/kg
+H2=Hg+(Cps*(Tsup-Tsat));
+Q=Ms*(H2-H1);
+printf('Amount of heat required: %2.2f kJ',Q);
+printf('\n');
diff --git a/2339/CH4/EX4.4.1/Ex4_4.sce b/2339/CH4/EX4.4.1/Ex4_4.sce new file mode 100755 index 000000000..a0af08a86 --- /dev/null +++ b/2339/CH4/EX4.4.1/Ex4_4.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+//Condition at pressure 15 bars
+P=15; //in bars
+Hf=844.9; // in kJ/kg
+Hfg=1947.3; //in kJ/kg
+Vg=0.132; //in m^3/kg
+x=0.9; //Dryness fraction
+
+W=P*100*x*Vg;
+printf('External Work Done: %3.2f kJ/kg',W);
+printf('\n');
+H=Hf+(x*Hfg);
+U=H-W;
+printf('Internal Energy: %3.1f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH4/EX4.40.1/Ex4_40.sce b/2339/CH4/EX4.40.1/Ex4_40.sce new file mode 100755 index 000000000..5bb0a7912 --- /dev/null +++ b/2339/CH4/EX4.40.1/Ex4_40.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+Ms=3; //in kg
+Tf=30; //in C
+P=8; //in bar
+Tsup=210+273; //in K
+Cps=2.1; //in kJ/kg K
+Cpw=4.186; //in kJ/kg K
+
+H1=Cpw*Tf;
+
+//At 8 bar pressure
+Tsat=170.4+273; //in K
+Hg=2769.1; //in kJ/kg
+H2=Hg+(Cps*(Tsup-Tsat));
+Q=Ms*(H2-H1);
+printf('Amount of heat required: %2.2f kJ',Q);
+printf('\n');
diff --git a/2339/CH4/EX4.41.1/Ex4_41.sce b/2339/CH4/EX4.41.1/Ex4_41.sce new file mode 100755 index 000000000..33d5a42ef --- /dev/null +++ b/2339/CH4/EX4.41.1/Ex4_41.sce @@ -0,0 +1,36 @@ +clc
+clear
+
+//At 7 bar pressure
+P1=7; //in bar
+P2=1; //in bar
+n=1.1;
+//Now according to law of expansion P(V)^1.1= Constant
+
+Vg1=0.273; //in m^3/kg
+V1=Vg1;
+V2=((P1/P2)^(1/n))*V1;
+
+W=((P1*100*V1)-(P2*100*V2))/(n-1);
+printf('Work Done: %3.1f kJ/kg',W);
+printf('\n');
+
+Hg=2763.5; //in kJ/kg
+H1=Hg;
+Vg=1.694;
+//At 1 bar, Vg=1.694 and as V2<Vg steam is wet
+x=V2/Vg;
+
+Hf=417.5; //in kJ/kg
+Hfg=2258; //in kJ/kg
+H2=Hf+(x*Hfg);
+
+U2=H2-(P2*100*V2);
+U1=H1-(P1*100*V1);
+U=U2-U1;
+printf('Change in Internal Energy: %3.2f kJ/kg',U);
+printf('\n');
+
+Q=U+W;
+printf('Heat transferred during the process: %3.2f kJ/kg',Q);
+printf('\n');
diff --git a/2339/CH4/EX4.5.1/Ex4_5.sce b/2339/CH4/EX4.5.1/Ex4_5.sce new file mode 100755 index 000000000..8b8dc9b62 --- /dev/null +++ b/2339/CH4/EX4.5.1/Ex4_5.sce @@ -0,0 +1,37 @@ +clc
+clear
+
+x=0.9; //Dryness Fraction
+m=1.5; //mass in kg
+Cps=2.1;
+//Condition at 10 bars
+P=10;
+Tsat=179.9; //in Celcius
+T=250; //in Celcius
+Hg=2778.1; //in kJ/kg
+Vg=0.194; //in m^3/kg
+Cps=2.1;
+H1=Hg+(Cps*(T-Tsat));
+Vsup=((T+273)/(Tsat+273))*Vg;
+U1=H1-(P*100*Vsup);
+Sf=2.139; //in kJ/kg K
+Sfg=4.448; //in kJ/kg K
+Sg=6.623; //in kJ/kg K
+S1=Sg+(Cps*log((T+273)/(Tsat+273)));
+
+//Conditions at 2.8 bars
+P2=2.8;
+Hf=551.4; //in kJ/kg
+Hfg=2170.7; //in kJ/kg
+Vg=0.646; //in m^3/kg
+H2=Hf+(x*Hfg);
+U2=H2-(P2*100*x*Vg);
+Sf=1.647; //in kJ/kg K
+Sfg=5.368; //in kJ/kg K
+S2=Sf+(x*Sfg);
+U=m*(U2-U1);
+printf('The change in internal energy: %3.1f kJ/kg',U);
+printf('\n');
+S=S2-S1;
+printf('The change in Entropy: %3.4f kJ/kg K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.6.1/Ex4_6.sce b/2339/CH4/EX4.6.1/Ex4_6.sce new file mode 100755 index 000000000..0497cde26 --- /dev/null +++ b/2339/CH4/EX4.6.1/Ex4_6.sce @@ -0,0 +1,34 @@ +clc
+clear
+
+//Conditions at 8 bar
+P=8; //Pressure in bar
+x=0.9; //dryness fraction
+Hf=721.1; //in kJ/kg
+Hfg=2048.0; //in kJ/kg
+Vg=0.240; //in m^3/kg
+H1=Hf+(x*Hfg);
+V1=x*Vg;
+
+//Enthalpy of superheated steam at 8 bar and 200 Celcius
+Hg=2769.1;
+Cps=2.1;
+Tsup=200+273; //in Celcius
+Tsat=170.4+273; //in Celcius
+H2=Hg+(Cps*(Tsup-Tsat));
+V2=(Vg*Tsup)/Tsat;
+H=H2-H1;
+printf('Heat supplied: %3.1f kJ/kg',H);
+printf('\n');
+W=P*100*(V2-V1);
+printf('Work Done: %3.3f kJ/kg',W);
+printf('\n');
+//At 8 bar
+Sf=2.046; //in kJ/kg K
+Sfg=4.617; //in kJ/kg K
+Sg=6.663; //in kJ/kg K
+S1=Sf+(x*Sfg);
+S2=Sg+(Cps*(log(Tsup/Tsat)));
+S=S2-S1;
+printf('The Enthalpy change during process: %3.1f kJ/kg K',S);
+printf('\n');
diff --git a/2339/CH4/EX4.7.1/Ex4_7.sce b/2339/CH4/EX4.7.1/Ex4_7.sce new file mode 100755 index 000000000..c06cd192c --- /dev/null +++ b/2339/CH4/EX4.7.1/Ex4_7.sce @@ -0,0 +1,30 @@ +clc
+clear
+
+//Conditions at 10 bar
+P1=10; //in bars
+Hg=2778.1; //in kJ/kg
+Tsat=179.9+273; //Temperature in K
+Vg=0.194; //in m^3/kg
+
+//Conditions at 10 bar and 300 Celcius
+Cps=2.1;
+Tsup=300+273;
+H1=Hg+(Cps*(Tsup-Tsat));
+V1=Vg*(Tsup/Tsat);
+U1=H1-(P1*100*V1);
+printf('The Internal energy: %3.1f kJ/kg',U1);
+printf('\n');
+
+//At 1.4 bar and other conditions
+P2=1.4; //in bars
+x=0.8; //Dryness Fraction
+Hf=458.4; //in kJ/kg
+Hfg=2232.0; //in kJ/kg
+Vg=1.237; //in m^3/kg
+H2=Hf+(x*Hfg);
+V2=x*Vg;
+U2=H2-(P2*100*V2);
+U=U2-U1;
+printf('The change in internal energy: %3.1f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH4/EX4.8.1/Ex4_8.sce b/2339/CH4/EX4.8.1/Ex4_8.sce new file mode 100755 index 000000000..9bdb825ae --- /dev/null +++ b/2339/CH4/EX4.8.1/Ex4_8.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+//Conditions at 8 bar
+P=8; //in bars
+x=0.8; //Dryness Fraction
+Hf=721.1; //in kJ/kg
+Hfg=2048.0; //in kJ/kg
+H1=Hf+(x*Hfg);
+H2=H1+410; //After adding 410 kJ of heat
+Hg=2769.1; //in kJ/kg
+printf('The Enthalpy of steam: %3.1f kJ/kg',H2);
+printf('\n');
+printf('The steam is superheated')
+printf('\n');
+V2=0.240; //in m^3/kg
+Vg=V2;
+Den=1/Vg;
+printf('The Density of steam: %3.3f kg/m^3',Den);
+printf('\n');
diff --git a/2339/CH4/EX4.9.1/Ex4_9.sce b/2339/CH4/EX4.9.1/Ex4_9.sce new file mode 100755 index 000000000..85d41dae4 --- /dev/null +++ b/2339/CH4/EX4.9.1/Ex4_9.sce @@ -0,0 +1,14 @@ +clc
+clear
+
+//For throttling H1=H2
+
+//At 11 bar
+Hf=781.3; //in kJ/kg
+Hfg=2000.4; //in kJ/kg
+
+//At 1 bar
+Hg=2675.5; //in kJ/kg
+x=(Hg-Hf)/Hfg;
+printf('The Dryness Fraction: %3.3f kJ/kg',x);
+printf('\n');
|