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

Last change on this file since 51 was 32, checked in by Argimiro Resende Secchi, 16 years ago

clean up.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 1.6 KB
Line 
1#*---------------------------------------------------------------------
2* This file is property of the author and cannot be used, copyed
3* or modified without permission.
4*
5* Copyright (C) 2002-2004  the author
6*----------------------------------------------------------------------
7* Author: Rafael de Pelegrini Soares
8* $Id: sample_bio.mso 32 2006-09-23 00:16:00Z arge $
9*----------------------------------------------------------------------
10*
11* EMSO sample file.
12*
13* FlowSheet with a bio-reactor based on Monod model.
14*
15*-------------------------------------------------------------------*#
16
17using "types";
18
19FlowSheet bio
20        PARAMETERS
21        miMax   as frequency(Brief="Monod parameter", Unit="1/h", Default=0.53);
22        D       as frequency(Brief="Dilution rate F/V", Unit="1/h");
23        Km      as conc_mass(Brief="Monod parameter", Default=0.12);
24        K1      as inv_conc_mass(Brief="Monod parameter", Default=0.4545);
25        Y       as coefficient(Brief="Yield biomass/substrate", Default=0.4);
26        x2f     as conc_mass(Brief="Substrate feed concentration", Default=4);
27 
28        VARIABLES
29        biomass   as conc_mass (Brief="Biomass concentration");
30        substrate as conc_mass (Brief="Substrate concentration");
31        mi        as frequency (Brief="Specific growth rate", Unit="1/h");
32       
33        EQUATIONS
34        "Specific growth rate"
35        mi = miMax*substrate/(Km+substrate+(K1*substrate^2));
36       
37        "Biomass production"
38        diff(biomass)/biomass = mi-D;
39       
40        "Substrate consumption"
41        diff(substrate) = D*(x2f-substrate)-(biomass*mi/Y);
42       
43        SET
44        D      = 0.3  *  "1/h";
45       
46        INITIAL
47        biomass = 1  *  "kg/m^3";
48        substrate = 0.5  *  "kg/m^3";
49 
50        OPTIONS
51        time = [0:0.5:10]*"h";
52        LASolver = "dense"; # dense or sparse
53end
Note: See TracBrowser for help on using the repository browser.