#*----------------------------------------------
* FlowSheet generated automaticaly by EMSO-GUI
*----------------------------------------------*#

using "streams";
using "heat_exchangers/PHE";

FlowSheet PHE_diagram
	PARAMETERS
	PP as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = ["toluene","benzene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp as Integer;
	SET
	NComp = PP.NumberOfComponents;

	DEVICES
	simple_sink_1 as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS
	simple_sink_1.Inlet.T = 350 * 'K';
	simple_sink_1.Inlet.P = 2 * 'atm';

	DEVICES
	simple_sink_2 as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS
	simple_sink_2.Inlet.T = 300 * 'K';
	simple_sink_2.Inlet.P = 2 * 'atm';

	DEVICES
	PHE_1 as PHE;

	SET
	PHE_1.ChevronAngle = "A30_Deg";
	PHE_1.SideOne = "cold";
	PHE_1.Geometry.Lv = 155 * 'cm';
	PHE_1.Geometry.Nplates = 129;
	PHE_1.Geometry.NpassHot = 2;
	PHE_1.Geometry.NpassCold = 2;
	PHE_1.Geometry.Dports = 15 * 'cm';
	PHE_1.Geometry.Lw = 70 * 'cm';
	PHE_1.Geometry.pitch = 3.7 * 'mm';
	PHE_1.Geometry.pt = 0.70 * 'mm';
	PHE_1.Geometry.Kwall = 17 * 'W/m/K';
	PHE_1.Geometry.Rfh = 0 * 'm^2*K/kW';
	PHE_1.Geometry.Rfc = 0 * 'm^2*K/kW';
	PHE_1.Geometry.PhiFactor = 1.17;

	SPECIFY

	INITIAL

	GUESS
	PHE_1.InletHot.z(1) = 1;
	PHE_1.InletHot.z(2) = 0;
	PHE_1.InletCold.z(1) = 0;
	PHE_1.InletCold.z(2) = 1;

	DEVICES
	source_1 as source;

	SET
	source_1.ValidPhases = "Vapour-Liquid";

	SPECIFY
	source_1.Composition(1) = 0;
	source_1.Composition(2) = 1;
	source_1.F = 290 * 'kmol/h';
	source_1.T = 293 * 'K';
	source_1.P = 2 * 'atm';

	INITIAL

	GUESS

	DEVICES
	source_2 as source;

	SET
	source_2.ValidPhases = "Vapour-Liquid";

	SPECIFY
	source_2.Composition(1) = 1;
	source_2.Composition(2) = 0;
	source_2.F = 82.5 * 'kmol/h';
	source_2.T = 360 * 'K';
	source_2.P = 2 * 'atm';

	INITIAL

	GUESS

	CONNECTIONS
	PHE_1.OutletCold to simple_sink_2.Inlet;
	PHE_1.OutletHot to simple_sink_1.Inlet;
	source_1.Outlet to PHE_1.InletCold;
	source_2.Outlet to PHE_1.InletHot;

	OPTIONS
	Dynamic = false;
	Integration = "index1";
	NLASolver(
		File = "sundials",
		RelativeAccuracy = 1e-3,
		AbsoluteAccuracy = 1e-6,
		MaxIterations = 200
	);
	DAESolver(
		File = "",
		RelativeAccuracy = 1e-3,
		AbsoluteAccuracy = 1e-6,
		EventAccuracy = 1e-2
	);
end