| Description |
|---|
| Field Noyes Version of Belousov- Zhabotinsky Reaction. BrO3 is held constant; HOBr is typically ignored, and can be replaced by an empty- set. The stoichiometry f is typically taken as 1/ 2 or 1 (denominator 1 or 2 in SBML) . In reaction5 it may be necessary to manually replace the denominator= "1/(f)" with denominator="2" or denominator= "1" to get f=1/2 or 1, respectively, because in standard SBML these should be integers and not variables. Some simulators may allow variable denominators, however. |
| Rate constant | Reaction |
|---|---|
| k1 = 1.3 | Br + BrO3 -> HBrO2 + HOBr |
| k2 = 2000000 | Br + HBrO2 -> HOBr^2 |
| k3 = 34 | BrO3 + HBrO2 -> Ce^2 + HBrO2^2 |
| k4 = 3000 | HBrO2^2 -> BrO3 + HOBr |
| k5 = 0.02 | Ce -> Br^f |
| Variable | IC | ODE |
|---|---|---|
| Br | 0.003 | Br'[t] == -(k1*Br[t]*BrO3[t]) + f*k5*Ce[t] - k2*Br[t]*HBrO2[t] |
| Ce | 0.05 | Ce'[t] == -(k5*Ce[t]) + 2*k3*BrO3[t]*HBrO2[t] |
| HBrO2 | 0.001 | HBrO2'[t] == k1*Br[t]*BrO3[t] - k2*Br[t]*HBrO2[t] + k3*BrO3[t]*HBrO2[t] - k4*HBrO2[t]^2 |
| HOBr | 0 | HOBr'[t] == k1*Br[t]*BrO3[t] + 2*k2*Br[t]*HBrO2[t] + k4*HBrO2[t]^2 |