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

Appendix 551

% problem table
fprintf("SN, t_h, dH_c, dH_h, D, I, O, H_h, t_c, H_c\n");
for i=1:1: N_SN
if (i < j1 jj i > j2)
tth = 0;
else
tth = t(i);
end
if (i < k1 jj i > k2)
ttc = 0;
else
ttc = t(i) - dtm;
end
fprintf("%d, %f, %f, %f, %f, %f, %f, %f, %f, %f\n", ...
i, tth, dH_c(i), dH_h(i), D(i), I(i), O(i), H_h(i), ttc, H_c(i));
end
if (N_SN +1<j1 jj N_SN +1>j2)
tth = 0;
else
tth = t(N_SN + 1);
end
if (N_SN +1<k1 jj N_SN +1>k2)
ttc = 0;
else
ttc = t(N_SN + 1) - dtm;
end
fprintf("%d, %f, –, –, –, –, –, –, %f, –\n", N_SN + 1, tth, ttc);
fprintf("t_h_pinch =%f, Q_HU_min =%f, Q_CU_min =%f\n", ...
t_h_pinch, Q_HU_min, Q_CU_min);

% cost calculation (isothermal mixing)

% heat exchanger costs
AE = zeros(6, 1);
QE = zeros(6, 1);
ChE = zeros(6, 1);
CcE = zeros(6, 1);
thE_in = zeros(6, 1);
thE_out = zeros(6, 1);