#*-------------------------------------------------------------------
* 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 86 2006-12-08 20:58:15Z paula $
*--------------------------------------------------------------------*#

using "stage_separators/batch_dist";

FlowSheet TesteBatch
	
	PARAMETERS
	PP	as CalcObject(Brief="Physical Properties",File="vrpp");
	NComp	as Integer;

	SET
	PP.Components = [ "isobutane", "benzene" ];
	PP.LiquidModel = "PR";
	PP.VapourModel = "PR";
	NComp = PP.NumberOfComponents;
	
	DEVICES
	batch 	as Diff_Dist;
	reflux 	as stream_therm;
	feed	as stream_therm;
	
	CONNECTIONS
	reflux to batch.InletL;
	feed to batch.Inlet;

	EQUATIONS
	if batch.Level > 1E-3 then
		batch.OutletV.F = 150 * "kmol/h";
		batch.Q = 3.7743e6 * "kJ/h";
	else
		batch.OutletV.F = 0 * "kmol/h";
		batch.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];
	reflux.v = 0;
	
	feed.F = 0 * "kmol/h";
	feed.T = 328 * "K";
	feed.P = 180 * "kPa";
	feed.z = [0.5, 0.5];
	feed.v = 0;
	
	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
	relativeAccuracy = 1e-3;
	time = [0:0.1:20]*"min";
end