initial commit / add all books

1 files changed, 35 insertions, 0 deletions

diff --git a/3685/CH9/EX9.5/Ex9_5.sce b/3685/CH9/EX9.5/Ex9_5.sce
new file mode 100644
index 000000000..4daefea05
--- /dev/null
+++ b/3685/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,35 @@
+clc

+Psat = 3.973 // Saturation pressure in MPa

+vf = 0.0012512 // specific volume of fluid in m^3/kg

+vg = 0.05013 // Specific volume of gas in m^3/kg

+hf = 1085.36  // Specific enthalpy of fluid in kJ/kg

+hfg = 1716.2 // Latent heat of vaporization in kJ/kg

+sf = 2.7927 // Specific entropy of fluid in kJ/kgK

+sfg = 3.2802 // Entropy change due to vaporization in kJ/kgK

+mf = 9 // Mass of liquid in kg

+V = 0.04 // Volume of vessel in m^3

+// at T = 250

+uf = 1080.39 //Specific internal energy in kJ/kg 

+ufg = 1522// Change in internal energy due to vaporization in kJ/kg

+

+printf("\n Example 9.5")

+Vf = mf*vf // volume of fluid

+Vg = V-Vf // volume of gas

+mg = Vg/vg // mass of gas

+m = mf+mg // mass if mixture

+x = mg/m // quality of steam

+v = vf+x*(vg-vf) // specific volume of mixture

+h = hf+x*hfg // enthalpy of mixture

+s = sf+(x*sfg) // entropy of mixture

+u = h-Psat*1e6*v*1e-03 // Internal energy of mixture

+u_ = uf+x*ufg // Internal energy at 250 degree Celsius

+printf("\n The pressure is %f MPa",Psat)

+printf("\n The total mass of mixture is %f kg",m)

+printf("\n Specific volume is %f m3/kg",v)

+printf("\n Enthalpy is is %f kJ/kg",h)

+printf("\n The entropy is %f kJ/kg K",s)

+printf("\n The internal energy is %f kJ/kg",u)

+printf("\n At 250 degree Celsius, internal energy is %fkJ/kg",u_) //The answer provided in the textbook is wrong

+

+//The answers vary due to round off error

+