#*-------------------------------------------------------------------
* 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 column model
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO(www.vrtech.com.br) to run
* SectionColumn_Test and needs VRTherm full to run the distillation
* column flowsheet.
*
*----------------------------------------------------------------------
* Author: Paula B. Staudt
* $Id: sample_column.mso 580 2008-07-26 19:34:58Z bicca $
*--------------------------------------------------------------------*#

using "stage_separators/columnTeste";

FlowSheet SectionColumn_Test_with2tray

PARAMETERS
	PP	as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = [ "isobutane", "benzene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp	as Integer;

DEVICES
	sec 	as Section_ColumnTeste;
	feed 	as liquid_stream;
	reb 	as vapour_stream;
	cond as liquid_stream;
	
SET
	NComp = PP.NumberOfComponents;
	sec.FeedTray=1;
	sec.NumberOfTrays = 2;
	#COLUMN
	sec.V = 4 * 'ft^3';
	sec.Ah = 0.394 * 'ft^2';
	sec.lw = 20.94 * 'in';
	sec.hw = 0.125 * 'ft';
	sec.Q = 0 * 'kW';
	sec.beta = 0.6;
	sec.alfa = 4;
	sec.Ap = 3.94 * 'ft^2';
	
CONNECTIONS
	feed to sec.Inlet;
	reb to sec.Vin;
	cond to sec.Lin;
	
SPECIFY
	feed.F = 113.4 * 'kmol/h';
	feed.T = 291 * 'K';
	feed.P = 168.3 * 'kPa';
	feed.z = [0.5, 0.5];

	cond.F = 68 * 'kmol/h';
	cond.P = 150 * 'kPa';
	cond.T = 281.75 * 'K';
	cond.z = [0.6664, 0.3336];

	reb.P = 185 * 'kPa';
	reb.T = 328.12 * 'K';
	reb.z = [0.001848, 0.9982];
	
	sec.trays.Emv = 1;
	sec.trays(1).OutletV.F = 150 * 'kmol/h';

INITIAL
	sec.trays.OutletL.T = 270 *'K';
	sec.trays.Level = 0.9 * sec.hw;
	sec.trays.OutletL.z(1) = 0.5;
	
	OPTIONS
	TimeStep = 10;
	TimeEnd = 1000;
	
end