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

using "streams.mso";
using "heat_exchangers\PHE.mso";

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;
	simple_sink_2 as simple_sink;
	PHE_1 as PHE;
	source_1 as simple_source;
	source_2 as simple_source;

	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;
#	source_1.Outlet to simple_sink_2.Inlet;
#	source_2.Outlet to simple_sink_1.Inlet;

	SET
	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 ;
	PHE_1.ChevronAngle = "A30_Deg" ;
	PHE_1.SideOne = "cold" ;
	source_1.ValidPhases = "Liquid-Only" ;
	source_2.ValidPhases = "Liquid-Only" ;

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

	INITIAL

	GUESS
	simple_sink_1.Inlet.T = 350 * 'K' ;
	simple_sink_1.Inlet.P = 2 * 'atm' ;
	simple_sink_2.Inlet.T = 300 * 'K' ;
	simple_sink_2.Inlet.P = 2 * 'atm' ;
	PHE_1.InletHot.z(1) = 1 ;
	PHE_1.InletHot.z(2) = 0 ;
	PHE_1.InletCold.z(1) = 0 ;
	PHE_1.InletCold.z(2) = 1 ;
	PHE_1.OutletCold.T = 300 *'K';
	PHE_1.OutletHot.T = 350 *'K';
	PHE_1.OutletCold.P = 2 *'atm';
	PHE_1.OutletHot.P = 2 *'atm';
	PHE_1.OutletHot.z(1) = 1 ;
	PHE_1.OutletHot.z(2) = 0 ;
	PHE_1.OutletCold.z(1) = 0 ;
	PHE_1.OutletCold.z(2) = 1 ;
	PHE_1.ColdSide.Properties.Outlet.rho=900*'kg/m^3';
	PHE_1.HotSide.Properties.Outlet.rho=800*'kg/m^3';
	PHE_1.ColdSide.Properties.Average.rho=900*'kg/m^3';
	PHE_1.HotSide.Properties.Average.rho=800*'kg/m^3';

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

end