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

using "water_steam/power_plant.mso";

FlowSheet planta_cogeracao
	PARAMETERS
	PP as Plugin(Brief="Physical Properties",
		Type="PP"
		,Components = ["water"]
		,LiquidModel = "IdealLiquid"
		,VapourModel = "Ideal"
	);
	NComp as Integer;
	SET
	NComp = PP.NumberOfComponents;

	DEVICES
	GV as Gerador_VaporMod;
	TAP as Turbina_sangra;
	TMP1 as Turbina_sangra;
	TMP2 as Turbina_sangra;
	TBP as Turbina;
	COND as Condensador;
	ATEMP as Tanque2;
	TD1 as Turbina;
	T4P as Tanque4perdas;
	DEAERA as Desaerador5;
	B4 as Bomba;
	T2P as Tanque2;
	B3 as Bomba;
	B5 as Bomba;
	B2 as Bomba;
	CONDP as Condensador;
	B1 as Bomba;
	CONDP6 as Condensador;
	FLASH as Flash;
	SP4 as Splitter4;
	Torre_1 as Torre;
	SUM as SumidouroQ;
	SumidouroQ_1 as SumidouroQ;
	SUM2 as Sumidouro;
	ETA_CICLO_1 as ETA_CICLO;
	GE as Gerador_Eletrico;
	SUMM as SumidouroQ;
	Fonte2_1 as Fonte2;
	Splitter_1 as Splitter;

	CONNECTIONS
	GV.Fout_sa to TAP.Fin;
	TAP.Fouts to ATEMP.Fin1;
	TAP.Fout to TMP1.Fin;
	ATEMP.Fout to TD1.Fin;
	TMP1.Fouts to T2P.Fin1;
	TMP1.Fout to TMP2.Fin;
	B4.Fout to ATEMP.Fin2;
	TD1.Fout to T4P.Fin2;
	B3.Fout to T2P.Fin2;
	TMP2.Fout to TBP.Fin;
	B5.Fout to GV.Fin_a;
	B2.Fout to T4P.Fin4;
	T4P.Fout to CONDP.Fin;
	CONDP.Fout to DEAERA.Fin1;
	COND.Fout to B1.Fin;
	B1.Fout to DEAERA.Fin2;
	TBP.Fout to COND.Fin;
	T2P.Fout to CONDP6.Fin;
	CONDP6.Fout to FLASH.Fin;
	FLASH.FoutL to DEAERA.Fin3;
	FLASH.FoutV to T4P.Fin3;
	DEAERA.Fout to SP4.Fin;
	SP4.Fout1 to B2.Fin;
	SP4.Fout2 to B3.Fin;
	SP4.Fout3 to B4.Fin;
	SP4.Fout4 to B5.Fin;
	COND.Q_COND to Torre_1.Qin;
	CONDP6.Q_COND to SUM.Qin;
	CONDP.Q_COND to SumidouroQ_1.Qin;
	T4P.Fperda to SUM2.Fin;
	GV.Q_GV to ETA_CICLO_1.POT_GV;
	B5.POT_BMB to ETA_CICLO_1.POT_BMB5;
	B4.POT_BMB to ETA_CICLO_1.POT_BMB1;
	B3.POT_BMB to ETA_CICLO_1.POT_BMB2;
	B2.POT_BMB to ETA_CICLO_1.POT_BMB3;
	B1.POT_BMB to ETA_CICLO_1.POT_BMB4;
	TAP.POT_TURB to GE.POT_TURB2;
	TMP1.POT_TURB to GE.POT_TURB1;
	TMP2.POT_TURB to GE.POT_TURB3;
	TBP.POT_TURB to GE.POT_TURB4;
	GE.POT_GE to ETA_CICLO_1.POT_GE;
	TD1.POT_TURB to SUMM.Qin;
	Fonte2_1.Fout to DEAERA.Fin4;
	TMP2.Fouts to Splitter_1.Fin;
	Splitter_1.Fouts to T4P.Fin1;
	Splitter_1.Fout to DEAERA.Fin5;

	SET

	SPECIFY
	GV.EF_GV = 0.8 ;
	GV.Fin_a.F = 81.3 * 'kg/s' ;
	GV.Feco.P = 8 * 'MPa' ;
	GV.Fvap.P = 8 * 'MPa' ;
	GV.Fvap.T = 568 * 'K' ;
	GV.Qpre = 1 * 'kW' ;
	GV.Fout_sa.P = 8 * 'MPa' ;
	GV.Fout_sa.T = 753 * 'K' ;
	GV.Feco.T = 563 * 'K' ;
	TAP.EF_T = 0.85 ;
	TAP.Fout.P = 1.9 * 'MPa' ;
	TAP.y = 0.3 ;
	TMP1.EF_T = 0.8 ;
	TMP1.Fout.P = 0.6 * 'MPa' ;
	TMP1.y = 0.05 ;
	TMP2.EF_T = 0.8 ;
	TMP2.Fout.P = 0.25 * 'MPa' ;
	TMP2.y = 0.05 ;
	TBP.EF_T = 0.7 ;
	TBP.Fout.P = 0.007 * 'MPa' ;
	COND.G_S = 2 * 'K' ;
	ATEMP.Fout.P = 1.9 * 'MPa' ;
	TD1.EF_T = 0.55 ;
	TD1.Fout.P = 0.25 * 'MPa' ;
	T4P.Fout.P = 0.25 * 'MPa' ;
	T4P.Fperda.F = 2.6 * 'kg/s' ;
	DEAERA.Fout.P = 0.16 * 'MPa' ;
	B4.EF_B = 0.95 ;
	B4.Fout.P = 1.9 * 'MPa' ;
	T2P.Fout.P = 0.6 * 'MPa' ;
	B3.EF_B = 0.95 ;
	B3.Fout.P = 0.6 * 'MPa' ;
	B5.EF_B = 0.95 ;
	B5.Fout.P = 8 * 'MPa' ;
	B2.EF_B = 0.95 ;
	B2.Fout.P = 0.25 * 'MPa' ;
	CONDP.G_S = 2 * 'K' ;
	B1.EF_B = 0.95 ;
	B1.Fout.P = 0.16 * 'MPa' ;
	CONDP6.G_S = 2 * 'K' ;
	FLASH.FoutL.P = 0.25 * 'MPa' ;
	SP4.y(1) = 0.001 ;
	SP4.y(2) = 0.001 ;
	SP4.y(3) = 0.001 ;
	Torre_1.F = 7000 * 'kg/s' ;
	Torre_1.Th = 308 * 'K' ;
	Torre_1.Tar_c = 293 * 'K' ;
	Torre_1.DTar = 10 * 'K' ;
	Fonte2_1.Fout.F = 2.6 * 'kg/s' ;
	Fonte2_1.Fout.P = 0.1013 * 'MPa' ;
	Fonte2_1.Fout.T = 293 * 'K' ;
	Splitter_1.y = 0.015 ;

	INITIAL

	GUESS

	OPTIONS
	Dynamic = false;
	Integration = "original";
	NLASolver(
		File = "sundials",
		RelativeAccuracy = 1e-5,
		AbsoluteAccuracy = 1e-6,
		MaxIterations = 100
	);
	DAESolver(
		File = "dassl",
		RelativeAccuracy = 1e-3,
		AbsoluteAccuracy = 1e-6,
		EventAccuracy = 1e-2
	);

#*
This tab can be used by users to insert small code
structures in the present simulation.

For example: Steps, Ramps, Sinusoidals, ... *#

PARAMETERS
PP2 as Plugin(Type="water");

#EQUATIONS

#Enter your code here.

#Example:
#*if time > 1000 * 's' then
	source_1.F = 100 * 'kmol/h';
else
	source_1.F = 200 * 'kmol/h';
end*#
end