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

506 Appendix

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(nx1 -i+1, ny1);
end
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 56
% analytical
nub = special_function_nu(bA, aA);
nu2ba = special_function_nu(2 ∗ bA, aA);
mub = special_function_mu(bA, aA);
epsilon_h = (nub - nu2ba + mub) / Rh(1);
% 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