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 | |
---|
12 | using "types"; |
---|
13 | |
---|
14 | FlowSheet 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; |
---|
77 | end |
---|