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

using "streams";
using "stage_separators/column";
using "pressure_changers/pump";

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

	DEVICES
	feed as simple_source;

	SET

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

	INITIAL

	GUESS

	DEVICES
	vapourOutlet as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	liquidOut as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	LiquidBottom as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	ToAnotherColumn as simple_sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	energy_source_1 as energy_source;

	SET

	SPECIFY
	energy_source_1.OutletQ = 3.77 * 'kW';

	INITIAL

	GUESS

	DEVICES
	energy_source_2 as energy_source;

	SET

	SPECIFY
	energy_source_2.OutletQ = -4E6 * 'kJ/h';

	INITIAL

	GUESS

	DEVICES
	Distillation_kettle_cond_1 as Distillation_kettle_cond;

	SET
	Distillation_kettle_cond_1.VolumeOfTray = 0.113267 * 'm^3';
	Distillation_kettle_cond_1.HeatSupply = 0 * 'kW';
	Distillation_kettle_cond_1.PlateArea = 0.366038 * 'm^2';
	Distillation_kettle_cond_1.HolesArea = 0.0366038 * 'm^2';
	Distillation_kettle_cond_1.WeirLength = 0.531876 * 'm';
	Distillation_kettle_cond_1.WeirHeight = 0.0381 * 'm';
	Distillation_kettle_cond_1.beta = 0.6;
	Distillation_kettle_cond_1.alfa = 4;
	Distillation_kettle_cond_1.CONDENSER.V = 2 * 'm^3';
	Distillation_kettle_cond_1.CONDENSER.Across = 1 * 'm^2';
	Distillation_kettle_cond_1.CONDENSER.Initial_Level = 1 * 'm';
	Distillation_kettle_cond_1.CONDENSER.Initial_Temperature = 260 * 'K';

	SPECIFY
	Distillation_kettle_cond_1.VapourDrawOffFlow = 0 * 'kmol/h';
	Distillation_kettle_cond_1.LiquidDrawOffFlow = 1 * 'kmol/h';
	Distillation_kettle_cond_1.MurphreeEff = 1;
	Distillation_kettle_cond_1.alfaTopo = 2;
	Distillation_kettle_cond_1.CONDENSER.OutletVapour.F = 0 * 'kmol/h';
	Distillation_kettle_cond_1.REBOILER.OutletLiquid.F = 28 * 'kmol/h';
	Distillation_kettle_cond_1.SPLITTER.FlowRatios(1) = 0.4;
	Distillation_kettle_cond_1.SPLITTER.Outlet2.F = 85 * 'kmol/h';
	Distillation_kettle_cond_1.PUMP.dP = 16 * 'kPa';

	INITIAL

	GUESS

	DEVICES
	Pump_Device as centrifugal_pump;

	SET
	Pump_Device.PumpEfficiency = 0.80;
	Pump_Device.MechanicalEff = 0.95;
	Pump_Device.NPSH_Options = "Default";
	Pump_Device.SuctionArea = 0.01 * 'm^2';

	SPECIFY
	Pump_Device.Pratio = 7;
	Pump_Device.StaticHead = 12 * 'm';

	INITIAL

	GUESS

	DEVICES
	work_inlet_1 as work_source;

	SET

	SPECIFY

	INITIAL

	GUESS

	CONNECTIONS
	feed.Outlet to Distillation_kettle_cond_1.FeedTray;
	energy_source_1.OutletQ to Distillation_kettle_cond_1.HeatToReboiler;
	Distillation_kettle_cond_1.BottomProduct to LiquidBottom.Inlet;
	Distillation_kettle_cond_1.LiquidDistillate to liquidOut.Inlet;
	energy_source_2.OutletQ to Distillation_kettle_cond_1.HeatToCondenser;
	Distillation_kettle_cond_1.VapourDistillate to vapourOutlet.Inlet;
	Distillation_kettle_cond_1.LiquidDrawOff to Pump_Device.Inlet;
	Pump_Device.Outlet to ToAnotherColumn.Inlet;
	work_inlet_1.Work to Pump_Device.WorkIn;

	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