#*-------------------------------------------------------------------
* 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 Diff_Dist
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Maurício Carvalho Maciel
* $Id: sample_batch_dist.mso 318 2007-07-19 04:23:31Z arge $
*--------------------------------------------------------------------*#

using "stage_separators/batch_dist";

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

	SET
	NComp = PP.NumberOfComponents;
	
	DEVICES
	batch 	as Diff_Dist;
	reflux 	as liquid_stream;
	feed	as source;
	Q		as energy_source;
	
	CONNECTIONS
	reflux to batch.InletL;
	feed.Outlet to batch.Inlet;
	Q.OutletQ to batch.InletQ;

	EQUATIONS
	if batch.Level > 1E-3 then
		batch.OutletV.F = 150 * 'kmol/h';
		Q.OutletQ.Q = 3.7743e6 * 'kJ/h';
	else
		batch.OutletV.F = 0 * 'kmol/h';
		Q.OutletQ.Q = 0 * 'kJ/h';
	end
	
	SPECIFY
	reflux.F = 0 * 'kmol/h';
	reflux.T = 328 * 'K';
	reflux.P = 180 * 'kPa';
	reflux.z = [0.5, 0.5];
	
	feed.Outlet.F = 0 * 'kmol/h';
	feed.Outlet.T = 328 * 'K';
	feed.Outlet.P = 180 * 'kPa';
	feed.Outlet.z = [0.5, 0.5];
	
	SET
	batch.V = 3 * 'm^3';
	batch.Across = 1 * 'm^2';
	
	INITIAL
	batch.T = 298 *'K';
	batch.Level = 2.5 * 'm';
	batch.x(2) = 0.2;
	
	OPTIONS
	TimeEnd = 20;
	TimeStep = 0.1;
	TimeUnit = 'min';
end