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

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

FlowSheet ColumnPacked_Diagram
	PARAMETERS
	PP as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = ["isobutane","n-pentane","propylene","benzene","isobutene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp as Integer;

DEVICES
	feed 				as simple_source;
	Packed_Section 		as Packed_Section_Column;
	Vapour_Out 			as simple_sink;
	Liquid_Out 			as simple_sink;
	condenser_device 	as simple_source;
	reboiler_device 	as simple_source;
	
SET
	NComp = PP.NumberOfComponents;
	feed.ValidPhases = "Liquid-Only";

SPECIFY
	feed.MolarComposition(1) = 0.2;
	feed.MolarComposition(2) = 0.2;
	feed.MolarComposition(3) = 0.2;
	feed.MolarComposition(4) = 0.2;
	feed.MolarComposition(5) = 0.2;
	feed.F = 113.4 * 'kmol/h';
	feed.T = 291 * 'K';
	feed.P = 1.66 * 'atm';

SET
	Packed_Section.NumberOfStages = 8;
	Packed_Section.FeedStageLocation = 5;
	Packed_Section.PackingHeight = 4 * 'm';
	Packed_Section.HeatSupply = 0 * 'kW';
	Packed_Section.ColumnDiameter = 1.009 * 'm';
	Packed_Section.VoidFraction = 0.951;
	Packed_Section.ResistanceCoeff = 0.1;
	Packed_Section.AreaPerPackingVol = 112.6 * 'm^2/m^3';
	Packed_Section.INITIALIZATION.TopStageTemperature = 300 * 'K';
	Packed_Section.INITIALIZATION.BottomStageTemperature = 320 * 'K';
	Packed_Section.INITIALIZATION.TopStageComposition(1) = 0.2;
	Packed_Section.INITIALIZATION.TopStageComposition(2) = 0.2;
	Packed_Section.INITIALIZATION.TopStageComposition(3) = 0.2;
	Packed_Section.INITIALIZATION.TopStageComposition(4) = 0.2;
	Packed_Section.INITIALIZATION.TopStageComposition(5) = 0.2;
	Packed_Section.INITIALIZATION.BottomStageComposition(1) = 0.2;
	Packed_Section.INITIALIZATION.BottomStageComposition(2) = 0.2;
	Packed_Section.INITIALIZATION.BottomStageComposition(3) = 0.2;
	Packed_Section.INITIALIZATION.BottomStageComposition(4) = 0.2;
	Packed_Section.INITIALIZATION.BottomStageComposition(5) = 0.2;
	Packed_Section.INITIALIZATION.ColumnHoldup = 0.01;

	SPECIFY
	Packed_Section.PressureDrop = 0.008 * 'atm';


SET
	condenser_device.ValidPhases = "Liquid-Only";

SPECIFY
	condenser_device.MolarComposition(1) = 0.599;
	condenser_device.MolarComposition(2) = 0.044;
	condenser_device.MolarComposition(3) = 0.035;
	condenser_device.MolarComposition(4) = 0.007;
	condenser_device.MolarComposition(5) = 0.315;
	condenser_device.F = 85 * 'kmol/h';
	condenser_device.T = 283.5 * 'K';
#	condenser_device.P = 2.33 * 'atm';

SET
	reboiler_device.ValidPhases = "Vapour-Only";

SPECIFY
	reboiler_device.MolarComposition(1) = 0.16;
	reboiler_device.MolarComposition(2) = 0.542;
	reboiler_device.MolarComposition(3) = 0.013;
	reboiler_device.MolarComposition(4) = 0.008;
	reboiler_device.MolarComposition(5) = 0.277;
	reboiler_device.F = 137.57 * 'kmol/h';
	reboiler_device.T = 325 * 'K';
	reboiler_device.P = 2.46 * 'atm';

CONNECTIONS
	feed.Outlet to Packed_Section.FeedStage;
	Packed_Section.VapourOutlet to Vapour_Out.Inlet;
	Packed_Section.LiquidOutlet to Liquid_Out.Inlet;
	condenser_device.Outlet to Packed_Section.LiquidInlet;
	reboiler_device.Outlet to Packed_Section.VapourInlet;

OPTIONS
	Dynamic = false;
	TimeStep = 0.1;
	TimeEnd = 1;
	TimeUnit = 'h';
	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