-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathtemp_K.m
151 lines (130 loc) · 5.47 KB
/
temp_K.m
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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
% Kandlikar temperature computation
function [Tp_k Tc_k h]=temp_K
global Tsat Psat rhol_R245FA rhog_R245FA Hl_R245FA Hg_R245FA Cpl_R245FA Cpg_R245FA mul_R245FA mug_R245FA lambdal_R245FA lambdag_R245FA sigma_R245FA
global R_contact
global rhol rhog mul mug hlg lambda Cpl Cpg D G Tf lmt
global x Zz P ggeom qgeom N Q_geom L_geom l_geom Np_geom
Tp_k=ones(length(x),1);
Tc_k=ones(length(x),1);
qq=qgeom(1);
gg=abs(ggeom(1));
for j=1:N(1)
for i=2:length(Psat)
if (P(j)<=Psat(i)) && (P(j)>Psat(i-1))
rhol=(rhol_R245FA(i)-rhol_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+rhol_R245FA(i);
rhog=(rhog_R245FA(i)-rhog_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+rhog_R245FA(i);
mul=(mul_R245FA(i)-mul_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+mul_R245FA(i);
mug=(mug_R245FA(i)-mug_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+mug_R245FA(i);
Hl=(Hl_R245FA(i)-Hl_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Hl_R245FA(i);
Hg=(Hg_R245FA(i)-Hg_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Hg_R245FA(i);
lambda=(lambdal_R245FA(i)-lambdal_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+lambdal_R245FA(i);
Cpl=(Cpl_R245FA(i)-Cpl_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Cpl_R245FA(i);
Cpg=(Cpg_R245FA(i)-Cpg_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Cpg_R245FA(i);
Tf=(Tsat(i)-Tsat(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Tsat(i);
hlg=Hg-Hl;
end
end
if x(j)<lmt
C0=((1-x(j))/x(j))^0.8*sqrt(rhog/rhol);
if (C0<0.65)
C1=1.1360;
C2=-0.9;
C3=667.2;
C4=0.7;
C5=0.3;
else
C1=0.6683;
C2=-0.2;
C3=1058;
C4=0.7;
C5=0.3;
end
Bo=qq/(G*hlg);
Fr=G^2/(rhol^2*gg*D);
if Fr>0.04
C5=0;
end
Pr=Cpl*mul/lambda;
Fk=1.4;
hl=lambda/D*0.023*(G*(1-x(j))*D/mul)^0.8*Pr^(1/3);
h(j)=hl*(C1*C0^C2*(25*Fr)^C5+C3*Bo^C4*Fk);
Tp_k(j)=Tf+qq/(h(j));
Tc_k(j)=Tp_k(j)+1/Conductance(h(j))*Q_geom(1)/(L_geom(1)*l_geom(1))*pi*D*3*Np_geom(1)/l_geom(1)+1/R_contact*Q_geom(1)/(L_geom(1)*l_geom(1));
if qq==0
h(j)=0;
Tc_k(j)=Tf;
end
else
Pr=Cpl*mul/lambda;
hl=0.023*lambda/D*(G*D/mul)^0.8*Pr^(1/3);
h(j)=hl;
Tp_k(j)=Tf+qq*(4/(D*G*Cpg)*(Zz(j)-Zz(1))+1/hl);
Tc_k(j)=Tp_k(j)+1/Conductance(h(j))*Q_geom(1)/(L_geom(1)*l_geom(1))*pi*D*3*Np_geom(1)/l_geom(1)+1/R_contact*Q_geom(1)/(L_geom(1)*l_geom(1));
if qq==0
h(j)=0;
Tc_k(j)=Tf;
end
end
end
for k=2:length(N)
qq=qgeom(k);
gg=abs(ggeom(k));
for j=sum(N(1:k-1))+1:sum(N(1:k))
for i=2:length(Psat)
if (P(j)<=Psat(i)) && (P(j)>Psat(i-1))
rhol=(rhol_R245FA(i)-rhol_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+rhol_R245FA(i);
rhog=(rhog_R245FA(i)-rhog_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+rhog_R245FA(i);
mul=(mul_R245FA(i)-mul_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+mul_R245FA(i);
mug=(mug_R245FA(i)-mug_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+mug_R245FA(i);
Hl=(Hl_R245FA(i)-Hl_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Hl_R245FA(i);
Hg=(Hg_R245FA(i)-Hg_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Hg_R245FA(i);
lambda=(lambdal_R245FA(i)-lambdal_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+lambdal_R245FA(i);
Cpl=(Cpl_R245FA(i)-Cpl_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Cpl_R245FA(i);
Cpg=(Cpg_R245FA(i)-Cpg_R245FA(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Cpg_R245FA(i);
Tf=(Tsat(i)-Tsat(i-1))/(Psat(i)-Psat(i-1))*(P(j)-Psat(i))+Tsat(i);
hlg=Hg-Hl;
end
end
if x(j)<lmt
C0=(((1-x(j))/x(j))^0.8)*sqrt(rhog/rhol);
if (C0<0.65)
C1=1.1360;
C2=-0.9;
C3=667.2;
C4=0.7;
C5=0.3;
else
C1=0.6683;
C2=-0.2;
C3=1058;
C4=0.7;
C5=0.3;
end
Bo=qq/(G*hlg);
Fr=G^2/(rhol^2*gg*D);
if Fr>0.04
C5=0;
end
Pr=Cpl*mul/lambda;
Fk=1.4;
hl=lambda/D*0.023*(G*(1-x(j))*D/mul)^0.8*Pr^(1/3);
h(j)=hl*(C1*(C0^C2)*(25*Fr)^C5+C3*(Bo^C4)*Fk);
Tp_k(j)=Tf+qq/(h(j));
Tc_k(j)=Tp_k(j)+1/Conductance(h(j))*Q_geom(1)/(L_geom(k)*l_geom(k))*pi*D*3*Np_geom(k)/l_geom(k)+1/R_contact*Q_geom(k)/(L_geom(k)*l_geom(k));
if qq==0
h(j)=0;
Tc_k(j)=Tf;
end
else
Pr=Cpl*mul/lambda;
hl=0.023*lambda/D*(G*D/mul)^0.8*Pr^(1/3);
h(j)=hl;
Tp_k(j)=Tf+qq*(4/(D*G*Cpg)*(Zz(j)-Zz(sum(N(1:k-1))))+1/hl);
Tc_k(j)=Tp_k(j)+1/Conductance(h(j))*Q_geom(1)/(L_geom(k)*l_geom(k))*pi*D*3*Np_geom(k)/l_geom(k)+1/R_contact*Q_geom(k)/(L_geom(k)*l_geom(k));
if qq==0
h(j)=0;
Tc_k(j)=Tf;
end
end
end
end