source: trunk/sample/estimation/Bio.mso @ 690

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

Updating samples.
File size: 3.6 KB
Line 
1
2#*
3    Considerações:
4        1- Biomassa X é composta por 2 componentes
5                O componente com atividade catalitica consiste em acido nucleicos e proteinas -> R
6                componente inerte P.
7    2- Nitrogenio, N é o substrato limitante q afeta a cinética de crescimento.
8  -----------------------------  *#
9Model Biop_NE_mod5t
10
11        PARAMETERS
12
13   
14        Kss     as Real(Brief="Constante de saturação do modelo de Monod para sacarose",Unit='kg/m^3'); 
15        Ksn as Real(Brief="Constante de saturação do modelo de Monod para nitrogênio",Unit='kg/m^3');
16        mim     as Real(Brief="taxa de crescimento máximo específico",Unit='1/s');
17        alfa    as Real; 
18        gama    as Real(Unit='1/s');
19        K1              as Real(Unit='kg/kg');#(g PHB / g biomassa) -> constante de formação de produto associado ao crescimento de biomassa
20        K2              as Real(Unit='1/s');#(g PHB / g biomassa*h) -> constante de formação de produto não associado ao crescimento de biomassa
21        Yn              as Real(Brief="1/rendimento de biomassa baseado no nitrogenio",Unit='kg/kg');#(g biomassa/g nitrogenio  )
22        kd              as Real(Unit='1/s');
23        #kc             as Real(Unit='1/s');
24 
25        VARIABLES
26        X       as Real(Unit='kg/m^3');
27        R       as Real(Unit='kg/m^3');
28        P       as Real(Unit='kg/m^3');
29        S       as Real(Unit='kg/m^3');
30        N       as Real(Unit='kg/m^3');
31        PHB as Real;
32        mi  as Real(Unit='1/s');
33       
34        EQUATIONS
35
36        "Produção de Biomassa total"
37                X = R + P;
38
39        "Produção de Biomassa residual"
40       
41        diff(R) = mi * R - kd * R;
42       
43        "Produção de Biopolímero"
44                diff(P) = (K1*mi+ K2)*R;
45       
46        "Consumo de Nitrogenio"
47                diff(N) = -(mi/Yn)*R;
48       
49        if (S > 0) then
50        "Consumo de Açúcar"     
51                diff(S) = -(alfa*mi + gama)*R;
52
53        "Taxa de crescimento específica  Monod 2 substratos"
54                mi = mim * (S/(S + Kss))*(N/(N + Ksn));
55
56        else
57                diff(S) = 0 * 'kg/m^3/s';
58        "Taxa de crescimento específica  Monod 2 substratos"
59                mi = 0/'s';
60        end
61       
62        #alfa =(1/Yrs1 + K1/Yps1), onde Yrs1 é o rendimento de biomassa baseado na frutose (g biomassa/g frutose),
63        #gama = (ms1 + K2/Yps1) ,onde ms1 é o coeficiente de manutenção da biomassa baseado na frutose (g frutose/g biomassa*h)
64
65        "Percentual de acumulo"
66        PHB * X = P * 100;
67end
68
69FlowSheet Biop_NE_process5t as Biop_NE_mod5t
70        SET
71        Kss = 0.1928 * 'kg/m^3'; 
72        Ksn = 0.0447 * 'kg/m^3';
73        mim = 0.7979 * '1/s';
74        alfa =2.0606; 
75        gama = 0.0849 * '1/s';
76        K1 = 0.4053 * 'kg/kg';
77        K2 = -0.0079 * '1/s';
78        Yn = 10.6314;
79        kd = 0.0073 * '1/s';
80        #kc = 0.01 * '1/s';
81
82        #SPECIFY
83       
84       
85        INITIAL
86        R = 0.017 * 'kg/m^3';
87        P = 0.007 * 'kg/m^3';
88        S = 17.519 * 'kg/m^3';
89        N = 0.425 * 'kg/m^3'; # x 0,212 conversão de sulfato de amonio para N
90
91       
92        OPTIONS
93       
94        TimeStep = 0.1;
95        TimeEnd  = 24;
96        TimeUnit = 's';
97#       DAESolver(File="dassl");   
98end
99
100Estimation Biop_NE_Estt5 as Biop_NE_process5t
101        ESTIMATE
102        # PARAMETER  START    LOWER     UPPER     UNIT
103
104
105        Kss                     0.2009          0.004   7               'kg/m^3'; 
106        Ksn                     0.0446          0.005   1               'kg/m^3';
107        mim                     0.7979      0.1         0.8             '1/s';
108        alfa            2.0293          1               5                       ; 
109        gama            0.08502         0.05            5               '1/s';
110        K1                      0.4059          0.1             3               'kg/kg';
111        K2                      -0.00795        -1              3               '1/s';
112        Yn                      10.62           0.1             18                       ;
113        kd                      0.007           0.0005  1               '1/s';
114        #k3                     0.0131          0.0001  1               '1/s';
115        #K2                     0.00527         0.001   3               '1/s';
116       
117        #*Kss           1.5679          0.005           7               'kg/m^3'; 
118        Kss                     0.0487          0.004   7               'kg/m^3'; 
119        Ksn                     0.049           0.005   1               'kg/m^3';
120        mim                     0.799       0.1         0.8             '1/s';
121        alfa            2.0458          1               5                       ; 
122        gama            0.086           0.05            5               '1/s';
123        K1                      0.404           0.1             3               'kg/kg';
124        K2                      0.0049          0.001   3               '1/s';
125        Yn                      10.69           0.5             18                       ;
126        kd                      0.0078          0.001   1               '1/s';*#
127        #k3                     0.01274         0.0001  1               '1/s';
128       
129        EXPERIMENTS
130        # DATA FILE    WEIGTH
131        "Bio.dat"      1;
132       
133        OPTIONS
134        Statistics(
135                Fit=true,
136                Parameter=false,
137                Prediction=false
138        );
139       
140        NLPSolver(
141                MaxIterations = 1000,
142                File = "complex"
143                #File = "ipopt_emso"
144        );
145       
146end
Note: See TracBrowser for help on using the repository browser.