initial commit / add all books

1 files changed, 26 insertions, 0 deletions

diff --git a/1328/CH16/EX16.1/16_1.sce b/1328/CH16/EX16.1/16_1.sce
new file mode 100644
index 000000000..94d8f7383
--- /dev/null
+++ b/1328/CH16/EX16.1/16_1.sce
@@ -0,0 +1,26 @@
+printf("\t example 16.1 \n");

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

+Af=(20*0.75*12*2)/(144);

+Ao=((3.14*1.25)-(20*0.035))*(12/144);

+printf("\t fin surface is : %.1f ft^2/lin ft \n",Af);

+printf("\t bare tube surface is : %.3f ft^2/lin ft \n",Ao);

+A=(Af+Ao);

+printf("\t total outside surface : %.2f ft^2/lin ft \n",A);

+Ai=(3.14*1.06*12)/(144);

+printf("\t total inside surface : %.3f ft^2/lin ft \n",Ai);

+printf("\t fin efficiencies \n");

+b=0.0625; // ft

+hf=4; // from table in solution

+m=(5.24*(hf^(1/2))); // m=((hf*P)/(Kax))^(1/2), eq 16.8

+n=(tanh(m*b))/(m*b); // efficiency , eq 16.26

+printf("\n hf      m      n \n "+string(hf)+"      "+string(m)+"      "+string(n)+" \n");

+// similarly efficiencies values are calculated at different hf values

+printf("\t weighted efficiency curve \n");

+hfi=((n*Af)+(Ao))*(hf/Ai); // eq 16.34

+printf("\n hf      hfi \n "+string(hf)+"      "+string(hfi)+" \n");

+// similarly efficiencies values are calculated at different hf values

+hf=[4 16 36 100 400 625 900]; // from 2nd table in the solution

+hfi=[35.4 110.8 193.5 370 935 1295 1700]; // from 2nd table in the solution

+plot2d("oll",hf,hfi);

+xtitle("weighted fin efficiency curve","heat transfer coefficient to fin,Btu/(ft^2)*(hr)","coefficient hf referred to the tube ID");

+//end