Skip to content

aerationcontrolinit_bsm1

Initialisation file for all states and parameters related to the aeration control system in reactor 3 to 5

All parameters and specifications are based on BSM1 model.

This file will be executed when running bsm1_cl.py.

CARB1_MAX

CARB1_MAX = 5

Maximum possible external carbon flow rate to reactor 1 [m³ ⋅ d⁻¹].

CARB2_MAX

CARB2_MAX = 5

Maximum possible external carbon flow rate to reactor 2 [m³ ⋅ d⁻¹].

CARB3_MAX

CARB3_MAX = 5

Maximum possible external carbon flow rate to reactor 3 [m³ ⋅ d⁻¹].

CARB4_MAX

CARB4_MAX = 5

Maximum possible external carbon flow rate to reactor 4 [m³ ⋅ d⁻¹].

CARB5_MAX

CARB5_MAX = 5

Maximum possible external carbon flow rate to reactor 5 [m³ ⋅ d⁻¹].

QINTR_MAX

QINTR_MAX = 5 * QIN0

Maximum pump capacity for internal recirculation rate Q_intr [m³ ⋅ d⁻¹].

QINTR_MIN

QINTR_MIN = 0

Minimum pump capacity for internal recirculation rate Q_intr [m³ ⋅ d⁻¹].

QW_MAX

QW_MAX = 0.1 * QIN0

Maximum pump capacity for waste sludge flow rate Q_w [m³ ⋅ d⁻¹].

QR_MAX

QR_MAX = 2 * QIN0

Maximum pump capacity for return sludge flow rate Q_r [m³ ⋅ d⁻¹].

QSTORAGE_MAX

QSTORAGE_MAX = 1500

Maximum pump capacity for outlet flow rate of the storage tank Q_st,out [m³ ⋅ d⁻¹].

QWT

QWT = 1 / 60 / 24 * 10

Time delay for artificial Q_w actuator (acts as first-order filter) [d].

T90_SNO2

T90_SNO2 = 10

Response time t_r for nitrate sensor [min].

MIN_SNO2

MIN_SNO2 = 0

Lower measuring limit of the nitrate sensor [g(N) ⋅ m⁻³].

MAX_SNO2

MAX_SNO2 = 10

Upper measuring limit of the nitrate sensor [g(N) ⋅ m⁻³].

T_SNO2

T_SNO2 = T90_SNO2 / 60 * 24 / 11.7724

Integral part time constant τ of transfer function for nitrate sensor [d].

STD_SNO2

STD_SNO2 = 0.025

Standard deviation for adding measurement noise for nitrate sensor [-].

KQINTR

KQINTR = 10000

Amplification constant for Q_intr controller [-].

TIQINTR

TIQINTR = 0.025

Integral part time constant τ for Q_intr controller [d].

TTQINTR

TTQINTR = 0.015

Integral part time constant τ of 'antiwindup' for Q_intr controller [d].

TDQINTR

TDQINTR = 0

Differential part time constant τ for Q_intr controller [d].

QINTRINTSTATE

QINTRINTSTATE = 0

Initial integration value for nitrate concentration of Q_intr controller [g(N) ⋅ m⁻³].

SNO2AWSTATE

SNO2AWSTATE = 0

Initial integration value for 'antiwindup' nitrate concentration of Q_intr controller [g(N) ⋅ m⁻³].

SNO2REF

SNO2REF = 1

Set point for nitrate concentration controller [g(N) ⋅ m⁻³].

QINTROFFSET

QINTROFFSET = 18500

Reasonable offset value for Q_intr controller around S_NO2^ref [m³ ⋅ d⁻¹].

USEANTIWINDUPSNO2

USEANTIWINDUPSNO2 = True

Boolean value to use 'antiwindup' for QINTR controller.

KFEEDFORWARD

KFEEDFORWARD = 0

Feedforward gain of QIN to QINTR [-].

T90_SO3

T90_SO3 = 1

Response time t_r for oxygen sensor 3 [min].

MIN_SO3

MIN_SO3 = 0

Lower measuring limit of the oxygen sensor 3 [g(O₂) ⋅ m⁻³].

MAX_SO3

MAX_SO3 = 10

"Upper measuring limit of the oxygen sensor 3 [g(O₂) ⋅ m⁻³].

T_SO3

T_SO3 = T90_SO3 / 60 * 24 / 3.89

Integral part time constant τ of transfer function for sensor 3 [d].

STD_SO3

STD_SO3 = 0.025

Standard deviation for adding measurement noise for sensor 3 [-].

KLA3_INIT

KLA3_INIT = 201.3015557168598

Initial value for PI controller 3 [d⁻¹].

KLA3_MIN

KLA3_MIN = 0

Lower limit of the adjustable KLa value for PI controller 3 [d⁻¹].

KLA3_MAX

KLA3_MAX = 360

Upper limit of the adjustable KLa value for PI controller 3 [d⁻¹].

KSO3

KSO3 = 25

Amplification constant for PI controller 3 [-].

TISO3

TISO3 = 0.002

Integral part time constant τ for PI controller 3 [d].

TTSO3

TTSO3 = 0.001

Integral part time constant τ of 'antiwindup' for PI controller 3 [d].

TDSO3

TDSO3 = 0

Differential part time constant τ for PI controller 3 [d].

SO3INTSTATE

SO3INTSTATE = 0

Initial integration value for saturated oxygen concentration of PI controller 3 [g(O₂) ⋅ m⁻³].

SO3AWSTATE

SO3AWSTATE = 0

Initial integration value for 'antiwindup' saturated oxygen concentration of PI controller 3
[g(O₂) ⋅ m⁻³].

SO3REF

SO3REF = 2

Set point for oxygen concentration controller 3 [g(O₂) ⋅ m⁻³].

KLA3OFFSET

KLA3OFFSET = 144

Reasonable offset value for PI controller 3 around S_O^ref [d⁻¹].

KLA3_LIM

KLA3_LIM = 0

KLa value after adjusting to upper and lower limit for PI controller 3 [d⁻¹].

KLA3_CALC

KLA3_CALC = 0

KLa value calculated from PI control [d⁻¹].

USEANTIWINDUPSO3

USEANTIWINDUPSO3 = True

Boolean value to use 'antiwindup' for PI controller 3.

T90_KLA3

T90_KLA3 = 4

Response time t_r for KLa actuator 3 [min].

T_KLA3

T_KLA3 = T90_KLA3 / 60 * 24 / 3.89

Integral part time constant τ for KLa actuator 3 [d].

T90_SO4

T90_SO4 = 1

Response time t_r for oxygen sensor 4 [min].

MIN_SO4

MIN_SO4 = 0

Lower measuring limit of the oxygen sensor 4 [g(O₂) ⋅ m⁻³].

MAX_SO4

MAX_SO4 = 10

"Upper measuring limit of the oxygen sensor 4 [g(O₂) ⋅ m⁻³].

T_SO4

T_SO4 = T90_SO4 / 60 * 24 / 3.89

Integral part time constant τ of transfer function for sensor 4 [d].

STD_SO4

STD_SO4 = 0.025

Standard deviation for adding measurement noise for sensor 4 [-].

KLA4_INIT

KLA4_INIT = 201.3015557168598

Initial value for PI controller 4 [d⁻¹].

KLA4_MIN

KLA4_MIN = 0

Lower limit of the adjustable KLa value for PI controller 4 [d⁻¹].

KLA4_MAX

KLA4_MAX = 360

Upper limit of the adjustable KLa value for PI controller 4 [d⁻¹].

KSO4

KSO4 = 25

Amplification constant for PI controller 4 [-].

TISO4

TISO4 = 0.002

Integral part time constant τ for PI controller 4 [d].

TTSO4

TTSO4 = 0.001

Integral part time constant τ of 'antiwindup' for PI controller 4 [d].

TDSO4

TDSO4 = 0

Differential part time constant τ for PI controller 4 [d].

SO4INTSTATE

SO4INTSTATE = 0

Initial integration value for saturated oxygen concentration of PI controller 4 [g(O₂) ⋅ m⁻³].

SO4AWSTATE

SO4AWSTATE = 0

Initial integration value for 'antiwindup' saturated oxygen concentration of PI controller 4
[g(O₂) ⋅ m⁻³].

SO4REF

SO4REF = 2

Set point for oxygen concentration controller 4 [g(O₂) ⋅ m⁻³].

KLA4OFFSET

KLA4OFFSET = 144

Reasonable offset value for PI controller 4 around S_O^ref [d⁻¹].

KLA4_LIM

KLA4_LIM = 0

KLa value after adjusting to upper and lower limit for PI controller 4 [d⁻¹].

KLA4_CALC

KLA4_CALC = 0

KLa value calculated from PI control [d⁻¹].

USEANTIWINDUPSO4

USEANTIWINDUPSO4 = True

Boolean value to use 'antiwindup' for PI controller 4.

T90_KLA4

T90_KLA4 = 4

Response time t_r for KLa actuator 4 [min].

T_KLA4

T_KLA4 = T90_KLA4 / 60 * 24 / 3.89

Integral part time constant τ for KLa actuator 4 [d].

T90_SO5

T90_SO5 = 1

Response time t_r for oxygen sensor 5 [min].

MIN_SO5

MIN_SO5 = 0

Lower measuring limit of the oxygen sensor 5 [g(O₂) ⋅ m⁻³].

MAX_SO5

MAX_SO5 = 10

"Upper measuring limit of the oxygen sensor 5 [g(O₂) ⋅ m⁻³].

T_SO5

T_SO5 = T90_SO5 / 60 * 24 / 3.89

Integral part time constant τ of transfer function for sensor 5 [d].

STD_SO5

STD_SO5 = 0.025

Standard deviation for adding measurement noise for sensor 5 [-].

KLA5_INIT

KLA5_INIT = 201.3015557168598

Initial value for PI controller 5 [d⁻¹].

KLA5_MIN

KLA5_MIN = 0

Lower limit of the adjustable KLa value for PI controller 5 [d⁻¹].

KLA5_MAX

KLA5_MAX = 360

Upper limit of the adjustable KLa value for PI controller 5 [d⁻¹].

KSO5

KSO5 = 25

Amplification constant for PI controller 5 [-].

TISO5

TISO5 = 0.002

Integral part time constant τ for PI controller 5 [d].

TTSO5

TTSO5 = 0.001

Integral part time constant τ of 'antiwindup' for PI controller 5 [d].

TDSO5

TDSO5 = 0

Differential part time constant τ for PI controller 5 [d].

SO5INTSTATE

SO5INTSTATE = 0

Initial integration value for saturated oxygen concentration of PI controller 5 [g(O₂) ⋅ m⁻³].

SO5AWSTATE

SO5AWSTATE = 0

Initial integration value for 'antiwindup' saturated oxygen concentration of PI controller 5
[g(O₂) ⋅ m⁻³].

SO5REF

SO5REF = 2

Set point for oxygen concentration controller 5 [g(O₂) ⋅ m⁻³].

KLA5OFFSET

KLA5OFFSET = 144

Reasonable offset value for PI controller 5 around S_O^ref [d⁻¹].

KLA5_LIM

KLA5_LIM = 0

KLa value after adjusting to upper and lower limit for PI controller 5 [d⁻¹].

KLA5_CALC

KLA5_CALC = 0

KLa value calculated from PI control [d⁻¹].

USEANTIWINDUPSO5

USEANTIWINDUPSO5 = True

Boolean value to use 'antiwindup' for PI controller 5.

T90_KLA5

T90_KLA5 = 4

Response time t_r for KLa actuator 5 [min].

T_KLA5

T_KLA5 = T90_KLA5 / 60 * 24 / 3.89

Integral part time constant τ for KLa actuator 5 [d].