93 lines
4.8 KiB
JavaScript
93 lines
4.8 KiB
JavaScript
class ASM3 {
|
|
|
|
kin_params = {
|
|
// Kinetic parameters (20 C for now)
|
|
|
|
// Hydrolysis
|
|
k_H: 3., // hydrolysis rate constant [g X_S g-1 X_H d-1]
|
|
K_X: 1., // hydrolysis saturation constant [g X_S g-1 X_H]
|
|
// Heterotrophs
|
|
k_STO: 5., // storage rate constant [g S_S g-1 X_H d-1]
|
|
nu_NO: 0.6, // anoxic reduction factor [-]
|
|
K_O: 0.2, // saturation constant S_0 [g O2 m-3]
|
|
K_NO: 0.5, // saturation constant S_NO [g NO3-N m-3]
|
|
K_S: 2., // saturation constant S_s [g COD m-3]
|
|
K_STO: 1., // saturation constant X_STO [g X_STO g-1 X_H]
|
|
mu_H_max: 2., // maximum specific growth rate [d-1]
|
|
K_NH: 0.01, // saturation constant S_NH3 [g NH3-N m-3]
|
|
K_HCO: 0.1, // saturation constant S_HCO [mole HCO3 m-3]
|
|
b_H_O: 0.2, // aerobic respiration rate [d-1]
|
|
b_H_NO: 0.1, // anoxic respiration rate [d-1]
|
|
b_STO_O: 0.2, // aerobic respitation rate X_STO [d-1]
|
|
b_STO_NO: 0.1, // anoxic respitation rate X_STO [d-1]
|
|
// Autotrophs
|
|
mu_A_max: 1.0, // maximum specific growth rate [d-1]
|
|
K_A_NH: 1., // saturation constant S_NH3 [g NH3-N m-3]
|
|
K_A_O: 0.5, // saturation constant S_0 [g O2 m-3]
|
|
K_A_HCO: 0.5, // saturation constant S_HCO [mole HCO3 m-3]
|
|
b_A_O: 0.15, // aerobic respiration rate [d-1]
|
|
b_A_NO: 0.05 // anoxic respiration rate [d-1]
|
|
}
|
|
|
|
stoi_params = {
|
|
// Stoichiometric and composition parameters
|
|
|
|
f_S_I: 0., // fraction S_I from hydrolysis [g S_I g-1 X_S]
|
|
// Yields
|
|
Y_STO_O: 0.85, // Aerobic yield X_STO per S_S [g X_STO g-1 S_S]
|
|
Y_STO_NO: 0.80, // Anoxic yield X_STO per S_S [g X_STO g-1 S_S]
|
|
Y_H_O: 0.63, // Aerobic yield X_H per X_STO [g X_H g-1 X_STO]
|
|
Y_H_NO: 0.54, // Anoxic yield X_H per X_STO [g X_H g-1 X_STO]
|
|
Y_A: 0.24, // anoxic yield X_A per S_NO [g X_A g-1 NO3-N]
|
|
// Composition (nitrogen)
|
|
i_NSI: 0.01, // nitrogen content S_I [g N g-1 S_I]
|
|
i_NSS: 0.03, // nitrogen content S_S [g N g-1 S_S]
|
|
i_NXI: 0.02, // nitrogen content X_I [g N g-1 X_I]
|
|
i_NXS: 0.04, // nitrogen content X_S [g N g-1 X_S]
|
|
i_NBM: 0.07, // nitrogen content X_H / X_A [g N g-1 X_H / X_A]
|
|
// Composition (TSS)
|
|
i_TS: 0.75, // TSS content X_I [g TS g-1 X_I]
|
|
i_TS: 0.75, // TSS content X_S [g TS g-1 X_S]
|
|
i_TS: 0.90, // TSS content X_H / X_A [g TS g-1 X_H / X_A]
|
|
i_TS: 0.60 // TSS content X_STO (PHB based) [g TS g-1 X_STO]
|
|
}
|
|
// S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
|
|
state = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
|
|
|
|
constructor(state) {
|
|
this.state = state;
|
|
}
|
|
|
|
_monod(c, K){
|
|
return c / (K + c);
|
|
}
|
|
|
|
_inv_monod(c, K){
|
|
return K / (K + c);
|
|
}
|
|
|
|
compute_rates(state) {
|
|
const rates = new Array(12);
|
|
const [S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS] = state;
|
|
|
|
// Hydrolysis
|
|
rates[0] = this.kin_params.k_H * this._monod(X_S / X_H, this.kin_params.K_X) * X_H;
|
|
|
|
// Heterotrophs
|
|
rates[1] = this.kin_params.k_STO * this._monod(S_O, this.kin_params.K_O) * this._monod(S_S, this.kin_params.K_S) * X_H;
|
|
rates[2] = this.kin_params.k_STO * this.kin_params.nu_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * this._monod(S_S, this.kin_params.K_S) * X_H;
|
|
rates[3] = this.kin_params.mu_H_max * this._monod(S_O, this.kin_params.K_O) * this._monod(S_NH, this.kin_params.K_NH) * this._monod(S_HCO, this.kin_params.K_HCO) * this._monod(X_STO/X_H, K_STO) * X_H;
|
|
rates[4] = this.kin_params.mu_H_max * this.kin_params.nu_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.NO) * this._monod(S_NH, this.kin_params.K_NH) * this._monod(S_HCO, this.kin_params.K_HCO) * this._monod(X_STO/X_H, this.kin_params.K_STO) * X_H;
|
|
rates[5] = this.kin_params.b_H_O * this._monod(S_O, this.kin_params.K_O) * X_H;
|
|
rates[6] = this.kin_params.b_H_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * X_H;
|
|
rates[7] = this.kin_params.b_STO_O * this._monod(S_O, this.kin_params.K_O) * X_H;
|
|
rates[8] = this.kin_params.b_STO_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * X_STO;
|
|
|
|
// Autotrophs
|
|
rates[9] = this.kin_params.mu_A_max * this._monod(S_O, this.kin_params.K_A_O) * this._monod(S_NH, this.kin_params.K_A_NH) * this._monod(S_HCO, this.kin_params.K_A_HCO) * X_A;
|
|
rates[10] = this.kin_params.b_A_O * this._monod(S_O, this.kin_params.K_O) * X_A;
|
|
rates[11] = this.kin_params.b_A_NO * this._inv_monod(S_O, this.kin_params.K_A_O) * this._monod(S_NO, this.kin_params.K_NO) * X_A;
|
|
|
|
return rates;
|
|
}
|
|
} |