source: branches/gui/sample/reactors/sample_gibbs.mso @ 832

#*---------------------------------------------------------------------
* 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 of a Gibbs reactor
*----------------------------------------------------------------------
*
* This sample file needs VRTherm (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
*
*       This is a sample based on a problem of the book "Introduction to
* Chemical Engineering Thermodynamics" by Smith, Van Ness and Abbott,
* 5th ed., 1996.
*
*----------------------------------------------------------------------
* Author: Rodolfo Rodrigues
* $Id$
*--------------------------------------------------------------------*#

using "reactors/gibbs";


#*---------------------------------------------------------------------
* Problem 15.13 - p.543 (Brazilian edition)
*
* System in gasosous fase contents CH4, H2O, CO, CO2, and H2
*
* Initial conditions
*   1bar, 1000K, 2mol of CH4 and 3mol of H2O
*
* Specify
*       Outlet.T
*       Outlet.P
*
* Return
*       Outlet.z
---------------------------------------------------------------------*#

FlowSheet sample_gibbs_vap
        PARAMETERS
        PP              as Plugin(Brief="External physical properties", Type="PP",
                        Components = ["methane","water","carbon monoxide","carbon dioxide","hydrogen"],
                        LiquidModel = "IdealLiquid",
                        VapourModel = "Ideal");
        NComp   as Integer;
        NElem   as Integer;
        
        DEVICES
        Fin     as source;
        R               as gibbs_vap; # gibbs reactor

        CONNECTIONS
        Fin.Outlet      to R.Inlet;
        
        SET
        NComp = PP.NumberOfComponents; # CH4, H2O, CO, CO2, H2
        NElem = 3; # C, O, H
        
        R.na(:,1) = [1.0, 0.0, 4.0]; # CH4
        R.na(:,2) = [0.0, 1.0, 2.0]; # H2O
        R.na(:,3) = [1.0, 1.0, 0.0]; # CO
        R.na(:,4) = [1.0, 2.0, 0.0]; # CO2
        R.na(:,5) = [0.0, 0.0, 2.0]; # H2

        SPECIFY
        Fin.Outlet.F = 5*'kmol/h';
        Fin.Outlet.z = [0.4, 0.6, 0.0, 0.0, 0.0];
        Fin.Outlet.P = 1*'atm';
        Fin.Outlet.T = 1000*'K';
        
        R.Outlet.F = 5*'kmol/h'*sqrt(R.Tank.Level/'m');
        R.Outlet.P = 1*'atm';
        R.Outlet.T = Fin.Outlet.T;

        R.Tank.Across = 5*'m^2';
        R.Tank.L = 1.5*'m';
        
        INITIAL
        R.Outletm.T = 1000*'K';
        R.M = [2, 3, 0, 0, 0]*'kmol';
        
        OPTIONS
        Dynamic = false;
        NLASolver(File="sundials",
                          RelativeAccuracy=1e-6,
                          AbsoluteAccuracy=1e-8);
end