Model Mathematics

Component: environment

Component: E

ddtime E = k2 ⁢ ENa + k4 ⁢ EH + k6 ⁢ ENH4 - k1 ⁢ Na ⁢ E + k3 ⁢ H ⁢ E + k5 ⁢ NH4 ⁢ E

Component: ENa

ddtime ENa = k1 ⁢ E ⁢ Na + k14 ⁢ ENa_ - k2 ⁢ ENa + k13 ⁢ ENa

Component: EH

ddtime EH = k3 ⁢ H ⁢ E + k16 ⁢ EH_ - k4 ⁢ EH + k15 ⁢ EH

Component: ENH4

ddtime ENH4 = k5 ⁢ NH4 ⁢ E + k18 ⁢ ENH4_ - k6 ⁢ ENH4 + k17 ⁢ ENH4

Component: E_

ddtime E_ = k8 ⁢ ENa_ + k10 ⁢ EH_ + k12 ⁢ ENH4_ - k7 ⁢ Na_ ⁢ E_ + k9 ⁢ H_ ⁢ E_ + k11 ⁢ NH4_ ⁢ E_

Component: ENa_

ddtime ENa_ = k7 ⁢ E_ ⁢ Na_ + k13 ⁢ ENa - k8 ⁢ ENa_ + k14 ⁢ ENa_

Component: EH_

ddtime EH_ = k9 ⁢ H_ ⁢ E_ + k15 ⁢ EH - k10 ⁢ EH_ + k16 ⁢ EH_

Component: ENH4_

ddtime ENH4_ = k11 ⁢ NH4_ ⁢ E_ + k17 ⁢ ENH4 - k12 ⁢ ENH4_ + k18 ⁢ ENH4_

Component: reaction_constants

W = H_ H_ + KH k13 = 2.41E0 ⁢ W k14 = 2.41E0 ⁢ W k15 = 4.8E-1 ⁢ W k16 = 4.8E-1 ⁢ W k17 = 2.41E0 ⁢ W k18 = 2.41E0 ⁢ W