From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2006/CH3/EX3.1/ex3_1.sce | 22 ++++++++++++++++++++++
 2006/CH3/EX3.2/ex3_2.sce | 26 ++++++++++++++++++++++++++
 2006/CH3/EX3.3/ex3_3.sce |  7 +++++++
 2006/CH3/EX3.4/ex3_4.sce | 10 ++++++++++
 2006/CH3/EX3.5/ex3_5.sce | 12 ++++++++++++
 2006/CH3/EX3.7/ex3_7.sce | 33 +++++++++++++++++++++++++++++++++
 2006/CH3/EX3.8/ex3_8.sce | 15 +++++++++++++++
 7 files changed, 125 insertions(+)
 create mode 100755 2006/CH3/EX3.1/ex3_1.sce
 create mode 100755 2006/CH3/EX3.2/ex3_2.sce
 create mode 100755 2006/CH3/EX3.3/ex3_3.sce
 create mode 100755 2006/CH3/EX3.4/ex3_4.sce
 create mode 100755 2006/CH3/EX3.5/ex3_5.sce
 create mode 100755 2006/CH3/EX3.7/ex3_7.sce
 create mode 100755 2006/CH3/EX3.8/ex3_8.sce

(limited to '2006/CH3')

diff --git a/2006/CH3/EX3.1/ex3_1.sce b/2006/CH3/EX3.1/ex3_1.sce
new file mode 100755
index 000000000..d13b237f7
--- /dev/null
+++ b/2006/CH3/EX3.1/ex3_1.sce
@@ -0,0 +1,22 @@
+clc;
+V=0.01; // Volume of water in a rigid vessel in m^3
+m=4.5; // Mass of water+ steam in a rigid vessel in kg
+T=35; // Temperature of water in a rigid vessel in degree celcius
+// (a)      
+v=V/m; // specific volume of water
+// From steam table
+vf=0.001006; vg=25.22; // specific volume in m^3/kg
+x=(v-vf)/(vg-vf); // Quality of steam
+x1=1-x; // Quality of water
+mg=x*m; // Mass of steam
+mf=x1*m; // Mass of water
+disp ("kg",mf,"Mass of water in a rigid vessel = ","kg",mg,"Mass of steam in a rigid vessel = ",x1,"Quality of water in a rigid vessel = ",x,"Quality of steam in a rigid vessel = ","  (a)   ");
+// (b)          
+vc=0.003155; //  Crictical volume for water in m^3/kg
+disp ("The level of liquid water will rise in the vessel. Since v < vc and refer figure 3.21","  (b)   ");
+// (c)  
+disp ("The final temperature after heating is 370.04 oC.  Because it is constant volume process and refer figure 3.21","  (c)  ");
+//  (d)   
+m1=0.45; // Mass of water in kg
+v1=V/m; // specific volume of water
+disp ("Level of liquid drops to bottom (v1 > vc).  Temperature on reaching saturation state is 298.5 oC and refer figure 3.21", "  (d)  ");
diff --git a/2006/CH3/EX3.2/ex3_2.sce b/2006/CH3/EX3.2/ex3_2.sce
new file mode 100755
index 000000000..d33501fba
--- /dev/null
+++ b/2006/CH3/EX3.2/ex3_2.sce
@@ -0,0 +1,26 @@
+clc;
+ // (a) Ammonia 26 oC and 0.074 m^3/kg
+// From saturation table of ammonia at 26 oC
+v=0.074; // specific volume of ammonia in m^3/kg
+vf=0.001663; vg=0.1245; // specific volume  of ammonia in m^3/kg
+x=(v-vf)/(vg-vf); // Quality of vapour  since v<vg
+disp (x,"The Quality of ammonia = ","(a) Ammonia 26 oC and 0.074 m^3/kg");
+// (b).Ammonia 550kPa and 0.31m^3/kg
+// From saturation table of ammonia at 550 kPa
+v=0.31; // specific volume of ammonia in m^3/kg
+vg=0.23; // specific volume  of ammonia in m^3/kg
+// v > vg . Since from superheated table by interpolation for 550kPa and v
+T=82.1; // Temperature of ammonia in degree celcius
+disp ("oC",T,"Temperature of ammonia = ","(b).Ammonia 550kPa and 0.31m^3/kg");
+// (c).Freon 12, 0.35MPa and 0.036 m^3/kg
+// From saturation table of Freon 12 at 0.35MPa
+v=0.036; // specific volume of Freon 12 in m^3/kg
+vf=0.000722; vg=0.049329; // specific volume  of Freon 12 in m^3/kg
+x=(v-vf)/(vg-vf); // Quality of vapour since v<vg
+disp (x,"The Quality of Freon 12  = ","(c).Freon 12, 0.35MPa and 0.036 m^3/kg");
+// (d).Methane 0.5MPa and 1.0 m^3/kmol
+v=1; // specific volume of Methane in m^3/kmol
+// From table at 0.5 MPa molar values are
+vf=0.04153; vg=2.007; // specific volume  of Methane in m^3/kmol
+x=(v-vf)/(vg-vf); // Quality of vapour since v<vg
+disp (x,"The Quality of Methane  = ","(d).Methane 0.5MPa and 1.0 m^3/kmol");
diff --git a/2006/CH3/EX3.3/ex3_3.sce b/2006/CH3/EX3.3/ex3_3.sce
new file mode 100755
index 000000000..39c305578
--- /dev/null
+++ b/2006/CH3/EX3.3/ex3_3.sce
@@ -0,0 +1,7 @@
+clc;
+V=300; // Volume of air in the room in m^3
+p=1; // Atmospheric pressure in bar
+T=25; // Temperature of air in Degree Celcius
+R=287; // Characteristic constant of Air in J/kg k
+m=(p*10^5*V)/(R*(T+273)); // Ideal gas equation
+disp ("kg",m,"Mass of air in room");
diff --git a/2006/CH3/EX3.4/ex3_4.sce b/2006/CH3/EX3.4/ex3_4.sce
new file mode 100755
index 000000000..aa9e3e181
--- /dev/null
+++ b/2006/CH3/EX3.4/ex3_4.sce
@@ -0,0 +1,10 @@
+clc;
+D=20; // Diameter of the sphere in cm
+m=2.54; // Mass of gas filled in sphere in gram
+p=10; // Pressure of gas in bar
+T=25; // Temperature of gas in Degree Celcius
+R=8.3144*10^3; // Universal gas constant in J/kmol K
+V=(3.14*(D*10^-2)^3)/16; // Volume of das in sphere in m^3
+M=(m*10^-3*R*(T+273))/(p*10^5*V); // Molecular weight of the gas
+disp (M,"Molecular weight of the gas");
+disp ("Therefore gas in sphere is Helium (unless mixture of two or more gases)");
diff --git a/2006/CH3/EX3.5/ex3_5.sce b/2006/CH3/EX3.5/ex3_5.sce
new file mode 100755
index 000000000..199a4afb1
--- /dev/null
+++ b/2006/CH3/EX3.5/ex3_5.sce
@@ -0,0 +1,12 @@
+clc;
+p2=2.5; // Pressure of air in the cylinder in bar
+T1=430; // Temperature of air in cylinder in Degree Celcius
+V1=1.2; // Volume of cylinder in m^3
+V2=0.6; // Volume of cylinder upto end stops in m^3
+// (a) Temperature of air when the piston reaches the stops
+T2=(T1+273)*(V2/V1); // constant pressure process
+disp ("K",T2,"Temperature of air when the piston reaches the stops");
+// (b) The pressure of air when its temperature equals to 25 oC
+T3=25; //Room temperature in Degree Celcius
+p3=p2*((T3+273)/T2); // constant volume process
+disp ("bar",p3,"The pressure of air when its temperature equals to 25 oC");
diff --git a/2006/CH3/EX3.7/ex3_7.sce b/2006/CH3/EX3.7/ex3_7.sce
new file mode 100755
index 000000000..e385cc9f8
--- /dev/null
+++ b/2006/CH3/EX3.7/ex3_7.sce
@@ -0,0 +1,33 @@
+clc;
+p=6000; // Pressure of nitrogen gas in kPa
+T=150; // Temperature of nitrogen gas in kelvin
+V=250; // Volume of tank in litres
+R_1=8.3143; // Universal gas constant in kJ/kmol K
+M=28.1013; // Molecular mass
+// (a).Beattie - Bridgeman equation of state
+// Constants for nitrogen gas
+c=4.2*10^4; Ao=136.2315; a=0.02617; Bo=0.05046; b=-0.00691;
+// By substituting these values in the following equation 
+// p=(R_1*T/v^2)*(1-(c/(vT^3)))*(v+Bo*(1-(b/v)))-(Ao/v^2*(1-(a/v)))
+// By trial and error we get
+v=0.1222; // specific volume in m^3/kmol
+m=(M*V/1000)/v; // Mass of nitrogen gas
+disp ("m^3/kmol",v,"specific volume of nitrogen gas = ","kg",m,"Mass of nitrogen gas = ","(a).Beattie - Bridgeman equation of state");
+// (b).Nitrogen tables
+// From property table of nitrogen fas
+v=0.004413;  // specific volume in m^3/kg
+m=(V/1000)/v; //  Mass of nitrogen gas
+disp ("m^3/kg",v,"specific volume of nitrogen gas = ","kg",m,"Mass of nitrogen gas = ","(b).Nitrogen tables");
+// (c).Ideal gas equation of state
+m=(p*V/1000)/(R_1*T/M); //Mass of nitrogen gas
+disp ("kg",m,"Mass of nitrogen gas = ","(c).Ideal gas equation of state");
+// (d).Generalized compressibility chart
+// The crictical properties for nitrogen gas 
+Tc=126.2; // Temperature in kelvin
+Pc=3.349; // Pressure in MPa
+// Reduced properties are
+Pr=p/Pc; Tr=T/Tc;
+z=0.6; // From chart
+m=(p*V/1000)/(z*R_1*T/M); //Mass of nitrogen gas
+disp ("kg",m,"Mass of nitrogen gas = ","(d).Generalized compressibility chart");
+disp ("Ideal gas equation of state","Generalized compressibility chart","Beattie - Bridgeman equation of state","Nitrogen tables"," (e).Arrangement the methods in order of percentage error : ");
diff --git a/2006/CH3/EX3.8/ex3_8.sce b/2006/CH3/EX3.8/ex3_8.sce
new file mode 100755
index 000000000..34e2fff9f
--- /dev/null
+++ b/2006/CH3/EX3.8/ex3_8.sce
@@ -0,0 +1,15 @@
+clc;
+T=-58.7; //Normal boling point of CF3Br in Degree Celcius
+Tc=340.9; // Crictical temperature of CF3Br in K
+pc=4.05; // Crictical pressure of CF3Br in MPa
+M=148.9; // Moleclar mass of CF3Br
+p=1.01325*10^5; // Atmospheric pressure in N/m^2
+R1=8314.4; // Universal gas constant in J/kmol K
+R=R1/M; // Gas constant of CF3Br
+a=(0.42748*R^2*Tc^2.5)/(pc*10^6); // Constant of Redlich-Kwong equation of state
+b=(0.08664*R*Tc)/(pc*10^6); // Constant of Redlich-Kwong equation of state
+vi=(R*(T+273))/p; // Ideal gas volume for assigning initial value
+// By substituting these values in the Redlich-Kwong equation of state 
+ // vi_1=(R*(T+273)/p)+b-((a/(p*(273+T)^0.5*vi))) and and solving it by trial and error method we get
+vi_1=0.11443; // in m^3/kg
+disp ("m^3/kg",vi_1,"Saturated vapour volume");
-- 
cgit