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

478 Appendix

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
t1(3, 1, :) = thm;
for i=1 :1:nx1
t2(1, i, 1) = t2(2, i, ny1);
end

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

[th, tc, thm, ] = crossflow_unmixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc);
for i=1 :1:nx1
for j=1:1:ny1
s = th(i, j) - t1(3, i, j);
if (abs(s) > abs(smax))
smax = s;
end
t1(3, i, j) = t1(3, i, j) + e ∗ s;
s = tc(i, j) - t2(3, i, j);
if (abs(s) > abs(smax))
smax = s;
end
t2(3, i, j) = t2(3, i, j) + e ∗ s;
end
end