diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3685/CH9/EX9.13/Ex9_13.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3685/CH9/EX9.13/Ex9_13.sce')
| -rw-r--r-- | 3685/CH9/EX9.13/Ex9_13.sce | 45 |
1 files changed, 45 insertions, 0 deletions
diff --git a/3685/CH9/EX9.13/Ex9_13.sce b/3685/CH9/EX9.13/Ex9_13.sce new file mode 100644 index 000000000..7f5729a63 --- /dev/null +++ b/3685/CH9/EX9.13/Ex9_13.sce @@ -0,0 +1,45 @@ +clc
+x1 = 1 // Steam quality in first vessel
+x2 = 0.8 // Steam quality in second vessel
+// at 0.2MPa
+vg = 0.8857 // Specific volume of gas in m^3/kg
+h1 = 2706.7 // Enthalpy in first vessel in kJ/kg
+v1 = vg // Specific volume of gas in first vessel in m^3/kg
+hg = h1 // Enthalpy in first vessel 1 in kJ/kg
+m1 = 5 // mass in first vessel in kg
+V1 = m1*v1 // Volume of first vessel in m^3
+// at 0.5MPa
+m2 = 10 // mass in second vessel in kg
+hf = 640.23 // Enthalpy in second vessel in kJ/kg
+hfg = 2108.5 // Latent heat of vaporization in kJ/kg
+vf = 0.001093 // Specific volume of fluid in second vessel in m^3/kg
+vfg = 0.3749 // Change in specific volume in second vessel due to evaporation of gas in m^3/kg
+v2 = vf+(x2*vfg) // Specific volume of gas in second vessel
+V2 = m2*v2 // Volume of second vessel in m^3
+//
+Vm = V1+V2 // Total volume
+m = m1+m2 // Total mass
+vm = Vm/m // net specific volume
+u1 = h1 // Internal energy
+h2 = hf+(x2*hfg) // Enthalpy calculation
+u2 = h2 // Internal energy calculation
+m3 = m // Net mass calculation
+h3 = ((m1*u1)+(m2*u2))/m3 // Resultant enthalpy calculation
+u3 = h3 // Resultant internal energy calculation
+v3 = vm // resultant specific volume calculation
+// From Mollier diagram
+x3 = 0.870 // Steam quality
+p3 = 3.5 // Pressure in MPa
+s3 = 6.29 // Entropy at state 3 in kJ/kgK
+s1 = 7.1271 // Entropy at state 1 in kJ/kgK
+sf = 1.8607 // Entropy in liquid state in kJ/kgK
+sfg = 4.9606 // Entropy change due to vaporization in kJ/kgK
+s2 = sf+(x2*sfg) // Entropy calculation
+E = m3*s3-((m1*s1)+(m2*s2)) // Entropy change during process
+
+printf("\n Example 9.13")
+printf("\n Final pressure is %f bar",p3)
+printf("\n Steam quality is %f ",x3)
+printf("\n Entropy change during the process is %f kJ/K",E)
+//The answers vary due to round off error
+
|