summaryrefslogtreecommitdiff
path: root/2510/CH15/EX15.11
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /2510/CH15/EX15.11
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '2510/CH15/EX15.11')
-rwxr-xr-x2510/CH15/EX15.11/Ex15_11.sce39
1 files changed, 39 insertions, 0 deletions
diff --git a/2510/CH15/EX15.11/Ex15_11.sce b/2510/CH15/EX15.11/Ex15_11.sce
new file mode 100755
index 000000000..2b49cce78
--- /dev/null
+++ b/2510/CH15/EX15.11/Ex15_11.sce
@@ -0,0 +1,39 @@
+
+//Variable declaration:
+h1 = 1200.0 //Hot film coefficient (Btu/h.ft^2..)
+h2 = 1175.0 //Cold film coefficient (Btu/h.ft^2..)
+L = 200.0 //Length of pipe (ft)
+MC = 30000.0
+mc = 22300.0
+T1 = 300.0 //Inlet temperature of hot fluid in pipe (.)
+t1 = 60.0 //Inlet temperature of cold fluid in pipe (.)
+syms T2 //Outlet temperature of hot fluid .
+syms t2 //Outlet temperature of cold fluid .
+//From table 6.2:
+ID = 2.067 //Inside diameter of pipe (in)
+OD = 2.375 //Outside diameter of pipe (in)
+Dx = 0.154 //Thickness of pipe (in)
+Ai = 0.541 //Inside sectional area of pipe (ft^2/ft)
+k = 25.0 //Thermal conductivity of pipe (Btu/h)
+
+//Calculation:
+Ui = 1.0/((1.0/h1) +(Dx/(k*12.0))+(1.0/(h2*(OD/ID)))) //Overall heat transfer coefficient (Btu/h.ft^2..)
+Ai1 = Ai*L //Inside area of pipe (ft^3/ft)
+QH = MC*(T1-T2) //Heat transfer rate of hot fluid (Btu/h)
+QC = mc*(t2-t1) //Heat transfer rate of cold fluid (Btu/h)
+t2ht = 195 //t2 by hit and trial
+[x] = fsolve(T2,QC-QH)
+T2 = x(1)
+DTlm = (T1-t1-T2+t2)/log((T1-t1)/(T2-t2)) //Log mean temperature difference (.)
+Q = Ui*Ai1*subst(t2ht,t2,DTlm) //Total heat transfer rate (Btu/h)
+
+//Result:
+disp("T2 :")
+disp(subst(t2ht,t2,T2))
+
+disp("t2 :")
+disp(subst(t2ht,t2,t2))
+
+disp("Qdot :")
+disp(Q/10**6)
+disp("x 10**6 Btu/h")