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

using "streams";
using "stage_separators/column";
using "stage_separators/reboiler";

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

	DEVICES
	feed as simple_source;

	SET

	SPECIFY
	feed.MolarComposition(1) = 0.5;
	feed.MolarComposition(2) = 0.5;
	feed.F = 113.4 * 'kmol/h';
	feed.T = 291 * 'K';
	feed.P = 168.30 * 'kPa';

	INITIAL

	GUESS

	DEVICES
	SectionBottom as Section_Column;

	SET
	SectionBottom.NumberOfTrays = 8;
	SectionBottom.FeedTrayLocation = 5;
	SectionBottom.VolumeOfTray = 4 * 'ft^3';
	SectionBottom.HeatSupply = 0 * 'kW';
	SectionBottom.PlateArea = 3.94 * 'ft^2';
	SectionBottom.HolesArea = 0.394 * 'ft^2';
	SectionBottom.WeirLength = 20.94 * 'in';
	SectionBottom.WeirHeight = 0.125 * 'ft';
	SectionBottom.beta = 0.6;
	SectionBottom.alfa = 4;

	SPECIFY
	SectionBottom.VapourDrawOffFlow = 0 * 'kmol/h';
	SectionBottom.LiquidDrawOffFlow = 0 * 'kmol/h';
	SectionBottom.MurphreeEff = 1;

	INITIAL

	GUESS

	DEVICES
	HeatSupplied as energy_source;

	SET

	SPECIFY
	HeatSupplied.OutletQ = 3.774 * 'kJ/h';

	INITIAL

	GUESS

	DEVICES
	LiquidOut as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	SectionTop as Section_Column;

	SET
	SectionTop.NumberOfTrays = 4;
	SectionTop.VolumeOfTray = 4 * 'ft^3';
	SectionTop.HeatSupply = 0 * 'kW';
	SectionTop.PlateArea = 3.94 * 'ft^2';
	SectionTop.HolesArea = 0.394 * 'ft^2';
	SectionTop.WeirLength = 21 * 'in';
	SectionTop.WeirHeight = 0.13 * 'ft';
	SectionTop.beta = 0.6;
	SectionTop.alfa = 4;
	SectionTop.INITIALIZATION.TopTemperature = 300 * 'K';
	SectionTop.INITIALIZATION.BottomTemperature = 330 * 'K';
	SectionTop.INITIALIZATION.LevelFraction = 0.5;
	SectionTop.INITIALIZATION.TopComposition(1) = 0.5;
	SectionTop.INITIALIZATION.TopComposition(2) = 0.5;
	SectionTop.INITIALIZATION.BottomComposition(1) = 0.5;
	SectionTop.INITIALIZATION.BottomComposition(2) = 0.5;

	SPECIFY
	SectionTop.VapourDrawOffFlow = 0 * 'kmol/h';
	SectionTop.LiquidDrawOffFlow = 0 * 'kmol/h';
	SectionTop.MurphreeEff = 1;
	SectionTop.VapourOutlet.F = 12 * 'kmol/h';

	INITIAL

	GUESS

	DEVICES
	vaporoutlet as sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	CondenserDevice as simple_source;

	SET

	SPECIFY
	CondenserDevice.MolarComposition(1) = 0.66;
	CondenserDevice.MolarComposition(2) = 0.33;
	CondenserDevice.F = 68 * 'kmol/h';
	CondenserDevice.T = 280 * 'K';
	CondenserDevice.P = 150 * 'kPa';

	INITIAL

	GUESS

	DEVICES
	NoFlow as sourceNoFlow;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	reboilerDevice as reboiler;

	SET
	reboilerDevice.Across = 1 * 'm^2';
	reboilerDevice.V = 2 * 'm^3';
	reboilerDevice.Initial_Level = 1 * 'm';
	reboilerDevice.Initial_Temperature = 350 * 'K';
	reboilerDevice.Initial_Composition(1) = 0.5;
	reboilerDevice.Initial_Composition(2) = 0.5;

	SPECIFY
	reboilerDevice.OutletLiquid.F = 28.4 * 'kmol/h';

	INITIAL

	GUESS

	CONNECTIONS
	feed.Outlet to SectionBottom.FeedTray;
	SectionBottom.VapourOutlet to SectionTop.VapourInlet;
	SectionTop.LiquidOutlet to SectionBottom.LiquidInlet;
	SectionTop.VapourOutlet to vaporoutlet.Inlet;
	CondenserDevice.Outlet to SectionTop.LiquidInlet;
	NoFlow.Outlet to SectionTop.FeedTray;
	HeatSupplied.OutletQ to reboilerDevice.InletQ;
	reboilerDevice.OutletVapour to SectionBottom.VapourInlet;
	SectionBottom.LiquidOutlet to reboilerDevice.InletLiquid;
	reboilerDevice.OutletLiquid to LiquidOut.Inlet;

	OPTIONS
	Dynamic = true;
	TimeStep = 0.1;
	TimeEnd = 100;
	TimeUnit = 's';
	Integration = "original";
	NLASolver(
		File = "sundials",
		RelativeAccuracy = 1e-3,
		AbsoluteAccuracy = 1e-6,
		MaxIterations = 100
	);
	DAESolver(
		File = "sundials",
		RelativeAccuracy = 1e-3,
		AbsoluteAccuracy = 1e-6,
		EventAccuracy = 1e-2
	);
end