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

474 Appendix

C(i) = C_C2 / M_C2;
G1(i, 3) = 1;
G2(3, i) = 1 / M_C2;
end
else % odd channel number for hot fluid
x_in(i) = L;
x_out(i) = 0;
C(i) = C_H1 / M_H1;
G1(i, 1) = 1;
G2(1, i) = 1 / M_H1;
end
if (i < M)
A(i, i + 1) = U_L / C(i);
A(i, i) = A(i, i) - A(i, i + 1);
end
if (i > 1)
A(i, i - 1) = U_L / C(i);
A(i, i) = A(i, i) - A(i, i - 1);
end
end
V_in = zeros(M, M);
V_out = zeros(M, M);
[A_V, A_D] = eig(A);
for i=1:M
for j=1:M
V_in(i, j) = A_V(i, j) ∗ exp(A_D(j, j) ∗ x_in(i));
V_out(i, j) = A_V(i, j) ∗ exp(A_D(j, j) ∗ x_out(i));
end
end
T_out = G2 ∗ V_out / (V_in - G ∗ V_out) ∗ G1 ∗ T_in;
f = x(1) + x(2) + 1E5 ∗ (max(T_out(1) - T_tar(1), 0) ...
+ max(T_tar(2) - T_out(2), 0) + max(T_tar(3) - T_out(3), 0));
if (f < f0)
f0 = f;
fprintf('f=%f,x=%f,%f,t_out=%f,%f,%f,dt_out=%f,%f,%f\n',...
f, x(1), x(2), T_out(1), T_out(2), T_out(3), ...
T_out(1)-T_tar(1),T_tar(2)-T_out(2),T_tar(3)-T_out(3));
end
end