initial commit / add all books

7 files changed, 344 insertions, 0 deletions

diff --git a/3516/CH17/EX17.1/Ex17_1.sce b/3516/CH17/EX17.1/Ex17_1.sce
new file mode 100644
index 000000000..a8e5bc6be
--- /dev/null
+++ b/3516/CH17/EX17.1/Ex17_1.sce
@@ -0,0 +1,11 @@
+printf("\t example 17.1 \n");

+pw=0.4298; // psia, at 75F, table 7

+pt=14.696; // psia

+t=75;

+Mw=18;

+Ma=29;

+X=(pw/(pt-pw))*(Mw/Ma);

+printf("\t humidity is : %.4f lb water/lb air \n",X);

+H=(X*t)+(1051.5*X)+(0.24*t); // eq 17.54

+printf("\t enthalpy at 75F is : %.1f Btu/lb dry air \n",H);

+// end

diff --git a/3516/CH17/EX17.2/Ex17_2.sce b/3516/CH17/EX17.2/Ex17_2.sce
new file mode 100644
index 000000000..e4283e5ac
--- /dev/null
+++ b/3516/CH17/EX17.2/Ex17_2.sce
@@ -0,0 +1,40 @@
+printf("\t example 17.2 \n");

+printf("\t approximate values are mentioned in the book \n");

+printf("\t by numerical integration \n");

+T1=85;

+T2=120;

+A=576; // ground area, from fig 17.12

+L=1500*(500/576);

+G=1400;

+R=(L/G);

+printf("\t R is : %.2f \n",R);

+H1=39.1; // fig 17.12

+H2=H1+(R*(T2-T1));

+printf("\t H2 is : %.1f Btu \n",H2);

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=85F

+Hs=50; // fig 17.12

+d1=(Hs-H1);

+printf("\t difference is : %.1f \n",d1);

+//at t=90

+Hs=56.7; // fig 17.12

+H=43.7; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.1f \n",d);

+dT=5; // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.70;

+printf("\t number of diffusing units : %.2f \n",nd);

+printf("\t log mean enthalpy difference \n");

+dt=49.9; // diff. of enthalpies at top of the tower, from table in solution

+db=10.9; // diff of enthalpies at bottom of the tower,from table in solution

+LME=(dt-db)/(2.3*log10(dt/db));

+printf("\t log mean of enthalpy : %.1f Btu/lb \n",LME);

+nd=(T2-T1)/(LME);

+printf("\t number of diffusing units are : %.2f \n",nd);

+// The error is naturally larger the greater the range

+//end

diff --git a/3516/CH17/EX17.3/Ex17_3.sce b/3516/CH17/EX17.3/Ex17_3.sce
new file mode 100644
index 000000000..6b8d77ef0
--- /dev/null
+++ b/3516/CH17/EX17.3/Ex17_3.sce
@@ -0,0 +1,11 @@
+printf("\t example 17.3 \n");

+printf("\t approximate values are mentioned in the book \n");

+// Since the loading is based on 1 ft2 of ground area

+nd=1.7;

+L=1302;

+Kxa=115;

+Z=(nd*L)/(Kxa);

+printf("\t Z is : %.1f ft \n",Z);

+HDU=(Z/nd);

+printf("\t height of diffusion unit : %.1ff ft \n",HDU);

+// end

diff --git a/3516/CH17/EX17.4/Ex17_4.sce b/3516/CH17/EX17.4/Ex17_4.sce
new file mode 100644
index 000000000..e39675bc4
--- /dev/null
+++ b/3516/CH17/EX17.4/Ex17_4.sce
@@ -0,0 +1,22 @@
+printf("\t example 17.4 \n");

+printf("\t approximate values are mentioned in the book \n");

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=79.3F

+Hs=43.4; // fig 17.12

+H=30.4; // fig 17.12

+d1=(Hs-H);

+printf("\t difference is : %.1f \n",d1);

+//at t=85

+Hs=50; // fig 17.12

+H=35.7; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.2f \n",d);

+dT=(85-79.3); // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.72;

+printf("\t number of diffusing units : %.2f \n",nd);

+// end

diff --git a/3516/CH17/EX17.5/Ex17_5.sce b/3516/CH17/EX17.5/Ex17_5.sce
new file mode 100644
index 000000000..5b5774818
--- /dev/null
+++ b/3516/CH17/EX17.5/Ex17_5.sce
@@ -0,0 +1,110 @@
+printf("\t example 17.5 \n");

+printf("\t approximate values are mentioned in the book \n");

+T1=85;

+T2=120;

+R=0.93; // R=(L/G), for 1500 gpm

+printf("\t for 120percent of design \n");

+R1=1.2*R;

+printf("\t R is : %.3f \n",R1);

+H1=39.1; // at 87.2F

+H2=H1+(R1*(T2-T1)); 

+printf("\t H2 is : %.1f Btu \n",H2);

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=87.2F

+Hs=53.1; // from table in the solution

+d1=(Hs-H1);

+printf("\t difference is : %.1f \n",d1);

+//at t=90

+Hs=56.7; // fig 17.12

+H=42; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.1f \n",d);

+dT=(90-87.2); // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.53;

+printf("\t number of diffusing units : %.2f \n",nd);

+printf("\t for 80 percent of design \n");

+R2=0.8*R;

+printf("\t R is : %.3f \n",R2);

+H1=39.1; // at 87.2F

+H2=H1+(R2*(T2-T1)); 

+printf("\t H2 is : %.0f Btu \n",H2);

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=82.5F

+Hs=47.2; // from table in the solution

+d1=(Hs-H1);

+printf("\t difference is : %.1f \n",d1);

+//at t=85

+Hs=50; // fig 17.12

+H=40.8; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.1f \n",d);

+dT=(85-82.5); // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.92;

+printf("\t number of diffusing units : %.2f \n",nd);

+X=[1.115 0.93 0.74];

+Y=[1.53 1.70 1.92];

+plot2d(X,Y,style=3,rect=[0.7,1.4,1.3,2]);

+xtitle("KxaV/L vs L/G","L/G","nd");

+printf("\t trial 1 \n");

+R3=1.1;

+printf("\t R is : %.3f \n",R3);

+H1=34.5; // at 87.2F

+H2=H1+(R3*(T2-T1)); 

+printf("\t H2 is : %.0f Btu \n",H2);

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=85F

+Hs=50; // from table in the solution

+d1=(Hs-H1);

+printf("\t difference is : %.1f \n",d1);

+//at t=90

+Hs=56.7; // fig 17.12

+H=40; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.1f \n",d);

+dT=(90-85); // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.48;

+printf("\t number of diffusing units : %.2f \n",nd);

+R3=1.19; // from fig 17.14

+printf("\t L/G is : %.2f \n",R3);

+printf("\t trial 2 \n");

+R4=1.2;

+printf("\t R4 is : %.3f \n",R4);

+H1=34.5; // at 87.2F

+H2=H1+(R4*(T2-T1)); 

+printf("\t H2 is : %.1f Btu \n",H2);

+// The area between the saturation line and the operating line represents the potential for heat transfer

+// at T=85F

+Hs=50; // from table in the solution

+d1=(Hs-H1);

+printf("\t difference is : %.1f \n",d1);

+//at t=90

+Hs=56.7; // fig 17.12

+H=40.5; // fig 17.12

+d2=Hs-H;

+printf("\t difference is : %.1f \n",d2);

+d=(d1+d2)/(2);

+printf("\t average of difference is : %.1f \n",d);

+dT=(90-85); // F

+nd1=(dT/d);

+printf("\t nd1 is : %.3f \n",nd1);

+// similarly calculating nd at each temperature and adding them will give you total nd value

+nd=1.56;

+printf("\t number of diffusing units : %.2f \n",nd);

+R3=1.08; // from fig 17.14

+printf("\t L/G is : %.2f \n",R3);

+// end

diff --git a/3516/CH17/EX17.6/Ex17_6.sce b/3516/CH17/EX17.6/Ex17_6.sce
new file mode 100644
index 000000000..4f8e7d2ee
--- /dev/null
+++ b/3516/CH17/EX17.6/Ex17_6.sce
@@ -0,0 +1,113 @@
+printf("\t example 17.6 \n");

+printf("\t approximate values are mentioned in the book \n");

+// basis 1ft^2 ground area

+//Assumption: 20 per cent of the initial vapor content of the gas enters the water body

+X1=(1.69/(14.7-1.69))*(18/29);

+printf("\t X1 : %.4f lb/lb \n",X1);

+G=1500;

+w1=G*X1;

+printf("\t total water in inlet gas : %.2f lb/hr \n",w1);

+// The inlet gas is at 300F and a 120F dew point. Use 0.25 Btu/(lb)(°F) for the specific heat of nitrogen

+H1=(0.0807*120)+(0.0807*1025.8)+(0.45*0.0807*(300-120))+(0.25*300); // eq 17.55

+printf("\t H1 : %.0f Btu/lb dry air \n",H1);

+X2=(w1*(1-.2)/G);

+printf("\t outlet gas humidity : %.5f lb/lb \n",X2);

+pw=(X2*29*14.7/18)/(1+(X2*29/18));

+printf("\t pw : %.3f psia \n",pw);

+Tw=112.9; // F, from table 7 for above pw

+// The outlet gas has a temperature of 200°F and a 112.9°F dew point

+H2=(X2*Tw)+(X2*1029.8)+(X2*0.45*(200-Tw))+(0.25*200); // eq 17.55

+printf("\t H2 : %.1f Btu/lb dry air \n",H2);

+q=G*(H1-H2);

+printf("\t total heat load : %.2e Btu/hr \n",q);

+w2=q/(120-85);

+printf("\t water loading : %.2e lb/hr \n",w2);

+printf("\t interval 1 \n");

+// (Kxa*delV/L)= 0 t0 0.05

+nd=0.05; // nd=Kxa*V/L

+Le=0.93; // fig 17.4 at 300F

+C=(0.25)+(0.45*X1);

+printf("\t C : %.3f Btu/(lb)*(F) \n",C);

+haV=(nd*w2*Le*C);

+printf("\t haV : %.1f Btu/(hr)*(F) \n",haV);

+qc=(haV*(300-120));

+printf("\t qc : %.2e Btu/hr \n",qc);

+delT=(qc/(C*G));

+printf("\t delT : %.1f F \n",delT);

+T1=(300-delT);

+printf("\t T(0.05) : %.1f F \n",T1);

+delt=(qc/w2);

+printf("\t delt : %.2f F \n",delt);

+t1=(120-delt);

+printf("\t t(0.05) : %.1f F \n",t1);

+printf("\t interval 2 \n");

+// (Kxa*delV/L)= 0.05 to 0.15

+nd1=0.1;

+haV1=(nd1*w2*Le*C);

+printf("\t haV1 : %.1f Btu/(hr)*(F) \n",haV1);

+qc1=(haV1*(T1-t1));

+printf("\t qc1 : %.1e Btu/hr \n",qc1);

+delT1=(qc1/(C*G));

+printf("\t delT1 : %.1f F \n",delT1);

+T2=(T1-delT1);

+printf("\t T(0.15) : %.2f F \n",T2);

+X3=0.0748; // at 117.6F

+w3=(nd1*w2*(0.0807-X3));

+printf("\t water diffused during interval : %.3f lb/hr \n",w3);

+w4=(w1-w3);

+printf("\t water remaining : %.2f lb/hr \n",w4);

+l1=1027; // Btu/lb, l1= lamda at 117.6F

+qd=(w3*l1);

+printf("\t qd : %.0f Btu/hr \n",qd);

+q1=(qd+qc1);

+printf("\t q1 : %.0f Btu/hr \n",q1);

+delt1=(q1/w2);

+printf("\t delt1 : %.2f F \n",delt1);

+t2=(t1-delt1);

+printf("\t t(0.15) : %.1f F \n",t2);

+X4=0.0640; // at 112.5

+X5=(w4/G);

+printf("\t X(112.5F) : %.4f lb/lb \n",X5);

+printf("\t interval 3 \n");

+// (Kxa*delV/L)= 0.15 to 0.25

+nd1=0.1;

+haV1=(nd1*w2*Le*C);

+printf("\t haV1 : %.1f Btu/(hr)*(F) \n",haV1);

+qc2=(haV1*(T2-t2));

+printf("\t qc2 : %.2e Btu/hr \n",qc2);

+delT2=(qc2/(C*G));

+printf("\t delT2 : %.1f F \n",delT2);

+T3=(T2-delT2);

+printf("\t T(0.25) : %.1f F \n",T3);

+w5=(nd1*w2*(X5-X4));

+printf("\t water diffused during interval : %.3f lb/hr \n",w5);

+w6=(w4-w5);

+printf("\t water remaining : %.2f lb/hr \n",w6);

+l2=1030; // Btu/lb, l1= lamda at 112.5F

+qd1=(w5*l2);

+printf("\t qd1 : %.2e Btu/hr \n",qd1);

+q2=(qd1+qc2);

+printf("\t q2 : %.3e Btu/hr \n",q2);

+delt2=(q2/w2);

+printf("\t delt2 : %.2f F \n",delt2);

+t3=(t2-delt2);

+printf("\t t(0.25) : %.1f F \n",t3);

+X6=0.0533; // at 106.5

+X7=(w6/G);

+printf("\t X(106.5F) : %.4f lb/lb \n",X7);

+// The calculations of the remaining intervals until a. gas temperature of 200°F is reached are shown in Fig. 17.17

+w7=21.92; // total water diffused from table in solution

+d=(w7/w1)*100;

+printf("\t calculated diffusion : %.0f \n",d);

+printf("\t Using some standard low-pressure-drop data \n");

+// For G = 1500, extrapolate to L = 2040 on logarithmic coordinates. Kxa = 510.

+ndt=.54; // from 1st table in solution

+Kxa=510; // from 2nd table in solution

+Z=(ndt*w2/Kxa);

+printf("\t tower height : %.2f ft \n",Z);

+A=(50000/G);

+printf("\t cross section : %.1f ft^2 \n",A);

+// end

+

+

+

diff --git a/3516/CH17/EX17.7/Ex17_7.sce b/3516/CH17/EX17.7/Ex17_7.sce
new file mode 100644
index 000000000..a75eaabbb
--- /dev/null
+++ b/3516/CH17/EX17.7/Ex17_7.sce
@@ -0,0 +1,37 @@
+printf("\t example 17.7 \n");

+printf("\t approximate values are mentioned in the book \n");

+C=0.28; // assumption

+w=50000; // lb/hr

+G=1500;

+Qs=(w*C*(500-200));

+Qd=(w/G)*(22685); // qd=22685, from previous prblm

+printf("\t sensible heat : %.1e Btu/hr \n",Qs);

+printf("\t approximate diffusion : %.2e Btu/hr \n",Qd);

+Q=(Qs+Qd);

+printf("\t total heat : %.3e Btu/hr \n",Q);

+// an allowance as high as 30 per cent of the sensible load can be made and the excess water compensated for by throttling when the tower is in operation

+w1=(Q/(120-85));

+printf("\t total water quantity : %.2e lb/hr \n",w1);

+// If the maximum liquid loading is taken as 2040 lb/(hr)(ft'!), the required tower cross section

+A=(w1/2040);

+printf("\t tower cross section : %.1f ft^2 \n",A);

+w3=(w/A);

+printf("\t new gas rate : %.0f lb/(hr)(ft^2) \n",w3);

+// The two terminal temperature differences are (200 - 85) and (500 - 120).

+LMTD=((500-120)-(200-85))/(log((500-120)/(200-85)));

+printf("\t LMTD : %.0f \n",LMTD);

+dt=35;

+N=(dt/LMTD); // eq 17.88

+printf("\t haV/L : %.2f \n",N);

+Le=0.93;

+nd=(N/(C*Le));

+printf("\t number diffusion units : %.2f \n",nd);

+// By extrapolation for G = 718 and L = 2040,Kxa=215

+L=2040;

+Kxa=215;

+Z=(nd*L/Kxa); // calculation mistake

+printf("\t height of tower : %.1f ft \n",Z);

+di=(A)^(1/2);

+printf(" ground dimensions : %.1f ft \n",di);

+// ground dimensions are 5.8*8.3*8.3 ft

+// end