F_c_max = 10;        % maximum cold fluid mass flow, kg/s
F_h_max = 3;         % maximum hot fluid mass flow, kg/s
u10=0.5;             % initial opening for cold valve, [-]
u20= 1;              % initial opening for hot valve, [-]
rho_c = 980 ;        % cold fluid density, kg/m3
rho_h = 890;         % hot fluid density,  kg/m3
U = 150*1e-3;        % heat transfer coefficient,   kW/m2/K
A = 90/10;           % heat transfer area, m2
cp_c = 1.5;          % specific heat for cold fluid, kJ/kg/K
cp_h = 1.2;          % specific heat for hot fluid,  kJ/kg/K
Tmax = (70+273.15);  % hot flid inlet temp = maximum temperature. will be used for scalling all other states
Tinh = (70+273.15)/Tmax;    % parameter
Tinc0 = (20 + 273.15)/Tmax; % cold fluid inlet temperature
b = 0.5;             % volume fraction,taken as equal for both shell and tube sides
Aspec = 100;         % specific area to volume, taken from mathisen.phd
V_c = A*b/Aspec;     % volume cold side, from mathisen.phd
V_h = A*b/Aspec;     % volume hot side, from mathisen.phd