source: trunk/sample/stage_separators/sample_batch_dist.mso @ 623

#*-------------------------------------------------------------------
* 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