initial commit / add all books

1 files changed, 35 insertions, 0 deletions

diff --git a/534/CH9/EX9.3/9_3_Rectangular_Duct.sce b/534/CH9/EX9.3/9_3_Rectangular_Duct.sce
new file mode 100644
index 000000000..587ded94a
--- /dev/null
+++ b/534/CH9/EX9.3/9_3_Rectangular_Duct.sce
@@ -0,0 +1,35 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.3   Page 577 \n'); //Example 9.3

+// Heat Loss from duct per meter of length

+

+//Operating Conditions

+Ts = 45+273;    //[K] Surface Temperature

+Tsurr = 15+273    ;//[K] Surrounding Temperature

+H = .3            ;//[m] Height 

+w = .75           ;//[m] Width

+

+//Table A.4 Air Properties T = 303 K

+k = 26.5*10^-3            ;//[W/m.K]

+uv = 16.2*10^-6         ;//[m^2/s] Kinematic Viscosity

+al = 22.9*10^-6        ;//[m^2/s] alpha

+be = 3.3*10^-3          ;//[K^-1]  Tf^-1

+Pr = .71               ;// Prandtl number 

+g = 9.81                ;//[m^2/s] gravitational constt

+

+Ra = g*be*(Ts-Tsurr)/al*H^3/uv;    //Length = Height

+//From equatiom 9.27

+Nu = [.68 + .67*Ra^.25/[1+(.492/Pr)^(9/16)]^(4/9)];

+//for Sides

+hs = Nu*k/H;

+

+Ra2 = g*be*(Ts-Tsurr)/al*(w/2)^3/uv;       //Length = w/2

+//For top eq 9.31

+ht = [k/(w/2)]*.15*Ra2^.3334;

+//For bottom Eq 9.32

+hb = [k/(w/2)]*.27*Ra2^.25;

+

+q = (2*hs*H+ht*w+hb*w)*(Ts-Tsurr);

+

+printf("\n Rate of heat loss per unit length of duct is %i W/m",q);

+//END
\ No newline at end of file