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

Appendix 547

N_SN = 2 ∗ (N_h + N_c);
for i=1:1: N_SN - 1
for j =i+1:1: N_SN
if (t(j) == t(i))
for k = j : 1 : N_SN - 1
t(k) = t(k + 1);
end
N_SN = N_SN - 1;
end
end
end
for i=1:1: N_SN - 1
for j =i+1:1: N_SN
if (t(j) > t(i))
s = t(i);
t(i) = t(j);
t(j) = s;
end
end
end
N_SN = N_SN - 1;

% total enthalpy change
dH_h = zeros(N_SN, 1);
dH_c = zeros(N_SN, 1);
for i=1:1: N_SN
C_h_i = 0;
for j =1:1:N_h
if (th_out(j) <= t(i + 1) && t(i) <= th_in(j))
C_h_i = C_h_i + C_h(j);
end
end
dH_h(i) = (t(i) - t(i + 1)) ∗ C_h_i;
C_c_i = 0;
for j =1:1:N_c
if (tc_s_in(j) <= t(i + 1) && t(i) <= tc_s_out(j))
C_c_i = C_c_i + C_c(j);
end
end