source: mso/sample/miscellaneous/sample_bio.mso @ 51

#*---------------------------------------------------------------------
* 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: sample_bio.mso 32 2006-09-23 00:16:00Z arge $
*----------------------------------------------------------------------
*
* EMSO sample file.
*
* FlowSheet with a bio-reactor based on Monod model.
*
*-------------------------------------------------------------------*#

using "types";

FlowSheet bio
        PARAMETERS
        miMax   as frequency(Brief="Monod parameter", Unit="1/h", Default=0.53);
        D       as frequency(Brief="Dilution rate F/V", Unit="1/h");
        Km      as conc_mass(Brief="Monod parameter", Default=0.12);
        K1      as inv_conc_mass(Brief="Monod parameter", Default=0.4545);
        Y       as coefficient(Brief="Yield biomass/substrate", Default=0.4);
        x2f     as conc_mass(Brief="Substrate feed concentration", Default=4);
 
        VARIABLES
        biomass   as conc_mass (Brief="Biomass concentration");
        substrate as conc_mass (Brief="Substrate concentration");
        mi        as frequency (Brief="Specific growth rate", Unit="1/h");
        
        EQUATIONS
        "Specific growth rate"
        mi = miMax*substrate/(Km+substrate+(K1*substrate^2));
        
        "Biomass production"
        diff(biomass)/biomass = mi-D;
        
        "Substrate consumption"
        diff(substrate) = D*(x2f-substrate)-(biomass*mi/Y);
        
        SET
        D      = 0.3  *  "1/h";
        
        INITIAL
        biomass = 1  *  "kg/m^3";
        substrate = 0.5  *  "kg/m^3";
 
        OPTIONS
        time = [0:0.5:10]*"h";
        LASolver = "dense"; # dense or sparse
end