#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*----------------------------------------------------------------------
* Author: Maurício Carvalho Maciel, Paula B. Staudt, Rafael P. Soares
* $Id: splitter.mso 910 2010-02-20 03:24:41Z arge $
*--------------------------------------------------------------------*#


using "streams";

Model splitter_n
	ATTRIBUTES
	Pallete 	= false;
	Icon 		= "icon/splitter_n"; 
	Brief 		= "Model of a splitter (NOT Handled by the GUI)";
	Info 		=
"== Assumptions ==
* thermodynamics equilibrium
* adiabatic
			
== Specify ==
* the inlet stream
* (Noutlet - 1) fraction of split of the outlet streams:

	frac(i) = (Mole Flow of the outlet stream i / 
				Mole Flow of the inlet stream)
				where i = 1, 2,...,Noutlet
";
	
	PARAMETERS
	NOutlet as Integer (Brief = "Number of Outlet Streams", Lower = 1);
	
	VARIABLES
in	Inlet   as stream (Brief = "Inlet stream", PosX=0, PosY=0.5001, Symbol="_{in}");
out Outlet(NOutlet)  as stream (Brief = "Outlet streams", PosX=1, PosY=0.5, Symbol="_{out}");
	frac(NOutlet) as fraction (Brief = "Distribution of Outlets", Default=0.5, Symbol="\phi");

	EQUATIONS
	
	sum(frac) = 1;
	
	for i in [1:NOutlet] do
		
		"Flow"
		Outlet(i).F = Inlet.F*frac(i);
		
		"Composition"
		Outlet(i).z = Inlet.z;
	
		"Pressure"
		Outlet(i).P = Inlet.P;
	
		"Enthalpy"
		Outlet(i).h = Inlet.h;
	
		"Temperature"	
		Outlet(i).T = Inlet.T;
	 
		"Vapourisation Fraction"
		Outlet(i).v = Inlet.v;
	end

end

Model splitter2
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/splitter_column"; 
	Brief 		= "Splitter with 2 outlet streams";
	Info 		=
"== Assumptions ==
*Thermodynamics equilibrium
*Adiabatic
			
== Specify ==
* The inlet stream
* One FlowRatios of split of the outlet streams:

	FlowRatios(i) = (Mole Flow of the outlet stream i / 
				Mole Flow of the inlet stream)
				where i = 1, 2
";

VARIABLES
in  Inlet   as stream (Brief = "Inlet stream", PosX=0.5, PosY=0, Symbol="_{in}");
out Outlet1 as stream (Brief = "Outlet stream 1", PosX=0.25, PosY=1, Symbol="_{out1}");
out Outlet2 as stream (Brief = "Outlet stream 2", PosX=0.75, PosY=1, Symbol="_{out2}");
	FlowRatios(2) 	as fraction 	(Brief = "Distribution of Outlets", Default=0.33, Symbol="\phi");

EQUATIONS

"Normalize Flow Ratios"	
	sum(FlowRatios) = 1;
	
"Flow"
	Outlet1.F = Inlet.F * FlowRatios(1);
	Outlet1.F + Outlet2.F = Inlet.F;

"Composition"
	Outlet1.z = Inlet.z;
	Outlet2.z = Inlet.z;

"Pressure"
	Outlet1.P = Inlet.P;
	Outlet2.P = Inlet.P;

"Enthalpy"
	Outlet1.h = Inlet.h;
	Outlet2.h = Inlet.h;

"Temperature"
	Outlet1.T = Inlet.T;
	Outlet2.T = Inlet.T;

"Vapourisation Fraction"
	Outlet1.v = Inlet.v;
	Outlet2.v = Inlet.v;

end

Model splitter3
	ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/splitter"; 
	Brief 		= "Model of a splitter with 3 outlet streams";
	Info 		=
"== Assumptions ==
* Thermodynamics equilibrium
* Adiabatic
			
== Specify ==
*The inlet stream
*Two FlowRatios of split of the outlet streams:

	FlowRatios(i) = (Mole Flow of the outlet stream i / 
				Mole Flow of the inlet stream)
				where i = 1, 2, 3
";
	
VARIABLES

in	Inlet   			as stream 	(Brief = "Inlet stream", PosX=0, PosY=0.5001, Symbol="_{in}");
out  Outlet1 	as stream 	(Brief = "Outlet stream 1", PosX=1, PosY=0.25, Symbol="_{Out1}");
out  Outlet2 	as stream 	(Brief = "Outlet stream 2", PosX=1, PosY=0.5059, Symbol="_{Out2}");
out  Outlet3 	as stream 	(Brief = "Outlet stream 3", PosX=1, PosY=0.75, Symbol="_{Out3}");

	FlowRatios(3) 	as fraction 	(Brief = "Distribution of Outlets", Default=0.33, Symbol="\phi");

EQUATIONS
	
"Normalize Flow Ratios"	
	sum(FlowRatios) = 1;
	
"Outlet1 Flow"
	Outlet1.F = Inlet.F*FlowRatios(1);
	
"Outlet2 Flow"
	Outlet2.F = Inlet.F*FlowRatios(2);
	
"Outlet3 Flow"
	Outlet3.F = Inlet.F*FlowRatios(3);
	
"Outlet1 Composition"
	Outlet1.z = Inlet.z;

"Outlet2 Composition"
	Outlet2.z = Inlet.z;

"Outlet3 Composition"
	Outlet3.z = Inlet.z;

"Outlet1 Pressure"
	Outlet1.P = Inlet.P;

"Outlet2 Pressure"
	Outlet2.P = Inlet.P;

"Outlet3 Pressure"
	Outlet3.P = Inlet.P;

"Outlet1 Enthalpy"
	Outlet1.h = Inlet.h;

"Outlet2 Enthalpy"
	Outlet2.h = Inlet.h;

"Outlet3 Enthalpy"
	Outlet3.h = Inlet.h;
	
"Outlet1 Temperature"	
	Outlet1.T = Inlet.T;

"Outlet2 Temperature"	
	Outlet2.T = Inlet.T;

"Outlet3 Temperature"	
	Outlet3.T = Inlet.T;

"Outlet1 Vapourisation Fraction"
	Outlet1.v = Inlet.v;

"Outlet2 Vapourisation Fraction"
	Outlet2.v = Inlet.v;
	
"Outlet3 Vapourisation Fraction"
	Outlet3.v = Inlet.v;

end

Model splitter_column

ATTRIBUTES
	Pallete 	= true;
	Icon 		= "icon/splitter_column"; 
	Brief 		= "Splitter with 2 outlet streams to be used with column section model";
	Info 		=
"== Assumptions ==
*Thermodynamics equilibrium
*Adiabatic

== Specify ==
* The inlet stream
* One Split Flow of the outlet stream:
**SplitFlow_RefluxStream: split of the outlet Reflux stream
**SplitFlow_DistillateStream: split of the outlet Distillate stream

OR

* The inlet stream and the RefluxRatio
";

VARIABLES
	in  Inlet   	as stream (Brief = "Inlet stream", PosX=0.5, PosY=0, Symbol="_{in}");
	out Reflux 		as stream (Brief = "Reflux stream", PosX=0.25, PosY=1, Symbol="_{reflux}");
	out Distillate 	as stream (Brief = "Distillate stream", PosX=0.75, PosY=1, Symbol="_{distillate}");

	SplitFlow_RefluxStream 	as fraction 	(Brief = "split of the outlet Reflux stream", Default=0.33);
	SplitFlow_DistillateStream 	as fraction 	(Brief = "split of the outlet Distillate stream", Default=0.33);
	RefluxRatio		as positive 	(Brief = "Reflux Ratio", Default=10, Lower = 0.05);

EQUATIONS

"Reflux Ratio"	
	#Reflux.F*RefluxRatio = SplitFlow_RefluxStream*(Reflux.F*RefluxRatio + Reflux.F);
	SplitFlow_RefluxStream*(1 + RefluxRatio) = RefluxRatio;
	
"Normalize Flow Ratios"	
	SplitFlow_RefluxStream + SplitFlow_DistillateStream = 1;
	
"Flow"
	Reflux.F = Inlet.F * SplitFlow_RefluxStream;
	Inlet.F = Reflux.F + Distillate.F;

"Composition"
	Reflux.z = Inlet.z;
	Distillate.z = Inlet.z;

"Pressure"
	Reflux.P = Inlet.P;
	Distillate.P = Inlet.P;

"Enthalpy"
	Reflux.h = Inlet.h;
	Distillate.h = Inlet.h;

"Temperature"
	Reflux.T = Inlet.T;
	Distillate.T = Inlet.T;

"Vapourisation Fraction"
	Reflux.v = Inlet.v;
	Distillate.v = Inlet.v;

end