diff options
Diffstat (limited to '965/CH9/EX9.18/18.sci')
| -rw-r--r-- | 965/CH9/EX9.18/18.sci | 22 |
1 files changed, 22 insertions, 0 deletions
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 =")
|