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

554 Appendix

Example 6.5 Pinch method for H2C2_150 (MatLab code)

% Example 6.5 Pinch method for H2C2_150

clc
clear

% n_opt = 0: rating
% 1: fmincon
n_opt = 0;
% minimum temperature difference
dtm = 10;
lb = 0.1;
ub = 100;
% subject to the linear inequalities A∗x <= b. If no inequalities exist,
% set A = [] and b = [].
A = [];
b = [];
% subject to the linear equalities Aeq∗x = beq. If no equalities exist,
% set Aeq = [] and beq = [].
Aeq = [];
beq = [];

switch n_opt
case 1
dtm = fmincon(@(dtm)pinch_method_for_H2C2_150(dtm), ...
dtm, A, b, Aeq, beq, lb, ub);
otherwise
pinch_method_for_H2C2_150(dtm);
end
function f = pinch_method_for_H2C2_150 (dtm)
% supply and target temperature and thermal capacity rates
N_h = 2;
th_in = [150, 170];
th_out = [ 50, 40];
C_h = [200, 100];
N_c = 2;
tc_in = [ 50, 80];
tc_out = [120, 110];
C_c = [300, 500];