#*-------------------------------------------------------------------
* Biorrefinaria Petrobras
*--------------------------------------------------------------------
* Nome do arquivo: heater_m.mso
* Projeto: Modelo integrado de producao de etanol 1G/2G
* Conteudo: trocador de calor 
*--------------------------------------------------------------------*#

#*-------------------------------------------------------------------
*
* Versao 2.2
* Data:    03/2016
* Autores:   Anderson R. A. Lino e Gabriel C. Fonseca
* 
*--------------------------------------------------------------------
* Descricao: modelo simplificado do trocador de calor que sera  
* empregado na biorrefinaria
*--------------------------------------------------------------------

*--------------------------------------------------------------------
*Hipoteses assumidas: 1 - operacao em estado estacionario
*                     2 - sem perdas de calor para o ambiente
*--------------------------------------------------------------------

*--------------------------------------------------------------------
*Notas: Foi feito o flowsheet teste para averiguar o modelo.
*As siglas R, L, U, D sao indicativos da posicao da corrente de saida
*sendo R=right, L=legth, U=up e D=down
*--------------------------------------------------------------------*#

using "energy_stream";
using "main_stream";

Model heatex_mL
	ATTRIBUTES
	Pallete = false;
	Brief = "Model of a Heat Exchanger";
	Info =
"== GENERAL ==
* Model of a simplified heat exchanger;
* The model performes mass and energy balances only.
	
== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.
";


#*-------------------------------------------------------------------
* Declaracao de variaveis
*--------------------------------------------------------------------*#
	
VARIABLES

in  Inlet		as main_stream 		(Brief="Inlet Stream", PosX=1.0, PosY=0.5, Protected = false, Symbol="_{in}");
out Outlet		as main_stream_PH  	(Brief="Outlet Stream", PosX=0.0, PosY=0.5, Protected = false, Symbol="_{out}");
	Q	   		as power			(Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
	
	Pdrop 		as press_delta		(Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
	lmtd		as temp_delta		(Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
	U	 		as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
	A 			as area				(Brief = "Exchange Surface Area", Lower = 1e-6);
	

#*-------------------------------------------------------------------
* Equacoes do modelo
*--------------------------------------------------------------------*#	

	EQUATIONS

	"Energy Balance"
		Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);

	"Molar Balance (Fluid Phase)"
		Inlet.Fluid.F  = Outlet.Fluid.F;
		
	"Molar Balance (Solid Phase)"
		Inlet.Solid.F = Outlet.Solid.F;

	"Fraction Constraint (Fluid Phase)"
		Outlet.Fluid.z = Inlet.Fluid.z;
		
	"Fraction Constraint (Solid Phase)"
		Outlet.Solid.z = Inlet.Solid.z;
		
	"Pressure Drop"
		Outlet.P = Inlet.P - Pdrop;
	
end

Model heater_mLD as heatex_mL
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterD";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";

	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{in}");

	EQUATIONS
	
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mLD as heatex_mL
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterD";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heater_mLU as heatex_mL
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterU";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mLU as heatex_mL
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterU";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heatex_mR

	ATTRIBUTES
	Pallete 	= false;
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
* Model of a simplified heat exchanger;
* The model performes mass and energy balances only.
	
== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.
";

#*-------------------------------------------------------------------
* Declaracao de variaveis
*--------------------------------------------------------------------*#
	
VARIABLES

in  Inlet		as main_stream 		(Brief="Inlet Stream", PosX=0.0, PosY=0.5, Protected = false, Symbol="_{in}");
out Outlet		as main_stream_PH  	(Brief="Outlet Stream", PosX=1.0, PosY=0.5, Protected = false, Symbol="_{out}");
	Q	   		as power			(Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
	
	Pdrop 		as press_delta		(Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
	lmtd		as temp_delta		(Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
	U	 		as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
	A 			as area				(Brief = "Exchange Surface Area", Lower = 1e-6);
	
#*-------------------------------------------------------------------
* Equacoes do modelo
*--------------------------------------------------------------------*#	

	EQUATIONS

	"Energy Balance"
		Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);

	"Molar Balance (Fluid Phase)"
		Inlet.Fluid.F  = Outlet.Fluid.F;
		
	"Molar Balance (Solid Phase)"
		Inlet.Solid.F = Outlet.Solid.F;

	"Fraction Constraint (Fluid Phase)"
		Outlet.Fluid.z	= Inlet.Fluid.z;
		
	"Fraction Constraint (Solid Phase)"
		Outlet.Solid.z = Inlet.Solid.z;
		
	"Pressure Drop"
		Outlet.P = Inlet.P - Pdrop;
end

Model heater_mRD as heatex_mR
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterD";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mRD as heatex_mR
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterD";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heater_mRU as heatex_mR
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterU";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mRU as heatex_mR
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterU";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heatex_mU

	ATTRIBUTES
	Pallete 	= false;
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
* Model of a simplified heat exchanger;
* The model performes mass and energy balances only.
	
== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.
";

#*-------------------------------------------------------------------
* Declaracao de variaveis
*--------------------------------------------------------------------*#
	
VARIABLES

in  Inlet		as main_stream 		(Brief="Inlet Stream", PosX=0.5, PosY=1.0, Protected = false, Symbol="_{in}");
out Outlet		as main_stream_PH  	(Brief="Outlet Stream", PosX=0.5, PosY=0.0, Protected = false, Symbol="_{out}");
	Q	   		as power			(Brief="Heat Duty", Default=7000, Lower=-1e10, Upper=1e10);
	
	Pdrop 		as press_delta		(Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
	lmtd		as temp_delta		(Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
	U	 		as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
	A 			as area				(Brief = "Exchange Surface Area", Lower = 1e-6);
	
#*-------------------------------------------------------------------
* Equacoes do modelo
*--------------------------------------------------------------------*#	

	EQUATIONS

	"Energy Balance"
		Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);

	"Molar Balance (Fluid Phase)"
		Inlet.Fluid.F  = Outlet.Fluid.F;
		
	"Molar Balance (Solid Phase)"
		Inlet.Solid.F = Outlet.Solid.F;

	"Fraction Constraint (Fluid Phase)"
		Outlet.Fluid.z = Inlet.Fluid.z;
		
	"Fraction Constraint (Solid Phase)"
		Outlet.Solid.z = Inlet.Solid.z;
		
	"Pressure Drop"
		Outlet.P = Inlet.P - Pdrop;	
end

Model heater_mUR as heatex_mU
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterR";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";

	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mUR as heatex_mU
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterR";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";

	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heater_mUL as heatex_mU
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterL";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mUL as heatex_mU
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterL";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heatex_mD

	ATTRIBUTES
	Pallete 	= false;
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
* Model of a simplified heat exchanger;
* The model performes mass and energy balances only.
	
== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.
";

#*-------------------------------------------------------------------
* Declaracao de variaveis
*--------------------------------------------------------------------*#
	
VARIABLES

in  Inlet		as main_stream 		(Brief="Inlet Stream", PosX=0.5, PosY=0.0, Protected = false, Symbol="_{in}");
out Outlet		as main_stream_PH  	(Brief="Outlet Stream", PosX=0.5, PosY=1.0, Protected = false, Symbol="_{out}");
	Q	   		as power			(Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
	
	Pdrop 		as press_delta		(Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
	lmtd		as temp_delta		(Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
	U	 		as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
	A 			as area				(Brief = "Exchange Surface Area", Lower = 1e-6);
	
#*-------------------------------------------------------------------
* Equacoes do modelo
*--------------------------------------------------------------------*#	

	EQUATIONS

	"Energy Balance"
		Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);

	"Molar Balance (Fluid Phase)"
		Inlet.Fluid.F  = Outlet.Fluid.F;
		
	"Molar Balance (Solid Phase)"
		Inlet.Solid.F = Outlet.Solid.F;

	"Fraction Constraint (Fluid Phase)"
		Outlet.Fluid.z = Inlet.Fluid.z;
		
	"Fraction Constraint (Solid Phase)"
		Outlet.Solid.z = Inlet.Solid.z;
		
	"Pressure Drop"
		Outlet.P = Inlet.P - Pdrop;	
end

Model heater_mDR as heatex_mD
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterR";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mDR as heatex_mD
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterL";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

Model heater_mDL as heatex_mD
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterR";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet stream temperature;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
in	Inlet_q	as heat_stream	(Brief = "Inlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{in}");

	EQUATIONS
	"Equate Heat Stream"
	Inlet_q.Q = Q;
	
	"General equation of heat exchange"
	Inlet_q.Q = U * A * lmtd;
	
end

Model cooler_mDL as heatex_mD
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/heaterR";
	Brief		= "Model of a Heat Exchanger";
	Info = 
"== GENERAL ==
*Model of a simplified heat exchanger;
*The model performes mass and energy balances only.

== ASSUMPTIONS ==
* Steady-State operation;
* No heat loss to the surroundings.

== SPECIFY ==
* The inlet streams:
  flow rate
  temperature
  pressure
  stream composition;
* Outlet vapour fraction;
* Overall Heat Transfer Coefficient;
* Exchange Surface Area;
* Pressure drop.

== SET ==
* Phase of the fluid entering the heater;
* Number of stream components(Ncomp/NcompS);
* Basic composition (mass or molar).
";
	
	VARIABLES
out	Outlet_q	as heat_stream	(Brief = "Outlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{out}");

	EQUATIONS
	"Equate Heat Stream"
	Outlet_q.Q = -Q;
	
	"General equation of heat exchange"
	Outlet_q.Q = U * A * lmtd;
	
end

FlowSheet teste_heater_m
	
#*-------------------------------------------------------------------
* Declaracao de dispositivos (ou blocos contendo o modelo)
*--------------------------------------------------------------------*#	
	
	DEVICES
	SHot as main_sourceR;
	H as heater_mLD;
	SQ as heat_sourceR;
	
#*-------------------------------------------------------------------
* Especifica as conexoes entre os modelos
*--------------------------------------------------------------------*#	
	
	CONNECTIONS
	SHot.Outlet to H.Inlet;
	SQ.Outlet_q to H.Inlet_q;
	
#*-------------------------------------------------------------------
* Especifica variaveis definidas no modelo
*--------------------------------------------------------------------*#
	
	SPECIFY

	SHot.Fluid.Fw = 52200 * 'kg/h';
	SHot.Solid.Fw = 48800 * 'kg/h';
	SHot.T = 340 * 'K';
	SHot.P = 1 * 'atm';
	SHot.CompositionOfSolid = [0.5110, 0.0657, 0.2190, 0.0001, 0, 0, 0, 0.0004, 0];
	SHot.CompositionOfFluid(1:11) = [0.982744, 0.000223825, 0.000463338, 0.00172757, 0, 0.00068594, 0.000144428, 0, 0, 0.000383508, 0.0130415];
	SHot.CompositionOfFluid(12:NComp) = 0;
	H.Pdrop = 0 * 'atm';
	H.Outlet.T = 350.0 * 'K';
	H.U = 0.69445 * 'kW/m^2/K';
	H.lmtd = 10 * 'K';



#*-------------------------------------------------------------------
#Parametros
*--------------------------------------------------------------------*#
	 
	PARAMETERS
	PP as Plugin	(Brief = "External Physical Properties", 
		Type="PP",
		Project = "../Flowsheets/v2_2/Fluid_v2_2.vrtherm" 
	);
	PPS as Plugin	(Brief = "External Physical Properties", 
		Type="PP",
		Project = "../Flowsheets/v2_2/Solid_v2_2.vrtherm"
	);
	
	NComp 	as Integer 	(Brief = "Number of chemical components in the fluid phase");
	NCompS 	as Integer 	(Brief = "Number of chemical components in the solid phase");
	
#*-------------------------------------------------------------------
* Define o valor dos parametros declarados no modelo
*--------------------------------------------------------------------*#
	
	SET
	NComp = PP.NumberOfComponents();
	NCompS = PPS.NumberOfComponents();
	SHot.ValidPhases = "Liquid-Only";
	SHot.CompositionBasis = "Molar";
#*-------------------------------------------------------------------
* Opcoes de Solver
*--------------------------------------------------------------------*#
	
	OPTIONS
	Dynamic = false;
	#GuessFile = "teste3.rlt";
	
end

FlowSheet teste_cooler_m
	
#*-------------------------------------------------------------------
* Declaracao de dispositivos (ou blocos contendo o modelo)
*--------------------------------------------------------------------*#	
	
	DEVICES
	Scold as main_sourceR;
	H as cooler_mLD;
	
	
#*-------------------------------------------------------------------
* Especifica as conexoes entre os modelos
*--------------------------------------------------------------------*#	
	
	CONNECTIONS
	Scold.Outlet to H.Inlet;
	
#*-------------------------------------------------------------------
* Especifica variaveis definidas no modelo
*--------------------------------------------------------------------*#
	
	SPECIFY

	Scold.Fluid.Fw = 52200 * 'kg/h';
	Scold.Solid.Fw = 4800 * 'kg/h';
	Scold.T = 345 * 'K';
	Scold.P = 1 * 'atm';
	Scold.CompositionOfSolid = [0.5110, 0.0657, 0.2190, 0.0001, 0, 0, 0, 0.0004, 0];
	Scold.CompositionOfFluid(1:11) = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0];
	Scold.CompositionOfFluid(12:NComp) = 0;
	H.Pdrop = 0 * 'atm';
	H.Outlet.T = 335 * 'K';
	H.U = 0.69445 * 'kW/m^2/K';
	H.lmtd = 10 * 'K';



#*-------------------------------------------------------------------
#Parametros
*--------------------------------------------------------------------*#
	 
	PARAMETERS
	PP as Plugin	(Brief = "External Physical Properties", 
		Type="PP",
		Project = "../Flowsheets/v2_2/Fluid_v2_2.vrtherm" 
	);
	PPS as Plugin	(Brief = "External Physical Properties", 
		Type="PP",
		Project = "../Flowsheets/v2_2/Solid_v2_2.vrtherm"
	);
	
	NComp 	as Integer 	(Brief = "Number of chemical components in the fluid phase");
	NCompS 	as Integer 	(Brief = "Number of chemical components in the solid phase");
	
#*-------------------------------------------------------------------
* Define o valor dos parametros declarados no modelo
*--------------------------------------------------------------------*#
	
	
	SET
	NComp = PP.NumberOfComponents();
	NCompS = PPS.NumberOfComponents();
	Scold.ValidPhases = "Liquid-Only";
	Scold.CompositionBasis = "Molar";
	
#*-------------------------------------------------------------------
* Opcoes de Solver
*--------------------------------------------------------------------*#
	
	OPTIONS
	Dynamic = false;

	
end