blob: b35245a0455a63d11b51d1f56db940c305547ea5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
//Variable declaration:
T1 = 99.0 //Mean film temperature (°C)
T2 = 98.0 //Plate surface temperature (°C)
g = 9.807 //Gravitational acceleration (m/s^2)
//From Appendix:
T3 = 100.0 //Saturation temperatre (°C)
h_vap1 = 970.3 //Latent heat of steam in Btu/lb (Btu/lb)
h_vap2 = 2.255*10**6 //Latent heat of steam in J/kg (J/kg)
p_v = 0.577 //Density of steam (kg/m^3)
p_l = 960.0 //Density of liquid water condensate (kg/m^3)
mu_l = 2.82*10**-4 //Absolute viscosity of liquid water condensate (kg/m.s)
k = 0.68 //Thermal conductivity of water (W/m.K)
//From table 12.2
Z = 0.4 //Height of rectangular plate (m)
Pw = 0.2 //Wetted perimeter of rectangular plate (m)
syms h //Average heat transfer coefficient (W/m^2.K)
//Calculation:
A = Z*Pw //Heat transfer area of plate (m^2)
R = A/Pw //Ratio A/Pw (m)
v_l = mu_l/p_l //Kinematic viscosity of liquid water condensate (m^2/s)
Co1 = (h/k)*(v_l**2/g/(1-p_v/p_l))**(1/3) //Condensation number (in terms of the average heat transfer coefficient)
Re = 4*h*Z*(T3-T2)/(mu_l*h_vap2) //Reynolds number in terms of the average heat transfer coefficient
//From equation 12.14:
CO1 = 0.0077*Re**Z //Co in terms of Reynolds number for flow type 1
x1 = solve(h,Co1-CO1) //Solving heat transfer coefficient (W/m^2.K)
h1 =x1(2); //Average heat transfer coefficient for flow type 1 (W/m^2.K)
Re1 = subst(h1,h,Re) //Reynolds number for flow type 1
CO2 = 1.874*Re**(-1/3) //Co in terms of Reynolds number for flow tupe 2
x2 = solve(Co1-CO2,h) //Solving average heat transfer coefficient for flow type 2 (W/m^2.K)
h2 = x2(1); //Average heat transfer coefficient for flow type 2 (W/m^2.K)
Re2 = subst(h2,h,Re) //Reynolds number for flow type 2
h2 = round(h2*10**-1)/10**-1
//Result:
printf("The type of condensation flow type 2 is laminar.")
disp("And the condensation heat transfer coefficient is : ")
disp(h2)
disp("W/m^2.K.")
|