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

480 Appendix

s = tc(i, j) - t2(1, i, j);
if (abs(s) > abs(smax))
smax = s;
end
t2(1, i, j) = t2(1, i, j) + e ∗ s;
end
end
t1(2, 1, :) = thm;

ntuh = NTUh(2);
ntuc = NTUc(2);
th(:, :) = t1(2, :, :);
tc(:, :) = t2(2, :, :);

th(1, :) = thm;
tc(:, 1) = t2(2, :, 1);
[th, tc, thm, ] = crossflow_mixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc, A, T, D);
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);