Page 508 - Design and Operation of Heat Exchangers and their Networks

Appendix 491

for i=1:1: nx1
for j=1:1:ny1
s = th(i, j) - t1(2, i, j);
if (abs(s) > abs(smax))
smax = s;
end
t1(2, i, j) = t1(2, i, j) + e ∗ s;
s = tc(i, j) - t2(2, i, j);
if (abs(s) > abs(smax))
smax = s;
end
t2(2, i, j) = t2(2, i, j) + e ∗ s;
end
end
for i=1:1: nx1
t2(1, i, 1) = t2(2, nx1 - i + 1, ny1);
end
n=n+ 1;
fprintf("%d, %e\n", n, smax);
end

fprintf("smax =%e, thm =%f, tcm =%f, 1-thm =%f, e_h =%f\n", ...
smax, thm, tcm, 1 - thm, epsilon_h);

case-46 % BA_4,6
% analytical
nu1 = special_function_nu(aA, bA);
nu2 = special_function_nu(aB, bB);
mu2 = special_function_mu2(aA, bA, bB, phi);
epsilon_h = 1 - Rh(1) ∗ nu1 ∗ nu2 ...
/ (1 - Rh(1) ∗ (1 - nu1 - nu2) - mu2);
% numerical
t1(:, :, :) = 1;
t2(:, :, :) = 0;
smax = 100;
n=0;
while (abs(smax) > eps && n < 100)
smax = 0;
ntuh = NTUh(1);