initial commit / add all books

15 files changed, 250 insertions, 0 deletions

diff --git a/2024/CH8/EX8.1/8_1.sce b/2024/CH8/EX8.1/8_1.sce
new file mode 100755
index 000000000..8862c59ad
--- /dev/null
+++ b/2024/CH8/EX8.1/8_1.sce
@@ -0,0 +1,9 @@
+clc

+//Initialization of variables

+R=1.986 //B/mol R

+t2=1100 //R

+t1=1000 //R

+//calculations

+wrev=R*(t2-t1)

+//results

+printf("work done = %.1f Btu/mol",wrev)

diff --git a/2024/CH8/EX8.10/8_10.sce b/2024/CH8/EX8.10/8_10.sce
new file mode 100755
index 000000000..cca0ee27c
--- /dev/null
+++ b/2024/CH8/EX8.10/8_10.sce
@@ -0,0 +1,13 @@
+clc

+//Initialization of variablesk=1.38

+R=1.986

+T1=900 //R

+M=29

+pr=0.1

+Wrev=50 //Btu/lbm

+cp=0.245

+k=1.3

+//calculations

+KE = -cp*T1*(pr^((k-1)/k) -1) -Wrev

+//results

+printf("Change in kinetic energy = %d Btu/lbm",KE)

diff --git a/2024/CH8/EX8.11/8_11.sce b/2024/CH8/EX8.11/8_11.sce
new file mode 100755
index 000000000..aa11aa69c
--- /dev/null
+++ b/2024/CH8/EX8.11/8_11.sce
@@ -0,0 +1,21 @@
+clc

+//Initialization of variables

+n=1.3

+p2=125 //psia

+m=1 //lbm

+c=0.04

+cv=0.171

+k=1.4

+p1=14.7 //psia

+T2=852 //R

+T1=520 //R

+//calculations

+eta=1+c-c*((p2/p1)^(1/n)) 

+md=m/eta

+m12=md*(c+1)

+m34=m12-m

+Q12=m12*cv*((k-n)/(1-n))*(T2-T1)

+Q34=m34*cv*((k-n)/(1-n))*(T1-T2)

+Q=Q12+Q34

+//results

+printf("Net heat transfer from air = %.1f Btu/lbm ",Q)

diff --git a/2024/CH8/EX8.12/8_12.sce b/2024/CH8/EX8.12/8_12.sce
new file mode 100755
index 000000000..948c141ee
--- /dev/null
+++ b/2024/CH8/EX8.12/8_12.sce
@@ -0,0 +1,31 @@
+clc

+//Initialization of variables

+stroke=14 //in

+n=1.3

+rpm=130

+pa=14.7 //psia

+ta=80+460 //R

+c=0.03

+p1=pa

+p4=200 //psia

+R=53.3

+cap=400 //cfm

+m=29.4

+cp=0.24

+//calculations

+pi=sqrt(pa*p4)

+nv=1+c-c*(pi/p1)^(1/n)

+Dl=cap/nv

+vdl=Dl/(2*rpm)

+dl=sqrt(vdl*1728*4/(%pi*stroke))

+vh=cap*p1/pi

+Dh=vh/nv

+vdh=Dh/(2*rpm)

+dh=sqrt(vdh*1728*4/(%pi*stroke))

+m=p1*144*cap/R/ta

+T2=ta*(pi/p1)^((n-1)/n)

+Q=m*cp*(T2-ta)

+//results

+printf("diameter of cylinder 1 = %d in",dl)

+printf("\n diameter of cylinder 2 = %.2f in",dh)

+printf("\n Heat transferred = %d Btu/min",Q)

diff --git a/2024/CH8/EX8.13/8_13.sce b/2024/CH8/EX8.13/8_13.sce
new file mode 100755
index 000000000..8fe9ed49f
--- /dev/null
+++ b/2024/CH8/EX8.13/8_13.sce
@@ -0,0 +1,14 @@
+clc

+//Initialization of variables

+Pr=10

+n=1.3

+T1=900 //R

+W=50 //Btu/lbm

+//calculations

+T2=T1/Pr^((n-1)/n)

+h1=120.86

+h2=30.69

+dh=h2-h1

+ke=-dh-W

+//results

+printf("Change in kinetic energy = %.2f Btu/lbm",ke)

diff --git a/2024/CH8/EX8.14/8_14.sce b/2024/CH8/EX8.14/8_14.sce
new file mode 100755
index 000000000..ef14a6f43
--- /dev/null
+++ b/2024/CH8/EX8.14/8_14.sce
@@ -0,0 +1,16 @@
+clc

+//Initialization of variables

+T1=900 //R

+p1=100 //psia

+p2=10 //psia

+w=50 //Btu/lbm

+//calculations

+h1=120.86 //Btu/lbm

+pr1=17.374

+pr2=pr1*p2/p1

+disp("From equilibrium charts,")

+T2=468 //R

+h2=16.3 //Btu/lbm

+ke=h1-h2-w

+//results

+printf("Change in kinetic energy = %.2f Btu/lbm",ke)

diff --git a/2024/CH8/EX8.15/8_15.sce b/2024/CH8/EX8.15/8_15.sce
new file mode 100755
index 000000000..699c63a1e
--- /dev/null
+++ b/2024/CH8/EX8.15/8_15.sce
@@ -0,0 +1,21 @@
+clc

+//Initialization of variables

+p1=100 //psi

+p2=10 //psia

+pa=14.7 //psi

+T2=468 //R

+T2r=528 //R

+R=1.986

+M=29

+//calculations

+disp("From air tables,")

+phi1=0.06657

+phi2=0.03762

+ds=phi1-phi2

+phi11=0.19569

+phi12=0.06657

+ds2= phi12-phi11 - R/M *log(p2/p1)

+//results

+printf("In case 1, change in entropy = %.5f Btu/lbm R",ds)

+printf("\n In case 2, change in entropy = %.5f btu/lbm R",ds2)

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

diff --git a/2024/CH8/EX8.16/8_16.sce b/2024/CH8/EX8.16/8_16.sce
new file mode 100755
index 000000000..efccee032
--- /dev/null
+++ b/2024/CH8/EX8.16/8_16.sce
@@ -0,0 +1,30 @@
+clc

+//Initialization of variables

+T1=520 //R

+disp("From air table,")

+vr1=5192

+u1=-6.87 //Btu/lbm

+pr1=2.504

+vrat=6

+p1=14.7

+R=1.986

+M=29

+//calculations

+vr2=vr1/vrat

+T2=1050 //R

+u2=86.1 //Btu/lbm

+pr2=30.35 

+p2=p1*pr2/pr1

+W=u1-u2

+k=1.39

+p22=p1*vrat^(k)

+T22=T1*(vrat)^(k-1)

+W2=R*(T22-T1)/(1-k)/M

+//results

+printf("in case 1, Final pressure = %d psia",p2)

+printf("\n in case 1, final temperature = %d R",T2)

+printf("\n in case 1, work done = %.2f Btu/lbm",W)

+printf("\n in case 2, Final pressure = %d psia",p22)

+printf("\n in case 2, final temperature = %d R",T22)

+printf("\n in case 2, work done = %.2f Btu/lbm",W2)

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

diff --git a/2024/CH8/EX8.2/8_2.sce b/2024/CH8/EX8.2/8_2.sce
new file mode 100755
index 000000000..c7c6de662
--- /dev/null
+++ b/2024/CH8/EX8.2/8_2.sce
@@ -0,0 +1,16 @@
+clc

+//Initialization of variables

+p1=20 //psia

+p2=40 //psia

+t1= 460+40 //R

+m=28

+R=1.986

+cp=0.246 //B/lbm R

+//calculations

+t2= t1*p2/p1

+cv=cp- (R/m)

+Qrev=cv*(t2-t1)

+ds= cv*log(t2/t1)

+//results

+printf("heat transferred = %.1f Btu/lbm",Qrev)

+printf("\n change in entropy = %.3f Btu/lbm R",ds)

diff --git a/2024/CH8/EX8.3/8_3.sce b/2024/CH8/EX8.3/8_3.sce
new file mode 100755
index 000000000..a2ee79fa7
--- /dev/null
+++ b/2024/CH8/EX8.3/8_3.sce
@@ -0,0 +1,16 @@
+clc

+//Initialization of variables

+t1=500 //R

+t2=1000 //R

+//calculations

+function y = cp1(t)

+    y= 7.484 - 3.47*10^3 /t + 1.16*10^6 /t^2

+endfunction

+function y = cp2(t)

+    y = 7.484/t - 3.47*10^3 /t^2 + 1.16*10^6 /t^3

+endfunction

+Q=intg(t1,t2,cp1)

+ds=intg(t1,t2,cp2)

+//results

+printf("heat transferred = %d Btu/mole",Q)

+printf("\n change in entropy = %.3f Btu/mole R",ds)

diff --git a/2024/CH8/EX8.4/8_4.sce b/2024/CH8/EX8.4/8_4.sce
new file mode 100755
index 000000000..6f465d288
--- /dev/null
+++ b/2024/CH8/EX8.4/8_4.sce
@@ -0,0 +1,12 @@
+clc

+//Initialization of variables

+v1=20.9 //ft^3/mol

+v2=23.2 //ft^3/mol

+p=500 //psia

+w1=198.6 //Btu/mol

+//calculations

+w=p*(v2-v1)*144/100*0.1285

+err = (w-w1)/w

+//results

+printf("Work done in this case = %d Btu/mol",w)

+printf("\n error = %.2f percent",err*100)

diff --git a/2024/CH8/EX8.5/8_5.sce b/2024/CH8/EX8.5/8_5.sce
new file mode 100755
index 000000000..b354c71c7
--- /dev/null
+++ b/2024/CH8/EX8.5/8_5.sce
@@ -0,0 +1,9 @@
+clc

+//Initialization of variables

+R=1.986

+T=1000 //R

+vr=2

+//calculations

+Q= R*T*log(vr)

+//results

+printf("heat transferred = work = %d Btu/mol",Q)

diff --git a/2024/CH8/EX8.7/8_7.sce b/2024/CH8/EX8.7/8_7.sce
new file mode 100755
index 000000000..db8a6b5ce
--- /dev/null
+++ b/2024/CH8/EX8.7/8_7.sce
@@ -0,0 +1,13 @@
+clc

+//Initialization of variables

+v2=41.8 //ft^3/mol

+v1=20.9 //ft^3/mol

+b=0.685 //ft^3/mol

+R=0.73 //atm ft^3 / R mol

+a=924.2

+T=1000 //R

+//calculations

+vr= log((v2-b)/(v1-b))

+W= R*T*vr + a*(1/v2 - 1/v1)

+//results

+printf("Work done = %.1f atm ft^3/mol",W)

diff --git a/2024/CH8/EX8.8/8_8.sce b/2024/CH8/EX8.8/8_8.sce
new file mode 100755
index 000000000..1fbbb22ee
--- /dev/null
+++ b/2024/CH8/EX8.8/8_8.sce
@@ -0,0 +1,17 @@
+clc

+//Initialization of variables

+R=1545 

+n=1.3

+T1=520 //R

+p2=125 //psia

+p1=14.7 //psia

+M=29

+cv=0.171

+k=1.4

+//calculations

+Wrev= R*T1/M/(1-n) *((p2/p1)^((n-1)/n) -1)

+T2= T1*(p2/p1)^((n-1)/n)

+Qrev= cv*((k-n)/(1-n))*(T2-T1)

+//results

+printf("Work done = %d ft lbf/lbm",Wrev)

+printf("\n Heat transferred = %.1f Btu/lbm",Qrev)

diff --git a/2024/CH8/EX8.9/8_9.sce b/2024/CH8/EX8.9/8_9.sce
new file mode 100755
index 000000000..f9a34fb67
--- /dev/null
+++ b/2024/CH8/EX8.9/8_9.sce
@@ -0,0 +1,12 @@
+clc

+//Initialization of variables

+k=1.38

+R=1.986

+T1=900 //R

+M=29

+pr=0.1

+Wrev=50 //Btu/lbm

+//calculations

+KE = k*R*T1/M/(1-k) *(pr^((k-1)/k) -1) -Wrev

+//results

+printf("Change in kinetic energy = %d Btu/lbm",KE)