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

500 Appendix

for i=1 :1:nx1
t2(1, i, 1) = t2(2, i, 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 15
% analytical
nu_2ab = special_function_nu_(2 ∗ aA, bA);
K2 = special_function_K(bB);
epsilon_h = (1 - bA ∗ K2 / 2) ∗ ( 1 - nu_2ab);

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

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

ntuh = NTUh(2);
ntuc = NTUc(2);
th(:, :) = thm;
for i=1 :1:nx1
tc(i, 1) = tc(nx1 -i+1, ny1);
end
[th, tc, thm, tcm] = crossflow_mixed_unmixed ...
(ntuh, ntuc, nx, ny, th, tc, A, T, D);

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

case 16
% analytical
nu_1 = special_function_nu_(aA, bA);