| 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 |
|---|
