From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 965/CH9/EX9.18/18.sci | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)
 create mode 100644 965/CH9/EX9.18/18.sci

(limited to '965/CH9/EX9.18/18.sci')

diff --git a/965/CH9/EX9.18/18.sci b/965/CH9/EX9.18/18.sci
new file mode 100644
index 000000000..2a9409f45
--- /dev/null
+++ b/965/CH9/EX9.18/18.sci
@@ -0,0 +1,22 @@
+clc;
+clear all;
+disp("heat transfer rate")
+L=3.2;//m
+d=0.006;//m diameter
+ts=54;// degree C
+tsat=100;//degree C
+rhol=973.7;//kg/m^3
+mu=365*10^(-6);// Ns/m^2
+k=0.668;// W/m.C
+rhov=0.596;// kg/m^3
+hfg=2257*10^3;// J/kg
+g=9.81;//m/s
+disp("h=0.0077*(rhol*(rhol-rhov)*k^2*g/(mu^2))^0.333*Re^0.4")
+disp("Eliminating h from euqation we get the condition that the flow will be turbulent if ")
+disp("0.00296*((rhol*(rhol-rhov)*k^3*g*(tsat-ts)^3*L^3/(mu^5*hfg^3))^(5/9)>1800")
+x=0.00296*(rhol*(rhol-rhov)*k^3*g*(tsat-ts)^3*L^3/(mu^5*hfg^3))^(5/9);
+if(x>1800)
+Re=x
+h=0.0077*(rhol*(rhol-rhov)*k^2*g/(mu^2))^0.333*Re^0.4
+Q=h*L*1*(tsat-ts);
+disp("kW/m",Q/1000,"hear transfer rate per unit width =")
-- 
cgit