summaryrefslogtreecommitdiff
path: root/2024/CH9
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /2024/CH9
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '2024/CH9')
-rwxr-xr-x2024/CH9/EX9.1/9_1.sce22
-rwxr-xr-x2024/CH9/EX9.10/9_10.sce16
-rwxr-xr-x2024/CH9/EX9.11/9_11.sce14
-rwxr-xr-x2024/CH9/EX9.12/9_12.sce20
-rwxr-xr-x2024/CH9/EX9.13/9_13.sce19
-rwxr-xr-x2024/CH9/EX9.14/9_14.sce12
-rwxr-xr-x2024/CH9/EX9.15/9_15.sce18
-rwxr-xr-x2024/CH9/EX9.16/9_16.sce21
-rwxr-xr-x2024/CH9/EX9.17/9_17.sce13
-rwxr-xr-x2024/CH9/EX9.2/9_2.sce24
-rwxr-xr-x2024/CH9/EX9.3/9_3.sce18
-rwxr-xr-x2024/CH9/EX9.4/9_4.sce26
-rwxr-xr-x2024/CH9/EX9.5/9_5.sce34
-rwxr-xr-x2024/CH9/EX9.6/9_6.sce10
-rwxr-xr-x2024/CH9/EX9.8/9_8.sce17
-rwxr-xr-x2024/CH9/EX9.9/9_9.sce23
16 files changed, 307 insertions, 0 deletions
diff --git a/2024/CH9/EX9.1/9_1.sce b/2024/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..64adf4c4c
--- /dev/null
+++ b/2024/CH9/EX9.1/9_1.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+h1=1329.1 //Btu/lbm
+v1=6.218 //ft^3/lbm
+J=778
+g=32.174
+m=1
+//calculations
+p=[80 60 54.6 40 20]
+h=[ 1304.1 1273.8 1265 1234.2 1174.8]
+v=[ 7.384 9.208 9.844 12.554 21.279]
+Fc=1
+V2=Fc*sqrt(2*J*g*(h1-h))
+A=m*v ./V2
+V2=[0 V2]
+A=[0 A]
+//results
+disp('velocity = ')
+disp(V2 )
+disp('Area = ')
+disp(A)
+//The initial values of velocity and area are 0 and infinity respectively
diff --git a/2024/CH9/EX9.10/9_10.sce b/2024/CH9/EX9.10/9_10.sce
new file mode 100755
index 000000000..938d620ee
--- /dev/null
+++ b/2024/CH9/EX9.10/9_10.sce
@@ -0,0 +1,16 @@
+clc
+//Initialization of variables
+K=0.6007
+Y1=0.91
+D1=6.065
+D2=1.820
+rho1=0.156
+p1=30
+p2=20.18
+//calculations
+bet=D2/D1
+m=0.525*K*Y1^2 *D2^2 *sqrt(rho1*(p1-p2))
+C=K*sqrt(1-bet^4)
+//results
+printf("mass flow rate = %.2f lbm/sec",m)
+printf("\n Coefficient of discharge = %.3f",C)
diff --git a/2024/CH9/EX9.11/9_11.sce b/2024/CH9/EX9.11/9_11.sce
new file mode 100755
index 000000000..9100238cc
--- /dev/null
+++ b/2024/CH9/EX9.11/9_11.sce
@@ -0,0 +1,14 @@
+clc
+//Initialization of variables
+zm=0.216
+pm=62.3 //lbm/ft^2
+p1=0.0736 //lbm/ft^2
+g=32.2
+d=4
+//calculations
+H=zm*(pm-p1)/12/p1
+V=sqrt(2*g*H)
+m=%pi/4 *d^2 *V*p1
+//results
+printf("average velocity = %.1f ft/sec",V)
+printf("\n mass flow rate = %.1f lbm/sec",m)
diff --git a/2024/CH9/EX9.12/9_12.sce b/2024/CH9/EX9.12/9_12.sce
new file mode 100755
index 000000000..334a3419e
--- /dev/null
+++ b/2024/CH9/EX9.12/9_12.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+p1=50 //psia
+pr=0.58
+//calculations
+p=p1*pr
+s1=1.6585
+h1=1174.1 //Btu/lbm
+sf=0.3680
+sfg=1.3313
+hfg=945.3
+vg=13.746
+hf=218.82
+x= (s1-sf)/sfg
+v2=vg*x
+h2=hf+x*hfg
+V2rev=223.77*sqrt(h1-h2)
+m=%pi/4 *1/144 *V2rev/v2
+//results
+printf("mass flow rate = %.3f lbm/sec",m)
diff --git a/2024/CH9/EX9.13/9_13.sce b/2024/CH9/EX9.13/9_13.sce
new file mode 100755
index 000000000..550bd3cc0
--- /dev/null
+++ b/2024/CH9/EX9.13/9_13.sce
@@ -0,0 +1,19 @@
+clc
+//Initialization of variables
+k=1.31
+p1=7200 //lbf/ft^2
+v1=8.515 //ft^3/lbm
+pr=0.6
+m1=0.574
+T1=741 //R
+//calculations
+V2rev=8.02*sqrt(k/(k-1) *p1*v1*(1- (pr)^((k-1)/k)))
+v2=v1*(1/pr)^(1/k)
+m=%pi/4 *1/144 *V2rev/v2
+C=m/m1
+T2=T1*(0.887)
+t=250+460 //R
+dt=t-T2
+//results
+printf("Mass flow rate = %.3f lbm/sec",m)
+printf("\n Meta stable under cooling = %d F",dt)
diff --git a/2024/CH9/EX9.14/9_14.sce b/2024/CH9/EX9.14/9_14.sce
new file mode 100755
index 000000000..71eb7b621
--- /dev/null
+++ b/2024/CH9/EX9.14/9_14.sce
@@ -0,0 +1,12 @@
+clc
+//Initialization of variables
+C=0.98
+m=1
+v=12.55 //ft^3/lbm
+V=1372 //ft/s
+//calculations
+A=m*v/(C*V) *144
+D=sqrt(A*4/%pi)
+//results
+printf("Area = %.3f in^2",A)
+printf("\n diameter = %.2f in",D)
diff --git a/2024/CH9/EX9.15/9_15.sce b/2024/CH9/EX9.15/9_15.sce
new file mode 100755
index 000000000..0b5c59afd
--- /dev/null
+++ b/2024/CH9/EX9.15/9_15.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+nn=0.95
+p1=50 //psia
+p2=30 //psia
+v1=8.515
+m=1 //lbm
+//calculations
+cv=sqrt(nn)
+V2rev=1372
+V2act=cv*V2rev
+n=1.283
+v2=v1*(p1/p2)^(1/n)
+A=m*v2/V2act *144
+D=sqrt(A*4/%pi)
+//results
+printf("Area = %.2f in^2",A)
+printf("\n diameter = %.3f in",D)
diff --git a/2024/CH9/EX9.16/9_16.sce b/2024/CH9/EX9.16/9_16.sce
new file mode 100755
index 000000000..e2acc3321
--- /dev/null
+++ b/2024/CH9/EX9.16/9_16.sce
@@ -0,0 +1,21 @@
+clc
+//Initialization of variables
+p1=100 //psia
+p2=14.7 //psia
+k=1.4
+T1=700 //R
+R=10.73/29
+V=50
+cv=0.171
+cp=0.24
+R2=1.986/29
+//calculations
+T2=T1/ (p1/p2)^((k-1)/k)
+T2=358 //R
+m1=p1*V/(R*T1)
+m2=p2*V/(R*T2)
+Wrev= cv*(m1*T1 - m2*T2) - (m1-m2)*(T2)*cp
+Wrev2=m1*cv*(T1-T2) - m1*R2*(T2-p2/p1*T1)
+//results
+printf("Work done in case 1 = %d Btu",Wrev)
+printf("\n Work done in case 2 = %d Btu",Wrev2)
diff --git a/2024/CH9/EX9.17/9_17.sce b/2024/CH9/EX9.17/9_17.sce
new file mode 100755
index 000000000..114f85e13
--- /dev/null
+++ b/2024/CH9/EX9.17/9_17.sce
@@ -0,0 +1,13 @@
+clc
+//Initialization of variables
+hf=1187.2 //Btu/lbm
+p2=100 //psia
+//calculations
+t=328 //F
+u2=hf
+disp("from steam table,")
+t2=540 //F
+p2=100 //psia
+dt=t2-t
+//results
+printf("Rise in temperature = %d F",dt)
diff --git a/2024/CH9/EX9.2/9_2.sce b/2024/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..e728b6a32
--- /dev/null
+++ b/2024/CH9/EX9.2/9_2.sce
@@ -0,0 +1,24 @@
+clc
+//Initialization of variables
+n=1.4
+p1=50 //psia
+J=778
+cp=0.24
+T1=520 //R
+k=n
+R=1545/29
+m=1
+p2=10 //psia
+//calculations
+rpt=(2/(n+1))^(n/(n-1))
+pt=p1*rpt
+Vtrev=223.77*sqrt(cp*T1*(1- rpt^((k-1)/k)))
+v1=R*T1/p1/144
+vt=v1*(p1/pt)^(1/k)
+At=m*vt/Vtrev
+V2rev=223.77*sqrt(cp*T1*(1-(p2/p1)^((k-1)/k)))
+v2=v1*(p1/p2)^(1/k)
+A2=m*v2/V2rev
+//results
+printf("Area required = %.5f ft^2",At)
+printf("\n Area in case 2 = %.5f ft^2",A2)
diff --git a/2024/CH9/EX9.3/9_3.sce b/2024/CH9/EX9.3/9_3.sce
new file mode 100755
index 000000000..ec7847d49
--- /dev/null
+++ b/2024/CH9/EX9.3/9_3.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+rpt=0.569
+b=0.8
+p1=50 //psia
+cp=0.24
+T1=520 //R
+k=1.4
+v2=12.2
+v1=3.86
+m=1
+//calculations
+pt=p1*rpt
+Vtrev=223.77*sqrt(cp*T1*(1- rpt^((k-1)/k)) /(1- b^4 *(v1/v2)^2))
+vt=(p1/pt)^(1/k) *v1
+At=m*vt/Vtrev
+//results
+printf("Area of throat = %.5f ft^2",At)
diff --git a/2024/CH9/EX9.4/9_4.sce b/2024/CH9/EX9.4/9_4.sce
new file mode 100755
index 000000000..c1eca187c
--- /dev/null
+++ b/2024/CH9/EX9.4/9_4.sce
@@ -0,0 +1,26 @@
+clc
+//Initialization of variables
+J=778
+g=32.2
+pc=54.6 //psia
+h1=1329.1 //Btu/lbm
+h2=1265 //btu/lbm
+V2rev=1790 //ft/s
+cv=0.99
+m=1 //lbm
+cv2=0.96
+//calculations
+V2d=cv*V2rev
+hd=cv^2 *(h1-h2)
+h2d=h1-hd
+v2d=9.946
+A2d=m*v2d/V2d
+dh=-154.3
+V3=2775 //ft/s
+V3d=cv2*V3
+h3d= h1+ cv2^2 *dh
+v3d=22.05 //ft^3/lbm
+A3d=m*v3d/V3d
+//results
+printf("Throat area in case 2 = %.4f ft^2",A2d)
+printf("\n Throat area in case 3 = %.5f ft^2",A3d)
diff --git a/2024/CH9/EX9.5/9_5.sce b/2024/CH9/EX9.5/9_5.sce
new file mode 100755
index 000000000..c66cee6ed
--- /dev/null
+++ b/2024/CH9/EX9.5/9_5.sce
@@ -0,0 +1,34 @@
+clc
+//Initialization of variables
+p2=26.4 //psia
+p1=50 //psia
+p3=10 //psia
+V2rev=1017 //ft/s
+cv=0.99
+J=778
+g=32.2
+cp=0.24
+T1=460+60 //R
+k=1.4
+R=1545/29
+m=1
+cv2=0.92
+//calculations
+V2d=cv*V2rev
+dhr= (V2rev^2 - V2d^2)/(2*g*J)
+dtr=dhr/cp
+T2=T1*(p2/p1)^((k-1)/k)
+T2d=T2+dtr
+v2d=R*T2d/(p2*144)
+A2d= m*v2d/V2d
+V3=1515
+V3d=V3*cv2^0.5
+T3=T1*(p3/p1)^((k-1)/k)
+dhr2=(V3^2 - V3d^2)/(2*J*g)
+dtr2=dhr2/cp
+T3d=T3+dtr2
+v3d=R*T3d/(p3*144)
+A3d=m*v3d/V3d
+//results
+printf("Area = %.5f ft^2",A2d)
+printf("\n Area in case 2= %.5f ft^2",A3d)
diff --git a/2024/CH9/EX9.6/9_6.sce b/2024/CH9/EX9.6/9_6.sce
new file mode 100755
index 000000000..c544f9c51
--- /dev/null
+++ b/2024/CH9/EX9.6/9_6.sce
@@ -0,0 +1,10 @@
+clc
+//Initialization of variables
+mum=0.0000121 //lbm/ft sec
+D=1.820 //in
+m=1.173 //lbm/sec
+//calculations
+Re=1.27*m*12/(D*mum)
+//results
+printf("Reynolds number = %d ",Re)
+//The answer is a bit different due to rounding off error in textbook
diff --git a/2024/CH9/EX9.8/9_8.sce b/2024/CH9/EX9.8/9_8.sce
new file mode 100755
index 000000000..d4956cd84
--- /dev/null
+++ b/2024/CH9/EX9.8/9_8.sce
@@ -0,0 +1,17 @@
+clc
+//Initialization of variables
+pt=54.6 //psia
+dh=64.1 //Btu/lbm
+dh2=154.3 //Btu/lbm
+vt=9.844 //ft^3/lbm
+vt2=21.279 //ft^3/lbm
+C=0.98
+J=778
+g=32.17
+m=1 //lbm/sec
+//calculations
+At=m/C *sqrt(vt^2 /(2*g*J*dh))
+A2=m/C *sqrt(vt2^2 /(2*g*J*dh2))
+//results
+printf("Throat area = %.4f ft^2",At)
+printf("\n Area at part 2 = %.5f ft^2",A2)
diff --git a/2024/CH9/EX9.9/9_9.sce b/2024/CH9/EX9.9/9_9.sce
new file mode 100755
index 000000000..fd7d32a36
--- /dev/null
+++ b/2024/CH9/EX9.9/9_9.sce
@@ -0,0 +1,23 @@
+clc
+//Initialization of variables
+k=1.4
+D2=1.820
+D1=6.065
+p2=20.18 //psia
+p1=30 //psia
+g=32.2
+G=13.59
+zm=20
+R=1545/29
+C=0.68
+T=520 //R
+//calculations
+dp=0.03609*G*zm
+rp=p2/p1
+bet=D2/D1
+v1=R*T/(p1*144)
+A2=%pi*D2^2 /(4*144)
+m=C*A2/sqrt(1 - bet^4 *rp^(1.43)) *sqrt(2*g*k*p1*144/(k-1)/v1 *(rp^(2/k) - rp^((k+1)/k)))
+//results
+printf("Mass flow rate = %.3f lbm/sec",m)
+//The answer is a bit different due to rounding off error in textbook
generated by cgit (git 2.34.1) at 2024-12-20 05:26:10 +0000