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

Appendix 507

for j=1:1:ny1
th(1, j) = th(nx1, ny1 -j+1);
end
[th, tc, thm, tcm] = crossflow_unmixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc);

fprintf("thm =%f, tcm =%f, 1-thm =%f, e_h =%f\n", ...
thm, tcm, 1 - thm, epsilon_h);
case 66
% analytical
nu1 = special_function_nu(aA, bA);
nu2ab = special_function_nu(2 ∗ aA, bA);
nua2b = special_function_nu(aA, 2 ∗ bA);
epsilon_h = 1 - nu2ab - 2 ∗ (nua2b - nu1);

% numerical
th(:,:) = 1;
tc(:,:) = 0;

ntuh = NTUh(1);
ntuc = NTUc(1);
[th, tc, thm, tcm] = crossflow_unmixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc);
fprintf("thm_i =%f, tcm_i =%f\n", thm, tcm);

ntuh = NTUh(2);
ntuc = NTUc(2);
for i=1:1:nx1
tc(i, 1) = tc(i, ny1);
end
for j=1:1:ny1
th(1, j) = th(nx1, j);
end
[th, tc, thm, tcm] = crossflow_unmixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc);

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

otherwise
fprintf("not available\n");
end