using "stage_separators/flash";

# exemplo baseado nos dados do artigo do Gani 1985.
FlowSheet flash_Test
	PARAMETERS
	PP as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = ["1,3-butadiene", "isobutene", "n-pentane",
			"1-pentene", "1-hexene", "benzene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp as Integer;

	VARIABLES
	Q	as energy_source (Brief="Heat supplied");
	
	SET
	NComp = PP.NumberOfComponents;
	
	DEVICES
	fl as flash_steady;
	s1 as source;
	
	CONNECTIONS
	s1.Outlet to fl.Inlet;
	Q.OutletQ to fl.InletQ;
	
	SPECIFY
	s1.Outlet.F = 496.3 * 'kmol/h';
	s1.Outlet.T = 338 * 'K';
	s1.Outlet.P = 507.1 * 'kPa';
	s1.Outlet.z = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283];
	
	fl.OutletL.P = 2.5 * 'atm';

	#Q.OutletQ.Q = 0 * 'kJ/h';
	fl.OutletL.T = 315.06 * 'K';
	
	OPTIONS
	Dynamic = false;
end

Reconciliation flash_rec as flash_Test

	RECONCILE
	s1.Outlet.F; fl.OutletL.F; fl.OutletV.F;
	
	FREE
	s1.Outlet.F;

	EXPERIMENTS
	# FILE            WEIGTH
	"flash.dat"       1;

	OPTIONS
	Filter = "mean";
	Significance = 0.95;
	ObjectiveFunction = "Normal";
	
	GrossErrorTests(
		Global = true,
		Nodal = true,
		Measurements = true
	);
	
	NLPSolver(
		MaxIterations=1000,
		File = "complex"
		#File = "ipopt_emso"
	);
	Dynamic = false;
end