Initialization file for all states and parameters related to the AS systems (reactors 1-5).
All parameters and specifications are based on BSM1 model. This file will be executed when running bsm2_cl.py, bsm2_ol.py or bsm2_olem.py.
QIN0
Flow rate of influent [m³ ⋅ d⁻¹].
QIN
Flow rate of influent [m³ ⋅ d⁻¹].
QINTR
Flow rate of internal recirculation [m³ ⋅ d⁻¹].
QR
Flow rate of sludge return [m³ ⋅ d⁻¹].
QW
Flow rate of waste sludge [m³ ⋅ d⁻¹].
S_I1
Soluble inert organic matter (reactor 1) [g(COD) ⋅ m⁻³].
S_S1
Readily biodegradable substrate (reactor 1) [g(COD) ⋅ m⁻³].
X_I1
Particulate inert organic matter (reactor 1) [g(COD) ⋅ m⁻³].
X_S1
Slowly biodegradable substrate (reactor 1) [g(COD) ⋅ m⁻³].
X_BH1
Active heterotrophic biomass (reactor 1) [g(COD) ⋅ m⁻³].
X_BA1
Active autotrophic biomass (reactor 1) [g(COD) ⋅ m⁻³].
X_P1
Particulate products arising from biomass decay (reactor 1) [g(COD) ⋅ m⁻³].
S_O1
Dissolved oxygen (reactor 1) [g(O₂) ⋅ m⁻³].
S_NO1
Nitrate and nitrite (reactor 1) [g(N) ⋅ m⁻³].
S_NH1
Ammonium plus ammonia nitrogen (reactor 1) [g(N) ⋅ m⁻³].
S_ND1
Soluble biodegradable organic nitrogen (reactor 1) [g(N) ⋅ m⁻³].
X_ND1
Particulate biodegradable organic nitrogen (reactor 1) [g(N) ⋅ m⁻³].
S_ALK1
Alkalinity (reactor 1) [mol(HCO₃⁻) ⋅ m⁻³].
TSS1
Total suspended solids (reactor 1) [g(TSS) ⋅ m⁻³].
Q1
Flow rate (reactor 1) [m³ ⋅ d⁻¹].
T1
Temperature (reactor 1) [°C].
S_D1_1
Dummy state 1 (reactor 1) [-].
S_D2_1
Dummy state 2 (reactor 1) [-].
S_D3_1
Dummy state 3 (reactor 1) [-].
X_D4_1
Dummy state 4 (reactor 1) [-].
X_D5_1
Dummy state 5 (reactor 1) [-].
S_I2
Soluble inert organic matter (reactor 2) [g(COD) ⋅ m⁻³].
S_S2
Readily biodegradable substrate (reactor 2) [g(COD) ⋅ m⁻³].
X_I2
Particulate inert organic matter (reactor 2) [g(COD) ⋅ m⁻³].
X_S2
Slowly biodegradable substrate (reactor 2) [g(COD) ⋅ m⁻³].
X_BH2
Active heterotrophic biomass (reactor 2) [g(COD) ⋅ m⁻³].
X_BA2
Active autotrophic biomass (reactor 2) [g(COD) ⋅ m⁻³].
X_P2
Particulate products arising from biomass decay (reactor 2) [g(COD) ⋅ m⁻³].
S_O2
Dissolved oxygen (reactor 2) [g(O₂) ⋅ m⁻³].
S_NO2
Nitrate and nitrite (reactor 2) [g(N) ⋅ m⁻³].
S_NH2
Ammonium plus ammonia nitrogen (reactor 2) [g(N) ⋅ m⁻³].
S_ND2
Soluble biodegradable organic nitrogen (reactor 2) [g(N) ⋅ m⁻³].
X_ND2
Particulate biodegradable organic nitrogen (reactor 2) [g(N) ⋅ m⁻³].
S_ALK2
Alkalinity (reactor 2) [mol(HCO₃⁻) ⋅ m⁻³].
TSS2
Total suspended solids (reactor 2) [g(TSS) ⋅ m⁻³].
Q2
Flow rate (reactor 2) [m³ ⋅ d⁻¹].
T2
Temperature (reactor 2) [°C].
S_D1_2
Dummy state 1 (reactor 2) [-].
S_D2_2
Dummy state 2 (reactor 2) [-].
S_D3_2
Dummy state 3 (reactor 2) [-].
X_D4_2
Dummy state 4 (reactor 2) [-].
X_D5_2
Dummy state 5 (reactor 2) [-].
S_I3
Soluble inert organic matter (reactor 3) [g(COD) ⋅ m⁻³].
S_S3
Readily biodegradable substrate (reactor 3) [g(COD) ⋅ m⁻³].
X_I3
Particulate inert organic matter (reactor 3) [g(COD) ⋅ m⁻³].
X_S3
Slowly biodegradable substrate (reactor 3) [g(COD) ⋅ m⁻³].
X_BH3
Active heterotrophic biomass (reactor 3) [g(COD) ⋅ m⁻³].
X_BA3
Active autotrophic biomass (reactor 3) [g(COD) ⋅ m⁻³].
X_P3
Particulate products arising from biomass decay (reactor 3) [g(COD) ⋅ m⁻³].
S_O3
Dissolved oxygen (reactor 3) [g(O₂) ⋅ m⁻³].
S_NO3
Nitrate and nitrite (reactor 3) [g(N) ⋅ m⁻³].
S_NH3
Ammonium plus ammonia nitrogen (reactor 3) [g(N) ⋅ m⁻³].
S_ND3
Soluble biodegradable organic nitrogen (reactor 3) [g(N) ⋅ m⁻³].
X_ND3
Particulate biodegradable organic nitrogen (reactor 3) [g(N) ⋅ m⁻³].
S_ALK3
Alkalinity (reactor 3) [mol(HCO₃⁻) ⋅ m⁻³].
TSS3
Total suspended solids (reactor 3) [g(TSS) ⋅ m⁻³].
Q3
Flow rate (reactor 3) [m³ ⋅ d⁻¹].
T3
Temperature (reactor 3) [°C].
S_D1_3
Dummy state 1 (reactor 3) [-].
S_D2_3
Dummy state 2 (reactor 3) [-].
S_D3_3
Dummy state 3 (reactor 3) [-].
X_D4_3
Dummy state 4 (reactor 3) [-].
X_D5_3
Dummy state 5 (reactor 3) [-].
S_I4
Soluble inert organic matter (reactor 4) [g(COD) ⋅ m⁻³].
S_S4
Readily biodegradable substrate (reactor 4) [g(COD) ⋅ m⁻³].
X_I4
Particulate inert organic matter (reactor 4) [g(COD) ⋅ m⁻³].
X_S4
Slowly biodegradable substrate (reactor 4) [g(COD) ⋅ m⁻³].
X_BH4
Active heterotrophic biomass (reactor 4) [g(COD) ⋅ m⁻³].
X_BA4
Active autotrophic biomass (reactor 4) [g(COD) ⋅ m⁻³].
X_P4
Particulate products arising from biomass decay (reactor 4) [g(COD) ⋅ m⁻³].
S_O4
Dissolved oxygen (reactor 4) [g(O₂) ⋅ m⁻³].
S_NO4
Nitrate and nitrite (reactor 4) [g(N) ⋅ m⁻³].
S_NH4
Ammonium plus ammonia nitrogen (reactor 4) [g(N) ⋅ m⁻³].
S_ND4
Soluble biodegradable organic nitrogen (reactor 4) [g(N) ⋅ m⁻³].
X_ND4
Particulate biodegradable organic nitrogen (reactor 4) [g(N) ⋅ m⁻³].
S_ALK4
Alkalinity (reactor 4) [mol(HCO₃⁻) ⋅ m⁻³].
TSS4
Total suspended solids (reactor 4) [g(TSS) ⋅ m⁻³].
Q4
Flow rate (reactor 4) [m³ ⋅ d⁻¹].
T4
Temperature (reactor 4) [°C].
S_D1_4
Dummy state 1 (reactor 4) [-].
S_D2_4
Dummy state 2 (reactor 4) [-].
S_D3_4
Dummy state 3 (reactor 4) [-].
X_D4_4
Dummy state 4 (reactor 4) [-].
X_D5_4
Dummy state 5 (reactor 4) [-].
S_I5
Soluble inert organic matter (reactor 5) [g(COD) ⋅ m⁻³].
S_S5
Readily biodegradable substrate (reactor 5) [g(COD) ⋅ m⁻³].
X_I5
Particulate inert organic matter (reactor 5) [g(COD) ⋅ m⁻³].
X_S5
Slowly biodegradable substrate (reactor 5) [g(COD) ⋅ m⁻³].
X_BH5
Active heterotrophic biomass (reactor 5) [g(COD) ⋅ m⁻³].
X_BA5
Active autotrophic biomass (reactor 5) [g(COD) ⋅ m⁻³].
X_P5
Particulate products arising from biomass decay (reactor 5) [g(COD) ⋅ m⁻³].
S_O5
Dissolved oxygen (reactor 5) [g(O₂) ⋅ m⁻³].
S_NO5
Nitrate and nitrite (reactor 5) [g(N) ⋅ m⁻³].
S_NH5
Ammonium plus ammonia nitrogen (reactor 5) [g(N) ⋅ m⁻³].
S_ND5
Soluble biodegradable organic nitrogen (reactor 5) [g(N) ⋅ m⁻³].
X_ND5
Particulate biodegradable organic nitrogen (reactor 5) [g(N) ⋅ m⁻³].
S_ALK5
Alkalinity (reactor 5) [mol(HCO₃⁻) ⋅ m⁻³].
TSS5
Total suspended solids (reactor 5) [g(TSS) ⋅ m⁻³].
Q5
Flow rate (reactor 5) [m³ ⋅ d⁻¹].
T5
Temperature (reactor 5) [°C].
S_D1_5
Dummy state 1 (reactor 5) [-].
S_D2_5
Dummy state 2 (reactor 5) [-].
S_D3_5
Dummy state 3 (reactor 5) [-].
X_D4_5
Dummy state 4 (reactor 5) [-].
X_D5_5
Dummy state 5 (reactor 5) [-].
YINIT1
YINIT1 = array([S_I1, S_S1, X_I1, X_S1, X_BH1, X_BA1, X_P1, S_O1, S_NO1, S_NH1, S_ND1, X_ND1, S_ALK1, TSS1, Q1, T1, S_D1_1, S_D2_1, S_D3_1, X_D4_1, X_D5_1])
Initial concentrations for the activated sludge reactor 1.
YINIT2
YINIT2 = array([S_I2, S_S2, X_I2, X_S2, X_BH2, X_BA2, X_P2, S_O2, S_NO2, S_NH2, S_ND2, X_ND2, S_ALK2, TSS2, Q2, T2, S_D1_2, S_D2_2, S_D3_2, X_D4_2, X_D5_2])
Initial concentrations for the activated sludge reactor 2.
YINIT3
YINIT3 = array([S_I3, S_S3, X_I3, X_S3, X_BH3, X_BA3, X_P3, S_O3, S_NO3, S_NH3, S_ND3, X_ND3, S_ALK3, TSS3, Q3, T3, S_D1_3, S_D2_3, S_D3_3, X_D4_3, X_D5_3])
Initial concentrations for the activated sludge reactor 3.
YINIT4
YINIT4 = array([S_I4, S_S4, X_I4, X_S4, X_BH4, X_BA4, X_P4, S_O4, S_NO4, S_NH4, S_ND4, X_ND4, S_ALK4, TSS4, Q4, T4, S_D1_4, S_D2_4, S_D3_4, X_D4_4, X_D5_4])
Initial concentrations for the activated sludge reactor 4.
YINIT5
YINIT5 = array([S_I5, S_S5, X_I5, X_S5, X_BH5, X_BA5, X_P5, S_O5, S_NO5, S_NH5, S_ND5, X_ND5, S_ALK5, TSS5, Q5, T5, S_D1_5, S_D2_5, S_D3_5, X_D4_5, X_D5_5])
Initial concentrations for the activated sludge reactor 5.
MU_H
Maximum heterotrophic growth rate [d⁻¹].
K_S
Substrate half-saturation coefficient for heterotrophic growth [g(COD) ⋅ m⁻³].
K_OH
Oxygen half-saturation coefficient for heterotrophic growth [g(O₂) ⋅ m⁻³].
K_NO
Nitrate half-saturation coefficient for anoxic heterotrophic growth [g(N) ⋅ m⁻³].
B_H
Heterotrophic decay rate [d⁻¹].
MU_A
Maximum autotrophic growth rate [d⁻¹].
K_NH
Ammonia half-saturation coefficient for autotrophic growth [g(N) ⋅ m⁻³].
K_OA
Oxygen half-saturation coefficient for autotrophic growth [g(O₂) ⋅ m⁻³].
B_A
Autotrophic decay rate [d⁻¹].
NY_G
Anoxic growth rate correction factor [-].
K_A
Ammonification rate [m³ ⋅ (g(COD) ⋅ d)⁻¹].
K_H
Maximum specific hydrolysis rate [g(COD) ⋅ (g(COD) ⋅ d)⁻¹].
K_X
Particulate substrate half-saturation coefficient for hydrolysis [g(COD) ⋅ g(COD)⁻¹].
NY_H
Anoxic hydrolysis rate correction factor [-].
Y_H
Heterotrophic yield [g(COD) ⋅ g(COD)⁻¹].
Y_A
Autotrophic yield [g(COD) ⋅ g(N)⁻¹].
F_P
Fraction of biomass leading to particulate inert products [-].
I_XB
Fraction of nitrogen in biomass [g(N) ⋅ g(COD)⁻¹].
I_XP
Fraction of nitrogen in organic particulate inerts [g(N) ⋅ g(COD)⁻¹].
X_I2TSS
Conversion factor for particulate inert organic matter to TSS [-].
X_S2TSS
Conversion factor for readily biodegradable substrate to TSS [-].
X_BH2TSS
Conversion factor for heterotrophic biomass to TSS [-].
X_BA2TSS
Conversion factor for autotrophic biomass to TSS [-].
X_P2TSS
Conversion factor for particulate products to TSS [-].
PAR1
PAR1 = array([MU_H, K_S, K_OH, K_NO, B_H, MU_A, K_NH, K_OA, B_A, NY_G, K_A, K_H, K_X, NY_H, Y_H, Y_A, F_P, I_XB, I_XP, X_I2TSS, X_S2TSS, X_BH2TSS, X_BA2TSS, X_P2TSS])
Parameters for the activated sludge reactor 1 at 15 °C, based on Alex et al (2018) (BSM1).
PAR2
Parameters for the activated sludge reactor 2 at 15 °C, based on Alex et al (2018) (BSM1).
PAR3
Parameters for the activated sludge reactor 3 at 15 °C, based on Alex et al (2018) (BSM1).
PAR4
Parameters for the activated sludge reactor 4 at 15 °C, based on Alex et al (2018) (BSM1).
PAR5
Parameters for the activated sludge reactor 5 at 15 °C, based on Alex et al (2018) (BSM1).
VOL1
Volume of reactor 1 [m³].
VOL2
Volume of reactor 2 [m³].
VOL3
Volume of reactor 3 [m³].
VOL4
Volume of reactor 4 [m³].
VOL5
Volume of reactor 5 [m³].
SOSAT1
Oxygen saturation concentration at 15 °C in reactor 1 [g(O₂) ⋅ m⁻³].
SOSAT2
Oxygen saturation concentration at 15 °C in reactor 2 [g(O₂) ⋅ m⁻³].
SOSAT3
Oxygen saturation concentration at 15 °C in reactor 3 [g(O₂) ⋅ m⁻³].
SOSAT4
Oxygen saturation concentration at 15 °C in reactor 4 [g(O₂) ⋅ m⁻³].
SOSAT5
Oxygen saturation concentration at 15 °C in reactor 5 [g(O₂) ⋅ m⁻³].