source: trunk/tutorial/ThreeTank3.mso @ 988

#*---------------------------------------------------------------------
* This file is property of the author and cannot be used, copyed
* or modified without permission.
*
* Copyright (C) 2002-2004  the author
*----------------------------------------------------------------------
* Author: Rafael de Pelegrini Soares
* $Id: ThreeTank3.mso 843 2009-09-03 18:36:47Z bicca $
*----------------------------------------------------------------------
*
* Arquivo de exemplo do EMSO
*
*-------------------------------------------------------------------*#

Model tank_simplified

PARAMETERS
        pi  as Real (Brief="pi number",Default=3.141592);
        k       as Real (Brief="Valve Constant", Unit = 'm^2.5/h', Default=4);
        Dh  as Real     (Brief="Hydraulic diameter", Unit = 'm', Default=2);
        A       as Real (Brief="Tank area", Unit = 'm^2');

SET
        pi  = 3.141592;
        A   = (pi*Dh^2)/4;

VARIABLES
        V               as Real (Brief="Tank Volume", Lower = 0, Unit='m^3');
        Level   as Real (Brief="Tank level", Lower = 0, Unit='m');
in      Fin     as Real (Brief="Input flow", Lower = 0, Unit='m^3/h');
out     Fout    as Real (Brief="Output flow",Lower = 0, Unit='m^3/h');

EQUATIONS

"Mass balance"
        diff(V) = Fin - Fout;

"Valve equation"
        Fout = k*sqrt(Level);

"Tank volume"
        V = A*Level;
end


FlowSheet ThreeTank

DEVICES
        Feed    as Real (Brief="Feed Stream", Unit='m^3/h');
        
        Tank1   as tank_simplified (Brief="First Tank");
        Tank2   as tank_simplified (Brief="Second Tank");
        Tank3   as tank_simplified (Brief="Third Tank");

CONNECTIONS
        Feed        to  Tank1.Fin;
        Tank1.Fout      to      Tank2.Fin;
        Tank2.Fout      to      Tank3.Fin;

SPECIFY
        Feed = 10 * 'm^3/h';
        
INITIAL
        Tank1.Level = 1 * 'm';
        Tank2.Level = 2 * 'm';
        Tank3.Level = 1 * 'm';
        
SET
        Tank2.k = 8 * 'm^2.5/h';
        Tank2.Dh = 4 * 'm';
        
OPTIONS
        TimeStep = 0.1;
        TimeEnd = 10;
        TimeUnit = 'h';

end