diff options
Diffstat (limited to '1073/CH3/EX3.44/3_44.sce')
| -rwxr-xr-x | 1073/CH3/EX3.44/3_44.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/1073/CH3/EX3.44/3_44.sce b/1073/CH3/EX3.44/3_44.sce new file mode 100755 index 000000000..b4036adc6 --- /dev/null +++ b/1073/CH3/EX3.44/3_44.sce @@ -0,0 +1,26 @@ +clc;
+clear;
+//Example 3.44
+rho_l=957.9 //[kg/m^3]
+lambda=2257 //[kJ/kg]
+lambda=lambda*10^3 //[J/kg]
+rho_v=31.54 //[kg/m^3]
+Cpv=4.64 //[kJ/kg.K]
+Cpv=Cpv*10^3 //[J/kg.K]
+kv=58.3*10^-3//[W/(m.K)]
+g=9.81 //[m/s^2]
+mu_v=18.6*10^-6 //[kg/(m.s)]
+e=1.0 //Emissivity
+sigma=5.67*10^-8;
+Ts=373 //[K]
+Tw=628 //[K]
+dT=Tw-Ts //[K]
+D=1.6*10^-3 //[m]
+T=(Tw+Ts)/2 //[K]
+hc=0.62*((kv^3)*rho_v*(rho_l-rho_v)*g*(lambda+0.40*Cpv*dT)/(D*mu_v*dT))^(1.0/4.0)//Convective heat transfer coeff [W/sq m.K]
+hr=e*sigma*(Tw^4-Ts^4)/(Tw-Ts) //Radiation heat transfer coeff in [W/sq m.K]
+h=hc+(3/4)*hr //Total heat transfer coefficient W/(sq m.K)
+Q_by_l=h*%pi*D*dT //Heat dissipation rate per unit length in [kW/m]
+printf("\n Stable film boiling point heat transfer coefficient is %f W/(sq m.K)",h);
+Q_by_l=Q_by_l/1000 //[kW/m]
+printf("\n Heat dissipated per unit length of the heater is %f kW/m",Q_by_l);
|