#*-------------------------------------------------------------------
* 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.
*
*--------------------------------------------------------------------
* Sample file for model streams
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Rafael de Pelegrini Soares
* $Id: sample_streams.mso 218 2007-03-16 02:55:04Z arge $
*--------------------------------------------------------------------*#

using "streams";

FlowSheet SampleStreams
	PARAMETERS
	NComp as Integer(Default=5);

	DEVICES
	s1 as stream;

	SPECIFY
	s1.F = 10  * 'kmol/h';
	s1.T = 300 * 'K';
	s1.P = 1   * 'atm';
	s1.z = 1/NComp;
	s1.h = 5   * 'J/kmol';
	s1.v = 0.0;

	OPTIONS
	Dynamic = false;
end

FlowSheet SampleStreamsTherm
	PARAMETERS
	PP   as Plugin(Brief="Physical Properties", Type="PP",
		Components = [ "isobutane", "benzene", "ethanol" ],
		LiquidModel = "IdealLiquid",
		VapourModel = "Ideal");
	NComp as Integer;

	SET
	NComp = PP.NumberOfComponents;

	DEVICES
	sl as source;
	sv as source;
	
	LiqOut  as sink;
	VapOut as sink;
	
	CONNECTIONS
	
	sl.Outlet to LiqOut.Inlet;
	sv.Outlet to VapOut.Inlet;

	
	SPECIFY
	sl.Outlet.F = 10  * 'kmol/h';
	sl.Outlet.T = 300 * 'K';
	sl.Outlet.P = 100 * 'kPa';
	sl.Outlet.z =  [0.24596, 0.370278, 0.383762];
#	sl.Outlet.v = 0.0; #liquid stream

	
	sv.Outlet.F = 10  * 'kmol/h';
	#sv.Outlet.h = 20248 * 'J/mol';
	sv.Outlet.T = 300 * 'K';
	sv.Outlet.P = 100 * 'kPa';
	#sv.Outlet.z = [0.824212, 0.12577, 0.0500176];
#	sv.Outlet.v = 1.0; #vapor stream

	EQUATIONS
	"Liquid-Vapour equilibrium"
	PP.VapourFugacityCoefficient(sv.Outlet.T, sv.Outlet.P, sv.Outlet.z) * sv.Outlet.z =
		PP.LiquidFugacityCoefficient(sl.Outlet.T, sl.Outlet.P, sl.Outlet.z) * sl.Outlet.z;
	
	OPTIONS
	Dynamic= true;
end