9 files changed, 134 insertions, 0 deletions
diff --git a/1460/CH6/EX6.1/6_1.sce b/1460/CH6/EX6.1/6_1.sce
new file mode 100755
index 000000000..7646de828
--- /dev/null
+++ b/1460/CH6/EX6.1/6_1.sce
@@ -0,0 +1,11 @@
+clc
+//intialization of variables
+T1=40+460 //R
+T2=340+460 //R
+//calculations
+function [cv] = Cv(T)
+	cv=0.162+0.00046*T
+endfunction
+du=intg(T1,T2,Cv)
+//results
+printf("Change in specific internal energy = %.1f B/lbm",du)
diff --git a/1460/CH6/EX6.2/6_2.sce b/1460/CH6/EX6.2/6_2.sce
new file mode 100755
index 000000000..f6d52ec44
--- /dev/null
+++ b/1460/CH6/EX6.2/6_2.sce
@@ -0,0 +1,8 @@
+clc
+//Initialization of variables
+cp=0.24 //B/lbm F
+R=53.3 //ft-lb/lbm F
+//calculations
+cv=cp-R/778
+//results
+printf("Specific heat at constant volume = %.3f B/lbm F",cv)
diff --git a/1460/CH6/EX6.3/6_3.sce b/1460/CH6/EX6.3/6_3.sce
new file mode 100755
index 000000000..281c71126
--- /dev/null
+++ b/1460/CH6/EX6.3/6_3.sce
@@ -0,0 +1,11 @@
+clc
+//Initialization of variables
+T1=1400+460 //R
+T2=1200+460 //R
+//calculations
+function [cp] = Cp(T)
+	cp=0.317- 1.2*100/T +4*10^4 /T^2
+endfunction
+dh=intg(T1,T2,Cp)
+//results
+printf("Change in stagnation enthalpy = %.1f B/lbm",dh)
diff --git a/1460/CH6/EX6.4/6_4.sce b/1460/CH6/EX6.4/6_4.sce
new file mode 100755
index 000000000..c3ae4836b
--- /dev/null
+++ b/1460/CH6/EX6.4/6_4.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+T1=100+460 //R
+T2=300+460 //R
+P1=15 //lb/in^2
+P2=30 //lb/in^2
+Cp=0.3 //B/lbm F
+R=40 //ft-lb/lbm R
+//calculations
+function [s] = fun(f)
+	s=Cp/f
+endfunction
+function [s2] = fun1(f)
+	s2=R/(f*778)
+endfunction
+ds=intg(T1,T2,fun) - intg(P1,P2,fun1)
+//results
+printf("Change in entropy = %.4f B/lbm R",ds)
diff --git a/1460/CH6/EX6.5/6_5.sce b/1460/CH6/EX6.5/6_5.sce
new file mode 100755
index 000000000..45136c8ae
--- /dev/null
+++ b/1460/CH6/EX6.5/6_5.sce
@@ -0,0 +1,12 @@
+clc
+//Initialization of variables
+T1=40+460 //R
+T2=340+460 //R
+P1=15 //lb/in^2
+cp=0.24
+cv=0.171
+//calculations
+gamma=cp/cv
+P2=P1 *(T2/T1)^(gamma/(gamma-1))
+//results
+printf("Final pressure = %.1f lb/in^2",P2)
diff --git a/1460/CH6/EX6.6/6_6.sce b/1460/CH6/EX6.6/6_6.sce
new file mode 100755
index 000000000..b402c9e04
--- /dev/null
+++ b/1460/CH6/EX6.6/6_6.sce
@@ -0,0 +1,16 @@
+clc
+//Initialization of variables
+P1=16 //lb/in^2
+P2=14 //lb/in^2
+Tt=83+460 //R
+gamma=1.4
+cp=0.24 //B/lbm F
+//calculations
+T=Tt *(P2/P1)^((gamma-1)/gamma)
+dh=cp*(Tt-T)
+V=sqrt(2*32.2*778*dh)
+//results
+printf("Actual temperature in the flow = %d R",T)
+printf("\n Flow velocity = %d ft/sec",V)
+
+
diff --git a/1460/CH6/EX6.7/6_7.sce b/1460/CH6/EX6.7/6_7.sce
new file mode 100755
index 000000000..ccdbaeca6
--- /dev/null
+++ b/1460/CH6/EX6.7/6_7.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+T1=400+460 //R
+P1=100 //lb/in^2
+P2=20 //lb/in^2
+T2=140+460 //R
+Cp=50
+//calculations
+Pratio=P1/P2
+Tratio=T1/T2
+C=log(Tratio) /log(Pratio)
+n=1/(1-C)
+v1=Cp*T1/(144*P1)
+v2=Cp*T2/(144*P2)
+w=144*P1*v1^n
+function[p]=fun(v)
+	p=w/v^n
+endfunction
+Work=intg(v1,v2,fun)
+//results
+printf("Work done = %.d ft-lb/lbm",Work)
+//The answers in the textbook varies a bit due to rounding off errors
diff --git a/1460/CH6/EX6.8/6_8.sce b/1460/CH6/EX6.8/6_8.sce
new file mode 100755
index 000000000..95b9256ae
--- /dev/null
+++ b/1460/CH6/EX6.8/6_8.sce
@@ -0,0 +1,13 @@
+clc
+//Initialization of variables
+P1=15 //lb/in^2
+P2=20 //lb/in^2
+T1=40+460 //R
+T2=540+460 //R
+//calculations
+disp("From table 6 at the two temperatures")
+phi1=0.58233
+phi2=0.75042
+ds=phi2-phi1-53.3*log(P2/P1) /778
+//results
+printf("Entropy change = %.5f B/lbm R",ds)
diff --git a/1460/CH6/EX6.9/6_9.sce b/1460/CH6/EX6.9/6_9.sce
new file mode 100755
index 000000000..bc0a15e80
--- /dev/null
+++ b/1460/CH6/EX6.9/6_9.sce
@@ -0,0 +1,23 @@
+clc
+//Initialization of variables
+T1=1440+460 //R
+T2=1000+460 //R
+n=1.4 
+//calculations
+Pratio=(T2/T1)^(n/(n-1))
+Vratio=(T1/T2)^(1/(n-1))
+disp("From table 6")
+Pr1=141.51
+Pr2=50.34
+vr1=4.974
+vr2=10.743
+Pratio2=Pr2/Pr1
+Vratio2=vr2/vr1
+//results
+//The answer in the textbook given for Vratio is wrong.
+printf("Case 1")
+printf("\n Pressure ratio = %.1f",Pratio+0.1)
+printf("\n Volume ratio = %.2f",Vratio)
+printf("\n Case 2")
+printf("\n Pressure ratio = %.3f",Pratio2)
+printf("\n Volume ratio = %.2f",Vratio2)