summaryrefslogtreecommitdiff
path: root/1709/CH12
diff options
context:
space:
mode:
Diffstat (limited to '1709/CH12')
-rwxr-xr-x1709/CH12/EX12.1/12_1.sce15
-rwxr-xr-x1709/CH12/EX12.10/12_10.sce16
-rwxr-xr-x1709/CH12/EX12.11/12_11.sce26
-rwxr-xr-x1709/CH12/EX12.12/12_12.sce20
-rwxr-xr-x1709/CH12/EX12.2/12_2.sce15
-rwxr-xr-x1709/CH12/EX12.3/12_3.sce26
-rwxr-xr-x1709/CH12/EX12.4/12_4.sce18
-rwxr-xr-x1709/CH12/EX12.5/12_5.sce22
-rwxr-xr-x1709/CH12/EX12.6/12_6.sce22
-rwxr-xr-x1709/CH12/EX12.7/12_7.sce18
-rwxr-xr-x1709/CH12/EX12.8/12_8.sce20
-rwxr-xr-x1709/CH12/EX12.9/12_9.sce21
12 files changed, 239 insertions, 0 deletions
diff --git a/1709/CH12/EX12.1/12_1.sce b/1709/CH12/EX12.1/12_1.sce
new file mode 100755
index 000000000..307bff068
--- /dev/null
+++ b/1709/CH12/EX12.1/12_1.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+disp("From Mollier diagram,")
+h1=1357 //500 psia, 700 F
+h2=935 //P2=2 psia
+h3=93.99 //sat liq at 2 psia
+vf=0.01613
+P4=500 //psia
+P3=2 //psia
+//calculations
+dh4=vf*(P4-P3)*144/778
+h4=h3+dh4
+eta= ((h1-h2)-(h4-h3))/(h1-h4)
+//results
+printf("Thermal efficiency = %.1f percent ",eta*100)
diff --git a/1709/CH12/EX12.10/12_10.sce b/1709/CH12/EX12.10/12_10.sce
new file mode 100755
index 000000000..2a2b538e2
--- /dev/null
+++ b/1709/CH12/EX12.10/12_10.sce
@@ -0,0 +1,16 @@
+clc
+//Initialization of variables
+e=0.83
+//calculations
+T1=530 //R
+T2d=838 //R
+T6d=T2d
+T3=1960 //R
+T4d=1375 //R
+T5d=T4d
+T5=e*(T5d-T2d) +T2d
+W=0.24*((T3-T4d)- (T2d-T1))
+Q=0.24*(T3-T5)
+eta=W/Q
+//results
+printf("Thermal efficiency = %d percent",eta*100+1)
diff --git a/1709/CH12/EX12.11/12_11.sce b/1709/CH12/EX12.11/12_11.sce
new file mode 100755
index 000000000..eb55fa123
--- /dev/null
+++ b/1709/CH12/EX12.11/12_11.sce
@@ -0,0 +1,26 @@
+clc
+//Initialization of variables
+T1=420 //R
+T11=530 //R
+T3=2460 //R
+V1=300 //ft/sec
+P1=5 //psia
+P5=P1
+P2=50 //psia
+P3=5 //psia
+P4=50 //psia
+g=1.4
+cp=0.24
+m=1
+//calculations
+T2=T1*(P2/P1)^((g-1)/g)
+T4=T3-T2+T11
+T5=T3*(P3/P4)^((g-1)/g)
+V5=sqrt(2*32.2*cp*(T4-T5)*778)
+T=m*(V1-V5)/32.2
+Qh=cp*(T3-T2)
+P=-T*V1
+//results
+printf("Thrust = %.1f lbf",T)
+printf("\n Heat input = %d B/lbm",Qh)
+printf("\n Power = %d ft-lbf /sec",P)
diff --git a/1709/CH12/EX12.12/12_12.sce b/1709/CH12/EX12.12/12_12.sce
new file mode 100755
index 000000000..9a9790807
--- /dev/null
+++ b/1709/CH12/EX12.12/12_12.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+h1=80.419 //B/lbm
+h3=36.013 //B/lbm
+h4=h3
+P3=172.35 //psia
+P2=P3
+m=5 //tons
+Q=12000
+//calculations
+h2=91.5 //B/lbm
+disp("From superheated steam tables,")
+COP=(h1-h4)/(h2-h1)
+W=h2-h1
+md=m*Q/(h1-h4)
+Wt=md*(h2-h1)
+Wt2=Wt/2545
+//results
+printf("Coefficient of performance = %.1f",COP)
+printf("\n Input work = %.1f hp",Wt2)
diff --git a/1709/CH12/EX12.2/12_2.sce b/1709/CH12/EX12.2/12_2.sce
new file mode 100755
index 000000000..685d6ad0a
--- /dev/null
+++ b/1709/CH12/EX12.2/12_2.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+disp("From molier diagram,")
+h1=1357 //500 psia 700F
+h2=1194 //P2=100 psia
+h3=1379 //100 psia, 700 F
+h4=1047 //p4=2 psia
+h5=93.99 //sat liq at 2 psia
+h6=95.02 //example 12.1
+//calculations
+W=h1-h2+h3-h4-(h6-h5)
+Q=(h1-h6)+(h3-h2)
+eta=W/Q
+//results
+printf("Thermal efficiency = %.2f percent",eta*100)
diff --git a/1709/CH12/EX12.3/12_3.sce b/1709/CH12/EX12.3/12_3.sce
new file mode 100755
index 000000000..c08607c7e
--- /dev/null
+++ b/1709/CH12/EX12.3/12_3.sce
@@ -0,0 +1,26 @@
+clc
+//Initialization of variables
+P=100 //psia
+//calculations
+disp("From mollier diagram,")
+h1=1357 //500 psia, 700F
+h2=1194 //100 psia
+h3=935//2 psia
+h4=93.99 //sat liq at 2 psia
+vf=0.01613
+vf2=0.01774
+P5=100 //psia
+P4=2 //psia
+dh4=vf*(P5-P4)*144/778
+h5=h4+dh4
+h6=298.4
+P7=500 //psia
+P6=100 //psia
+dh6=vf2*(P7-P6)*144/778
+h7=dh6+h6
+m=(h6-h5)/(h2-h5)
+W=h1-h2 + (1-m)*(h2-h3) - (1-m)*(h5-h4) -(h7-h6)
+Q=h1-h7
+etath=W/Q
+//results
+printf("Thermal efficiency = %.1f percent",etath*100)
diff --git a/1709/CH12/EX12.4/12_4.sce b/1709/CH12/EX12.4/12_4.sce
new file mode 100755
index 000000000..6fa8899bf
--- /dev/null
+++ b/1709/CH12/EX12.4/12_4.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+x=0.8
+//calculations
+disp("From molier diagram,")
+h1=1357 //500 psia 700F
+h2=1194 //P2=100 psia
+h3=1379 //100 psia, 700 F
+h4=1047 //p4=2 psia
+h5=93.99 //sat liq at 2 psia
+h6=95.02 //example 12.1
+h2d=h1- x*(h1-h2)
+h4d=h3- x*(h3-h4)
+W=(h1-h2d) +(h3-h4d) - (h6-h5)
+Q=(h1-h6) + (h3-h2d)
+eta=W/Q
+//results
+printf("Thermal efficiency = %d percent",eta*100+1)
diff --git a/1709/CH12/EX12.5/12_5.sce b/1709/CH12/EX12.5/12_5.sce
new file mode 100755
index 000000000..acd7db3f7
--- /dev/null
+++ b/1709/CH12/EX12.5/12_5.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+P4=50 //psia
+P1=14.7 //psia
+P3=50 //psia
+P2=14.7 //psia
+g=1.4
+//calculations
+V1r=(P4/P1)^(1/g)
+V2r=(P3/P2)^(1/g)
+//After solving,
+V4=5.38 //ft^3/min
+V1=12.9 //ft^3/min
+V2=112.9 //ft^3/min
+PD=V2-V4
+etavol=(V2-V1)/(V2-V4)
+W32=g*P2*144*V2*((P3/P2)^((g-1)/g) -1 ) /(1-g)
+W41=g*P4*144*V4*((P1/P4)^((g-1)/g) -1 ) /(1-g)
+Wt=W32+W41
+//results
+printf("Total work = %.2e ft-lbf /min",Wt)
+//The answer given in textbook is wrong . please verify it using a calculator
diff --git a/1709/CH12/EX12.6/12_6.sce b/1709/CH12/EX12.6/12_6.sce
new file mode 100755
index 000000000..e9a328de7
--- /dev/null
+++ b/1709/CH12/EX12.6/12_6.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+P1=14.7 //psia
+P4=100 //psia
+T1=530 //R
+T3=T1
+g=1.4
+m=10 //lbm
+cp=0.24
+//calculations
+P2=sqrt(P1*P4)
+T2=T1*(P2/P1)^((g-1)/g)
+T4=T2
+W=2*cp*(T2-T1)
+Wt=W*m
+hp=Wt*60/2545
+Q=m*cp*(T2-T3)
+T4=T1*(P4/P1)^((g-1)/g)
+W2=m*cp*(T4-T1)
+//results
+printf("Work required in case 1 = %d Btu/min",Wt+1)
+printf("\n Work required in case 2 = %d Btu/min",W2+1)
diff --git a/1709/CH12/EX12.7/12_7.sce b/1709/CH12/EX12.7/12_7.sce
new file mode 100755
index 000000000..f4c8e7546
--- /dev/null
+++ b/1709/CH12/EX12.7/12_7.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+g=1.4
+r1=10
+r2=12
+r3=15
+Tl=530 //R
+Th=1960 //R
+//calculations
+eta1=1- (r1)^(1-g)
+eta2=1- (r2)^(1-g)
+eta3=1- (r3)^(1-g)
+etac=1-Tl/Th
+//results
+printf("Efficiency in case 1 = %.1f percent",eta1*100)
+printf("\n Efficiency in case 2 = %.1f percent",eta2*100)
+printf("\n Efficiency in case 3 = %.1f percent",eta3*100)
+printf("\n Carnot efficiency = %.2f percent",etac*100)
diff --git a/1709/CH12/EX12.8/12_8.sce b/1709/CH12/EX12.8/12_8.sce
new file mode 100755
index 000000000..deda55bbd
--- /dev/null
+++ b/1709/CH12/EX12.8/12_8.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+T1=70+460 //R
+P1=14.7 //psia
+g=1.4
+r=15
+rc=2
+cp=0.24
+cp2=0.1715
+//calculations
+T2=T1*(r)^(g-1)
+T3=rc*T2
+T4=T3*(rc/r)^(g-1)
+Qh=cp*(T3-T2)
+Ql=cp2*(T4-T1)
+W=Qh-Ql
+eta=W/Qh
+//results
+printf("Work output = %d B/lbm",W)
+printf("\n Efficiency = %.1f percent",eta*100)
diff --git a/1709/CH12/EX12.9/12_9.sce b/1709/CH12/EX12.9/12_9.sce
new file mode 100755
index 000000000..341624e9c
--- /dev/null
+++ b/1709/CH12/EX12.9/12_9.sce
@@ -0,0 +1,21 @@
+clc
+//Initialization of variables
+P1=14.7 //psia
+P4=14.7 //psia
+T1=530 //R
+T3=1960 //R
+P2=60 //psia
+P3=P2
+g=1.4
+eta1=0.85
+eta2=0.9
+//calculations
+T2=T1*(P2/P1)^((g-1)/g)
+T4=T3*(P4/P3)^((g-1)/g)
+T2d=(T2-T1)/eta1 + T1
+T4d=-eta2*(T3-T4) +T3
+Wact=0.24*(T3-T4d - (T2d-T1))
+Qh=0.24*(T3-T2d)
+etath=Wact/Qh
+//results
+printf("Thermal efficiency = %.1f percent",etath*100)
generated by cgit (git 2.34.1) at 2024-12-20 04:05:30 +0000