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 | |
| 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')
207 files changed, 5223 insertions, 0 deletions
diff --git a/2339/CH1/EX1.1.1/Ex1_1.sce b/2339/CH1/EX1.1.1/Ex1_1.sce new file mode 100755 index 000000000..4ae97ea33 --- /dev/null +++ b/2339/CH1/EX1.1.1/Ex1_1.sce @@ -0,0 +1,6 @@ +clc +clear + +//Solving the problem +Q=84-8.4-21+4.2; +printf('The Net Work Done= %2.1f kJ',Q); //Displaying result diff --git a/2339/CH1/EX1.1.2/Ex1_2.sce b/2339/CH1/EX1.1.2/Ex1_2.sce new file mode 100755 index 000000000..3e1f36005 --- /dev/null +++ b/2339/CH1/EX1.1.2/Ex1_2.sce @@ -0,0 +1,10 @@ +clc +clear + +//Declaring values +Q=-700; +W=-3000; +m=5; +U=Q-W; +Us=U/m; +printf('Change in Specific Energy= %3.0f J/kg',Us); //displaying result diff --git a/2339/CH1/EX1.1.3/Ex1_3.sce b/2339/CH1/EX1.1.3/Ex1_3.sce new file mode 100755 index 000000000..8ef58d5e5 --- /dev/null +++ b/2339/CH1/EX1.1.3/Ex1_3.sce @@ -0,0 +1,8 @@ +clc +clear + +//Declaring values +Q=50; +W=40; +U=Q-W; +printf('Change in Internal Energy= %2.0f kJ',U); diff --git a/2339/CH1/EX1.1.4/Ex1_4.sce b/2339/CH1/EX1.1.4/Ex1_4.sce new file mode 100755 index 000000000..cc03dbed9 --- /dev/null +++ b/2339/CH1/EX1.1.4/Ex1_4.sce @@ -0,0 +1,11 @@ +clc +clear + +m=3000; //mass in kg +P=736; //Power in kW +t=5*3600; //Time in seconds +HV=27170 //Heating value in kJ/kg + +E=P/((m/t)*HV); +Eff=E*100; +printf('Thermal Efficiency= %2.2f percent',Eff); diff --git a/2339/CH1/EX1.1.5/Ex1_5.sce b/2339/CH1/EX1.1.5/Ex1_5.sce new file mode 100755 index 000000000..a41386bd1 --- /dev/null +++ b/2339/CH1/EX1.1.5/Ex1_5.sce @@ -0,0 +1,12 @@ +clc +clear + +U=22; //Internal Energy in kJ/s +P2=0.95*1000; //Pressure in kPa +V2=0.09; //Volume in m^3/s; +P1=0.5*1000; +V1=0.15; +X=(P2*V2)-(P1*V1); +H=U+X; + +printf('Change in Enthalpy: %2.1f kJ/s',H); diff --git a/2339/CH1/EX1.1.6/Ex1_6.sce b/2339/CH1/EX1.1.6/Ex1_6.sce new file mode 100755 index 000000000..b1ec7dc7e --- /dev/null +++ b/2339/CH1/EX1.1.6/Ex1_6.sce @@ -0,0 +1,15 @@ +clc +clear + +Th=0.22; //Thermal Efficiency +Hr=1260; //Heat Rejected in MJ/hr +CV=42; //Calorific Value of Coal +X=1-Th; +HI=Hr/X; //Heat Input in MJ/hr + +O=((HI-Hr)*1000)/3600; //Output +Mf=HI/CV; //Mass of Fuel Used + +printf('Power Output is %2.2f kW',O); +printf('\n'); +printf('Mass of Fuel used per hour: %2.1f kg/hr',Mf); diff --git a/2339/CH1/EX1.1.7/Ex1_7.sce b/2339/CH1/EX1.1.7/Ex1_7.sce new file mode 100755 index 000000000..61849dd79 --- /dev/null +++ b/2339/CH1/EX1.1.7/Ex1_7.sce @@ -0,0 +1,11 @@ +clc +clear + +m=2; //mass in kg +T1=30+273; //Temperature in K +T2=60+273; +Cp=4.187; +T=T2/T1; +X=double(log(T)); +S=m*Cp*X; +printf('Entropy Change of Water: %1.4f kJ/K',S); diff --git a/2339/CH1/EX1.1.8/Ex1_8.sce b/2339/CH1/EX1.1.8/Ex1_8.sce new file mode 100755 index 000000000..3ae88ac7f --- /dev/null +++ b/2339/CH1/EX1.1.8/Ex1_8.sce @@ -0,0 +1,16 @@ +clc +clear + +//Declaring Values +m=600; //Mass in kg +z=50000; //Distance in meters +V=2500000; //Velocity in m/hr +g=7.9; //Gravitational Field in m/s^2 +Vel=V/3600; +KE=(0.5*m*Vel*Vel)/1000000; //Kinetic Energy in MJ +PE=(m*g*z)/1000000; //Potential Energy in MJ + +//Displaying Results +printf('The Kinetic Energy is %3.2f MJ',KE); +printf('\n') +printf('The Potential Energy is %3.2f MJ',PE); diff --git a/2339/CH13/EX13.1.1/Ex13_1.sce b/2339/CH13/EX13.1.1/Ex13_1.sce new file mode 100755 index 000000000..74da46e46 --- /dev/null +++ b/2339/CH13/EX13.1.1/Ex13_1.sce @@ -0,0 +1,10 @@ +clc +clear + +N1=250; +D1=53; +D2=32; + +N2=N1*(D1/D2); +printf('Speed of shaft: %2.2f RPM',N2); +printf('\n'); diff --git a/2339/CH13/EX13.2.1/Ex13_2.sce b/2339/CH13/EX13.2.1/Ex13_2.sce new file mode 100755 index 000000000..a9ff24b37 --- /dev/null +++ b/2339/CH13/EX13.2.1/Ex13_2.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+D1=600;
+D2=300;
+N1=100;
+VR=D1/D2;
+N2=VR*N1;
+
+printf('Case One \n');
+printf('Velocity Ratio= %2.2f',VR);
+printf('\n');
+printf('Speed of driven shaft= %2.2f RPM',N2);
+printf('\n\n');
+
+printf('Case Two \n');
+VR=(D1+5)/(D2+5);
+N2=VR*N1;
+printf('Velocity Ratio= %2.2f',VR);
+printf('\n');
+printf('Speed of driven shaft= %2.2f RPM',N2);
+printf('\n\n');
+
+printf('Case Three \n');
+S=4;
+VR=[(D1+5)/(D2+5)]*[(100-S)/100];
+N2=VR*N1;
+printf('Velocity Ratio= %2.2f',VR);
+printf('\n');
+printf('Speed of driven shaft= %2.2f RPM',N2);
+printf('\n\n');
diff --git a/2339/CH13/EX13.3.1/Ex13_3.sce b/2339/CH13/EX13.3.1/Ex13_3.sce new file mode 100755 index 000000000..0104a021e --- /dev/null +++ b/2339/CH13/EX13.3.1/Ex13_3.sce @@ -0,0 +1,13 @@ +clc
+clear
+
+D1=0.3;
+D2=0.2;
+C=3;
+
+L1=[(22/7)*(1/2)*(D1+D2)]+[((D1+D2)^2)/(4*C)]+(2*C);
+L2=[(22/7)*(1/2)*(D1+D2)]+[((D1-D2)^2)/(4*C)]+(2*C);
+
+L=L2-L1;
+printf('The belt length is to be reduced by %2.4f mm',(0-L)*1000);
+printf('\n');
diff --git a/2339/CH13/EX13.4.1/Ex13_4.sce b/2339/CH13/EX13.4.1/Ex13_4.sce new file mode 100755 index 000000000..fc61f49ff --- /dev/null +++ b/2339/CH13/EX13.4.1/Ex13_4.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+D=1;
+P=5000;
+N=250;
+Mew=0.25;
+PP=20;
+Theta=170*(22/7)*(1/180);
+V=((22/7)*D*N)/60;
+
+T12=exp(Mew*Theta)-1;
+T2=(P/(V*T12));
+T1=(T12+1)*T2;
+W=T1/PP;
+
+printf('Width of belt= %2.2f mm',W);
+printf('\n');
diff --git a/2339/CH13/EX13.5.1/Ex13_5.sce b/2339/CH13/EX13.5.1/Ex13_5.sce new file mode 100755 index 000000000..dec97ea9b --- /dev/null +++ b/2339/CH13/EX13.5.1/Ex13_5.sce @@ -0,0 +1,15 @@ +clc
+clear
+
+N1=1000;
+Z1=30;
+Z2=45;
+Z3=75;
+
+N13=Z3/Z1;
+N3=N1/N13;
+
+printf('Velocity Ratio of gear train= %2.1f ',N13);
+printf('\n');
+printf('N3= %2.1f RPM',N3);
+printf('\n');
diff --git a/2339/CH13/EX13.6.1/Ex13_6.sce b/2339/CH13/EX13.6.1/Ex13_6.sce new file mode 100755 index 000000000..29af54562 --- /dev/null +++ b/2339/CH13/EX13.6.1/Ex13_6.sce @@ -0,0 +1,13 @@ +clc
+clear
+
+Na=600;
+Za=25;
+Zb=50;
+Zc=20;
+Zd=40;
+Nad=(Zb/Za)*(Zd/Zc);
+Nd=Na/Nad;
+
+printf('Speed of Output Shaft= %2.1f RPM',Nd);
+printf('\n');
diff --git a/2339/CH3/EX3.1.1/Ex3_1.sce b/2339/CH3/EX3.1.1/Ex3_1.sce new file mode 100755 index 000000000..bf26befad --- /dev/null +++ b/2339/CH3/EX3.1.1/Ex3_1.sce @@ -0,0 +1,24 @@ +clc +clear + +//Declaring Values +V=3; //Volume in m^3 +P1=2500; //Pressure in kilobar +P2=1500; +T2=21+273; //Temperature in Kelvin +T1=(T2*P1)/P2; +Cp=1.005; +Cv=0.718; +R=Cp-Cv; //Universal Gas Constant +m=(P1*V)/(R*T1); //Calculating mass +H=m*Cp*(T2-T1); +U=m*Cv*(T2-T1); +Q=U; //Since Constant Volume Process: Work Done=0 +//Displaying Results +printf('Change in Enthalpy: %5.2f kJ',H); +printf('\n'); +printf('Change in Internal Energy: %5.2f kJ',U); +printf('\n'); +printf('Heat Transfer: %4.2f kJ',Q); +printf('\n') +printf('As Answer is negative, system rejects heat'); diff --git a/2339/CH3/EX3.13.1/Ex3_13.sce b/2339/CH3/EX3.13.1/Ex3_13.sce new file mode 100755 index 000000000..aec25a480 --- /dev/null +++ b/2339/CH3/EX3.13.1/Ex3_13.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=1;
+P2=25;
+P1=1;
+pV=260;
+T1=17+273;
+T2=T1;
+V1=(pV*T1)/(P1*100000);
+printf('As process is Isothermal, Initial and Final Temperatures are same \n');
+printf('The Final Temperature: %3.0f K',T1);
+printf('\n');
+V2=(pV*T2)/(P2*100000);
+printf('The Final Volume: %3.5f m^3',V2);
+printf('\n');
+CR=P2/P1;
+printf('The Compression Ratio: %3.0f ',CR);
+printf('\n');
+printf('Change in Enthalpy is zero as it is Isothermal process \n');
+W=P1*100*V1*(log(P1/P2));
+printf('Work Done is: %3.1f kJ',W);
+printf('\n');
diff --git a/2339/CH3/EX3.14.1/Ex3_14.sce b/2339/CH3/EX3.14.1/Ex3_14.sce new file mode 100755 index 000000000..62a3005d5 --- /dev/null +++ b/2339/CH3/EX3.14.1/Ex3_14.sce @@ -0,0 +1,36 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P2=6;
+Cp=1.75;
+P1=1;
+n=1.3;
+T1=30+273;
+M=30;
+m=2;
+Ro=8314.4;
+R=(Ro/M)/1000;
+printf('The Gas Constant: %3.3f kJ/kg K', R);
+printf('\n');
+Cv=Cp-R;
+G=Cp/Cv;
+printf('The value of Gamma: %1.2f ',G);
+printf('\n');
+T2=(T1)*((P2/P1)^((n-1)/n));
+printf('Final Temperature: %3.2f K',T2);
+printf('\n');
+W=(m*R*(T1-T2))/(n-1);
+printf('The work done on the gas: %3.2f kJ',W);
+printf('\n');
+Q=((G-n)/(G-1))*W;
+printf('The Heat Transfer is %3.2f kJ',Q);
+printf('\n');
+U=m*Cv*(T2-T1);
+printf('The change in Internal Energy is %3.2f kJ',U);
+printf('\n');
diff --git a/2339/CH3/EX3.15.1/Ex3_15.sce b/2339/CH3/EX3.15.1/Ex3_15.sce new file mode 100755 index 000000000..a2b7b8e42 --- /dev/null +++ b/2339/CH3/EX3.15.1/Ex3_15.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P1=350;
+P2=130;
+R=0.287;
+T1=450;
+G=1.4;
+m=1;
+T2=T1*((P2/P1)^((G-1)/G));
+W=(m*R*(T1-T2))/(G-1);
+printf('Amount of External Work done: %3.2f kJ/kg',W);
+printf('\n');
+U=-W;
+printf('Change in Internal Energy: %3.2f kJ/kg',U);
+printf('\n');
diff --git a/2339/CH3/EX3.16.1/Ex3_16.sce b/2339/CH3/EX3.16.1/Ex3_16.sce new file mode 100755 index 000000000..2348c61e4 --- /dev/null +++ b/2339/CH3/EX3.16.1/Ex3_16.sce @@ -0,0 +1,18 @@ +clc +clear + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K +P1=1.5; +V1=0.1; +V2=0.04; +P2=7.5; +T1=30+273; + +//Using ideal gas equation +T2=(P2*V2*T1)/(P1*V1); +printf('The Value of Temperature of gas: %3.0f K',T2); diff --git a/2339/CH3/EX3.17.1/Ex3_17.sce b/2339/CH3/EX3.17.1/Ex3_17.sce new file mode 100755 index 000000000..3ceec93cc --- /dev/null +++ b/2339/CH3/EX3.17.1/Ex3_17.sce @@ -0,0 +1,32 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P1=1.5;
+V1=3;
+T1=27+273;
+P2=30;
+T2=60+273;
+R=0.287;
+
+m1=(P1*100*V1)/(R*T1);
+m2=(P2*100*V1)/(R*T2);
+
+m=m2-m1;
+printf('The mass pumped: %2.2f kg',m);
+printf('\n');
+V=(m*R*(17+273))/(1*100);
+
+printf('Volume: %2.2f m^3',V);
+printf('\n');
+
+T3=27+273;
+P3=(T3*P2)/T2;
+
+printf('Final air pressure in the vessel: %2.2f bar',P3);
+printf('\n');
diff --git a/2339/CH3/EX3.18.1/Ex3_18.sce b/2339/CH3/EX3.18.1/Ex3_18.sce new file mode 100755 index 000000000..97010da21 --- /dev/null +++ b/2339/CH3/EX3.18.1/Ex3_18.sce @@ -0,0 +1,26 @@ +clc +clear + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K +V1=1.5; +m=2; +T1=27+273; +T2=207+273; +V2=V1; +M=28; +Ro=8.314; + +R=Ro/M; + +P1=(m*R*T1)/V1; +printf('The initial pressure of gas: %3.3f bar',P1/100); +printf('\n'); + +P2=(P1*T2)/T1; +printf('The final pressure of gas: %3.3f bar',P2/100); +printf('\n'); diff --git a/2339/CH3/EX3.19.1/Ex3_19.sce b/2339/CH3/EX3.19.1/Ex3_19.sce new file mode 100755 index 000000000..c1fadba4f --- /dev/null +++ b/2339/CH3/EX3.19.1/Ex3_19.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+T1=27+273;
+V1=0.06;
+P1=150;
+Ro=8.314;
+M=28;
+
+R=Ro/M;
+m=(P1*100*V1)/(R*T1);
+printf('Mass of gas at design condition: %2.1f kg',m);
+printf('\n');
+
+P2=170;
+T2=(T1*P2)/P1;
+printf('Fusible plug should melt at: %3.0f K',T2);
+printf('\n');
diff --git a/2339/CH3/EX3.2.1/Ex3_2.sce b/2339/CH3/EX3.2.1/Ex3_2.sce new file mode 100755 index 000000000..9ff0e53e9 --- /dev/null +++ b/2339/CH3/EX3.2.1/Ex3_2.sce @@ -0,0 +1,28 @@ +clc +clear + +//Inputs +m=1; //Mass in kg +P1=7; //Pressure in bar +T1=90+273; //Temperature in K +P2=1.4; +R=0.287; +G=1.4; //Gamma for air +//Calculations according to data required +x=P2/P1; +y=0.1/1.1; +z=x^y; +T2=T1*z; //calculating T2 +printf('Final Temperature is: %3.1f K',T2); +printf('\n'); +W=(m*R*(T1-T2))/(1.1-1); +printf('Work Done is: %3.1f kJ',W); +printf('\n'); +Cv=(R)/(G-1); +Cp=R+Cv; +CI=m*Cv*(T2-T1); +printf('Change in Internal Energy is: %3.2f kJ',CI); +printf('\n'); +Q=CI+W; +printf('Heat Transfer is: %3.2f kJ',Q); +printf('\n'); diff --git a/2339/CH3/EX3.20.1/Ex3_20.sce b/2339/CH3/EX3.20.1/Ex3_20.sce new file mode 100755 index 000000000..f7c830c37 --- /dev/null +++ b/2339/CH3/EX3.20.1/Ex3_20.sce @@ -0,0 +1,36 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P1=7;
+m=3.7;
+V1=1.5;
+P2=1.4;
+V2=4.5;
+U=648;
+Cv=1.05;
+
+RT1=(P1*100*V1)/(m);
+RT2=(P2*100*V2)/(m);
+
+RT=RT1-RT2;
+T=(U)/(m*Cv);
+R=(RT/T)
+printf('The value of R: %1.3f kJ/kg K',R);
+printf('\n');
+
+Cp=Cv+R;
+H=m*Cp*(-T);
+printf('The change in enthalpy is: %4.2f kJ',H);
+printf('\n');
+
+T1=RT1/R;
+T2=RT2/R;
+
+printf('Initial Temperature: %3.2f K \n',T1);
+printf('Final Temperature: %3.2f K \n',T2);
diff --git a/2339/CH3/EX3.21.1/Ex3_21.sce b/2339/CH3/EX3.21.1/Ex3_21.sce new file mode 100755 index 000000000..ff2653d4c --- /dev/null +++ b/2339/CH3/EX3.21.1/Ex3_21.sce @@ -0,0 +1,23 @@ +clc +clear + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K +V=1.6; +P=1; +m=2; +T=17+273; +G=1.4; + +R=(P*100*V)/(m*T); +Cv=(R)/(G-1); +printf('The Value of Cv: %1.2f kJ/kg K',Cv); +printf('\n'); + +Cp=Cv+R; +printf('The Value of Cp: %1.3f kJ/kg K',Cp); +printf('\n') diff --git a/2339/CH3/EX3.22.1/Ex3_22.sce b/2339/CH3/EX3.22.1/Ex3_22.sce new file mode 100755 index 000000000..e0fd6156f --- /dev/null +++ b/2339/CH3/EX3.22.1/Ex3_22.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.091;
+P1=2.73;
+T1=187+273;
+T2=27+273;
+Cp=1.005;
+Cv=0.718;
+
+R=Cp-Cv;
+
+m=(P1*100*V1)/(R*T1);
+Q=m*Cp*(T2-T1);
+printf('The Value of heat transferred: %1.2f kJ',Q);
+printf('\n')
+
+V2=(T2*V1)/T1;
+W=P1*100*(V2-V1);
+printf('The Value of Work done: %1.2f kJ',W);
+printf('\n')
diff --git a/2339/CH3/EX3.23.1/Ex3_23.sce b/2339/CH3/EX3.23.1/Ex3_23.sce new file mode 100755 index 000000000..7f01ed0aa --- /dev/null +++ b/2339/CH3/EX3.23.1/Ex3_23.sce @@ -0,0 +1,42 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=28;
+V1=3;
+T1=100+273;
+T2=37+273;
+G=1.4;
+Ro=8.314;
+
+v=V1/m;
+R=Ro/m;
+
+P1=(m*R*T1)/V1;
+
+printf('The Specific Volume: %1.3f m^3/kg',v);
+printf('\n')
+
+printf('The Initial Pressure: %1.2f kPa',P1);
+printf('\n')
+
+P2=(P1*T2)/T1;
+printf('The Final Pressure: %1.2f kPa',P2);
+printf('\n')
+
+Cv=(R)/(G-1);
+Cp=Cv*G;
+U=m*Cv*(T2-T1);
+H=m*Cp*(T2-T1);
+
+printf('Change in Internal Energy: %1.2f kJ',U);
+printf('\n');
+
+
+printf('Change in Heat energy: %1.2f kJ',H);
+printf('\n')
diff --git a/2339/CH3/EX3.24.1/Ex3_24.sce b/2339/CH3/EX3.24.1/Ex3_24.sce new file mode 100755 index 000000000..2e1ee7d26 --- /dev/null +++ b/2339/CH3/EX3.24.1/Ex3_24.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=3;
+V1=V1*100; //In kPa
+P1=2;
+T1=73+273;
+P2=7;
+R=0.287;
+Cv=0.718;
+Cp=1.005;
+
+m=(P1*V1)/(R*T1);
+T2=(P2*T1)/P1;
+
+U=m*Cv*(T2-T1);
+H=m*Cp*(T2-T1);
+
+printf('Change in Internal Energy: %1.2f kJ',U);
+printf('\n');
+
+printf('Change in heat Energy: %1.2f kJ',H);
+printf('\n')
diff --git a/2339/CH3/EX3.25.1/Ex3_25.sce b/2339/CH3/EX3.25.1/Ex3_25.sce new file mode 100755 index 000000000..a8a4e70a9 --- /dev/null +++ b/2339/CH3/EX3.25.1/Ex3_25.sce @@ -0,0 +1,34 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=1;
+T1=27+273;
+T2=197+273;
+V1=2.1;
+R=0.287;
+Cp=1.005;
+
+W=m*R*(T2-T1);
+printf('Work Done: %1.2f kJ',W);
+printf('\n');
+
+Q=m*Cp*(T2-T1);
+U=Q-W;
+
+printf('Change in Heat Energy: %1.2f kJ',Q);
+printf('\n');
+printf('Change in Internal Energy: %1.2f kJ',U);
+printf('\n');
+
+P=(m*R*T1)/(V1);
+V2=(V1*T2)/(T1);
+printf('Pressure: %1.2f kPa',P);
+printf('\n');
+printf('Final Volume: %1.1f m^3',V2);
+printf('\n');
diff --git a/2339/CH3/EX3.26.1/Ex3_26.sce b/2339/CH3/EX3.26.1/Ex3_26.sce new file mode 100755 index 000000000..d391ecdb7 --- /dev/null +++ b/2339/CH3/EX3.26.1/Ex3_26.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.5;
+P1=0.3;
+V2=0.1;
+P2=(P1*V1)/(V2);
+printf('Final Pressure: %1.2f bar',P2);
+printf('\n');
+
+W=(P1*100*V1)*(log(V2/V1));
+printf('Work Done: %1.2f kJ',W);
+printf('\n');
+
+printf('Change in Internal Energy: 0 kJ as it is Isothermal Process');
+printf('\n');
+
+printf('Change in Heat Energy: %1.2f kJ',W);
+printf('\n');
diff --git a/2339/CH3/EX3.27.1/Ex3_27.sce b/2339/CH3/EX3.27.1/Ex3_27.sce new file mode 100755 index 000000000..4b8841f25 --- /dev/null +++ b/2339/CH3/EX3.27.1/Ex3_27.sce @@ -0,0 +1,49 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=0.19;
+P1=1;
+V1=0.16;
+T1=17+273;
+P2=4.1;
+V2=0.046;
+Ro=8.314;
+
+R=(P1*100*V1)/(m*T1);
+printf('Gas Constant: %1.2f kJ/kg K',R);
+printf('\n');
+
+M=Ro/R;
+printf('Molecular Mass: %1.2f kg/kg mole',M);
+printf('\n');
+
+G=(log(P1/P2))/(log(V2/V1));
+printf('Ratio of Specific Heats: %1.2f ',G);
+printf('\n');
+
+Cv=(R)/(G-1);
+printf('Value of Cv: %1.2f kJ/kg K',Cv);
+printf('\n');
+
+Cp=G*Cv;
+printf('Value of Cp: %1.2f kJ/kg K',Cp);
+printf('\n');
+
+T2=(P2*100*V2)/(m*R);
+U=m*Cv*(T2-T1);
+printf('Change in Internal Energy: %1.2f kJ',U);
+printf('\n');
+
+H=m*Cp*(T2-T1);
+printf('Heat Transfer: %1.2f kJ',H);
+printf('\n');
+
+W=((P1*100*V1)-(P2*100*V2))/(G-1);
+printf('Work Done: %1.2f kJ',W);
+printf('\n');
diff --git a/2339/CH3/EX3.28.1/Ex3_28.sce b/2339/CH3/EX3.28.1/Ex3_28.sce new file mode 100755 index 000000000..2c497a32c --- /dev/null +++ b/2339/CH3/EX3.28.1/Ex3_28.sce @@ -0,0 +1,32 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.19;
+P1=5;
+T1=190+273;
+P2=1;
+H=100;
+G=1.4;
+R=0.287;
+Cp=1.005;
+
+V2=V1*((P1/P2)^(1/G));
+W=((P1*100*V1)-(P2*100*V2))/(G-1);
+m=(P1*100*V1)/(R*T1);
+T2=T1*((P2/P1)^((G-1)/G))
+x=H/(m*Cp);
+T3=x+T2;
+
+V3=(V2*T3)/T2;
+Wo=P2*100*(V3-V2);
+Wf=W+Wo;
+printf('Total Work Done: %1.2f kJ',Wf);
+printf('\n');
+
+
diff --git a/2339/CH3/EX3.29.1/Ex3_29.sce b/2339/CH3/EX3.29.1/Ex3_29.sce new file mode 100755 index 000000000..c334f4e4b --- /dev/null +++ b/2339/CH3/EX3.29.1/Ex3_29.sce @@ -0,0 +1,44 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.1;
+V3=V1;
+P1=10;
+T1=200+273;
+P2=3;
+R=0.287;
+G=1.4;
+Cv=0.718;
+
+m=(P1*100*V1)/(R*T1);
+T2=T1*((P2/P1)^((G-1)/G));
+V2=V1*((P1/P2)^((1)/G));
+T3=T2;
+P3=(P2*V2)/V3;
+printf('Pressure after Isothermal Compression: %1.2f bar',P3);
+printf('\n');
+printf('Temperature after isothermal compression: %1.2f K',T2);
+printf('\n');
+
+W1=((P1*100*V1)-(P2*100*V2))/(G-1);
+printf('Work Developed during adiabatic expansion: %2.0f kJ',W1);
+printf('\n');
+
+W2=(P2*100*V2)*log(V3/V2);
+printf('Work of Compression: %1.2f kJ',W2);
+printf('\n');
+
+Q=m*Cv*(T1-T3);
+printf('Heat supplied in 3rd Process: %1.2f kJ',Q);
+printf('\n');
+
+U=m*Cv*(T2-T1);
+printf('Change in Internal Energy: %1.2f kJ',U);
+printf('\n');
+
diff --git a/2339/CH3/EX3.3.1/Ex3_3.sce b/2339/CH3/EX3.3.1/Ex3_3.sce new file mode 100755 index 000000000..06e8dd6a3 --- /dev/null +++ b/2339/CH3/EX3.3.1/Ex3_3.sce @@ -0,0 +1,22 @@ +clear +clc + +P1=2.75*100; //Pressure in kPa +V1=0.09 //Volume in m^3 +T1=185+273; //Temperature in Kelvin +T2=15+273; +R=0.29; +Cp=1.005; +Cv=0.715; + +m=(P1*V1)/(R*T1); +V2=(T2/T1)*V1; +W=P1*(V2-V1); +printf('The Work Done: %2.3f kJ',W); +printf('\n'); +Q=m*Cp*(T2-T1); +printf('The Heat Transfer: %2.3f kJ',Q); +printf('\n'); +U=Q-W; +printf('The change in Internal Energy: %2.3f kJ',U); +printf('\n'); diff --git a/2339/CH3/EX3.30.1/Ex3_30.sce b/2339/CH3/EX3.30.1/Ex3_30.sce new file mode 100755 index 000000000..87212cfdf --- /dev/null +++ b/2339/CH3/EX3.30.1/Ex3_30.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+//Inputs
+V1=0.028;
+P1=1;
+T1=27+273;
+n=1.3;
+V2=0.0046;
+T3=T1;
+
+T2=T1*((V1/V2)^(n-1));
+printf('Temperature after compression: %1.2f K',T2);
+printf('\n');
+
+P2=P1*((V1/V2)^n);
+W=((P1*100*V1)-(P2*100*V2))/(n-1);
+printf('Work Done: %1.2f kJ',W);
+printf('\n');
+
+P3=(T3*P2)/T2;
+printf('Final Pressure: %1.2f bar',P3);
+printf('\n');
+
+
+
diff --git a/2339/CH3/EX3.31.1/Ex3_31.sce b/2339/CH3/EX3.31.1/Ex3_31.sce new file mode 100755 index 000000000..bdae96ae0 --- /dev/null +++ b/2339/CH3/EX3.31.1/Ex3_31.sce @@ -0,0 +1,36 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.15;
+P1=900;
+T1=300+273;
+T3=T1;
+V2=3*V1;
+R=0.287;
+Cp=1.005;
+G=1.4;
+n=1.5;
+Cv=0.718;
+
+//Calculations
+m=(P1*V1)/(R*T1);
+T2=(V2*T1)/V1;
+Q1=m*Cp*(T2-T1);
+printf('Heat Received: %3.2f kJ',Q1);
+printf('\n');
+
+Q2=(m*Cv)*((n-G)/(n-1))*(T3-T2);
+Q3=m*R*T3*(log(1/27));
+Qr=0-(Q2+Q3);
+printf('Heat Rejected: %3.2f kJ',Qr);
+printf('\n');
+
+Eff=(1-(Qr/Q1))*100;
+printf('Efficiency: %3.2f percent',Eff);
+printf('\n');
diff --git a/2339/CH3/EX3.32.1/Ex3_32.sce b/2339/CH3/EX3.32.1/Ex3_32.sce new file mode 100755 index 000000000..5365650c3 --- /dev/null +++ b/2339/CH3/EX3.32.1/Ex3_32.sce @@ -0,0 +1,33 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+M=27;
+P1=1;
+T1=60+273;
+n=1.3;
+Cvm=21;
+Ro=8.314;
+R=Ro/M;
+
+//Calculations
+V1=(R*T1)/(P1*100);
+V2=V1/12;
+x=V1/V2;
+P2=P1*(x^n);
+W=((P1*100*V1)-(P2*100*V2))/(n-1);
+printf('Work Done: %3.2f kJ/kg',W);
+printf('\n');
+
+Cv=Cvm/M;
+Cp=Cv+R;
+G=Cp/Cv;
+
+Q=((G-n)/(G-1))*W;
+printf('Heat Transfer during the process: %3.2f kJ/kg',Q);
+printf('\n');
diff --git a/2339/CH3/EX3.33.1/Ex3_33.sce b/2339/CH3/EX3.33.1/Ex3_33.sce new file mode 100755 index 000000000..c403e7013 --- /dev/null +++ b/2339/CH3/EX3.33.1/Ex3_33.sce @@ -0,0 +1,35 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+D=0.550;
+L=0.740;
+r=12;
+P1=100;
+T1=27+273;
+n=1.32;
+R=0.287;
+G=1.4;
+V=((22/7)/4)*D*D*L;
+V2=V/11;
+V1=V2+V;
+P2=P1*((r)^n);
+T2=T1*((r)^(n-1));
+printf('The Pressure at end of Compression: %3.3f kPa',P2);
+printf('\n');
+printf('The Temperature at end of Compression: %3.3f K',T2);
+printf('\n');
+m=(P1*V1)/(R*T1);
+printf('The Mass in the cylinder: %3.3f kg',m);
+printf('\n');
+W=((P1*V1)-(P2*V2))/(n-1);
+printf('The Work Done: %3.3f kJ',W);
+printf('\n');
+Q=((G-n)/(G-1))*W;
+printf('The Heat Transfer: %3.3f kJ',Q);
+printf('\n');
diff --git a/2339/CH3/EX3.34.1/Ex3_34.sce b/2339/CH3/EX3.34.1/Ex3_34.sce new file mode 100755 index 000000000..b279f7fa6 --- /dev/null +++ b/2339/CH3/EX3.34.1/Ex3_34.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=1;
+P1=10;
+T1=337+273;
+P2=1;
+V=6;
+R=0.287;
+G=1.4;
+x=log(P2/P1);
+y=log(1/V);
+n=x/y;
+printf('The Value of n: %3.3f ',n);
+printf('\n');
+V1=(m*R*T1)/(P1*100);
+V2=V1*6;
+W=((P1*100*V1)-(P2*100*V2))/(n-1);
+printf('The Work Done: %3.1f kJ',W);
+printf('\n');
+Q=((G-n)/(G-1))*W;
+printf('The Heat Transfer: %3.2f kJ',Q);
+printf('\n');
diff --git a/2339/CH3/EX3.37.1/Ex3_37.sce b/2339/CH3/EX3.37.1/Ex3_37.sce new file mode 100755 index 000000000..6ad4fb836 --- /dev/null +++ b/2339/CH3/EX3.37.1/Ex3_37.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+T1=430;
+T2=289.25;
+P2=100;
+P1=400;
+G=1.4;
+V1=0.2;
+R=287;
+Q=60;
+Cp=1.005;
+T2=T1*((P2/P1)^((G-1)/G));
+V2=V1*((P1/P2)^(1/G));
+m=(P1*1000*V1)/(R*T1);
+W1=(m*(R/1000)*(T1-T2))/(G-1);
+T3=(Q/(m*Cp))+T2;
+V3=(V2*T3)/T2;
+W2=P2*(V3-V2);
+W=W1+W2;
+printf('The Net Work Done: %3.3f kJ',W);
+printf('\n');
+n=((m*(R/1000)*(T1-T3))/W)+1;
+printf('The value of n: %3.2f ',n);
+printf('\n');
diff --git a/2339/CH3/EX3.38.1/Ex3_38.sce b/2339/CH3/EX3.38.1/Ex3_38.sce new file mode 100755 index 000000000..05ee24bce --- /dev/null +++ b/2339/CH3/EX3.38.1/Ex3_38.sce @@ -0,0 +1,30 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg Km=1;
+P1=6;
+V1=0.01;
+V2=0.05;
+P2=2;
+W1=(((P1+P2)/2)*100)*(V2-V1);
+printf('The Work done for first cycle: %3.1f kJ',W1);
+printf('\n');
+P3=P2;
+V3=(P1*V1)/P3;
+W2=P2*100*(V3-V2);
+printf('The Work done for second cycle: %3.1f kJ',W2);
+printf('\n');
+W3=(P3*100*V3)*(log(V1/V3));
+printf('The Work done for third cycle: %3.2f kJ',W3);
+printf('\n');
+W=W1+W2+W3;
+printf('The net Work done: %3.2f kJ',W);
+printf('\n');
+Q=W; //As process is cyclic
+printf('The Heat Transfer: %3.2f kJ',Q);
+printf('\n');
diff --git a/2339/CH3/EX3.39.1/Ex3_39.sce b/2339/CH3/EX3.39.1/Ex3_39.sce new file mode 100755 index 000000000..43ec4bdcd --- /dev/null +++ b/2339/CH3/EX3.39.1/Ex3_39.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.6;
+P1=1;
+T1=90+273;
+V2=0.18;
+P2=5;
+R=0.287;
+G=1.4;
+
+m=(P1*100*V1)/(R*T1);
+printf('The mass of Gas: %3.4f kg',m);
+printf('\n');
+n=(log(P2/P1))/(log(V1/V2));
+printf('The value of n: %3.3f ',n);
+printf('\n');
+Cv=R/(G-1);
+T2=((P2*V2)/(P1*V1))*T1;
+U=m*Cv*(T2-T1);
+printf('The change in Internal Energy: %3.3f kJ',U);
+printf('\n');
diff --git a/2339/CH3/EX3.4.1/Ex3_4.sce b/2339/CH3/EX3.4.1/Ex3_4.sce new file mode 100755 index 000000000..ce9fd918f --- /dev/null +++ b/2339/CH3/EX3.4.1/Ex3_4.sce @@ -0,0 +1,30 @@ +clc +clear + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K + +m=0.67; +P1=14; +T1=290+273; +R=287; +V1=(m*R*T1)/(14*(10^5)); +printf('The Volume: %2.3f kJ',V1); +printf('\n'); +V2=4*V1; +printf('The Final Volume: %2.3f kJ',V2); +printf('\n'); +x=V1/V2; +y=x^1.3; +P2=P1*y; +printf('The Final Pressure: %2.3f bar',P2); +printf('\n'); +x=V1/V2; +y=x^0.3; +T2=T1*y; +printf('The Final Temperature: %2.3f K',T2); +printf('\n'); diff --git a/2339/CH3/EX3.40.1/Ex3_40.sce b/2339/CH3/EX3.40.1/Ex3_40.sce new file mode 100755 index 000000000..908a28816 --- /dev/null +++ b/2339/CH3/EX3.40.1/Ex3_40.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+R=0.29;
+Cp=1.005;
+P1=2.75;
+P2=P1;
+V1=0.09;
+T1=185+273;
+T2=15+273;
+
+//Calculations
+V2=(V1*T2)/T1;
+m=(P1*100*V1)/(R*T1);
+Q=m*Cp*(T2-T1);
+printf('The Heat Transfer: %3.3f kJ',Q);
+printf('\n');
+W=P1*100*(V2-V1);
+printf('The Work done: %3.3f kJ',W);
+printf('\n');
diff --git a/2339/CH3/EX3.41.1/Ex3_41.sce b/2339/CH3/EX3.41.1/Ex3_41.sce new file mode 100755 index 000000000..af0722512 --- /dev/null +++ b/2339/CH3/EX3.41.1/Ex3_41.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+T1=25+273;
+T2=145+273;
+m=2;
+R=267;
+G=1.4;
+Cv=R/(G-1);
+printf('The value of Cv: %3.1f J/kg K',Cv);
+printf('\n');
+Cp=G*Cv;
+printf('The value of Cp: %3.1f J/kg K',Cp);
+printf('\n');
+U=m*Cv*(T2-T1)*(1/1000);
+printf('The change in Internal Energy: %3.1f kJ',U);
+printf('\n');
+H=m*Cp*(T2-T1)*(1/1000);
+printf('The Heat Transfer: %3.1f kJ',H);
+printf('\n');
diff --git a/2339/CH3/EX3.42.1/Ex3_42.sce b/2339/CH3/EX3.42.1/Ex3_42.sce new file mode 100755 index 000000000..8735be51d --- /dev/null +++ b/2339/CH3/EX3.42.1/Ex3_42.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+D=1;
+h=4;
+P1=100;
+T1=27+273;
+P2=125;
+Cp=14.307;
+Cv=10.183;
+V1=(22/7)*(1/4)*(D*D*h);
+R=Cp-Cv;
+m=(P1*V1)/(R*T1);
+T2=(P2*T1)/P1;
+Q=m*Cv*(T2-T1);
+printf('The Heat Transfer: %3.0f kJ',Q);
+printf('\n');
diff --git a/2339/CH3/EX3.43.1/Ex3_43.sce b/2339/CH3/EX3.43.1/Ex3_43.sce new file mode 100755 index 000000000..6346d8f27 --- /dev/null +++ b/2339/CH3/EX3.43.1/Ex3_43.sce @@ -0,0 +1,35 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+V1=0.15;
+V2=3*V1;
+P1=900;
+P2=P1;
+T1=300+273;
+T3=T1;
+n=1.5;
+G=1.4;
+R=0.287;
+Cp=1.005;
+m=(P1*V1)/(R*T1);
+T2=(T1*V2)/V1;
+Q1=m*Cp*(T2-T1);
+W1=m*R*(T2-T3)/(n-1);
+Q2=(G-n)*W1/(G-1);
+P3=P2*((T3/T2)^(G/(G-1)));
+Q3=m*R*T1*log(P3/P1);
+H_rec=Q1;
+printf('The Heat received: %3.1f kJ',H_rec);
+printf('\n');
+H_rej=0-(Q2+Q3);
+printf('The Heat Rejected: %3.1f kJ',H_rej);
+printf('\n');
+eff=100*(1-(H_rej/H_rec));
+printf('Efficiency: %3.2f percent',eff);
+printf('\n');
diff --git a/2339/CH3/EX3.44.1/Ex3_44.sce b/2339/CH3/EX3.44.1/Ex3_44.sce new file mode 100755 index 000000000..b8073b567 --- /dev/null +++ b/2339/CH3/EX3.44.1/Ex3_44.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg Km=1;
+V1=0.15;
+P1=1;
+V2=0.05;
+G=1.4;
+
+//Calculations
+P2=(V1*P1)/V2;
+W_it=P1*100*V1*log(P1/P2);
+printf('Work done in Isothermal process: %2.2f kJ',W_it);
+printf('\n');
+P2=P1*((V1/V2)^G);
+W_ad=((P1*100*V1)-(P2*100*V2))/(G-1);
+printf('Work done in Adiabatic process: %2.2f kJ',W_ad);
+printf('\n');
diff --git a/2339/CH3/EX3.45.1/Ex3_45.sce b/2339/CH3/EX3.45.1/Ex3_45.sce new file mode 100755 index 000000000..d4d764c2f --- /dev/null +++ b/2339/CH3/EX3.45.1/Ex3_45.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=1;
+Cp=1.005;
+P1=100;
+T1=17+273;
+T2=T1;
+P2=2500;
+printf('Final Temperature: %2.2f K',T2);
+printf('\n');
+
+V1=(260*T1)/(P1*1000);
+V2=(P1*V1)/P2;
+printf('Final Volume: %2.5f m^3',V2);
+printf('\n');
+n=P2/P1;
+printf('Compression ratio: %2.0f ',n);
+printf('\n');
+H=m*Cp*(T2-T1);
+printf('Change in Enthalpy: %2.2f kJ',H);
+printf('\n');
+W=P1*V1*log(P1/P2);
+printf('Work done: %2.2f kJ/kg',W);
+printf('\n');
diff --git a/2339/CH3/EX3.46.1/Ex3_46.sce b/2339/CH3/EX3.46.1/Ex3_46.sce new file mode 100755 index 000000000..9fa62652a --- /dev/null +++ b/2339/CH3/EX3.46.1/Ex3_46.sce @@ -0,0 +1,35 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P1=150;
+T1=17+273;
+P2=750;
+n=1.3;
+m=1;
+R=0.287;
+Cp=1.001;
+
+//Calculations
+T2=T1*((P2/P1)^((n-1)/n));
+printf('The final temperature: %2.1f K',T2);
+printf('\n');
+W=m*R*(T1-T2)/(n-1);
+printf('Work done: %2.2f kJ/kg',W);
+printf('\n');
+Cv=Cp-R;
+U=m*Cv*(T2-T1);
+printf('Change in internal energy: %2.2f kJ/kg',U);
+printf('\n');
+G=Cp/Cv;
+Q=((G-n)/(G-1))*W;
+printf('Amount of heat transfer: %2.2f kJ/kg',Q);
+printf('\n');
+H=m*Cp*(T2-T1);
+printf('Change in enthalpy: %2.2f kJ/kg',H);
+printf('\n');
diff --git a/2339/CH3/EX3.5.1/Ex3_5.sce b/2339/CH3/EX3.5.1/Ex3_5.sce new file mode 100755 index 000000000..f0e0616ba --- /dev/null +++ b/2339/CH3/EX3.5.1/Ex3_5.sce @@ -0,0 +1,25 @@ +clc +clear + + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K +P1=510; +V1=0.142; +P2=170; +V2=0.275; +H=-65; +Cv=0.718; +X=(P2*V2)-(P1*V1); +U=H-X; +printf('The Change in Internal Energy: %2.2f kJ',U); +printf('\n'); +G=H/U; +Cp=G*Cv; +R=Cp-Cv; +printf('The Value of R: %2.3f kJ/kg K',R); +printf('\n'); diff --git a/2339/CH3/EX3.6.1/Ex3_6.sce b/2339/CH3/EX3.6.1/Ex3_6.sce new file mode 100755 index 000000000..d7f7e0ede --- /dev/null +++ b/2339/CH3/EX3.6.1/Ex3_6.sce @@ -0,0 +1,26 @@ +clc +clear + +//Inputs +//The Values in the program are as follows: +//Temperature in Celcius converted to Kelvin(by adding 273) +//Pressure in bar converted to kPa (by multiplying 100) +//Volume in m^3 +//Value of R,Cp and Cv in kJ/kg K +P1=25; +T1=27+273; +P2=5; +T2=20+273; +V1=0.7; +Et=1.43; +Pn=101.325; +Tn=273; + +//Calculations +R=(Pn)/(Et*Tn); +m1=(Pn*V1)/(R*Tn); +V2=(m1*R*T1)/(P1*100); +m2=(P2*100*V2)/(R*T2); +mf=m1-m2; +printf('The mass of Oxygen used: %3.3f kg',mf); +printf('\n'); diff --git a/2339/CH3/EX3.7.1/Ex3_7.sce b/2339/CH3/EX3.7.1/Ex3_7.sce new file mode 100755 index 000000000..861165d33 --- /dev/null +++ b/2339/CH3/EX3.7.1/Ex3_7.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+n=1.3;
+m=1;
+T1=300;
+T2=200;
+W=90;
+Ro=8.3143;
+R=((n-m)*W)/((T1-T2)*m);
+M=Ro/R;
+printf('The molecular mass of gas is: %3.1f kg/kg mole',M);
+printf('\n');
diff --git a/2339/CH3/EX3.8.1/Ex3_8.sce b/2339/CH3/EX3.8.1/Ex3_8.sce new file mode 100755 index 000000000..5784e1404 --- /dev/null +++ b/2339/CH3/EX3.8.1/Ex3_8.sce @@ -0,0 +1,37 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+m=0.18;
+V1=0.15;
+T1=15+273;
+P1=100;
+V2=0.056;
+P2=400;
+R=(P1*V1)/(m*T1);
+printf('The Gas Constant: %3.3f kJ/kg K',R);
+printf('\n');
+M=8.3141/R;
+printf('The Molecular Mass of Gas: %3.2f kg/kg mole',M);
+printf('\n');
+x=log(P2/P1);
+y=log(V2/V1);
+G=-(x/y);
+Cv=R/(G-1);
+printf('The Cv: %3.2f kJ/kg K',Cv);
+printf('\n');
+Cp=Cv+R;
+printf('The Cp: %3.2f kJ/kg K',Cp);
+printf('\n');
+x=(G-1)/G;
+y=P2/P1;
+z=y^x;
+T2=T1*z;
+U=m*Cv*(T2-T1);
+printf('The change in Internal Energy: %3.2f kJ',U);
+printf('\n');
diff --git a/2339/CH3/EX3.9.1/Ex3_9.sce b/2339/CH3/EX3.9.1/Ex3_9.sce new file mode 100755 index 000000000..0fcf773a5 --- /dev/null +++ b/2339/CH3/EX3.9.1/Ex3_9.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+//Inputs
+//The Values in the program are as follows:
+//Temperature in Celcius converted to Kelvin(by adding 273)
+//Pressure in bar converted to kPa (by multiplying 100)
+//Volume in m^3
+//Value of R,Cp and Cv in kJ/kg K
+P1=100;
+V1=0.25;
+T1=100+273;
+V2=0.05;
+P2=750;
+G=1.4;
+R=0.298;
+n=(log(P2/P1))/(log(V1/V2));
+printf('The Index n: %1.2f ',n);
+printf('\n');
+T2=T1*((P2/P1)^((n-1)/n));
+Cv=R/(G-1);
+Cp=R+Cv;
+m=(P1*V1)/(R*T1);
+W=(m*R*(T1-T2))/(n-1);
+Q=((G-n)/(G-1))*W;
+printf('The Heat change: %2.2f kJ',Q);
+printf('\n');
+U=m*Cv*(T2-T1);
+printf('The change in Internal Energy: %2.2f kJ',U);
+printf('\n');
+
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');
diff --git a/2339/CH5/EX5.1.1/Ex5_1.sce b/2339/CH5/EX5.1.1/Ex5_1.sce new file mode 100755 index 000000000..c8dc96a30 --- /dev/null +++ b/2339/CH5/EX5.1.1/Ex5_1.sce @@ -0,0 +1,27 @@ +clc +clear + +Mf=1300; //in kg +Ma=13000; //in kg +P=7; //in bar +Cpw=4.187; //in kJ/kg K +CV=30000; //in kJ/kg +x=0.95; //Dryness Fraction +Tfw=40; //in C + +Hfw=Tfw*Cpw; + +//At 7 bar +Hf=697.2; //in kJ/kg +Hfg=2066.3; //in kJ/kg + +H=Hf+(x*Hfg); +Ms=Ma/Mf; + +Me=(Ms*(H-Hfw))/(2257); +printf('Equivalent evaporation: %3.2f kg/kg of coal',Me); +printf('\n'); + +Eff=100*(Ma*(H-Hfw))/(Mf*CV); +printf('Boiler Efficiency: %3.1f Percent',Eff); +printf('\n'); diff --git a/2339/CH5/EX5.10.1/Ex5_10.sce b/2339/CH5/EX5.10.1/Ex5_10.sce new file mode 100755 index 000000000..ba0d98e5a --- /dev/null +++ b/2339/CH5/EX5.10.1/Ex5_10.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+//At 10 bar pressure
+Tsat=179.9;
+Tsup=250;
+Cps=2.1; //in kJ/kg K
+
+Hg=2778.1; //in kJ/kg
+Ms=10; //in kg/kg of coal
+Hsup=Hg+(Cps*(Tsup-Tsat));
+
+Hfw=155;
+Me=(Ms*(Hsup-Hfw))/2257;
+
+FOE=Me/Ms; //Factor of Evaporation
+BP=(Me*370)/21.296;
+printf('Equivalent Evaporation: %3.1f kg/kg of coal',Me);
+printf('\n');
+printf('Boiler Power: %3.1f kW',BP);
+printf('\n');
diff --git a/2339/CH5/EX5.11.1/Ex5_11.sce b/2339/CH5/EX5.11.1/Ex5_11.sce new file mode 100755 index 000000000..864b8a225 --- /dev/null +++ b/2339/CH5/EX5.11.1/Ex5_11.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+Ma=1100; //in kg/hr
+CV=33000; //in kJ/kg
+Tfw=46; //in C
+P=10; //in bar
+x=0.9; //Dryness Fraction
+Eff=0.81; //Efficiency
+
+Hf=762.8;
+Hfg=2015.3;
+H=Hf+(x*Hfg);
+Hfw=192.6;
+
+Mf=(Ma*(H-Hfw))/(CV*Eff);
+printf('Amount of Coal Consumed per hour: %3.1f kg',Mf);
+printf('\n');
diff --git a/2339/CH5/EX5.12.1/Ex5_12.sce b/2339/CH5/EX5.12.1/Ex5_12.sce new file mode 100755 index 000000000..e1f88b177 --- /dev/null +++ b/2339/CH5/EX5.12.1/Ex5_12.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+Ms=7.3; //kg/kg of fuel
+Tfw=46; //in C
+P=10; //in bar
+FOE=1.17; //Factor of Evaporation
+Eff=0.79;
+Me=FOE*Ms;
+printf('Equivalent Evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
+
+Hfw=192.6; //in kJ/kg
+Hg=2778.1; //in kJ/kg
+Tsat=179.9; //in C
+Cps=2.1; //in kJ/kg K
+H=(2257*FOE)+Hfw;
+Tsup=((H-Hg)/Cps)+Tsat;
+printf('Temperature of Superheated Steam: %3.1f C',Tsup);
+printf('\n');
+
+CV=(Ms*(H-Hfw))/Eff;
+printf('Calorific Value: %3.1f kJ/kg',CV);
+printf('\n');
diff --git a/2339/CH5/EX5.13.1/Ex5_13.sce b/2339/CH5/EX5.13.1/Ex5_13.sce new file mode 100755 index 000000000..caa43538d --- /dev/null +++ b/2339/CH5/EX5.13.1/Ex5_13.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+Ma=18000; //in kg/hr
+P=10; //in bar
+x=0.97; //Dryness Fraction
+Tfw=40; //in C
+Mf=2050; //in kg/hr
+CV=28000; //kJ/kg
+
+//At 10 bar
+Hf1=762.8;
+Hfg1=2015.3;
+H=Hf1+(x*Hfg1);
+
+Hfw=167.6;
+
+Eff=(Ma*100*(H-Hfw))/(Mf*CV);
+printf('Boiler efficiency: %3.2f Percent',Eff);
+printf('\n');
+
+EA=((Ma/Mf)*(H-Hfw))/2257;
+printf('Equivalent Evaporation: %3.2f kg/kg of coal',EA);
+printf('\n');
diff --git a/2339/CH5/EX5.14.1/Ex5_14.sce b/2339/CH5/EX5.14.1/Ex5_14.sce new file mode 100755 index 000000000..3f4632214 --- /dev/null +++ b/2339/CH5/EX5.14.1/Ex5_14.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+Ma=18000; //in kg/hr
+P=12; //in bar
+x=0.97; //Dryness Fraction
+CV=27400; //in kJ/kg
+Mf=2050; //in kg>hr
+
+Qs=Mf*CV;
+printf('Heat Supplied per hour: %3.1f kJ/hr',Qs);
+printf('\n');
+
+//At 12 bar
+Hf=798.6; //in kJ/kg
+Hfg=1986.2; //in kJ/kg
+H1=Hf+(x*Hfg);
+
+//At 105 C
+Hfw=438.9; //in kJ/kg
+Eff=(Ma*100*(H1-Hfw))/Qs;
+printf('Thermal Efficiency: %3.2f Percent',Eff);
+printf('\n');
+
+Ms=Ma/Mf;
+printf('Factor of Evaporation: %3.2f ',Ms);
+printf('\n');
diff --git a/2339/CH5/EX5.15.1/Ex5_15.sce b/2339/CH5/EX5.15.1/Ex5_15.sce new file mode 100755 index 000000000..f6fa12ce2 --- /dev/null +++ b/2339/CH5/EX5.15.1/Ex5_15.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+Ms=7.5; //kg/kg of coal
+P=11; //in bar
+Tf=70; //in C
+Eff=0.75; //Efficiency
+FOE=1.15; //Factor of Evaporation
+Cps=2.1; //in kJ/kg K
+Hfw=293; //in kJ/kg
+H=(FOE*2257)+Hfw;
+
+//At 11 bar
+Hg=2781.7; //in kJ/kg
+Tsat=184.1; //in C
+Tsup=((H-Hg)/Cps)+Tsat;
+DOS=Tsup-Tsat; //Degree of Superheat
+printf('Degree of Superheat: %3.1f C',DOS);
+printf('\n');
+
+Me=(Ms*(H-Hfw))/2257;
+printf('Equivalent evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
+
+CV=(Ms*(H-Hfw))/Eff;
+printf('Calorific value of Boiler: %3.2f kJ/kg ',CV);
+printf('\n');
diff --git a/2339/CH5/EX5.16.1/Ex5_16.sce b/2339/CH5/EX5.16.1/Ex5_16.sce new file mode 100755 index 000000000..ccc45980c --- /dev/null +++ b/2339/CH5/EX5.16.1/Ex5_16.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+Ma=17000; //in kg/hr
+P=14; //in bar
+x=0.95; //Dryness Fraction
+Tfw=102; //in C
+Mf=2050; //in kg/hr
+CV=27400; //Calorific Value
+HS=Mf*CV;
+printf('Heat Supplied per hour: %3.2f kJ',HS);
+printf('\n');
+
+Hf=830.3; //in kJ/kg
+Hfg=1959.7; //in kJ/kg
+Hfw=427.5; //in kJ/kg
+
+H=Hf+(x*Hfg);
+Eff=(Ma*100*(H-Hfw))/(Mf*CV);
+printf('Efficiency of Boiler: %3.2f Percent',Eff);
+printf('\n');
+
+Ms=Ma/Mf;
+Me=(Ms*(H-Hfw))/2257;
+printf('Equivalent evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
diff --git a/2339/CH5/EX5.17.1/Ex5_17.sce b/2339/CH5/EX5.17.1/Ex5_17.sce new file mode 100755 index 000000000..0c55d621e --- /dev/null +++ b/2339/CH5/EX5.17.1/Ex5_17.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+Ma=1800; //kg/hr
+P=12; //in bar
+x=0.97; //Dryness Fraction
+Tfw=105; //in C
+Mf=2050; //in kg/hr
+CV=27400; //in kJ/kg
+
+Q=Mf*CV;
+printf('Heat Supplied: %3.2f kJ',Q);
+printf('\n');
+
+//At 12 bar pressure
+Hf=798.6; //in kJ/kg
+Hfg=1986.2; //in kJ/kg
+H=Hf+(x*Hfg);
+Hfw=4.187*Tfw;
+
+Me=(Ma*(H-Hfw))/(2257*Mf);
+printf('Equivalent Evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
+
+Eff=(Ma*100*(H-Hfw))/(CV*Mf);
+printf('Efficiency of boiler: %3.2f Percent',Eff);
+printf('\n');
diff --git a/2339/CH5/EX5.18.1/Ex5_18.sce b/2339/CH5/EX5.18.1/Ex5_18.sce new file mode 100755 index 000000000..d08a6fc0d --- /dev/null +++ b/2339/CH5/EX5.18.1/Ex5_18.sce @@ -0,0 +1,11 @@ +clc
+clear
+
+Me=10; //kg/kg
+CV=34000 //kJ/kg
+
+x=Me*2257;
+
+Eff=100*x/CV;
+printf('Efficiency of Boiler: %3.2f Percent',Eff);
+printf('\n');
diff --git a/2339/CH5/EX5.19.1/Ex5_19.sce b/2339/CH5/EX5.19.1/Ex5_19.sce new file mode 100755 index 000000000..c8e4a4461 --- /dev/null +++ b/2339/CH5/EX5.19.1/Ex5_19.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+Ma=5500; //kg/hr
+P=1; //bar
+x=0.94; //Dryness Fraction
+Tfw=40; //in C
+Mf=600; //kg/hr
+CV=32000; //kJ/kg
+Hfw=Tfw*4.187;
+
+//At 1 bar pressure
+Hf=417.5; //kJ/kg
+Hfg=2258; //kJ/kg
+H=Hf+(x*Hfg);
+Ms=Ma/Mf;
+
+Me=(Ms*(H-Hfw))/2257;
+printf('Equivalent Evaporation: %3.3f kg/kg of coal',Me);
+printf('\n');
+
+Eff=(Ms*100*(H-Hfw))/CV;
+printf('Efficiency: %3.2f percent',Eff);
+printf('\n');
diff --git a/2339/CH5/EX5.2.1/Ex5_2.sce b/2339/CH5/EX5.2.1/Ex5_2.sce new file mode 100755 index 000000000..451bf7068 --- /dev/null +++ b/2339/CH5/EX5.2.1/Ex5_2.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+Ma=5400; //in kg/hr
+Tfw=42; //in C
+P=7.6; //in bar
+Mf=670; //in kg/hr
+x=0.98; //Dryness Fraction
+CV=31000; //kJ/kg
+Ms=Ma/Mf;
+Hf=175.81; //in kJ/kg
+Hfw=Hf;
+
+//Now at 7.6 bar pressure
+Hf=711.8; //in kJ/kg
+Hfg=2055.2; //in kJ/kg
+
+H=Hf+(x*Hfg);
+Eff=100*(Ma*(H-Hfw))/(Mf*CV);
+printf('Boiler Efficiency %3.1f percent',Eff);
+printf('\n');
+
+Me=(Ms*(H-Hfw))/(2257);
+printf('Equivalent evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
diff --git a/2339/CH5/EX5.3.1/Ex5_3.sce b/2339/CH5/EX5.3.1/Ex5_3.sce new file mode 100755 index 000000000..75e3cc2f6 --- /dev/null +++ b/2339/CH5/EX5.3.1/Ex5_3.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+P=12; //in bar
+CV=34000; //in kJ/kg
+T=250; //in C
+Ms=10; //in kg/kg of coal
+Tfw=36; //in C
+Hfw=150.74; //in kJ/kg
+Hg=2784.8; //in kJ/kg
+Tsup=T;
+Tsat=188; //in C
+Cps=2.1; //in kJ/kg K
+H=Hg+(Cps*(Tsup-Tsat));
+
+Me=(Ms*(H-Hfw))/2257;
+printf('Equivalent evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
+
+Eff=(Me*250)/21.296;
+printf('Boiler Power: %3.2f kW',Eff);
+printf('\n');
diff --git a/2339/CH5/EX5.4.1/Ex5_4.sce b/2339/CH5/EX5.4.1/Ex5_4.sce new file mode 100755 index 000000000..6dedfc8b2 --- /dev/null +++ b/2339/CH5/EX5.4.1/Ex5_4.sce @@ -0,0 +1,30 @@ +clc
+clear
+
+Ma=35500; //kg of steam
+Mf=3460;
+CV=39500;
+Ms=Ma/Mf;
+
+
+Hfw2=313.9; //in kJ/kg
+Hfw1=71.4; //in kJ/kg
+
+Q=Ma*(Hfw2-Hfw1); //Heat added in economizer
+H=2915.0; //in kJ/kg
+
+Me=(Ms*(H-Hfw2))/2257;
+printf('Equivalent evaporation: %3.2f kg/kg of Oil',Me);
+printf('\n');
+
+Eff1=(Ma*100*(H-Hfw2))/(Mf*CV);
+printf('Thermal Efficiency of boiler: %3.1f Percent',Eff1);
+printf('\n');
+
+Eff2=(Ma*100*(H-Hfw1))/(Mf*CV);
+printf('Thermal Efficiency of Boiler plant: %3.1f Percent',Eff2);
+printf('\n');
+
+HU=860875000/(Mf*CV);
+printf('Heat Utilized by Economizer: %3.1f Percent',HU);
+printf('\n');
diff --git a/2339/CH5/EX5.5.1/Ex5_5.sce b/2339/CH5/EX5.5.1/Ex5_5.sce new file mode 100755 index 000000000..8c34472d8 --- /dev/null +++ b/2339/CH5/EX5.5.1/Ex5_5.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+Ma=10000; //in kg/hr
+P=7; //in bar
+
+Tfw=40; //in C
+Hfw=167.6; //in kJ/kg
+H=2763.5; //in kJ/kg
+
+Q=Ma*(H-Hfw)/60; //Heat per minute
+SA=Q/2720; //Heating surface area required
+printf('Heating surface area required: %3.1f m^2',SA);
+printf('\n');
+
+GA=SA/25;
+printf('Grate area required: %3.1f m^2',GA);
+printf('\n');
diff --git a/2339/CH5/EX5.6.1/Ex5_6.sce b/2339/CH5/EX5.6.1/Ex5_6.sce new file mode 100755 index 000000000..d0e4643c9 --- /dev/null +++ b/2339/CH5/EX5.6.1/Ex5_6.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+Ma=2400; //in kg
+Mf=240; //in kg
+P=12; //in bar
+CV=33500; //in kJ/kg
+Tfw=120; //in C
+Cpw=4.187;
+Hfw=Cpw*Tfw;
+H=2784.8; //in kJ/kg
+Mfa=Mf-(0.1*Mf);
+Eff=(Ma*100*(H-Hfw))/(Mfa*CV);
+printf('Thermal Efficiency: %3.1f percent',Eff);
+printf('\n');
+
+Eff1=(Ma*100*(H-Hfw))/(Mf*CV);
+printf('Thermal Efficiency of boiler and grate: %3.1f percent',Eff1);
+printf('\n');
diff --git a/2339/CH5/EX5.7.1/Ex5_7.sce b/2339/CH5/EX5.7.1/Ex5_7.sce new file mode 100755 index 000000000..b6f28b908 --- /dev/null +++ b/2339/CH5/EX5.7.1/Ex5_7.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+Mf=255; //in kg
+x=0.94; //Dryness Fraction
+CV=30100; //in kJ/kg
+P=11.5; //in bar
+Ma=2100; //in kg
+Tfw=25; //in C
+Ms=Ma/Mf;
+
+Hfw=104.9; //in kJ/kg
+Hf=790.1; //in kJ/kg
+Hfg=1993.2; //in kJ/kg
+H=Hf+(x*Hfg);
+
+Me=(Ms*(H-Hfw))/2257;
+Eff=(Ma*100*(H-Hfw))/(Mf*CV);
+printf('Equivalent Evaporation: %3.2f kg/kg of coal \n',Me)
+printf('Thermal Efficiency: %3.1f percent',Eff);
+printf('\n');
diff --git a/2339/CH5/EX5.8.1/Ex5_8.sce b/2339/CH5/EX5.8.1/Ex5_8.sce new file mode 100755 index 000000000..d9bdeed6e --- /dev/null +++ b/2339/CH5/EX5.8.1/Ex5_8.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+Hf=762.8; //in kJ/kg
+Hfg=2015.3; //in kJ/kg
+x=0.95; //Dryness Fraction
+Ma=1000;
+Eff=0.75;
+CV=31000;
+
+H=Hf+(x*Hfg);
+Cpw=4.187;
+T=50;
+
+Hfw=Cpw*T;
+Q=Ma*(H-Hfw);
+
+Mf=Q/(Eff*CV);
+y=Mf/0.9;
+
+Eff1=(Q*100)/(y*CV);
+printf('Efficiency of Boiler and grate: %3.1f percent',Eff1);
+printf('\n');
diff --git a/2339/CH5/EX5.9.1/Ex5_9.sce b/2339/CH5/EX5.9.1/Ex5_9.sce new file mode 100755 index 000000000..ca466b32d --- /dev/null +++ b/2339/CH5/EX5.9.1/Ex5_9.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+//At 10 bar
+Hg=2778.1; //in kJ/kg
+Cp=2.1; //in kJ/kg K
+T=50;
+CV=30000; //in kJ/kg
+
+H=Hg+(Cp*T);
+C=4.187;
+Tf=30;
+Hfw=C*Tf;
+
+Ms=800/100;
+
+Me=(Ms*(H-Hfw))/2257;
+printf('Equivalent Evaporation: %3.2f kg/kg of coal',Me);
+printf('\n');
+
+Eff=(Ms*100*(H-Hfw))/CV;
+printf('Efficiency of Boiler and grate: %3.1f percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.1.1/Ex6_1.sce b/2339/CH6/EX6.1.1/Ex6_1.sce new file mode 100755 index 000000000..91f770d90 --- /dev/null +++ b/2339/CH6/EX6.1.1/Ex6_1.sce @@ -0,0 +1,35 @@ +clc +clear + +m=1; //in kg +R=0.287; //Universal Gas Constant +r=7; //Compression Ratio +P1=1; //in bar +T1=24+273; //in K +T3=2000; //in K +G=1.4; //Gamma + +ASE=(1-(1/(r)^(G-1)))*100; +printf('Air Standard Efficiency is %3.1f Percent',ASE); +printf('\n'); + +P2=P1*(r)^G; +printf('Pressure at end of Compression is %3.2f Bar',P2); +printf('\n'); + +T2=T1*((r)^(G-1)); +printf('Temperature at end of Compression is %3.2f K',T2); +printf('\n'); + +Cv=0.718; +Q=Cv*(T3-T2); +printf('Heat Supplied is %3.2f kJ/kg',Q); +printf('\n'); + +W=ASE*Q/100; +V1=(m*R*T1)/(P1*100); +V2=V1/r; +V=V1-V2; +Pm=W/V; +printf('Mean Effective Pressure is %3.2f kPa',Pm); +printf('\n'); diff --git a/2339/CH6/EX6.10.1/Ex6_10.sce b/2339/CH6/EX6.10.1/Ex6_10.sce new file mode 100755 index 000000000..1edc9d786 --- /dev/null +++ b/2339/CH6/EX6.10.1/Ex6_10.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+T3=1350+273; //in K
+T1=30+273; //in K
+Qs=750; //in kJ/kg
+Cv=0.718;
+G=1.4;
+
+//For Process 2-3
+T2=T3-(Qs/Cv);
+r=(T2/T1)^(1/(G-1));
+printf('The compression Ratio is %3.2f ',r);
+printf('\n');
+
+Eff=100*(1-(1/(r^(G-1))));
+printf('The Efficiency is %3.1f Percent',Eff);
+printf('\n');
+
+W=Eff*Qs/100;
+printf('The Work Output is %3.0f kJ/kg',W);
+printf('\n');
+
+P21=(r^G);
+P32=T3/T2;
+P31=P21*P32;
+printf('Ratio of maximum to minimum pressure is %3.2f ',P31);
+printf('\n');
diff --git a/2339/CH6/EX6.11.1/Ex6_11.sce b/2339/CH6/EX6.11.1/Ex6_11.sce new file mode 100755 index 000000000..7f8d640ad --- /dev/null +++ b/2339/CH6/EX6.11.1/Ex6_11.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+Vs=500; //in cm^3
+Vc=55; //in cm^3
+T1=30+273; //in K
+P1=1; //in bar
+T3=1450+273; //in K
+G=1.4;
+R=0.287;
+Cv=0.718;
+r=(Vs+Vc)/Vc;
+Eff=100*(1-(1/(r^(G-1))));
+printf('The standard air Efficiency is %3.2f Percent',Eff);
+printf('\n');
+
+T2=T1*(r^(G-1));
+Qs=Cv*(T3-T2);
+W=Eff*Qs;
+
+V1=Vc+Vs;
+m=(P1*100*V1*(10^-6))/(R*T1);
+Pm=(W*m)/(100*(Vs*(10^-6)));
+printf('The Mean Effective Pressure is %3.1f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.12.1/Ex6_12.sce b/2339/CH6/EX6.12.1/Ex6_12.sce new file mode 100755 index 000000000..805bdc838 --- /dev/null +++ b/2339/CH6/EX6.12.1/Ex6_12.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+r=6; //Compression Ratio
+T1=20+273; //in K
+G=1.4;
+Cv=0.718;
+Qs=1900;
+Eff=100*(1-(1/(r^(G-1))));
+printf('The Efficiency is %3.2f Percent',Eff);
+printf('\n');
+
+T2=T1*(r^(G-1));
+printf('The value of T2 is %3.0f K',T2);
+printf('\n');
+
+T3=(Qs/Cv)+T2;
+printf('The value of T3 is %3.0f K',T3);
+printf('\n');
+
+T4=T3/(r^(G-1));
+printf('The value of T4 is %3.0f K',T4);
+printf('\n');
+
+W=Qs*Eff/100;
+printf('The Work Output is %3.0f kJ/kg',W);
+printf('\n');
diff --git a/2339/CH6/EX6.13.1/Ex6_13.sce b/2339/CH6/EX6.13.1/Ex6_13.sce new file mode 100755 index 000000000..bd69f8318 --- /dev/null +++ b/2339/CH6/EX6.13.1/Ex6_13.sce @@ -0,0 +1,56 @@ +clc
+clear
+
+D=0.1; //Diameter
+L=0.12; //Length
+V=(22/7)*(1/4)*D*D*L;
+T1=19+273;
+r=6.5; //Compression ratio
+P1=1; //in bar
+G=1.4; //Gamma
+Vs=9.425*(10^-4);
+Vc=Vs/(r-1);
+V1=Vc+Vs;
+printf('V1= %3.5f m^3',V1);
+printf('\n');
+
+V2=Vc;
+V3=Vc;
+P2=P1*(r^G);
+printf('P2= %3.1f bar',P2);
+printf('\n');
+
+T2=T1*(r^(G-1));
+printf('T2= %3.1f K',T2);
+printf('\n');
+
+//For process 2-3
+Qs=1900;
+Cv=0.718;
+T3=(Qs/Cv)+T2;
+printf('T3= %3.1f K',T3);
+printf('\n');
+
+P3=P2*(T3/T2);
+printf('P3= %3.1f bar',P3);
+printf('\n');
+
+//For process 4-1
+V4=V1;
+P4=P3*((V3/V4)^G);
+printf('P4= %3.1f bar',P4);
+printf('\n');
+
+T4=T1*(P4/P1);
+printf('T4= %3.1f K',T4);
+printf('\n');
+
+Eff=100*(1-(1/(r^(G-1))));
+printf('Efficiency= %3.1f Percent',Eff);
+printf('\n');
+
+R=0.287;
+m=(P1*100*V1)/(R*T1);
+Pm=(Eff*Qs*m)/(10000*Vs);
+printf('Mean Effective Pressure= %3.1f bar',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.14.1/Ex6_14.sce b/2339/CH6/EX6.14.1/Ex6_14.sce new file mode 100755 index 000000000..d47470e85 --- /dev/null +++ b/2339/CH6/EX6.14.1/Ex6_14.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+P1=1; //in bar
+T1=20+273; //in K
+P2=39; //in bar
+P3=P2;
+T3=1100+273; //in K
+G=1.4;
+
+//For reversible Adiabatic Process 1-2
+T2=T1*((P2/P1)^((G-1)/G));
+r=(P2/P1)^(1/G);
+
+Z=T3/T2;
+
+Eff=100*(1-((1/(r^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency: %2.2f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.15.1/Ex6_15.sce b/2339/CH6/EX6.15.1/Ex6_15.sce new file mode 100755 index 000000000..c772e30de --- /dev/null +++ b/2339/CH6/EX6.15.1/Ex6_15.sce @@ -0,0 +1,62 @@ +clc
+clear
+
+r=16; //Compression Ratio
+D=0.21; //Diameter
+L=0.3; //Length
+P1=1; //in bar
+G=1.4;
+T1=17+273; //in K
+Z=[0.1*(r-1)+1];
+Vs=(22/7)*(1/4)*D*D*L;
+Vc=Vs/15;
+V2=Vc;
+printf('Vc= V2= %2.6f m^3',Vc);
+printf('\n');
+V1=Vc+Vs;
+printf('V1= %2.4f m^3',V1);
+printf('\n');
+V3=(0.1*(Vs))+(V2);
+printf('V3= %2.4f m^3',V3);
+printf('\n');
+
+//For Process 1-2
+P2=r^G;
+printf('P2= %2.1f bar',P2);
+printf('\n');
+
+T2=T1*(r^(G-1));
+printf('T2= %2.1f K',T2);
+printf('\n');
+
+T3=Z*T2;
+printf('T3= %2.1f K',T3);
+printf('\n');
+
+P3=P2;
+P4=P3*((V3/V1)^G);
+printf('P4= %2.1f bar',P4);
+printf('\n');
+
+T4=T3*((V3/V1)^(G-1));
+printf('T4= %2.1f K',T4);
+printf('\n');
+
+Cv=0.718;
+Cp=1.005;
+
+Eff=100*(1-((Cv*(T4-T1))/(Cp*(T3-T2))));
+printf('Efficiency: %2.1f Percent',Eff);
+printf('\n');
+
+R=0.287;
+m=(P1*100*V1)/(R*T1);
+Pm=(m*((Cp*(T3-T2))-(Cv*(T4-T1))))/(Vs);
+printf('Mean Effective Pressure= %2.1f kPa',Pm);
+printf('\n');
+
+N=300; //Cycles per minute
+W=10.41;
+EP=W*(N/60);
+printf('Engine Power= %2.2f kW',EP);
+printf('\n');
diff --git a/2339/CH6/EX6.16.1/Ex6_16.sce b/2339/CH6/EX6.16.1/Ex6_16.sce new file mode 100755 index 000000000..239e1c995 --- /dev/null +++ b/2339/CH6/EX6.16.1/Ex6_16.sce @@ -0,0 +1,66 @@ +clc
+clear
+
+r=19; //Compression Ratio
+P1=1; //in bar
+T1=17+273; //in K
+Qs=730; //in kJ/cycle
+G=1.4;
+
+//For process 1-2
+m=1; //in kg
+R=0.287; //Universal Gas Constant
+V1=(m*R*T1)/(P1*100);
+printf('V1= %2.4f m^3/kg',V1);
+printf('\n');
+
+V2=V1/r;
+printf('V2= %2.4f m^3/kg',V2);
+printf('\n');
+
+P2=P1*(r^G);
+printf('P2= %2.1f bar',P2);
+printf('\n');
+
+T2=T1*(r^(G-1));
+printf('T2= %2.1f K',T2);
+printf('\n');
+
+//For Process 2-3
+Cv=0.718;
+T3=(Qs/(Cv*m))+T2;
+printf('T3= %2.1f K',T3);
+printf('\n');
+
+P3=P2;
+printf('P3= %2.1f bar',P3);
+printf('\n');
+
+//As pressure is constant
+V3=(T3/T2)*V2;
+printf('V3= %2.4f m^3/kg',V3);
+printf('\n');
+
+//For process 3-4
+V4=V1;
+T4=T3*((V3/V4)^(G-1));
+printf('T4= %2.1f K',T4);
+printf('\n');
+
+P4=P3*((V3/V4)^G);
+printf('P4= %2.2f bar',P4);
+printf('\n');
+
+Cp=1.005;
+
+W=((Cp)*(T3-T2))-((Cv*(T4-T1)));
+printf('Work Done= %2.1f kJ/kg',W);
+printf('\n');
+
+Eff=100*(W/(Cp*(T3-T2)));
+printf('Efficiency= %2.2f Percent',Eff);
+printf('\n');
+
+Pm=W/(V1-V2);
+printf('Mean Effective Pressure= %2.2f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.17.1/Ex6_17.sce b/2339/CH6/EX6.17.1/Ex6_17.sce new file mode 100755 index 000000000..acc7b5af3 --- /dev/null +++ b/2339/CH6/EX6.17.1/Ex6_17.sce @@ -0,0 +1,10 @@ +clc
+clear
+
+r=19; //Compression Ratio
+Re=9.1; //Expansion Ratio
+Z=r/Re;
+G=1.4;
+Eff=100*(1-((1/(r^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency: %2.2f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.18.1/Ex6_18.sce b/2339/CH6/EX6.18.1/Ex6_18.sce new file mode 100755 index 000000000..8450c59a6 --- /dev/null +++ b/2339/CH6/EX6.18.1/Ex6_18.sce @@ -0,0 +1,37 @@ +clc
+clear
+
+D=16; //in cm
+L=24; //in cm
+Vc=340;
+V2=Vc;
+G=1.4;
+
+Vs=(22/7)*(1/4)*D*D*L;
+V1=Vs+Vc;
+r=V1/V2;
+
+//Cut-off is 6% of the stroke
+Co1=0.06;
+
+V3=(Co1*(V1-V2))+V2;
+Z=V3/V2;
+x=(Z^G)-1;
+y=(r^(G-1))*(G)*(Z-1);
+Eff1=100*(1-((x)/(y)));
+
+
+
+//Cut-off is 10% of the stroke
+Co2=0.10;
+
+V3=(Co2*(V1-V2))+V2;
+Z=V3/V2;
+x=(Z^G)-1;
+y=(r^(G-1))*(G)*(Z-1);
+Eff2=100*(1-((x)/(y)));
+
+Loss=((Eff1-Eff2)*100)/Eff1;
+
+printf('Loss: %2.2f Percent',r);
+printf('\n');
diff --git a/2339/CH6/EX6.2.1/Ex6_2.sce b/2339/CH6/EX6.2.1/Ex6_2.sce new file mode 100755 index 000000000..59e7cc8f9 --- /dev/null +++ b/2339/CH6/EX6.2.1/Ex6_2.sce @@ -0,0 +1,15 @@ +clc
+clear
+
+T1=323; //in K
+T2=673; //in K
+G=1.4;
+
+r_G=T2/T1;
+r=(r_G)^(1/(G-1));
+printf('Compression Ratio is %2.2f ',r);
+printf('\n');
+
+Eff=100*(1-(1/(r^(G-1))));
+printf('Air Standard Efficiency is %2.0f Percent ',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.20.1/Ex6_20.sce b/2339/CH6/EX6.20.1/Ex6_20.sce new file mode 100755 index 000000000..182d4b26d --- /dev/null +++ b/2339/CH6/EX6.20.1/Ex6_20.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+T3=1000+273; //in K
+T1=27+273; //in K
+G=1.25;
+
+r=(T3/T1)^G;
+printf('Compression Ratio: %2.1f ',r);
+printf('\n');
+
+T2=sqrt(T1*T3);
+T4=T2;
+printf('T2=T4= %2.0f K',T2);
+printf('\n');
+
+Cv=0.718;
+W=Cv*[(sqrt(T3))-(sqrt(T1))]^2;
+printf('Maximum Work Done: %2.0f kJ/kg',W);
+printf('\n');
diff --git a/2339/CH6/EX6.21.1/Ex6_21.sce b/2339/CH6/EX6.21.1/Ex6_21.sce new file mode 100755 index 000000000..32f83276e --- /dev/null +++ b/2339/CH6/EX6.21.1/Ex6_21.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+r=6; //Compression Ratio
+G=1.4;
+
+Eff=100*(1-(1/(r^(G-1))));
+printf('Efficiency: %2.2f Percent',Eff);
+printf('\n');
+
+m=1; //in kg
+R=0.287; //Universal Gas Constant
+T1=27+273; //in K
+P1=1; //in bar
+
+V1=(m*R*T1)/(P1*100);
+V2=V1/r;
+Vc=V2;
+Vs=V1-Vc;
+
+T2=T1*(r^(G-1));
+Cv=0.718;
+Qs=1046;
+T3=(Qs/Cv)+T2;
+T4=T3/(r^(G-1));
+W=Qs-(Cv*(T4-T1));
+Pm=W/Vs;
+printf('Effective Mean Pressure: %2.2f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.22.1/Ex6_22.sce b/2339/CH6/EX6.22.1/Ex6_22.sce new file mode 100755 index 000000000..33adb2003 --- /dev/null +++ b/2339/CH6/EX6.22.1/Ex6_22.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+T1=87+273; //in K
+r=14; //Compression Ratio
+T3=1795+273; //in K
+T4=677+273; //in K
+G=1.4;
+T2=T1*(r^(G-1));
+printf('T2= %2.1f K',T2);
+printf('\n');
+
+Cp=1.005;
+Cv=0.718;
+W=[Cp*(T3-T2)]-[Cv*(T4-T1)];
+Qs=Cp*(T3-T2);
+Eff=(W*100)/Qs;
+printf('Efficiency: %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.23.1/Ex6_23.sce b/2339/CH6/EX6.23.1/Ex6_23.sce new file mode 100755 index 000000000..8fcea1f9b --- /dev/null +++ b/2339/CH6/EX6.23.1/Ex6_23.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+r=16; //Compression Ratio
+P1=1; //in bar
+T1=20+273; //in K
+T3=1431+273; //in K
+G=1.4;
+T2=T1*[r^(G-1)];
+m=1;
+R=0.287;
+V1=(m*R*T1)/(P1*100);
+V2=V1/r;
+
+//For Constant Pressure Process 2-3
+V3=V2*(T3/T2);
+Z=V3/V2;
+Eff=100*(1-((1/(r^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+Cp=1.005;
+Qs=Cp*(T3-T2);
+W=Qs*Eff/100;
+Vs=V1-V2;
+Pm=W/Vs;
+printf('Effective Mean Pressure %2.1f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.24.1/Ex6_24.sce b/2339/CH6/EX6.24.1/Ex6_24.sce new file mode 100755 index 000000000..f4315b50a --- /dev/null +++ b/2339/CH6/EX6.24.1/Ex6_24.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+r=8;
+T1=310; //in K
+T3=1600; //in K
+G=1.4;
+Cv=0.717;
+
+//For process 1-2
+T2=T1*(r^(G-1));
+
+//Now Heat Supplied
+Qs=Cv*(T3-T2);
+printf('Heat Supplied= %2.1f kJ/kg',Qs);
+printf('\n');
+
+//Efficiency of Cycle
+Eff=100*[1-(1/(r^(G-1)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.25.1/Ex6_25.sce b/2339/CH6/EX6.25.1/Ex6_25.sce new file mode 100755 index 000000000..c1446301e --- /dev/null +++ b/2339/CH6/EX6.25.1/Ex6_25.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+r=15; //Compression Ratio
+P1=100; //in kPa
+T1=27+273;
+Cp=1.006;
+Cv=0.717;
+G=1.4;
+
+//Cut off takes place at 12% of Working Stroke
+T2=T1*(r^(G-1));
+printf('T2= %2.1f K',T2);
+printf('\n');
+
+P2=P1*(r^G);
+printf('P2= %2.1f kPa',P2);
+printf('\n');
+
+Z=(0.12*(r-1))+1;
+Eff=100*(1-((1/(r^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.26.1/Ex6_26.sce b/2339/CH6/EX6.26.1/Ex6_26.sce new file mode 100755 index 000000000..323713a9a --- /dev/null +++ b/2339/CH6/EX6.26.1/Ex6_26.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+T1=288; //in K
+T3=1673; //in K
+Qs=800; //in kJ/kg
+G=1.4;
+Cv=0.718;
+R=0.287;
+P1=1;
+
+Cp=Cv*G;
+T2=T3-(Qs/Cp);
+
+x=T2/T1;
+r=x^(1/(G-1));
+printf('Compression Ratio %2.1f ',r);
+printf('\n');
+
+Eff=100*[1-(1/(r^(G-1)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+P3=r*T3*P1/T1;
+printf('P3= %2.1f bar',P3);
+printf('\n');
diff --git a/2339/CH6/EX6.27.1/Ex6_27.sce b/2339/CH6/EX6.27.1/Ex6_27.sce new file mode 100755 index 000000000..9566f58fb --- /dev/null +++ b/2339/CH6/EX6.27.1/Ex6_27.sce @@ -0,0 +1,8 @@ +clc
+clear
+
+T2=293; //in K
+Eff=0.7;
+T1=T2/(1-Eff);
+printf('T1= %2.1f K',T1);
+printf('\n');
diff --git a/2339/CH6/EX6.28.1/Ex6_28.sce b/2339/CH6/EX6.28.1/Ex6_28.sce new file mode 100755 index 000000000..10d3b1c92 --- /dev/null +++ b/2339/CH6/EX6.28.1/Ex6_28.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+T1=330; //in K
+T2=876; //in K
+T3=2223; //in K
+T4=1143; //in K
+P1=1; //in bar
+G=1.4;
+
+Cv=0.718;
+Cp=1.005;
+Eff=100*[1-((Cv*(T4-T1))/(Cp*(T3-T2)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+//For Process 1-2
+P2=P1*[(T2/T1)^(G/(G-1))];
+printf('Maximum Pressure %2.1f bar',P2);
+printf('\n');
diff --git a/2339/CH6/EX6.29.1/Ex6_29.sce b/2339/CH6/EX6.29.1/Ex6_29.sce new file mode 100755 index 000000000..a60d539e4 --- /dev/null +++ b/2339/CH6/EX6.29.1/Ex6_29.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+T1=25+273; //in K
+T3=1500+273; //in K
+Qa=900; //in kJ/kg
+Cv=0.718;
+G=1.4;
+
+T2=T3-(Qa/Cv);
+r=(T2/T1)^(1/(G-1));
+printf('Compression Ratio is %2.1f ',r);
+printf('\n');
+
+Eff=100*[1-(1/(r^(G-1)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+Px=r^G; //Max Pressure
+Py=T3/T2; //1/Min Pressure
+P=Px*Py;
+printf('Pressure Ratio %2.1f ',P);
+printf('\n');
diff --git a/2339/CH6/EX6.3.1/Ex6_3.sce b/2339/CH6/EX6.3.1/Ex6_3.sce new file mode 100755 index 000000000..ed5c5fa27 --- /dev/null +++ b/2339/CH6/EX6.3.1/Ex6_3.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+P1=97; //in kPa
+T1=323; //in K
+r=5; //Compression Ratio
+Q=930; //in kJ/kg
+G=1.4;
+Cv=0.718;
+T2=T1*(r^(G-1));
+T3=(Q/Cv)+T2;
+printf('Maximum Temperature Attained is %2.2f K',T3);
+printf('\n');
+
+Eff=100*(1-(1/(r)^(G-1)));
+printf('Thermal Efficiency of cycle is %2.1f Percent',Eff);
+printf('\n');
+
+W=Eff*Q/100;
+printf('Work Done is %2.2f kJ/kg',W);
+printf('\n');
diff --git a/2339/CH6/EX6.30.1/Ex6_30.sce b/2339/CH6/EX6.30.1/Ex6_30.sce new file mode 100755 index 000000000..97078eed0 --- /dev/null +++ b/2339/CH6/EX6.30.1/Ex6_30.sce @@ -0,0 +1,29 @@ +clc
+clear
+
+P1=1; //in bar
+T1=15+273; //in K
+P2=15; //in bar
+P3=40; //in bar
+G=1.4;
+Cv=0.718;
+
+r=(P2/P1)^(1/G);
+printf('Compression Ratio is %2.1f ',r);
+printf('\n');
+
+Eff=100*[1-(1/r^(G-1))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+T2=T1*(r^(G-1));
+T3=T2*(P3/P2);
+T4=T3/(r^(G-1));
+W=Cv*[T3-T2+(T1-T4)];
+
+R=0.287;
+V1=(R*T1)/P1;
+V2=V1/r;
+
+Pm=W/(V1-V2);
+printf('Mean Effective Pressure %2.1f bar',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.31.1/Ex6_31.sce b/2339/CH6/EX6.31.1/Ex6_31.sce new file mode 100755 index 000000000..5aec29766 --- /dev/null +++ b/2339/CH6/EX6.31.1/Ex6_31.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+P2=44; //in bar
+P3=P2;
+T3=1600+273; //in K
+P1=1; //in bar
+T1=27+273; //in K
+G=1.4;
+
+T2=T1*[(P2/P1)^((G-1)/G)];
+R=0.287;
+
+V1=(R*T1)/(P1*100);
+
+
+r=(P2/P1)^(1/G);
+Z=T3/T2;
+Eff=100*(1-((1/(r^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.32.1/Ex6_32.sce b/2339/CH6/EX6.32.1/Ex6_32.sce new file mode 100755 index 000000000..ba36fe74d --- /dev/null +++ b/2339/CH6/EX6.32.1/Ex6_32.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+r=16; //Compression Ratio
+P1=1; //in bar
+T1=20+273;
+T3=1431+273; //in K
+G=1.4;
+
+T2=T1*(r^(G-1));
+Z=T3/T2;
+T4=(Z^G)*T1;
+Eff=100*[1-((T4-T1)/(G*(T3-T2)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+Cp=1.005;
+Qs=Cp*(T3-T2);
+W=Eff*(Qs/100);
+R=0.287;
+V1=(R*T1)/(P1*100);
+V2=V1/r;
+V=V1-V2;
+
+Pm=W/(V);
+printf('Mean Effective Pressure %2.1f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.33.1/Ex6_33.sce b/2339/CH6/EX6.33.1/Ex6_33.sce new file mode 100755 index 000000000..1d23d007d --- /dev/null +++ b/2339/CH6/EX6.33.1/Ex6_33.sce @@ -0,0 +1,30 @@ +clc
+clear
+
+P1=1; //in bar
+T1=15+273; //in K
+P2=15; //in bar
+P3=40; //in bar
+G=1.4;
+
+r=(P2/P1)^(1/G);
+Eff=100*[1-(1/(r^(G-1)))];
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
+
+T2=T1*[(P2/P1)^((G-1)/G)];
+T3=T2*(P3/P2);
+Cv=0.718;
+
+Qs=Cv*(T3-T2);
+W=Eff*Qs;
+R=0.287;
+
+V1=(R*T1)/(P1*100);
+V2=V1/r;
+
+Vs=V1-V2;
+Pm=W/(Vs*100);
+
+printf('Mean Effective Pressure is %2.1f kPa',Pm);
+printf('\n');
diff --git a/2339/CH6/EX6.4.1/Ex6_4.sce b/2339/CH6/EX6.4.1/Ex6_4.sce new file mode 100755 index 000000000..e562cd038 --- /dev/null +++ b/2339/CH6/EX6.4.1/Ex6_4.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+T1=57+273; //in K
+T2=603+273; //in K
+T3=1950+273; //in K
+T4=870+273; //in K
+G=1.4;
+P1=1; //in bar
+Cp=1.005;
+Cv=0.718;
+
+P2=P1*((T2/T1)^((G)/(G-1)));
+printf('Maximum Pressure attained is %2.1f bar',P2);
+printf('\n');
+
+Qs=Cp*(T3-T2); //Heat Supplied
+Qr=Cv*(T4-T1); //Heat Rejected
+Eff=100*(1-(Qr/Qs));
+printf('Efficiency is %2.0f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.5.1/Ex6_5.sce b/2339/CH6/EX6.5.1/Ex6_5.sce new file mode 100755 index 000000000..0f03d8940 --- /dev/null +++ b/2339/CH6/EX6.5.1/Ex6_5.sce @@ -0,0 +1,15 @@ +clc
+clear
+
+V2=0.2; //in cm^3
+V3=V2;
+Vc=V2;
+
+Vs=1.2; //in cm^3
+V1=V2+Vs;
+G=1.4;
+
+r=(V1/V2);
+Eff=100*(1-(1/(r^(G-1))));
+printf('Efficiency of Engine is %2.0f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.6.1/Ex6_6.sce b/2339/CH6/EX6.6.1/Ex6_6.sce new file mode 100755 index 000000000..c5ee0b4b9 --- /dev/null +++ b/2339/CH6/EX6.6.1/Ex6_6.sce @@ -0,0 +1,12 @@ +clc
+clear
+
+rc=14; //Adiabatic Compression Ratio
+ra=8; //Adiabatic Expansion Ratio
+G=1.4;
+Z=rc/ra; //Cutoff Ratio
+
+//It is a diesel Cycle
+Eff=100*(1-((1/(rc^(G-1)))*(1/G)*((Z^G)-1)/(Z-1)));
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.7.1/Ex6_7.sce b/2339/CH6/EX6.7.1/Ex6_7.sce new file mode 100755 index 000000000..19941559e --- /dev/null +++ b/2339/CH6/EX6.7.1/Ex6_7.sce @@ -0,0 +1,8 @@ +clc
+clear
+
+Eff=0.6; //Efficiency
+T2=283; //in K
+T1=T2/(1-Eff);
+printf('Initial Temperature is %2.1f K',T1);
+printf('\n');
diff --git a/2339/CH6/EX6.8.1/Ex6_8.sce b/2339/CH6/EX6.8.1/Ex6_8.sce new file mode 100755 index 000000000..c9e11bf4d --- /dev/null +++ b/2339/CH6/EX6.8.1/Ex6_8.sce @@ -0,0 +1,15 @@ +clc
+clear
+
+D=10; //Diameter in cm
+L=15; //Length in cm
+Vs=(22/7)*(1/4)*D*D*L; //in cm^3
+Vc=250; //in cm^3
+V2=Vc;
+V1=Vs+Vc;
+r=V1/V2;
+G=1.4;
+
+Eff=100*(1-(1/(r^(G-1))));
+printf('Efficiency is %2.1f Percent',Eff);
+printf('\n');
diff --git a/2339/CH6/EX6.9.1/Ex6_9.sce b/2339/CH6/EX6.9.1/Ex6_9.sce new file mode 100755 index 000000000..78229e7d2 --- /dev/null +++ b/2339/CH6/EX6.9.1/Ex6_9.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+T3=15+273; //in K
+T4=T3;
+P3=1.1; //in bar
+P4=4; //in bar
+P1=12; //in bar
+N=150; //in rpm
+G=1.4;
+
+T1=T4*((P1/P4)^((G-1)/G));
+Eff=100*(1-(T4/T1));
+printf('The Efficiency is %3.2f Percent',Eff);
+printf('\n');
+
+r=P4/P3;
+R=0.287;
+m=1;
+
+W=m*R*(T1-T3)*(log(r));
+P=W*(N/60);
+printf('The Power is %3.1f kW',P);
+printf('\n');
diff --git a/2339/CH7/EX7.1.1/Ex7_1.sce b/2339/CH7/EX7.1.1/Ex7_1.sce new file mode 100755 index 000000000..d61598608 --- /dev/null +++ b/2339/CH7/EX7.1.1/Ex7_1.sce @@ -0,0 +1,10 @@ +clc +clear + +Vs=0.01; //in m^3 +Pm=600; //in kPa +N=300; //in rpm +n=N/2; +IP=(Vs*Pm*n)/60; +printf('Indicated Power= %2.0f kW',IP); +printf('\n'); diff --git a/2339/CH7/EX7.10.1/Ex7_10.sce b/2339/CH7/EX7.10.1/Ex7_10.sce new file mode 100755 index 000000000..abc12a26a --- /dev/null +++ b/2339/CH7/EX7.10.1/Ex7_10.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+D=16;
+L=19;
+Vc=700;
+Pm=5;
+N=1000;
+Eith=0.32;
+Vs=(22/7)*D*D*L*(1/4);
+Vc=700;
+G=1.4;
+r=(Vs+Vc)/Vc;
+Ea=[1-(1/(r^(G-1)))];
+Er=Eith/Ea;
+printf('Relative Efficiency= %3.2f Percent',Er*100);
+printf('\n');
+
+IP=(Pm*100*Vs*(10^-6)*N)/60;
+printf('IP= %3.2f KW',IP);
+printf('\n');
diff --git a/2339/CH7/EX7.11.1/Ex7_11.sce b/2339/CH7/EX7.11.1/Ex7_11.sce new file mode 100755 index 000000000..7a74ac5be --- /dev/null +++ b/2339/CH7/EX7.11.1/Ex7_11.sce @@ -0,0 +1,14 @@ +clc
+clear
+
+T=50;
+Vst=870;
+N=300;
+Pm=10;
+n=N/2;
+
+BP=(2*(22/7)*N*T)/(60*1000);
+IP=(Pm*100*Vst*(10^-6)*N)/(60*2);
+Em=BP/IP;
+printf('Mechanical Efficiency= %3.2f Percent',Em*100);
+printf('\n');
diff --git a/2339/CH7/EX7.12.1/Ex7_12.sce b/2339/CH7/EX7.12.1/Ex7_12.sce new file mode 100755 index 000000000..24e58f365 --- /dev/null +++ b/2339/CH7/EX7.12.1/Ex7_12.sce @@ -0,0 +1,11 @@ +clc
+clear
+
+Pm=7;
+A=(22/7)*(1/4)*((0.15/1.25)^2);
+n=900;
+L=0.15;
+N=2
+IP=(Pm*100*A*L*n*N)/(60*2);
+printf('IP= %3.2f kW',IP);
+printf('\n');
diff --git a/2339/CH7/EX7.13.1/Ex7_13.sce b/2339/CH7/EX7.13.1/Ex7_13.sce new file mode 100755 index 000000000..0d7b242df --- /dev/null +++ b/2339/CH7/EX7.13.1/Ex7_13.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+N=900;
+D=0.1;
+L=0.14;
+Mf=2.1;
+CV=42000;
+Pm=7.5;
+Vc=0.15;
+G=1.4;
+A=(22/7)*(1/4)*D*D;
+IP=(Pm*100*A*L*N*2)/(60*2);
+Eith=(IP*3600)/(Mf*CV);
+printf('Eith= %3.1f Percent',Eith*100);
+printf('\n');
+
+r=(1+0.15)/(0.15);
+Ea=1-[1/(r^(G-1))];
+Er=Eith/Ea;
+printf('Relative Efficiency= %3.2f Percent',Er*100);
+printf('\n');
diff --git a/2339/CH7/EX7.14.1/Ex7_14.sce b/2339/CH7/EX7.14.1/Ex7_14.sce new file mode 100755 index 000000000..f6062bb83 --- /dev/null +++ b/2339/CH7/EX7.14.1/Ex7_14.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+NOC=6;
+N=820;
+n=N/2;
+IP=90;
+LD=1.4;
+Pbm=5;
+Em=0.79;
+BP=IP*Em;
+D=[[IP*60*2]/[Pbm*100*(22/7)*(1/4)*LD*N*NOC]]^(1/3);
+printf('D= %3.0f mm',D*1000);
+printf('\n');
+L=LD*D;
+printf('L= %3.0f mm',L*1000);
+printf('\n');
diff --git a/2339/CH7/EX7.15.1/Ex7_15.sce b/2339/CH7/EX7.15.1/Ex7_15.sce new file mode 100755 index 000000000..d4f324092 --- /dev/null +++ b/2339/CH7/EX7.15.1/Ex7_15.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+NOC=4;
+N=2500;
+n=N/2;
+BP=200;
+LD=1.2;
+Pm=10;
+Em=0.81;
+Mf=65;
+CV=42000;
+IP=BP/Em;
+D=[[IP*60*2*4]/[Pm*100*(22/7)*(1.2*(N)*NOC)]]^(1/3);
+printf('D= %3.0f mm',D*1000);
+printf('\n');
+
+L=LD*D;
+printf('L= %3.0f mm',L*1000);
+printf('\n');
+
+Eith=(IP*3600)/(Mf*CV);
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
+
+Ebth=Eith*Em;
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
diff --git a/2339/CH7/EX7.16.1/Ex7_16.sce b/2339/CH7/EX7.16.1/Ex7_16.sce new file mode 100755 index 000000000..de70d669e --- /dev/null +++ b/2339/CH7/EX7.16.1/Ex7_16.sce @@ -0,0 +1,23 @@ +clc
+clear
+
+IP=42;
+FP=7;
+ES=1800;
+
+BP=IP-FP;
+
+Em=BP/IP;
+printf('Mechanical Efficiency= %3.0f Percent',Em*100);
+printf('\n');
+
+BSFC=0.3;
+CV=43000;
+
+Ebth=3600/(BSFC*CV);
+printf('Brake Thermal Efficiency= %3.0f Percent',Ebth*100);
+printf('\n');
+
+Eith=Ebth/Em;
+printf('Indicated Thermal Efficiency= %3.0f Percent',Eith*100);
+printf('\n');
diff --git a/2339/CH7/EX7.17.1/Ex7_17.sce b/2339/CH7/EX7.17.1/Ex7_17.sce new file mode 100755 index 000000000..c395ec1dd --- /dev/null +++ b/2339/CH7/EX7.17.1/Ex7_17.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+D=0.3;
+L=0.45;
+N=300;
+Pimep=6;
+F=1.5;
+Reff=(180+4)/2;
+
+IP=(Pimep*100*L*(22/7)*(1/4)*(D*D)*N)/(2*60);
+printf('Indicated Power= %3.2f kW',IP);
+printf('\n');
+
+BP=(2*(22/7)*N*F*Reff)/6000;
+printf('Brake Power= %3.2f kW',BP);
+printf('\n');
+
+Em=BP/IP;
+printf('Mechanical Efficiency= %3.2f Percent',Em*100);
+printf('\n');
diff --git a/2339/CH7/EX7.18.1/Ex7_18.sce b/2339/CH7/EX7.18.1/Ex7_18.sce new file mode 100755 index 000000000..e9b45be9d --- /dev/null +++ b/2339/CH7/EX7.18.1/Ex7_18.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+D=0.27;
+L=0.38;
+Pmep=6;
+N=250;
+F=1000;
+Reff=0.75;
+Mf=10;
+CV=44400;
+
+BP=(2*(22/7)*N*(F*Reff))/60;
+printf('Brake Power= %3.2f kW',BP/1000);
+printf('\n');
+
+A=(22/7)*(1/4)*(D*D);
+IP=[Pmep*100*L*A*N]/(2*60);
+printf('Indicated Power= %3.2f kW',IP);
+printf('\n');
+
+Em=BP/(IP*1000);
+printf('Mechanical Efficiency= %3.2f Percent',Em*100);
+printf('\n');
+
+Eith=(IP*3600)/(Mf*CV);
+printf('Indicated Thermal Power= %3.2f Percent',Eith*100);
+printf('\n');
diff --git a/2339/CH7/EX7.19.1/Ex7_19.sce b/2339/CH7/EX7.19.1/Ex7_19.sce new file mode 100755 index 000000000..a9288591c --- /dev/null +++ b/2339/CH7/EX7.19.1/Ex7_19.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+NOC=6;
+IP=89.5;
+N=800;
+LD=1.25;
+Em=0.8;
+Pbemp=5;
+Em=0.8;
+Pimep=Pbemp/0.8;
+
+D3=(IP*2*60*4)/(Pimep*100*LD*(22/7)*N*NOC);
+D=D3^(1/3);
+L=LD*D;
+printf('L= %3.0f mm',L*1000);
+printf('\n');
+printf('D= %3.0f mm',D*1000);
+printf('\n');
diff --git a/2339/CH7/EX7.2.1/Ex7_2.sce b/2339/CH7/EX7.2.1/Ex7_2.sce new file mode 100755 index 000000000..fe08340a3 --- /dev/null +++ b/2339/CH7/EX7.2.1/Ex7_2.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+n=6; //Number of Cylinders
+IP=90; //Indicated Power in kW
+Eff=0.85; //Mechanical Efficiency
+Pmb=5; //in bar
+LD=1.5;
+Pm=Pmb/Eff;
+N=800;
+nx=N/2;
+
+//Length=1.5*D
+D=[[IP*60*4]/[Pm*100*(22/7)*LD*nx*n]]^(1/3);
+printf('D= %3.4f mm',D*100);
+printf('\n');
+L=D*LD;
+printf('L= %3.4f mm',L*100);
+printf('\n');
diff --git a/2339/CH7/EX7.20.1/Ex7_20.sce b/2339/CH7/EX7.20.1/Ex7_20.sce new file mode 100755 index 000000000..0154e5dd2 --- /dev/null +++ b/2339/CH7/EX7.20.1/Ex7_20.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+D=0.25;
+L=0.4;
+Pm=6.5;
+N=250;
+W=1080;
+Ddrum=1.5;
+Mf=10;
+CV=44300;
+
+A=(22/7)*(1/4)*D*D;
+IP=(Pm*100*A*L*N)/(60*2);
+printf('Indicated Power= %3.2f kW',IP);
+printf('\n');
+
+Reff=Ddrum/2;
+W=1.08;
+
+BP=[2*(22/7)*N*W*Reff]/60;
+printf('Brake Power= %3.2f kW',BP);
+printf('\n');
+
+Em=BP/IP;
+Eith=(IP*3600)/(Mf*CV);
+printf('Em= %3.2f Percent',Em*100);
+printf('\n');
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
+
diff --git a/2339/CH7/EX7.21.1/Ex7_21.sce b/2339/CH7/EX7.21.1/Ex7_21.sce new file mode 100755 index 000000000..7f437fbbb --- /dev/null +++ b/2339/CH7/EX7.21.1/Ex7_21.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+W=50;
+S=7;
+D=1.25;
+N=450;
+Mf=4;
+CV=43000;
+Em=0.7;
+Reff=9.81*(D/2);
+
+BP=[2*(22/7)*N*(W-S)*Reff]/(60*1000);
+Ebth=(BP*3600)/(Mf*CV);
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
+
+Eith=Ebth/Em;
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
diff --git a/2339/CH7/EX7.22.1/Ex7_22.sce b/2339/CH7/EX7.22.1/Ex7_22.sce new file mode 100755 index 000000000..1c8cca110 --- /dev/null +++ b/2339/CH7/EX7.22.1/Ex7_22.sce @@ -0,0 +1,33 @@ +clc
+clear
+
+T=640;
+D=0.21;
+N=350;
+L=0.28;
+Pm=5.6;
+Mf=8.16;
+CV=42705;
+
+BP=[2*(22/7)*N*T]/60000;
+printf('Brake Power= %3.2f kW',BP/1000);
+printf('\n');
+
+A=(22/7)*(1/4)*D*D;
+IP=(Pm*100*A*L*N)/60;
+
+Em=BP/IP;
+printf('Em= %3.2f Percent',Em*100);
+printf('\n');
+
+Eith=(IP*3600)/(Mf*CV);
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
+
+Ebth=(BP*3600)/(Mf*CV);
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
+
+BSFC=Mf/BP;
+printf('BSFC= %3.2f kg/kWh',BSFC);
+printf('\n');
diff --git a/2339/CH7/EX7.23.1/Ex7_23.sce b/2339/CH7/EX7.23.1/Ex7_23.sce new file mode 100755 index 000000000..40a62a471 --- /dev/null +++ b/2339/CH7/EX7.23.1/Ex7_23.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+IP=37;
+FP=6;
+BSFC=0.28;
+CV=44300;
+
+BP=IP-FP;
+Em=(IP-FP)/IP;
+printf('Em= %3.2f Percent',Em*100);
+printf('\n');
+
+Mf=BSFC*BP;
+Ebth=(BP*3600)/(Mf*CV);
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
+
+Eith=Ebth/Em;
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
diff --git a/2339/CH7/EX7.24.1/Ex7_24.sce b/2339/CH7/EX7.24.1/Ex7_24.sce new file mode 100755 index 000000000..f9137ca24 --- /dev/null +++ b/2339/CH7/EX7.24.1/Ex7_24.sce @@ -0,0 +1,44 @@ +clc
+clear
+
+D=0.1;
+L=0.125;
+Pm=2.6;
+W=60;
+S=19;
+Reff=0.4;
+r=6;
+Mf=1;
+CV=42000;
+N=2000;
+
+A=(22/7)*(1/4)*D*D;
+
+IP=(Pm*100*A*L*N)/(60*2);
+printf('indicated Power= %3.2f kW',IP);
+printf('\n');
+
+BP=(2*(22/7)*N*(W-S)*Reff)/60000;
+printf('Brake Power= %3.2f kW',BP);
+printf('\n');
+
+Em=BP/IP;
+printf('Em= %3.2f Percent',Em*100);
+printf('\n');
+
+Ebth=(BP*3600)/(Mf*CV);
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
+
+Eith=Ebth/Em;
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
+
+G=1.4;
+Ea=1-[1/(r^(G-1))];
+printf('Ea= %3.2f Percent',Ea*100);
+printf('\n');
+
+Er=Ebth/Ea;
+printf('Er= %3.2f Percent',Er*100);
+printf('\n');
diff --git a/2339/CH7/EX7.25.1/Ex7_25.sce b/2339/CH7/EX7.25.1/Ex7_25.sce new file mode 100755 index 000000000..314f9febd --- /dev/null +++ b/2339/CH7/EX7.25.1/Ex7_25.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+IP=30;
+N=2500;
+Pm=800;
+Em=0.8;
+LD=1.5;
+Ebth=0.28;
+CV=44000;
+
+BP=IP*Em;
+printf('Brake Power= %3.2f kW',BP);
+printf('\n');
+
+Mf=(BP/(Ebth*CV));
+printf('Mass Flow Rate= %3.2f kg/hr',Mf*3600);
+printf('\n');
diff --git a/2339/CH7/EX7.3.1/Ex7_3.sce b/2339/CH7/EX7.3.1/Ex7_3.sce new file mode 100755 index 000000000..33ed305bb --- /dev/null +++ b/2339/CH7/EX7.3.1/Ex7_3.sce @@ -0,0 +1,19 @@ +clc
+clear
+
+BP=22; //Brake Power
+Eff=0.85; //Mechanical Efficiency
+IP=BP/Eff;
+mf=6.5;
+CV=30000; //Calorific Value
+Ebth=BP/((mf/3600)*CV);
+printf('Brake Thermal Eff= %3.1f Percent',Ebth*100);
+printf('\n');
+
+Eith=IP/((mf/3600)*CV);
+printf('Indicated Thermal Eff= %3.1f Percent',Eith*100);
+printf('\n');
+
+BSFC=mf/BP;
+printf('BSFC= %3.1f kg/kWh',BSFC);
+printf('\n');
diff --git a/2339/CH7/EX7.4.1/Ex7_4.sce b/2339/CH7/EX7.4.1/Ex7_4.sce new file mode 100755 index 000000000..dc830d3c2 --- /dev/null +++ b/2339/CH7/EX7.4.1/Ex7_4.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+BP=185; //Brake Power
+Eff=0.75;
+IP=BP/Eff;
+LD=1.5;
+N=35;
+n=N/2;
+nx=4;
+Pm=830; //in kPa
+D=[[IP*4]/[Pm*(22/7)*LD*nx*n]]^(1/3);
+printf('D= %3.0f mm',D*1000);
+printf('\n');
+L=D*LD;
+printf('L= %3.0f mm',L*1000);
+printf('\n');
diff --git a/2339/CH7/EX7.5.1/Ex7_5.sce b/2339/CH7/EX7.5.1/Ex7_5.sce new file mode 100755 index 000000000..ff9d8119b --- /dev/null +++ b/2339/CH7/EX7.5.1/Ex7_5.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+Vc=5*(10^-4);
+D=0.15;
+L=0.2;
+Vs=(22/7)*D*D*L*(1/4);
+r=(Vc+Vs)/Vc;
+G=1.4;
+Ea=[1-(1/(r^(G-1)))];
+Eith=0.3;
+Erel=Eith/Ea;
+printf('Erel= %3.2f Percent',Erel*100);
+printf('\n');
+
+Pm=500; //in kPa
+n=1000/2;
+IP=(Pm*Vs*n)/60;
+printf('IP= %3.2f kW',IP);
+printf('\n');
diff --git a/2339/CH7/EX7.6.1/Ex7_6.sce b/2339/CH7/EX7.6.1/Ex7_6.sce new file mode 100755 index 000000000..91ab3f3f8 --- /dev/null +++ b/2339/CH7/EX7.6.1/Ex7_6.sce @@ -0,0 +1,11 @@ +clc
+clear
+
+Pm=600;
+A=(22/7)*(1/4)*0.11*0.11*0.14;
+n=1000;
+IP=(Pm*A*n)/60;
+Em=0.8;
+BP=Em*IP;
+printf('BP= %3.2f kW',BP);
+printf('\n');
diff --git a/2339/CH7/EX7.7.1/Ex7_7.sce b/2339/CH7/EX7.7.1/Ex7_7.sce new file mode 100755 index 000000000..51e235dfc --- /dev/null +++ b/2339/CH7/EX7.7.1/Ex7_7.sce @@ -0,0 +1,12 @@ +clc
+clear
+
+r=6;
+G=1.4;
+Ea=100*[1-(1/(r^(G-1)))];
+Ebt=Ea/2;
+CV=41500;
+BP=15;
+Mf=BP/(CV*(Ebt/100));
+printf('Mf= %3.2f kg/hr',Mf*3600);
+printf('\n');
diff --git a/2339/CH7/EX7.8.1/Ex7_8.sce b/2339/CH7/EX7.8.1/Ex7_8.sce new file mode 100755 index 000000000..9d3843bae --- /dev/null +++ b/2339/CH7/EX7.8.1/Ex7_8.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+n=4;
+;
+
+DL=1.2;
+BP=32;
+N=2500;
+Pm=9;
+Em=0.86;
+Mf=9;
+CV=43000;
+
+IP=BP/Em;
+D=[[IP*60*4]/[Pm*100*(22/7)*DL*N*n]]^(1/3);
+printf('D= %3.0f mm',D*1000);
+printf('\n');
+
+L=DL*D;
+printf('L= %3.0f mm',L*1000);
+printf('\n');
+
+Ebth=BP/(Mf*CV/3600);
+printf('Ebth= %3.2f Percent',Ebth*100);
+printf('\n');
+
+Eith=Ebth/Em;
+printf('Eith= %3.2f Percent',Eith*100);
+printf('\n');
+
diff --git a/2339/CH7/EX7.9.1/Ex7_9.sce b/2339/CH7/EX7.9.1/Ex7_9.sce new file mode 100755 index 000000000..3f54ff3ba --- /dev/null +++ b/2339/CH7/EX7.9.1/Ex7_9.sce @@ -0,0 +1,14 @@ +clc
+clear
+
+Eith=0.29;
+Em=0.77;
+BP=5.5;
+SG=0.87;
+CV=43000;
+Ebth=Em*Eith;
+Mf=(BP*3600)/(Ebth*CV);
+D=SG*1000;
+Mff=(Mf*1000)/D
+printf('Mf= %3.2f litre/hr',Mff);
+printf('\n');
diff --git a/2339/CH8/EX8.1.1/Ex8_1.sce b/2339/CH8/EX8.1.1/Ex8_1.sce new file mode 100755 index 000000000..0bc9df6a8 --- /dev/null +++ b/2339/CH8/EX8.1.1/Ex8_1.sce @@ -0,0 +1,12 @@ +clc +clear + +P1=100; +T1=300; +P2=650; +n=1.25; +r=0.05; + +Ev=1-[r*(((P2/P1)^(1/n))-1)]; +printf('Volumetric Efficiency= %2.2f Percent',Ev*100); +printf('\n'); diff --git a/2339/CH8/EX8.10.1/Ex8_10.sce b/2339/CH8/EX8.10.1/Ex8_10.sce new file mode 100755 index 000000000..a17fbf221 --- /dev/null +++ b/2339/CH8/EX8.10.1/Ex8_10.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+IP=15;
+n=1.2;
+P1=100;
+P2=700;
+x=[(P2/P1)^((n-1)/n)]-1;
+V1N=[IP*(n-1)*60]/[n*P1*x*2];
+LN=150/2;
+D2=V1N*4/[(22/7)*LN];
+D=D2^0.5;
+L=D*1.5;
+printf('D= %2.0f mm',D*1000);
+printf('\n');
+printf('L= %2.0f mm',L*1000);
+printf('\n');
diff --git a/2339/CH8/EX8.11.1/Ex8_11.sce b/2339/CH8/EX8.11.1/Ex8_11.sce new file mode 100755 index 000000000..66628804c --- /dev/null +++ b/2339/CH8/EX8.11.1/Ex8_11.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+P1=1;
+P2=16;
+n=1.3;
+LN=100;
+N=350;
+IP=30;
+Ev=0.95;
+
+L=LN/N;
+x=[((P2/P1)^((n-1)/n))-1];
+V14=[IP*(n-1)*60]/[n*P1*100*x*N];
+Vs=V14/Ev;
+D2=Vs*4/[(22/7)*L];
+D=D2^0.5;
+printf('D= %2.0f mm',D*1000);
+printf('\n');
+printf('L= %2.0f mm',L*1000);
+printf('\n');
diff --git a/2339/CH8/EX8.12.1/Ex8_12.sce b/2339/CH8/EX8.12.1/Ex8_12.sce new file mode 100755 index 000000000..8c7252acd --- /dev/null +++ b/2339/CH8/EX8.12.1/Ex8_12.sce @@ -0,0 +1,62 @@ +clc
+clear
+
+D=0.2;
+L=0.3;
+Vs=(22/7)*(1/4)*D*D*L;
+Vc=0.04*Vs;
+
+V3=Vc;
+P2=9;
+P1=1;
+n=1.3;
+V4=V3*[(P2/P1)^(1/n)];
+V1=Vs+Vc;
+W=[n/(n-1)]*[P1*100]*[V1-V4]*[((P2/P1)^((n-1)/n))-1];
+R=0.287;
+T1=15+273;
+Md=[P1*(V1-V4)*100]/[R*T1];
+Wpkg=W/Md;
+printf('Word done per kg: %3.2f kJ/kg of air',Wpkg);
+printf('\n');
+
+T2=T1*[(P2/P1)^((n-1)/n)];
+G=1.4;
+Q=[(G-n)/(G-1)]*[(R*(T1-T2))/(n-1)];
+printf('Heat Transfereed: %3.2f kJ/kg',Q);
+printf('\n');
+
+Pm=W/Vs;
+printf('Mean Effective Pressure: %3.2f kPa',Pm);
+printf('\n');
+
+Mac=V1/(V1-V4);
+printf('Mass of air compressed to delivered: %3.2f ',Mac);
+printf('\n');
+
+Tf=T1;
+Pf=101.325;
+Vf=[P1*100*(V1-V4)*Tf]/[Pf*T1];
+RPM=500;
+Vf=Vf*RPM;
+printf('FAD at standard condition: %3.2f m^3/min',Vf);
+printf('\n');
+
+IP=[W*RPM]/60;
+Etrans=0.92;
+Emech=0.85;
+Emotor=0.75;
+MP=IP/[Etrans*Emech*Emotor];
+printf('Motor Power: %3.2f kW',MP);
+printf('\n');
+
+MAC=Md*RPM;
+printf('Mass of air compressed: %3.2f kg/min',MAC);
+printf('\n');
+
+ACC=MAC*Mac;
+printf('Air compressed in cylinder: %3.2f kg/min',ACC);
+printf('\n');
+
+printf('End Temperature: %3.2f K',T2);
+printf('\n');
diff --git a/2339/CH8/EX8.13.1/Ex8_13.sce b/2339/CH8/EX8.13.1/Ex8_13.sce new file mode 100755 index 000000000..ec05a56da --- /dev/null +++ b/2339/CH8/EX8.13.1/Ex8_13.sce @@ -0,0 +1,27 @@ +clc
+clear
+
+Vs=0.015;
+Vc=0.06*Vs;
+V3=Vc;
+n=1.3;
+P2=10;
+P1=1;
+N=280;
+
+V4=V3*[(P2/P1)^(1/n)];
+printf('V4: %3.4f m^3/cycle',V4);
+printf('\n');
+printf('V3: %3.4f m^3/cycle',V3);
+printf('\n');
+
+V1=Vs+Vc;
+printf('V1: %3.4f m^3/cycle',V1);
+printf('\n');
+
+V14=V1-V4; //Suction Volume
+V2=V1*[(P1/P2)^(1/n)];
+IP=[n/(n-1)]*[P1*100*(V14)]*[((P2/P1)^((n-1)/n))-1]*[N/60];
+printf('IP: %3.0f kW',IP);
+printf('\n');
+
diff --git a/2339/CH8/EX8.14.1/Ex8_14.sce b/2339/CH8/EX8.14.1/Ex8_14.sce new file mode 100755 index 000000000..0648bd3e1 --- /dev/null +++ b/2339/CH8/EX8.14.1/Ex8_14.sce @@ -0,0 +1,35 @@ +clc
+clear
+
+P2=6;
+P1=0.96;
+n=1.3;
+CV=0.04;
+
+Ev=[1-[CV*[((P2/P1)^(1/n))-1]]]*100;
+printf('Clearance Volumetric Efficiency: %3.1f Percent',Ev);
+printf('\n');
+
+D=0.09;
+L=0.1;
+Vs=(22/7)*(1/4)*(D*D*L);
+Vc=0.04*Vs;
+V4=Vc*[(P2/P1)^(1/n)];
+V1=Vc+Vs;
+EDV=V1-V4;
+printf('Effective Displacement Volume: %3.5f m^3',EDV);
+printf('\n');
+
+T1=313;
+Tf=293;
+Pf=1;
+Vf=[(P1*(EDV)*Tf)]/[T1*Pf];
+N=410;
+FAD=Vf*N*2*60;
+printf('Free air delivered: %3.2f m^3',FAD);
+printf('\n');
+
+W=[n/(n-1)]*[P1*100*(V1-V4)]*[((P2/P1)^((n-1)/n))-1];
+IP=W*2*N/60;
+printf('Indicated Power: %3.2f kW',IP);
+printf('\n');
diff --git a/2339/CH8/EX8.15.1/Ex8_15.sce b/2339/CH8/EX8.15.1/Ex8_15.sce new file mode 100755 index 000000000..eca769cea --- /dev/null +++ b/2339/CH8/EX8.15.1/Ex8_15.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+P1=1;
+P2=5;
+T1=27+273;
+m=1;
+R=0.287;
+
+W1=m*R*T1*(log(P2/P1));
+printf('Work in isothermal process: %3.1f kJ',W1);
+printf('\n');
+
+G=1.4;
+W2=[G/(G-1)]*[m*R*T1]*[((P2/P1)^((G-1)/G))-1];
+printf('Work in isentropic process: %3.0f kJ',W2);
+printf('\n');
+
+n=1.25;
+W3=[n/(n-1)]*[m*R*T1]*[((P2/P1)^((n-1)/n))-1];
+printf('Work in polytropic process: %3.1f kJ',W3);
+printf('\n');
diff --git a/2339/CH8/EX8.16.1/Ex8_16.sce b/2339/CH8/EX8.16.1/Ex8_16.sce new file mode 100755 index 000000000..d6562b287 --- /dev/null +++ b/2339/CH8/EX8.16.1/Ex8_16.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+IP=41;
+P1=1;
+T1=17+273;
+P2=7;
+N=100;
+n=1.2;
+
+L=150/[2*N];
+V1=(22/7)*(1/4)*(L); //Along with D^2
+W=[n/(n-1)]*[P1*100*V1]*[((P2/P1)^((n-1)/n))-1];
+D2=[IP*60]/[W*2*N];
+D=sqrt(D2);
+printf('D: %3.3f m',D);
+printf('\n');
diff --git a/2339/CH8/EX8.17.1/Ex8_17.sce b/2339/CH8/EX8.17.1/Ex8_17.sce new file mode 100755 index 000000000..7f8d4660c --- /dev/null +++ b/2339/CH8/EX8.17.1/Ex8_17.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+D=0.15;
+L=0.2;
+P1=1;
+T1=17+273;
+P2=7;
+N=100;
+R=0.287;
+V1=(22/7)*(1/4)*D*D*L;
+m=[P1*100*V1]/[R*T1];
+Mpm=m*N;
+n=1.25;
+
+IP=[n/(n-1)]*[P1*100*V1]*[((P2/P1)^((n-1)/n))-1]*[N/60];
+printf('Mass/min: %3.1f Mpm',Mpm);
+printf('\n');
+
+printf('Indicated Power: %3.1f kW',IP);
+printf('\n');
+
+T2=T1*[(P2/P1)^((n-1)/n)];
+printf('T2: %3.1f K',T2);
+printf('\n');
diff --git a/2339/CH8/EX8.18.1/Ex8_18.sce b/2339/CH8/EX8.18.1/Ex8_18.sce new file mode 100755 index 000000000..197075b5c --- /dev/null +++ b/2339/CH8/EX8.18.1/Ex8_18.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+D=0.15;
+N=100;
+L=0.2;
+P1=1;
+T1=27+273;
+P2=6;
+n=1.25;
+
+Vs=(22/7)*(1/4)*D*D*L;
+Vc=0.05*Vs;
+V1=Vs+Vc;
+V4=Vc*[(P2/P1)^(1/n)];
+
+IP=[n/(n-1)]*[P1*100*(V1-V4)]*[((P2/P1)^((n-1)/n))-1];
+IPf=IP**(N/60)
+printf('IP: %3.2f kJ',IPf);
+printf('\n');
+
+Pm=IP/Vs;
+printf('Mean Effective Pressure: %3.2f kN/m^2',Pm);
+printf('\n');
diff --git a/2339/CH8/EX8.19.1/Ex8_19.sce b/2339/CH8/EX8.19.1/Ex8_19.sce new file mode 100755 index 000000000..0cfc6cc44 --- /dev/null +++ b/2339/CH8/EX8.19.1/Ex8_19.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+n=1.2;
+m=5;
+R=0.287;
+T2=107+273;
+T1=27+273;
+IP=[n/(n-1)]*[m/60]*[R*(T2-T1)];
+printf('Air Power: %3.2f kW',IP);
+printf('\n');
+
+BP=14;
+Em=IP*100/BP;
+printf('Mechanical Efficiency: %3.0f Percent',Em);
+printf('\n');
diff --git a/2339/CH8/EX8.2.1/Ex8_2.sce b/2339/CH8/EX8.2.1/Ex8_2.sce new file mode 100755 index 000000000..d612dc4b6 --- /dev/null +++ b/2339/CH8/EX8.2.1/Ex8_2.sce @@ -0,0 +1,14 @@ +clc
+clear
+
+D=0.24;
+LN=5/6;
+P1=100;
+P2=1000;
+n=1.35;
+
+A=(22/7)*(1/4)*D*D;
+
+IP=[n/(n-1)]*[P1*A*LN]*[((P2/P1)^((n-1)/n))-1];
+printf('Indicated Power= %2.2f kW',IP);
+printf('\n');
diff --git a/2339/CH8/EX8.20.1/Ex8_20.sce b/2339/CH8/EX8.20.1/Ex8_20.sce new file mode 100755 index 000000000..e92014092 --- /dev/null +++ b/2339/CH8/EX8.20.1/Ex8_20.sce @@ -0,0 +1,21 @@ +clc
+clear
+
+V1=50;
+P1=1;
+P2=5.5;
+n=1.3;
+Em=0.82;
+
+IP=[n/(n-1)]*[P1*100*V1]*[((P2/P1)^((n-1)/n))-1]*[1/60];
+BP=IP/Em;
+
+printf('IP: %3.1f kW',IP);
+printf('\n');
+printf('BP: %3.1f kW',BP);
+printf('\n');
+
+IsoP=P1*100*V1*[log(P2/P1)]*(100/60);
+Eo=IsoP/BP;
+printf('Isothermal Efficiecy: %3.1f Percent',Eo);
+printf('\n');
diff --git a/2339/CH8/EX8.21.1/Ex8_21.sce b/2339/CH8/EX8.21.1/Ex8_21.sce new file mode 100755 index 000000000..152469966 --- /dev/null +++ b/2339/CH8/EX8.21.1/Ex8_21.sce @@ -0,0 +1,14 @@ +clc
+clear
+
+P1=1;
+P2=5.5;
+T1=27+273;
+Pa=1.01325;
+Ta=17+273;
+C=0.06;
+n=1.3;
+
+Ev=[(P1*Ta)/(Pa*T1)]*[1+C-(C*((P2/P1)^(1/n)))];
+printf('Volumetric Efficiency: %3.0f Percent',Ev*100);
+printf('\n');
diff --git a/2339/CH8/EX8.22.1/Ex8_22.sce b/2339/CH8/EX8.22.1/Ex8_22.sce new file mode 100755 index 000000000..c70d2ccc5 --- /dev/null +++ b/2339/CH8/EX8.22.1/Ex8_22.sce @@ -0,0 +1,31 @@ +clc
+clear
+
+V14=7.5;
+P1=1;
+T1=27+273;
+P2=5.5;
+n=1.3;
+C=0.06;
+
+T2=T1*[(P2/P1)^((n-1)/n)];
+printf('T2: %3.1f K',T2);
+printf('\n');
+
+Ev=1+C-[C*((P2/P1)^((1/n)))];
+printf('Vol Eff: %3.1f Percent',Ev*100);
+printf('\n');
+
+AP=[n/(n-1)]*[P1*100*V14/60]*[((P2/P1)^((n-1)/n))-1];
+printf('Air Power: %3.1f kW',AP);
+printf('\n');
+
+Em=0.9;
+BP=AP/Em;
+printf('BP: %3.1f kW',BP);
+printf('\n');
+
+Emot=0.96;
+EMC=BP/Emot;
+printf('Electric Motor Capacity: %3.1f kW',EMC);
+printf('\n');
diff --git a/2339/CH8/EX8.23.1/Ex8_23.sce b/2339/CH8/EX8.23.1/Ex8_23.sce new file mode 100755 index 000000000..69c308dc5 --- /dev/null +++ b/2339/CH8/EX8.23.1/Ex8_23.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+V1=5;
+P1=1;
+P2=5;
+n=1.25;
+Em=0.9;
+IP=[n/(n-1)]*[P1*100*V1/60]*[((P2/P1)^((n-1)/n))-1];
+SP=IP/Em;
+printf('Shaft Power: %3.1f kW',SP);
+printf('\n');
+
+IsoP=P1*100*V1*(log(P2/P1))*(1/60);
+Eo=IsoP/SP;
+printf('Overall Efficiency: %3.0f Percent',Eo*100);
+printf('\n');
diff --git a/2339/CH8/EX8.24.1/Ex8_24.sce b/2339/CH8/EX8.24.1/Ex8_24.sce new file mode 100755 index 000000000..9947d2239 --- /dev/null +++ b/2339/CH8/EX8.24.1/Ex8_24.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+V1=25;
+P1=1;
+P2=7;
+N=460;
+Em=0.8;
+Ev=0.76;
+Ei=0.81;
+
+IsoP=P1*100*V1*(log(P2/P1));
+IndP=IsoP/Ei;
+Vs=V1/Ev;
+Pm=IndP/Vs;
+BP=IndP/(3600*Em);
+
+printf('Mean Effective Pressure: %3.2f bar',Pm/100);
+printf('\n');
+
+printf('BP: %3.2f kW',BP);
+printf('\n');
diff --git a/2339/CH8/EX8.25.1/Ex8_25.sce b/2339/CH8/EX8.25.1/Ex8_25.sce new file mode 100755 index 000000000..4fe171a19 --- /dev/null +++ b/2339/CH8/EX8.25.1/Ex8_25.sce @@ -0,0 +1,32 @@ +clc
+clear
+
+Va=3;
+Pa=1;
+Ta=17+273;
+P2=8.2;
+N=300;
+n=1.35
+LD=1.2;
+Em=0.9;
+C=0.05;
+
+P1=Pa-0.05;
+T1=Ta+10;
+V14=[Pa*Va*T1]/[P1*Ta];
+
+IP=[n/(n-1)]*[P1*100*V14/60]*[((P2/P1)^((n-1)/n))-1];
+BP=IP/Em;
+printf('BP: %3.1f kW',BP);
+printf('\n');
+
+Ev=1+C-[C*((P2/P1)^(1/n))];
+printf('Volumetric Efficiency: %3.1f Percent',Ev*100);
+printf('\n');
+
+Vs=(22/7)*(1/4)*LD;
+VsMin=Vs*2*N;
+D3=V14/[VsMin*Ev];
+D=D3^(1/3);
+printf('Cylinder Diameter: %3.0f mm',D*1000);
+printf('\n');
diff --git a/2339/CH8/EX8.26.1/Ex8_26.sce b/2339/CH8/EX8.26.1/Ex8_26.sce new file mode 100755 index 000000000..01962a3da --- /dev/null +++ b/2339/CH8/EX8.26.1/Ex8_26.sce @@ -0,0 +1,26 @@ +clc
+clear
+
+V1=1;
+P1=1.013;
+T1=15+273;
+P2=7;
+R=0.287;
+n=1.35;
+
+m=[P1*100*V1]/[R*T1];
+printf('Mass of air per minute: %3.1f kg',m);
+printf('\n');
+
+T2=T1*[((P2/P1)^((n-1)/n))];
+printf('T2: %3.1f K',T2);
+printf('\n');
+
+IP=[n/(n-1)]*[P1*100*V1/60]*[((P2/P1)^((n-1)/n))-1];
+printf('IP: %3.1f kW',IP);
+printf('\n');
+
+IsoP=P1*100*V1*(1/60)*log(P2/P1);
+Ei=IsoP/IP;
+printf('Isothermal Efficiency: %3.0f Percent',Ei*100);
+printf('\n');
diff --git a/2339/CH8/EX8.27.1/Ex8_27.sce b/2339/CH8/EX8.27.1/Ex8_27.sce new file mode 100755 index 000000000..d50752580 --- /dev/null +++ b/2339/CH8/EX8.27.1/Ex8_27.sce @@ -0,0 +1,36 @@ +clc
+clear
+
+P1=1.013;
+T1=15+273;
+P2=7;
+FAD=0.3;
+G=1.4;
+
+IP=[G/(G-1)]*[P1*100*FAD/60]*[((P2/P1)^((G-1)/G))-1];
+printf('For Isentropic process \n');
+printf('IP: %3.1f kW',IP);
+printf('\n');
+
+T2=T1*[(P2/P1)^((G-1)/G)];
+printf('T2: %3.0f K',T2);
+printf('\n');
+
+printf('For Reversible Isothermal process \n');
+IP=P1*100*FAD*(1/60)*[log(P2/P1)];
+printf('IP: %3.3f kW',IP);
+printf('\n');
+
+T2=T1;
+printf('T2: %3.0f K',T2);
+printf('\n');
+
+printf('For Polytropic process \n');
+n=1.25
+IP=[n/(n-1)]*[P1*100*FAD/60]*[((P2/P1)^((n-1)/n))-1];
+printf('IP: %3.3f kW',IP);
+printf('\n');
+
+T2=T1*[(P2/P1)^((n-1)/n)];
+printf('T2: %3.2f K',T2);
+printf('\n');
diff --git a/2339/CH8/EX8.28.1/Ex8_28.sce b/2339/CH8/EX8.28.1/Ex8_28.sce new file mode 100755 index 000000000..f307af183 --- /dev/null +++ b/2339/CH8/EX8.28.1/Ex8_28.sce @@ -0,0 +1,53 @@ +clc
+clear
+
+V1=94;
+P1=1;
+T1=25+273;
+P2=9;
+
+printf('For isothermal process \n');
+T2=T1;
+printf('T2: %3.0f K',T2);
+printf('\n');
+
+P=P1*100*V1*log(P2/P1);
+printf('Power required: %3.0f kW',P/60);
+printf('\n');
+
+Q=P;
+printf('Heat Rejected: %3.0f KW',Q/60);
+printf('\n');
+
+printf('\n For adiabatic process \n');
+G=1.4;
+T2=T1*[(P2/P1)^((G-1)/G)];
+printf('T2: %3.0f K',T2);
+printf('\n');
+
+P=[G/(G-1)]*[P1*100*V1/60]*[((P2/P1)^((G-1)/G))-1];
+printf('Power required: %3.0f kW',P);
+printf('\n');
+
+Q=0;
+printf('Heat Rejected: %3.0f kW',Q);
+printf('\n');
+
+printf('\n For adiabatic process \n');
+n=1.25;
+T2=T1*[(P2/P1)^((n-1)/n)];
+printf('T2: %3.0f K',T2);
+printf('\n');
+
+P=[n/(n-1)]*[P1*100*V1/60]*[((P2/P1)^((n-1)/n))-1];
+printf('Power required: %3.0f kW',P);
+printf('\n');
+
+R=0.287;
+Cp=1.005;
+
+m=[P1*100*V1]/[R*T1];
+H=m*(1/60)*Cp*(T2-T1);
+Q=H-P;
+printf('Heat Rejected: %3.0f kW',Q);
+printf('\n');
diff --git a/2339/CH8/EX8.29.1/Ex8_29.sce b/2339/CH8/EX8.29.1/Ex8_29.sce new file mode 100755 index 000000000..b1126097d --- /dev/null +++ b/2339/CH8/EX8.29.1/Ex8_29.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+P1=1;
+P2=12;
+n=1.3;
+N=350;
+L=180/(2*N);
+IP=30;
+Ev=0.92;
+
+W=[n/(n-1)]*[P1*100]*[((P2/P1)^((n-1)/n))-1]; //with (V1-V4)
+V14=[IP*60]/[N*W];
+Vs=V14/Ev;
+D2=Vs*4/[(22/7)*L];
+D=sqrt(D2);
+printf('D: %3.3f m',D);
+printf('\n');
+printf('L: %3.3f m',L);
+printf('\n');
diff --git a/2339/CH8/EX8.3.1/Ex8_3.sce b/2339/CH8/EX8.3.1/Ex8_3.sce new file mode 100755 index 000000000..217eee7fa --- /dev/null +++ b/2339/CH8/EX8.3.1/Ex8_3.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+N=300;
+D=0.2;
+L=0.24;
+P1=1.01325;
+P2=8*1.01325;
+n=1.35;
+Et=0.96;
+Em=0.85;
+Vs=(22/7)*(1/4)*D*D*L;
+
+IP=[n/(n-1)]*[P1*Vs]*[N/60]*[((P2/P1)^((n-1)/n))-1];
+printf('Indicated Power= %2.1f kW',IP*100);
+printf('\n');
+
+BP=IP/(Et*Em);
+printf('Brake Power= %2.1f kW',BP*100);
+printf('\n');
diff --git a/2339/CH8/EX8.30.1/Ex8_30.sce b/2339/CH8/EX8.30.1/Ex8_30.sce new file mode 100755 index 000000000..df80e36c9 --- /dev/null +++ b/2339/CH8/EX8.30.1/Ex8_30.sce @@ -0,0 +1,22 @@ +clc
+clear
+
+m=1;
+P1=1;
+P2=5;
+T1=27+273;
+n=1.25;
+R=0.287;
+
+W=m*R*T1*log(P2/P1);
+printf('Work Done for Isothermal Process: %3.2f kJ/kg \n\n',W);
+printf('\n');
+
+G=1.4;
+W=[G/(G-1)]*[m*R*T1]*[((P2/P1)^((G-1)/G))-1];
+printf('Work Done for Isentropic Process: %3.2f kJ/kg \n\n',W);
+printf('\n');
+
+W=[n/(n-1)]*[m*R*T1]*[((P2/P1)^((n-1)/n))-1];
+printf('Work Done for Polytropic Process: %3.2f kJ/kg \n\n',W);
+printf('\n');
diff --git a/2339/CH8/EX8.31.1/Ex8_31.sce b/2339/CH8/EX8.31.1/Ex8_31.sce new file mode 100755 index 000000000..ac48187d0 --- /dev/null +++ b/2339/CH8/EX8.31.1/Ex8_31.sce @@ -0,0 +1,24 @@ +clc
+clear
+
+D=0.15;
+L=0.3;
+P1=1;
+T1=27+273;
+P2=8;
+N=120;
+G=1.4;
+R=0.287;
+Vs=(22/7)*(1/4)*D*D*L;
+
+m=[P1*100*Vs]/[R*T1];
+printf('Mass of air compressed per cycle: %3.4f kJ/cycle ',m);
+printf('\n');
+
+W=[G/(G-1)]*[P1*100*Vs]*[((P2/P1)^((G-1)/G))-1];
+printf('Work required per cycle: %3.3f kJ/cycle ',W);
+printf('\n');
+
+P=(W*N)/60;
+printf('Power required to drive compressor: %3.2f kJ/cycle ',P);
+printf('\n');
diff --git a/2339/CH8/EX8.4.1/Ex8_4.sce b/2339/CH8/EX8.4.1/Ex8_4.sce new file mode 100755 index 000000000..fe049f670 --- /dev/null +++ b/2339/CH8/EX8.4.1/Ex8_4.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+n=1.35;
+P1=1.013;
+V1=1/60;
+P2=7;
+Et=0.85;
+Em=0.9;
+
+IP=[n/(n-1)]*[P1*100*V1]*[((P2/P1)^((n-1)/n))-1];
+printf('Indicated Power= %2.1f kW',IP);
+printf('\n');
+
+BP=IP/(Et*Em);
+printf('Brake Power= %2.1f kW',BP);
+printf('\n');
diff --git a/2339/CH8/EX8.5.1/Ex8_5.sce b/2339/CH8/EX8.5.1/Ex8_5.sce new file mode 100755 index 000000000..5ed415f52 --- /dev/null +++ b/2339/CH8/EX8.5.1/Ex8_5.sce @@ -0,0 +1,15 @@ +clc
+clear
+
+n=1.2;
+P1=1;
+P2=6;
+Vs=1.5/60;
+
+IP=[n/(n-1)]*[P1*100*Vs]*[((P2/P1)^((n-1)/n))-1];
+printf('Indicated Power= %2.1f kW',IP);
+printf('\n');
+MP=6.55;
+Em=IP/MP;
+printf('Mechanical Efficiency= %2.1f Percent',Em*100);
+printf('\n');
diff --git a/2339/CH8/EX8.6.1/Ex8_6.sce b/2339/CH8/EX8.6.1/Ex8_6.sce new file mode 100755 index 000000000..4e7a91298 --- /dev/null +++ b/2339/CH8/EX8.6.1/Ex8_6.sce @@ -0,0 +1,20 @@ +clc
+clear
+
+N=300;
+V14=14/(2*N);
+Vs=0.023/(1.05-0.22);
+n=1.3;
+P1=1.013;
+P2=7;
+IP=[n/(n-1)]*[P1*100*V14]*[((P2/P1)^((n-1)/n))-1]*[2*N/60];
+printf('Indicated Power= %2.1f kW',IP);
+printf('\n');
+
+T1=288;
+T2=T1*[(P2/P1)^((n-1)/n)];
+printf('Delivery Temperature= %2.0f K',T2);
+printf('\n');
+
+printf('Swept Volume= %2.4f m^3',Vs);
+printf('\n');
diff --git a/2339/CH8/EX8.7.1/Ex8_7.sce b/2339/CH8/EX8.7.1/Ex8_7.sce new file mode 100755 index 000000000..e4a95e05a --- /dev/null +++ b/2339/CH8/EX8.7.1/Ex8_7.sce @@ -0,0 +1,30 @@ +clc
+clear
+
+P1=1;
+P2=10;
+Vs=0.015;
+FAD=3;
+Vc=Vs*0.06;
+n=1.3;
+T1=20+273;
+
+IP=[n/(n-1)]*[P1*100*3]*[((P2/P1)^((n-1)/n))-1];
+printf('Indicated Power= %2.1f kW',IP/60);
+printf('\n');
+
+V4=Vc*[(P2/P1)^(1/n)];
+V1=Vs+Vc;
+V14=0.0107;
+RS=3/V14;
+printf('Rotation Speed= %2.0f RPM',RS);
+printf('\n');
+Tf=288;
+Pf=101.325;
+Vf=[P1*100*(FAD)*Tf]/[T1*Pf];
+printf('Vf= %2.4f m^3/min',Vf);
+printf('\n');
+
+Mcd=V1/(V14);
+printf('Mcd= %2.1f ',Mcd);
+printf('\n');
diff --git a/2339/CH8/EX8.8.1/Ex8_8.sce b/2339/CH8/EX8.8.1/Ex8_8.sce new file mode 100755 index 000000000..2c434614f --- /dev/null +++ b/2339/CH8/EX8.8.1/Ex8_8.sce @@ -0,0 +1,17 @@ +clc
+clear
+
+P1=1;
+P2=10;
+Vs=0.014;
+n=1.3;
+V1=3;
+FAD=3;
+
+W=[n/(n-1)]*[P1*100*V1/60]*[((P2/P1)^((n-1)/n))-1];
+printf('Power required= %2.1f kW',W);
+printf('\n');
+
+RPM=FAD/Vs;
+printf('Rotational Speed= %2.0f rpm',RPM);
+printf('\n');
diff --git a/2339/CH8/EX8.9.1/Ex8_9.sce b/2339/CH8/EX8.9.1/Ex8_9.sce new file mode 100755 index 000000000..984387457 --- /dev/null +++ b/2339/CH8/EX8.9.1/Ex8_9.sce @@ -0,0 +1,28 @@ +clc
+clear
+
+Vs=5.665/600;
+Vc=0.04*Vs;
+V3=Vc;
+n=1.3;
+P3=5.6;
+P2=0.97;
+V4=V3*[(P3/P2)^(1/n)];
+V1=Vs+Vc;
+Vd=V1-V4;
+T1=300;
+Tf=288;
+P1=0.96;
+Pf=1.01325;
+Vf=[Tf*P1*Vd]/[Pf*T1];
+Mcd=V1/(Vd);
+printf('Vf= %2.4f m^3/cycle',Vf);
+printf('\n');
+printf('Mc/Md= %2.2f ',Mcd);
+printf('\n');
+
+N=600;
+W=[n/(n-1)]*[P1*100*Vd]*[((P3/P1)^((n-1)/n))-1];
+IP=W*N/60;
+printf('Indicated Power= %2.2f kW',IP);
+printf('\n');
diff --git a/2339/CH9/EX9.1.1/Ex9_1.sce b/2339/CH9/EX9.1.1/Ex9_1.sce new file mode 100755 index 000000000..8afe6593f --- /dev/null +++ b/2339/CH9/EX9.1.1/Ex9_1.sce @@ -0,0 +1,33 @@ +clc +clear + +D=0.3; +L=0.6; +N=60; +Hs=5; +Hd=10; +Ep=0.8; +Qa=0.075; + +A=(22/7)*(1/4)*D*D; +Rho=1000; +g=9.81; + +F1=Rho*g*Hs*A; +F2=Rho*g*Hd*A; + +TF=F1+F2; +printf('Total Force Required: %2.2f kN',TF/1000); +printf('\n'); + +Q=(2*L*A*N)/60; +Qa=0.075; +Slip=(Q-Qa)/Q; +printf('Percentage Slip: %2.2f Percent',Slip*100); +printf('\n'); + +Cd=Qa/Q; + +P=(Rho*g*Qa*(Hs+Hd))/Ep; +printf('Power input: %2.2f kW',P/1000); +printf('\n'); diff --git a/2339/CH9/EX9.2.1/Ex9_2.sce b/2339/CH9/EX9.2.1/Ex9_2.sce new file mode 100755 index 000000000..21edafb40 --- /dev/null +++ b/2339/CH9/EX9.2.1/Ex9_2.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+Qa=0.025;
+Hm=20;
+L=0.4;
+D=0.3;
+
+A=(22/7)*(1/4)*D*D;
+Slip=0.02;
+Q=25/[1000*(1-Slip)];
+
+N=(Q*60)/(L*A);
+
+printf('Speed of Pump: %2.2f RPM',N);
+printf('\n');
diff --git a/2339/CH9/EX9.3.1/Ex9_3.sce b/2339/CH9/EX9.3.1/Ex9_3.sce new file mode 100755 index 000000000..ac7aa0cbc --- /dev/null +++ b/2339/CH9/EX9.3.1/Ex9_3.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+Hs=32;
+N=1450;
+Eff=0.85;
+Q=0.05;
+Hfs=1;
+Hfd=6;
+Hm=Hs+Hfd+Hfs;
+Rho=1000;
+g=9.81;
+
+P=[Rho*g*Q*Hm]/Eff;
+printf('Power Consumed: %2.2f kW',P/1000);
+printf('\n');
diff --git a/2339/CH9/EX9.4.1/Ex9_4.sce b/2339/CH9/EX9.4.1/Ex9_4.sce new file mode 100755 index 000000000..ae212d734 --- /dev/null +++ b/2339/CH9/EX9.4.1/Ex9_4.sce @@ -0,0 +1,18 @@ +clc
+clear
+
+Pm=25;
+Em=0.9;
+Q=0.063;
+Hs=4;
+Hd=25;
+Rho=1000;
+Hm=Hs+Hd;
+g=9.81;
+
+Ph=Rho*g*Q*Hm/1000;
+Ps=Em*Pm;
+Ep=Ph/Ps;
+
+printf('Efficiency of Pump: %2.2f Percent',Ep*100);
+printf('\n');
|