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

using "pressure_changers/valve";
using "streams";
using "controllers/PIDs";
using "stage_separators/flash";

FlowSheet sample_flash
	PARAMETERS
	PP as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = ["n-butane","benzene","n-octane"],
		LiquidModel = "SRK",
		VapourModel = "SRK"
	);
	NComp as Integer;
	SET
	NComp = PP.NumberOfComponents;

	DEVICES
	LCV as valve_flow;

	SET
	LCV.MinFlow = 0 * 'kmol/h';
	LCV.MaxFlow = 700 * 'kmol/h';

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	Feed as simple_source;

	SET

	SPECIFY
	Feed.MolarComposition(1) = 0.3;
	Feed.MolarComposition(2) = 0.3;
	Feed.MolarComposition(3) = 0.4;
	Feed.T = 338 * 'K';
	Feed.P = 5 * 'atm';

	INITIAL

	GUESS

	DEVICES
	PCV as valve_flow;

	SET
	PCV.MinFlow = 0 * 'kmol/h';
	PCV.MaxFlow = 500 * 'kmol/h';

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	VaporProduct as sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	LiquidProduct as sink;

	SET

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	FeedValve as valve_flow;

	SET
	FeedValve.MinFlow = 0 * 'kmol/h';
	FeedValve.MaxFlow = 1000 * 'kmol/h';

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	FeedStep as StepSignal;

	SET
	FeedStep.StepTime = 1.5 * 'h';
	FeedStep.StartValue = 0.5;
	FeedStep.FinalValue = 0.6;

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	PID_PC as PID;

	SET
	PID_PC.PID_Select = "Parallel_AWBT";
	PID_PC.Action = "Direct";
	PID_PC.Mode = "Automatic";
	PID_PC.Clip = "Clipped";
	PID_PC.alpha = 1;
	PID_PC.beta = 1;
	PID_PC.bias = 0.5;
	PID_PC.derivTime = 18 * 's';
	PID_PC.intTime = 72 * 's';
	PID_PC.gain = 0.9275;
	PID_PC.gamma = 1;
	PID_PC.tau = 1 * 's';
	PID_PC.tauSet = 1 * 's';
	PID_PC.MinInput = 0;
	PID_PC.MaxInput = 6;
	PID_PC.MinOutput = 0;
	PID_PC.MaxOutput = 1;

	SPECIFY
	PID_PC.SetPoint = 0.6;

	INITIAL

	GUESS

	DEVICES
	PID_LC as PID;

	SET
	PID_LC.PID_Select = "Parallel_AWBT";
	PID_LC.Action = "Direct";
	PID_LC.Mode = "Automatic";
	PID_LC.Clip = "Clipped";
	PID_LC.alpha = 1;
	PID_LC.beta = 1;
	PID_LC.bias = 0.5;
	PID_LC.derivTime = 19 * 's';
	PID_LC.intTime = 76 * 's';
	PID_LC.gain = 1.1776;
	PID_LC.gamma = 1;
	PID_LC.tau = 1 * 's';
	PID_LC.tauSet = 1 * 's';
	PID_LC.MinInput = 0;
	PID_LC.MaxInput = 2;
	PID_LC.MinOutput = 0;
	PID_LC.MaxOutput = 1;

	SPECIFY
	PID_LC.SetPoint = 0.7;

	INITIAL

	GUESS

	DEVICES
	F101 as flash;

	SET
	F101.Orientation = "vertical";
	F101.Heads = "hemispherical";
	F101.Diameter = 2 * 'm';
	F101.Lenght = 5 * 'm';
	F101.Levelpercent_Initial = 0.80;
	F101.Temperature_Initial = 373 * 'K';
	F101.Composition_Initial(1) = 0.30;
	F101.Composition_Initial(2) = 0.40;
	F101.Composition_Initial(3) = 0.30;

	SPECIFY

	INITIAL

	GUESS

	DEVICES
	heat_duty as energy_source2;

	SET

	SPECIFY
	heat_duty.OutletQ = 1026.32 * 'kW';

	INITIAL

	GUESS

	CONNECTIONS
	PCV.Outlet to VaporProduct.Inlet;
	LCV.Outlet to LiquidProduct.Inlet;
	Feed.Outlet to FeedValve.Inlet;
	FeedStep.OutSignal to FeedValve.FlowFraction;
	PID_LC.Output to LCV.FlowFraction;
	PID_PC.Output to PCV.FlowFraction;
	FeedValve.Outlet to F101.Inlet;
	heat_duty.OutletQ to F101.InletQ;
	F101.OutletLiquid to LCV.Inlet;
	F101.OutletVapour to PCV.Inlet;
	F101.LI to PID_LC.Input;
	F101.PI to PID_PC.Input;

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