6 files changed, 112 insertions, 0 deletions
diff --git a/1709/CH3/EX3.1/3_1.sce b/1709/CH3/EX3.1/3_1.sce
new file mode 100755
index 000000000..cfe39b906
--- /dev/null
+++ b/1709/CH3/EX3.1/3_1.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+V=1 //ft^3
+m=30 //lbm
+//calculations
+v=V/m
+vf1=0.01665
+vfg1=32.38 //ft^3/lbm
+x1=0.000515
+uf1=169.92
+ufg1=904.8
+u1=uf1+x1*ufg1
+vfg=0.0216
+vfg2=0.4240
+v2=v
+x2=0.0277
+uf2=538.4
+ufg2=571
+u2=uf2+x2*ufg2
+Q=m*(u2-u1)
+//results
+printf("Heat transfer = %d Btu",Q)
diff --git a/1709/CH3/EX3.2/3_2.sce b/1709/CH3/EX3.2/3_2.sce
new file mode 100755
index 000000000..35704d586
--- /dev/null
+++ b/1709/CH3/EX3.2/3_2.sce
@@ -0,0 +1,23 @@
+clc
+//Initialization of variables
+V2=2.5 //ft^3
+V1=0.5 //ft^3
+P=100 //psia
+x1=0.5
+//calculations
+W=-P*144*(V2-V1)
+vf1=0.01774
+vfg1=4.414
+v1=vf1+x1*vfg1
+m=V1/v1
+v2=V2/m
+disp("From tables ,")
+uf1=298.08
+ufg1=807.1
+u1=uf1+x1*ufg1
+h2=1747.9
+u2=h2-P*144*v2/778
+Q=m*(u2-u1)
+//results
+printf("Amount of heat = %d Btu",Q)
+//The answer for u2 is given wrong in the textbook. Please use a calculator to find it
diff --git a/1709/CH3/EX3.3/3_3.sce b/1709/CH3/EX3.3/3_3.sce
new file mode 100755
index 000000000..b8ab8b1ce
--- /dev/null
+++ b/1709/CH3/EX3.3/3_3.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+V1=1.735*10^-4 //ft^3
+v1=0.016080 //ft^3/lbm
+h1=70.61 //B/lbm
+P1=100 //psia
+V2=1 //ft^3
+//calculations
+u1=h1-P1*v1*144/778
+m=V1/v1
+v2=V2/m
+vf2=0.01613
+vfg2=350.3
+x2=(v2-vf2)/vfg2
+hf2=67.97
+hfg2=1037.2
+h2=hf2+x2*hfg2
+P2=0.9492
+u2=h2- P2*144*v2/778
+Q=m*(u2-u1)
+//results
+printf("Enthalpy change = %.2f Btu",Q)
diff --git a/1709/CH3/EX3.4/3_4.sce b/1709/CH3/EX3.4/3_4.sce
new file mode 100755
index 000000000..977971fdb
--- /dev/null
+++ b/1709/CH3/EX3.4/3_4.sce
@@ -0,0 +1,13 @@
+clc
+//Initialization of variables
+P=20 //psia
+V=1 //ft^3
+T=560 //R
+cv=0.1715
+Q=10//Btu
+//calculations
+m=P*144*V/(53.35*T)
+T2=Q/(m*cv) +T
+P2=m*53.35*T2/V
+//results
+printf("Fina pressure = %d lbf/ft^2",P2)
diff --git a/1709/CH3/EX3.5/3_5.sce b/1709/CH3/EX3.5/3_5.sce
new file mode 100755
index 000000000..fc0294ddd
--- /dev/null
+++ b/1709/CH3/EX3.5/3_5.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+T1=560 //R
+T2=3460  //R
+m=28.02 //lb
+cv=0.248
+//calculations
+function [q]=fun(T)
+    q=9.47 - 3.29*10^3 /T +1.07*10^6 /T^2
+endfunction
+Q1=intg(T1,T2,fun)
+Q2=m*cv*(T2-T1)
+Error=(Q1-Q2)/Q1
+//results
+printf("Percentage error = %.1f percent",Error*100)
diff --git a/1709/CH3/EX3.6/3_6.sce b/1709/CH3/EX3.6/3_6.sce
new file mode 100755
index 000000000..a1a7c6635
--- /dev/null
+++ b/1709/CH3/EX3.6/3_6.sce
@@ -0,0 +1,17 @@
+clc
+//Initialization of variables
+rate=20 //gal/min
+P1=20 //psia
+P2=1000 //psia
+T=100+460 //R
+//calculations
+vf=0.01613
+disp("From table A-8")
+dv=-5.2*10^-5 //ft^3/lbm
+K=-dv/(vf*P2*144)
+wt=K*vf*(P2^2 - P1^2)*144*144*10^4 /2
+m=rate*8.33
+Wt=wt*m
+Wthp=Wt/33000
+//results
+printf("Pump power required = %d hp",Wthp)