source: trunk/sample/auto/ethanol_bif.mso @ 463

Last change on this file since 463 was 463, checked in by Argimiro Resende Secchi, 15 years ago

Adding new example using AUTO.

File size: 2.3 KB
Line 
1#*
2 Jobses' Model
3 Ref: Andres Mahecha-Botero, Parag Garhyan, S.S.E.H. Elnashaie
4          "Non-linear characteristics of a membrane fermentor
5                for ethanol production and their implications"
6          Nonlinear Analysis: RealWorld Applications 7 (2006) 432 – 457.
7*----------------------------------------------------------------------
8* Author: Fábio Diehl and Argimiro R. Secchi
9* $Id: ethanol_bif.mso 1 2008-02-19 12:17:36Z arge $
10*--------------------------------------------------------------------*#
11
12using "types";
13
14FlowSheet Ethanol_bif
15        PARAMETERS
16        Cs0             as conc_mass (Brief = "Substrate Initial Concentration");
17        D               as frequency (Brief = "Dilution Rate", DisplayUnit = '1/h');
18        Ks              as conc_mass (Brief = "Monod Constant");
19        Ysx             as positive  (Brief = "Yield factor based on substrate (kg/kg)");
20        Ypx             as positive  (Brief = "Yield factor based on product (kg/kg)");
21        ms              as frequency (Brief = "Maintenance factor based on substrate (kg/kg h)");
22        mp              as frequency (Brief = "Maintenance factor based on product (kg/kg h)");
23        k1              as frequency (Brief = "Empirical constant");
24        k2              as positive (Brief = "Empirical constant", Unit = 'm^3/kg/h');
25        k3              as positive (Brief = "Empirical constant", Unit = 'm^6/kg^2/h');
26       
27        VARIABLES
28        Cs              as conc_mass (Brief = "Substrate Concentration");
29        Ce              as conc_mass (Brief = "Key-Product Concentration");
30        Cx              as conc_mass (Brief = "Biomass Concentration");
31        Cp              as conc_mass (Brief = "Ethanol Concentration");
32
33        EQUATIONS
34        diff(Cs) = -Cs*Ce/(Ks+Cs)/(Ysx*'h') - ms*Cx + D*(Cs0-Cs);
35        diff(Ce) = (k1-k2*Cp+k3*Cp^2)*Cs*Ce/(Ks+Cs) - D*Ce;
36        diff(Cx) = Cs*Ce/(Ks+Cs)/'h' - D*Cx;
37        diff(Cp) = Cs*Ce/(Ks+Cs)/(Ypx*'h') + mp*Cx - D*Cp;
38
39        SET
40        Cs0 = 200 * 'kg/m^3';
41        D = 1/'h';
42        Ks = 0.5*'kg/m^3';
43        Ysx = 1/40.9;
44        Ypx = 1/19;
45        ms = 2.16/'h';
46        mp = 1.1/'h';
47        k1 = 16/'h';
48        k2 = 0.497*'m^3/kg/h';
49        k3 = 0.00383*'m^6/kg^2/h';
50
51        GUESS
52#*
53        Cs = 40 * 'kg/m^3';
54        Ce = 10 * 'kg/m^3';
55        Cx = 10 * 'kg/m^3';
56        Cp = 90 * 'kg/m^3';
57*#
58        Cs = 20 * 'kg/m^3';
59        Ce = 1 * 'kg/m^3';
60        Cx = 10 * 'kg/m^3';
61        Cp = 40 * 'kg/m^3';
62               
63        INITIAL
64        Cs = 40 * 'kg/m^3';
65        Ce = 10 * 'kg/m^3';
66        Cx = 10 * 'kg/m^3';
67#       Cp = 90 * 'kg/m^3';
68        Cp = 40 * 'kg/m^3';
69
70        OPTIONS
71        TimeStart = 0;
72        TimeStep = 0.5;
73        TimeEnd = 30;
74        TimeUnit = 'h';
75        DAESolver(File = "dassl");
76#       Dynamic = false;
77end
Note: See TracBrowser for help on using the repository browser.