8 files changed, 128 insertions, 0 deletions
diff --git a/1709/CH7/EX7.1/7_1.sce b/1709/CH7/EX7.1/7_1.sce
new file mode 100755
index 000000000..6343ae798
--- /dev/null
+++ b/1709/CH7/EX7.1/7_1.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+disp("Using gas tables,")
+T1=1160 //R
+h1=281.14 //B/lbm
+Pr1=21.18 
+P2=30 //psia
+P1=100 //psia
+//calculations
+Pr2=Pr1*P2/P1
+T2=833 //R
+h2=199.45 //B/lbm
+dh=h2-h1
+//results
+printf("Change in enthalpy = %.2f B/lbm",dh)
diff --git a/1709/CH7/EX7.2/7_2.sce b/1709/CH7/EX7.2/7_2.sce
new file mode 100755
index 000000000..ac4d93362
--- /dev/null
+++ b/1709/CH7/EX7.2/7_2.sce
@@ -0,0 +1,11 @@
+clc
+//Initialization of variables
+T2=860 //R
+phi1=0.78767
+phi2=0.71323
+P2=30 //psia
+P1=100 //psia
+//calculations
+dS=phi2-phi1- 53.35/778 *log(P2/P1)
+//results
+printf("Net change of entropy = %.5f B/lbm R",dS)
diff --git a/1709/CH7/EX7.3/7_3.sce b/1709/CH7/EX7.3/7_3.sce
new file mode 100755
index 000000000..267e1addb
--- /dev/null
+++ b/1709/CH7/EX7.3/7_3.sce
@@ -0,0 +1,13 @@
+clc
+//Initialization of variables
+T1=540 //R
+T2=960 //R
+disp("From gas tables,")
+h2=231.06 //B/lbm
+h1=129.06 //B/lbm
+cp=0.24
+//calculations
+W=h2-h1
+dh=cp*(T2-T1)
+//results
+printf("Change in enthalpy = %.1f B/lbm",dh)
diff --git a/1709/CH7/EX7.4/7_4.sce b/1709/CH7/EX7.4/7_4.sce
new file mode 100755
index 000000000..de9ee49b3
--- /dev/null
+++ b/1709/CH7/EX7.4/7_4.sce
@@ -0,0 +1,12 @@
+clc
+//Initialization of variables
+T1=420 //R
+T2=380 //R
+hig=1221.2
+P1=0.0019
+//calculations
+lnp=hig*778*(1/T1 - 1/T2)/85.6
+pra=exp(lnp)
+P2=pra*P1
+//results
+printf("Final pressure = %.3e psia",P2)
diff --git a/1709/CH7/EX7.5/7_5.sce b/1709/CH7/EX7.5/7_5.sce
new file mode 100755
index 000000000..e045d1a5d
--- /dev/null
+++ b/1709/CH7/EX7.5/7_5.sce
@@ -0,0 +1,17 @@
+clc
+//Initialization of variables
+disp("from critical constant tables")
+pc=482//psia
+Tc=227 //R
+vc=1.44 //ft^3/lbm mol
+P=600 //psia
+T=310 //R
+//calculations
+Pr=P/pc
+Tr=T/Tc
+disp("From Z tables,")
+Z=0.83
+v=Z*55.12*T/(P*144)
+rho=1/v
+//results
+printf("Density = %.1f lbm/ft^3",rho)
diff --git a/1709/CH7/EX7.6/7_6.sce b/1709/CH7/EX7.6/7_6.sce
new file mode 100755
index 000000000..e7e49cb63
--- /dev/null
+++ b/1709/CH7/EX7.6/7_6.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+T=-150+460 //R
+v=0.6 //ft^3/lbm
+vc=1.44
+Tc=227 //R
+Pc=482 //psia
+//calculations
+disp("From tables of z")
+vr=v/vc
+Tr=T/Tc
+Pr=1.75
+P=Pr*Pc
+//results
+printf("Final pressure = %d psia",P)
diff --git a/1709/CH7/EX7.7/7_7.sce b/1709/CH7/EX7.7/7_7.sce
new file mode 100755
index 000000000..60b59d465
--- /dev/null
+++ b/1709/CH7/EX7.7/7_7.sce
@@ -0,0 +1,21 @@
+clc
+//Initialization of variables
+disp("Critical tables suggest,")
+Tc=344 //R
+Pc=673 //psia
+P1=20 //psia
+P2=500 //psia
+M=16
+T=560 //R
+//calculations
+pr1=P1/Pc
+pr2=P2/Pc
+Tr=T/Tc
+dh2=0.65*Tc
+dsp=0.35 //B/lbm mol R
+dsp2=0.018-dsp- 1545/778 *log(P2/P1)
+W=dh2-dsp2*T
+W2=W/M
+//results
+printf("Work per pound mass = %d B/lbm",W2)
+//The answer is a bit different due to rounding off error
diff --git a/1709/CH7/EX7.8/7_8.sce b/1709/CH7/EX7.8/7_8.sce
new file mode 100755
index 000000000..5c48db17a
--- /dev/null
+++ b/1709/CH7/EX7.8/7_8.sce
@@ -0,0 +1,24 @@
+clc
+//Initialization of variables
+P=1000 //psia
+T1=100 + 460 //R
+T2=800+460 //R
+//calculations
+pc=1070 //psia
+Tc=548 //R
+pr1=P/pc
+Tr1=T1/Tc
+Tr2=T2/Tc
+M=44
+disp("from fig 7.7")
+h1=4235.8 //B/lbm mol
+h2=11661 //B/lbm mol
+h2bar=3.5 //B/lbm mol
+h1bar=0.48 //B/lbm mol
+dhbar=Tc*(h2bar-h1bar) + h2-h1
+Q=dhbar/M
+cp=0.202 //B/lbm F
+Q2=cp*(T2-T1)
+Error=(Q-Q2)/Q
+//results
+printf("Error in calculation = %d percent",Error*100)