initial commit / add all books

10 files changed, 158 insertions, 0 deletions

diff --git a/1460/CH11/EX11.1/11_1.sce b/1460/CH11/EX11.1/11_1.sce
new file mode 100755
index 000000000..fc82c6ff5
--- /dev/null
+++ b/1460/CH11/EX11.1/11_1.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of variables

+P1=100 //lb/in^2

+P2=14.7 //lb/in^2

+T1=600+460 //R

+T2=300+460 //R

+area=1 //in^2

+//calculations

+disp("From steam tables")

+ht1=1329.1 

+h2=1192.8

+v2=30.53

+Vel=sqrt(2*32.2*778*(ht1-h2))

+wdot=area*Vel/(144*v2)

+//results

+printf("Exit velocity = %d ft/sec",Vel)

+printf("\n Mass flow rate = %.3f lbm/sec",wdot)

diff --git a/1460/CH11/EX11.10/11_10.sce b/1460/CH11/EX11.10/11_10.sce
new file mode 100755
index 000000000..8e70ca44c
--- /dev/null
+++ b/1460/CH11/EX11.10/11_10.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of variables

+P=300 //lb/in^2

+P2=14.7 //lb/in^2

+T=4540+460 //R

+w=100 //lbm/sec

+g=1.25

+MW=30

+R=1544

+//calculations

+R=R/MW

+Pratio=P2/P

+V4=sqrt(2*g*32.2*51.5*T*(1-(Pratio)^((g-1)/g))/(g-1)) 

+Fn=w*V4/32.2

+//results

+printf("Thrust = %d lb",Fn)

+//The answer in the textbook is a bit different due to rounding off error.

diff --git a/1460/CH11/EX11.2/11_2.sce b/1460/CH11/EX11.2/11_2.sce
new file mode 100755
index 000000000..c3cb5deee
--- /dev/null
+++ b/1460/CH11/EX11.2/11_2.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of variables

+Pt1=100 //lb/in^2

+P2=15 //lb/in^2

+A=1 //in^2

+T=500+460 //F

+gamma=1.4

+//calculations

+Pratio=P2/Pt1

+r1=(P2/Pt1)^((gamma-1)/gamma)

+r2=(P2/Pt1)^(2/gamma)

+r3=(P2/Pt1)^((gamma+1)/gamma)

+V2=sqrt(2*gamma*32.2*53.3*T*(1-r1)/(gamma-1))

+wdot=A*Pt1*sqrt(2*gamma*(r2-r3)/(gamma-1)) /(sqrt(53.3*T/32.2))

+//results

+printf("Exit velocity = %d ft/sec",V2)

+printf("\n Mass flow rate = %.2f lbm/sec",wdot)

diff --git a/1460/CH11/EX11.3/11_3.sce b/1460/CH11/EX11.3/11_3.sce
new file mode 100755
index 000000000..5d3bcecb9
--- /dev/null
+++ b/1460/CH11/EX11.3/11_3.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of variables

+Pt1=100 //lb/in^2

+Tt1=960 //RP2=15 //lb/in^2

+wdot=1.13 //lbm/sec

+gamma=1.4

+//calculations

+Pstar=Pt1*(2/(1+gamma))^(gamma/(gamma-1))

+Tstar=Tt1*(2/(1+gamma))

+Vstar=sqrt(gamma*32.2*53.3*Tstar)

+vstar=53.3*Tstar/(144*Pstar)

+Astar=wdot*vstar*144/Vstar

+//results

+printf("Ideal throat area = %.3f in^2",Astar)

+printf("\n Ideal pressure = %.1f lb/in^2",Pstar)

+printf("\n Ideal temperature = %d R",Tstar)

+printf("\n Ideal throat specific volume = %.1f ft^3/lbm",vstar)

diff --git a/1460/CH11/EX11.4/11_4.sce b/1460/CH11/EX11.4/11_4.sce
new file mode 100755
index 000000000..29c1e716d
--- /dev/null
+++ b/1460/CH11/EX11.4/11_4.sce
@@ -0,0 +1,10 @@
+clc

+//initialization of variables

+ht1=1329.1

+st1=1.7581

+h2s=1151.4

+s2s=1.7581

+//calculations

+eta=sqrt((ht1-1192.8)/(ht1-h2s))

+//results

+printf("\n efficiency = %.2f",eta)

diff --git a/1460/CH11/EX11.5/11_5.sce b/1460/CH11/EX11.5/11_5.sce
new file mode 100755
index 000000000..bc3e2eee0
--- /dev/null
+++ b/1460/CH11/EX11.5/11_5.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of variables

+P1=100 //lb/in^2

+T1=500+460 //R

+P2=15 //lb/in^2

+eta=0.95

+A=1 //in^2

+gamma=1.4

+//calculations

+Ve=2200 //ft/sec

+V2=eta*Ve

+T2=T1*(1-eta*(1-(P2/P1)^((gamma-1)/gamma)))

+vexit=53.3*T2/(144*P2)

+wdot=A*V2/(144*vexit)

+//results

+printf("Exit velocity = %.1f ft^3/lbm",vexit)

+printf("\n Mass flow = %.3f lbm/sec",wdot)

diff --git a/1460/CH11/EX11.6/11_6.sce b/1460/CH11/EX11.6/11_6.sce
new file mode 100755
index 000000000..b86495d01
--- /dev/null
+++ b/1460/CH11/EX11.6/11_6.sce
@@ -0,0 +1,14 @@
+clc

+//initialization of variables

+v=500 //ft/sec

+P=14.7 //lb/in^2

+T=60+460 //R

+eta=0.85

+cp=0.24

+gamma=1.4

+//calculations

+Pt2=eta*P*(1+ (gamma-1)*v^2 /(2*gamma*32.2*53.3*T))^(gamma/(gamma-1))

+Tratio=1+ (gamma-1)*v*v/(2*gamma*32.2*53.3*T)

+Tt2=T*Tratio

+//results

+printf("Exit stagnation temperature  = %d R",Tt2+1)

diff --git a/1460/CH11/EX11.7/11_7.sce b/1460/CH11/EX11.7/11_7.sce
new file mode 100755
index 000000000..c6c6aa050
--- /dev/null
+++ b/1460/CH11/EX11.7/11_7.sce
@@ -0,0 +1,16 @@
+clc

+//initialization of variables

+P=30 //lb/in^2

+T=1000+460 //R

+Pd=14.7 //lb/in^2

+w=60 //lbm/sec

+eta=0.95 //percent

+R=53.3

+gamma=1.35

+cp=0.264

+//calculations

+V2s=sqrt(2*gamma*32.2*53.3*T*(1-(Pd/P)^(0.259))/(gamma-1))

+V2=eta*V2s

+Fn=w*(V2)/32.2

+//results

+printf("Thrust of the engine = %d ft/sec",Fn)

diff --git a/1460/CH11/EX11.8/11_8.sce b/1460/CH11/EX11.8/11_8.sce
new file mode 100755
index 000000000..8a38f16ae
--- /dev/null
+++ b/1460/CH11/EX11.8/11_8.sce
@@ -0,0 +1,21 @@
+clc

+//initialization of variables

+v=600 //ft/sec

+T=60+460 //R

+P=14.7 //lb/in^2

+Pratio=6

+Tin=1540+460 //R

+cp=0.264

+cpratio=1.35

+//calculations

+Pt2byP1=(1+ (cpratio-1)*v^2 /(cpratio*2*32.2*53.3*T))^(3.86)

+Pt3byP1=Pt2byP1*Pratio

+eta=1- 1/(Pt3byP1)^0.259

+Tt3=T*(Pt3byP1)^((cpratio-1)/cpratio)

+Q=cp*(Tin-Tt3)

+V6=sqrt(eta*2*32.2*778*Q + v^2)

+Fn=(V6-v)/32.2

+//resullts

+printf("Thermal efficiency = %.2f ",eta)

+printf("\n thrust per pound of air per sec = %.1f lb-sec/lbm",Fn)

+//The answers are a bit different due to rounding off error in textbook

diff --git a/1460/CH11/EX11.9/11_9.sce b/1460/CH11/EX11.9/11_9.sce
new file mode 100755
index 000000000..c2a1995f9
--- /dev/null
+++ b/1460/CH11/EX11.9/11_9.sce
@@ -0,0 +1,12 @@
+clc

+//initialization of variables

+V=1000 //mph

+P=14.7 //lb/in^2

+T=60 //F

+g=1.4

+//calculations

+V1=V*(88/T)

+Pratio=(1+ (g-1)*V1^2 /(2*g*32.2*53.3*(T+460)))^(g/(g-1))

+eta=1-1/(Pratio)^0.286

+//results

+printf("Theoretical cycle efficiency = %.3f",eta)