From c2122e95372bb268f2902d42bf5323e306b4346a Mon Sep 17 00:00:00 2001 From: "p.vanderwilt" Date: Wed, 4 Jun 2025 17:28:52 +0200 Subject: [PATCH] Implemented rates in rates function --- dependencies/asm3_class.js | 42 ++++++++++++++++++++++++++++++++------ 1 file changed, 36 insertions(+), 6 deletions(-) diff --git a/dependencies/asm3_class.js b/dependencies/asm3_class.js index 2b85fe5..1b214db 100644 --- a/dependencies/asm3_class.js +++ b/dependencies/asm3_class.js @@ -5,7 +5,7 @@ class ASM3 { // 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] + 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 [-] @@ -16,16 +16,16 @@ class ASM3 { 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_O2: 0.2, // aerobic respiration rate [d-1] + b_H_O: 0.2, // aerobic respiration rate [d-1] b_H_NO: 0.1, // anoxic respiration rate [d-1] - b_STO_O2: 0.2, // aerobic respitation rate X_STO [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_O2: 0.15, // aerobic respiration rate [d-1] + b_A_O: 0.15, // aerobic respiration rate [d-1] b_A_NO: 0.05 // anoxic respiration rate [d-1] } @@ -51,13 +51,43 @@ class ASM3 { 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() { + 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); - rates[0] = this.parameters[] + 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; } } \ No newline at end of file